Pyörivän keskiökärjen rakenteen kehittäminen työstövärähtelyjen välttämiseksi sorvauksessa

Tässä diplomityössä käsitellään sorvauksen työstövärähtelyjen ja sorvin keskiökärjen rakenteen yhteyttä. Työ on osa Lappeenrannan teknillisen yliopiston VMAX-projektia, ja sen taustalla on pyrkimys uudenlaisen, sorvin kärkipylkän puristusvoiman ajonaikaiseen säätämiseen perustuvan työstövärähtelyjen välttämismenetelmän kehittämiseen. Tämän menetelmän toiminnan todentaminen oli työn ensimmäinen tavoite. Menetelmän toteuttaminen asettaa kuitenkin käytetyn keskiökärjen rakenteelle tiettyjä vaatimuksia. Työn toisena tavoitteena oli nämä vaatimukset täyttävän keskiökärjen prototyypin kehittäminen. Tutkimus eteni seuraavasti. Ensimmäiseksi ongelma määriteltiin tutustumalla työn teoreettiseen taustaan ja aiheeseen liittyvään tutkimukseen Lappeenrannan teknillisestä yliopistosta ja muualta. Myös keskiökärkiä valmistavien yritysten tuotekatalogeja tarkasteltiin. Seuraavaksi siirryttiin alustavaan suunnitteluvaiheeseen, jossa verifioitiin menetelmän toiminta ja luotiin konsepteja keskiökärjen rakenteen kehittämistä varten. Tämän alustavan vaiheen jälkeen suoritettiin suunnitteluprosessi keskiökärjen prototyypille. Lopuksi, suunnitellun prototyypin rakenteen käyttäytymistä arvioitiin tietokonemallinnuksen avulla. Lisätuloksena tutkimuksen aikana johdettiin yksinkertaistettu elementtimenetelmään perustuva laskentamalli järjestelmän ominaistaajuuksien selvittämiseksi. Laskentamallin tarkkuutta arvoitiin. Suunnitteluprosessin tuloksena saatiin kaikki menetelmän toiminnan sekä normaalin käytön asettamat vaatimukset täyttävä rakenne keskiökärjen prototyypille. Myös johdetun laskentamallin tulokset ovat varsin lähellä 3D-elementtimallinnuksen antamia tuloksia. Tutkimuksen tavoitteiden voidaan siis sanoa toteutuneen. Koska prototyyppiä ja laskentamallia ei kuitenkaan ole vielä kokeellisesti verifioitu, tämä ei ole täysin varmaa.

Elementtikoon vaikutus väsymismenetelmien tarkkuuteen

Kandidaatintyössä tutkittiin kolmen hitsatun palstalevyliitoksen väsymiskestävyyttä äärellisten elementtien menetelmällä. Tutkimuksen pääpaino oli 2D-mallin antamien tulosten suppenevuudessa kolmella eri elementtikoolla ja kolmessa eri mittasuhdetapauksessa. Tutkimustuloksissa havaittiin noin 5 % muutos harvimman ja tiheimmän elementtiverkon tulosten välillä. Tutkimuksessa havaittiin myös että muutos harvimman ja tiheimmän elementtiverkon tuloksissa pieneni hiukan mitä suurempia kappaleet ja niiden hitsausliitokset olivat.

Saline distribution during multicomponent salting in pre-cooked quail eggs

The relationship of NaCl with problems of arterial hypertension has led to a reduction in the levels of this salt in food production. KCl has been used as a partial substitute for NaCl since it cannot be completely substituted without affecting the acceptability of the end product. In this study, the diffusion that occurs during quail egg salting in static and stirred brine was simulated. The mathematical model used was based on a generalization of the Fick's 2nd law, and the COMSOL Multiphysics software was used to simulate the diffusion in the NaCl-KCl-water system. The deviations in the simulated data and experimental data were 2.50% for NaCl and 6.98% for KCl in static brine, while in the stirred brine they were 3.48% for NaCl and 4.72% for KCl. The simulation results presented good agreement with the experimental values and validated the predictive capacity of the model.

