Modelagem e simulação da desidratação osmótica em pedaços de abacaxi utilizando o método de elementos finitos

Water loss and sugar gain were modelling during the osmotic dehydration process of pieces of pineaplle. The transfer of solute to the fruit and the water to the solution was based on Fick's 2nd law. The three dimensional model was solved by the finite element method with the usage of the software COMSOL Multiphysics 3.2. The main and cross diffusion coefficients and the Biot number were determined on the simulation and the deviation between the experimental and the simulated data were 4,28% to sucrose and 1,66 to the water.

Kestomagneettigeneraattorin staattorin mekaanisen rakenteen analysointi ja kehitys

Työssä tutkittiin kestomagneettigeneraattorin staattoriytimen aksiaalisen puristamisen ja lämpötilan vaikutuksia staattorin jännityksiin ja muodonmuutoksiin. Lisäksi tutkittiin aksiaalisen puristusjännityksen vaikutusta staattorin aksiaalisuuntaisen ominaismuodon taajuuteen. Rakenteesta tehtiin elementtimalli, jonka avulla laskettiin rakenteen siirtymät, jännitykset ja ominaistaajuudet. Rakenteen jännityksiä ja siirtymiä tutkittiin myös yksinkertaisella jousikuvauksella, jonka tuloksia verrattiin elementtimallilla saatuihin tuloksiin. Työn tavoitteena oli kehittää uusi staattorirakenne. Elementtimenetelmän tulosten perusteella staattorirakennetta voidaan yksinkertaistaa poistamalla nyt käytetty puristusrengas, jolloin rakenteesta tulee yksinkertaisempi ja kustannustehokkaampi.

Polymeeripinnoitteen aiheuttamat itseherätteiset värähtelyt nipillisessä telasysteemissä

Lukuisissa teollisuussovelluksissa materiaalien, kuten paperin ja teräslevyjen, muokkaamiseen käytettävät pyörivät nippitelat kärsivät aina erilaisten herätteiden synnyttämistä mekaanisista värähtelyistä, jotka voivat aiheuttaa virheitä valmistettaviin tuotteisiin. Tässä työssä tutkittiin viskoelastisia polymeerejä ja polymeeripinnoitteen nipilliseen telasysteemiin synnyttämiä haitallisia itseherätteisiä värähtelyjä. Työn polymeerejä käsittelevässä kirjallisuusosassa luotiin katsaus amorfisten polymeerien fysikaalisiin ominaisuuksiin. Kokeellisessa osuudessa tutkittiin tarkemmin kahden amorfisen telapinnoitepolymeerin termoreologisia ja mekaanisia ominaisuuksia suoritettujen DMTA-mittausten perusteella. Sovittamalla toisen polymeerin master-käyrään yleistetty lineaarisen standardiaineen malli saatiin selville polymeerin mekaaniset parametrit ja approksimaatio sen relaksaatiospektrille. Telapinnoitteen nipilliseen systeemiin synnyttämiä itseherätteisiä värähtelyjä ja niiden seurauksia tarkasteltiin kahdelle telalle ja polymeeripinnoitteelle kehitetyn analyyttisen mallin ja numeeristen laskujen avulla. Pinnoite mallinnettiin lineaarisen standardiaineen mukaisesti. Telasysteemin parametrit määritettiin DMTA-mittaustuloksista ja systeemiä vastaavasta koelaitteesta kokeellisella moodianalyysillä ja elementtimenetelmällä. Numeerisesta stabiilisuusanalyysistä ja liikeyhtälöiden integroinneista saadut tulokset kertovat telapinnoitteen aaltomaisista deformaatiomuodoista ja niiden synnyttämistä taajuusalueittain esiintyvistä epästabiileista värähtelyistä. Telasysteemi on epästabiili pinnoitedeformaatiokuvion systeemiin aiheuttaman herätevoiman taajuuden ollessa lähellä systeemin korkeampaa ominaistaajuutta. Numeerisista tuloksista voitiin ennustaa nopean ja hitaan barringin olemassaolo.

