Assessment of meshing techniques for fatigue analysis of load carrying welded joints using the effective notch approach

The effective notch stress approach for the fatigue strength assessment of welded structures as included in the Fatigue Design Recommendation of the IIW requires the numerical analysis of the elastic notch stress in the weld toe and weld root which is fictitiously rounded with a radius of 1mm. The goal of this thesis work was to consider alternate meshing strategies when using the effective notch stress approach to assess the fatigue strength of load carrying partial penetration fillet-welded cruciform joints. In order to establish guidelines for modeling the joint and evaluating the results, various two-dimensional (2D) finite element analyses were carried out by systematically varying the thickness of the plates, the weld throat thickness, the degree of bending, and the shape and location of the modeled effective notch. To extend the scope of this work, studies were also carried out on the influence of

Application of capillary electrophoretic methods for determining carbohydrates and aliphatic carboxylic acids formed during wood processing

Knowledge of the behaviour of cellulose, hemicelluloses, and lignin during wood and pulp processing is essential for understanding and controlling the processes. Determination of monosaccharide composition gives information about the structural polysaccharide composition of wood material and helps when determining the quality of fibrous products. In addition, monitoring of the acidic degradation products gives information of the extent of degradation of lignin and polysaccharides. This work describes two capillary electrophoretic methods developed for the analysis of monosaccharides and for the determination of aliphatic carboxylic acids from alkaline oxidation solutions of lignin and wood. Capillary electrophoresis (CE), in its many variants is an alternative separation technique to chromatographic methods. In capillary zone electrophoresis (CZE) the fused silica capillary is filled with an electrolyte solution. An applied voltage generates a field across the capillary. The movement of the ions under electric field is based on the charge and hydrodynamic radius of ions. Carbohydrates contain hydroxyl groups that are ionised only in strongly alkaline conditions. After ionisation, the structures are suitable for electrophoretic analysis and identification through either indirect UV detection or electrochemical detection. The current work presents a new capillary zone electrophoretic method, relying on in-capillary reaction and direct UV detection at the wavelength of 270 nm. The method has been used for the simultaneous separation of neutral carbohydrates, including mono- and disaccharides and sugar alcohols. The in-capillary reaction produces negatively charged and UV-absorbing compounds. The optimised method was applied to real samples. The methodology is fast since no other sample preparation, except dilution, is required. A new method for aliphatic carboxylic acids in highly alkaline process liquids was developed. The goal was to develop a method for the simultaneous analysis of the dicarboxylic acids, hydroxy acids and volatile acids that are oxidation and degradation products of lignin and wood polysaccharides. The CZE method was applied to three process cases. First, the fate of lignin under alkaline oxidation conditions was monitored by determining the level of carboxylic acids from process solutions. In the second application, the degradation of spruce wood using alkaline and catalysed alkaline oxidation were compared by determining carboxylic acids from the process solutions. In addition, the effectiveness of membrane filtration and preparative liquid chromatography in the enrichment of hydroxy acids from black liquor was evaluated, by analysing the effluents with capillary electrophoresis.

Qualitative Characteristics and Quantitative Measures of Solution's Reliability in Discrete Optimization: Traditional Analytical Approaches, Innovative Computational Methods and Applicability

