Design Aspects of Megawatt-Range Direct-Driven Permanent MagnetWind Generators

Direct-driven permanent magnet synchronous generator is one of the most promising topologies for megawatt-range wind power applications. The rotational speed of the direct-driven generator is very low compared with the traditional electrical machines. The low rotational speed requires high torque to produce megawatt-range power. The special features of the direct-driven generators caused by the low speed and high torque are discussed in this doctoral thesis. Low speed and high torque set high demands on the torque quality. The cogging torque and the load torque ripple must be as low as possible to prevent mechanical failures. In this doctoral thesis, various methods to improve the torque quality are compared with each other. The rotor surface shaping, magnet skew, magnet shaping, and the asymmetrical placement of magnets and stator slots are studied not only by means of torque quality, but also the effects on the electromagnetic performance and manufacturability of the machine are discussed. The heat transfer of the direct-driven generator must be designed to handle the copper losses of the stator winding carrying high current density and to keep the temperature of the magnets low enough. The cooling system of the direct-driven generator applying the doubly radial air cooling with numerous radial cooling ducts was modeled with a lumped-parameter-based thermal network. The performance of the cooling system was discussed during the steady and transient states. The effect of the number and width of radial cooling ducts was explored. The large number of radial cooling ducts drastically increases the impact of the stack end area effects, because the stator stack consists of numerous substacks. The effects of the radial cooling ducts on the effective axial length of the machine were studied by analyzing the crosssection of the machine in the axial direction. The method to compensate the magnet end area leakage was considered. The effect of the cooling ducts and the stack end area effects on the no-load voltages and inductances of the machine were explored by using numerical analysis tools based on the three-dimensional finite element method. The electrical efficiency of the permanent magnet machine with different control methods was estimated analytically over the whole speed and torque range. The electrical efficiencies achieved with the most common control methods were compared with each other. The stator voltage increase caused by the armature reaction was analyzed. The effect of inductance saturation as a function of load current was implemented to the analytical efficiency calculation.

Pienahitsin sijainnin vaikutus kuormaa kantamattoman T-liitoksen väsymiskestävyyteen

Hitsatuissa rakenteissa on paljon hitsien aiheuttamia epäjatkuvuuskohtia, jotka ovat alttiita vaurioitumaan väsyttävän kuormituksen alaisina. Standardit ja normit auttavat suunnittelua rakenneyksityiskohtien analysoinnissa, mutta ohjeiden arvot vaihtelevat paljon. Onkin tär-keää päättää, mitä standardia tai normia käytetään ja saada tutkimustuloksia tukemaan ohjeen käyttöä. Tämä työ on tehty Konecranes Finland Oy:lle. Työn tarkoituksena on etsiä käytössä olevi-en ohjeiden joukosta suunnittelun kannalta parhaiten todellista rakennetta kuvaava. Työssä esitellään käytössä olevia standardeja ja normeja sekä vertaillaan niiden välisiä eroja rakenneyksityiskohdan väsymiskestävyyden määrittämisessä. Työssä tarkastellaan myös hitsatun rakenteen väsymiseen vaikuttavia tekijöitä sekä väsymisanalyysissä käytet-täviä menetelmiä. Tutkittavaan liitokseen valittiin kolme erilaista hitsiä, yksi –ja kaksipuolinen pienahitsi sekä katkopienahitsi. Elementtimenetelmän avulla luotiin malleja koekappaleista ja laskettiin niille väsymiskestävyydet. Tämä lisäksi väsymistä tarkasteltiin väsytyskokeilla, jotka suoritti Lappeenrannan Teknillisen Yliopiston Teräsrakenteiden laboratorio. Lopputuloksena työstä saatiin tietoa hitsien väsymiskestävyyksistä. Liitosten väsymiskes-tävyydet olivat hyviä ja varsinkin katkopienahitsin tulokset olivat odotettua paremmat. Kat-kopienahitsin väsymiskestävyys vaatii tämän työn tulosten perusteella lisää tutkimista.

High-speed solid-rotor induction machine – calculation program

