Methods and Models for Accelerating Dynamic Simulation of Fluid Power Circuits

The objective of this dissertation is to improve the dynamic simulation of fluid power circuits. A fluid power circuit is a typical way to implement power transmission in mobile working machines, e.g. cranes, excavators etc. Dynamic simulation is an essential tool in developing controllability and energy-efficient solutions for mobile machines. Efficient dynamic simulation is the basic requirement for the real-time simulation. In the real-time simulation of fluid power circuits there exist numerical problems due to the software and methods used for modelling and integration. A simulation model of a fluid power circuit is typically created using differential and algebraic equations. Efficient numerical methods are required since differential equations must be solved in real time. Unfortunately, simulation software packages offer only a limited selection of numerical solvers. Numerical problems cause noise to the results, which in many cases leads the simulation run to fail. Mathematically the fluid power circuit models are stiff systems of ordinary differential equations. Numerical solution of the stiff systems can be improved by two alternative approaches. The first is to develop numerical solvers suitable for solving stiff systems. The second is to decrease the model stiffness itself by introducing models and algorithms that either decrease the highest eigenvalues or neglect them by introducing steady-state solutions of the stiff parts of the models. The thesis proposes novel methods using the latter approach. The study aims to develop practical methods usable in dynamic simulation of fluid power circuits using explicit fixed-step integration algorithms. In this thesis, twomechanisms whichmake the systemstiff are studied. These are the pressure drop approaching zero in the turbulent orifice model and the volume approaching zero in the equation of pressure build-up. These are the critical areas to which alternative methods for modelling and numerical simulation are proposed. Generally, in hydraulic power transmission systems the orifice flow is clearly in the turbulent area. The flow becomes laminar as the pressure drop over the orifice approaches zero only in rare situations. These are e.g. when a valve is closed, or an actuator is driven against an end stopper, or external force makes actuator to switch its direction during operation. This means that in terms of accuracy, the description of laminar flow is not necessary. But, unfortunately, when a purely turbulent description of the orifice is used, numerical problems occur when the pressure drop comes close to zero since the first derivative of flow with respect to the pressure drop approaches infinity when the pressure drop approaches zero. Furthermore, the second derivative becomes discontinuous, which causes numerical noise and an infinitely small integration step when a variable step integrator is used. A numerically efficient model for the orifice flow is proposed using a cubic spline function to describe the flow in the laminar and transition areas. Parameters for the cubic spline function are selected such that its first derivative is equal to the first derivative of the pure turbulent orifice flow model in the boundary condition. In the dynamic simulation of fluid power circuits, a tradeoff exists between accuracy and calculation speed. This investigation is made for the two-regime flow orifice model. Especially inside of many types of valves, as well as between them, there exist very small volumes. The integration of pressures in small fluid volumes causes numerical problems in fluid power circuit simulation. Particularly in realtime simulation, these numerical problems are a great weakness. The system stiffness approaches infinity as the fluid volume approaches zero. If fixed step explicit algorithms for solving ordinary differential equations (ODE) are used, the system stability would easily be lost when integrating pressures in small volumes. To solve the problem caused by small fluid volumes, a pseudo-dynamic solver is proposed. Instead of integration of the pressure in a small volume, the pressure is solved as a steady-state pressure created in a separate cascade loop by numerical integration. The hydraulic capacitance V/Be of the parts of the circuit whose pressures are solved by the pseudo-dynamic method should be orders of magnitude smaller than that of those partswhose pressures are integrated. The key advantage of this novel method is that the numerical problems caused by the small volumes are completely avoided. Also, the method is freely applicable regardless of the integration routine applied. The superiority of both above-mentioned methods is that they are suited for use together with the semi-empirical modelling method which necessarily does not require any geometrical data of the valves and actuators to be modelled. In this modelling method, most of the needed component information can be taken from the manufacturer’s nominal graphs. This thesis introduces the methods and shows several numerical examples to demonstrate how the proposed methods improve the dynamic simulation of various hydraulic circuits.