The purpose of this thesis is twofold. The first and major part is devoted to sensitivity analysis of various discrete optimization problems while the second part addresses methods applied for calculating measures of solution stability and solving multicriteria discrete optimization problems. Despite numerous approaches to stability analysis of discrete optimization problems two major directions can be single out: quantitative and qualitative. Qualitative sensitivity analysis is conducted for multicriteria discrete optimization problems with minisum, minimax and minimin partial criteria. The main results obtained here are necessary and sufficient conditions for different stability types of optimal solutions (or a set of optimal solutions) of the considered problems. Within the framework of quantitative direction various measures of solution stability are investigated. A formula for a quantitative characteristic called stability radius is obtained for the generalized equilibrium situation invariant to changes of game parameters in the case of the H¨older metric. Quality of the problem solution can also be described in terms of robustness analysis. In this work the concepts of accuracy and robustness tolerances are presented for a strategic game with a finite number of players where initial coefficients (costs) of linear payoff functions are subject to perturbations. Investigation of stability radius also aims to devise methods for its calculation. A new metaheuristic approach is derived for calculation of stability radius of an optimal solution to the shortest path problem. The main advantage of the developed method is that it can be potentially applicable for calculating stability radii of NP-hard problems. The last chapter of the thesis focuses on deriving innovative methods based on interactive optimization approach for solving multicriteria combinatorial optimization problems. The key idea of the proposed approach is to utilize a parameterized achievement scalarizing function for solution calculation and to direct interactive procedure by changing weighting coefficients of this function. In order to illustrate the introduced ideas a decision making process is simulated for three objective median location problem. The concepts, models, and ideas collected and analyzed in this thesis create a good and relevant grounds for developing more complicated and integrated models of postoptimal analysis and solving the most computationally challenging problems related to it.

Quantum transport phenomena and shallow impurity states in CdSb

This work is dedicated to investigation of the energy spectrum of one of the most anisotropic narrow-gap semiconductors, CdSb. At the beginning of the present studies even the model of its energy band structure was not clear. Measurements of galvanomagnetic effects in wide temperature range (1.6 - 300 K) and in magnetic fields up to 30 T were chosen for clarifying of the energy spectrum in the intentionally undoped CdSb single crystals and doped with shallow impurities (In, Ag). Detection of the Shubnikov - de Haas oscillations allowed estimating the fundamental energy spectrum parameters. The shapes of the Fermi surfaces of electrons (sphere) and holes (ellipsoid), the number of the equivalent extremums for valence band (2) and their positions in the Brillouin zone were determined for the first time in this work. Also anisotropy coefficients, components of the tensor of effective masses of carriers, effective masses of density of states, nonparabolicity of the conduction and valence bands, g-factor and its anisotropy for n- and p-CdSb were estimated for the first time during these studies. All the results obtained are compared with the cyclotron resonance data and the corresponding theoretical calculations for p-CdSb. This is basic information for the analyses of the complex transport properties of CdSb and for working out the energy spectrum model of the shallow energy levels of defects and impurities in this semiconductor. It was found out existence of different mechanisms of hopping conductivity in the presence of metal - insulator transition induced by magnetic field in n- and p-CdSb. Quite unusual feature opened in CdSb is that different types of hopping conductivity may take place in the same crystal depending on temperature, magnetic field or even orientation of crystal in magnetic field. Transport properties of undoped p-CdSb samples show that the anisotropy of the resistivity in weak and strong magnetic fields is determined completely by the anisotropy of the effective mass of the holes. Temperature and magnetic field dependence of the Hall coefficient and magnetoresistance is attributed to presence of two groups of holes with different concentrations and mobilities. The analysis demonstrates that below Tcr ~ 20 K and down to ~ 6 - 7 K the low-mobile carriers are itinerant holes with energy E2 ≈ 6 meV. The high-mobile carriers, at all temperatures T < Tcr, are holes activated thermally from a deeper acceptor band to itinerant states of a shallower acceptor band with energy E1 ≈ 3 meV. Analysis of temperature dependences of mobilities confirms the existence of the heavy-hole band or a non-equivalent maximum and two equivalent maxima of the light-hole valence band. Galvanomagnetic effects in n-CdSb reveal the existence of two groups of carriers. These are the electrons of a single minimum in isotropic conduction band and the itinerant electrons of the narrow impurity band, having at low temperatures the energies above the bottom of the conduction band. It is found that above this impurity band exists second impurity band of only localized states and the energy of both impurity bands depend on temperature so that they sink into the band gap when temperature is increased. The bands are splitted by the spin, and in strong magnetic fields the energy difference between them decreases and redistribution of the electrons between the two impurity bands takes place. Mobility of the conduction band carriers demonstrates that scattering in n-CdSb at low temperatures is strongly anisotropic. This is because of domination from scattering on the neutral impurity centers and increasing of the contribution to mobility from scattering by acoustic phonons when temperature increases. Metallic conductivity in zero or weak magnetic field is changed to activated conductivity with increasing of magnetic field. This exhibits a metal-insulator transition (MIT) induced by the magnetic field due to shift of the Fermi level from the interval of extended states to that of the localized states of the electron spectrum near the edge of the conduction band. The Mott variablerange hopping conductivity is observed in the low- and high-field intervals on the insulating side of the MIT. The results yield information about the density of states, the localization radius of the resonant impurity band with completely localized states and about the donor band. In high magnetic fields this band is separated from the conduction band and lies below the resonant impurity bands.

