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.

Development of beam and plate finite elements based on the absolute nodal coordinate formulation

The focus of this dissertation is to develop finite elements based on the absolute nodal coordinate formulation. The absolute nodal coordinate formulation is a nonlinear finite element formulation, which is introduced for special requirements in the field of flexible multibody dynamics. In this formulation, a special definition for the rotation of elements is employed to ensure the formulation will not suffer from singularities due to large rotations. The absolute nodal coordinate formulation can be used for analyzing the dynamics of beam, plate and shell type structures. The improvements of the formulation are mainly concentrated towards the description of transverse shear deformation. Additionally, the formulation is verified by using conventional iso-parametric solid finite element and geometrically exact beam theory. Previous claims about especially high eigenfrequencies are studied by introducing beam elements based on the absolute nodal coordinate formulation in the framework of the large rotation vector approach. Additionally, the same high eigenfrequency problem is studied by using constraints for transverse deformation. It was determined that the improvements for shear deformation in the transverse direction lead to clear improvements in computational efficiency. This was especially true when comparative stress must be defined, for example when using elasto-plastic material. Furthermore, the developed plate element can be used to avoid certain numerical problems, such as shear and curvature lockings. In addition, it was shown that when compared to conventional solid elements, or elements based on nonlinear beam theory, elements based on the absolute nodal coordinate formulation do not lead to an especially stiff system for the equations of motion.

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.

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.

Dynamic Model of a Small Once-Through Boiler

In this study, a model for the unsteady dynamic behaviour of a once-through counter flow boiler that uses an organic working fluid is presented. The boiler is a compact waste-heat boiler without a furnace and it has a preheater, a vaporiser and a superheater. The relative lengths of the boiler parts vary with the operating conditions since they are all parts of a single tube. The present research is a part of a study on the unsteady dynamics of an organic Rankine cycle power plant and it will be a part of a dynamic process model. The boiler model is presented using a selected example case that uses toluene as the process fluid and flue gas from natural gas combustion as the heat source. The dynamic behaviour of the boiler means transition from the steady initial state towards another steady state that corresponds to the changed process conditions. The solution method chosen was to find such a pressure of the process fluid that the mass of the process fluid in the boiler equals the mass calculated using the mass flows into and out of the boiler during a time step, using the finite difference method. A special method of fast calculation of the thermal properties has been used, because most of the calculation time is spent in calculating the fluid properties. The boiler was divided into elements. The values of the thermodynamic properties and mass flows were calculated in the nodes that connect the elements. Dynamic behaviour was limited to the process fluid and tube wall, and the heat source was regarded as to be steady. The elements that connect the preheater to thevaporiser and the vaporiser to the superheater were treated in a special way that takes into account a flexible change from one part to the other. The model consists of the calculation of the steady state initial distribution of the variables in the nodes, and the calculation of these nodal values in a dynamic state. The initial state of the boiler was received from a steady process model that isnot a part of the boiler model. The known boundary values that may vary during the dynamic calculation were the inlet temperature and mass flow rates of both the heat source and the process fluid. A brief examination of the oscillation around a steady state, the so-called Ledinegg instability, was done. This examination showed that the pressure drop in the boiler is a third degree polynomial of the mass flow rate, and the stability criterion is a second degree polynomial of the enthalpy change in the preheater. The numerical examination showed that oscillations did not exist in the example case. The dynamic boiler model was analysed for linear and step changes of the entering fluid temperatures and flow rates.The problem for verifying the correctness of the achieved results was that there was no possibility o compare them with measurements. This is why the only way was to determine whether the obtained results were intuitively reasonable and the results changed logically when the boundary conditions were changed. The numerical stability was checked in a test run in which there was no change in input values. The differences compared with the initial values were so small that the effects of numerical oscillations were negligible. The heat source side tests showed that the model gives results that are logical in the directions of the changes, and the order of magnitude of the timescale of changes is also as expected. The results of the tests on the process fluid side showed that the model gives reasonable results both on the temperature changes that cause small alterations in the process state and on mass flow rate changes causing very great alterations. The test runs showed that the dynamic model has no problems in calculating cases in which temperature of the entering heat source suddenly goes below that of the tube wall or the process fluid.

Hydrostaattisen kuormituselementin parametrien vaikutus tehontarpeeseen

Työssä tutkitaan paperiteollisuuden käyttämän taipumakompensoidun telan hydrostaattisen kuormituselementin parametrien vaikutusta tehontarpeeseen ja kuormituselementin käyttäytymistä eri parametreilla. Työssä käsitellään hydrostaattisen kuormituselementin FE-laskennassa käytetyn mallinosien rakentaminen, kokoonpano ja teoria. Lisäksi analyyttinen laskenta esitetään lyhyesti. FE-laskennassa käytetään virtaus-rakenne -vuorovaikutusanalyysiä. Laskenta suoritetaan ANSYS/Flotran ohjelmistolla (versio 5.5). Hydrostaattisen kuormituselementin tehontarpeen vähentämiseen löytyi kolme tekijää. Suurimpana yksittäisenä tekijänä on elementin kitkapinta-alan pienentäminen, jolla saavutetaan 44 % kitkatehontarpeen pieneneminen. Seuraavaksi eniten vaikuttaa kitkatehontarpeeseen kuormituselementin muotoilu, jolla saavutetaan 20 % kitkatehontarpeen pieneneminen. Kolmantena tekijänä on viskositeetin pienentäminen, jolla suurimmillaan saavutetaan 15 % kitkatehontarpeen pieneneminen.