Chassis frame design of a mobile platform equipped with robotic arms

The overall objective of the thesis is to design a robot chassis frame which is a bearing structure of a vehicle supporting all mechanical components and providing structure and stability. Various techniques and scientific principles were used to design a chassis frame.Design principles were applied throughout the process. By using Solid-Works software,virtual models was made for chassis frame. Chassis frame of overall dimension 1597* 800*950 mm3 was designed. Center of mass lieson 1/3 of the length from front wheel at height 338mm in the symmetry plane. Overall weight of the chassis frame is 80.12kg. Manufacturing drawing is also provided. Additionally,structural analysis was done in FEMAP which gives the busting result for chassis design by taking into consideration stress and deflection on different kind of loading resembling real life case. On the basis of simulated result, selected material was verified. Resulting design is expected to perform its intended function without failure. As a suggestion for further research, additional fatigue analysis and proper dynamic analysis can be conducted to make the study more robust.

Sellumassan varastosäiliön rakenneanalyysi

Tämä työ käsittelee valkaistun ja valkaisemattoman sellumassan varastosäiliön tärinää ja värähtelyä. Värähtelyn seurauksena säiliön seinämän jäykisterenkaan hitsausliitokseen on syntynyt särö. Työn tavoitteena on selvittää, mikä johtaa särön syntyyn ja miten kestäväm-piä varastosäiliöitä voidaan rakentaa materiaali-, valmistus- tai rakennemuutoksien avulla. Työ alkoi tutkimalla rakennemateriaalina olevan duplex-teräksen mikrorakenteen ominai-suuksia, sekä sen hitsattavuutta ja seostamista kirjallisuustutkimuksena. Kirjallisuustutki-musta jatkettiin selvittämällä mahdollisia vaurion syntymekanismeja seinämän särölle. Työssä analysoitiin myös tehtaalla mitattuja värähtelyarvoja. Lopuksi laskettiin FE-analyysillä tyhjän varastosäiliön ominaismuodot ja -taajuudet moodianalyysillä, sekä selvi-tettiin harmonisella analyysillä pinnankorkeuden vaihtelun vaikutus siirtymävasteeseen ja kriittisiin värähtelytaajuuksiin. Varastosäiliöön kohdistuvaa värähtelyä ei ole mahdollista poistaa kokonaan, mutta väräh-telyn aiheuttamia seurauksia kyetään rajaamaan useilla keinoilla. Toimenpiteinä voivat olla ainakin seinämän materiaalin paksuuden lisääminen, jäykisteripojen lisääminen ja kriittisten sellun pinnankorkeuksien välttäminen. Kriittiseksi pinnankorkeudeksi havaittiin 40–45 %:n täyttöaste ja turvalliseksi korkeudeksi 35–38 %:n täyttöaste. Varastosäiliölle kriittisen ominaistaajuuden katsotaan syntyvän taajuuksilla 3,3–3,8 Hz ja 5,8–6,4 Hz. Sellumassa putoaa varastosäiliöön noin 2 Hz taajuudella.

Evaluation of the Efficiency of Line-Start Permanent-Magnet Machines as a Function of the Operating Temperature

The standard squirrel-cage induction machine has nearly reached its maximum efficiency. In order to further increase the energy efficiency of electrical machines, the use of permanent magnets in combination with the robust design and the line start capability of the induction machine is extensively investigated. Many experimental designs have been suggested in literature, but recently, these line-start permanent-magnet machines (LSPMMs) have become off-the-shelf products available in a power range up to 7.5 kW. The permanent magnet flux density is a function of the operating temperature. Consequently, the temperature will affect almost every electrical quantity of the machine, including current, torque, and efficiency. In this paper, the efficiency of an off-the-shelf 4-kW three-phase LSPMM is evaluated as a function of the temperature by both finite-element modeling and by practical measurements. In order to obtain stator, rotor, and permanent magnet temperatures, lumped thermal modeling is used.

