Numerical simulation of two-chamber plasma sources using finite element method

Numerical simulation of plasma sources is very important. Such models allows to vary different plasma parameters with high degree of accuracy. Moreover, they allow to conduct measurements not disturbing system balance.Recently, the scientific and practical interest increased in so-called two-chamber plasma sources. In one of them (small or discharge chamber) an external power source is embedded. In that chamber plasma forms. In another (large or diffusion chamber) plasma exists due to the transport of particles and energy through the boundary between chambers.In this particular work two-chamber plasma sources with argon and oxygen as active mediums were onstructed. This models give interesting results in electric field profiles and, as a consequence, in density profiles of charged particles.

Simulation and material calibration of ultra high strength steel (UHSS) S960 welded joint under static tensile test utilizing finite element method (FEM)

The aim of this work was to calibrate the material properties including strength and strain values for different material zones of ultra-high strength steel (UHSS) welded joints under monotonic static loading. The UHSS is heat sensitive and softens by heat due to welding, the affected zone is heat affected zone (HAZ). In this regard, cylindrical specimens were cut out from welded joints of Strenx® 960 MC and Strenx® Tube 960 MH, were examined by tensile test. The hardness values of specimens’ cross section were measured. Using correlations between hardness and strength, initial material properties were obtained. The same size specimen with different zones of material same as real specimen were created and defined in finite element method (FEM) software with commercial brand Abaqus 6.14-1. The loading and boundary conditions were defined considering tensile test values. Using initial material properties made of hardness-strength correlations (true stress-strain values) as Abaqus main input, FEM is utilized to simulate the tensile test process. By comparing FEM Abaqus results with measured results of tensile test, initial material properties will be revised and reused as software input to be fully calibrated in such a way that FEM results and tensile test results deviate minimum. Two type of different S960 were used including 960 MC plates, and structural hollow section 960 MH X-joint. The joint is welded by BöhlerTM X96 filler material. In welded joints, typically the following zones appear: Weld (WEL), Heat affected zone (HAZ) coarse grained (HCG) and fine grained (HFG), annealed zone, and base material (BaM). Results showed that: The HAZ zone is softened due to heat input while welding. For all the specimens, the softened zone’s strength is decreased and makes it a weakest zone where fracture happens while loading. Stress concentration of a notched specimen can represent the properties of notched zone. The load-displacement diagram from FEM modeling matches with the experiments by the calibrated material properties by compromising two correlations of hardness and strength.

Effective notch method assessment of misalignment for fatigue loaded welds

 The objective of this thesis work was to assess axial misalignment in fatigue loaded welds using the effective notch method. As a result, the fatigue behaviour of non-load carrying cruciform fillet welded joint under cyclic tensile loading has been studied. Various degrees of axial misalignment have been found in one series of non-load carrying cruciform fillet welded joints used in a laboratory investigation. As a result, it was important to carry out a comprehensive investigation since axial misalignment forms part of thequality of fatigue loaded structure and can reduce the fatigue strength. To extend the study, the correlation between fatigue strength and stress ratio, as well as stress concentration factor, were also studied. Moreover, a closer investigation of place of crack initiation and its dependence on weld sequence and imperfections of test specimen (angular distortion) was studied. For the fatigue class calculations, FEM (finite element method) and the effectivenotch approach are used. The addressed variable is the axial misalignment whichis introduce by modeling the entire joint. Fracture mechanics based calculations are also used and quantitatively compared with effective notch and experimental results.

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.

Analysis of Industrial Granular Flow Applications by Using Advanced Collision Models

