MODELIRANJE ELEKTROMAGNETNEGA SPROŽNIKA S KRATKOSTIČNIM OBROČKOM

Size: px
Start display at page:

Download "MODELIRANJE ELEKTROMAGNETNEGA SPROŽNIKA S KRATKOSTIČNIM OBROČKOM"

Transcription

1 UNIVERZA V LJUBLJANI Fakulteta za elektrotehniko Borut Drnovšek MODELIRANJE ELEKTROMAGNETNEGA SPROŽNIKA S KRATKOSTIČNIM OBROČKOM MAGISTRSKO DELO Mentor: izr. prof. dr. Dejan Križaj Ljubljana, 2013

2 II ZAHVALA Zahvaljujem se mentorju izr. prof. dr. Dejanu Križaju, ki me je s svojim znanjem in nasveti usmerjal in me podpiral pri izdelavi magistrske naloge. Zahvaljujem se tudi Topčagić Zumretu, ki mi je s svojimi nasveti pomagal, da sem prišel do hitrejših rezultatov simulacij. Zahvala velja tudi mojim staršem, ki so mi omogočili študij in me skozi vsa leta moralno in drugače podpirali. Posebna zahvala pa gre Špeli in mojemu sinčku Žanu, ki sta me navdihnila, da sem končno uspel dokončati magistrsko delo. Borut Drnovšek

3 III VSEBINA 1 POVZETEK UVOD SPLOŠNO O INŠTALACIJSKEM ODKLOPNIKU ZGRADBA INŠTALACIJSKEGA ODKLOPNIKA SELEKTIVNI ODKLOPNIK RAZRED SELEKTIVNOSTI INŠTALACIJSKIH ODKLOPNIKOV Popolna selektivnost Delna selektivnost ZGRADBA IN DELOVANJE SELEKTIVNEGA ODKLOPNIKA Zgradba selektivnega odklopnika Delovanje selektivnega odklopnika ZGRADBA IN DELOVANJE ELEKTROMAGNETNEGA SPROŽNIKA Zgradba elektromagnetnega sprožnika Delovanje elektromagnetnega sprožnika METODE ZA NUMERIČNI IZRAČUN MAGNETNEGA POLJA IN SIL V ELEKTROMAGNETNEM SPROŽNIKU OSNOVNE ENAČBE ELEKTROMAGNETNEGA POLJA ROBNI POGOJI MAXWELLOVA ENAČBA V KROŽNO VALJNIH KOORDINATAH SILA NA KOTVO ELEKTROMAGNETNEGA SPROŽNIKA NUMERIČNE SIMULACIJE ELEKTROMAGNETNEGA SPROŽNIKA SELEKTIVNEGA INŠTALACIJSKEGA ODKLOPNIKA NUMERIČNE SIMULACIJE ELEKTROMAGNETNEGA SPROŽNIKA BREZ KRATKOSTIČNEGA OBROČKA IN S KONSTANTNIM R Izmenični tok Vpliv zračne reže med kotvo in jedrom na gostoto magnetnega pretoka Velikost sile med jedrom in kotvo v odvisnosti od velikosti zračne reže med njima ELEKTROMAGNETNI SPROŽNIK BREZ KRATKOSTIČNEGA OBROČKA IN Z UPOŠTEVANJEM RELATIVNE PERMEABILNOSTI ŽELEZA Vpliv toka na gostoto magnetnega pretoka v sprožniku Vpliv toka na velikost magnetne sile med kotvo in jedrom ELEKTROMAGNETNI SPROŽNIK Z BAKRENIM KRATKOSTIČNIM OBROČKOM IN Z UPOŠTEVANJEM RELATIVNE PERMEABILNOSTI ŽELEZA Vpliv toka na gostoto magnetnega pretoka v sprožniku Vpliv vzbujalnega toka na inducirani tok v bakrenem obročku Vpliv toka na velikost magnetne sile med kotvo in jedrom z vgrajenim kratkostičnim bakrenim obročkom v jedru sprožnika Vpliv položaja vgradnje bakrenega kratkostičnega obročka na velikost sile med kotvo in jedrom Položaj obročka na zunanji strani jedra Položaj obročka na notranji strani jedra ELEKTROMAGNETNI SPROŽNIK Z BAKRENIM KRATKOSTIČNIM OBROČKOM, Z UPOŠTEVANJEM RELATIVNE PERMEABILNOSTI ŽELEZA, VENDAR BREZ MAGNETNEGA JARMA ZAKLJUČEK LITERATURA...71

4 4 Uporabljeni simboli I N - nazivni tok H - magnetna poljska jakost B - gostota magnetnega pretoka - permeabilnost - frekvenca - dielektričnost - prostorska gostota elektrine σ - specifična električna prevodnost N - število ovojev tuljave - magnetni pretok F - magnetna sila G električna prevodnost J gostota toka E električna poljska jakost v hitrost naelektrenih delcev K ploskovni tok A vektorski magnetni potencial D gostota električnega polja V električni potencial

5 5 1 Povzetek Namen magistrske naloge je numerična simulacija magnetnega polja in sil elektromagnetnega sprožnika, ki ima vgrajen kratkostični bakreni obroček. Ta elektromagnetni sprožnik je sestavni del selektivnega inštalacijskega odklopnika, ki mora zagotoviti zakasnjen izklop v primerjavi z naknadno vgrajenimi inštalacijskimi odklopniki, s čimer preprečimo izpad celotnega sistema. Selektivnost v kratkostičnem delu izklopne krivulje se zagotovi z vgrajenim preduporom, ki omeji tok skozi selektivni inštalacijski odklopnik in posledično zakasni izklop aparata. S selektivnim inštalacijskim odklopnikom lahko nadomestimo klasične talilne varovalke tipa NV, ki se uporabljajo v priključnih omaricah. Vse simulacije s programom COMSOL, ki deluje po metodi končnih elementov, so bile narejene v osno-simetričnem sistemu, kjer smo strukturo rotirali okoli»z«osi. Uporaba osne simetrije nam je izredno olajšala izdelavo same strukture ter povečala hitrost samega reševanja problema. Osno-simetričen sistem smo lahko uporabili, ker je elektromagnetni sprožnik skoraj v celoti osno-simetričen, saj le magnetni jarem izstopa iz osne simetrije, kar pa smo rešili s prilagoditvijo preseka samega jarma. Vse simulacije so bile narejene pri vzbujanju s harmoničnim signalom. Enostavnejša je geometrija strukture, lažji in hitrejši je izračun s programom COMSOL. Z večanjem števila končnih elementov se izboljšujejo rezultati analize, vendar smo tu omejeni z zmogljivostjo računalnika. Pri zelo velikem številu elementov se pogosto tudi oteži konvergiranje numeričnega izračuna k rešitvi. Poleg gostote mreže pa moramo biti previdni tudi na velikost okolice, saj z zunanjimi mejami okolice definiramo področje magnetnega polja. Najprej je v nalogi obdelan model elektromagnetnega sprožnika brez kratkostičnega obročka ter s konstantno relativno permeabilnostjo železa. V tem delu naloge smo s pomočjo simulacij ugotavljali vpliv velikosti zračne reže med kotvo in jedrom na gostoto magnetnega pretoka v elektromagnetnem sprožniku, posledično pa tudi na velikost sile med kotvo in jedrom. V drugem delu so prikazane simulacije na elektromagnetnem sprožniku, kjer smo upoštevali magnetilno krivuljo feromagnetnih elementov elektromagnetnega sprožnika. Še vedno pa je bila to struktura brez bakrenega kratkostičnega obročka. V tem delu naloge smo preverjali vpliv vzbujalnega toka na gostoto magnetnega pretoka v elektromagnetnem sprožniku ter

6 6 posledično velikost sile med kotvo in jarmom. Hkrati je v tem delu naloge prikazan tudi prehod sestavnih delov sprožnika v področje nasičenja, ko je vzbujalni tok dovolj velik. V tretjem delu naloge je prikazan vpliv bakrenega kratkostičnega obročka, ki je vgrajen v jedro elektromagnetnega sprožnika. Vpliv obročka je največji v času, ko je vrednost vzbujalnega sinusnega signala toka minimalna, saj je takrat pri sinusnem vzbujanju največja sprememba fluksa skozi obroček in s tem posledično največja inducirana napetost v obročku. Zaradi induktivnih in uporovnih lastnosti obročka je tok v obročku fazno zamaknjen glede na vzbujalni tok. Zaradi tega pojava se pojavi dodatna sila med kotvo in jarmom, ki je različna od nič tudi v času, ko je vzbujalni tok enak nič. Na ta način nam bakreni obroček preprečuje vibriranje gibljivega kontakta, saj nam zagotavlja stalno silo med jedrom in kotvo. Ker pa je tako velikost kot oblika induciranega toka odvisna od položaja obročka, je v tem delu naloge prikazana tudi razlika med obema skrajnima legama vgradnje bakrenega kratkostičnega obročka. Zaradi problemov s prostorom v aparatu je v zadnjem delu naloge prikazana še zgradba elektromagnetnega sprožnika, kjer je odstranjen magnetni jarem, s čimer smo preverili možne smeri optimizacije elektromagnetnega sprožnika. Ključne besede: selektivni odklopnik, elektromagnetni sprožnik s kratkostičnim obročkom, magnetni pretok, magnetna sila, metoda končnih elementov, COMSOL, histerezna zanka. Abstract The objective of the thesis is numerical simulation of magnetic field and electromagnetic forces in electromagnetically actuated magnet tripping unit with a built-in copper short-circuit ring. This tripping unit is an integral part of a Selective circuit breaker (SMCB) which should provide a time delayed switch off, compared to subsequently build in MCB in order to prevent breakage of the entire system. Selectivity in the short-circuit breaking part of the current signal is provided by the built-in resistor which limits the current through the selective circuit breaker and consequently delays the switch-off of the selective circuit breaker. With a

7 7 selective MCB we could replace a conventional fuse type NV that are nowadays used in the distribution box. The simulations have been performed using a numerical simulation program COMSOL Multiphysics which is based on a finite element method. Due to largely (but not completely) axial symmetric shape of the circuit breaker we have decided to design a simulation structure assuming completely axially symmetrical structure. The use of axial symmetry significantly simplifies the simulation process as it enables use of 2D simulation instead of the threedimensional one. The nonsymmetrical magnetic yoke has been modeled by adjusting the material parameters of the simulated yoke to approximate the magnetic resistance of the yoke. All simulations were performed with harmonic signals. This could in principal lead to another simplifications of usage of complex description of Ampere s law, however, for more realistic modeling we took into account also the magnetic nonlinearities which necessitates use of time domain modeling. By increasing the number of finite elements we get more accurate analysis results, but eventually we are limited by the memory capacity of the computer. Besides the density of the mesh we need to be careful also on the size of the complete simulation area, since the boundary conditions on external borders cannot be absolutely accurately defined. The Results section (Chapter 5) is separated to several subchapters. First we analyzed the model of an electromagnetic tripping unit without a short-circuit ring and with constant relative permeability of the iron parts. This enabled basic study of the influence of the size of the air gap between the anchor and the core on the magnetic flux density in the electromagnetic tripping unit and consequently on the magnitude of the force between the anchor and the core. As expected the force is following the sinusoidal excitation current with doubled frequency and is reaching zero at zero flux (zero current excitation). In the second part we analyzed the influence of the magnetization curve on the development of the force in the tripping unit. The magnetization of the iron core significantly influences the magnitude as well as the shape of the force between the anchor and the yoke in particular at larger current excitations. In the third part we analyzed the influence of a short-circuit copper ring which is built into the core of the electromagnetic trigger. During device operation the time variation of the flux through the ring results in induced voltage in the ring which drives the (short-circuit) current through the ring. This current builds its own magnetic field around the ring that superimposes onto the primary one. The influence of the ring is largest when the value of the current is zero as at that moment the change of the flux and the induced voltage through the

8 8 ring is largest. Due to inductive and resistive properties of the ring the induced current in the ring is phase shifted in comparison with the driving current through the coil. This results in an additional force that acts between the core and the yoke also during the time of zero excitation current. In this way the copper ring reduces vibrations of the movable contact since it ensures a non-zero force between the core and the anchor also during the time the excitation signal is zero. Simulations reveal that the size and the shape of the induced current depend on the position of the ring. In the last part of the work we analyzed the possibility of removal of the yoke as in some cases it would be advantageous to get additional space around the trigger. The simulations show that in that case the magnitude of the force between the anchor and the yoke would be significantly reduced. We have shown that numerical simulation can be a valuable tool for analyzing behavior of electromagnetic structures such as an electromagnetic tripping unit. Simulations can be used for improved understanding of device operation but also for verification of some improvements of design and even optimizations of device operation. Key-words: selective miniature circuit breaker, magnetic tripping unit with short circuit ring, magnetic flux, magnetic force, finite element method, COMSOL, hysteresis loop.

9 9 2 Uvod Pri uporabi električne energije je pomembno varovanje električnih strojev, vodnikov in uporabnikov. Zaradi napak, ki se pojavljajo v električnih strojih ali vodnikih, zaradi udara strele kamorkoli v prenos električne energije se pojavljajo večji tokovi, kot pa jih zahteva porabnik. Posledice teh napak so uničeni vodniki, pokvarjeni ali uničeni električni stroji in, kar je najpomembneje, ogroženo je zdravje in življenje ljudi in živali. Da se izognemo tem posledicam, vključujemo v električne tokokroge električne varovalke kot zaščito pred preobremenitvami. Električne varovalke so namenoma oslabljena mesta v tokokrogu, kjer se tokokrog prekine, kadar tok v tem tokokrogu prekorači določeno vrednost v določenem času. Električne varovalke, ki se uporabljajo v industriji in gospodinjstvu, se delijo na dve vrsti: - varovalke s talilnimi vložki (talilne varovalke), - inštalacijski odklopniki oz. avtomatske varovalke. Telo talilnega vložka je iz kvalitetnega steatita, zelo odpornega proti temperaturnim preobremenitvam. Pokrova sta iz aluminija in sta odporna proti učinkom korozivne atmosfere. V notranjem delu keramičnega telesa je nameščen bakreni talilni element, ki je točkasto privarjen na posebno oblikovan notranji del kontaktnega noža. Preostanek notranjosti je zasut s kremenčevim peskom točno določene granulacije in sestave. Na mestu prekinitve se pojavi oblok, ki upari bakreni talilni element in tali pesek. Hkrati se dviguje tlak v obločnem kanalu. Visok tlak požene bakrene pare v okoliški pesek, tako da v obločnem kanalu ni več kovinskih delcev, stene obločnega kanala pa so izključno iz neprevodnih materialov. Taljenje peska hkrati tudi ohlaja obločni plamen. Omenjena pojava dvigujeta obločno napetost. Ko obločna napetost preseže pritisnjeno napetost, je izpolnjen pogoj za uspešno omejevanje in posledično tudi prekinitev toka [1, stran 3]. 2.1 Splošno o inštalacijskem odklopniku Inštalacijski odklopnik ali MCB (ang. Miniature Circuit Breaker) je naprava, ki služi za zaščito pred preobremenitvami in kratkimi stiki. Večina inštalacijskih odklopnikov je zgrajena tako, da v preobremenitvenem delu izklaplja bimetal, saj morajo biti časi izklopa v tem primeru napake dokaj dolgi. Časi se gibljejo od nekaj sekund pa vse do 1 ure.

10 10 V primeru pojava kratkega stika pa izklaplja elektromagnetni sprožnik, ki zagotavlja izredno hiter izklop aparata, saj so kratkostični tokovi zelo veliki in lahko dosežejo tudi nekaj ka. Inštalacijski odklopniki se ločijo glede na izklopno karakteristiko (karakteristika proženja), kjer je definiran čas, v katerem mora aparat pri določeni vrednosti toka izklopiti. Glede na izklopno karakteristiko razlikujemo naslednje tipe odklopnikov: - B (za zaščito električnih vodnikov, v gospodinjstvu, ), - C (za zaščito naprav, ki imajo večje zagonske tokove, npr.: elektromotorji), - D (povsod tam, kjer se pojavljajo izjemno velike kratkotrajne tokovne konice ob vklopih (npr.: določeni motorji, transformatorji, halogenske razsvetljave), da ob vklopih takih naprav ne pride do nepotrebnih izklopov). Glavna razlika med temi tremi karakteristikami je v začetku proženja elektromagneta, saj pri B tipu elektromagnet začne delovati v področju med 3- in 5-kratnikom nazivnega toka, pri C tipu začne elektromagnet delovati v področju med 5- in 10-kratnikom nazivnega toka, pri D karakteristiki pa med 10- in 20-kratnikom nazivnega toka. Obstaja še karakteristika K, ki pa se uporablja zelo redko, in sicer je ta karakteristika nekakšen približek motorskemu zaščitnemu stikalu, saj začne elektromagnet delovati med 8- in 12-kratnikom nazivnega toka. Največ se uporabljajo odklopniki tipa B in C. 2.2 Zgradba inštalacijskega odklopnika Na sliki 2.1 so razvidni sestavni deli, ki določajo delovanje inštalacijskega odklopnika in s tem njegove lastnosti. Ti deli so: a elektromagnetni sprožnik: zagotavlja zahtevan prag proženja in trenutno delovanje do nazivne kratkostične zmogljivosti, z njim se definira izklopna karakteristika aparata; b bimetalni sprožnik: zagotavlja izklop v področju preobremenitev, to je od minimalnega toka delovanja do praga proženja elektromagnetnega sprožnika; c stikalni mehanizem: poskrbi, da se delovanje elektromagnetnega in bimetalnega sprožnika prenese na kontaktni sklop ter omogoča ročni vklop in izklop; d kontaktni sklop: sestavljen je iz fiksnega in gibljivega kontakta;