Convenience of scallops in plate structures

In this thesis was performed comprehensive study about the convenience of scallops in plate structures. A literature review was performed and lack of knowledge was fulfilled with fatigue tests performed in the laboratory of Steel Structures at the Lappeenranta University of Technology and with finite element method. The aim of this thesis was to produce design guidance for the use of scallops for different structural details and different loading conditions. An additional aim was to include more precise instructions for scallop design to produce good fatigue resistance and appropriate manufacturing quality. The literature review was performed searching bridge engineering and maritime standards and design guides and studies from scientific databases and reference lists from the literature of this field. Fatigue tests were used to research the effect of using scallops or not using scallops to fatigue strength of bracket specimen. Tests were performed on three specimens with different scallop radii and to five specimens without scallops with different weld penetration depths. Finite element method using solid elements, symmetry and submodels was used to determine stress concentration factors for I-beams with scallops. Stresses were defined with hot spot stress method. Choosing to use a scallop or not in the structure is affected by many factors, such as structural and loading conditions and manufacturability. As a rule of thumb, scallops should be avoided because those cause stress concentration points to the structure and take a lot of time to manufacture. When scallops are not used, good quality welding should be provided and full weld penetration is recommended to be used in load carrying corner weld areas. In some cases, it is advisable to use scallops. In that case, circular scallops are recommended to be used and radius should be chosen from fatigue strength or manufacturing point of view.

Investigation of the fatigue performance of welded duplex stainless steel

Diplomityössä tutkitaan hitsatun duplex-teräksen, laatu: EN 1.4462 (Outokumpu laatu 2205) väsymislujuutta. Tutkimusmetodologia noudattaa sekä kokeellisia että laskennallisia menetelmiä. Kokeelliset menetelmät sisältävät hitsatun teräksen väsytystestaukset laboratoriossa, hitsausten jälkikäsittelyt (HiFIT) sekä perusaineelle ja hitseille tehtävät metallurgiset tutkimukset. Väsytyskokeista saatavia tuloksia verrataan kansainvälisen hitsausinstituutin (IIW) vahvistamiin rakennekohtaisiin standardeihin sekä kirjallisuudessa esiintyviin tutkimustuloksiin. Laskennalliset menetelmät sisältävät vertailulaskelmia tehollisen lovijännityksen (ENS) menetelmää hyödyntäen. Tehollisen lovijännityksen menetelmässä liitoksissa vaikuttavat teholliset lovijännitykset selvitetään elementtimenetelmän (FEM) avulla. Tulokset vahvistavat, että hitsauksella ja hitsausten jälkikäsittelyllä on suuri merkitys rakenteen kestoikään. Suurin osa väsytyskokeiden tuloksista osoitti parempia väsymiskestävyyden arvoja kuin rakennekohtaiset standardit, mutta liitosten liitosvirheiden todettiin heikentävän väsytyskestävyyttä. Jälkikäsittelyiden todettiin parantavan liitosten väsymiskestävyyden tuloksia ja todettiin tulosten olevan hyödynnettävissä mitoituksessa.

Mechanical simulation of the endovascular repair of abdominal aortic aneurysms