Kokeellisen moodianalyysin jälkikäsittely ominaismuotojen visuaaliseksi tarkastelemiseksi

Kokeellinen moodianalyysi on kokeellinen menetelmä, jolla voidaan selvittää rakenteille ominaista värähtelyä. Työssä oli kolme tavoitetta, jotka pyrittiin saavuttamaan. Ensimmäinen tavoite oli ohjeistuksen luominen rakenteiden ominaismuotojen visuaaliseksi tarkastelemiseksi Lappeenrannan teknillisellä yliopistolla käytettävissä olevilla laitteistoilla ja ohjelmilla. Ohjeistuksen perustana on ajatus siitä, että kokeellinen moodianalyysi saataisiin tehokkaampaan käyttöön koneensuunnittelun opetuksessa. Ohjeistus tehtiin pääasiassa kuvien ja kuvia tukevien selitysten avulla. Työn toisena tavoitteena oli verrata rakenteen ominaisvärähtelyä, kun se oli tuettu vapaasti ja kun se oli tuettu kiinteästi sen luonnolliseen ympäristöön. Taajuuksia verrattaessa huomattiin, että eri tavoin tuetut rakenteet käyttäytyvät eri tavalla, mikä on otettava huomioon, kun tutkitaan kriittisiä värähtelyjä. Ominaisvärähtelyjä voidaan selvittää myös matemaattisesti esimerkiksi äärellisten elementtien menetelmällä. Työn kolmantena tavoitteena oli verrata elementtimenetelmällä ja kokeellisesti saatuja ominaisvärähtelyn arvoja. Elementtimallia pyrittiin tarkentamaan eri parametrejä muuttamalla niin, että ominaistaajuuksien arvot vastaavat mahdollisimman hyvin toisiaan. Tavoite saavutettiin.

Induction Motor Drive Energy Efficiency - Simulation and Analysis

A coupled system simulator, based on analytical circuit equations and a finite element method (FEM) model of the motor has been developed and it is used to analyse a frequency-converterfed industrial squirrel-cage induction motor. Two control systems that emulate the behaviour of commercial direct-torque-controlled (DTC) and vector-controlled industrial frequency converters have been studied, implemented in the simulation software and verified by extensive laboratory tests. Numerous factors that affect the operation of a variable speed drive (VSD) and its energy efficiency have been investigated, and their significance in the simulation of the VSD results has been studied. The dependency of the frequency converter, induction motor and system losses on the switching frequency is investigated by simulations and measurements at different speeds for both the vector control and the DTC. Intensive laboratory measurements have been carried out to verify the simulation results.

Miehistönkuljetusajoneuvon rungon suunnittelu

Tämän diplomityön tavoitteena oli suunnitella miehistönkuljetusajoneuvon runko. Rungosta suunniteltiin mahdollisimman hyvin energiaa absorboiva. Rakenne toteutettiin kennora-kenteena. Suunnittelussa sovellettiin koneensuunnittelun periaatteiden lisäksi energiaa ab-sorboivien rakenteiden suunnittelun periaatteita. Myös valmistustekniset näkökohdat otet-tiin huomioon. Rakenteessa hyödynnettiin Ruukki Oy:n Ramor 500 suojausterästä sekä OPTIM 500 MC terästä. Lisäksi erilaisten täyteaineiden käyttöä tutkittiin. Suunnittelun työkaluna käytettiin epälineaarista elementtimenetelmää, koska energiaa ab-sorboivien rakenteiden suunnittelussa on otettava huomioon materiaalien epälineaarinen käyttäytyminen. Rakenteen suunnittelu jakaantui viiteen vaiheeseen. Aluksi rakenteeseen kohdistuvat kuormitukset laskettiin elementtimenetelmän avulla. Esisuunnittelussa lasket-tiin plastisuusteorian avulla alustavasti tarvittavat materiaalipaksuudet. Tämän jälkeen ra-kenteen ydingeometria optimoitiin mahdollisimman hyvin energiaa absorboivaksi. Opti-moinnissa hyödynnettiin elementtimenetelmää. Seuraavassa vaiheessa varmistettiin raken-teen globaalit ominaisuudet. Lopuksi rakenteen kestävyyttä tarkasteltiin elementtimene-telmällä. Runko ei mallien mukaan kestänyt siltä vaadittuja kuormitustapauksia. Mallin kaikki ole-tukset pidettiin varmalla puolella. Reunaehdot oletettiin todellisuutta jäykemmiksi. Myös-kään materiaalin venymänopeudesta johtuvaa lujittumista ei otettu huomioon. Koska mii-naräjähdys on monimutkainen tapahtuma, rungon todellinen kestävyys joudutaan ar-viomaan räjähdystesteillä. Elementtimallien perusteella voidaan kuitenkin sanoa, että ener-giaa absorboiva ajoneuvon runko on mahdollista toteuttaa kennorakenteena. Lisäksi voi-daan todeta, että elementtimenetelmää sopii työvälineeksi tämän tyyppisten rakenteiden suunnitteluun.