11 11 e obločni kanal: po njem električni oblok potuje do gasilne komore; f gasilna komora: pogasi oblok, ki nastane med kratkim stikom. Slika 2.1: Inštalacijski odklopnik

12 12 3 Selektivni odklopnik Inštalacijski odklopnik (MCB) ima kar nekaj prednosti v primerjavi s talilnimi varovalkami. Prednosti inštalacijskih odklopnikov so: - enostaven vklop, - možnost večkratnega vklopa, - večpolni izklop pri pojavu napake le v enem polu, - potreben manjši prostor za vgradnjo. Zaradi teh prednosti se je v podjetju ETI d.d. porodila zamisel, da bi tudi v priključne omarice namesto dosedanjih talilnih varovalk tipa NV vgrajevali inštalacijske odklopnike, za kar pa je potrebno razviti povsem nov selektivni odklopnik, v nadaljevanju SO. Napetostno neodvisen selektivni odklopnik že izdeluje podjetje ABB, napetostno odvisnega pa izdeluje podjetje HAGER. Oba obstoječa selektivna odklopnika sta dokaj velikih dimenzij, približno štirikratne velikosti klasičnega odklopnika. SO podjetja ETI d.d. pa bi bil velikosti 2M (dveh modulov), kar je dvakratna velikost klasičnih inštalacijskih odklopnikov. Le-ti pa zasedejo manj prostora kot pa klasične talilne varovalke tipa NV, s tem pa bi pridobili tudi nekaj prostora v priključnih omaricah. Poleg prednosti v velikosti aparata pa bi imel ta selektivni inštalacijski odklopnik tudi manjšo notranjo upornost, zaradi česar bi bile njegove lastne izgube manjše. 3.1 Razred selektivnosti inštalacijskih odklopnikov Razred selektivnosti inštalacijskega odklopnika govori o uspešnosti njegovega delovanja pri kratkih stikih. Starejši tipi odklopnikov so prekinjali kratkostični tok šele v trenutku, ko je sinusni potek toka dosegel vrednost nič. Novejši inštalacijski odklopnik pa kratkostični tok tudi omeji in ga začne prekinjati že prej, preden le-ta doseže svojo maksimalno vrednost. Tako pri pričakovanem toku kratkega stika 10 ka (efektivna vrednost) odklopnik prekine tok že pri vrednosti 5 do 5,5 ka. Zaradi te njegove lastnosti je tudi energija, ki jo prepusti, bistveno manjša. Velikost te prepuščene energije pa je osnova za razvrstitev odklopnikov v razrede selektivnosti. Najslabši je razred 1, najboljši pa razred 3. Odklopniki, ki spadajo v razred 3, zelo dobro omejujejo kratkostični tok in prepustijo zelo malo energije. Prepuščena energija se sprosti in porabi na samem inštalacijskem odklopniku ter na električnih inštalacijah in napravah, ki jih ta odklopnik ščiti. Odklopnik s slabšim razredom

13 13 selektivnosti prepušča bistveno več energije, ki veliko hitreje uničuje sam inštalacijski odklopnik, ki bo zato hitro dokončno odpovedal in ga bo potrebno zamenjati. Hkrati pa tak odklopnik tudi slabše opravlja svojo osnovno funkcijo (zaščita električnih inštalacij). Prihaja do večjega obremenjevanja varovane inštalacije (segrevanja), hitrejšega staranja inštalacij in tako tudi do bistveno povečane nevarnosti požara [3]. Kvaliteta odklopnikov ima za uporabnike bistveno vlogo, kajti če: 1. odklopnik izklopi prezgodaj, se po nepotrebnem prekine delovni proces in lahko nastane tudi posredna škoda, 2. pri trenutnem izklopu ni izklapljanja, pomeni, da je odklopnik nekaj sekund izpostavljen relativno visoki termični obremenitvi, ki pospešuje njegovo staranje in krajša življenjsko dobo, 3. je nizka kratkostična zmogljivost, je potrebno tak odklopnik hitro zamenjati z novim. Prevelika prepuščena energija obremenjuje celotno inštalacijo, ki jo odklopnik ščiti, zato prihaja do močnejšega segrevanja, hitrejšega staranja in nevarnosti požara [1]. Obstajata dva nivoja selektivnosti v zaščiti električnih sistemov: popolna selektivnost, delna selektivnost Popolna selektivnost Popolno selektivnost imamo v primeru, ko odklopnik bližje napaki izklopi, odklopnik bližje napajanju pa ostane vklopljen, s čimer zagotovimo maksimalno razpoložljivost sistema. Popolna selektivnost mora biti dosežena v obeh področjih delovanja, in sicer tako v preobremenitvenem (t > 0,1s) kot tudi v kratkostičnem področju (t < 0,1s) delovanja odklopnika. Za boljšo predstavo si poglejmo sliko 3.1, kjer je prikazana popolna selektivnost med dvema inštalacijskima odklopnikoma.

14 14 Slika 3.1: Prikaz popolne selektivnosti med dvema inštalacijskima odklopnikoma Delna selektivnost V primeru, da pogoj popolne selektivnosti ni izpolnjen do pričakovanega kratkostičnega toka, govorimo o delni selektivnosti. Poglejmo si sliko 3.2, kjer je prikazana delna selektivnost dveh inštalacijskih odklopnikov. Slika 3.2: Prikaz delne selektivnosti med dvema inštalacijskima odklopnikoma Na sliki 3.2 se lahko vidi, da je selektivnost izpolnjena le v preobremenitvenem delu obeh odklopnikov, v kratkostičnem delu pa je čas izklopa obeh odklopnikov enak, zato pogoj selektivnosti ni izpolnjen. Do tega pride, ker aparat 2 nima časovne zakasnitve trenutnega delovanja pri kratkostičnem toku.

15 Zgradba in delovanje selektivnega odklopnika Zgradba selektivnega odklopnika V moji diplomski nalogi je že bilo opisano delovanje selektivnega odklopnika, vendar je bila zgradba in princip delovanja aparata takrat malo drugačen. Ker smo z meritvami prišli do zaključka, da dvokotveni elektromagnetni sprožnik ne deluje dovolj zanesljivo, smo naredili rekonstrukcijo celotnega aparata ter seveda tudi elektromagnetnega sprožnika. Zgradba selektivnega odklopnika je zelo podobna zgradbi klasičnega inštalacijskega odklopnika. Za boljše razumevanje delovanja selektivnega odklopnika si na tem mestu poglejmo sliko 3.3, kjer je prikazana njegova zgradba, hkrati pa je opisana tudi razlaga zgradbe aparata. Začetni sestavni del selektivnega odklopnika ter vseh inštalacijskih odklopnikov sta ohišje 1 ter pokrov. Ohišje in pokrov morata biti narejena iz električno neprevodnega materiala in morata omogočati nizanje stikal v baterijo na zbiralni letvi. Slika 3.3: Selektivni odklopnik

16 16 V ohišju sta vstavljeni dve sponki, in sicer dovodna sponka 2 in odvodna sponka 3. Nato je v ohišju nameščen elektromagnetni sprožnik 4, katerega podrobnejšo sestavo si bomo ogledali v enem od naslednjih poglavij. Glavni sestavni del, ki pa ga v klasičnem inštalacijskem odklopniku ni, je predupor 5, ki skrbi za selektivni odklop aparata. Pomemben del selektivnega odklopnika je tudi sekundarni bimetal 17, ki je zaporedno vezan na predupor. Povezava med fiksnim kontaktom 6 ter priključnimi sponkami 2 je izvedena preko gibljivega kontakta 7. Na fiksni kontakt 6 je privarjena tuljava elektromagnetnega sprožnika 16, ki je potem na svojem drugem koncu privarjena na odvodno sponko 3. Za izredno hiter izklop aparata skrbi stikalni mehanizem, ki ga sestavljajo: stikalni vzvod 8, kljukica 9, stikalni locen 10 ter gumb 11. Obločni sklop selektivnega odklopnika pa sestavljajo: ploščica obločnega kanala 12, gasilna komora 13 ter obločno gasilni kanal 14. Poleg vseh teh sestavnih delov pa potrebujemo še zaskočnik 15, s katerim lahko aparat pritrdimo na DIN letev Delovanje selektivnega odklopnika V primeru normalnega obratovanja selektivnega odklopnika, to je področje do 1,2-kratnik nazivnega toka, teče tok čez dovodno sponko preko gibljivega kontakta na fiksni kontakt in naprej preko tuljave na odvodno sponko. Aparat v tem področju delovanja ne sme izklopiti, saj je tok manjši od preobremenitvenega toka. V primeru pojava preobremenitvenega toka, to je tok od 1,2-kratnika pa do 6,25-kratnika nazivnega toka, je potek oz. pot toka popolnoma ista, razlika je v tem, da bo v tem področju izklopil bimetalni sprožnik, ki je indirektno ogrevan s pomočjo tuljave elektromagnetnega sprožnika. Ko se pojavi kratkostični tok, to je tok višji od 6,25-kratnika nazivnega toka, pa začne delovati elektromagnetni sprožnik. Ko ta tok steče skozi tuljavo elektromagnetnega sprožnika, povzroči osni pomik kotve do jedra. S tem se igla pomakne in povzroči razklenitev glavnega kontakta, zaradi česar steče potem tok skozi pomožni kontakt preko bimetalnega sprožnika in

17 17 vgrajenega predupora na odvodno sponko. Ob razkleniti glavnega kontakta se pojavi električni oblok, ki se pogasi preko gasilne komore. Če je pojav kratkostičnega toka le kratkotrajen, se povečano elektromagnetno polje, ki ga ustvari tuljava elektromagnetnega sprožnika, zmanjša na prvotno vrednost in zato se glavni kontakt znova sklene. Če pa je pojav kratkostičnega toka dolgotrajnejši, se zaradi povečanega toka skozi pomožni kontakt in s tem tudi bimetalni sprožnik, bimetalni sprožnik začne upogibati in pritisne na stikalni mehanizem, zaradi česar aparat po določenem času izklopi. V tem primeru je povrnitev aparata v normalno obratovanje možno samo ročno preko gumba. 3.3 Zgradba in delovanje elektromagnetnega sprožnika Zgradba elektromagnetnega sprožnika Eden od najpomembnejših sestavnih delov tega selektivnega inštalacijskega odklopnika je elektromagnetni sprožnik, ki poleg predupora skrbi za selektivno izklapljanje v kratkostičnem področju delovanja selektivnega inštalacijskega odklopnika. Na sliki 3.4 je prikazana zgradba elektromagnetnega sprožnika. Slika 3.4: Skica sprožnika v osni simetriji Na sliki 3.4 je prikazana zgradba sprožnika. Slika je narisana v osni simetriji in je enaka, kot se bo v nadaljevanju uporabljala pri simulacijah. Glavni parameter je širina zračne reže med

18 18 malo kotvo in jedrom, ki v normalnem režimu obratovanja avtomata znaša 2,5 mm in se v primeru preobremenitve zmanjša oziroma je ni več, saj se kotva in jedro združita. Jedro je fiksni del magnetnega sprožnika, na katerem je pritrjen jarem, tako da je zračna reža med jedrom in jarmom minimalna. Se pa na drugem koncu jarma pojavi zračna reža med jarmom in kotvo, saj se med njima nahaja izolacijska cevka, ki služi za pozicioniranje in vodilo kotvi. Na aluminijasto cevko je navita tuljava. V našem primeru se na njej nahaja 8 ovojev izolirane in na koncu razmaščene žice. Sprožnik pa sestavlja še en pomemben del, in sicer je to udarna igla, ki je osno pomično nameščena v jedru in se na eni strani nalega na kotvo, na nasprotni strani pa prosto moli izven jedra. V primeru klasičnega inštalacijskega odklopnika je zgradba elektromagnetnega sprožnika popolnoma identična do sedaj opisani zgradbi. Razlika med selektivnim elektromagnetnim sprožnikom in klasičnim elektromagnetnim sprožnikom je v bakrenem obročku. Ta obroček skrbi, da selektivni inštalacijski odklopnik tudi v kratkostičnem področju obratovanja zagotovi selektivnost v primerjavi s klasičnim inštalacijskim odklopnikom Delovanje elektromagnetnega sprožnika Ko skozi tuljavo teče tok, se v jedru, kotvi in magnetnem jarmu vzpostavi elektromagnetno polje in posledično sila med kotvo in jedrom. Ko je sila med kotvo in jedrom tako velika, da preseže silo vzmeti med kotvo in jedrom, se začne kotva gibati proti jedru. Hkrati pa mora sila med kotvo in jedrom premagati tudi proti-silo vzmeti gibljivega kontakta. Zato mora biti ta sila večja od sile, s katero gibljivi kontakt pritiska na fiksni kontakt. Ena izmed poglavitnih zahtev je, da sprožnik začne delovati pri 6,25-kratniku nazivnega toka, ki je meja med preobremenitvenim in kratkostičnim tokom. Pri tem toku se mora kotva začeti premikati proti jedru, s čimer premakne udarno iglo. Ta nato premakne glavni gibljivi kontakt in s tem razklene glavni tokokrog. Ker tok v tem primeru ne more več teči po prvotni poti, steče skozi predupor in bimetal po pomožnem tokokrogu. Velikost toka je sedaj omejena z velikostjo vgrajenega predupora. Ker imamo vgrajen upor, katerega upornost znaša 600 m je maksimalen tok, ki steče skozi selektivni odklopnik, približno 400 A. Ker gibljivi kontakt primarnega tokokroga drži odprt elektromagnetni sprožnik, je potrebno zagotoviti dovolj veliko silo med kotvo in jedrom tudi v času na vsake 10 ms, ko gre sinusni tok skozi

19 19 nič, saj bi se v nasprotnem primeru pojavilo vibriranje gibljivega kontakta in s tem varjenje med gibljivim in fiksnim kontaktom. Da smo odpravili to neželeno vibriranje gibljivega kontakta, smo v jedro elektromagnetnega sprožnika vgradili bakreni obroček, ki v popolnosti odpravi to vibriranje. Po določenem času pojava kratkostičnega toka se bimetal segreje in posledično upogne, zaradi česar pritisne na stikalni mehanizem in izklopi aparat. Če se kratkostični tok v času, krajšem od 30 ms, zmanjša na normalno vrednost, ki je manjša oziroma enaka nazivni vrednosti selektivnega inštalacijskega odklopnika, se mora znova skleniti glavni kontakt, s čimer se znova vzpostavi primarni tokokrog.

20 20 4 Metode za numerični izračun magnetnega polja in sil v elektromagnetnem sprožniku 4.1 Osnovne enačbe elektromagnetnega polja Za vakuum in neferomagnetne materiale velja Ampéreov zakon (4.1) Enačba (4.1) nam pove povezavo med tokom I kot povzročiteljem magnetnega polja in vektorjem gostote magnetnega pretoka B. Diferencialna oblika enačbe (4.1) je rot B = 0 J, (4.2) kjer je J gostota konduktivnega toka. Poleg konduktivnega toka J pa lahko upoštevamo še t.i.»premikalni«tok. Ko upoštevamo še ta tok, pa dobi enačba (4.2) obliko. (4.3) Pri obravnavi magnetnih pojavov v prostoru, v katerem se nahajajo tudi feromagnetni pojavi, je bolj primerno uporabiti zapis Ampereovega zakona z jakostjo magnetnega polja H, (4.4) kjer je I J d A (4.5) tok, ki ga zaobjema zanka L, A pa je površina, ki jo popisuje ta zanka. A Za vakuum (in neferomagnetne materiale) velja zveza B = μ 0 H in torej enačba (4.3) dobi obliko. (4.6) Izraz predstavlja I. Maxwellovo diferencialno enačbo. Za kvazistatično polje lahko izpustimo v enačbi (4.6) še člen D in dobimo enačbo t. (4.7)

