Dynamic Simulations of Rotors during Drop on Retainer Bearings

The active magnetic bearings have recently been intensively developed because of noncontact support having several advantages compared to conventional bearings. Due to improved materials, strategies of control, and electrical components, the performance and reliability of the active magnetic bearings are improving. However, additional bearings, retainer bearings, still have a vital role in the applications of the active magnetic bearings. The most crucial moment when the retainer bearings are needed is when the rotor drops from the active magnetic bearings on the retainer bearings due to component or power failure. Without appropriate knowledge of the retainer bearings, there is a chance that an active magnetic bearing supported rotor system will be fatal in a drop-down situation. This study introduces a detailed simulation model of a rotor system in order to describe a rotor drop-down situation on the retainer bearings. The introduced simulation model couples a finite element model with component mode synthesis and detailed bearing models. In this study, electrical components and electromechanical forces are not in the focus. The research looks at the theoretical background of the finite element method with component mode synthesis that can be used in the dynamic analysis of flexible rotors. The retainer bearings are described by using two ball bearing models, which include damping and stiffness properties, oil film, inertia of rolling elements and friction between races and rolling elements. Thefirst bearing model assumes that the cage of the bearing is ideal and that the cage holds the balls in their predefined positions precisely. The second bearing model is an extension of the first model and describes the behavior of the cageless bearing. In the bearing model, each ball is described by using two degrees of freedom. The models introduced in this study are verified with a corresponding actual structure. By using verified bearing models, the effects of the parameters of the rotor system onits dynamics during emergency stops are examined. As shown in this study, the misalignment of the retainer bearings has a significant influence on the behavior of the rotor system in a drop-down situation. In this study, a stability map of the rotor system as a function of rotational speed of the rotor and the misalignment of the retainer bearings is presented. In addition, the effects of parameters of the simulation procedure and the rotor system on the dynamics of system are studied.

A novel procedure for establishing clinching parameters for high strength steel sheet

The building industry has a particular interest in using clinching as a joining method for frame constructions of light-frame housing. Normally many clinch joints are required in joining of frames.In order to maximise the strength of the complete assembly, each clinch joint must be as sound as possible. Experimental testing is the main means of optimising a particular clinch joint. This includes shear strength testing and visual observation of joint cross-sections. The manufacturers of clinching equipment normally perform such experimental trials. Finite element analysis can also be used to optimise the tool geometry and the process parameter, X, which represents the thickness of the base of the joint. However, such procedures require dedicated software, a skilled operator, and test specimens in order to verify the finite element model. In addition, when using current technology several hours' computing time may be necessary. The objective of the study was to develop a simple calculation procedure for rapidly establishing an optimum value for the parameter X for a given tool combination. It should be possible to use the procedure on a daily basis, without stringent demands on the skill of the operator or the equipment. It is also desirable that the procedure would significantly decrease thenumber of shear strength tests required for verification. The experimental workinvolved tests in order to obtain an understanding of the behaviour of the sheets during clinching. The most notable observation concerned the stage of the process in which the upper sheet was initially bent, after which the deformation mechanism changed to shearing and elongation. The amount of deformation was measured relative to the original location of the upper sheet, and characterised as the C-measure. By understanding in detail the behaviour of the upper sheet, it waspossible to estimate a bending line function for the surface of the upper sheet. A procedure was developed, which makes it possible to estimate the process parameter X for each tool combination with a fixed die. The procedure is based on equating the volume of material on the punch side with the volume of the die. Detailed information concerning the behaviour of material on the punch side is required, assuming that the volume of die does not change during the process. The procedure was applied to shear strength testing of a sample material. The sample material was continuously hot-dip zinc-coated high-strength constructional steel,with a nominal thickness of 1.0 mm. The minimum Rp0.2 proof stress was 637 N/mm2. Such material has not yet been used extensively in light-frame housing, and little has been published on clinching of the material. The performance of the material is therefore of particular interest. Companies that use clinching on a daily basis stand to gain the greatest benefit from the procedure. By understanding the behaviour of sheets in different cases, it is possible to use data at an early stage for adjusting and optimising the process. In particular, the functionality of common tools can be increased since it is possible to characterise the complete range of existing tools. The study increases and broadens the amount ofbasic information concerning the clinching process. New approaches and points of view are presented and used for generating new knowledge.

