Simulation of structural stress history based on dynamic analysis

Modern machine structures are often fabricated by welding. From a fatigue point of view, the structural details and especially, the welded details are the most prone to fatigue damage and failure. Design against fatigue requires information on the fatigue resistance of a structure’s critical details and the stress loads that act on each detail. Even though, dynamic simulation of flexible bodies is already current method for analyzing structures, obtaining the stress history of a structural detail during dynamic simulation is a challenging task; especially when the detail has a complex geometry. In particular, analyzing the stress history of every structural detail within a single finite element model can be overwhelming since the amount of nodal degrees of freedom needed in the model may require an impractical amount of computational effort. The purpose of computer simulation is to reduce amount of prototypes and speed up the product development process. Also, to take operator influence into account, real time models, i.e. simplified and computationally efficient models are required. This in turn, requires stress computation to be efficient if it will be performed during dynamic simulation. The research looks back at the theoretical background of multibody dynamic simulation and finite element method to find suitable parts to form a new approach for efficient stress calculation. This study proposes that, the problem of stress calculation during dynamic simulation can be greatly simplified by using a combination of floating frame of reference formulation with modal superposition and a sub-modeling approach. In practice, the proposed approach can be used to efficiently generate the relevant fatigue assessment stress history for a structural detail during or after dynamic simulation. In this work numerical examples are presented to demonstrate the proposed approach in practice. The results show that approach is applicable and can be used as proposed.

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.

Non-linear analysis of laminated metal matrix composites by an integrated micro/macro-mechanical model

The present paper describes an integrated micro/macro mechanical study of the elastic-viscoplastic behavior of unidirectional metal matrix composites (MMC). The micromechanical analysis of the elastic moduli is based on the Composites Cylinder Assemblage model (CCA) with comparisons also draw with a Representative Unit Cell (RUC) technique. These "homogenization" techniques are later incorporated into the Vanishing Fiber Diameter (VFD) model and a new formulation is proposed. The concept of a smeared element procedure is employed in conjunction with two different versions of the Bodner and Partom elastic-viscoplastic constitutive model for the associated macroscopic analysis. The formulations developed are also compared against experimental and analytical results available in the literature.

Object-oriented Programming Applied to the Development of Structural Analysis and Optimization Software

This paper presents the development of a two-dimensional interactive software environment for structural analysis and optimization based on object-oriented programming using the C++ language. The main feature of the software is the effective integration of several computational tools into graphical user interfaces implemented in the Windows-98 and Windows-NT operating systems. The interfaces simplify data specification in the simulation and optimization of two-dimensional linear elastic problems. NURBS have been used in the software modules to represent geometric and graphical data. Extensions to the analysis of three-dimensional problems have been implemented and are also discussed in this paper.

Numerical Simulation of the Cavitation in the Hydrodynamic Lubrication of Journal Bearings: a Parallel Algorithm

In this paper we present an algorithm for the numerical simulation of the cavitation in the hydrodynamic lubrication of journal bearings. Despite the fact that this physical process is usually modelled as a free boundary problem, we adopted the equivalent variational inequality formulation. We propose a two-level iterative algorithm, where the outer iteration is associated to the penalty method, used to transform the variational inequality into a variational equation, and the inner iteration is associated to the conjugate gradient method, used to solve the linear system generated by applying the finite element method to the variational equation. This inner part was implemented using the element by element strategy, which is easily parallelized. We analyse the behavior of two physical parameters and discuss some numerical results. Also, we analyse some results related to the performance of a parallel implementation of the algorithm.

A Design Methodology for the Compensation of Positioning Deviation in Gantry Manipulators

This work presents recent results concerning a design methodology used to estimate the positioning deviation for a gantry (Cartesian) manipulator, related mainly to structural elastic deformation of components during operational conditions. The case-study manipulator is classified as gantry type and its basic dimensions are 1,53m x 0,97m x 1,38m. The dimensions used for the calculation of effective workspace due to end-effector path displacement are: 1m x 0,5m x 0,5m. The manipulator is composed by four basic modules defined as module X, module Y, module Z and terminal arm, where is connected the end-effector. Each module controlled axis performs a linear-parabolic positioning movement. The planning path algorithm has the maximum velocity and the total distance as input parameters for a given task. The acceleration and deceleration times are the same. Denavit-Hartemberg parameterization method is used in the manipulator kinematics model. The gantry manipulator can be modeled as four rigid bodies with three degrees-of-freedom in translational movements, connected as an open kinematics chain. Dynamic analysis were performed considering inertial parameters specification such as component mass, inertia and center of gravity position of each module. These parameters are essential for a correct manipulator dynamic modelling, due to multiple possibilities of motion and manipulation of objects with different masses. The dynamic analysis consists of a mathematical modelling of the static and dynamic interactions among the modules. The computation of the structural deformations uses the finite element method (FEM).

Voimalaitoksen kanavien kannakkeiden rakennesuunnittelu

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

Tehollisen lovijännityksen menetelmän käytettävyys ultralujien terästen korkealaatuisten hitsien väsymismitoituksessa