21 21 Za FEM analizo se uporablja deskritiziran zapis enačb 4.6 in 4.7 v diferencialni obliki. Na trgu obstajajo različni simulacijski programi, s pomočjo katerih se da simulirati razmere v elektromagnetnem polju. Jaz sem izbral simulacijski program COMSOL Multyphisics, saj je nadgradnja programa FEMLAB, s pomočjo katerega sem že opravljal simulacije za potrebe diplomske naloge. Različni moduli simulacijskega programa COMSOL omogočajo reševanje raznovrstnih problemov kot so elektromagnetika, fluidika, termodinamika, kemijske sestavine, mehanika, V našem primeru smo problem reševali z AC/DC modulom, ki omogoča poleg ostalih tudi reševanje magnetostatičnih problemov. Enačba 4.7 je zelo podobna enačbi, ki jo uporablja program COMSOL za izračun magnetostatičnega polja v 2D osno simetričnem sistemu. Z upoštevanjem izraza H = B/μ in B = rota dobimo enačbo, ki jo uporablja COMSOL, če uporabljamo simulacijo v 2D sistemu z osno simetrijo: ( ) = J cond, (4.8) kjer lahko člen, ostaneta le še člena:, ki je posledica premikanja nabojev, zanemarimo. Tako nam v enačbi, ki predstavlja inducirane tokove ter ( ), ki je enak izrazu rot H Kot lahko vidimo je končna oblika enačbe, ki jo COMSOL uporablja za izračun magnetostatičnega polja zelo podobna enačbi (4.7), saj se glasi ( ) (4.9) Tu je zaenkrat omenjena samo prva Maxwellova diferencialna enačba, ker program COMSOL za svoje izračune uporablja prvo Maxwellovo enačbo v diferencialni obliki. Seveda pa so v splošnem za opis elektromagnetnih pojavov pomembne še ostale tri Maxwellove enačbe, ki skupaj s prvo sestavljajo osnovne enačbe elektromagnetnega polja.

22 22 Tok lahko v simulacijskem programu COMSOL podajamo kot vsiljen tok, ali pa kot tok, ki je posledica priključene napetosti. 4.2 Robni pogoji Magnetne silnice (gostotnice) prehajajo iz ene snovi v drugo (iz železa v zrak in obratno) po lomnem zakonu magnetnega polja [5, stran 69]. Kjer se permeabilnost snovi skokoma spremeni in če na meji ni nikakršnega ploskovnega toka (K = 0), kar za naš primer velja, magnetne silnice izstopajo iz snovi 1 pod izstopnim kotom 1 in vstopajo v snov 2 pod vpadnim kotom 2 (slika 4.1). Povezava med njimi je naslednja: iz pretočnega zakona izvedemo e n1 x (H 1 -H 2 ) = K H 1t - H 2t = K, zaradi brezizvornosti magnetnega polja je e n1 (B 1 - B 2 ) = 0 B 1n = B 2n oziroma 1 H 1n = 2 H 2n in tg 1 /tg 2 = 1 / 2. Slika 4.1: Lomni zakon v magnetiki Pri elektromagnetnem sprožniku imamo dve različni snovi, to je feromagnetik (železo) in zrak. Vzemimo, da je snov 1 zrak, snov 2 pa železo. Permeabilnost zraka je : 1 = 0 = V s/(a m). Permeabilnost železa je približno: 2 = r 0 = = 0,0012 V s/(a m). Zaradi velikega količnika 2 / 1 izstopajo gostotnice pretežno pravokotno iz površine železa.

23 23 Ker pa ima program COMSOL za osnovno veličino vektor magnetnega potenciala, si na tem mestu poglejmo, kako bi se enačba za ploskovni tok glasila, če bi jo napisali s pomočjo A. Ob upoštevanju nekaterih osnovnih zvez med A in H pridemo do enačbe 1 A 1 A n n K. (4.10) Pomembno je, da razumemo naravo magnetnega pretoka, ki si najraje utira pot tam, kjer je r čim večji. COMSOL sam zagotovi pravilne mejne pogoje znotraj elementov strukture, medtem ko je potrebno na zunanjih robovih okolice ročno nastaviti ustrezen mejni pogoj. Na zunanji meji nastavimo pogoj, (4.11) kar predstavlja magnetno izolacijo sistema. To je ustrezen pogoj le za neskončno oddaljenost, zato je v praksi potrebno izbrati dovolj veliko okolico. V našem primeru ta pogoj ne predstavlja večjega problema, saj se magnetno polje zaključuje po jarmu in je izven same strukture magnetno polje izredno majhno, kar pa nam pokažejo tudi simulacije Maxwellova enačba v krožno valjnih koordinatah Ker smo zaradi lažjega in hitrejšega načina reševanja enačb v simulacijah uporabili model reševanja s pomočjo osne simetrije, si poglejmo osnovne enačbe krožno valjnega koordinatnega sistema. Koordinatne ploskve so: koncentrični valji, polravnine in vzporedne ravnine. Točka T v prostoru je podana s presečiščem ploskve krožnega valja r = konst., polravnine φ = konst. in ravnine z = konst.. Spremenljivke (koordinate) so: q r, q j, q z. (4.12) 1 2 3

24 24 Enotski vektorji so: e = e, e = e, e = e. (4.13) 1 r 2 3 z Enotski vektorji e 1, e 2 in e 3 v smeri koordinat q 1, q 2 in q 3 so pravokotni drug na drugega. Poleg tega vektorski produkt dveh da tretjega. Element dolžine je: dl e dl e dl e dl (4.14) Pri tem je dl i (i = 1, 2, 3) element dolžine v smeri koordinate q i. Element dolžine dl i je proporcionalen spremembi koordinate dq i. Faktor sorazmernosti med elementom dolžine in spremembo (elementom) koordinate imenujemo metrični ali Lamejev koeficient in ga označujemo s h i (dl i = h i dq i ). Tako je element dolžine: d l e h dq e h dq e h dq (4.15) Lamejev koeficient je določen z h i x y z q q q i i i. (4.16) S pomočjo formule (4.16) dobimo: h h 1, h h r, h h 1. (4.17) 1 r 2 3 Povezavo med kartezijevimi in krožnimi valjnimi koordinatami lahko takoj zapišemo s pomočjo slike 4.2 : x r cos, y r sin, z z. (4.18) z Iz te skupine enačb lahko pridemo nato do enačb 2 y r x y 2, arctg, z z. (4.19) x Spremenljivke r, φ in z lahko zavzamejo vrednosti v mejah: 0 r, 0 2, - z. (4.20)

25 25 Slika 4.2: Krožno valjni koordinatni sistem Rotor v splošnih pravokotnih koordinatah: e h e h e h rot B. (4.21) h h h q q q h B h B h B Ker imamo v našem primeru zaradi rotiranja okoli osi z opraviti s krožno valjnim koordinatnim sistemom, si poglejmo vrednost rot B v krožno valjnih koordinatah: e e r e r 1 rot B. (4.22) r r z B rb B r Potem, ko smo pogledali, kako se izračuna rot B v krožno valjnih koordinatah, si poglejmo še 1. Maxwellovo enačbo v krožno valjnih koordinatah. Kot smo videli, se 1.Maxwellova enačba lahko zapiše: ( ). (4.23) z z

26 26 Če na levi strani enačbe predvidimo, da je v = 0 ter da nimamo induciranega toka, se enačba poenostavi in jo lahko zapišemo: r J rot rot A. (4.24) Ko s pomočjo enačbe (4.23) rešimo izraz rot( krožno valjnih koordinatah, ki se glasi: r rota), dobimo 1.Maxwellovo enačbo v u r 1 r 1 2 r u J u 0 z z, (4.25) kjer je»u«odvisna nenična komponenta magnetnega potenciala (A φ ), deljenega z radialno koordinato r. To transformacijo izvedemo, da se izognemo singularnosti na simetrični osi [6, stran 3-37 ]. 4.4 Sila na kotvo elektromagnetnega sprožnika Na meji dveh različnih permeabilnosti deluje sila, ki ima smer normale na mejno ploskev in je usmerjena v prostor z manjšo permeabilnostjo. To pomeni, da je sila usmerjena iz kotve v zračno režo in iz jedra v zračno režo. Kotva in jedro se skušata približati, vendar, ker je jedro pritrjeno na jarem, se premika le kotva. Koristna ploskovna sila, ki bo premaknila kotvo, deluje le na horizontalnih ploskvah kotve. Na to silo torej vplivajo le magnetne silnice oz. magnetno polje, ki izstopa iz čelne ploskve kotve. Silnice, ki izstopajo iz stranske ploskve kotve, določajo silo, ki skuša kotvo raztegniti, ne pripomore pa nič k premiku kotve. Koristna ploskovna sila v nekoliko poenostavljeni obliki [15, stran 382] se glasi: f B1 t B1 n (4.26) B 1n je normalna, B 1t pa tangencialna komponenta gostote magnetnega pretoka v prostoru z 1 (zrak), 2 pa je permeabilnost železa. Pri obravnavi lomnega zakona smo ugotovili, da vektor magnetnega pretoka vstopa pravokotno v prostor z manjšo permeabilnostjo (zrak), če je

27 27 2 >> 1. Zato lahko običajno zanemarimo tangencialno komponento gostote magnetnega pretoka (B 1t = 0) in dobimo poenostavljeno obliko enačbe (4.26) : 11 2 f B1 n. (4.27) 2 1 Integral ploskovne sile po mejni ploskvi A med obema permeabilnostima nam da celotno silo: F fda. (4.28) A Ker imamo v našem primeru krožno valjni koordinatni sistem, je da = 2πr dr. Sila na kotvo za obravnavani elektromagnetni sprožnik je: 2 r2 F f r d d r. (4.29) 0 r1 Nas zanima, kolikšna magnetna sila deluje na kotvo, oziroma pritiska kotvo na jedro. Če imamo podano porazdelitev gostote magnetnega pretoka na meji med kotvo in zračno režo, lahko magnetno silo izračunamo s pomočjo enačbe (4.29). Sedaj se magnetna sila glasi: F f r r. (4.30) To enačbo bi uporabili, če bi računali sile s pomočjo magnetnega vezja. Poslužil se bom enostavnejšega načina in bom računal silo s pomočjo računalniškega programa COMSOL, ki omogoča izračun sile s pomočjo Maxwellovega tenzorja napetosti (ang. Maxwell Stress Tensor) (COMSOL 4.2a, Dynamic Help) ( ) ( ), (4.31) ki predstavlja površinsko silo, ki jo je potrebno še integrirati, da dobimo sumarno silo. V našem primeru zaradi simulacije strukture v osni simetriji uporabimo enačbo (4.29).

28 28 5 Numerične simulacije elektromagnetnega sprožnika selektivnega inštalacijskega odklopnika V tem poglavju bomo opravili različne tipe simulacij, ki nam bodo potrdili smiselnost uporabe bakrenega kratkostičnega obročka v elektromagnetnem sprožniku. Začeli bomo z enostavnimi simulacijami, ki nam bodo pokazale osnovni princip delovanja elektromagnetnega sprožnika, kasneje pa bomo nadgrajevali naš model do končnega najzahtevnejšega modela sprožnika. Kot najenostavnejši primer bomo prikazali rezultate simulacij, kjer bomo vzeli konstanten r, hkrati pa tudi ne bo vgrajenega bakrenega kratkostičnega obročka. Nadaljevali bomo s simulacijami, kjer bomo upoštevali magnetilno krivuljo za železo in s tem dejstvo, da je r =f(b). Sledile bodo simulacije, kjer bomo poleg magnetilne krivulje železa upoštevali tudi funkcijo bakrenega obročka. Opravile se bodo simulacije z različnim položajem bakrenega obročka, s čimer bomo prikazali, da je prispevek sile, ki je posledica induciranega polja v bakrenem obročku, zelo odvisen od položaja vgrajenega bakrenega kratkostičnega obročka. Nazadnje pa bomo prikazali še vpliv magnetnega jarma na velikost sile med kotvo in jedrom. 5.1 Numerične simulacije elektromagnetnega sprožnika brez kratkostičnega obročka in s konstantnim r Najprej bom prikazal elektromagnetno polje v poenostavljenem elektromagnetnem sprožniku, kjer bo r konstanten in kjer ne bo vgrajenega bakrenega kratkostičnega obročka. Elektromagnetno polje bo posledica izmeničnega toka skozi tuljavo. V nadaljevanju pa bom potem prikazal potek sile med kotvo in jedrom v odvisnosti od velikosti zračne reže med kotvo in jedrom. Za začetek si poglejmo nastavitve v simulacijskem programu COMSOL, ki so bile potrebne, da smo prišli do pravilnih rezultatov. Kot prvo je potrebno seveda ustrezno narisati strukturo. Ker smo lahko zaradi oblike elektromagnetnega sprožnika uporabili osno simetrijo, nam je to zelo olajšalo celotno delo, saj nam ni bilo treba risati 3D strukture. Glede na to, da so se opravile simulacije s pomočjo osne simetrije, je bilo potrebno magnetni jarem prilagoditi, saj

29 29 je to edini element elektromagnetnega sprožnika, ki ne zadosti pogojem osne simetrije. V mojem primeru sem vzel manjšo debelino jarma, s čimer sem dobil približno enak presek, kot ga ima jarem dejansko. Lahko pa bi tudi prilagodil relativno permeabilnost železa za jarem, s čimer bi prišel do enakih rezultatov. Ko je narisana struktura v osni simetriji, pa je potrebno določiti materiale za posamezne dele te strukture. Materiale lahko izberemo iz knjižnice, ki jo ima COMSOL, lahko pa tudi ročno nastavimo določene parametre, s čimer definiramo materiale posameznim sestavnim delom strukture. Ker je knjižnica simulacijskega programa zelo bogata z materiali, sem materiale izbral kar iz te knjižnice. Izbrati je bilo potrebno sledeče materiale: 4. zrak (okolica, zračne reže), 5. baker (tuljava, kratkostični obroček v nadaljevanju), 6. železo (kotva, jedro, magnetni jarem), 7. plastika (udarna igla). Ko imamo izdelano geometrijo simulacijske strukture, ki je razdeljena na posamezne elemente, se tem elementom določijo materiali. Te materiale določimo v Materials, ki je podzavihek Model-a. Glej sliko 5.1. Slika 5.1: Definiranje materialov sestavnim delom strukture

30 30 Ko je narisana struktura ter določeni materiali, ki to strukturo sestavljajo, pa je potrebno določiti parametre, ki so potrebni za zagon simulacije. Prvi parameter, ki ga je bilo potrebno nastaviti, je sinusni tok, ki smo ga nastavili kot parameter zaradi lažjega nadaljnega spreminjanja. Definirali smo ga kot I v Single-Turn Coil Domain, ki je podzavihek Magnetic Fieldsa, kar je prikazano na sliki 5.2. Slika 5.2: Definiranje toka skozi navitje Ker pa je ta tok definiran kot parameter, pa je bilo potrebno ta parameter določiti. Določili smo ga v Global Definitions in njegovem podzavihku Variables 1a, kar je prikazano na sliki 5.3. Slika 5.3: Določitev parametra toka

31 31 Ker smo za začetek opravili simulacije s konstantno relativno permeabilnostjo železa, je bilo potrebno le-to nastaviti. Nastavi se v zavihku Amperes Law1, ki se nahaja v Magnetic Fields, kar je prikazano na sliki 5.4. Slika 5.4: Določitev konstantne relativne permeabilnosti železa Kot zadnje, pa je potrebno nastaviti še gostoto mreže oz. velikost elementov, saj je potrebno najti nek kompromis med gostoto mreže in hitrostjo oz. natančnostjo izračuna. Jaz sem s poizkušanjem ter lastnimi izkušnjami določil različne gostote mreže v različnih sestavnih delih sprožnika. Kje in kako se nastavi gostota mreže (Mesh), pa je prikazano na sliki 5.5. Kot lahko vidimo na tej sliki, sem za tiste dele strukture, kjer se lahko pojavijo stresanja polja ali pa kakšne simulacijske nesingularnosti, vzel bolj gosto mrežo, s čimer sem tudi pospešil računanje. Če je mreža preredka, se lahko hitro pojavi problem nekonvergiranja rešitve ter s

32 32 tem nepravilen rezultat same simulacije. Če vzamemo zelo gosto mrežo, pa se lahko hitro pojavi problem z notranjim pomnilnikom računalnika, oz. hitrostjo izračuna. Slika 5.5: Določitev gostote mreže strukturi Ko se določijo vsi parametri, pa se lahko zažene simulacijo. Pri vseh nadaljnjih simulacijah sem na podoben način določeval parametre, kot je to opisano v tem podpoglavju. Vsi rezultati, tako tabelarični kot tudi grafični, se pregledujejo v zavihku Results Izmenični tok Najprej si na sliki 5.6 poglejmo sinusni potek vzbujalnega toka, ko je bila efektivna vrednost izmeničnega toka skozi tuljavo približno 70 A. Ker je vsiljeni sinusni tok skozi tuljavo frekvence 50 Hz, je perioda signala 20 ms. Tok doseže svojo maksimalno vrednost po 5 ms in nato na vsakih 10 ms.

33 33 Slika 5.6: Potek vzbujalnega sinusnega toka efektivne vrednosti 70 A Ker se zaradi toka skozi tuljavo v železu pojavi magnetno polje, si poglejmo potek oz. porazdelitev gostote magnetnega pretoka v elektromagnetnem sprožniku ob različnih časih. Slika 5.7: Porazdelitev gostote magnetnega pretoka po 35 ms pri AC toku

34 34 Slika 5.8: Porazdelitev gostote magnetnega pretoka po 40 ms pri AC toku Kot lahko vidimo na slikah 5.7 in 5.8 se gostota magnetnega pretoka s časom pri izmeničnem vsiljenem toku spreminja. V trenutku, ko gre tok skozi 0, to je pri frekvenci signala 50 Hz na vsakih 10 ms, je tudi gostota magnetnega pretoka skoraj enaka 0. Za boljšo predstavo si poglejmo sliko 5.9, na kateri bosta prikazana tako tok kot tudi gostota magnetnega pretoka v odvisnosti od časa. Slika 5.9: Gostota magnetnega pretoka in potek vsiljenega toka v odvisnosti od časa