Electromagnetic Design of a Solid Steel Rotor Motor for Demanding Operation Environments

The solid-rotor induction motor provides a mechanically and thermally reliable solution for demanding environments where other rotor solutions are prohibited or questionable. Solid rotors, which are manufactured of single pieces of ferromagnetic material, are commonly used in motors in which the rotationspeeds exceed substantially the conventional speeds of laminated rotors with squirrel-cage. During the operation of a solid-rotor electrical machine, the rotor core forms a conductor for both the magnetic flux and the electrical current. This causes an increase in the rotor resistance and rotor leakage inductance, which essentially decreases the power factor and the efficiency of the machine. The electromagnetic problems related to the solid-rotor induction motor are mostly associated with the low performance of the rotor. Therefore, the main emphasis in this thesis is put on the solid steel rotor designs. The rotor designs studied in thisthesis are based on the fact that the rotor construction should be extremely robust and reliable to withstand the high mechanical stresses caused by the rotational velocity of the rotor. In addition, the demanding operation environment sets requirements for the applied materials because of the high temperatures and oxidizing acids, which may be present in the cooling fluid. Therefore, the solid rotors analyzed in this thesis are made of a single piece of ferromagnetic material without any additional parts, such as copper end-rings or a squirrel-cage. A pure solid rotor construction is rigid and able to keep its balance over a large speed range. It also may tolerate other environmental stresses such as corroding substances or abrasive particles. In this thesis, the main target is to improve the performance of an induction motor equipped with a solid steel rotor by traditional methods: by axial slitting of the rotor, by selecting a proper rotor core material and by coating the rotor with a high-resistive stainless ferromagnetic material. In the solid steel rotor calculation, the rotor end-effects have a significant effect on the rotor characteristics. Thus, the emphasis is also put on the comparison of different rotor endfactors. In addition, a corrective slip-dependent end-factor is proposed. The rotor designs covered in this thesis are the smooth solid rotor, the axially slitted solid rotor and the slitted rotor having a uniform ferromagnetic coating cylinder. The thesis aims at design rules for multi-megawatt machines. Typically, mega-watt-size solidrotor machines find their applications mainly in the field of electric-motor-gas-compression systems, in steam-turbine applications, and in various types of largepower pump applications, where high operational speeds are required. In this thesis, a 120 kW, 10 000 rpm solid-rotor induction motor is usedas a small-scale model for such megawatt-range solid-rotor machines. The performance of the 120 kW solid-rotor induction motors is determined by experimental measurements and finite element calculations.

Direct-on-line axial flux permanent magnet synchronous generator static and dynamic performance

In distributed energy production, permanent magnet synchronous generators (PMSG) are often connected to the grid via frequency converters, such as voltage source line converters. The price of the converter may constitute a large part of the costs of a generating set. Some of the permanent magnet synchronous generators with converters and traditional separately excited synchronous generators couldbe replaced by direct-on-line (DOL) non-controlled PMSGs. Small directly networkconnected generators are likely to have large markets in the area of distributed electric energy generation. Typical prime movers could be windmills, watermills and internal combustion engines. DOL PMSGs could also be applied in island networks, such as ships and oil platforms. Also various back-up power generating systems could be carried out with DOL PMSGs. The benefits would be a lower priceof the generating set and the robustness and easy use of the system. The performance of DOL PMSGs is analyzed. The electricity distribution companies have regulations that constrain the design of the generators being connected to the grid. The general guidelines and recommendations are applied in the analysis. By analyzing the results produced by the simulation model for the permanent magnet machine, the guidelines for efficient damper winding parameters for DOL PMSGs are presented. The simulation model is used to simulate grid connections and load transients. The damper winding parameters are calculated by the finite element method (FEM) and determined from experimental measurements. Three-dimensional finite element analysis (3D FEA) is carried out. The results from the simulation model and 3D FEA are compared with practical measurements from two prototype axial flux permanent magnet generators provided with damper windings. The dimensioning of the damper winding parameters is case specific. The damper winding should be dimensioned based on the moment of inertia of the generating set. It is shown that the damper winding has optimal values to reach synchronous operation in the shortest period of time after transient operation. With optimal dimensioning, interferenceon the grid is minimized.