Symmetristen komponenttien mallinnus

Tässä insinöörityössä suunniteltiin Helsingin ammattikorkeakoululle jakeluverkoissa tapahtuvien oikosulkujen symmetristen komponenttien laskutavan havainnollistamiseen sopiva sijaiskytkentä. Sijaiskytkennässä tärkeitä huomioitavia asioita olivat mm. jännitetaso, havainnollistavien muuntajien oikosulkukestoisuus, jatkokäyttö laboratoriotyönä ja yleinen havainnollistavuus. Työssä on aluksi perehdytty symmetristen komponenttien ja jakeluverkoissa tapahtuvien oikosulkujen teoriaan. Tämän jälkeen mitoitettiin tarvittavan kytkennän komponenttien jännite- ja virtakestoisuudet mahdolliset lisäkäytöt huomioiden. Näiden rajoitusten mukaan perusteella työtä ruvettiin toteuttamaan. Työssä tilattiin sähkön 40 V:n pääjännitetasolle alentava muuntaja syöttämään oikosulun kestävää muuntajaa, jolla simuloitiin jakeluverkon yleisimpiä vikatyyppejä. Jälkimmäiselle muuntajalle mitoitettiin ja hankittiin sisäistä impedanssia vastaava induktanssi. Tämän avulla rakennettiin kokonaisuus, jonka avulla voidaan simuloida kaikkia tapahtuvia oikosulkuja vastaavat sijaiskytkennät. Työhön jätettiin kehittämisvaraa ja muita laboratoriotyön rakentamismahdollisuuksia tulevien insinööritöiden tekijöille.

Re-usable ASIC integrationverification software architecture

Mikropiirien valmistus- ja suunnittelutekniikoiden kehittyminen mahdollistaa yhä monimutkaisempien mikropiirien valmistamisen. Piirien verifioinnista onkin tullut prosessin aikaa vievin osa,sillä kompleksisuuden kasvaessa kasvaa verifioinnin tarve eksponentiaalisesti. Vaikka erinäisiä strategioita piirien integroinnin verifiointiin on esitetty, mm. verifioinnin jakaminen koko suunnitteluprosessin ajalle, jopa yli puolet koko piirin suunnitteluun ja valmistukseen käytetystä työmäärästä kuluu verifiointiin. Uudelleenkäytettävät komponentit ovat pääosassa piirin suunnittelussa, mutta verifioinnissa uudelleenkäytettävyyttä ei ole otettu kunnolla käyttöön ainakaan verifiointiohjelmistojen osalta. Tämä diplomityö esittelee uudelleenkäytettävän mikropiirien verifiointiohjelmistoarkkitehtuurin, jolla saadaan verifiointitaakkaa vähennettyä poistamalla verifioinnissa käytettävien ohjelmistojen uudelleensuunnittelun ja toteuttamisen tarvetta.

Suurnopeushilaohjain IGB-transistoreille

Hilalta ohjattavien tehokytkimien, kuten IGBT tai MOSFET, ohjaaminen toteutetaan ohjausjännitettä muuttamalla. Hilaohjaimen tehtävä on ladata tehokytkimen hilakapasitanssia ohjattuun jännitteeseen. Hilaohjaimen siirtämän varauksen suuruuteen vaikuttaa hilakapasitanssin todellinen suuruus ja käytettävä ohjausjännite. Hilaohjaimen virrankäsittelykyky määrittää hilaohjaimen nopeuden siinä mielessä, että ohjattavan komponentin hilavaraus tulisi pystyä lataamaan ja purkamaan tietyssä suhteessa kytkentätaajuuteen. Kytkentätaajuuksien kasvaessa myös hilaohjaimen virrankäsittelykyky sekä hilaohjaimen teho joutuvat uudelleen arvioitaviksi. Työssä perehdytään IGBT:n toimintaan ja hilaohjainratkaisuihin suurilla kytkentätaajuuksilla. Työssä ontutkittu mahdollisuutta rakentaa IGBT:n hilaohjain jo markkinoilla olevista MOSFET:ien hilaohjaimista. Suurnopeushilaohjaimesta suunnitellaanja rakennetaan prototyyppi. Prototyypin suoriutumista tehtävistään tarkastellaan tekemällä mittauksia. Lopuksi arvioidaan asioita joiden tulee muuttua matkalla kohti megahertsien kytkentätaajuuksia.