35 35 Zaradi boljšega prikaza je gostota magnetnega pretoka pomnožena s faktorjem 80000, da dobimo isto vrednostno skalo. Vidimo lahko, da sta vsiljeni tok skozi tuljavo in gostota magnetnega pretoka v elektromagnetnem sprožniku popolnoma v fazi, tako je vrednost gostote magnetnega pretoka enaka 0, ko je vrednost sinusnega signala toka enaka 0. V simulacijskem programu COMSOL smo se na začetku soočili tudi s problemom induciranega toka v tuljavi, ki pa ga v dejanskem sprožniku ni, saj smo tuljavo simulirali z 10-imi posamično zaključenimi ovoji. Ta inducirani tok nam je povzročil določeno zamaknitev magnetnega polja, zato vsiljeni tok in magnetno polje nista bila sinhronizirana. Ta problem smo rešili na dokaj enostaven način, in sicer tako, da smo tuljavi določili dokaj visoko upornost, s čimer smo preprečili induciranje toka v tuljavi. Efektivno vrednost toka doseže signal po 2,5 ms ter nato na vsake 5 ms. Za primerjavo je na sliki 5.10 prikazana gostota magnetnega pretoka skozi elektromagnetni sprožnik v času 12,6 ms, to je v času, ko je vrednost izmeničnega toka enaka 70,7 A. Slika 5.10: Porazdelitev gostote magnetnega pretoka po 12,6 ms pri AC toku

36 Vpliv zračne reže med kotvo in jedrom na gostoto magnetnega pretoka Vse simulacije, ki so bile narejene v tem podpoglavju, so bile narejene pri efektivni vrednosti izmeničnega toka skozi tuljavo I ef = 70,7 A; hkrati pa je bila tudi relativna permeabilnost železa konstantna. Začetna zračna reža med jedrom in kotvo znaša 2,5 mm. S tako veliko zračno režo dosežemo dovolj velik hod udarne igle, ki mora v primeru pojava kratkostičnega toka odmakniti gibljivi kontakt od fiksnega kontakta za vsaj 5 mm. Večji je razmak med gibljivim in fiksnim kontaktom, boljše je gašenje kratkostičnega obloka med kontaktoma. Za začetek si poglejmo gostoto magnetnega pretoka v sprožniku pri zračni reži velikosti 2,5 mm med kotvo in jedrom, ki je prikazana na sliki 5.6. Slika 5.6: Porazdelitev gostote magnetnega pretoka pri zračni reži 2,5 mm

37 37 Na sliki 5.6 je lepo vidno, da prihaja do velikega stresanja magnetnega polja v zračni reži med jedrom in kotvo. Do tega pride, ker je zračna reža med njima relativno velika, zaradi česar je tudi gostota magnetnega pretoka v železu, in sicer tako v jedru kot v kotvi, dokaj majhna. Na sliki 5.7 je prikazano le polje v zračni reži med kotvo in jedrom, s čimer je bolje vidno stresanje magnetnega polja v tem področju. Slika 5.7: Podrobnejši prikaz porazdelitev gostote magnetnega pretoka pri zračni reži 2,5 mm Takoj, ko se kotva začne premikati proti jedru in se zračna reža med njima manjša, se tudi stresano polje zmanjšuje, posledično pa se veča gostota magnetnega pretoka v jedru in kotvi. Za primerjavo je na sliki 5.8 prikazana gostota magnetnega pretoka v sprožniku pri velikosti zračne reže, med jedrom in kotvo, 0,01 mm. Vidimo lahko, da v tem primeru stresanja polja skoraj ni, saj je zračna reža med kotvo in jedrom minimalna.

38 38 Slika 5.8: Porazdelitev gostote magnetnega pretoka pri zračni reži 0,01 mm Velikost sile med jedrom in kotvo v odvisnosti od velikosti zračne reže med njima V prejšnjem podpoglavju je bila prikazana odvisnost gostote magnetnega pretoka v sprožniku od velikosti zračne reže med kotvo in jedrom. Ker pa je sila kvadratično odvisna od gostote magnetnega pretoka, se zato tudi sila povečuje z zmanjšanjem zračne reže. Na grafu 5.1 je prikazan potek sile v odvisnosti od velikosti zračne reže (dx) med jedrom in kotvo.

39 sila [N] ,5 3 2,5 2 1,5 1 0,5 Electromagnetic force, z component (N) dx=2,5mm Electromagnetic force, z component (N) dx=2mm Electromagnetic force, z component (N) dx=1mm Electromagnetic force, z component (N) dx=0,5mm Electromagnetic force, z component (N) dx=0,1mm Electromagnetic force, z component (N) dx=0,01mm 0 0,0000 0,0100 0,0200 0,0300 0,0400 čas [s] Graf 5.1: Sila v odvisnosti od velikosti zračne reže med jedrom in kotvo Z zmanjševanjem zračne reže med kotvo in jedrom maksimalna sila med njima strmo narašča in pri minimalni zračni reži 0,01 mm znaša že okoli 3,9 N. Na tem mestu naj ponovno opomnim, da so bile te simulacije narejena na enostavnem modelu sprožnika, kjer je bila relativna permeabilnost železa konstantna.

40 B [T] Elektromagnetni sprožnik brez kratkostičnega obročka in z upoštevanjem relativne permeabilnosti železa V prejšnjem poglavju smo opravili simulacije na enostavnem modelu sprožnika, kjer je bila relativna permeabilnost železa konstantna, in sicer je znašala Vendar pa je potrebno za realne rezultate upoštevati tudi relativno permeabilnost železa ( r ), ki pa se z večanjem gostote magnetnega pretoka (B) manjša, se pravi, da je r = f(b). V našem primeru smo za naše simulacije vzeli kar železo, ki se nahaja v bazi programa COMSOL 4.2a. Na grafu 5.2 se nahaja magnetilna krivulja za železo oz. prikaz odvisnosti gostote magnetnega pretoka (B) od magnetne poljske jakosti (H) 2,5 2 1,5 1 0, H[A/m] Graf 5.2: Gostota magnetnega pretoka v odvisnosti od magnetne poljske jakosti Kot lahko vidimo na grafu, pri določeni gostoti magnetnega pretoka (B) železo preide v nasičenje. To se začne dogajati, ko je gostota magnetnega pretoka višja od 1,4 T, zato mu po

41 41 tej vrednosti gostote magnetnega pretoka začne relativna permeabilnost r naglo padati in v popolnem nasičenju doseže vrednost 1, kar je tudi relativna permeabilnost zraka. Za razliko od predhodnih simulacij je bilo potrebno sedaj nastaviti le še relativno permeabilnost kot funkcijo gostote magnetnega pretoka. Relativna permeabilnost, ki je odvisna od gostote magnetnega pretoka, se določi v Magnetic Fields, in sicer v zavihku Ampere's Law 2. V tem zavihku je potrebno določiti, za katere sestavne dele strukture bodo te nastavitve veljale, potem pa je potrebno izbrati pogoj H=f(B). Za boljšo predstavo si poglejmo sliko 5.9, kjer je prikazano okno, kjer se nastavijo parametri, da simulacijski program COMSOL operira z relativno permeabilnostjo železa, ki je odvisna od gostote magnetnega pretoka. Slika 5.9: Prikaz nastavitve parametrov za r =f(b) Ko imamo določen ta zelo pomemben parameter, pa se lahko lotimo simulacij, kjer bomo lahko videli, kako se spreminja magnetno polje in posledično sila med kotvo in jedrom v

42 42 odvisnosti od jakosti toka skozi tuljavo, oziroma v odvisnosti od velikosti zračne reže med kotvo in jedrom, kjer pa relativna permeabilnost železa ni več konstantna Vpliv toka na gostoto magnetnega pretoka v sprožniku S povečevanjem toka oz. amperovojev se povečuje tudi gostota magnetnega pretoka v elektromagnetnem sprožniku. Za boljšo predstavo si poglejmo rezultate simulacij, kjer je na slikah 5.10, 5.11 in 5.12 prikazana gostota magnetnega pretoka v sprožniku pri različnih izmeničnih tokovih skozi tuljavo elektromagnetnega sprožnika. V prvem primeru, glej sliko 5.10, smo opravili simulacijo, ko skozi navitje tuljave teče izmenični tok I ef =70,7 A. Slika 5.10: Porazdelitev gostote magnetnega pretoka pri zračni reži 0,01mm in I ef =70,7 A

43 43 Pri tej vrednosti toka je maksimalna gostota magnetnega pretoka, ki se pojavi na prehodu iz jedra na magnetni jarem, velikosti 1,14 T. Pri tej vrednosti gostote magnetnega pretoka železo še ni prešlo v nasičenje. Oba sestavna dela, jedro in kotva, med katerima se pojavi elektromagnetna sila, pa sta pri tej vrednosti toka še daleč od nasičenja. Na sliki 5.11 pa se že vidi, da določeni deli sprožnika prehajajo v nasičenje, predvsem jedro in magnetni jarem, ki sta volumsko manjša od kotve. Zaradi prehoda magnetnega jarma v nasičenje se del magnetnega polja sedaj zaključuje po zraku. V tem primeru je skozi navitje tuljave tekel tok I ef =283 A. Slika 5.11: Porazdelitev gostote magnetnega pretoka pri zračni reži 0,01mm in I ef =283 A Na sliki 5.12 pa je prikazana gostota magnetnega pretoka, ko skozi navitje tuljave teče tok I ef =1414 A. Na tej sliki se vidi, da je že večina sestavnih delov elektromagnetnega sprožnika, ki so narejeni iz železa, prešlo v nasičenje, zato se pojavlja dokaj veliko stresanje magnetnega polja in se zato posledično določena količina magnetnega polja zaključuje po zraku.

44 B [T] 44 Slika 5.12: Porazdelitev gostote magnetnega pretoka pri zračni reži 0,01mm in I ef =1414 A Na grafu 5.3 je prikazana razlika v gostoti magnetnega pretoka, če imamo konstantno relativno permeabilnost železa, ter če imamo relativno permeabilnost, ki je odvisna od magnetnega polja. 6,0 4,0 2,0 0,0-2,0 0 0,01 0,02 0,03 0,04 relativna permeabilnost=f(b) relativna permeabilnost=konst -4,0-6,0 t [s] Graf 5.3: Gostota magnetnega polja v odvisnosti od časa za različne oblike relatvne permeabilnosti

45 45 Presek, kjer smo računali povprečno gostoto magnetnega polja, je prikazan na sliki Slika 5.13: Presek, kjer smo računali povprečno vrednost gostote magnetnega polja Obe simulaciji sta bili narejeni pri velikosti izmeničnega pritisnjenega toka skozi tuljavo 400A. Vidimo lahko pomembno razliko, saj je v primeru, ko je relativna permeabilnost konstantna, maksimalna vrednost gostote magnetnega pretoka nekaj čez 5 T. Kot vemo, je to vrednost, ki je ne moremo doseči v nobenem materialu. Za razliko pa v primeru, ko uporabimo relativno permeabilnost, ki je funkcija gostote magnetnega polja, dosežemo maksimalno vrednost gostote magnetnega pretoka okoli 1,5 T. Pri tej vrednosti že večina magnetno prevodnih materialov preide v nasičenje in je njegova relativna permeabilnost zelo blizu vrednosti relatvne permeabilnosti zraka, ki je enaka 1.

46 F[N] Vpliv toka na velikost magnetne sile med kotvo in jedrom Z večanjem toka se veča tudi gostota magnetnega pretoka v sprožniku, vendar se ta gostota veča le do nasičenja železa. Kot smo videli v prejšnjem podpoglavju, preide železo pri dokaj velikih tokovih v nasičenje. Na grafu 5.4 je prikazan potek sile med jedrom in kotvo v odvisnosti od velikosti toka skozi tuljavo elektromagnetnega sprožnika. Kot smo že v prejšnjih podpoglavjih pokazali, je gostota magnetnega pretoka odvisna od velikosti zračne reže. Pri vseh simulacijah v tem podpoglavju je bila zračna reža med kotvo in jedrom minimalna, in sicer je znašala le 0,01 mm Electromagnetic force, z component (N) Ieff=141A Electromagnetic force, z component (N) Ieff=283A Electromagnetic force, z component (N) Ieff=495A Electromagnetic force, z component (N) Ieff=707A Electromagnetic force, z component (N) Ieff=1061A Electromagnetic force, z component (N) Ieff=1414A ,000 0,010 0,020 0,030 0,040 t[s] Graf 5.4: Velikost sile v odvisnosti od toka skozi tuljavo

47 47 Kot lahko vidimo na grafu 5.4, se pri relativno majhnih vrednostih izmeničnega toka vsako njegovo povečanje zelo pozna na velikosti sile med kotvo in jedrom. Ko enkrat železni sestavni deli elektromagnetnega sprožnika preidejo v nasičenje, pa veliko povečanje toka zelo malo doprinese k povečanju sile med kotvo in jedrom. Ker se v nasičenju gostota magnetnega pretoka skoraj ne povečuje več, se zaradi tega tudi sila, ki je kvadratično odvisna od gostote magnetnega pretoka, povečuje zelo malo. 5.3 Elektromagnetni sprožnik z bakrenim kratkostičnim obročkom in z upoštevanjem relativne permeabilnosti železa Pri pravilnem delovanju elektromagnetnega sprožnika se mora zračna reža med malo kotvo in jedrom iz začetnih 2,5 mm zmanjšati na 0 mm. Vse to naj bi se dogajalo, ko bi tok dosegel tako imenovani kratkostični tok, ki pri aparatih, v katerih je vgrajen tovrstni sprožnik, znaša 6,25-kratnik nazivnega toka. Dovolj veliko gostoto magnetnega pretoka v elektromagnetnem sprožniku zagotovimo s pravilnim številom ovojev tuljave in pravilno izbiro materiala, iz katerega so narejeni vsi elektromagnetno prevodni sestavni deli. Pri običajnem elektromagnetnem sprožniku se pojavi osciliranje elektromagnetne sile med jedrom in kotvo, s tem pa tudi osciliranje gibljivega kontakta, kar pa je neželeni pojav. Ob vsakem prehodu sinusnega vzbujalnega toka skozi ničlo sila med kotvo in jedrom pade na nič in zato se gibljivi kontakt odmakne od fiksnega kontakta. Zaradi vzbujalnega toka frekvence 50 Hz zato kontakt vibrira s frekvenco 100 Hz. To vibriranje je zelo problematično, saj se pri dokaj visokih tokovih ustvarijo idealni pogoji za obločno varjenje. Zaradi tega je bilo potrebno razmisliti, kako bi lahko odpravili vibriranje gibljivega kontakta. Za najenostavnejšo in najučinkovitejšo rešitev se je pokazala vgradnja bakrenega t.i. kratkostičnega obročka v jedro elektromagnetnega sprožnika. Bakreni obroček služi za sekundarni tokokrog, v katerem bi se, ob idealnih razmerah, induciral tok, ki bi bil v primerjavi s primarnim tokom zamaknjen za 90. Kako obroček dejansko vpliva na potek elektromagnetnega polja v elektromagnetnem sprožniku in posledično na medsebojno silo med kotvo in jedrom, pa je opisano v nadaljevanju.

48 48 Za razliko od prejšnjih simulacij je bilo potrebno sedaj le še določiti položaj bakrenega obročka in mu določiti pravilne materialne lastnosti. Vse to se je nastavilo v zavihku Models in nato podzavihku Materials. Pot za določitev lastnosti strukture je prikazana na sliki Slika 5.14: Določitev lastnosti strukture Vpliv toka na gostoto magnetnega pretoka v sprožniku Kot sem že prej zapisal, pade gostota magnetnega pretoka v železu na nič, ko gre sinusni signal toka skozi nič. V primeru uporabe bakrenega obročka, ki je vgrajen v jedro elektromagnetnega sprožnika, pa se v obročku inducira določen tok, ki pripomore, da gostota magnetnega pretoka v elektromagnetnem sprožniku nikoli ne pade na nič. Na sliki 5.15 je prikazana gostota magnetnega pretoka v elektromagnetnem sprožniku pri izmeničnem toku I ef =70,1 A. Gostota elektromagnetnega pretoka je prikazana v času 35 ms, to je čas, ko je gostota magnetnega pretoka najvišja, saj v tem trenutku sinusni signal toka doseže svojo maksimalno vrednost.

49 49 Slika 5.15: Porazdelitev gostote magnetnega pretoka pri zračni reži 0,01mm in I ef =70,7 A ter času 35 ms Pri tej vrednosti toka skozi tuljavo je maksimalna gostota magnetnega pretoka, ki se pojavi na prehodu iz jedra na magnetni jarem ter v področju jedra, kjer se nahaja kratkostični obroček, velikosti 1,2 T. Pri tej vrednosti gostote magnetnega pretoka železo še ni prešlo v nasičenje. Na sliki 5.16 je prikazana gostota magnetnega pretoka v elektromagnetnem sprožniku pri izmeničnem toku I ef =70,7 A, ter v času 40 ms. To je čas, ko je gostota magnetnega pretoka zaradi vzbujalnega toka najmanjša, saj gre v tem trenutku sinusni signal toka, ki teče skozi tuljavo, skozi nič.