Automated formulation of optimisation models for steel beam structures

Over 70% of the total costs of an end product are consequences of decisions that are made during the design process. A search for optimal cross-sections will often have only a marginal effect on the amount of material used if the geometry of a structure is fixed and if the cross-sectional characteristics of its elements are property designed by conventional methods. In recent years, optimalgeometry has become a central area of research in the automated design of structures. It is generally accepted that no single optimisation algorithm is suitable for all engineering design problems. An appropriate algorithm, therefore, mustbe selected individually for each optimisation situation. Modelling is the mosttime consuming phase in the optimisation of steel and metal structures. In thisresearch, the goal was to develop a method and computer program, which reduces the modelling and optimisation time for structural design. The program needed anoptimisation algorithm that is suitable for various engineering design problems. Because Finite Element modelling is commonly used in the design of steel and metal structures, the interaction between a finite element tool and optimisation tool needed a practical solution. The developed method and computer programs were tested with standard optimisation tests and practical design optimisation cases. Three generations of computer programs are developed. The programs combine anoptimisation problem modelling tool and FE-modelling program using three alternate methdos. The modelling and optimisation was demonstrated in the design of a new boom construction and steel structures of flat and ridge roofs. This thesis demonstrates that the most time consuming modelling time is significantly reduced. Modelling errors are reduced and the results are more reliable. A new selection rule for the evolution algorithm, which eliminates the need for constraint weight factors is tested with optimisation cases of the steel structures that include hundreds of constraints. It is seen that the tested algorithm can be used nearly as a black box without parameter settings and penalty factors of the constraints.

Theoretical Analysis and Simulations of Vertically Vibrated Granular Materials

The dynamical properties ofshaken granular materials are important in many industrial applications where the shaking is used to mix, segregate and transport them. In this work asystematic, large scale simulation study has been performed to investigate the rheology of dense granular media, in the presence of gas, in a three dimensional vertical cylinder filled with glass balls. The base wall of the cylinder is subjected to sinusoidal oscillation in the vertical direction. The viscoelastic behavior of glass balls during a collision, have been studied experimentally using a modified Newton's Cradle device. By analyzing the results of the measurements, using numerical model based on finite element method, the viscous damping coefficient was determinedfor the glass balls. To obtain detailed information about the interparticle interactions in a shaker, a simplified model for collision between particles of a granular material was proposed. In order to simulate the flow of surrounding gas, a formulation of the equations for fluid flow in a porous medium including particle forces was proposed. These equations are solved with Large Eddy Simulation (LES) technique using a subgrid-model originally proposed for compressible turbulent flows. For a pentagonal prism-shaped container under vertical vibrations, the results show that oscillon type structures were formed. Oscillons are highly localized particle-like excitations of the granular layer. This self-sustaining state was named by analogy with its closest large-scale analogy, the soliton, which was first documented by J.S. Russell in 1834. The results which has been reportedbyBordbar and Zamankhan(2005b)also show that slightly revised fluctuation-dissipation theorem might apply to shaken sand, which appears to be asystem far from equilibrium and could exhibit strong spatial and temporal variations in quantities such as density and local particle velocity. In this light, hydrodynamic type continuum equations were presented for describing the deformation and flow of dense gas-particle mixtures. The constitutive equation used for the stress tensor provides an effective viscosity with a liquid-like character at low shear rates and a gaseous-like behavior at high shear rates. The numerical solutions were obtained for the aforementioned hydrodynamic equations for predicting the flow dynamics ofdense mixture of gas and particles in vertical cylindrical containers. For a heptagonal prism shaped container under vertical vibrations, the model results were found to predict bubbling behavior analogous to those observed experimentally. This bubbling behavior may be explained by the unusual gas pressure distribution found in the bed. In addition, oscillon type structures were found to be formed using a vertically vibrated, pentagonal prism shaped container in agreement with computer simulation results. These observations suggest that the pressure distribution plays a key rolein deformation and flow of dense mixtures of gas and particles under vertical vibrations. The present models provide greater insight toward the explanation of poorly understood hydrodynamic phenomena in the field of granular flows and dense gas-particle mixtures. The models can be generalized to investigate the granular material-container wall interactions which would be an issue of high interests in the industrial applications. By following this approach ideal processing conditions and powder transport can be created in industrial systems.