Ce travail de thèse porte sur la simulation du déploiement des prothèses vasculaires de type stent-graft (SG) lors de la réparation endovasculaire (EVAR) des anévrismes de l’aorte abdominale (AAA). Cette étude se présente en trois parties: (i) tests mécaniques en flexion et compression de SG couramment utilisés (corps et jambage de marque Cook) ainsi que la simulation numérique desdits tests, (ii) développement d’un modèle numérique d’anévrisme, (iii) stratégie de simulation du déploiement des SG. La méthode numérique employée est celle des éléments finis. Dans un premier temps, une vérification du modèle éléments finis (MEF) des SG est realisée par comparaison des différents cas de charge avec leur pendant expérimental. Ensuite, le MEF vasculaire (AAA) est lui aussi vérifié lors d’une comparaison des niveaux de contraintes maximales principales dans la paroi avec des valeurs de la littérature. Enfin, le déploiement est abordé tout en intégrant les cathéters. Les tests mécaniques menés sur les SG ont été simulés avec une différence maximale de 5,93%, tout en tenant compte de la pré-charge des stents. Le MEF de la structure vasculaire a montré des contraintes maximales principales éloignées de 4,41% par rapport à un modèle similaire précédemment publié. Quant à la simulation du déploiement, un jeu complet de SG a pu être déployé avec un bon contrôle de la position relative et globale, dans un AAA spécifique pré-déformé, sans toutefois inclure de thrombus intra-luminal (TIL). La paroi du AAA a été modélisée avec une loi de comportement isotropique hyperélastique. Étant donné que la différence maximale tolérée en milieu clinique entre réalité et simulation est de 5%, notre approche semble acceptable et pourrait donner suite à de futurs développements. Cela dit, le petit nombre de SG testés justifie pleinement une vaste campagne de tests mécaniques et simulations supplémentaires à des fins de validation.

Développement d'une nouvelle méthode de caractérisation tissulaire basée sur l'élastographie ultrasonore : application pour le dépistage précoce du cancer du sein