Sarjakytkentäisen taajuusmuuttajan ohjauksen toteutus FPGA-piireillä

Tässä diplomityössä tutustutaan sarjakytkentäisen taajuusmuuttajan rakenteeseen ja ohjaukseen. Työssä tarkastellaan myös hajautetun järjestelmän testaamista ja ohjelmointia. Työssä toteutettiin sarjakytkentäisen taajuusmuuttajan modulaattori hajautetusti FPGA-piireille. Lisäksi diplomityössä kehitettiin ja toteutettiin optinen JTAG-rajapinta hajautetun järjestelmän testaukseen ja ohjelmointiin. Laboratoriokoelaitteisto koostui yhdeksästä taajuusmuuttajan ohjauskortista, LUT Master -ohjausyksiköstä ja dSPACE-ohjausjärjestelmästä. Laboratoriomittaukset tehtiin Lappeenrannan teknillisen yliopiston Säätötekniikan laboratoriossa.

Verkkovaihtosuuntaajan vektorisäädön toteutus FPGA-piirillä

Verkkovaihtosuuntaajalla pystytään muuntamaan tasajännite vaihtojännitteeksi ja päinvastoin. Verkkovaihtosuuntaajan toiminta perustuu tehokytkinten ohjaukseen ja sopivan modulointimenetelmän käyttöön. Vektorisäädössä vaihtosuuntaajanvirrat ja jännitteet esitetään kompleksitasossa, jolloin virta- ja jännitekomponentit voidaan esittää vektoreina. Vektorisäädössä verkkovaihtosuuntaajan ohjaustoteutetaan laskemalla kompleksitasossa vektoreille arvot, jotka tuottavat vaihtosuuntaajan lähtöön halutun vektorin. Koska FPGA-piirit mahdollistavat nopean rinnakkaisen laskennan, soveltuvat ne hyvin vektorisäädön toteuttamiseen. FPGA-piirien rakenteesta johtuen on säätöjärjestelmän suunnittelussa huomioitava kiinteän pilkun lukujen riittävä bittileveys ja järjestelmän diskretointiaika. Työssä suunnitellaan verkkovaihtosuuntaajan vektorisäätö ja tutkitaan bittileveyden vaikutusta säädön toteuttamiseen FPGA-piirillä. Bittileveyden tarkasteluun esitetään käytettäväksi tilastollisia menetelmiä. Työssä tarkastellaan kiinteän pilkun järjestelmän ja liukulukujärjestelmän erosuureen tilastollisia tunnusmerkkejä sekä histogrammia. Tarkasteluissa huomattiin, että maksimivirhe itsessään ei tarjoa riittävästi tietoa erosuureen jakautumisesta. Näin ollen maksimivirhe ei ole kaikissa tilanteissa sovelias menetelmä riittävän bittitarkkuuden määrittämiseen. Työssä esitetään riittävän bittitarkkuuden määrittelemiseen käytettäväksi otossuureista otosvarianssia, keskipoikkeamaa ja vaihteluväliä.

Hydraulisen kuristinmallin ja liikealustan ohjauksen kehittäminen reaaliaikasimulointiin