Galvanomagnetic effects in strongly doped p-Bi2Te3:Sn crystals

Magnetic field dependencies of Hall coefficient and magnetoresistivity are investigated in classical and quantizing magnetic fields in p-Bi2Te3 crystals heavily doped with Sn grown by Czochralsky method. Magnetic field was parallel to the trigonal axis C3. Shubnikov-de Haas effect and quantum oscillations of the Hall coefficient were measured at temperatures 4.2 K and 11 K. On the basis of the magnetic field dependence of the Hall coefficient a method of estimation of the Hall factor and Hall mobility using the Drabble- Wolf six ellipsoid model is proposed. Shubnikov-de Haas effect and quantum oscillations of the Hall coefficient were observed at 4.2 K and 11 K. New evidence for the existence of the narrow band of Sn impurity states was shown. This band is partly filled by electrons and it is overlapping with the valence states of the light holes. Parameters of the impurity states, their energy ESn - 15 meV, band broadening ¿<< k0T and localization radius of the impuritystate R - 30 Å were obtained.

Metsäpalveluyritysten menestyminen ja siihen vaikuttavat tekijät

Tutkimuksessa kartoitettiin metsäpalveluyritysten toimintaympäristöä sekä menestymistä ja siihen vaikuttavia tekijöitä. Tutkimusta varten haastateltiin 27 metsäpalveluyrittäjää kesäkuussa 1999. Heikommin menestyneiden ja Menestyneimpien metsäpalveluyritysten löytämiseksi luotiin kaksi menestyspistejärjestelmää: KASKARA (KASvu, KAnnattavuus & RAhoitus) -mittaristo ja Tasapainotettu (Balanced Scorecard) mittaristo. KASKARA-mittaristossa tarkastelun kohteena olivat yrityksen kasvu (yrityksen liikevaihdon kasvuprosentti ja investoinnit) ja kannattavuus sekä rahoituksellinen tila (maksuvalmius ja vakavaraisuus) yrittäjän itsensä arvioimana. Toisaalta käytettäessä Tasapainotettua mittaristoa KASKARA-mittaristoon liitettiin vielä neljä muuttujaa: yrityksen tärkeimmän asiakkaan antama palaute yrittäjälle, yrittäjän hankkima täydennyskoulutus, kovien kilpailijoiden lukumäärä yrityksen nykyisellä toiminta-alueella ja yrittäjän oma arvio yrityksensä paremmuudesta kilpailijoihin nähden. Tutkimuksen Menestyneimpiä yrityksiä voitiin kuvata seuraavasti: yrittäjä kertoi olleensa luottavainen yritystä perustaessaan; yhtenä tärkeimmistä motiiveista yrityksensä perustamiseen yrittäjä piti pyrkimystä ansiotulojen kasvattamiseen; hieman alle puolet yrityksistä oli hakkuuseen painottuneita yrityksiä - muut päätoimialat olivat puukauppa, metsätalouden suunnittelupalvelut sekä mittauspalvelut; yrittäjän lisäksi yritykseen on palkattu keskimäärin kaksi vakinaista työntekijää; yrityksen kirjanpidosta huolehtii joku muu (useimmiten tilitoimisto) kuin yrittäjä; yrityksen kaluston määrä on kasvanut yrityksen perustamisvuodesta, mutta yrityksen asiakaskunta on pysynyt lähes samana; yrityksen tärkeimmän asiakkaan osuus yrityksen liikevaihdosta on 90%; tulevaisuudessa yrittäjä haluaa ylläpitää nykyisiä asiakassuhteita; ja kartoitettaessa yrittäjän halukkuutta ryhtyä uudelleen yrittäjäksi saamien kokemustensa pohjalta valtaosin yrittäjät sanoivat varmasti ryhtyvänsä uudelleen metsäpalveluyrittäjäksi. Vastaavasti Heikommin menestyneitä yrityksiä ja niiden toimintaa voitiin luonnehtia seuraavien tekijöiden avulla: yrittäjä oli ollut iältään yli 40-vuotias, kun hän oli perustanut metsäpalveluyrityksensä; yrittäjä oli toiminut maanviljelijänä ennen metsäpalveluyrittäjäuraansa; lähes puolet yrittäjistä olivat olleet melko epävarmoja yritystä käynnistäessään; yrityksen päätoimiala on useimmiten hakkuu; yrityksessä ei ole yrittäjän lisäksi palkattua, vakinaista henkilöstöä; yrityksen toiminta-alueen säde on alle 40 kilometriä; valtaosin yrittäjät hoitavat itse yrityksensä kirjanpidon; yrityksen kaluston määrä ei ole kasvanut yrityksen perustamisvuodesta; yrityksellä on nykyisin alle viisi asiakasta; tärkeimmän asiakkaan osuus yrityksen liikevaihdosta on noin 60%; tulevaisuudessa yrittäjä suunnittelee hankkivansa uusia asiakkaita nykyisten asiakassuhteiden ylläpidon lisäksi; ja jos yrittäjä saisi tehdä uudelleen yrityksensä perustamispäätöksen, runsas kaksi kolmasosaa yrittäjistä kertoi harkitsevansa tarkkaan perustamispäätöstä, mutta perustavansa uudelleen metsäpalveluyrityksen.