Theoretical and numerical study of thermomagnetic convection in magnetic fluids

In this thesis, the magnetic field control of convection instabilities and heat and mass transfer processesin magnetic fluids have been investigated by numerical simulations and theoretical considerations. Simulation models based on finite element and finite volume methods have been developed. In addition to standard conservation equations, themagnetic field inside the simulation domain is calculated from Maxwell equations and the necessary terms to take into account for the magnetic body force and magnetic dissipation have been added to the equations governing the fluid motion.Numerical simulations of magnetic fluid convection near the threshold supportedexperimental observations qualitatively. Near the onset of convection the competitive action of thermal and concentration density gradients leads to mostly spatiotemporally chaotic convection with oscillatory and travelling wave regimes, previously observed in binary mixtures and nematic liquid crystals. In many applications of magnetic fluids, the heat and mass transfer processes including the effects of external magnetic fields are of great importance. In addition to magnetic fluids, the concepts and the simulation models used in this study may be applied also to the studies of convective instabilities in ordinary fluids as well as in other binary mixtures and complex fluids.

Aortic blood pressure measured via EIT: investigation of different measurement settings.

Electrical impedance tomography (EIT) allows the measurement of intra-thoracic impedance changes related to cardiovascular activity. As a safe and low-cost imaging modality, EIT is an appealing candidate for non-invasive and continuous haemodynamic monitoring. EIT has recently been shown to allow the assessment of aortic blood pressure via the estimation of the aortic pulse arrival time (PAT). However, finding the aortic signal within EIT image sequences is a challenging task: the signal has a small amplitude and is difficult to locate due to the small size of the aorta and the inherent low spatial resolution of EIT. In order to most reliably detect the aortic signal, our objective was to understand the effect of EIT measurement settings (electrode belt placement, reconstruction algorithm). This paper investigates the influence of three transversal belt placements and two commonly-used difference reconstruction algorithms (Gauss-Newton and GREIT) on the measurement of aortic signals in view of aortic blood pressure estimation via EIT. A magnetic resonance imaging based three-dimensional finite element model of the haemodynamic bio-impedance properties of the human thorax was created. Two simulation experiments were performed with the aim to (1) evaluate the timing error in aortic PAT estimation and (2) quantify the strength of the aortic signal in each pixel of the EIT image sequences. Both experiments reveal better performance for images reconstructed with Gauss-Newton (with a noise figure of 0.5 or above) and a belt placement at the height of the heart or higher. According to the noise-free scenarios simulated, the uncertainty in the analysis of the aortic EIT signal is expected to induce blood pressure errors of at least ± 1.4 mmHg.

Virtuaaliprototyypin käyttö roottorin dynamiikan analysoinnissa