Diplomityön ensimmäisessä vaiheessa tutkittiin hydraulisen kuristimen ominaisuuksia ja esiteltiin numeerisesti tehokas kuristinmalli käyttäenpolynomifunktiota virtauksen laminaarisen ja transitioalueen kuvaukseen. Puoliempiirisen mallin paremmuus tulee esille siinä, että kuristimen geometriatietoja ei tarvita laskettaessa virtausta paine-eron perusteella. Reaaliaikasimuloinnissa esiintyy kompromisseja tarkkuuden ja laskentanopeuden välillä. Tätä asiaa tutkittiin kahden virtausalueen kuristinmallilla. Transition paine-eron ja integrointiaika-askelen valinnan vaikutus tarkkuuteen ja laskentanopeuteen tutkittiin. Toisessa vaiheessa tutkittiin mahdollisimman hyvän liiketuntuman tuottamista liikealustalla ohjaussignaalia kehittämällä. Liikealustan liikeradan rajallisuudesta johtuen ohjauksessa on perinteisesti käytetty washout-suodatusta, joka erottelee simuloitavan järjestelmän kiihtyvyyssignaalista vain nopeatkiihtyvyydet. Tässä työssä tutkittiin hitaiden kiihtyvyyksien ottamista mukaan liikealustan ohjaukseen liikealustan liikeradan puitteissa. Tämä toteutettiin kuvaamalla hitaat kiihtyvyydet kallistamalla liikealustaa, jolloin käyttäjään kohdistuva voima saatiin kuvattua gravitaatiota hyväksi käyttäen.

Virranmittausanturin suunnittelu ja simulointi taajuusmuuttajakäyttöön

Virranmittausantureita tarvitaan monenlaisissa käyttökohteissa, joissa ne mittaavat sekä virran suuruuttaettä laatua ja toimivat osana niiden säätelyjärjestelmää. Virranmittausantureita tarvitaan myös vikatilanteiden määrittämiseen erilaisissa suojauspiireissä. Taajuusmuuttajissa virranmittaus on hyvin tärkeää ja suurista virroista sekä taajuuksista johtuen se täytyy suunnitella huolella. Tässä diplomityössä käsitellään ja tutkitaan eri virranmittausmenetelmiä, joiden avulla taajuusmuuttajan luotettava virranmittaus voidaan toteuttaa. Työssä tutkitaan eri menetelmiä virranmittauksen toteuttamiseksi, minkä jälkeen niistä valitaan sopiva menetelmä ja tutkitaan sen eri toteutusvaihtoehtoja. Sopivan toteutusvaihtoehdon valinnan jälkeen työssä suunnitellaan oma virranmittausanturi, joka sopii nimenomaisesti taajuusmuuttajakäyttöön. Suunnitellun anturin ominaisuuksia tutkitaan lopuksi simuloimalla, jonka jälkeen arvioidaan sen soveltumista käytännön sovelluksiin sekä arvioidaan erilaisia keinoja sen parantamiseksi.

FPGA-Based Control Design for Power Electronic Applications