Lasikuituputkien kelauslinjan automatisointi

Lasikuituputkien kelaamiseen tarkoitettu automaatio voidaan toteuttaa suorakulmaisen kolmion trigonometrian perusteella. Kelaaminen tapahtuu pyörittämällä muottia ja liikuttamalla samalla muottia sivuttain. Näin tehdessä saadaan kelattavat lasikuitulangat leviämään koko muotin leveydelle. Jos langat halutaan saada kelattua tiiviisti toistensa viereen, on muotin liikuttava tietty matka yhtä muotin pyörähdystä kohden. Matkaa kutsutaan nousuksi. Nousu riippuu lankojen muodostaman kimpun leveydestä, kelauskulmasta ja muotin säteestä. Kelauskulmalla voidaan vaikuttaa kelattavan tuotteen aksiaaliseen ja pitkittäiseen lujuuteen. Kelauskoneen automaatiojärjestelmä koostuu sähkökäytöistä, ohjaavasta yksiköstä, valvomosta ja erilaisista antureista. Eri osat yhdistetään toisiinsa väylätekniikalla. Sähkökäytöt ovat yleensä servokäyttöjä, koska niiltä vaaditaan tarkkuutta ja nopeutta. Ohjaavana yksikkönä toimii yleensä teollisuus logiikka. Anturointia tarvitaan esimerkiksi hartsin lämpötilan säätöön ja lankojen ratakireyden säätöön. Anturoinnilla ja säädöillä pyritään parantamaan kelattavan tuotteen laatua. Koneen ohjelmallinen toteutus vaatii kelauskoneen mekaniikan ja kelausprosessin ymmärtämistä. Lisäksi on huomioitava koneen käyttäjät ja heidän turvallisuus.