Työn tavoitteena oli selvittää kaupallisen dynamiikansimulointiohjelmiston soveltuvuus roottoridynamiikan analysointiin. Työssä keskityttiin erityisesti roottorin dynamiikkaan vaikuttavien epäideaalisuuksien mallintamiseen. Simulointitulosten tarkkuutta selvitettiin mittauksilla. Lisäksi vertailtiin yleiskäyttöisen dynamiikan simulointiohjelmiston ja roottoridynamiikan erikoisohjelmiston teoriaa. Tutkittava roottori oli paperikoneen putkitela. Telan joustavuus kuvattiin elementtimenetelmällä ratkaistujen moodien avulla. Elementtimallissa huomioitiin telan vaipan seinämänpaksuusvaihtelu, joka vaikuttaa telan massa- ja jäykkyysjakaumaan. Dynamiikkaohjelmistossa mallinnettiin telan tuennasta tulevat herätteet. Dynamiikkaohjelmistona käytettiin ADAMS:ia ja FEM-ohjelmana ANSYS:stä. Tuloksista havaittiin käytetyn menetelmän soveltuvan roottoridynamiikan analysointiin ja roottorin epäideaalisuuksien mallintamiseen. Simulointimallilla saatiin esille murtolukukriittiset pyörimisnopeudet ja telan kriittinen pyörimisnopeus vastasi hyvin mittaustuloksia.

Modelling the curing process in magneto-sensitive polymers: Rate-dependence and shrinkage

This paper deals with a phenomenologically motivated magneto-viscoelastic coupled finite strain framework for simulating the curing process of polymers under the application of a coupled magneto-mechanical road. Magneto-sensitive polymers are prepared by mixing micron-sized ferromagnetic particles in uncured polymers. Application of a magnetic field during the curing process causes the particles to align and form chain-like structures lending an overall anisotropy to the material. The polymer curing is a viscoelastic complex process where a transformation from fluid. to solid occurs in the course of time. During curing, volume shrinkage also occurs due to the packing of polymer chains by chemical reactions. Such reactions impart a continuous change of magneto-mechanical properties that can be modelled by an appropriate constitutive relation where the temporal evolution of material parameters is considered. To model the shrinkage during curing, a magnetic-induction-dependent approach is proposed which is based on a multiplicative decomposition of the deformation gradient into a mechanical and a magnetic-induction-dependent volume shrinkage part. The proposed model obeys the relevant laws of thermodynamics. Numerical examples, based on a generalised Mooney-Rivlin energy function, are presented to demonstrate the model capacity in the case of a magneto-viscoelastically coupled load.

Reluktanssiverkkomalli kestomagneettitahtikoneen staattisen ilmavälivuontiheyden mallinnukseen

Kaikissa pyörivissä sähkömoottoreissa vääntömomentin tuoton kannalta olennainen magneettivuo kulkee staattorin ja roottorin välillä ilmavälin kautta. Ilmaväli mallinnetaan koneensuunnittelun yhteydessä tämän vuoksi tarkasti. Elementtimenetelmällä voidaan analysoida moottoreita varsin tarkasti, mutta menetelmän käyttö vie paljon aikaa ja sovittaminen muihin laskentaympäristöihin on usein hankalaa. Tämän vuoksi voidaan käyttää riittävän tarkkuuden omaavia analyyttisiä laskentamenetelmiä, joiden sovittaminen muihin ohjelmaympäristöihin on helpompaa kuin elementtimenetelmää käytettäessä. Diplomityössä kehitetään reluktanssiverkkomalli kestomagneettien aikaansaaman ilmavälivuontiheyden mallintamiseen kestomagneettitahtikoneille, joissa on pinta-asennetut kestomagneetit. Kehitetyn reluktanssiverkkomallin toimivuutta vertaillaan muihin ilmavälivuontiheyden laskentamenetelmiin.

Hydrostaattisen kuormituselementin parametrien vaikutus tehontarpeeseen