50 50 Slika 5.16: Porazdelitev gostote magnetnega pretoka pri zračni reži 0,01mm in I ef =70,7 A ter času 40 ms Za primerjavo si na sliki 5.17 poglejmo gostoto magnetnega pretoka v elektromagnetnem sprožniku, ki nima vgrajenega kratkostičnega bakrenega obročka, oziroma smo za material elementa, ki predstavlja kratkostični obroček, določili železo. Na tej sliki lahko vidimo, da je v času prehoda sinusnega signala toka skozi nič gostota magnetnega pretoka v jedru in kotvi skoraj enaka nič. Za razliko pa se v istem času, 40 ms, na sliki 5.15 vidi, da se zaradi induciranega toka v obročku pojavi določeno magnetno polje v okolici bakrenega obročka.

51 51 Slika 5.17: Porazdelitev gostote magnetnega pretoka pri zračni reži 0,01mm in I ef =70,7 A ter času 40 ms; brez kratkostičnega obročka Kot lahko vidimo na sliki 5.15, se magnetno polje zaključuje okoli bakrenega kratkostičnega obročka in je v tem področju tudi najmočnejše. Iz tega lahko sklepamo, da je položaj bakrenega obročka optimalen, saj bi v primeru, ko bi ga vgradili v kakšen drug sestavni del, dobili zanemarljivo medsebojno silo med kotvo in jedrom. Sedaj pa si poglejmo še, kako velikost vzbujalnega sinusnega toka vpliva na gostoto magnetnega pretoka v času 40 ms. Na sliki 5.18 in sliki 5.19 je prikaz gostote magnetnega pretoka pri vzbujalnem toku 283 A skozi navitje tuljave. Na sliki 5.18 je prikazana gostota magnetnega pretoka v času 35 ms, na sliki 5.19 pa v času 40 ms.

52 52 Slika 5.18: Porazdelitev gostote magnetnega pretoka pri zračni reži 0,01mm in I ef =283 A ter času 35 ms Slika 5.19: Porazdelitev gostote magnetnega pretoka pri zračni reži 0,01mm in I ef =283 A ter času 40 ms

53 53 Kot lahko vidimo na sliki 5.19 je gostota magnetnega pretoka okoli bakrenega kratkostičnega obročka mnogo večja od gostote magnetnega pretoka, prikazanega na sliki Gostota magnetnega pretoka, ki je posledica induciranega toka v bakrenem kratkostičnem obročku, bi se poviševala do vrednosti, ko bi jedro oz. kotva prišla v nasičenje. Kako elementi elektromagnetnega sprožnika, ki so iz železa, prehajajo v nasičenje, se najlepše vidi preko spreminjanja relativne permeabilnosti železa. Za lažjo predstavo si najprej oglejmo razmere v elektromagnetnem sprožniku, ko skozi tuljavo»vsilimo«tok 100 A. To je prikazano na sliki Slika 5.20: Relativna permeabilnost pri toku skozi tuljavo 100 A Slika 5.20 prikazuje velikost relativne permeabilnosti v elektromagnetnem sprožniku v času 35 ms, to je čas, ko doseže sinusni signal toka skozi tuljavo maksimalno vrednost. Kot lahko vidimo, elementi sprožnika še niso prešli v nasičenje, saj je relativna permeabilnost železnih elementov v večini primerov povsod večja od 1000.

54 54 Da bi videli, kako določeni elementi elektromagnetnega sprožnika prehajajo v nasičenje, smo povečali tok skozi tuljavo na 400 A. Vrednosti relativne permeabilnosti v elektromagnetnem sprožniku pri tem toku prikazujeta sliki 5.21 in Na sliki 5.21 je prikazana vrednost relativne permeabilnosti po času 35 ms, to je času, ko je vrednost sinusnega signala toka skozi tuljavo maksimalna. Slika 5.21: Relativna permeabilnost pri toku skozi tuljavo 400A in času 35 ms Na sliki 5.21 lahko vidimo, da je večino jedra že v področju nasičenja, predvsem del, kjer je vgrajen bakreni kratkostični obroček, saj je tam presek jedra najmanjši. V nasičenju je tudi že del magnetnega jarma, vendar le tisti del, ki je bližje jedru. Za razliko od jedra pa je relativna permeabilnost v kotvi še vedno skoraj po celotnem preseku večja od 1000, kar nam nakazuje, da je še daleč od področja nasičenja.

55 55 Na sliki 5.22 pa je prikazana porazdelitev relativne permeabilnost v času 30 ms, to je v času, ko gre vrednost sinusnega toka skozi tuljavo, skozi 0, zato se v tem času pojavi le magnetno polje, ki je posledica induciranega toka v kratkostičnem bakrenem obročku. Slika 5.22: Relativna permeabilnost pri toku skozi tuljavo 400A in času 30 ms Za razliko od slike 5.21 lahko na sliki 5.22 vidimo, da zopet večina sestavnih delov elektromagnetnega sprožnika ni v nasičenju, ker se je jakost magnetnega polja zmanjšala pod mejo, ki povzroči nasičenje železa. Za boljšo predstavo, kdaj material preide v področje nasičenja, si poglejmo graf 5.2, kjer je podana magnetilna krivulja za železo. V področju, ko je vrednost gostote magnetnega pretoka med 1 in 1,5 T, to je področje, kjer se krivulja začne lomiti, prehaja material v področje nasičenja. Ko postane krivulja skoraj vzporedna z x osjo, pa je material v popolnem nasičenju, saj je v tem primeru njegova relativna permeabilnost enaka 1, kar je vrednost relativne permeabilnosti za zrak.

56 induciran tok [A] Vpliv vzbujalnega toka na inducirani tok v bakrenem obročku Ker je gostota magnetnega pretoka v času 40 ms v večji meri posledica induciranega toka v bakrenem kratkostičnem obročku, pa si na tem mestu poglejmo še odvisnost induciranega toka v obročku v odvisnosti od časa ter velikosti vzbujalnega toka skozi tuljavo. Na grafu 5.5 je prikazan časovni potek induciranega toka v bakrenem kratkostičnem obročku za različne vrednosti vzbujalnega toka skozi tuljavo. Iind (Ivzb= 70,3A) Iind (Ivzb= 40A) Iind (Ivzb= 100A) Iind (Ivzb= 150A) Iind (Ivzb= 200A) Iind (Ivzb= 250A) Iind (Ivzb= 300A) Iind (Ivzb= 400A) 200 Iind (Ivzb= 500A) Iind (Ivzb= 700A) Iind (Ivzb= 1000A) ,01 0,02 0,03 0, čas [s] Graf 5.5: Velikost induciranega toka v KS obročku pri različnih velikostih vzbujalnega toka Zračna reža med jedrom in kotvo je bila neka minimalna zračna reža in je znašala 0,01mm. Kot lahko vidimo se induciran tok v bakrenem kratkostičnem obročku povečuje z velikostjo vzbujalnega toka skozi navitje tuljave. Velja še omeniti, da bi bila v idealnih razmerah ( r =konst.) oblika induciranega toka sinusne oblike. Ker pa jedro pri večjih vrednostih

57 57 vzbujalnega toka preide v nasičenje, tudi fluks skozi obroček ni več sinusne oblike, posledično pa tudi inducirani tok. Za lepši prikaz fazne zakasnitve med vzbujalnim in induciranim tokom si poglejmo graf 5.6, na katerem sta prikazana oba toka pri efektivni vrednosti vzbujalnega toka 250 A. Graf 5.6: Potek vzbujalnega in induciranega toka pri I eff =250 A Na grafu 5.6 lahko vidimo, da fazni zamik med vzbujalnim in induciranim tokom ni idealnih 90, ampak je tam nekje okoli 120. Vidimo lahko tudi, da je induciran tok popačene sinusne oblike, saj ga požene inducirana napetost, ki se inducira v bakrenem obročku, ki pa je posledica spremembe magnetnega polja po času. Fazni kot je tako odvisen od lastne induktivnosti zanke kot tudi od upornosti zanke, torej od induktivnih in uporovnih lastnosti obročka. Poleg vzbujalnega in induciranega toka je na grafu prikazana tudi gostota magnetnega pretoka po preseku na prehodu med kotvo in jedrom. Kot lahko vidimo, vrednost gostote magnetnega

58 58 pretoka nikoli ne pade na vrednost nič, kar je posledica induciranega toka v bakrenem kratkostičnem obročku, ki se nahaja v jedru kombiniranega elektromagnetnega sprožnika. Če gledamo graf 5.5, lahko vidimo, da se fazni zamik med vzbujalnim in induciranim tokom z večanjem vzbujalnega toka manjša in znaša okoli 110 pri efektivni vrednosti vzbujalnega toka A Vpliv toka na velikost magnetne sile med kotvo in jedrom z vgrajenim kratkostičnim bakrenim obročkom v jedru sprožnika Magnetni sprožnik mora odkloniti gibljivi kontakt, ko tok doseže tok»kratkega stika«, ki znaša 6,25-kratnik nazivnega toka. Pri tem toku se mora kotva premakniti proti jedru. Da pa se kotva začne premikati, mora biti magnetna sila med kotvo in jedrom dovolj velika. Sila med jedrom in kotvo se povečuje z zmanjšanjem zračne reže med kotvo in jedrom, saj se povečuje tudi jakost magnetnega polja v jedru in kotvi, sila pa je posledično odvisna od magnetne poljske jakosti. Ko se zračna reža zmanjša na minimum, mora biti stalna sila med kotvo in jedrom večja od sile, s katero vzmet deluje na gibljivi kontakt v smeri fiksnega kontakta. Ta sila ne sme nikoli pasti pod vrednostjo sile vzmeti kontakta, saj bi v tem primeru začel gibljivi kontakt vibrirati, kar pa lahko pripelje do pojava varjenja med gibljivim in fiksnim kontaktom. Na grafu 5.4 lahko vidimo, da v primeru, ko v jedru ni vgrajenega kratkostičnega bakrenega obročka, pade sila med kotvo in jedrom na vrednost nič vsakič, ko gre sinusni signal vzbujalnega toka skozi nič, se pravi, da gibljivi kontakt vibrira s frekvenco 100 Hz. Na grafu 5.7 pa lahko vidimo, kako je sila odvisna od velikosti vzbujalnega toka v primeru, ko imamo v jedru vgrajen bakreni kratkostični obroček. V tem primeru sila med kotvo in jedrom nikoli ne pade na nič. Iz grafa lahko vidimo, da se sila med kotvo in jedrom pri dokaj nizkih efektivnih vrednostih vzbujalnega toka z vsakim povečanjem le-tega procentualno zelo poveča. Ko pa kotva in jedro preideta v nasičenje, je pa procentualno povečevanje sile zelo majhno, saj je gostota

59 sila [N] 59 magnetnega polja v obeh sestavnih delih že dosegla maksimum. Relativna permeabilnost železa je v tem primeru 1 in se železo obnaša kot zrak, zato so izgube v tem primeru zelo velike I=70A I=100A I=150A I=200A I=300A I=500A I=750A I=1000A ,005 0,01 0,015 0,02 0,025 0,03 0,035 0,04 čas [s] Graf 5.7: Potek sile pri različnih tokovih skozi tuljavo Kot je bilo že omenjeno, ne sme prihajati do vibriranja gibljivega kontakta med preobremenitvijo, saj v tem primeru lahko pride do varjenja, kar za posledico prinese nedelovanje aparata. Zaradi tega mora biti tudi minimalna sila med jedrom in kotvo večja od proti-sile vzmeti. To proti-silo vzmeti sestavljata dve vzmeti, in sicer vzmet, ki je vgrajena v elektromagnetnem sprožniku in skrbi, da v normalnem področju obratovanja, to so tokovi do nazivne obremenitve aparata, ne prihaja do neželenega vibriranja kotve in s tem neželenih šumov v aparatu. Druga vzmet pa je kontaktna vzmet, ki skrbi za dovolj veliko kontaktno silo med gibljivim in fiksnim kontaktom. Če je ta sila premajhna, lahko prihaja do prevelikega segrevanja aparata.

60 F [N] 60 Iz prakse velja»nenapisano pravilo«, da mora biti kontaktna sila med gibljivim in fiksnim kontaktom za vsak nazivni amper aparata 0,1 N, iz česar sledi, da je potrebna kontaktna sila za aparat z nazivnim tokom 32 A približno 3 N. Potek minimalne in maksimalne sile v odvisnosti od vzbujalnega toka si lahko pogledamo na grafu Fmin [N] Fmax [N] I vzb [A] Graf 5.8: Potek minimalne in maksimalne sile med jedrom in kotvo pri različnih vzbujalnih tokovih skozi tuljavo Graf 5.8 nam pokaže, da smo z določeno strukturo dobili v vsakem področju oziroma ob vsakem času dovolj veliko silo, da se ne pojavi neželeno vibriranje gibljivega kontakta.

61 Vpliv položaja vgradnje bakrenega kratkostičnega obročka na velikost sile med kotvo in jedrom Kljub temu, da smo z obstoječo geometrijo pridobili dovolj veliko silo med kotvo in jedrom tudi v času, ko gre sinusni signal vsiljenega toka skozi nič, nas zanima, kako položaj bakrenega kratkostičnega obročka vpliva na velikost sile med kotvo in jedrom. Ker je obstoječo geometrijo jedra izredno težko izdelovati v velikoserijski proizvodnji, si poglejmo, kaj se zgodi, če si za položaj bakrenega obročka izberemo dva skrajna položaja, ki sta iz vidika velikoserijske proizvodnje veliko enostavnejša za montažo Položaj obročka na zunanji strani jedra V prvem primeru si poglejmo, kaj se zgodi s silo, če vgradimo bakreni kratkostični obroček povsem na zunanjo stran jedra. Na sliki 5.23 je prikazana nova oblika strukture, na kateri so bile opravljene simulacije. Slika 5.23: Struktura z bakrenim obročkom na zunanji strani jedra Presek obročka je v tem primeru ostal isti, da smo lahko dobili primerjavo v velikosti sile med kotvo in jedrom, če le spremenimo položaj kratkostičnega bakrenega obročka. Na grafu 5.9 je prikazan potek induciranga toka skozi bakreni obroček. Zaradi lažje primerjave s predhodnimi rezultati sem znova vzel iste vrednosti vzbujalnega toka, kot so bile pri grafu 5.6.

62 sila [N] 62 Graf 5.9: Potek vzbujalnega in induciranega toka pri I eff =250 A Na grafu 5.9 lahko vidimo, da se tudi v tem primeru, ko imamo vgrajen bakreni kratkostični obroček na zunanji strani jedra, inducira tok v tem bakrenem obročku, vendar dosti manjši kot če imamo vgrajen obroček v notranjosti jedra. Kot pa lahko vidimo na naslednjem grafu, grafu 5.10, pa ta induciran tok nima nobenega vpliva na medsebojno silo med kotvo in jedrom I=70A I=100A I=150A I=200A I=300A I=500A I=750A I=1000A 0 0 0,01 0,02 0,03 0,04 čas [s] Graf 5.10: Potek sile pri različnih tokovih skozi tuljavo

63 63 Na grafu 5.10 vidimo, da kljub induciranem toku v bakrenem obročku pade sila med kotvo in jedrom na nič. Glede na prikazani graf 5.10 lahko pridemo do zaključka, da obroček na skrajni zunanji strani jedra nima pozitivnega učinka na medsebojno silo med kotvo in jedrom Položaj obročka na notranji strani jedra V tem primeru pa si poglejmo, kaj se zgodi s silo, če vgradimo bakreni kratkostični obroček povsem na notranjo stran jedra. Na sliki 5.24 je prikazana nova oblika strukture, na kateri so bile opravljene simulacije. Slika 5.24: Struktura z bakrenim obročkom na notranji strani jedra Tudi v tem primeru je ostal presek obročka isti, da smo lahko dobili primerjavo v velikosti sile med kotvo in jedrom, če le spremenimo položaj kratkostičnega bakrenega obročka. Na grafu 5.11 je prikazan potek induciranga toka skozi bakreni obroček. Zaradi lažje primerjave s predhodnimi rezultati sem znova vzel iste vrednosti vzbujalnega toka, kot so bili pri grafu 5.6 in 5.9.

64 64 Graf 5.11: Potek vzbujalnega in induciranega toka pri I eff =250 A Iz grafa 5.11 lahko razberemo, da je ta položaj bakrenega kratkostičnega obročka popolnoma nekoristen, saj je induciran tok v obročku minimalen, in sicer kar približno 100-krat manjši od induciranega toka, ko je bil bakreni obroček vgrajen v notranjosti jedra. S pomočjo gornjih grafov lahko pridemo do zaključka, da mora biti bakreni kratkostični obroček vgrajen v jedru na mestu, kjer bo iz obeh strani obdan z železom. Glede na rezultate lahko z gotovostjo trdim, da mora biti obroček vgrajen v jedru tako, da se nahaja v zunanji polovici jedra, saj se je pokazalo, da se v primeru, ko je obroček popolnoma na zunanji strani, inducira v njem mnogo večji tok, kot pa če je obroček vgrajen popolnoma na notranji strani.