The Use of Radiostereometric Analysis in Fractures of the Distal Radius: From Phantom Models to Clinical Application

Radiostereometric analysis (RSA) is a highly accurate method for the measurement of in vivo micromotion of orthopaedic implants. Validation of the RSA method is a prerequisite for performing clinical RSA studies. Only a limited number of studies have utilised the RSA method in the evaluation of migration and inducible micromotion during fracture healing. Volar plate fixation of distal radial fractures has increased in popularity. There is still very little prospective randomised evidence supporting the use of these implants over other treatments. The aim of this study was to investigate the precision, accuracy, and feasibility of using RSA in the evaluation of healing in distal radius fractures treated with a volar fixed-angle plate. A physical phantom model was used to validate the RSA method for simple distal radius fractures. A computer simulation model was then used to validate the RSA method for more complex interfragmentary motion in intra-articular fractures. A separate pre-clinical investigation was performed in order to evaluate the possibility of using novel resorbable markers for RSA. Based on the validation studies, a prospective RSA cohort study of fifteen patients with plated AO type-C distal radius fractures with a 1-year follow-up was performed. RSA was shown to be highly accurate and precise in the measurement of fracture micromotion using both physical and computer simulated models of distal radius fractures. Resorbable RSA markers demonstrated potential for use in RSA. The RSA method was found to have a high clinical precision. The fractures underwent significant translational and rotational migration during the first two weeks after surgery, but not thereafter. Maximal grip caused significant translational and rotational interfragmentary micromotion. This inducible micromotion was detectable up to eighteen weeks, even after the achievement of radiographic union. The application of RSA in the measurement of fracture fragment migration and inducible interfragmentary micromotion in AO type-C distal radius fractures is feasible but technically demanding. RSA may be a unique tool in defining the progress of fracture union.

Statistical variation of weld profiles and their expected influence on fatigue strength

The general objective of this study was to conduct astatistical analysis on the variation of the weld profiles and their influence on the fatigue strength of the joint. Weld quality with respect to its fatigue strength is of importance which is the main concept behind this thesis. The intention of this study was to establish the influence of weld geometric parameters on the weld quality and fatigue strength. The effect of local geometrical variations of non-load carrying cruciform fillet welded joint under tensile loading wasstudied in this thesis work. Linear Elastic Fracture Mechanics was used to calculate fatigue strength of the cruciform fillet welded joints in as-welded condition and under cyclic tensile loading, for a range of weld geometries. With extreme value statistical analysis and LEFM, an attempt was made to relate the variation of the cruciform weld profiles such as weld angle and weld toe radius to respective FAT classes.

Fatigue Analysis of Welded Cruciform Joints with Cold Laps

Työssä on tutkittu vetojännityskuormituksen alaisena olevien hitsattujen kuormaa kantamattomien X-liitosten hitsin paikallisen geometrian variaation vaikutusta väsymislujuuteen. Muuttujina olivat reunan pyöristyssäde, kylmäjuoksun suuruus ja kylkikulma. Geometristen muuttujien parametrinen riippuvuussuhde on analysoitu usealla elementtimallilla. Väsymistarkastelu on suoritettu käyttämällä lineaaris-elastista murtumismekaniikkaa (LEFM) tasovenymätilassa ja materiaalina terästä. Särönkasvun suunnan ennustamisessaon käytetty maksimipääjännityskriteeriä sekä jännitysintensiteettikertoimet on määritetty J-integraalilla. Särön ydintymisvaihetta ei ole otettu huomioon. Rakenteen on oletettu olevan hitsatussa tilassa ja jännitysheilahdus on kokonaan tehollinen. Särön kasvunopeuden ennustamiseen on käytetty Paris'n lakia. Väsymislujuustulokset on esitetty karakteristisina väsymisluokkina (FAT) ja sovitettu parametriseksi yhtälöksi. Lopuksi väsymisanalyysin ennustamia tuloksia on verrattu saatavilla oleviin väsytystestituloksiin.

