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.

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.

Development of finite elements for analysis of biomechanical structures using flexible multibody formulations

The absolute nodal coordinate formulation was originally developed for the analysis of structures undergoing large rotations and deformations. This dissertation proposes several enhancements to the absolute nodal coordinate formulation based finite beam and plate elements. The main scientific contribution of this thesis relies on the development of elements based on the absolute nodal coordinate formulation that do not suffer from commonly known numerical locking phenomena. These elements can be used in the future in a number of practical applications, for example, analysis of biomechanical soft tissues. This study presents several higher-order Euler–Bernoulli beam elements, a simple method to alleviate Poisson’s and transverse shear locking in gradient deficient plate elements, and a nearly locking free gradient deficient plate element. The absolute nodal coordinate formulation based gradient deficient plate elements developed in this dissertation describe most of the common numerical locking phenomena encountered in the formulation of a continuum mechanics based description of elastic energy. Thus, with these fairly straightforwardly formulated elements that are comprised only of the position and transverse direction gradient degrees of freedom, the pathologies and remedies for the numerical locking phenomena are presented in a clear and understandable manner. The analysis of the Euler–Bernoulli beam elements developed in this study show that the choice of higher gradient degrees of freedom as nodal degrees of freedom leads to a smoother strain field. This improves the rate of convergence.

Hygro-elasto-plastic behavior of planar orthotropic material

The mechanical and hygroscopic properties of paper and board are factors affecting the whole lifecycle of a product, including paper/board quality, production, converting, and material and energy savings. The progress of shrinkage profiles, loose edges of web, baggy web causing wrinkling and misregistration in printing are examples of factors affecting runnability and end product quality in the drying section and converting processes, where paper or board is treated as a moving web. The structural properties and internal stresses or plastic strain differences built up during production also cause the end-product defects related to distortion of the shape of the product such as sheet or box. The objective of this work was to construct a model capable of capturing the characteristic behavior of hygroscopic orthotropic material under moisture change, during different external in-plane stretch or stress conditions. Two independent experimental models were constructed: the elasto-plastic material model and the hygroexpansivity-shrinkage model. Both describe the structural properties of the sheet with a fiber orientation probability distribution, and both are functions of the dry solids content and fiber orientation anisotropy index. The anisotropy index, introduced in this work, simplifies the procedure of determining the constitutive parameters of the material model and the hygroexpansion coefficients in different in-plane directions of the orthotropic sheet. The mathematically consistent elasto-plastic material model and the dry solids content dependent hygroexpansivity have been constructed over the entire range from wet to dry. The presented elastoplastic and hygroexpansivity-shrinkage models can be used in an analytical approach to estimate the plastic strain and shrinkage in simple one-dimensional cases. For studies of the combined and more complicated effects of hygro-elasto-plastic behavior, both models were implemented in a finite element program for a numerical solution. The finite element approach also offered possibilities for studying different structural variations of orthotropic planar material, as well as local buckling behavior and internal stress situations of the sheet or web generated by local strain differences. A comparison of the simulation examples presented in this work to results published earlier confirms that the hygro-elasto-plastic model provides at least qualitatively reasonable estimates. The application potential of the hygro-elasto-plastic model is versatile, including several phenomena and defects appearing in the drying, converting and end-use conditions of the paper or board webs and products, or in other corresponding complex planar materials.

Inductance Calculation of Tooth-Coil Permanent-Magnet Synchronous Machines

Analytical calculation methods for all the major components of the synchronous inductance of tooth-coil permanentmagnet synchronous machines are reevaluated in this paper. The inductance estimation is different in the tooth-coil machine compared with the one in the traditional rotating field winding machine. The accuracy of the analytical torque calculation highly depends on the estimated synchronous inductance. Despite powerful finite element method (FEM) tools, an accurate and fast analytical method is required at an early design stage to find an initialmachine design structure with the desired performance. The results of the analytical inductance calculation are verified and assessed in terms of accuracy with the FEM simulation results and with the prototype measurement results.

The Fluid Structure Interaction Analysis of a Peristaltic Pump

The thesis work models the squeezing of the tube and computes the fluid motion of a peristaltic pump. The simulations have been conducted by using COMSOL Multiphysics FSI module. The model is setup in axis symmetric with several simulation cases to have a clear understanding of the results. The model captures total displacement of the tube, velocity magnitude, and average pressure fluctuation of the fluid motion. A clear understanding and review of many mathematical and physical concepts are also discussed with their applications in real field. In order to solve the problems and work around the resource constraints, a thorough understanding of mass balance and momentum equations, finite element concepts, arbitrary Lagrangian-Eulerian method, one-way coupling method, two-way coupling method, and COMSOL Multiphysics simulation setup are understood and briefly narrated.

Sähkökoneessa esiintyvien voimien yksityiskohtainen mittausjärjestelmä