Tässä diplomityössä on tutkittu tehollisen lovijännityksen menetelmän soveltuvuutta ultralujien terästen korkealaatuisten hitsien väsymismitoitukseen. International Institute of Welding suosittelee käyttämään elementtimenetelmässä hitsin rajaviivoilla sekä juuressa fiktiivistä 1 mm pyöristystä, jonka avulla tehollinen lovijännitys määritetään. Kaikille liitostyypeille sovelletaan samaa, kaltevuudeltaan m = 3 olevaa SN-käyrää, jolloin maksimipääjännitystä vastaava väsymisluokka FAT saa arvon 225. Nykyisiä mitoitusohjeita on pidetty kuitenkin liian konservatiivisina, etenkin jos kyseessä on suurilujuuksisesta teräksestä valmistettu korkealaatuinen hitsi. Rajaviivalla vaikuttavaa lovijännitystä on tutkittu mallintamalla liitokset FEMAP – elementtimenetelmäohjelmalla varioimalla rajaviivan pyöristystä. Elementtimenetelmän tuloksia on verrattu analyyttisiin loven muotoluvun laskentakaavoihin. Tutkittavana on ollut Ruukin Optim 960 QC sekä Optim 1100 QC – teräksistä valmistettuja koesauvoja. Koesauvat on valmistettu sekä koestettu pääasiassa Lappeenrannan teknillisen yliopiston teräsrakenteiden laboratoriossa. Tutkittavat koesauvat ovat olleet kuormaa kantamattomia ristiliitoksia sekä päittäisliitoksia. Suurin osa koesauvoista on väsytetty käyttämällä jännityssuhdetta R < 0,11. Koesauvat on jaoteltu jännityssuhteen sekä liitostyypin mukaan. Kaikkien koekappaleiden karakteristiseksi väsymisluokan arvoksi on määritetty FAT 200. Alle 0,11 jännityssuhteella väsytettyjen koekappaleiden karakteristinen väsymisluokka on FAT 230 ja isoilla jännityssuhteilla väsytettyjen FAT 126. Tulosten perusteella nykyiset mitoitusohjeet eivät ole liian konservatiivisia. Väsymisluokkaa FAT 225 voidaan käyttää väsymislaskennassa, mikäli rakenteen kuormitusten suhde on alle 0,1. Isoilla jännityssuhteilla koestettujen koekappaleiden lukumäärä on ollut pieni, joten niiden mitoitukselle ei voida antaa tarkkoja ohjeita.

Turveperävaunun väsymiskestävyyden parantaminen

Työssä tutkittiin turveperävaunun renkaan kiinnitysrakenteen väsymiskestävyyttä. Väsymiskestävyyden parantamiseksi renkaiden kiinnityskohdan muotoilu ja mitoitus suunniteltiin uudelleen. Suunnittelussa keskityttiin ensisijaisesti nimellisten ja rakenteellisten jännitysheilahdusten pienentämiseen. Sekä vanhan, että uuden rakenneratkaisun väsymiskestoikiä tarkasteltiin käyttämällä hot spot – menetelmää.

Value engineering of InnoTrackTM

Middle section module of InnoTrackTM moving walk was re-engineered according to value analysis process. Self-supporting steel structure for moving walk was created as a result of this process. Designed structure was verified and validated by prototype tests and finite element method calculations. Self-supporting steel structure replaces the original design of middle section module in InnoTrackTM. Designed structure provides higher satisfaction to customers’ needs and at the same time, it uses less resources. The redesigned middle section module provides higher value to the customer.

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

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

Vibrations in Rotating Machinery Arising from Minor Imperfections in Component Geometries

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

Trukkitraktorin rungon suunnittelu ja laskenta

Tässä työssä suunniteltiin lappeenrantalaisen Astex Oy:n tilauksesta liikkuvan työkoneen runkorakenne. Työ tehtiin Lappeenrannan teknillisen yliopiston teräsrakenteiden laboratoriossa. Suunniteltava rakenne kuului linkkuohjattuun, noin 5000 kg painoiseen trukkitraktoriin. Lähtökohtana rakenteen suunnittelulle olivat tilaajan asettamat rakenteen geometriaan ja suorituskykyyn liittyvät reunaehdot ja rajoitteet. Uuden rakenteen suunnittelussa hyödynnettiin myös tilaajatahon kehittelemää vastaavan tyyppistä prototyyppirakennetta. Rakenteen suunnitteluprosessi koostui neljästä eri vaiheesta. Suunniteltavan rakenteen lähtökohtana olleelle prototyyppirakenteelle suoritettiin koneen käytönaikaisia venymäliuskamittauksia rakenteen kuormitushistorian selvittämiseksi. Mittauksista saatujen tulosten perusteella määritettiin kuormitukset uudelle suunniteltavalle rakenteelle. Määritettyjä kuormituksia hyödyntäen ideoitiin, suunniteltiin ja mallinnettiin uusi tilaajan vaatimuksia vastaava rakenne. Uudelle rakenteelle suoritettiin lujuusanalyysit FE-analyysiä hyödyntäen. Uuden rakenteen suunnittelussa kiinnitettiin huomiota rakenteen hyvään valmistettavuuteen ja suunniteltiin rakenneratkaisut tilaajataholle mahdollisimman optimaalisiksi. Suunnittelu- ja mallinnustyö tehtiin Solidworks 2014 ohjelmistolla. Rakenteen lujuustekniset tarkastelut sisälsivät rakennedetaljien analyyttistä mitoitusta ja laskentaa. FE-laskennalla selvitettiin rakenteen ääri- ja väsymiskestävyys. Laskenta sisälsi koko rakenteen globaaleja tarkasteluja, sekä eri kriittisten rakennedetaljien yksityiskohtaisempia analyysejä. FE-laskennan pääpaino oli rakenteen väsymisanalyyseissä, jotka toteutettiin Hot-Spot- ja särönkasvumenetelmillä. Rakenteen FE-analyysien suorittamisessa käytettiin Femap, NxNastran ja Abaqus-ohjelmistoja.