PEAP 802.1x-tunnistuksen toteutus RADIUS-palvelimeen

Valimo iDServer -tunnistuspalvelin on ohjelmisto, joka tukee eri käyttäjientunnistusmenetelmiä, kuten tekstiviestillä lähetettävää kertakäyttösalasanaa tai normaaliakäyttäjätunnusta ja salasanaa. Tässä diplomityössä kuvataan, kuinka palvelimeenon lisätty tuki käyttäjien kirjautumiselle langattoman verkon tukiasemien ja virtuaalilähiverkkoa tukevien kytkinten kautta käyttäen normaaleja Windows-käyttöjärjestelmän mukana tulevia asiakasohjelmistoja. Työn ensimmäisessä vaiheessa kuvataan lähtökohdat ja vaatimukset tulevalle järjestelmälle. Työn osana käytännössä tehty kokonaisuus muodostuu useista eri määrityksistä koostuvista osista. Työn toisessa vaiheessa käydään läpi korkealla tasolla sovelluksen vaatimat protokollat. Osana näihin protokolliin kuului erilaisten avainten jatarkisteiden laskenta sekä salausmenetelmien käyttö, jotka myös kuvataan tässä työssä. Viimeisessä kappaleessa analysoidaan työn tuloksia jakäydään läpi toteutukseen ja itse sovelluksen toimintaan liittyvät ongelmat. Suurin osa havaituista ongelmista liittyi tilanteisiin, joihin itse palvelinsovelluksen toteutuksella ei voitu vaikuttaa. Eniten ongelmia aiheuttivat asiakasohjelmiston sekä verkkokorttien ja niiden ajureiden toiminta ongelmatilanteissa. Asiakasohjelmistoa ei selkeästi ole suunniteltu käytettäväksi kuin muuttumattomien salasanojen kanssa, koska käyttäjän näkökulmasta käyttökokemus ei ollut optimaalinen. Ongelmista huolimatta työn tuloksena saatiin asiakkaan vaatimukset täyttävä järjestelmä. Myös tuotekehitysnäkökulmasta projektia voitaneen pitää onnistuneena, koska nyt tehty sovellus luo pohjan uusien tunnistustapojen ja menetelmien toteuttamiselle tuotteen jatkokehitystä ajatellen.

Computational Analysis and Optimization of Real Gas Flow In Small Centrifugal Compressors