Tehoelektoniikkalaitteella tarkoitetaan ohjaus- ja säätöjärjestelmää, jolla sähköä muokataan saatavilla olevasta muodosta haluttuun uuteen muotoon ja samalla hallitaan sähköisen tehon virtausta lähteestä käyttökohteeseen. Tämä siis eroaa signaalielektroniikasta, jossa sähköllä tyypillisesti siirretään tietoa hyödyntäen eri tiloja. Tehoelektroniikkalaitteita vertailtaessa katsotaan yleensä niiden luotettavuutta, kokoa, tehokkuutta, säätötarkkuutta ja tietysti hintaa. Tyypillisiä tehoelektroniikkalaitteita ovat taajuudenmuuttajat, UPS (Uninterruptible Power Supply) -laitteet, hitsauskoneet, induktiokuumentimet sekä erilaiset teholähteet. Perinteisesti näiden laitteiden ohjaus toteutetaan käyttäen mikroprosessoreja, ASIC- (Application Specific Integrated Circuit) tai IC (Intergrated Circuit) -piirejä sekä analogisia säätimiä. Tässä tutkimuksessa on analysoitu FPGA (Field Programmable Gate Array) -piirien soveltuvuutta tehoelektroniikan ohjaukseen. FPGA-piirien rakenne muodostuu erilaisista loogisista elementeistä ja niiden välisistä yhdysjohdoista.Loogiset elementit ovat porttipiirejä ja kiikkuja. Yhdysjohdot ja loogiset elementit ovat piirissä kiinteitä eikä koostumusta tai lukumäärää voi jälkikäteen muuttaa. Ohjelmoitavuus syntyy elementtien välisistä liitännöistä. Piirissä on lukuisia, jopa miljoonia kytkimiä, joiden asento voidaan asettaa. Siten piirin peruselementeistä voidaan muodostaa lukematon määrä erilaisia toiminnallisia kokonaisuuksia. FPGA-piirejä on pitkään käytetty kommunikointialan tuotteissa ja siksi niiden kehitys on viime vuosina ollut nopeaa. Samalla hinnat ovat pudonneet. Tästä johtuen FPGA-piiristä on tullut kiinnostava vaihtoehto myös tehoelektroniikkalaitteiden ohjaukseen. Väitöstyössä FPGA-piirien käytön soveltuvuutta on tutkittu käyttäen kahta vaativaa ja erilaista käytännön tehoelektroniikkalaitetta: taajuudenmuuttajaa ja hitsauskonetta. Molempiin testikohteisiin rakennettiin alan suomalaisten teollisuusyritysten kanssa soveltuvat prototyypit,joiden ohjauselektroniikka muutettiin FPGA-pohjaiseksi. Lisäksi kehitettiin tätä uutta tekniikkaa hyödyntävät uudentyyppiset ohjausmenetelmät. Prototyyppien toimivuutta verrattiin vastaaviin perinteisillä menetelmillä ohjattuihin kaupallisiin tuotteisiin ja havaittiin FPGA-piirien mahdollistaman rinnakkaisen laskennantuomat edut molempien tehoelektroniikkalaitteiden toimivuudessa. Työssä on myösesitetty uusia menetelmiä ja työkaluja FPGA-pohjaisen säätöjärjestelmän kehitykseen ja testaukseen. Esitetyillä menetelmillä tuotteiden kehitys saadaan mahdollisimman nopeaksi ja tehokkaaksi. Lisäksi työssä on kehitetty FPGA:n sisäinen ohjaus- ja kommunikointiväylärakenne, joka palvelee tehoelektroniikkalaitteiden ohjaussovelluksia. Uusi kommunikointirakenne edistää lisäksi jo tehtyjen osajärjestelmien uudelleen käytettävyyttä tulevissa sovelluksissa ja tuotesukupolvissa.

Applicability of Power-Line Communications to Data Transfer of On-line Condition Monitoring of Electrical Drives

The research of power-line communications has been concentrated on home automation, broadband indoor communications and broadband data transfer in a low voltage distribution network between home andtransformer station. There has not been carried out much research work that is focused on the high frequency characteristics of industrial low voltage distribution networks. The industrial low voltage distribution network may be utilised as a communication channel to data transfer required by the on-line condition monitoring of electric motors. The advantage of using power-line data transfer is that it does not require the installing of new cables. In the first part of this work, the characteristics of industrial low voltage distribution network components and the pilot distribution network are measured and modelled with respect topower-line communications frequencies up to 30 MHz. The distributed inductances, capacitances and attenuation of MCMK type low voltage power cables are measured in the frequency band 100 kHz - 30 MHz and an attenuation formula for the cables is formed based on the measurements. The input impedances of electric motors (15-250 kW) are measured using several signal couplings and measurement based input impedance model for electric motor with a slotted stator is formed. The model is designed for the frequency band 10 kHz - 30 MHz. Next, the effect of DC (direct current) voltage link inverter on power line data transfer is briefly analysed. Finally, a pilot distribution network is formed and signal attenuation in communication channels in the pilot environment is measured. The results are compared with the simulations that are carried out utilising the developed models and measured parameters for cables and motors. In the second part of this work, a narrowband power-line data transfer system is developed for the data transfer ofon-line condition monitoring of electric motors. It is developed using standardintegrated circuits. The system is tested in the pilot environment and the applicability of the system for the data transfer required by the on-line condition monitoring of electric motors is analysed.