Broadband excitation in the system identification of active magnetic bearing rotor systems

One of the targets of the climate and energy package of the European Union is to increase the energy efficiency in order to achieve a 20 percent reduction in primary energy use compared with the projected level by 2020. The energy efficiency can be improved for example by increasing the rotational speed of large electrical drives, because this enables the elimination of gearboxes leading to a compact design with lower losses. The rotational speeds of traditional bearings, such as roller bearings, are limited by mechanical friction. Active magnetic bearings (AMBs), on the other hand, allow very high rotational speeds. Consequently, their use in large medium- and high-speed machines has rapidly increased. An active magnetic bearing rotor system is an inherently unstable, nonlinear multiple-input, multiple-output system. Model-based controller design of AMBs requires an accurate system model. Finite element modeling (FEM) together with the experimental modal analysis provides a very accurate model for the rotor, and a linearized model of the magneticactuators has proven to work well in normal conditions. However, the overall system may suffer from unmodeled dynamics, such as dynamics of foundation or shrink fits. This dynamics can be modeled by system identification. System identification can also be used for on-line diagnostics. In this study, broadband excitation signals are adopted to the identification of an active magnetic bearing rotor system. The broadband excitation enables faster frequency response function measurements when compared with the widely used stepped sine and swept sine excitations. Different broadband excitations are reviewed, and the random phase multisine excitation is chosen for further study. The measurement times using the multisine excitation and the stepped sine excitation are compared. An excitation signal design with an analysis of the harmonics produced by the nonlinear system is presented. The suitability of different frequency response function estimators for an AMB rotor system are also compared. Additionally, analytical modeling of an AMB rotor system, obtaining a parametric model from the nonparametric frequency response functions, and model updating are discussed in brief, as they are key elements in the modeling for a control design. Theoretical methods are tested with a laboratory test rig. The results conclude that an appropriately designed random phase multisine excitation is suitable for the identification of AMB rotor systems.

Application of Computational Fluid Dynamics in Modeling Blood Flow in Human Thoracic Aorta