Hissiteollisuudessa nostokoneistoina käytettyjen sähkömoottoreiden laatuvaatimukset ovat tiuken-tuneet viime vuosina. Erityisesti koneistojen tuottama ääni ja mekaaninen värähtely ovat olleet jat-kuvasti tiukentuneen tarkastelun alaisena. Hissikoriin ja hissiä ympäröiviin rakenteisiin välittyvästä värähtelystä johtuva ääni on yksi hissin laatuvaikutelmaan merkittävimmin vaikuttavia tekijöitä. Nostokoneisto on yksi tärkeimmistä äänen ja värähtelyn lähteistä hissijärjestelmässä. Koneiston suunnittelulla edellä mainittuja tekijöitä voidaan minimoida. Sähkökoneiden suunnittelussa finiit-tielementtimenetelmien (FEM) käyttö on vakiintunut haastavimmissa sovelluksissa. Kone Oyj:llä nostokoneistoina käytetään aksiaalivuokestomagneettitahtikoneita (AFPMSM), joiden FEM simu-lointiin käytetään yleisesti kolmea eri tapaa. Kukin näistä vaihtoehdoista pitää sisällään omat hyö-tynsä, että haittansa. Suunnittelun kannalta tärkeää on oikean menetelmän valinta ai-ka/informatiivisuus suhteen maksimoimiseksi. Erittäin tärkeää on myös saatujen tulosten oikeelli-suus. Tämän diplomityön tavoite on kehittää järjestelmä, jonka avulla AFPMS-koneen voimia voidaan mitata yksityiskohtaisella tasolla. Järjestelmän avulla voidaan tarkastella käytössä olevien FE-menetelmien tulosten oikeellisuutta sekä äänen että värähtelyn syntymekanismeja. Järjestelmän tarkoitus on myös syventää Kone Oyj tietotaitoa AFPMS-koneiden toiminnasta. Tässä työssä esitellään AFPMS-koneen epäideaalisuuksia, jotka voivat vaikuttaa mittajärjestelmän suunnitteluun. Myös koneen epäideaalisuuksiin lukeutuvaa ääntä on tarkasteltu tässä työssä. Jotta työn tavoitteiden mukaista FE-menetelmien vertailua ja tulosten oikeellisuuden tarkastelua voitai-siin tehdä, myös yleisimpiä AFPMS-koneen FE-menetelmiä tarkastellaan. Työn tuloksena on mittajärjestelmän suunnitelma, jonka avulla voidaan toteuttaa kuuden vapausas-teen voimamittaus jokaiselle koneistomagneetille alle 1N resoluutiolla. Suunnitellun järjestelmän toimivuutta on tarkasteltu FE-menetelmiä käyttäen ja järjestelmässä käytettävän voima-anturin ky-vykkyyttä on todennettu referenssimittauksin. Suunniteltu mittajärjestelmä mahdollistaa sähkömoottorin useiden eri epäideaalisuuksien tarkaste-lun yksityiskohtaisella tasolla. Mittausajatuksen soveltaminen myös muiden koneiden tutkimiseen tarjoaa mahdollisuuksia jatkotutkimuksille.

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.

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.

Virtuaaliprototyypin käyttö roottorin dynamiikan analysoinnissa

Työn tavoitteena oli selvittää kaupallisen dynamiikansimulointiohjelmiston so-veltuvuus roottoridynamiikan analysointiin. Työssä keskityttiin erityisesti rootto-rin dynamiikkaan vaikuttavien epäideaalisuuksien mallintamiseen. Simulointitu-losten tarkkuutta selvitettiin mittauksilla. Lisäksi vertailtiin yleiskäyttöisen dyna-miikan simulointiohjelmiston ja roottoridynamiikan erikoisohjelmiston teoriaa. Tutkittava roottori oli paperikoneen putkitela. Telan joustavuus kuvattiin ele-menttimenetelmällä ratkaistujen moodien avulla. Elementtimallissa huomioitiin telan vaipan seinämänpaksuusvaihtelu, joka vaikuttaa telan massa- ja jäykkyysja-kaumaan. Dynamiikkaohjelmistossa mallinnettiin telan tuennasta tulevat herätteet. Dynamiikkaohjelmistona käytettiin ADAMS:ia ja FEM-ohjelmana ANSYS:stä. Tuloksista havaittiin käytetyn menetelmän soveltuvan roottoridynamiikan ana-lysointiin ja roottorin epäideaalisuuksien mallintamiseen. Simulointimallilla saa-tiin esille murtolukukriittiset pyörimisnopeudet ja telan kriittinen pyörimisnopeus vastasi hyvin mittaustuloksia.

Field Equilibrium Finite Elements for Structural Analysis

Many finite elements used in structural analysis possess deficiencies like shear locking, incompressibility locking, poor stress predictions within the element domain, violent stress oscillation, poor convergence etc. An approach that can probably overcome many of these problems would be to consider elements in which the assumed displacement functions satisfy the equations of stress field equilibrium. In this method, the finite element will not only have nodal equilibrium of forces, but also have inner stress field equilibrium. The displacement interpolation functions inside each individual element are truncated polynomial solutions of differential equations. Such elements are likely to give better solutions than the existing elements.In this thesis, a new family of finite elements in which the assumed displacement function satisfies the differential equations of stress field equilibrium is proposed. A general procedure for constructing the displacement functions and use of these functions in the generation of elemental stiffness matrices has been developed. The approach to develop field equilibrium elements is quite general and various elements to analyse different types of structures can be formulated from corresponding stress field equilibrium equations. Using this procedure, a nine node quadrilateral element SFCNQ for plane stress analysis, a sixteen node solid element SFCSS for three dimensional stress analysis and a four node quadrilateral element SFCFP for plate bending problems have been formulated.For implementing these elements, computer programs based on modular concepts have been developed. Numerical investigations on the performance of these elements have been carried out through standard test problems for validation purpose. Comparisons involving theoretical closed form solutions as well as results obtained with existing finite elements have also been made. It is found that the new elements perform well in all the situations considered. Solutions in all the cases converge correctly to the exact values. In many cases, convergence is faster when compared with other existing finite elements. The behaviour of field consistent elements would definitely generate a great deal of interest amongst the users of the finite elements.

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.

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.