The aim of this thesis is to utilize the technology developed at LUT and to provide an easy tool for high-speed solid-rotor induction machine preliminary design. Computer aided design tool MathCAD has been chosen as the environment for realizing the calculation program. Four versions of the design program have been made depending on the motor rotor type. The first rotor type is an axially slitted solid-rotor with steel end rings. The next one is an axially slitted solid-rotor with copper end rings. The third machine type is a solid rotor with deep, rectangular copper bars and end rings (squirrel cage). And the last one is a solid-rotor with round copper bars and end rings (squirrel cage). Each type of rotor has its own specialties but a general thread of design is common. This paper follows the structure of the calculating program and explains some features and formulas. The attention is concentrated on the difference between laminated and solid-rotor machine design principles. There is no deep analysis of the calculation ways are presented. References for all solution methods appearing during the design procedure are given for more detailed studying. This thesis pays respect to the latest innovations in solid-rotor machines theory. Rotor ends’ analytical calculation follows the latest knowledge in this field. Correction factor for adjusting the rotor impedance is implemented. The purpose of the created design program is to calculate the preliminary dimensions of the machine according to initial data. Obtained results are not recommended for exact machine development. Further more detailed design should be done in a finite element method application. Hence, this thesis is a practical tool for the prior evaluating of the high-speed machine with different solid-rotor types parameters.

Soodakattiloiden ilmakanavien hyödyntäminen jäykistävänä rakenteena

Työssä tutkittiin soodakattiloiden ilmakanavien hyödyntämistä jäykistävänä rakenteena. Työssä käsiteltiin yksittäisiä jäykistämättömiä ja jäykistettyjä levykenttiä ja niiden lommahduskestävyyttä Eurokoodi standardin mukaisesti ja elementtimenetelmän avulla. Lisäksi käsiteltiin lommahduksen teoriaa ja levykenttien käyttäytymistä yleisellä tasolla erilaisilla kuormituksilla ja reunaehdoilla. Työn tavoitteena oli selvittää kuinka lommahdus tutkitaan Eurokoodin mukaisesti ja elementtimenetelmää hyödyntäen, kun levykentän kuormituksena on poikittainen kuormitus tason suuntaisen kuormituksen lisäksi. Työssä tutkittiin kahden eri elementtimenetelmään pohjautuvan ratkaisuvaihtoehdon käyttöä lommahduslaskennassa. Työssä kehitettiin Eurokoodin sovellettu yhteisvaikutuskaavan käyttö lineaarisen ominaisarvotehtävän ratkaisun lisänä, jossa otetaan huomioon painekuorman vaikutus levykentän lommahduksessa. Kehitettyä menetelmää sovellettiin ilmakanavan esimerkkirakenteen mitoituksessa.

Eri laskentamenetelmien vertailu hitsatun rakenteen väsymislaskennassa

Tässä työssä tutkittiin eri mitoitusmenetelmien soveltuvuutta hitsattujen rakenteiden vä-symislaskennassa. Käytetyt menetelmät olivat rakenteellinen jännityksen menetelmä, te-hollisen lovijännityksen menetelmä ja murtumismekaniikka. Lisäksi rakenteellisen jänni-tyksen määrittämiseksi käytettiin kolmea eri menetelmää. Menetelmät olivat pintaa pitkin ekstrapolointi, paksuuden yli linearisointi ja Dongin menetelmä. Väsymiskestävyys määritettiin kahdelle hitsiliitoksen yksityiskohdalle. Laskenta tehtiin käyttäen elementtimenetelmää rakenteen 3D-mallille. Tutkittavasta aggregaattirungosta oli olemassa FE-malli mutta alimallinnustekniikkaa hyödyntämällä pystyttiin yksityiskohtai-semmin tutkimaan vain pientä osaa koko rungon mallista. Rakenteellisen jännityksen menetelmä perustuu nimellisiin jännityksiin. Kyseinen mene-telmä ei vaadi geometrian muokkausta. Yleensä rakenteellisen jännityksen menetelmää käytetään hitsin rajaviivan väsymislaskennassa, mutta joissain tapauksissa sitä on käytetty juuren puolen laskennassa. Tässä työssä rakenteellisen jännityksen menetelmää käytettiin myös juuren puolen tutkimisessa. Tehollista lovijännitystä tutkitaan mallintamalla 1 mm fiktiiviset pyöristykset sekä rajaviivalle että juuren puolelle. Murtumismekaniikan so-veltuvuutta tutkittiin käyttämällä Franc2D särön kasvun simulointiohjelmaa. Väsymislaskennan tulokset eivät merkittävästi poikkea eri laskentamenetelmien välillä. Ainoastaan rakenteellisen jännityksen Dongin menetelmällä saadaan poikkeavia tuloksia. Tämä johtuu pääasiassa siitä, että menetelmän laskentaetäisyydestä ei ole tietoa. Raken-teellisen jännityksen menetelmällä, tehollisen lovijännityksen menetelmällä ja murtumis-mekaniikalla saadaan samansuuntaiset tulokset. Suurin ero menetelmien välillä on mal-linnuksen ja laskennan vaatima työmäärä.