Torque vibration model of axial-flux surface-mounted permanent magnet synchronous machine

In order that the radius and thus ununiform structure of the teeth and otherelectrical and magnetic parts of the machine may be taken into consideration the calculation of an axial flux permanent magnet machine is, conventionally, doneby means of 3D FEM-methods. This calculation procedure, however, requires a lotof time and computer recourses. This study proves that also analytical methods can be applied to perform the calculation successfully. The procedure of the analytical calculation can be summarized into following steps: first the magnet is divided into slices, which makes the calculation for each section individually, and then the parts are submitted to calculation of the final results. It is obvious that using this method can save a lot of designing and calculating time. Thecalculation program is designed to model the magnetic and electrical circuits of surface mounted axial flux permanent magnet synchronous machines in such a way, that it takes into account possible magnetic saturation of the iron parts. Theresult of the calculation is the torque of the motor including the vibrations. The motor geometry and the materials and either the torque or pole angle are defined and the motor can be fed with an arbitrary shape and amplitude of three-phase currents. There are no limits for the size and number of the pole pairs nor for many other factors. The calculation steps and the number of different sections of the magnet are selectable, but the calculation time is strongly depending on this. The results are compared to the measurements of real prototypes. The permanent magnet creates part of the flux in the magnetic circuit. The form and amplitude of the flux density in the air-gap depends on the geometry and material of the magnetic circuit, on the length of the air-gap and remanence flux density of the magnet. Slotting is taken into account by using the Carter factor in the slot opening area. The calculation is simple and fast if the shape of the magnetis a square and has no skew in relation to the stator slots. With a more complicated magnet shape the calculation has to be done in several sections. It is clear that according to the increasing number of sections also the result will become more accurate. In a radial flux motor all sections of the magnets create force with a same radius. In the case of an axial flux motor, each radial section creates force with a different radius and the torque is the sum of these. The magnetic circuit of the motor, consisting of the stator iron, rotor iron, air-gap, magnet and the slot, is modelled with a reluctance net, which considers the saturation of the iron. This means, that several iterations, in which the permeability is updated, has to be done in order to get final results. The motor torque is calculated using the instantaneous linkage flux and stator currents. Flux linkage is called the part of the flux that is created by the permanent magnets and the stator currents passing through the coils in stator teeth. The angle between this flux and the phase currents define the torque created by the magnetic circuit. Due to the winding structure of the stator and in order to limit the leakage flux the slot openings of the stator are normally not made of ferromagnetic material even though, in some cases, semimagnetic slot wedges are used. In the slot opening faces the flux enters the iron almost normally (tangentially with respect to the rotor flux) creating tangential forces in the rotor. This phenomenon iscalled cogging. The flux in the slot opening area on the different sides of theopening and in the different slot openings is not equal and so these forces do not compensate each other. In the calculation it is assumed that the flux entering the left side of the opening is the component left from the geometrical centre of the slot. This torque component together with the torque component calculated using the Lorenz force make the total torque of the motor. It is easy to assume that when all the magnet edges, where the derivative component of the magnet flux density is at its highest, enter the slot openings at the same time, this will have as a result a considerable cogging torque. To reduce the cogging torquethe magnet edges can be shaped so that they are not parallel to the stator slots, which is the common way to solve the problem. In doing so, the edge may be spread along the whole slot pitch and thus also the high derivative component willbe spread to occur equally along the rotation. Besides forming the magnets theymay also be placed somewhat asymmetric on the rotor surface. The asymmetric distribution can be made in many different ways. All the magnets may have a different deflection of the symmetrical centre point or they can be for example shiftedin pairs. There are some factors that limit the deflection. The first is that the magnets cannot overlap. The magnet shape and the relative width compared to the pole define the deflection in this case. The other factor is that a shifting of the poles limits the maximum torque of the motor. If the edges of adjacent magnets are very close to each other the leakage flux from one pole to the other increases reducing thus the air-gap magnetization. The asymmetric model needs some assumptions and simplifications in order to limit the size of the model and calculation time. The reluctance net is made for symmetric distribution. If the magnets are distributed asymmetrically the flux in the different pole pairs will not be exactly the same. Therefore, the assumption that the flux flows from the edges of the model to the next pole pairs, in the calculation model from one edgeto the other, is not correct. If it were wished for that this fact should be considered in multi-pole pair machines, this would mean that all the poles, in other words the whole machine, should be modelled in reluctance net. The error resulting from this wrong assumption is, nevertheless, irrelevant.