65 Elektromagnetni sprožnik z bakrenim kratkostičnim obročkom, z upoštevanjem relativne permeabilnosti železa, vendar brez magnetnega jarma Ker želimo vedno narediti pocenitve, oziroma narediti čim bolj enostavno strukturo, pa si poglejmo, kaj se zgodi z magnetnimi razmerami v sprožniku, če mu odstranimo magnetni jarem. Nova struktura elektromagnetnega sprožnika je prikazana na sliki Slika 5.25: Elektromagnetni sprožnik brez magnetnega jarma Tudi v tem primeru, ko imamo strukturo brez magnetnega jarma, nas zanima, ali se v bakrenemu obročku inducira dovolj toka, da preprečimo vibriranje gibljivega kontakta. Poleg tega pa nas zanima tudi, če takšna struktura sploh zagotovi dovolj veliko medsebojno silo med kotvo in jedrom, da premaga silo kontaktiranja. Najprej si na grafu 5.12 poglejmo velikosti induciranega toka v bakrenem kratkostičnem obročku pri različnih tokovih vzbujanja. Zaradi lažje primerjave s strukturo, ki je imela vgrajen bakreni obroček, bomo naredili simulacije z istimi vrednostmi vzbujalnega toka, kot je to prikazano na grafu 5.5.

Numerical simulation aided design of the selective electromagnetic trigger

Numerical simulation aided design of the selective electromagnetic trigger Elektrotehniški vestnik 74(5): 73-78, 7 Electrotechnical Review: Ljubljana, Slovenija Načrtovanje elektromagnetnega sprožnika s pomočjo numerične simulacije Borut Drnovšek, Dejan Križaj ETI Elektroelement

More information

Reševanje problemov in algoritmi

Reševanje problemov in algoritmi Reševanje problemov in algoritmi Vhod Algoritem Izhod Kaj bomo spoznali Zgodovina algoritmov. Primeri algoritmov. Algoritmi in programi. Kaj je algoritem? Algoritem je postopek, kako korak za korakom rešimo

More information

Operation of an Electromagnetic Trigger with a Short-circuit Ring

Operation of an Electromagnetic Trigger with a Short-circuit Ring Operation of an Electromagnetic Trigger with a Short-circuit Ring Dejan Križaj 1*, Zumret Topčagić 1, and Borut Drnovšek 1,2 1 Faculty of Electrical Engineering, University of Ljubljana, Ljubljana, Slovenia,

More information

TOPLJENEC ASOCIIRA LE V VODNI FAZI

TOPLJENEC ASOCIIRA LE V VODNI FAZI TOPLJENEC ASOCIIRA LE V VODNI FAZI V primeru asociacij molekul topljenca v vodni ali organski fazi eksperimentalno določeni navidezni porazdelitveni koeficient (P n ) v odvisnosti od koncentracije ni konstanten.

More information

ZASNOVA IN RAZVOJ DUŠILKE ZA ENERGETSKI TRANSFORMATOR

ZASNOVA IN RAZVOJ DUŠILKE ZA ENERGETSKI TRANSFORMATOR Univerza v Ljubljani Fakulteta za elektrotehniko Borut Prašnikar ZASNOVA IN RAZVOJ DUŠILKE ZA ENERGETSKI TRANSFORMATOR Magistrsko delo Mentor: prof. dr. Danjel Vončina, univ. dipl. inž. el. Ljubljana,

More information

ENAČBA STANJA VODE IN VODNE PARE

ENAČBA STANJA VODE IN VODNE PARE ENAČBA STANJA VODE IN VODNE PARE SEMINARSKA NALOGA PRI PREDMETU JEDRSKA TEHNIKA IN ENERGETIKA TAMARA STOJANOV MENTOR: IZRED. PROF. DR. IZTOK TISELJ NOVEMBER 2011 Enačba stanja idealni plin: pv = RT p tlak,

More information

Vrste kratkih stikov

Vrste kratkih stikov Vrste kratkih stikov Seminarska naloga pri predmetu Avtor:, dipl. inž. el. (UN) Mentor: prof. dr. Grega Bizjak, univ. dipl. inž. el. Ljubljana, študijsko leto 2016/2017 Kazalo: 1 Uvod... 3 2 Nastanek kratkega

More information

OA07 ANNEX 4: SCOPE OF ACCREDITATION IN CALIBRATION

OA07 ANNEX 4: SCOPE OF ACCREDITATION IN CALIBRATION OA07 ANNEX 4: SCOPE OF ACCREDITATION IN CALIBRATION Table of contents 1 TECHNICAL FIELDS... 2 2 PRESENTING THE SCOPE OF A CALIBRATION LABOORATORY... 2 3 CONSIDERING CHANGES TO SCOPES... 6 4 CHANGES WITH

More information

Izmenični signali moč (17)

Izmenični signali moč (17) Izenicni_signali_MOC(17c).doc 1/7 8.5.007 Izenični signali oč (17) Zania nas potek trenutne oči v linearne dvopolne (dve zunanji sponki) vezju, kjer je napetost na zunanjih sponkah enaka u = U sin( ωt),

More information

Državni izpitni center. Izpitna pola 2. Ponedeljek, 28. avgust 2017 / 90 minut

Državni izpitni center. Izpitna pola 2. Ponedeljek, 28. avgust 2017 / 90 minut Š i f r a k a n d i d a t a : Državni izpitni center *M177711* JESENSKI IZPITNI ROK Izpitna pola Ponedeljek, 8. avgust 017 / 90 minut Dovoljeno gradivo in pripomočki: Kandidat prinese nalivno pero ali

More information

DOKTORSKA DISERTACIJA

DOKTORSKA DISERTACIJA UNIVERZA V LJUBLJANI NARAVOSLOVNOTEHNIŠKA FAKULTETA DOKTORSKA DISERTACIJA GAŠPER NOVAK LJUBLJANA, 2015 UNIVERZA V LJUBLJANI NARAVOSLOVNOTEHNIŠKA FAKULTETA ODDELEK ZA MATERIALE IN METALURGIJO Načrtovanje

More information

ENERGY AND MASS SPECTROSCOPY OF IONS AND NEUTRALS IN COLD PLASMA

ENERGY AND MASS SPECTROSCOPY OF IONS AND NEUTRALS IN COLD PLASMA UDK621.3:(53+54+621 +66), ISSN0352-9045 Informaclje MIDEM 3~(~UU8)4, Ljubljana ENERGY AND MASS SPECTROSCOPY OF IONS AND NEUTRALS IN COLD PLASMA Marijan Macek 1,2* Miha Cekada 2 1 University of Ljubljana,

More information

Vpliv navitja na prostorske harmonske komponente enofaznega motorja z obratovalnim kondenzatorjem

Vpliv navitja na prostorske harmonske komponente enofaznega motorja z obratovalnim kondenzatorjem Elektrotehniški vestnik 69(3-4): 175 180, 00 Electrotechnical Review, Ljubljana, Slovenija Vpliv navitja na prostorske harmonske komponente enofaznega motorja z obratovalnim kondenzatorjem Ivan Zagradišnik,

More information

Simulacija dinamičnih sistemov s pomočjo osnovnih funkcij orodij MATLAB in Simulink

Simulacija dinamičnih sistemov s pomočjo osnovnih funkcij orodij MATLAB in Simulink Laboratorijske vaje Računalniška simulacija 2012/13 1. laboratorijska vaja Simulacija dinamičnih sistemov s pomočjo osnovnih funkcij orodij MATLAB in Simulink Pri tej laboratorijski vaji boste spoznali

More information

OPTIMIRANJE IZDELOVALNIH PROCESOV

OPTIMIRANJE IZDELOVALNIH PROCESOV OPTIMIRANJE IZDELOVALNIH PROCESOV asist. Damir GRGURAŠ, mag. inž. str izr. prof. dr. Davorin KRAMAR damir.grguras@fs.uni-lj.si Namen vaje: Ugotoviti/določiti optimalne parametre pri struženju za dosego

More information

Numerično modeliranje linearnih električnih motorjev z metodo končnih elementov

Numerično modeliranje linearnih električnih motorjev z metodo končnih elementov Univerza v Ljubljani Fakulteta za elektrotehniko Jaka Peternelj Numerično modeliranje linearnih električnih motorjev z metodo končnih elementov Magistrsko delo Mentor: doc. dr. Selma Čorović Ljubljana,

More information

GEOMETRIJSKE FAZE V KVANTNI MEHANIKI

GEOMETRIJSKE FAZE V KVANTNI MEHANIKI GEOMETRIJSKE FAZE V KVANTNI MEHANIKI LARA ULČAKAR Fakulteta za matematiko in fiziko Univerza v Ljubljani V članku so predstavljene geometrijske faze, ki nastopijo pri obravnavi kvantnih sistemov. Na začetku

More information

56 1 Upogib z osno silo

56 1 Upogib z osno silo 56 1 Upogib z osno silo PREGLEDNICA 1.5 (nadaljevanje): Upogibnice in notranje sile za nekatere nosilce d) Upogibnica prostoležečega nosilca obteženega s silo F Pomik in zasuk v polju 1: w 1 = F b x (L

More information

UNIVERZA V LJUBLJANI PEDAGOŠKA FAKULTETA POLONA ŠENKINC REŠEVANJE LINEARNIH DIFERENCIALNIH ENAČB DRUGEGA REDA S POMOČJO POTENČNIH VRST DIPLOMSKO DELO

UNIVERZA V LJUBLJANI PEDAGOŠKA FAKULTETA POLONA ŠENKINC REŠEVANJE LINEARNIH DIFERENCIALNIH ENAČB DRUGEGA REDA S POMOČJO POTENČNIH VRST DIPLOMSKO DELO UNIVERZA V LJUBLJANI PEDAGOŠKA FAKULTETA POLONA ŠENKINC REŠEVANJE LINEARNIH DIFERENCIALNIH ENAČB DRUGEGA REDA S POMOČJO POTENČNIH VRST DIPLOMSKO DELO LJUBLJANA, 2016 UNIVERZA V LJUBLJANI PEDAGOŠKA FAKULTETA

More information

Calculation of stress-strain dependence from tensile tests at high temperatures using final shapes of specimen s contours

Calculation of stress-strain dependence from tensile tests at high temperatures using final shapes of specimen s contours RMZ Materials and Geoenvironment, Vol. 59, No. 4, pp. 331 346, 2012 331 Calculation of stress-strain dependence from tensile tests at high temperatures using final shapes of specimen s contours Določitev

More information

Geometrijske faze v kvantni mehaniki

Geometrijske faze v kvantni mehaniki Seminar 1-1. letnik, 2. stopnja Geometrijske faze v kvantni mehaniki Avtor: Lara Ulčakar Mentor: prof. dr. Anton Ramšak Ljubljana, november 2014 Povzetek V seminarju so predstavljene geometrijske faze,

More information

Izračun magnetnega polja okrog dolgih prevodnih nemagnetnih zaslonov

Izračun magnetnega polja okrog dolgih prevodnih nemagnetnih zaslonov Elektrotehniški vestnik 761-2): 31 37, 29 Electrotechnical Review, Ljubljana, Slovenija Izračun magnetnega polja okrog dolgih prevodnih nemagnetnih zaslonov Edi Bulić Univerza v Ljubljani, Fakulteta za

More information

SIMETRIČNE KOMPONENTE

SIMETRIČNE KOMPONENTE Univerza v Ljubljani Fakulteta za elektrotehniko SIMETRIČNE KOMPONENTE Seminarska naloga pri predmetu Razdelilna in industrijska omrežja Poročilo izdelala: ELIZABETA STOJCHEVA Mentor: prof. dr. Grega Bizjak,

More information

Izgube moči sinhronskega reluktančnega motorja

Izgube moči sinhronskega reluktančnega motorja Elektrotehniški vestnik 70(5): 267 272, 2003 Electrotechnical Review, Ljubljana, Slovenija Izgube moči sinhronskega reluktančnega motorja Damijan Miljavec 1, Miralem Hadžiselimovič 2, Konrad Lenasi 1,

More information

ANALIZA MAGNETNIH SKLEPOV TRIFAZNEGA PREKLOPNO RELUKTANČNEGA MOTORJA

ANALIZA MAGNETNIH SKLEPOV TRIFAZNEGA PREKLOPNO RELUKTANČNEGA MOTORJA UNIVERZA V MARIBORU FAKULTETA ZA ELEKTROTEHNIKO, RAČUNALNIŠTVO IN INFORMATIKO Uroš Škoflek ANALIZA MAGNETNIH SKLEPOV TRIFAZNEGA PREKLOPNO RELUKTANČNEGA MOTORJA Diplomska naloga Maribor, december 2007 2

More information

Evolucija dinamike Zemljine precesije

Evolucija dinamike Zemljine precesije Univerza v Ljubljani Fakulteta za matematiko in fiziko oddelek za fiziko Evolucija dinamike Zemljine precesije Avtor: Ivo Krajnik Ljubljana, 15. marec 2011 Povzetek Bistvo tega seminarja je v sklopu klasične

More information

Iskanje najcenejše poti v grafih preko polkolobarjev

Iskanje najcenejše poti v grafih preko polkolobarjev Univerza v Ljubljani Fakulteta za računalništvo in informatiko Veronika Horvat Iskanje najcenejše poti v grafih preko polkolobarjev DIPLOMSKO DELO VISOKOŠOLSKI STROKOVNI ŠTUDIJSKI PROGRAM PRVE STOPNJE

More information

Multipla korelacija in regresija. Multipla regresija, multipla korelacija, statistično zaključevanje o multiplem R

Multipla korelacija in regresija. Multipla regresija, multipla korelacija, statistično zaključevanje o multiplem R Multipla koelacia in egesia Multipla egesia, multipla koelacia, statistično zaklučevane o multiplem Multipla egesia osnovni model in ačunane paametov Z multiplo egesio napoveduemo vednost kiteia (odvisne

More information

Optimizacija indukcijskega segrevanja z numeričnim modeliranjem in genetskim algoritmom

Optimizacija indukcijskega segrevanja z numeričnim modeliranjem in genetskim algoritmom Elektrotehniški vestnik 76(1-2): 63-68, 2009 Electrotechnical Review: Ljubljana, Slovenija Optimizacija indukcijskega segrevanja z numeričnim modeliranjem in genetskim algoritmom Matej Kranjc, Anže Županič,

More information

Attempt to prepare seasonal weather outlook for Slovenia

Attempt to prepare seasonal weather outlook for Slovenia Attempt to prepare seasonal weather outlook for Slovenia Main available sources (ECMWF, EUROSIP, IRI, CPC.NCEP.NOAA,..) Two parameters (T and RR anomally) Textual information ( Met Office like ) Issued

More information

Državni izpitni center. Izpitna pola 1. Četrtek, 4. junij 2015 / 90 minut

Državni izpitni center. Izpitna pola 1. Četrtek, 4. junij 2015 / 90 minut Š i f r a k a n d i d a t a : Državni izpitni center *M15177111* SPOMLADANSKI IZPITNI ROK Izpitna pola 1 Četrtek, 4. junij 015 / 90 minut Dovoljeno gradivo in pripomočki: Kandidat prinese nalivno pero

More information

Simulation of multilayer coating growth in an industrial magnetron sputtering system

Simulation of multilayer coating growth in an industrial magnetron sputtering system RMZ Materials and Geoenvironment, Vol. 57, No. 3, pp. 317 330, 2010 317 Simulation of multilayer coating growth in an industrial magnetron sputtering system Simulacija rasti večplastnih prevlek v industrijski

More information

Univerza v Ljubljani Fakulteta za matematiko in fiziko Oddelek za fiziko. Fizika RFID. Seminar iz uporabne fizike

Univerza v Ljubljani Fakulteta za matematiko in fiziko Oddelek za fiziko. Fizika RFID. Seminar iz uporabne fizike Univerza v Ljubljani Fakulteta za matematiko in fiziko Oddelek za fiziko Fizika RFID Seminar iz uporabne fizike Marko Mravlak Mentor: doc. dr. Primož Ziherl 28. maj 2008 Povzetek V seminarju bomo predstavili

More information

Determining the Leakage Flow through Water Turbines and Inlet- Water Gate in the Doblar 2 Hydro Power Plant

Determining the Leakage Flow through Water Turbines and Inlet- Water Gate in the Doblar 2 Hydro Power Plant Elektrotehniški vestnik 77(4): 39-44, 010 Electrotechnical Review: Ljubljana, Slovenija Določanje puščanja vodnih turbin in predturbinskih zapornic v hidroelektrarni Doblar Miha Leban 1, Rajko Volk 1,

More information

Seminar - 1. letnik bolonjske magistrske stopnje. O energijskih bilanci v fuzijskem reaktorju - Lawsonov kriterij. Avtor: Matic Kunšek

Seminar - 1. letnik bolonjske magistrske stopnje. O energijskih bilanci v fuzijskem reaktorju - Lawsonov kriterij. Avtor: Matic Kunšek Seminar - 1. letnik bolonjske magistrske stopnje O energijskih bilanci v fuzijskem reaktorju - Lawsonov kriterij Avtor: Matic Kunšek Mentor: dr. Tomaž Gyergyek Ljubljana, marec 2014 Povzetek: V tem seminarju

More information

1) V diagramu sta prikazana plazemska koncentracijska profila po večkratnem intravenskem odmerjanju učinkovine v dveh različnih primerih (1 in 2).