Ultralujasta teräksestä valmistetun henkilönostimen teleskooppipuomin kehittäminen

Tässä työssä kehitettiin palo- ja pelastuskäyttöön tarkoitettuun henkilönostimeen teleskooppipuomin profiilit. Profiilien valmistusmateriaalina oli kuumavalssattu, ultraluja säänkestävä rakenneteräs. Työssä kehitettiin standardien ja ohjeiden pohjalta laskentapohja, jolla voidaan tutkia teleskooppipuomin jaksojen tukireaktioita, taivutus- ja vääntömomentteja ja leikkaus ja normaalivoimia. Laskentapohjassa voidaan varioida eri kuormitusten suuntia, teleskooppipuomin sivusuuntaista ulottumaa ja nostokulmaa. Profiilien alustavassa mitoituksessa hyödynnettiin paikallisen lommahduksen huomioon ottavia standardeja ja suunnitteluohjeita. Eri poikkileikkausten ominaisuuksia verrattiin keskenään ja profiili valittiin yhdessä kohdeyrityksen kanssa. Alustavan mitoituksen yhteydessä muodostettiin apuohjelma valitulle poikkileikkaukselle, jolla voitiin tutkia profiilin eri muuttujien vaikutusta mm. paikalliseen lommahdukseen ja jäykkyyteen. Laskentapohjaan sisällytettiin myös optimointirutiini, jolla voitiin minimoida poikkileikkauksen pinta-ala ja tätä kautta profiilin massa. Lopullinen mitoitus suoritettiin elementtimenetelmällä. Mitoituksessa tutkittiin alustavasti mitoitettujen profiilien paikallista lommahdusta lineaarisen stabiilius- ja epälineaarisen analyysin pohjalta. Profiilien jännityksiä tutkittiin tarkemmin mm. varioimalla kuormituksia ja osittelemalla elementtien normaalijännityksiä. Diplomityössä kehitetyllä ja analysoidulla teleskooppipuomilla voitiin keventää jaksojen painoja 15-30 %. Sivusuuntainen ulottuma parani samalla lähes 20 % ja nimelliskuorma kasvoi 25 %.

Sekaelementit absoluuttisten solmukoordinaattien menetelmässä

Tässä lisensiaatintyössä käsitellään sekaelementtien sovellusmahdollisuuksia absoluuttisten solmukoordinaattien menetelmässä. Absoluuttisten solmukoordinaattien menetelmä on uudentyyppinen lähestymistapa elementtimenetelmän elementtien koordinaattien määrittämiseksi ja sen yhtenä tavoitteena on tehostaa suuria siirtymiä tai kiertymiä sisältävien elementtien laskentatehokkuutta. Tässä työssä absoluuttisten solmukoordinaattien menetelmä esitellään pääpiirteittäin sekä annetaan esimerkkejä muutamista tyypillisimmistä elementeistä lausuttuna edellä mainittujen koordinaattien perusteella. Sekaelementeiksi kutsutaan elementtityyppejä, missä tuntemattomien muuttujien joukkoja on aina enemmän kuin yksi. Sekaelementit erottavat redusoitumattomista elementeistä siirtymäkentän sisältyminen muuttujaryhmään ja hybridielementeistä muuttujien identtiset ulottuvuudet. Sekaelementtejä käytetään esimerkiksi kokoonpuristumattomien materiaalien rakenneanalyyseissä, alentamaan elementiltä vaadittavia jatkuvuusehtoja tai mallintamaan ilmiöitä, missä fysikaaliset ominaisuudet ovat jostain syystä voimakkaasti toisistaan riippuvaisia. Tämän lisensiaatintyön kirjoittamiseksi on tehty tutkimusta sekaelementtien mahdollisuuksista toimia absoluuttisten solmukoordinaattien menetelmässä. Tutkimuksen tuloksena on saatu aikaan kaksi tässä työssä esiteltävää, varsin rajatun toimintakyvyn omaavaa sekaelementtityyppiä, joiden siirtymäkentät on määritelty globaalien koordinaattien suhteen sisältäen myös orientaatiotermit. Tutkimusaihe vaatii kuitenkin vielä paljon lisätyötä, ennen kuin sekaelementtityyppejä voidaan kauttaaltaan soveltaa absoluuttisten solmukoordinaattien menetelmällä toteutetuissa rakenneanalyyseissä.

Sellumassan varastosäiliön rakenneanalyysi

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

Convenience of scallops in plate structures

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

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.