Design and implementation of FPGA-based LQ control of active magnetic bearings

The need for high performance, high precision, and energy saving in rotating machinery demands an alternative solution to traditional bearings. Because of the contactless operation principle, the rotating machines employing active magnetic bearings (AMBs) provide many advantages over the traditional ones. The advantages such as contamination-free operation, low maintenance costs, high rotational speeds, low parasitic losses, programmable stiffness and damping, and vibration insulation come at expense of high cost, and complex technical solution. All these properties make the use of AMBs appropriate primarily for specific and highly demanding applications. High performance and high precision control requires model-based control methods and accurate models of the flexible rotor. In turn, complex models lead to high-order controllers and feature considerable computational burden. Fortunately, in the last few years the advancements in signal processing devices provide new perspective on the real-time control of AMBs. The design and the real-time digital implementation of the high-order LQ controllers, which focus on fast execution times, are the subjects of this work. In particular, the control design and implementation in the field programmable gate array (FPGA) circuits are investigated. The optimal design is guided by the physical constraints of the system for selecting the optimal weighting matrices. The plant model is complemented by augmenting appropriate disturbance models. The compensation of the force-field nonlinearities is proposed for decreasing the uncertainty of the actuator. A disturbance-observer-based unbalance compensation for canceling the magnetic force vibrations or vibrations in the measured positions is presented. The theoretical studies are verified by the practical experiments utilizing a custom-built laboratory test rig. The test rig uses a prototyping control platform developed in the scope of this work. To sum up, the work makes a step in the direction of an embedded single-chip FPGA-based controller of AMBs.

Kiinnirullaimen dynamiikan ja värähtelyjen simulointi

Diplomityössä tutkittiin kaupallisen monikappaledynamiikkaohjelmiston soveltuvuutta kiinnirullaimen dynamiikan ja värähtelyjen tutkimiseen. Erityisen kiinnostuneita oltiin nipin kuvauksesta sekä nipissä tapahtuvista värähtelyistä. Tässä diplomityössä mallinnettiin kiinnirullaimen ensiö- ja toisiokäytöt sekä tampuuritela. Malli yhdistettiin myöhemmin Metso Paper Järvenpäässä rinnakkaisena diplomityönä tehtyyn malliin, joista muodostui kahteen ratkaisijaan perustuva simulointimalli. Simulointimalli rakennettiin käyttämään kahta erillistä ratkaisijaa, joista toinen on mekaniikkamallin rakentamisessa käytetty ADAMS-ohjelmisto ja toinen säätöjärjestelmää ja hydraulipiirejä kuvaava Simulink-malli. Nipin mallintamiseksi tampuuritela ja rullaussylinteri mallinnettiin joustaviksi käyttäen keskitettyjen massojen menetelmää. Siirtolaitteissa sekä runkorakenteissa tapahtuvat joustot kuvattiin yhden vapausasteen jousi-vaimennin voimilla kuvattuina järjestelminä. Tässä diplomityössä on myös keskitytty esittelemään ADAMS-ohjelmiston toimintaa ohjeistavasti sekä käsittelemään parametrisen mallintamisen etuja. Työssä havaittiin monikappaledynamiikan soveltuvuus kiinnirullaimen dynamiikan sekä dynaamisten voimien aiheuttamien värähtelyjen tutkimiseen. Suoritetuista värähtelymittauksista voitiin tehdä vain arvioita. Mallin havaittiin vaativan lisätutkimusta ja kehitystyötä