1) V diagramu sta prikazana plazemska koncentracijska profila po večkratnem intravenskem odmerjanju učinkovine v dveh različnih primerih (1 in 2). NALOGE ) V diagramu sta prikazana plazemska koncentracijska profila po večkratnem intravenskem odmerjanju učinkovine v dveh različnih primerih ( in ). 0.8 0.6 0.4 0. 0.0 0.08 0.06 0.04 0.0 0.00 0 0 0 30

More information

HIGGSOV MEHANIZEM MITJA FRIDMAN. Fakulteta za matematiko in fiziko Univerza v Ljubljani

HIGGSOV MEHANIZEM MITJA FRIDMAN. Fakulteta za matematiko in fiziko Univerza v Ljubljani HIGGSOV MEHANIZEM MITJA FRIDMAN Fakulteta za matematiko in fiziko Univerza v Ljubljani V članku je predstavljen Higgsov mehanizem, ki opisuje generiranje mase osnovnih delcev. Vpeljan je Lagrangeov formalizem,

More information

USING SIMULATED SPECTRA TO TEST THE EFFICIENCY OF SPECTRAL PROCESSING SOFTWARE IN REDUCING THE NOISE IN AUGER ELECTRON SPECTRA

USING SIMULATED SPECTRA TO TEST THE EFFICIENCY OF SPECTRAL PROCESSING SOFTWARE IN REDUCING THE NOISE IN AUGER ELECTRON SPECTRA UDK 543.428.2:544.171.7 ISSN 1580-2949 Original scientific article/izvirni znanstveni ~lanek MTAEC9, 49(3)435(2015) B. PONIKU et al.: USING SIMULATED SPECTRA TO TEST THE EFFICIENCY... USING SIMULATED SPECTRA

More information

11 Osnove elektrokardiografije

11 Osnove elektrokardiografije 11 Osnove elektrokardiografije Spoznali bomo lastnosti električnega dipola in se seznanili z opisom srca kot električnega dipola. Opisali bomo, kakšno električno polje ta ustvarja v telesu, kako ga merimo,

More information

Merjenje difuzije z magnetno resonanco. Avtor: Jasna Urbanija Mentor: doc.dr.igor Serša

Merjenje difuzije z magnetno resonanco. Avtor: Jasna Urbanija Mentor: doc.dr.igor Serša Merjenje difuzije z magnetno resonanco Avtor: Jasna Urbanija Mentor: doc.dr.igor Serša Univerza v Ljubljani Fakulteta za matematiko in fiziko Februar 2005 1 Povzetek Pojav jedrske magnetne resonance omogoča

More information

Državni izpitni center ELEKTROTEHNIKA. Izpitna pola 1. Četrtek, 1. junij 2017 / 90 minut

Državni izpitni center ELEKTROTEHNIKA. Izpitna pola 1. Četrtek, 1. junij 2017 / 90 minut Š i f r a k a n d i d a t a : Državni izpitni center *M17177111* SPOMLADANSKI IZPITNI ROK ELEKTROTEHNIKA Izpitna pola 1 Četrtek, 1. junij 017 / 90 minut Dovoljeno gradivo in pripomočki: Kandidat prinese

More information

Călugăreanu-White-Fullerjev teorem in topologija DNA

Călugăreanu-White-Fullerjev teorem in topologija DNA Univerza v Ljubljani Fakulteta za matematiko in fiziko Oddelek za fiziko Călugăreanu-White-Fullerjev teorem in topologija DNA Seminar Jure Aplinc, dipl. fiz. (UN) Mentor: prof. dr. Rudolf Podgornik 26.

More information

FOTONSKI POGON. Avtor: Črt Harej Mentor: prof. dr. Simon Širca. Ljubljana, Maj 2016

FOTONSKI POGON. Avtor: Črt Harej Mentor: prof. dr. Simon Širca. Ljubljana, Maj 2016 FOTONSKI POGON Seminar I b - 1. letnik, II. stopnja Avtor: Črt Harej Mentor: prof. dr. Simon Širca Ljubljana, Maj 2016 Povzetek Človeštvo že skoraj 60 let raziskuje in uresničuje vesoljske polete. V tem

More information

Termoelektrični pojav

Termoelektrični pojav Oddelek za fiziko Seminar 4. letnik Termoelektrični pojav Avtor: Marko Fajs Mentor: prof. dr. Janez Dolinšek Ljubljana, marec 2012 Povzetek Seminar govori o termoelektričnih pojavih. Koncentriran je predvsem

More information

MICROWAVE PLASMAS AT ATMOSPHERIC PRESSURE: NEW THEORETICAL DEVELOPMENTS AND APPLICATIONS IN SURFACE SCIENCE

MICROWAVE PLASMAS AT ATMOSPHERIC PRESSURE: NEW THEORETICAL DEVELOPMENTS AND APPLICATIONS IN SURFACE SCIENCE UDK621.3:(53+54+621 +66), ISSN0352-9045 Informacije MIDEM 38(2008)4, Ljubljana MICROWAVE PLASMAS AT ATMOSPHERIC PRESSURE: NEW THEORETICAL DEVELOPMENTS AND APPLICATIONS IN SURFACE SCIENCE T. 8elmonte*,

More information

Elektromagnetno polje: 1. vaje (4. in )