Työssä tutkitaan paperiteollisuuden käyttämän taipumakompensoidun telan hydrostaattisen kuormituselementin parametrien vaikutusta tehontarpeeseen ja kuormituselementin käyttäytymistä eri parametreilla. Työssä käsitellään hydrostaattisen kuormituselementin FE-laskennassa käytetyn mallinosien rakentaminen, kokoonpano ja teoria. Lisäksi analyyttinen laskenta esitetään lyhyesti. FE-laskennassa käytetään virtaus-rakenne -vuorovaikutusanalyysiä. Laskenta suoritetaan ANSYS/Flotran ohjelmistolla (versio 5.5). Hydrostaattisen kuormituselementin tehontarpeen vähentämiseen löytyi kolme tekijää. Suurimpana yksittäisenä tekijänä on elementin kitkapinta-alan pienentäminen, jolla saavutetaan 44 % kitkatehontarpeen pieneneminen. Seuraavaksi eniten vaikuttaa kitkatehontarpeeseen kuormituselementin muotoilu, jolla saavutetaan 20 % kitkatehontarpeen pieneneminen. Kolmantena tekijänä on viskositeetin pienentäminen, jolla suurimmillaan saavutetaan 15 % kitkatehontarpeen pieneneminen.

Materiaalimallin määrityksestä kylmämuovattujen putkipalkkien K-liitosten FE-analyysissä

Työssä on tutkittu kylmämuovattujen nelikulmaisten putkipalkkien K-liitosten mallinnusta epälineaarisella elementtimenetelmällä. Työn tärkeimpänä tavoitteena on ollut kehittää putkipalkin osien materiaalimalleja siten, että liitosten kestävyyttä voidaan tutkia laboratoriokokeiden ohella luotettavasti myös elementtimenetelmällä. Toisena tavoitteena on ollut tutkia, voidaanko putkipalkkien liitosten mitoitusohjeita turvallisesti soveltaa kylmämuovatuille putkipalkeille, joissa valmistusprosessi aiheuttaa muutoksia materiaaliominaisuuksiin, erityisesti muodonmuutoskykyyn. Työssä tehtyjen laboratoriokokeiden ja elementtianalyysien perusteella elementti-menetelmä on käyttökelpoinen työkalu putkipalkkiliitosten staattista kestävyyttä määritettäessä, kun materiaalimallit on määritetty oikein. Erityisesti liitoksen käyttö-rajatilan mukaisen kestävyyden laskennassa elementtimenetelmällä saadaan hyvin laboratoriokokeita vastaavia tuloksia. Tehdyt laboratoriokokeet osoittavat myös, että Eurocode 3:n mukaisia putkipalkkien liitosten mitoitusohjeita voi turvallisesti käyttää kylmämuovatuille putkipalkeille.

Tuulikuormien vaikutus teleskooppimaston toimintaan

Tämä diplomityö on tehty Patria Vehicles Oy:n toimeksiannosta. Patria Vehicles Oy:n tuotantoon kuuluvat vaativiin maasto-olosuhteisiin soveltuvat sotilasajoneuvot sekä teleskooppimastot. Tutkimuksen tarkoituksena oli mallintaa mastoperävaunusta joustava malli, johon vaikuttavat tuulikuormat. Mallin avulla voidaan tutkia maston siirtymiä, kallistumia sekä kiertymiä. Tutkimuksessa on käytetty ADAMS-simulointiohjelmistoa sekä I-DEAS- FEM ohjelmistoa. Dynaamisten ongelmien ratkaisemiseksi on ymmärrettävä rakenteiden käyttäytymistä. Tuulikuormien mallintamisen edellytyksenä on tuulikuormien syntymisen ymmärtäminen. Tämän työn peruslähtökohtana on mallintaa kaikki maston jäykkyyteen vaikuttavat komponentit joustavina FE-menetelmän avulla. Luodaan superelementit Craig-Bamptonin ominaismuotojen superponointimenetelmällä. Nämä superelementit liitetään toisiinsa ja asetetaan niille tuulikuormat. Luodaan kosketukset puomien, sekä maan ja maston välille. Pienennetään joustavien osien ominaismuotojen määrää, jotta saataisiin nopeammat analyysit. Parametrisoidaan malli, jolloin voidaan analysoida mallilla useampia tapauksia. Verifioidaan malli varmistaaksemme sen oikeellisuuden. Taulukoidaan tulokset.