Granular flow phenomena are frequently encountered in the design of process and industrial plants in the traditional fields of the chemical, nuclear and oil industries as well as in other activities such as food and materials handling. Multi-phase flow is one important branch of the granular flow. Granular materials have unusual kinds of behavior compared to normal materials, either solids or fluids. Although some of the characteristics are still not well-known yet, one thing is confirmed: the particle-particle interaction plays a key role in the dynamics of granular materials, especially for dense granular materials. At the beginning of this thesis, detailed illustration of developing two models for describing the interaction based on the results of finite-element simulation, dimension analysis and numerical simulation is presented. The first model is used to describing the normal collision of viscoelastic particles. Based on some existent models, more parameters are added to this model, which make the model predict the experimental results more accurately. The second model is used for oblique collision, which include the effects from tangential velocity, angular velocity and surface friction based on Coulomb's law. The theoretical predictions of this model are in agreement with those by finite-element simulation. I n the latter chapters of this thesis, the models are used to predict industrial granular flow and the agreement between the simulations and experiments also shows the validation of the new model. The first case presents the simulation of granular flow passing over a circular obstacle. The simulations successfully predict the existence of a parabolic steady layer and show how the characteristics of the particles, such as coefficients of restitution and surface friction affect the separation results. The second case is a spinning container filled with granular material. Employing the previous models, the simulation could also reproduce experimentally observed phenomena, such as a depression in the center of a high frequency rotation. The third application is about gas-solid mixed flow in a vertically vibrated device. Gas phase motion is added to coherence with the particle motion. The governing equations of the gas phase are solved by using the Large eddy simulation (LES) and particle motion is predicted by using the Lagrangian method. The simulation predicted some pattern formation reported by experiment.

The influence ofvariable amplitude loading on the fatigue strength of rhs corners

Rectangular hollow section (RHS) members are components widely used in engineering applications because of their good-looking, good properties in engineering areas and inexpensive cost comparing to members with other sections. The increasing use of RHS in load bearing structures makes it necessary to analyze the fatigue behavior of the RHS members. In this thesis, concentration will be given to the fatigue behavior of the RHS members under variable amplitude pure torsional loading. For the RHS members, failure will normally occur in the corner region if the welded regions are under full penetration. This is because of the complicated stress components' distributions at the RHScorners, where all of three fracture mechanics modes will happen. Mode I is mainly caused by the residual stresses that caused by the manufacturing process. Modes II and III are caused by the applied torsional loading. Stress based Findleymodel is also used to analyze the stress components. Constant amplitude fatigue tests have been done as well as variable amplitude fatigue tests. The specimens under variable amplitude loading gave longer fatigue lives than those under constant amplitude loading. Results from tests show an S-N curvewith slope around 5.

Fatigue Analysis of Welded Cruciform Joints with Cold Laps

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

Napakäämityn murtovakokestomagneettitahtikoneen suunnittelu

Diplomityössä tutustutaan napakäämittyihin murtovakokestomagneettitahtikoneisiin sekä syvennytään niiden suunnitteluun analyyttisen laskennanavulla. Analyyttisen laskennan haasteena on perinteisesti käytettyjen kestomagneettitahtikoneiden laskentayhtälöiden sovittaminen murtovakokäämityille kestomagneettitahtikoneille. Työn lopussa analyyttisen laskennan tuloksia verrataan kenttälaskentaohjelmalla saatuihin tuloksiin. Tärkeimmät analyyttiset yhtälöt sekä laskentamallin rakenne on esitetty työssä.

Rakenteellisen jouston kuvaus reaaliaikasimuloinnissa

Työn tavoitteena oli tuottaa rakenteellisen jouston huomioiva monikappaledynmiikan simulointiohjelma Matlab-ympäristöön. Rakenteellinen jousto huomioitiin kelluvan koordinaatiston menetelmällä ja joustavuutta kuvaavat muodot ratkaistiin elementtimenetelmällä. Tehdyn ohjelman avulla voidaan koostaa joustavista kappaleista koostuvia avaruusmekanismeja ja tutkia niiden dynaamista käyttäytymistä. Simulointitulosta verrattiin kaupallisen ohjelmiston tuottamaan tulokseen. Työssä havaittiin, että kelluvan koordinaatiston menetelmä on käyttökelpoinen reaaliaikaiseen simulointiin. Työssä toteutetun ohjelman tulokset vastasivat kaupallisen simulointiohjelman tuloksia.