Rullaussylinterin ja painotelan mallinnus kiinnirullaimen simulointiin

Diplomityössä tutkitaan kaupallisen simulointiohjelmiston soveltuvuutta nykyaikaisen kiinnirullaimen dynamiikan tutkimiseen. Kiinnostuksen kohteena on erityisesti kahden telan välinen nippi, sekä siinä tapahtuvat värähtelyt. Työssä mallinnetaan rullaussylinterin ja telapainolaitteen simulointimallit. Rullaussylinterin simulointimalli yhdistetään Lappeenrannan teknillisessä korkeakoulussa mallinnettuun tampuuritelan simulointimalliin, jolloin nippikontaktin tutkiminen on mahdollista. Simuloituja tuloksia verrataan todellisella laitteella tehtyihin mittauksiin sekä elementtimenetelmällä laskettuihin tuloksiin. Diplomityön mekaniikka mallinnetaan ADAMS-ohjelmistossa monikappaledynamiikan keinoin. Toimilaitteiden sekä säätöjärjestelmien kuvaukseen käytetään MATLAB Simulink-ohjelmistoa. Telojen joustavuuden mallinnuksessa käytetään hyväksi keskittyneiden massojen periaatetta. Järjestelmän hydraulipiirit mallinnetaan keskittyneiden paineiden teorian mukaisesti ja toimilaitteiden mallinnuksessa käytetään puoliempiiristä mallinnustekniikkaa. Työssä havaitaan monikappaledynamiikan soveltuvan kiinnirullaimen dynamiikan tutkimiseen. Kahden diplomityön tuloksena laaditun nippimallin avulla voidaan kuvata rullaustapahtumassa vaikuttavat voimat oikein. Värähtelymittausten perusteella voidaan tehdä karkeita johtopäätöksiä, mallin toimivuuden arvioimiseksi värähtelyjen kuvaamisessa, joskin mallin havaitaan vaativan lisätutkimusta ja kehitystyötä.

Anaerobic treatment of deinking pulp plant effluents

Anaerobic treatment as a first biological stage in wastewater treatment is nowadays a well-established technology in recycled paper processing mills using closed water circuits. Today further developed high-rate processes and especially high-tower reactors are also able to handle lower organic loads and become therefore feasible for deinking pulp plant effluents. The interest in the anaerobic method is based on a positive energy balance in form of biogas production and low biomass yield from the process. The anaerobic treatment method was researched and its suitability for the deinking pulp plant effluents was tested experimentally at Stora Enso Maxau mill. In the theory, the deinking pulp process is introduced and the effluents from the deinking process are characterized. The anaerobic treatment is brought up in depth in terms of its use for the deinking effluents, and different kind of reactor types are presented. In addition, other wastewater treatment methods are shortly introduced with the focus on tertiary treatment. Static biodegradability tests were carried out for the wastewaters both anaerobically and aerobically. Based on the results, the deinking effluents can be degraded anaerobically, and inhibition to the methanogenic bacteria was not noticed. In the aerobic static test a good performance of the existing wastewater treatment plant at Maxau mill was proved. Later on pilot trials with sequential anaerobic-aerobic treatment were carried out for the deinking effluents. The anaerobic reactor used was a so called internal circulation reactor. The results confirmed that the combination of the anaerobic treatment and the aerobic activated sludge process is a suitable method for deinking wastewaters with a COD reduction as good as with a two stage aerobic method. When combined with the outstanding quality of the produced biogas and the cost savings acquired from the lower sludge production, the anaerobic treatment was found to be an especially favorable treatment method.