Modeling of residual stresses and distortion due to welding in fillet welds

Welding has a growing role in modern world manufacturing. Welding joints are extensively used from pipes to aerospace industries. Prediction of welding residual stresses and distortions is necessary for accurate evaluation of fillet welds in relation to design and safety conditions. Residual stresses may be beneficial or detrimental, depending whether they are tensile or compressive and the loading. They directly affect the fatigue life of the weld by impacting crack growth rate. Beside theoretical background of residual stresses this study calculates residual stresses and deformations due to localized heating by welding process and subsequent rapid cooling in fillet welds. Validated methods are required for this purpose due to complexity of process, localized heating, temperature dependence of material properties and heat source. In this research both empirical and simulation methods were used for the analysis of welded joints. Finite element simulation has become a popular tool of prediction of welding residual stresses and distortion. Three different cases with and without preload have been modeled during this study. Thermal heat load set is used by calculating heat flux from the given heat input energy. First the linear and then nonlinear material behavior model is modeled for calculation of residual stresses. Experimental work is done to calculate the stresses empirically. The results from both the methods are compared to check their reliability. Residual stresses can have a significant effect on fatigue performance of the welded joints made of high strength steel. Both initial residual stress state and subsequent residual stress relaxation need to be considered for accurate description of fatigue behavior. Tensile residual stresses are detrimental and will reduce the fatigue life and compressive residual stresses will increase it. The residual stresses follow the yield strength of base or filler material and the components made of high strength steel are typically thin, where the role of distortion is emphasizing.

Capacity of an I-beam made of S960 QC based on Eurocode 3

The capacity of beams is a very important factor in the study of durability of structures and structural members. The capacity of a high-strength steel I-beam made of S960 QC was investigated in this study. The investigation included assessment of the service limits and ultimate limits of the steel beam. The thesis was done according to European standards for steel structures, Eurocode 3. An analytical method was used to determine the throat thickness, deformation, elastic and plastic moment capacities as well as the fatigue life of the beam. The results of the analytical method were compared with those obtained by Finite Element Analysis (FEA). Elastic moment capacity obtained by the analytical method was 172 kNm. FEA and the analytical method predicted the maximum lateral-torsional buckling (LTB) capacity in the range of 90-93 kNm and the probability of failure as a result of LTB is estimated to be 50%. The lateral buckling capacity meant that the I-beam can carry a safe load of 300 kN instead of the initial load of 600 kN. The beam is liable to fail shortly after exceeding the elastic moment capacity. Based on results in of the different approaches, it was noted that FEA predicted higher deformation values on the load-deformation curve than the analytical results. However, both FEA and the analytical methods predicted identical results for nominal stress range and moment capacities. Fatigue life was estimated to be in the range of 53000-64000 cycles for bending stress range using crack propagation equation and strength-life approach. As Eurocode 3 is limited to steel grades up to S690, results for S960 must be verified with experimental data and appropriate design rules.

BEM/FEM non-linear model applied to transient analysis with viscous damping

In this article a two-dimensional transient boundary element formulation based on the mass matrix approach is discussed. The implicit formulation of the method to deal with elastoplastic analysis is considered, as well as the way to deal with viscous damping effects. The time integration processes are based on the Newmark rhoand Houbolt methods, while the domain integrals for mass, elastoplastic and damping effects are carried out by the well known cell approximation technique. The boundary element algebraic relations are also coupled with finite element frame relations to solve stiffened domains. Some examples to illustrate the accuracy and efficiency of the proposed formulation are also presented.

On the numerical simulation of machining processes

Numerical simulation of machining processes can be traced back to the early seventies when finite element models for continuous chip formation were proposed. The advent of fast computers and development of new techniques to model large plastic deformations have favoured machining simulation. Relevant aspects of finite element simulation of machining processes are discussed in this paper, such as solution methods, material models, thermo-mechanical coupling, friction models, chip separation and breakage strategies and meshing/re-meshing strategies.

An Hp-Adaptive Hierarchical Formulation for the Boundary Element Method Applied to Elasticity in Two Dimensions