Deep-hole drilling process, monitoring system and sensors modifications

 In the drilling processes and especially deep-hole drilling process, the monitoring system and having control on mechanical parameters (e.g. Force, Torque,Vibration and Acoustic emission) are essential. The main focus of this thesis work is to study the characteristics of deep-hole drilling process, and optimize the monitoring system for controlling the process. The vibration is considered as a major defect area of the deep-hole drilling process which often leads to breakage of the drill, therefore by vibration analysis and optimizing the workpiecefixture, this area is studied by finite element method and the suggestions are explained. By study on a present monitoring system, and searching on the new sensor products, the modifications and recommendations are suggested for optimize the present monitoring system for excellent performance in deep-hole drilling process research and measurements.

Terässiilojen mitoitusohjelman kehittäminen

Siiloja käytetään jauhemaisten ja rakeisten aineiden varastointiin ja annosteluun maataloudessa sekä mm. lasi-, laasti- ja tasoiteteollisuudessa. Tällaiset laitokset toimitetaan usein kokonaistoimituksina, joissa tavanomaisten teräsrakenteiden kustannusvaikutus voi olla jopa 30 % koko toimituksen arvosta. Niinpä osaltaan siilojen tehokkaalla mitoituksella ja erilaisiarakenneratkaisuja vertailemalla voidaan alentaa toimitusten kokonaiskustannuksia. Tässä diplomityössä perehdyttiin siilojen rakenteisiin ja niiden lujuustekniseen mitoitukseen etenkin tuentaratkaisujen osalta. Mitoituksen osalta työssä perehdyttiin pääasiassa uusiin Eurooppalaisiin standardeihin SFS-EN 1991-4 ja SFS-EN 1993-4-1. Niiden mukaan määritetään mitoituskuormat sekä tehdään lujuustekninen mitoitus. Työn tuloksena laadittiin Microsoft Excel -pohjainen mitoitusohjelma, jolla voidaan nopeasti mitoittaa siilojen olennaisimmat teräsrakenteet. Ohjelmalla voidaan määrittää mitoitus koskien siilon lieriön ja kartion seinämiä sekä tuentaratkaisuja. Ohjelma nopeuttaa erityisesti tarjousvaiheen suunnittelua ja antaa pohjan varsinaiselle siilojen yksityiskohtaiselle mitoittamiselle ja suunnittelulle.

Dynamic Simulations of Rotors during Drop on Retainer Bearings

The active magnetic bearings have recently been intensively developed because of noncontact support having several advantages compared to conventional bearings. Due to improved materials, strategies of control, and electrical components, the performance and reliability of the active magnetic bearings are improving. However, additional bearings, retainer bearings, still have a vital role in the applications of the active magnetic bearings. The most crucial moment when the retainer bearings are needed is when the rotor drops from the active magnetic bearings on the retainer bearings due to component or power failure. Without appropriate knowledge of the retainer bearings, there is a chance that an active magnetic bearing supported rotor system will be fatal in a drop-down situation. This study introduces a detailed simulation model of a rotor system in order to describe a rotor drop-down situation on the retainer bearings. The introduced simulation model couples a finite element model with component mode synthesis and detailed bearing models. In this study, electrical components and electromechanical forces are not in the focus. The research looks at the theoretical background of the finite element method with component mode synthesis that can be used in the dynamic analysis of flexible rotors. The retainer bearings are described by using two ball bearing models, which include damping and stiffness properties, oil film, inertia of rolling elements and friction between races and rolling elements. Thefirst bearing model assumes that the cage of the bearing is ideal and that the cage holds the balls in their predefined positions precisely. The second bearing model is an extension of the first model and describes the behavior of the cageless bearing. In the bearing model, each ball is described by using two degrees of freedom. The models introduced in this study are verified with a corresponding actual structure. By using verified bearing models, the effects of the parameters of the rotor system onits dynamics during emergency stops are examined. As shown in this study, the misalignment of the retainer bearings has a significant influence on the behavior of the rotor system in a drop-down situation. In this study, a stability map of the rotor system as a function of rotational speed of the rotor and the misalignment of the retainer bearings is presented. In addition, the effects of parameters of the simulation procedure and the rotor system on the dynamics of system are studied.

A novel procedure for establishing clinching parameters for high strength steel sheet