The aim of this study was to simulate blood flow in thoracic human aorta and understand the role of flow dynamics in the initialization and localization of atherosclerotic plaque in human thoracic aorta. The blood flow dynamics in idealized and realistic models of human thoracic aorta were numerically simulated in three idealized and two realistic thoracic aorta models. The idealized models of thoracic aorta were reconstructed with measurements available from literature, and the realistic models of thoracic aorta were constructed by image processing Computed Tomographic (CT) images. The CT images were made available by South Karelia Central Hospital in Lappeenranta. The reconstruction of thoracic aorta consisted of operations, such as contrast adjustment, image segmentations, and 3D surface rendering. Additional design operations were performed to make the aorta model compatible for the numerical method based computer code. The image processing and design operations were performed with specialized medical image processing software. Pulsatile pressure and velocity boundary conditions were deployed as inlet boundary conditions. The blood flow was assumed homogeneous and incompressible. The blood was assumed to be a Newtonian fluid. The simulations with idealized models of thoracic aorta were carried out with Finite Element Method based computer code, while the simulations with realistic models of thoracic aorta were carried out with Finite Volume Method based computer code. Simulations were carried out for four cardiac cycles. The distribution of flow, pressure and Wall Shear Stress (WSS) observed during the fourth cardiac cycle were extensively analyzed. The aim of carrying out the simulations with idealized model was to get an estimate of flow dynamics in a realistic aorta model. The motive behind the choice of three aorta models with distinct features was to understand the dependence of flow dynamics on aorta anatomy. Highly disturbed and nonuniform distribution of velocity and WSS was observed in aortic arch, near brachiocephalic, left common artery, and left subclavian artery. On the other hand, the WSS profiles at the roots of branches show significant differences with geometry variation of aorta and branches. The comparison of instantaneous WSS profiles revealed that the model with straight branching arteries had relatively lower WSS compared to that in the aorta model with curved branches. In addition to this, significant differences were observed in the spatial and temporal profiles of WSS, flow, and pressure. The study with idealized model was extended to study blood flow in thoracic aorta under the effects of hypertension and hypotension. One of the idealized aorta models was modified along with the boundary conditions to mimic the thoracic aorta under the effects of hypertension and hypotension. The results of simulations with realistic models extracted from CT scans demonstrated more realistic flow dynamics than that in the idealized models. During systole, the velocity in ascending aorta was skewed towards the outer wall of aortic arch. The flow develops secondary flow patterns as it moves downstream towards aortic arch. Unlike idealized models, the distribution of flow was nonplanar and heavily guided by the artery anatomy. Flow cavitation was observed in the aorta model which was imaged giving longer branches. This could not be properly observed in the model with imaging containing a shorter length for aortic branches. The flow circulation was also observed in the inner wall of the aortic arch. However, during the diastole, the flow profiles were almost flat and regular due the acceleration of flow at the inlet. The flow profiles were weakly turbulent during the flow reversal. The complex flow patterns caused a non-uniform distribution of WSS. High WSS was distributed at the junction of branches and aortic arch. Low WSS was distributed at the proximal part of the junction, while intermedium WSS was distributed in the distal part of the junction. The pulsatile nature of the inflow caused oscillating WSS at the branch entry region and inner curvature of aortic arch. Based on the WSS distribution in the realistic model, one of the aorta models was altered to induce artificial atherosclerotic plaque at the branch entry region and inner curvature of aortic arch. Atherosclerotic plaque causing 50% blockage of lumen was introduced in brachiocephalic artery, common carotid artery, left subclavian artery, and aortic arch. The aim of this part of the study was first to study the effect of stenosis on flow and WSS distribution, understand the effect of shape of atherosclerotic plaque on flow and WSS distribution, and finally to investigate the effect of lumen blockage severity on flow and WSS distributions. The results revealed that the distribution of WSS is significantly affected by plaque with mere 50% stenosis. The asymmetric shape of stenosis causes higher WSS in branching arteries than in the cases with symmetric plaque. The flow dynamics within thoracic aorta models has been extensively studied and reported here. The effects of pressure and arterial anatomy on the flow dynamic were investigated. The distribution of complex flow and WSS is correlated with the localization of atherosclerosis. With the available results we can conclude that the thoracic aorta, with complex anatomy is the most vulnerable artery for the localization and development of atherosclerosis. The flow dynamics and arterial anatomy play a role in the localization of atherosclerosis. The patient specific image based models can be used to diagnose the locations in the aorta vulnerable to the development of arterial diseases such as atherosclerosis.

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.

Design and Control of Permanent Magnet Synchronous Machine for Hybrid Electric Vehicle