This paper presents an HP-Adaptive Procedure with Hierarchical formulation for the Boundary Element Method in 2-D Elasticity problems. Firstly, H, P and HP formulations are defined. Then, the hierarchical concept, which allows a substantial reduction in the dimension of equation system, is introduced. The error estimator used is based on the residual computation over each node inside an element. Finally, the HP strategy is defined and applied to two examples.

Voimalaitoksen kanavien kannakkeiden rakennesuunnittelu

Tässä diplomityössä esitetään voimalaitoksen kanavien kannakkeiden rakennesuunnittelussa tarvittavat laskentamenetelmät. Työssä rakenteiden suunnitteluun ja mitoitukseen käytetään pääasiassa Eurokoodi 3 teräsrakenteiden suunnittelustandardin mukaista rajatilamitoitusta. Lisäksi kehitetään mitoitustyökaluja tärkeimpien kanavakannakkeiden suunnitteluun. Toteutettujen mitoitustyökalujen toiminta verifioidaan lujuusopin elementtimenetelmällä tehtävin tarkistuslaskelmin. Laskentatyökalujen analyyttisen ratkaisun verifioitiin olevan varmalla puolella kaikissa tutkituissa ilmiöissä. Työssä verifioituja menetelmiä voidaan soveltaa myös muiden vastaavien rakenteiden mitoittamiseen. Työssä luotujen laskentatyökalujen sisältämät laskentamenetelmät mahdollistavat monenlaisten rakenteiden vaatimustenmukaisen suunnittelun.

Modeling and Parameter Estimation of Double-Star Permanent Magnet Synchronous Machines

The power rating of wind turbines is constantly increasing; however, keeping the voltage rating at the low-voltage level results in high kilo-ampere currents. An alternative for increasing the power levels without raising the voltage level is provided by multiphase machines. Multiphase machines are used for instance in ship propulsion systems, aerospace applications, electric vehicles, and in other high-power applications including wind energy conversion systems. A machine model in an appropriate reference frame is required in order to design an efficient control for the electric drive. Modeling of multiphase machines poses a challenge because of the mutual couplings between the phases. Mutual couplings degrade the drive performance unless they are properly considered. In certain multiphase machines there is also a problem of high current harmonics, which are easily generated because of the small current path impedance of the harmonic components. However, multiphase machines provide special characteristics compared with the three-phase counterparts: Multiphase machines have a better fault tolerance, and are thus more robust. In addition, the controlled power can be divided among more inverter legs by increasing the number of phases. Moreover, the torque pulsation can be decreased and the harmonic frequency of the torque ripple increased by an appropriate multiphase configuration. By increasing the number of phases it is also possible to obtain more torque per RMS ampere for the same volume, and thus, increase the power density. In this doctoral thesis, a decoupled d–q model of double-star permanent-magnet (PM) synchronous machines is derived based on the inductance matrix diagonalization. The double-star machine is a special type of multiphase machines. Its armature consists of two three-phase winding sets, which are commonly displaced by 30 electrical degrees. In this study, the displacement angle between the sets is considered a parameter. The diagonalization of the inductance matrix results in a simplified model structure, in which the mutual couplings between the reference frames are eliminated. Moreover, the current harmonics are mapped into a reference frame, in which they can be easily controlled. The work also presents methods to determine the machine inductances by a finite-element analysis and by voltage-source inverters on-site. The derived model is validated by experimental results obtained with an example double-star interior PM (IPM) synchronous machine having the sets displaced by 30 electrical degrees. The derived transformation, and consequently, the decoupled d–q machine model, are shown to model the behavior of an actual machine with an acceptable accuracy. Thus, the proposed model is suitable to be used for the model-based control design of electric drives consisting of double-star IPM synchronous machines.

Vääntökuormitetun koteloprofiilipuomin väsymiskestävyys

Tässä työssä tutkittiin vääntökuormitettua hitsattua koteloprofiilipuomia. Vääntökuormitus aiheuttaa koteloprofiiliin vinouttavan voimasysteemin, joka aiheuttaa kotelopalkkiin sekä poikittaisia taivutusjännityksiä että pitkittäisiä jännityksiä. Vinoutumisen aiheuttamia lisära-situksia on tutkittu analyyttisesti BEF-analogian avulla sekä elementtimenetelmää apuna käyttäen. Lisäksi vääntökuormitus aiheuttaa puomiin estetyn väännön jännityksiä. Väsymiskestoikää tutkittiin laboratoriossa suoritettujen väsytyskokeiden avulla sekä las-kennallisesti. Kestoikälaskennassa käytettiin tehollisen lovijännityksen menetelmää hitsin juuren puolen väsymisen arvioinnissa sekä hot spot- jännityksen menetelmää hitsin rajavii-van väsymisen arvioinnissa. Teholliset lovijännitykset sekä hot spot jännitykset ovat määri-tetty elementtimenetelmän avulla. Laboratoriokokeiden ja elementtimenetelmä laskennan perusteella saatiin rakenteen kes-toiän kannalta kriittiset detaljit määritettyä. Kriittiset detaljit sijaitsevat puomin päädyn rakenteissa. Tutkimuksessa saatiin selville, että kriittisten detaljien rakenteellisilla ja valmis-tusteknisillä ratkaisuilla on merkittävä vaikutus lopullisen tuotteen väsymiskestoikään.