The building industry has a particular interest in using clinching as a joining method for frame constructions of light-frame housing. Normally many clinch joints are required in joining of frames.In order to maximise the strength of the complete assembly, each clinch joint must be as sound as possible. Experimental testing is the main means of optimising a particular clinch joint. This includes shear strength testing and visual observation of joint cross-sections. The manufacturers of clinching equipment normally perform such experimental trials. Finite element analysis can also be used to optimise the tool geometry and the process parameter, X, which represents the thickness of the base of the joint. However, such procedures require dedicated software, a skilled operator, and test specimens in order to verify the finite element model. In addition, when using current technology several hours' computing time may be necessary. The objective of the study was to develop a simple calculation procedure for rapidly establishing an optimum value for the parameter X for a given tool combination. It should be possible to use the procedure on a daily basis, without stringent demands on the skill of the operator or the equipment. It is also desirable that the procedure would significantly decrease thenumber of shear strength tests required for verification. The experimental workinvolved tests in order to obtain an understanding of the behaviour of the sheets during clinching. The most notable observation concerned the stage of the process in which the upper sheet was initially bent, after which the deformation mechanism changed to shearing and elongation. The amount of deformation was measured relative to the original location of the upper sheet, and characterised as the C-measure. By understanding in detail the behaviour of the upper sheet, it waspossible to estimate a bending line function for the surface of the upper sheet. A procedure was developed, which makes it possible to estimate the process parameter X for each tool combination with a fixed die. The procedure is based on equating the volume of material on the punch side with the volume of the die. Detailed information concerning the behaviour of material on the punch side is required, assuming that the volume of die does not change during the process. The procedure was applied to shear strength testing of a sample material. The sample material was continuously hot-dip zinc-coated high-strength constructional steel,with a nominal thickness of 1.0 mm. The minimum Rp0.2 proof stress was 637 N/mm2. Such material has not yet been used extensively in light-frame housing, and little has been published on clinching of the material. The performance of the material is therefore of particular interest. Companies that use clinching on a daily basis stand to gain the greatest benefit from the procedure. By understanding the behaviour of sheets in different cases, it is possible to use data at an early stage for adjusting and optimising the process. In particular, the functionality of common tools can be increased since it is possible to characterise the complete range of existing tools. The study increases and broadens the amount ofbasic information concerning the clinching process. New approaches and points of view are presented and used for generating new knowledge.

Direct-on-line axial flux permanent magnet synchronous generator static and dynamic performance

In distributed energy production, permanent magnet synchronous generators (PMSG) are often connected to the grid via frequency converters, such as voltage source line converters. The price of the converter may constitute a large part of the costs of a generating set. Some of the permanent magnet synchronous generators with converters and traditional separately excited synchronous generators couldbe replaced by direct-on-line (DOL) non-controlled PMSGs. Small directly networkconnected generators are likely to have large markets in the area of distributed electric energy generation. Typical prime movers could be windmills, watermills and internal combustion engines. DOL PMSGs could also be applied in island networks, such as ships and oil platforms. Also various back-up power generating systems could be carried out with DOL PMSGs. The benefits would be a lower priceof the generating set and the robustness and easy use of the system. The performance of DOL PMSGs is analyzed. The electricity distribution companies have regulations that constrain the design of the generators being connected to the grid. The general guidelines and recommendations are applied in the analysis. By analyzing the results produced by the simulation model for the permanent magnet machine, the guidelines for efficient damper winding parameters for DOL PMSGs are presented. The simulation model is used to simulate grid connections and load transients. The damper winding parameters are calculated by the finite element method (FEM) and determined from experimental measurements. Three-dimensional finite element analysis (3D FEA) is carried out. The results from the simulation model and 3D FEA are compared with practical measurements from two prototype axial flux permanent magnet generators provided with damper windings. The dimensioning of the damper winding parameters is case specific. The damper winding should be dimensioned based on the moment of inertia of the generating set. It is shown that the damper winding has optimal values to reach synchronous operation in the shortest period of time after transient operation. With optimal dimensioning, interferenceon the grid is minimized.