Le cancer du sein est le cancer le plus fréquent chez la femme. Il demeure la cause de mortalité la plus importante chez les femmes âgées entre 35 et 55 ans. Au Canada, plus de 20 000 nouveaux cas sont diagnostiqués chaque année. Les études scientifiques démontrent que l'espérance de vie est étroitement liée à la précocité du diagnostic. Les moyens de diagnostic actuels comme la mammographie, l'échographie et la biopsie comportent certaines limitations. Par exemple, la mammographie permet de diagnostiquer la présence d’une masse suspecte dans le sein, mais ne peut en déterminer la nature (bénigne ou maligne). Les techniques d’imagerie complémentaires comme l'échographie ou l'imagerie par résonance magnétique (IRM) sont alors utilisées en complément, mais elles sont limitées quant à la sensibilité et la spécificité de leur diagnostic, principalement chez les jeunes femmes (< 50 ans) ou celles ayant un parenchyme dense. Par conséquent, nombreuses sont celles qui doivent subir une biopsie alors que leur lésions sont bénignes. Quelques voies de recherche sont privilégiées depuis peu pour réduire l`incertitude du diagnostic par imagerie ultrasonore. Dans ce contexte, l’élastographie dynamique est prometteuse. Cette technique est inspirée du geste médical de palpation et est basée sur la détermination de la rigidité des tissus, sachant que les lésions en général sont plus rigides que le tissu sain environnant. Le principe de cette technique est de générer des ondes de cisaillement et d'en étudier la propagation de ces ondes afin de remonter aux propriétés mécaniques du milieu via un problème inverse préétabli. Cette thèse vise le développement d'une nouvelle méthode d'élastographie dynamique pour le dépistage précoce des lésions mammaires. L'un des principaux problèmes des techniques d'élastographie dynamiques en utilisant la force de radiation est la forte atténuation des ondes de cisaillement. Après quelques longueurs d'onde de propagation, les amplitudes de déplacement diminuent considérablement et leur suivi devient difficile voir impossible. Ce problème affecte grandement la caractérisation des tissus biologiques. En outre, ces techniques ne donnent que l'information sur l'élasticité tandis que des études récentes montrent que certaines lésions bénignes ont les mêmes élasticités que des lésions malignes ce qui affecte la spécificité de ces techniques et motive la quantification de d'autres paramètres mécaniques (e.g.la viscosité). Le premier objectif de cette thèse consiste à optimiser la pression de radiation acoustique afin de rehausser l'amplitude des déplacements générés. Pour ce faire, un modèle analytique de prédiction de la fréquence de génération de la force de radiation a été développé. Une fois validé in vitro, ce modèle a servi pour la prédiction des fréquences optimales pour la génération de la force de radiation dans d'autres expérimentations in vitro et ex vivo sur des échantillons de tissu mammaire obtenus après mastectomie totale. Dans la continuité de ces travaux, un prototype de sonde ultrasonore conçu pour la génération d'un type spécifique d'ondes de cisaillement appelé ''onde de torsion'' a été développé. Le but est d'utiliser la force de radiation optimisée afin de générer des ondes de cisaillement adaptatives, et de monter leur utilité dans l'amélioration de l'amplitude des déplacements. Contrairement aux techniques élastographiques classiques, ce prototype permet la génération des ondes de cisaillement selon des parcours adaptatifs (e.g. circulaire, elliptique,…etc.) dépendamment de la forme de la lésion. L’optimisation des dépôts énergétiques induit une meilleure réponse mécanique du tissu et améliore le rapport signal sur bruit pour une meilleure quantification des paramètres viscoélastiques. Il est aussi question de consolider davantage les travaux de recherches antérieurs par un appui expérimental, et de prouver que ce type particulier d'onde de torsion peut mettre en résonance des structures. Ce phénomène de résonance des structures permet de rehausser davantage le contraste de déplacement entre les masses suspectes et le milieu environnant pour une meilleure détection. Enfin, dans le cadre de la quantification des paramètres viscoélastiques des tissus, la dernière étape consiste à développer un modèle inverse basé sur la propagation des ondes de cisaillement adaptatives pour l'estimation des paramètres viscoélastiques. L'estimation des paramètres viscoélastiques se fait via la résolution d'un problème inverse intégré dans un modèle numérique éléments finis. La robustesse de ce modèle a été étudiée afin de déterminer ces limites d'utilisation. Les résultats obtenus par ce modèle sont comparés à d'autres résultats (mêmes échantillons) obtenus par des méthodes de référence (e.g. Rheospectris) afin d'estimer la précision de la méthode développée. La quantification des paramètres mécaniques des lésions permet d'améliorer la sensibilité et la spécificité du diagnostic. La caractérisation tissulaire permet aussi une meilleure identification du type de lésion (malin ou bénin) ainsi que son évolution. Cette technique aide grandement les cliniciens dans le choix et la planification d'une prise en charge adaptée.

A Variational Calculation of TE and TM Cutoff Wavenumbers in Circular Eccentric Guides by Conformal Mapping

The cutoff wavenumbers of higher order modes in circular eccentric guides are computed with the variational analysis combined with a conformal mapping. A conformal mapping is applied to the variational formulation, and the variational equation is solved by the finite-element method. Numerical results for TE and TM cutoff wavenumbers are presented for different distances between the centers and ratio of the radii. Comparisons with numerical results found in the literature validate the presented method

Development and Analysis of a Novel Isosceles Trapezoidal Dielectric Resonator Antenna For Wireless Communication

This thesis describes the development and analysis of an Isosceles Trapezoidal Dielectric Resonator Antenna (ITDRA) by realizing different DR orientations with suitable feed configurations enabling it to be used as multiband, dual band dual polarized and wideband applications. The motivation for this work has been inspired by the need for compact, high efficient, low cost antenna suitable for multi band application, dual band dual polarized operation and broadband operation with the possibility of using with MICs, and to ensure less expensive, more efficient and quality wireless communication systems. To satisfy these challenging demands a novel shaped Dielectric Resonator (DR) is fabricated and investigated for the possibility of above required properties by trying out different orientations of the DR on a simple microstrip feed and with slotted ground plane as well. The thesis initially discusses and evaluates recent and past developments taken place within the microwave industry on this topic through a concise review of literature. Then the theoretical aspects of DRA and different feeding techniques are described. Following this, fabrication and characterization of DRA is explained. To achieve the desired requirements as above both simulations and experimental measurements were undertaken. A 3-D finite element method (FEM) electromagnetic simulation tool, HFSSTM by Agilent, is used to determine the optimum geometry of the dielectric resonator. It was found to be useful in producing approximate results although it had some limitations. A numerical analysis technique, finite difference time domain (FDTD) is used for validating the results of wide band design at the end. MATLAB is used for modeling the ITDR and implementing FDTD analysis. In conclusion this work offers a new, efficient and relatively simple alternative for antennas to be used for multiple requirements in the wireless communication system.

Investigations on Structural response of Water Backed Perforated plates

The study envisaged herein contains the numerical investigations on Perforated Plate (PP) as well as numerical and experimental investigations on Perforated Plate with Lining (PPL) which has a variety of applications in underwater engineering especially related to defence applications. Finite element method has been adopted as the tool for analysis of PP and PPL. The commercial software ANSYS has been used for static and free vibration response evaluation, whereas ANSYS LS-DYNA has been used for shock analysis. SHELL63, SHELL93, SOLID45, SOLSH190, BEAM188 and FLUID30 finite elements available in the ANSYS library as well as SHELL193 and SOLID194 available in the ANSYS LS-DYNA library have been made use of. Unit cell of the PP and PPL which is a miniature of the original plate with 16 perforations have been used. Based upon the convergence characteristics, the utility of SHELL63 element for the analysis of PP and PPL, and the required mesh density are brought out. The effect of perforation, geometry and orientation of perforation, boundary conditions and lining plate are investigated for various configurations. Stress concentration and deflection factor are also studied. Based on these investigations, stadium geometry perforation with horizontal orientation is recommended for further analysis.Linear and nonlinear static analysis of PP and PPL subjected to unit normal pressure has been carried out besides the free vibration analysis. Shock analysis has also been carried out on these structural components. The analytical model measures 0.9m x 0.9m with stiffener of 0.3m interval. The influence of finite element, boundary conditions, and lining plate on linear static response has been estimated and presented. Comparison of behavior of PP and PPL in the nonlinear strain regime has been made using geometric nonlinear analysis. Free vibration analysis of the PP and PPL has been carried out ‘in vacuum’ condition and in water backed condition, and the influence of water backed condition and effect of perforation on natural frequency have been investigated.Based upon the studies on the vibration characteristics of NPP, PP and PPL in water backed condition and ‘in vacuum’ condition, the reduction in the natural frequency of the plate in immersed condition has been rightly brought out. The necessity to introduce the effect of water medium in the analysis of water backed underwater structure has been highlighted.Shock analysis of PP and PPL for three explosives viz., PEK, TNT and C4 has been carried out and deflection and stresses on plate as well as free field pressure have been estimated using ANSYS LS-DYNA. The effect of perforations and the effect of lining plate have been predicted. Experimental investigations of the measurement of free field pressure using PPL have been conducted in a shock tank. Free field pressure has been measured and has been validated with finite element analysis results. Besides, an experiment has been carried out on PPL, for the comparison of the static deflection predicted by finite element analysis.The distribution of the free field pressure and the estimation of differential pressure from experimentation and the provision for treating the differential pressure as the resistance, as a part of the design load for PPL, has been brought out.

Performance Assessment of Sandwich Structures with Debonds and Dents

A sandwich construction is a special form of the laminated composite consisting of light weight core, sandwiched between two stiff thin face sheets. Due to high stiffness to weight ratio, sandwich construction is widely adopted in aerospace industries. As a process dependent bonded structure, the most severe defects associated with sandwich construction are debond (skin core bond failure) and dent (locally deformed skin associated with core crushing). Reasons for debond may be attributed to initial manufacturing flaws or in service loads and dent can be caused by tool drops or impacts by foreign objects. This paper presents an evaluation on the performance of honeycomb sandwich cantilever beam with the presence of debond or dent, using layered finite element models. Dent is idealized by accounting core crushing in the core thickness along with the eccentricity of the skin. Debond is idealized using multilaminate modeling at debond location with contact element between the laminates. Vibration and buckling behavior of metallic honeycomb sandwich beam with and without damage are carried out. Buckling load factor, natural frequency, mode shape and modal strain energy are evaluated using finite element package ANSYS 13.0. Study shows that debond affect the performance of the structure more severely than dent. Reduction in the fundamental frequencies due to the presence of dent or debond is not significant for the case considered. But the debond reduces the buckling load factor significantly. Dent of size 8-20% of core thickness shows 13% reduction in buckling load capacity of the sandwich column. But debond of the same size reduced the buckling load capacity by about 90%. This underscores the importance of detecting these damages in the initiation level itself to avoid catastrophic failures. Influence of the damages on fundamental frequencies, mode shape and modal strain energy are examined. Effectiveness of these parameters as a damage detection tool for sandwich structure is also assessed

FEM based Virtual Prototyping and Design of Third Harmonic Excitation System for Low Voltage Brushless Alternators

Salient pole brushless alternators coupled to IC engines are extensively used as stand-by power supply units for meeting in- dustrial power demands. Design of such generators demands high power to weight ratio, high e ciency and low cost per KVA out- put. Moreover, the performance characteristics of such machines like voltage regulation and short circuit ratio (SCR) are critical when these machines are put into parallel operation and alterna- tors for critical applications like defence and aerospace demand very low harmonic content in the output voltage. While designing such alternators, accurate prediction of machine characteristics, including total harmonic distortion (THD) is essential to mini- mize development cost and time. Total harmonic distortion in the output voltage of alternators should be as low as possible especially when powering very sophis- ticated and critical applications. The output voltage waveform of a practical AC generator is replica of the space distribution of the ux density in the air gap and several factors such as shape of the rotor pole face, core saturation, slotting and style of coil disposition make the realization of a sinusoidal air gap ux wave impossible. These ux harmonics introduce undesirable e ects on the alternator performance like high neutral current due to triplen harmonics, voltage distortion, noise, vibration, excessive heating and also extra losses resulting in poor e ciency, which in turn necessitate de-rating of the machine especially when connected to non-linear loads. As an important control unit of brushless alternator, the excitation system and its dynamic performance has a direct impact on alternator's stability and reliability. The thesis explores design and implementation of an excitation i system utilizing third harmonic ux in the air gap of brushless al- ternators, using an additional auxiliary winding, wound for 1=3rd pole pitch, embedded into the stator slots and electrically iso- lated from the main winding. In the third harmonic excitation system, the combined e ect of two auxiliary windings, one with 2=3rd pitch and another third harmonic winding with 1=3rd pitch, are used to ensure good voltage regulation without an electronic automatic voltage regulator (AVR) and also reduces the total harmonic content in the output voltage, cost e ectively. The design of the third harmonic winding by analytic methods demands accurate calculation of third harmonic ux density in the air gap of the machine. However, precise estimation of the amplitude of third harmonic ux in the air gap of a machine by conventional design procedures is di cult due to complex geome- try of the machine and non-linear characteristics of the magnetic materials. As such, prediction of the eld parameters by conven- tional design methods is unreliable and hence virtual prototyping of the machine is done to enable accurate design of the third har- monic excitation system. In the design and development cycle of electrical machines, it is recognized that the use of analytical and experimental methods followed by expensive and in exible prototyping is time consum- ing and no longer cost e ective. Due to advancements in com- putational capabilities over recent years, nite element method (FEM) based virtual prototyping has become an attractive al- ternative to well established semi-analytical and empirical design methods as well as to the still popular trial and error approach followed by the costly and time consuming prototyping. Hence, by virtually prototyping the alternator using FEM, the important performance characteristics of the machine are predicted. Design of third harmonic excitation system is done with the help of results obtained from virtual prototype of the machine. Third harmonic excitation (THE) system is implemented in a 45 KVA ii experimental machine and experiments are conducted to validate the simulation results. Simulation and experimental results show that by utilizing third harmonic ux in the air gap of the ma- chine for excitation purposes during loaded conditions, triplen harmonic content in the output phase voltage is signi cantly re- duced. The prototype machine with third harmonic excitation system designed and developed based on FEM analysis proved to be economical due to its simplicity and has the added advan- tage of reduced harmonics in the output phase voltage.