Modal testing and numerical modeling of the dynamic properties of layered sheet-steel structure

Recently, due to the increasing total construction and transportation cost and difficulties associated with handling massive structural components or assemblies, there has been increasing financial pressure to reduce structural weight. Furthermore, advances in material technology coupled with continuing advances in design tools and techniques have encouraged engineers to vary and combine materials, offering new opportunities to reduce the weight of mechanical structures. These new lower mass systems, however, are more susceptible to inherent imbalances, a weakness that can result in higher shock and harmonic resonances which leads to poor structural dynamic performances. The objective of this thesis is the modeling of layered sheet steel elements, to accurately predict dynamic performance. During the development of the layered sheet steel model, the numerical modeling approach, the Finite Element Analysis and the Experimental Modal Analysis are applied in building a modal model of the layered sheet steel elements. Furthermore, in view of getting a better understanding of the dynamic behavior of layered sheet steel, several binding methods have been studied to understand and demonstrate how a binding method affects the dynamic behavior of layered sheet steel elements when compared to single homogeneous steel plate. Based on the developed layered sheet steel model, the dynamic behavior of a lightweight wheel structure to be used as the structure for the stator of an outer rotor Direct-Drive Permanent Magnet Synchronous Generator designed for high-power wind turbines is studied.

Vibrations in Rotating Machinery Arising from Minor Imperfections in Component Geometries

Vibrations in machines can cause noise, decrease the performance, or even damage the machine. Vibrations appear if there is a source of vibration that excites the system. In the worst case scenario, the excitation frequency coincides with the natural frequency of the machine causing resonance. Rotating machines are a machine type, where the excitation arises from the machine itself. The excitation originates from the mass imbalance in the rotating shaft, which always exists in machines that are manufactured using conventional methods. The excitation has a frequency that is dependent on the rotational speed of the machine. The rotating machines in industrial use are usually designed to rotate at a constant rotational speed, the case where the resonances can be easily avoided. However, the machines that have a varying operational speed are more problematic due to a wider range of frequencies that have to be avoided. Vibrations, which frequencies equal to rotational speed frequency of the machine are widely studied and considered in the typical machine design process. This study concentrates on vibrations, which arise from the excitations having frequencies that are multiples of the rotational speed frequency. These vibrations take place when there are two or more excitation components in a revolution of a rotating shaft. The dissertation introduces four studies where three kinds of machines are experiencing vibrations caused by different excitations. The first studied case is a directly driven permanent magnet generator used in a wind power plant. The electromagnetic properties of the generator cause harmonic excitations in the system. The dynamic responses of the generator are studied using the multibody dynamics formulation. In another study, the finite element method is used to study the vibrations of a magnetic gear due to excitations, which frequencies equal to the rotational speed frequency. The objective is to study the effects of manufacturing and assembling inaccuracies. Particularly, the eccentricity of the rotating part with respect to non-rotating part is studied since the eccentric operation causes a force component in the direction of the shortest air gap. The third machine type is a tube roll of a paper machine, which is studied while the tube roll is supported using two different structures. These cases are studied using different formulations. In the first case, the tube roll is supported by spherical roller bearings, which have some wavinesses on the rolling surfaces. Wavinesses cause excitations to the tube roll, which starts to resonate at the frequency that is a half of the first natural frequency. The frequency is in the range where the machine normally operates. The tube roll is modeled using the finite element method and the bearings are modeled as nonlinear forces between the tube roll and the pedestals. In the second case studied, the tube roll is supported by freely rotating discs, which wavinesses are also measured. The above described phenomenon is captured as well in this case, but the simulation methodology is based on the flexible multibody dynamics formulation. The simulation models that are used in both of the last two cases studied are verified by measuring the actual devices and comparing the simulated and measured results. The results show good agreement.