This master’s thesis mainly focuses on the design requirements of an Electric drive for Hybrid car application and its control strategy to achieve a wide speed range. It also emphasises how the control and performance requirements are transformed into its design variables. A parallel hybrid topology is considered where an IC engine and an electric drive share a common crank shaft. A permanent magnet synchronous machine (PMSM) is used as an electric drive machine. Performance requirements are converted into Machine design variables using the vector model of PMSM. Main dimensions of the machine are arrived using analytical approach and Finite Element Analysis (FEA) is used to verify the design and performance. Vector control algorithm was used to control the machine. The control algorithm was tested in a low power PMSM using an embedded controller. A prototype of 10 kW PMSM was built according to the design values. The prototype was tested in the laboratory using a high power converter. Tests were carried out to verify different operating modes. The results were in agreement with the calculations.

Fatigue study of the fillet weld waving on weld toe and the penetration in the weld root based on the effective notch stress approach

The fatigue failure of structures under fluctuating loads in fillet weld joints raises a demand to determine the parameters related to this type of loading. In this study, the stress distribution in the susceptible area of weld toe and weld root in fillet welded models analyzed by finite element method applying FEMAP software. To avoid the geometrical singularity on the path of analytical stress analysis in the toe and root area of a weld model the effective notch stress approach applied by which a proper fictitious rounding that mostly depend on the material of structure is applied. The models with different weld toe waving width and radius are analyzed while the flank angle of weld varied in 45 and 30 degrees. The processed results shows that the waving compare to the straight weld toe makes differences in the value of stress and consequently the stress concentration factor between the tip and depth of the waves in the weld toe which helps to protect the crack of propagation and gives enough time and tools to be informed of the crack initiation in the structure during the periodical observation of structure. In the weld root study the analyses among the models with the welding penetration percentage from non-penetration to the full-penetration shows a slightly increase in the root area stress value which comparing with the stiffening effect of penetration conclude that the half-penetration can make an optimization between the stress increase and stiffening effect of deep penetration.

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ä.

Liikkuvan kivimurskaimen syötinosan rakenneanalyysi osarakennetekniikan avulla

Työn tarkoituksena oli selvittää, miten osarakennetekniikkaa voidaan soveltaa siirrettävän kivimurskaimen syötinosan simuloinnissa. Tätä tutkittiin luomalla kahdella eri ohjelmistolla simulaatiomalli syötinosasta ja mallintamalla syötinosan runko joustavaksi kappaleeksi osarakennetekniikan avulla. Luotujen simulointimallien tarkkuutta selvitettiin vertaamalla niistä saatuja rungon jännityksiä tutkittavan rakenteen rungosta mitattuihin jännityksiin. Työn tarkoituksena oli myös tutkia, miten hyvin simulaatiomallit soveltuvat käytettäväksi syötinosan tuotekehityksessä. Tässä työssä käytettiin syötinosan simulaatiomallin luomiseen ANSYS-ohjelmistoa ja ADAMS-ohjelmistoa. Simulaatiomalleihin lisättiin tutkittavasta järjestelmästä mitattu ohjaussignaali sekä syötinosan jousien arvot. Järjestelmän rakenneominaisuudet saatiin suoraan valmistajan luovuttamista tiedoista. ADAMS-ohjelmistolla mallinnetussa simulaatiomallissa runko mallinnettiin joustavaksi ANSYS-ohjelmistossa, josta se siirrettiin ADAMS-ohjelmistoon. Saaduista tuloksista kävi ilmi, että osarakennetekniikkaa voidaan hyödyntää syötinosan joustavan rungon simuloinnissa. Tutkittavasta järjestelmästä mitatuissa jännityksissä ja simulaatiomalleista saaduissa jännityksissä oli eroja, mutta jännityshistorian muodot ja suuruusluokat vastasivat pääosin toisiaan. Tulosten parantamiseksi tulee selvittää lisää alkuarvoja tutkittavasta järjestelmästä ja varmistua nyt saatujen jousiparametrien oikeellisuudesta.

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.

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.