Elektromagnetno polje: 1. vaje (4. in ) 1 Elektromagnetno polje: 1. vaje (4. in 5. 10. 2016) asistent: Martin Klanjšek (01 477 3866, martin.klanjsek@ijs.si) 0. Uvertura in ponovitev nekaj matematičnih pripomočkov [vektorska analiza (Gaussov

More information

Električne lastnosti vodov. Ohmske upornosti. Induktivnost vodov. Kapacitivnost vodov. Odvodnost vodov. Vod v svetlobi telegrafske enačbe.

Električne lastnosti vodov. Ohmske upornosti. Induktivnost vodov. Kapacitivnost vodov. Odvodnost vodov. Vod v svetlobi telegrafske enačbe. Električne lastnosti vodov Ohmske upornosti. Induktivnost vodov. Kapacitivnost vodov. Odvodnost vodov. Vod v svetlobi telegrafske enačbe. Primarne konstante vodov Če opazujemo električni vod iz istega

More information

UNIVERZA NA PRIMORSKEM FAKULTETA ZA MATEMATIKO, NARAVOSLOVJE IN INFORMACIJSKE TEHNOLOGIJE. O neeksaknotsti eksaktnega binomskega intervala zaupanja

UNIVERZA NA PRIMORSKEM FAKULTETA ZA MATEMATIKO, NARAVOSLOVJE IN INFORMACIJSKE TEHNOLOGIJE. O neeksaknotsti eksaktnega binomskega intervala zaupanja UNIVERZA NA PRIMORSKEM FAKULTETA ZA MATEMATIKO, NARAVOSLOVJE IN INFORMACIJSKE TEHNOLOGIJE Zaključna naloga (Final project paper) O neeksaknotsti eksaktnega binomskega intervala zaupanja (On the inexactness

More information

UNIVERZA NA PRIMORSKEM FAKULTETA ZA MATEMATIKO, NARAVOSLOVJE IN INFORMACIJSKE TEHNOLOGIJE

UNIVERZA NA PRIMORSKEM FAKULTETA ZA MATEMATIKO, NARAVOSLOVJE IN INFORMACIJSKE TEHNOLOGIJE UNIVERZA NA PRIMORSKEM FAKULTETA ZA MATEMATIKO, NARAVOSLOVJE IN INFORMACIJSKE TEHNOLOGIJE Zaključna naloga Uporaba logistične regresije za napovedovanje razreda, ko je število enot v preučevanih razredih

More information

UNIVERZA V LJUBLJANI PEDAGOŠKA FAKULTETA SAŠO ZUPANEC MAX-PLUS ALGEBRA DIPLOMSKO DELO

UNIVERZA V LJUBLJANI PEDAGOŠKA FAKULTETA SAŠO ZUPANEC MAX-PLUS ALGEBRA DIPLOMSKO DELO UNIVERZA V LJUBLJANI PEDAGOŠKA FAKULTETA SAŠO ZUPANEC MAX-PLUS ALGEBRA DIPLOMSKO DELO Ljubljana, 2013 UNIVERZA V LJUBLJANI PEDAGOŠKA FAKULTETA ODDELEK ZA MATEMATIKO IN RAČUNALNIŠTVO SAŠO ZUPANEC Mentor:

More information

Primer optimizacije obratovanja distribucijskega omrežja s prevezavami, obratovanjem v zanki in generiranjem jalove moči

Primer optimizacije obratovanja distribucijskega omrežja s prevezavami, obratovanjem v zanki in generiranjem jalove moči 26. MEDNARODNO POSVETOVANJE»KOMUNALNA ENERGETIKA 2017«J. Pihler Primer optimizacije obratovanja distribucijskega omrežja s prevezavami, obratovanjem v zanki in generiranjem jalove moči MATEJ PINTARIČ,

More information

UČNI NAČRT PREDMETA / COURSE SYLLABUS (leto / year 2017/18) Parcialne diferencialne enačbe Partial differential equations. Študijska smer Study field

UČNI NAČRT PREDMETA / COURSE SYLLABUS (leto / year 2017/18) Parcialne diferencialne enačbe Partial differential equations. Študijska smer Study field Predmet: Course title: UČNI NAČRT PREDMETA / COURSE SYLLABUS (leto / year 2017/18) Parcialne diferencialne enačbe Partial differential equations Študijski program in stopnja Study programme and level Magistrski

More information

1 Luna kot uniformni disk

1 Luna kot uniformni disk 1 Luna kot uniformni disk Temperatura lune se spreminja po površini diska v širokem razponu, ampak lahko luno prikažemo kot uniformni disk z povprečno temperaturo osvetlitve (brightness temperature) izraženo

More information

Cveto Trampuž PRIMERJAVA ANALIZE VEČRAZSEŽNIH TABEL Z RAZLIČNIMI MODELI REGRESIJSKE ANALIZE DIHOTOMNIH SPREMENLJIVK

Cveto Trampuž PRIMERJAVA ANALIZE VEČRAZSEŽNIH TABEL Z RAZLIČNIMI MODELI REGRESIJSKE ANALIZE DIHOTOMNIH SPREMENLJIVK Cveto Trampuž PRIMERJAVA ANALIZE VEČRAZSEŽNIH TABEL Z RAZLIČNIMI MODELI REGRESIJSKE ANALIZE DIHOTOMNIH SPREMENLJIVK POVZETEK. Namen tega dela je prikazati osnove razlik, ki lahko nastanejo pri interpretaciji

More information

ZDRAVLJENJE BOLNICE S VON WILLEBRANDOVO BOLEZNIJO TIPA 3 IN INHIBITORJI

ZDRAVLJENJE BOLNICE S VON WILLEBRANDOVO BOLEZNIJO TIPA 3 IN INHIBITORJI ZDRAVLJENJE BOLNICE S VON WILLEBRANDOVO BOLEZNIJO TIPA 3 IN INHIBITORJI B. Faganel Kotnik, L. Kitanovski, J. Jazbec, K. Strandberg, M. Debeljak, Bakija, M. Benedik Dolničar A. Trampuš Laško, 9. april 2016

More information

Meritve Casimirjevega efekta z nanomembranami

Meritve Casimirjevega efekta z nanomembranami Oddelek za fiziko Seminar a -. letnik, II. stopnja Meritve Casimirjevega efekta z nanomembranami avtor: Žiga Kos mentor: prof. dr. Rudolf Podgornik Ljubljana, 29. januar 203 Povzetek V tem seminarju bo

More information

Univerza v Ljubljani Fakulteta za matematiko in fiziko Oddelek za fiziko. Projekt ITER SEMINAR. Avtor: Jure Maglica Mentor: doc. dr.

Univerza v Ljubljani Fakulteta za matematiko in fiziko Oddelek za fiziko. Projekt ITER SEMINAR. Avtor: Jure Maglica Mentor: doc. dr. Univerza v Ljubljani Fakulteta za matematiko in fiziko Oddelek za fiziko Projekt ITER SEMINAR Avtor: Jure Maglica Mentor: doc. dr. Milan Čerček Ljubljana, April 005 POVZETEK V seminarju je opisan kratek

More information

JEDRSKA URA JAN JURKOVIČ. Fakulteta za matematiko in fiziko Univerza v Ljubljani

JEDRSKA URA JAN JURKOVIČ. Fakulteta za matematiko in fiziko Univerza v Ljubljani JEDRSKA URA JAN JURKOVIČ Fakulteta za matematiko in fiziko Univerza v Ljubljani Natančnost časa postaja vse bolj uporabna in pomembna, zato se rojevajo novi načini merjenja časa. Do danes najbolj natančnih

More information

Acta Chim. Slov. 2003, 50,

Acta Chim. Slov. 2003, 50, 771 IMPACT OF STRUCTURED PACKING ON BUBBE COUMN MASS TRANSFER CHARACTERISTICS EVAUATION. Part 3. Sensitivity of ADM Volumetric Mass Transfer Coefficient evaluation Ana akota Faculty of Chemistry and Chemical

More information

Državni izpitni center. Izpitna pola 1. Sobota, 27. avgust 2016 / 90 minut

Državni izpitni center. Izpitna pola 1. Sobota, 27. avgust 2016 / 90 minut Š i f r a k a n d i d a t a : Državni izpitni center *M1677111* JESENSKI IZPITNI OK Izpitna pola 1 Sobota, 7. avgust 016 / 90 minut Dovoljeno gradivo in pripomočki: Kandidat prinese nalivno pero ali kemični

More information

Zakasneli nevtroni v reaktorjih s tekočim gorivom

Zakasneli nevtroni v reaktorjih s tekočim gorivom Seminar Zakasneli nevtroni v reaktorjih s tekočim gorivom Avtor: Janez Kokalj januar, 2015 Mentor: Dr. Luka Snoj Povzetek Četrta generacija jedrskih reaktorjev, kamor spadajo tudi reaktorji na staljeno

More information

Verifikacija napovedi padavin

Verifikacija napovedi padavin Oddelek za Meteorologijo Seminar: 4. letnik - univerzitetni program Verifikacija napovedi padavin Avtor: Matic Šavli Mentor: doc. dr. Nedjeljka Žagar 26. februar 2012 Povzetek Pojem verifikacije je v meteorologiji

More information

NIKJER-NIČELNI PRETOKI

NIKJER-NIČELNI PRETOKI UNIVERZA V LJUBLJANI PEDAGOŠKA FAKULTETA ALJA ŠUBIC NIKJER-NIČELNI PRETOKI DIPLOMSKO DELO LJUBLJANA, 2016 UNIVERZA V LJUBLJANI PEDAGOŠKA FAKULTETA Dvopredmetni učitelj: matematika - računalništvo ALJA

More information

Modeli dinamičnega vzgona letalskih kril. Drugi del.

Modeli dinamičnega vzgona letalskih kril. Drugi del. Modeli dinamičnega vzgona letalskih kril. Drugi del. Sašo Knez in Rudolf Podgornik Oddelek za fiziko, Fakulteta za Matematiko in Fiziko Univerza v Ljubljani Povzetek V drugem delu tega članka se bova posvetila

More information

Testiranje programov za račun vodnega udara in uporaba na realnem primeru derivacijske hidroelektrarne

Testiranje programov za račun vodnega udara in uporaba na realnem primeru derivacijske hidroelektrarne Univerza v Ljubljani Fakulteta za gradbeništvo in geodezijo Jamova 2 1000 Ljubljana, Slovenija telefon (01) 47 68 500 faks (01) 42 50 681 fgg@fgg.uni-lj.si Univerzitetni program Gradbeništvo, Hidrotehniška

More information

IZRAČUN MEMBRANSKE RAZTEZNE POSODE - "MRP" za HLADNOVODNE SISTEME (DIN 4807/2)

IZRAČUN MEMBRANSKE RAZTEZNE POSODE - MRP za HLADNOVODNE SISTEME (DIN 4807/2) IZPIS IZRAČUN MEMBRANSKE RAZTEZNE POSODE - "MRP" za HLADNOVODNE SISTEME Izhodiščni podatki: Objkt : Vrtc Kamnitnik Projkt : PZI Uporaba MRP : Črpalna vrtina Datum : 30.8.2017 Obdlal : Zupan Skupna hladilna

More information

Univerza na Primorskem. Fakulteta za matematiko, naravoslovje in informacijske tehnologije. Zaznavanje gibov. Zaključna naloga

Univerza na Primorskem. Fakulteta za matematiko, naravoslovje in informacijske tehnologije. Zaznavanje gibov. Zaključna naloga Univerza na Primorskem Fakulteta za matematiko, naravoslovje in informacijske tehnologije Boštjan Markežič Zaznavanje gibov Zaključna naloga Koper, september 2011 Mentor: doc. dr. Peter Rogelj Kazalo Slovarček

More information

USING THE DIRECTION OF THE SHOULDER S ROTATION ANGLE AS AN ABSCISSA AXIS IN COMPARATIVE SHOT PUT ANALYSIS. Matej Supej* Milan Čoh

USING THE DIRECTION OF THE SHOULDER S ROTATION ANGLE AS AN ABSCISSA AXIS IN COMPARATIVE SHOT PUT ANALYSIS. Matej Supej* Milan Čoh Kinesiologia Slovenica, 14, 3, 5 14 (28) Faculty of Sport, University of Ljubljana, ISSN 1318-2269 5 Matej Supej* Milan Čoh USING THE DIRECTION OF THE SHOULDER S ROTATION ANGLE AS AN ABSCISSA AXIS IN COMPARATIVE

More information

Different Techniques for Calculating Apparent and Incremental Inductances using Finite Element Method

Different Techniques for Calculating Apparent and Incremental Inductances using Finite Element Method Different Techniques for Calculating Apparent and Incremental Inductances using Finite Element Method Dr. Amer Mejbel Ali Electrical Engineering Department Al-Mustansiriyah University Baghdad, Iraq amerman67@yahoo.com

More information

MECHANICAL EFFICIENCY, WORK AND HEAT OUTPUT IN RUNNING UPHILL OR DOWNHILL

MECHANICAL EFFICIENCY, WORK AND HEAT OUTPUT IN RUNNING UPHILL OR DOWNHILL original scientific article UDC: 796.4 received: 2011-05-03 MECHANICAL EFFICIENCY, WORK AND HEAT OUTPUT IN RUNNING UPHILL OR DOWNHILL Pietro Enrico DI PRAMPERO University of Udine, Department of Biomedical

More information

UNIVERZA V MARIBORU FAKULTETA ZA STROJNIŠTVO

UNIVERZA V MARIBORU FAKULTETA ZA STROJNIŠTVO UNIVERZA V MARIBORU FAKULTETA ZA STROJNIŠTVO Jan TIBAUT RAČUNSKA ANALIZA OBTEKANJA LOPATICE LOPATIČNE REŠETKE univerzitetnega študijskega programa 1. stopnje Strojništvo Maribor, september 2012 1 Fakulteta

More information

REGULACIJA ULTRASENZITIVNOSTI LINEARNO SKLOPLJENIH PROTEINSKIH KASKAD

REGULACIJA ULTRASENZITIVNOSTI LINEARNO SKLOPLJENIH PROTEINSKIH KASKAD REGULACIJA ULTRASENZITIVNOSTI LINEARNO SKLOPLJENIH PROTEINSKIH KASKAD Seminar iz fizike na dvopredmetnem študijskem programu Fizika (stari program) Aleš Vunjak Mentor: asist. dr. Rene Markovič Maribor,

More information

Acta Chim. Slov. 2000, 47, Macroion-macroion correlations in the presence of divalent counterions. Effects of a simple electrolyte B. Hrib

Acta Chim. Slov. 2000, 47, Macroion-macroion correlations in the presence of divalent counterions. Effects of a simple electrolyte B. Hrib Acta Chim. Slov. 2000, 47, 123-131 123 Macroion-macroion correlations in the presence of divalent counterions. Effects of a simple electrolyte B. Hribar and V. Vlachy Faculty of Chemistry and Chemical

More information

UNIVERZA NA PRIMORSKEM FAKULTETA ZA MATEMATIKO, NARAVOSLOVJE IN INFORMACIJSKE TEHNOLOGIJE. Verjetnostni algoritmi za testiranje praštevilskosti

UNIVERZA NA PRIMORSKEM FAKULTETA ZA MATEMATIKO, NARAVOSLOVJE IN INFORMACIJSKE TEHNOLOGIJE. Verjetnostni algoritmi za testiranje praštevilskosti UNIVERZA NA PRIMORSKEM FAKULTETA ZA MATEMATIKO, NARAVOSLOVJE IN INFORMACIJSKE TEHNOLOGIJE Zaključna naloga Verjetnostni algoritmi za testiranje praštevilskosti (Algorithms for testing primality) Ime in

More information

Computing the steady-state response of nonlinear circuits by means of the ǫ-algorithm

Computing the steady-state response of nonlinear circuits by means of the ǫ-algorithm Elektrotehniški vestnik XX(Y): 6, YEAR Electrotechnical Review, Ljubljana, Slovenija Computing the steady-state response of nonlinear circuits by means of the ǫ-algorithm Borut Wagner, Árpád Bűrmen, Janez

More information

Measurements of the magnetic field uniformity in a solenoid multilayer

Measurements of the magnetic field uniformity in a solenoid multilayer Elektrotehniški vestnik 74(5): 33-38, 27 Electrotechnical Review: Ljubljana, Slovenija Merjenje homogenosti magnetnega polja v večplastni solenoidni tuljavi Gregor Geršak Univerza v Ljubljani, Fakulteta

More information

Predmet: Seminar Avtor: Matic Pirc Mentor: prof. dr. Rudolf Podgornik

Predmet: Seminar Avtor: Matic Pirc Mentor: prof. dr. Rudolf Podgornik Fakulteta za matematiko in fiziko Univerza v Ljubljani MAVRICA Predmet: Seminar 2011 Avtor: Matic Pirc Mentor: prof. dr. Rudolf Podgornik Profesorja: dr. Martin Čopič in dr. Igor Poberaj Brežice, 29.4.2011

More information

1. UVOD UPORABA PROGRAMSKIH ORODIJ ZA DOLOČITEV NAPETOSTI KORAKA. Peter KITAK POVZETEK

1. UVOD UPORABA PROGRAMSKIH ORODIJ ZA DOLOČITEV NAPETOSTI KORAKA. Peter KITAK POVZETEK 5. posvetovanje "KOMUNALNA ENERGETIKA / POWER ENGINEERING", Maribor, 016 1 UPORABA PROGRAMSKIH ORODIJ ZA DOLOČITEV NAPETOSTI KORAKA Peter KITAK POVZETEK Na primeru paličnega ozemljila so v članku prikazani

More information

Makroekonomija 1: 4. vaje. Igor Feketija

Makroekonomija 1: 4. vaje. Igor Feketija Makroekonomija 1: 4. vaje Igor Feketija Teorija agregatnega povpraševanja AD = C + I + G + nx padajoča krivulja AD (v modelu AS-AD) učinek ponudbe denarja premiki vzdolž krivulje in premiki krivulje mikro

More information

UČNI NAČRT PREDMETA / COURSE SYLLABUS (leto / year 2017/18) Predmet: Analiza 3 Course title: Analysis 3. Študijska smer Study field ECTS

UČNI NAČRT PREDMETA / COURSE SYLLABUS (leto / year 2017/18) Predmet: Analiza 3 Course title: Analysis 3. Študijska smer Study field ECTS UČNI NAČRT PREDMETA / COURSE SYLLABUS (leto / year 2017/18) Predmet: Analiza 3 Course title: Analysis 3 Študijski program in stopnja Study programme and level Univerzitetni študijski program Matematika

More information

FIZIKA VIRUSOV. Avtor: Miran Dragar Mentor: prof. dr. Rudolf Podgornik. Maj Povzetek

FIZIKA VIRUSOV. Avtor: Miran Dragar Mentor: prof. dr. Rudolf Podgornik. Maj Povzetek UNIVERZA V LJUBLJANI Fakulteta za matematiko in fiziko Oddelek za fiziko FIZIKA VIRUSOV Avtor: Miran Dragar Mentor: prof. dr. Rudolf Podgornik Maj 2007 Povzetek V seminarju bo predstavljen preprost model,

More information

Aging Of Overvoltage Protection Elements Caused By Past Activations. Staranje elementov za prenaptostno zaščito zaradi premostitev

Aging Of Overvoltage Protection Elements Caused By Past Activations. Staranje elementov za prenaptostno zaščito zaradi premostitev Professional article Aging Of Overvoltage Protection Elements Caused By Past Activations Aleksandra Vasić 1, Miloš Vujisić 2, Koviljka Stanković 2, Predrag Osmokrović 2 * 1 Faculty of Mechanical Engineering,

More information

Univerza v Ljubljani, Fakulteta za matematiko in fiziko. Seminar

Univerza v Ljubljani, Fakulteta za matematiko in fiziko. Seminar Univerza v Ljubljani, Fakulteta za matematiko in fiziko Seminar Disperzijski modeli za modeliranje izpustov Avtor: Maruška Mole Mentor: asist. Rahela Žabkar Ljubljana, februar 2009 Povzetek Seminar predstavi

More information

Praktikum II. M. Zgonik, M. Horvat. 4. november 2013

Praktikum II. M. Zgonik, M. Horvat. 4. november 2013 M. Zgonik, M. Horvat 4. november 2013 Predgovor 1. Predgovor Trenutna skripta navodil k nalogam iz Praktikuma 2 je nastala po predlogi navodil prejšnjih let, od 1985 do 2002, ki jih je spisal Ivan Zupančič.

More information

UNIVERZA V MARIBORU FAKULTETA ZA NARAVOSLOVJE IN MATEMATIKO. Oddelek za matematiko in računalništvo DIPLOMSKO DELO.

UNIVERZA V MARIBORU FAKULTETA ZA NARAVOSLOVJE IN MATEMATIKO. Oddelek za matematiko in računalništvo DIPLOMSKO DELO. UNIVERZA V MARIBORU FAKULTETA ZA NARAVOSLOVJE IN MATEMATIKO Oddelek za matematiko in računalništvo DIPLOMSKO DELO Sabina Skornšek Maribor, 2012 UNIVERZA V MARIBORU FAKULTETA ZA NARAVOSLOVJE IN MATEMATIKO

More information

Statistika 2 z računalniško analizo podatkov. Neizpolnjevanje predpostavk regresijskega modela

Statistika 2 z računalniško analizo podatkov. Neizpolnjevanje predpostavk regresijskega modela Statistika 2 z računalniško analizo podatkov Neizpolnjevanje predpostavk regresijskega modela 1 Predpostavke regresijskega modela (ponovitev) V regresijskem modelu navadno privzamemo naslednje pogoje:

More information

EVA MARKELJ RAČUNALNIŠKO SIMULIRANJE SIPANJA SVETLOBE V ATMOSFERI

EVA MARKELJ RAČUNALNIŠKO SIMULIRANJE SIPANJA SVETLOBE V ATMOSFERI UNIVERZA V LJUBLJANI PEDAGOŠKA FAKULTETA EVA MARKELJ RAČUNALNIŠKO SIMULIRANJE SIPANJA SVETLOBE V ATMOSFERI DIPLOMSKO DELO LJUBLJANA, 2016 UNIVERZA V LJUBLJANI PEDAGOŠKA FAKULTETA DVOPREDMETNI UČITELJ:

More information

UNIVERZA V LJUBLJANI FAKULTETA ZA MATEMATIKO IN FIZIKO ODDELEK ZA FIZIKO SEMINAR. Pulzni eksperiment

UNIVERZA V LJUBLJANI FAKULTETA ZA MATEMATIKO IN FIZIKO ODDELEK ZA FIZIKO SEMINAR. Pulzni eksperiment UNIVERZA V LJUBLJANI FAKULTETA ZA MATEMATIKO IN FIZIKO ODDELEK ZA FIZIKO SEMINAR Pulzni eksperiment AVTOR: Andraž Petrović MENTOR: prof. Matjaž Ravnik Ljubljana, Maj 2004 POVZETEK: V seminarju bom opisal

More information

Osnove numerične matematike

Osnove numerične matematike Univerza v Ljubljani Fakulteta za računalništvo in informatiko Osnove numerične matematike Bojan Orel Ljubljana, 2004 Kazalo 1 Uvod 1 1.1 Zakaj numerične metode..................... 1 1.2 Napake in numerično

More information

UNIVERZA V MARIBORU FAKULTETA ZA NARAVOSLOVJE IN MATEMATIKO. Oddelek za matematiko in računalništvo MAGISTRSKA NALOGA. Tina Lešnik

UNIVERZA V MARIBORU FAKULTETA ZA NARAVOSLOVJE IN MATEMATIKO. Oddelek za matematiko in računalništvo MAGISTRSKA NALOGA. Tina Lešnik UNIVERZA V MARIBORU FAKULTETA ZA NARAVOSLOVJE IN MATEMATIKO Oddelek za matematiko in računalništvo MAGISTRSKA NALOGA Tina Lešnik Maribor, 2014 UNIVERZA V MARIBORU FAKULTETA ZA NARAVOSLOVJE IN MATEMATIKO

More information

2A skupina zemeljskoalkalijske kovine

2A skupina zemeljskoalkalijske kovine 1. NALOGA: V ČEM SE RAZLIKUJETA BeO IN MgO? 1. NALOGA: ODGOVOR Elementi 2. periode (od Li do F) se po fizikalnih in kemijskih lastnostih (diagonalne lastnosti) znatno razlikujejo od elementov, ki so v

More information

NUMERICAL SIMULATION OF THE PROGRESSIVE DAMAGE TO FRC PANELS DUE TO SHOCK LOADING

NUMERICAL SIMULATION OF THE PROGRESSIVE DAMAGE TO FRC PANELS DUE TO SHOCK LOADING UDK 539.4:519.6 ISSN 1580-949 Izvirni znanstveni ~lanek MTAEC9, 39(3)77(005) NUMERICAL SIMULATION OF THE PROGRESSIVE DAMAGE TO FRC PANELS DUE TO SHOCK LOADING NUMERI^NI MODEL NARA[^AJO^E PO[KODBE FRC-PANELOV

More information

Katastrofalno zaporedje okvar v medsebojno odvisnih omrežjih

Katastrofalno zaporedje okvar v medsebojno odvisnih omrežjih Katastrofalno zaporedje okvar v medsebojno odvisnih omrežjih Daniel Grošelj Mentor: Prof. Dr. Rudi Podgornik 2. marec 2011 Kazalo 1 Uvod 2 2 Nekaj osnovnih pojmov pri teoriji omrežij 3 2.1 Matrika sosednosti.......................................

More information

Assessment of surface deformation with simultaneous adjustment with several epochs of leveling networks by using nd relative pedaloid

Assessment of surface deformation with simultaneous adjustment with several epochs of leveling networks by using nd relative pedaloid RMZ - Materials and Geoenvironment, Vol. 53, No. 3, pp. 315-321, 2006 315 Assessment of surface deformation with simultaneous adjustment with several epochs of leveling networks by using nd relative pedaloid

More information

Dušan Čalić. Meritev pozitivnega temperaturnega koeficienta reaktivnosti v reaktorju TRIGA SEMINAR

Dušan Čalić. Meritev pozitivnega temperaturnega koeficienta reaktivnosti v reaktorju TRIGA SEMINAR UNIVERZA V LJUBLJANI FAKULTETA ZA MATEMATIKO IN FIZIKO Oddelek za fiziko Dušan Čalić Meritev pozitivnega temperaturnega koeficienta reaktivnosti v reaktorju TRIGA SEMINAR Mentor: prof. dr. Matjaž Ravnik

More information

Vpliv delovanja napetostnega stabilizatorja MAGTECH na NN distribucijsko omrežje

Vpliv delovanja napetostnega stabilizatorja MAGTECH na NN distribucijsko omrežje 9. KONFERENCA SLOVENSKIH ELEKROENERGEIKOV Kranjska Gora 9 Vpliv delovanja napetostnega stabilizatorja MAGECH na NN distribucijsko omrežje Miran Rošer Elektro Celje d.d. Vrunčeva a, Celje E-mail: miran.roser@elektro-celje.si,

More information

MODELI CESTNEGA PROMETA

MODELI CESTNEGA PROMETA MODELI CESTNEGA PROMETA LUKA ŠEPEC Fakulteta za matematiko in fiziko Univerza v Ljubljani V članku so predstavljeni različni pristopi k modeliranju cestnega prometa. Najprej so predstavljene empirične

More information

UPORABA SIMOVERT VC MASTERJA ZA POGON TRANSPORTNEGA TRAKU

UPORABA SIMOVERT VC MASTERJA ZA POGON TRANSPORTNEGA TRAKU Alen Bračič UPORABA SIMOVERT VC MASTERJA ZA POGON TRANSPORTNEGA TRAKU Diplomsko delo Maribor, april 2010 II III Diplomsko delo visokošolskega strokovnega študijskega programa Študent: Študijski program:

More information

Računalniško simuliranje dinamike rotorjev Computer Simulation of the Dynamics of Rotors

Računalniško simuliranje dinamike rotorjev Computer Simulation of the Dynamics of Rotors STROJNIŠKI VESTNIK - JOURNAL OF MECHANICAL ENGINEERING, LJUBLJANA (42) 1996/9 10 1 Računalniško simuliranje dinamike rotorjev Computer Simulation of the Dynamics of Rotors Robert Cokan, Miha Boltežar,

More information

SOLARNI DIMNIK. Zaključni seminar na študijskem programu 1.stopnje Fizika MARTIN KRANER. Maribor, 2010

SOLARNI DIMNIK. Zaključni seminar na študijskem programu 1.stopnje Fizika MARTIN KRANER. Maribor, 2010 SOLARNI DIMNIK Zaključni seminar na študijskem programu 1.stopnje Fizika MARTIN KRANER V seminarju je predstavljen solarni dimnik, njegovi sestavni deli in delovanje. Kot primer ţe delujočega solarnega

More information