Small centrifugal compressors are more and more widely used in many industrialsystems because of their higher efficiency and better off-design performance comparing to piston and scroll compressors as while as higher work coefficient perstage than in axial compressors. Higher efficiency is always the aim of the designer of compressors. In the present work, the influence of four partsof a small centrifugal compressor that compresses heavy molecular weight real gas has been investigated in order to achieve higher efficiency. Two parts concern the impeller: tip clearance and the circumferential position of the splitter blade. The other two parts concern the diffuser: the pinch shape and vane shape. Computational fluid dynamics is applied in this study. The Reynolds averaged Navier-Stokes flow solver Finflo is used. The quasi-steady approach is utilized. Chien's k-e turbulence model is used to model the turbulence. A new practical real gas model is presented in this study. The real gas model is easily generated, accuracy controllable and fairly fast. The numerical results and measurements show good agreement. The influence of tip clearance on the performance of a small compressor is obvious. The pressure ratio and efficiency are decreased as the size of tip clearance is increased, while the total enthalpy rise keeps almost constant. The decrement of the pressure ratio and efficiency is larger at higher mass flow rates and smaller at lower mass flow rates. The flow angles at the inlet and outlet of the impeller are increased as the size of tip clearance is increased. The results of the detailed flow field show that leakingflow is the main reason for the performance drop. The secondary flow region becomes larger as the size of tip clearance is increased and the area of the main flow is compressed. The flow uniformity is then decreased. A detailed study shows that the leaking flow rate is higher near the exit of the impeller than that near the inlet of the impeller. Based on this phenomenon, a new partiallyshrouded impeller is used. The impeller is shrouded near the exit of the impeller. The results show that the flow field near the exit of the impeller is greatly changed by the partially shrouded impeller, and better performance is achievedthan with the unshrouded impeller. The loading distribution on the impeller blade and the flow fields in the impeller is changed by moving the splitter of the impeller in circumferential direction. Moving the splitter slightly to the suction side of the long blade can improve the performance of the compressor. The total enthalpy rise is reduced if only the leading edge of the splitter ismoved to the suction side of the long blade. The performance of the compressor is decreased if the blade is bended from the radius direction at the leading edge of the splitter. The total pressure rise and the enthalpy rise of thecompressor are increased if pinch is used at the diffuser inlet. Among the fivedifferent pinch shape configurations, at design and lower mass flow rates the efficiency of a straight line pinch is the highest, while at higher mass flow rate, the efficiency of a concave pinch is the highest. The sharp corner of the pinch is the main reason for the decrease of efficiency and should be avoided. The variation of the flow angles entering the diffuser in spanwise direction is decreased if pinch is applied. A three-dimensional low solidity twisted vaned diffuser is designed to match the flow angles entering the diffuser. The numerical results show that the pressure recovery in the twisted diffuser is higher than in a conventional low solidity vaned diffuser, which also leads to higher efficiency of the twisted diffuser. Investigation of the detailed flow fields shows that the separation at lower mass flow rate in the twisted diffuser is later than in the conventional low solidity vaned diffuser, which leads to a possible wider flow range of the twisted diffuser.

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.

Analytical Study on the Influence of Geometrical Parameters on the Fatigue Strength of Welds with Undercuts

The purpose of this study is to analyze the influence of geometric parameters on the fatigue strength of the welded joint. The thesis includes the analytical study pertaining to variation of undercut and the weld toe radius. The theoretical part includes the concepts of fracture mechanics, different stages of the crack propagation and finally the results. The numerical results are obtained from FRANC2D/L. The results show that the fatigue strength of the welded structure depends on the parameters of the welds such as the weld toe radius, weld angle, height of the undercut and plate thickness. The results show that the fatigue life can be improved by having the undercut as minimum as possible. However there has been significant variation on the fatigue life of the structure when the weld toe radius is increased.

Calculation of the critical crack size for cold form rectangular hollow section tube (CRHS) under bending load at -40º C temperature

To predict the capacity of the structure or the point which is followed by instability, calculation of the critical crack size is important. Structures usually contain several cracks but not necessarily all of these cracks lead to failure or reach the critical size. So, defining the harmful cracks or the crack size which is the most leading one to failure provides criteria for structure’s capacity at elevated temperature. The scope of this thesis was to calculate fracture parameters like stress intensity factor, the J integral and plastic and ultimate capacity of the structure to estimate critical crack size for this specific structure. Several three dimensional (3D) simulations using finite element method by Ansys program and boundary element method by Frank 3D program were carried out to calculate fracture parameters and results with the aid of laboratory tests (loaddisplacement curve, the J resistance curve and yield or ultimate stress) leaded to extract critical size of the crack. Two types of the fracture which is usually affected by temperature, Elastic and Elasti-Plastic fractures were simulated by performing several linear elastic and nonlinear elastic analyses. Geometry details of the weldment; flank angle and toe radius were also studied independently to estimate the location of crack initiation and simulate stress field in early stages of crack extension in structure. In this work also overview of the structure’s capacity in room temperature (20 ºC) was studied. Comparison of the results in different temperature (20 ºC and -40 ºC) provides a threshold of the structure’s behavior within the defined range.