Investigation of Two Bridge Alternatives for Low Volume Roads Volume 2 of 2 Concept 2: Beam In Slab Bridge, HR-382, 1997

Recent reports have indicated that 23.5% of the nation's highway bridges are structurally deficient and 17.7% are functionally obsolete. A significant number of these bridges are on the Iowa secondary road system where over 86% of the rural bridge management responsibilities are assigned to the counties. Some of the bridges can be strengthened or otherwise rehabilitated, but many more are in need of immediate replacement. In a recent investigation (HR-365 "Evaluation of Bridge Replacement Alternatives for the County Bridge System") several types of replacement bridges that are currently being used on low volume roads were identified. It was also determined that a large number of counties (69%) have the ability and are interested in utilizing their own forces to design and construct short span bridges. After reviewing the results from HR-365, the research team developed one "new" bridge replacement concept and a modification of a replacement system currently being used. Both of these bridge replacement alternatives were investigated in this study, the results of which are presented in two volumes. This volume (Volume 2) presents the results of Concept 2 - Modification of the Beam-in-Slab Bridge. Concept 1 - Steel Beam Precast Units is presented in Volume 1. Concept 2 involves various laboratory tests of the Beam-in-Slab bridge (BISB) currently being used by Benton County and several other Iowa counties. In this investigation, the behavior and strength of the BISB were determined; a new method of obtaining composite action between the steel beams and concrete was also tested. Since the Concept 2 bridge is primarily intended for use on low-volume roads, the system can be constructed with new or used beams. In the experimental part of the investigation, there were three types of laboratory tests: push-out tests, service and ultimate load tests of models of the BISB, and composite beam tests utilizing the newly developed shear connection. In addition to the laboratory tests, there was a field test in which an existing BISB was service load tested. An equation was developed for predicting the strength of the shear connection investigated; in addition, a finite element model for analyzing the BISB was also developed. Push-out tests were completed to determine the strength of the recently developed shear connector. A total of 36 specimens were tested, with variables such as hole diameter, hole spacing, presence of reinforcement, etc. being investigated. In the model tests of the BISB, two and four beam specimens [L=9,140 mm (30 ft)] were service load tested for behavior and load distribution data. Upon completion of these tests, both specimens were loaded to failure. In the composite beam tests, four beams, one with standard shear studs and three using the shear connection developed, were tested. Upon completion of the service load tests, all four beams were loaded to failure. The strength and behavior of the beams with the new shear connection were found to be essentially the same as that of the specimen with standard shear studs.

Field and Laboratory Evaluation of Precast Concrete Bridges, TR-440, 2001

Recent data compiled by the National Bridge Inventory revealed 29% of Iowa's approximate 24,600 bridges were either structurally deficient or functionally obsolete. This large number of deficient bridges and the high cost of needed repairs create unique problems for Iowa and many other states. The research objective of this project was to determine the load capacity of a particular type of deteriorating bridge – the precast concrete deck bridge – which is commonly found on Iowa's secondary roads. The number of these precast concrete structures requiring load postings and/or replacement can be significantly reduced if the deteriorated structures are found to have adequate load capacity or can be reliably evaluated. Approximately 600 precast concrete deck bridges (PCDBs) exist in Iowa. A typical PCDB span is 19 to 36 ft long and consists of eight to ten simply supported precast panels. Bolts and either a pipe shear key or a grouted shear key are used to join adjacent panels. The panels resemble a steel channel in cross-section; the web is orientated horizontally and forms the roadway deck and the legs act as shallow beams. The primary longitudinal reinforcing steel bundled in each of the legs frequently corrodes and causes longitudinal cracks in the concrete and spalling. The research team performed service load tests on four deteriorated PCDBs; two with shear keys in place and two without. Conventional strain gages were used to measure strains in both the steel and concrete, and transducers were used to measure vertical deflections. Based on the field results, it was determined that these bridges have sufficient lateral load distribution and adequate strength when shear keys are properly installed between adjacent panels. The measured lateral load distribution factors are larger than AASHTO values when shear keys were not installed. Since some of the reinforcement had hooks, deterioration of the reinforcement has a minimal affect on the service level performance of the bridges when there is minimal loss of cross-sectional area. Laboratory tests were performed on the PCDB panels obtained from three bridge replacement projects. Twelve deteriorated panels were loaded to failure in a four point bending arrangement. Although the panels had significant deflections prior to failure, the experimental capacity of eleven panels exceeded the theoretical capacity. Experimental capacity of the twelfth panel, an extremely distressed panel, was only slightly below the theoretical capacity. Service tests and an ultimate strength test were performed on a laboratory bridge model consisting of four joined panels to determine the effect of various shear connection configurations. These data were used to validate a PCDB finite element model that can provide more accurate live load distribution factors for use in rating calculations. Finally, a strengthening system was developed and tested for use in situations where one or more panels of an existing PCDB need strengthening.

Evaluation of the Need for Longitudinal Median Joints in Bridge Decks on Dual Structures, TR-661, 2015

The primary objective of this project was to determine the effect of bridge width on deck cracking in bridges. Other parameters, such as bridge skew, girder spacing and type, abutment type, pier type, and number of bridge spans, were also studied. To achieve the above objectives, one bridge was selected for live-load and long-term testing. The data obtained from both field tests were used to calibrate a three-dimensional (3D) finite element model (FEM). Three different types of loading—live loading, thermal loading, and shrinkage loading—were applied. The predicted crack pattern from the FEM was compared to the crack pattern from bridge inspection results. A parametric study was conducted using the calibrated FEM. The general conclusions/recommendations are as follows: -- Longitudinal and diagonal cracking in the deck near the abutment on an integral abutment bridge is due to the temperature differences between the abutment and the deck. Although not likely to induce cracking, shrinkage of the deck concrete may further exacerbate cracks developed from thermal effects. -- Based upon a limited review of bridges in the Iowa DOT inventory, it appears that, regardless of bridge width, longitudinal and diagonal cracks are prevalent in integral abutment bridges but not in bridges with stub abutments. -- The parametric study results show that bridge width and skew have minimal effect on the strain in the deck bridge resulting from restrained thermal expansion. -- Pier type, girder type, girder spacing, and number of spans also appear to have no influence on the level of restrained thermal expansion strain in the deck near the abutment.

Use of anisotropic modelling in electrical impedance tomography: description of method and preliminary assessment of utility in imaging brain function in the adult human head.

Electrical Impedance Tomography (EIT) is an imaging method which enables a volume conductivity map of a subject to be produced from multiple impedance measurements. It has the potential to become a portable non-invasive imaging technique of particular use in imaging brain function. Accurate numerical forward models may be used to improve image reconstruction but, until now, have employed an assumption of isotropic tissue conductivity. This may be expected to introduce inaccuracy, as body tissues, especially those such as white matter and the skull in head imaging, are highly anisotropic. The purpose of this study was, for the first time, to develop a method for incorporating anisotropy in a forward numerical model for EIT of the head and assess the resulting improvement in image quality in the case of linear reconstruction of one example of the human head. A realistic Finite Element Model (FEM) of an adult human head with segments for the scalp, skull, CSF, and brain was produced from a structural MRI. Anisotropy of the brain was estimated from a diffusion tensor-MRI of the same subject and anisotropy of the skull was approximated from the structural information. A method for incorporation of anisotropy in the forward model and its use in image reconstruction was produced. The improvement in reconstructed image quality was assessed in computer simulation by producing forward data, and then linear reconstruction using a sensitivity matrix approach. The mean boundary data difference between anisotropic and isotropic forward models for a reference conductivity was 50%. Use of the correct anisotropic FEM in image reconstruction, as opposed to an isotropic one, corrected an error of 24 mm in imaging a 10% conductivity decrease located in the hippocampus, improved localisation for conductivity changes deep in the brain and due to epilepsy by 4-17 mm, and, overall, led to a substantial improvement on image quality. This suggests that incorporation of anisotropy in numerical models used for image reconstruction is likely to improve EIT image quality.

Siltanosturin teräsrakenteen suunnittelun avuksi luodun ohjelman jatkokehitys

Diplomityössä tehdään jatkokehitystä KCI Konecranes yrityksen siltanosturin laskentaohjelmaan. Ohjelman tärkeimmät jatkokehityskohteet kartoitettiin käyttäjäkyselyn avulla ja niistä valittiin toivotuimmat, sekä diplomityön lujuusopilliseen aihepiiriin parhaiten soveltuvat. Työhön valitut kaksi aihetta ovat koteloprofiilin kaksiosaisen uuman lujuuslaskennan selvittäminen ja siltanosturin kahdeksanpyöräisenpäätykannattajan elementtimallin suunnittelu. Diplomityössä selvitetään jatkokehityskohteisiin liittyvä teoria, mutta varsinainen ohjelmointi jätetään työn ulkopuolelle. Kaksiosaisella uumalla varustetussa koteloprofiilissa nostovaunun kulkukiskon alla olevan uuman yläosa tehdään paksummaksi, jotta uuma kestäisi nostovaunun pyöräkuormasta aiheutuvan paikallisen jännityksen, eliniin sanotun rusennusjännityksen. Rusennusjännityksen määrittäminen uumalevyissä on kaksiosaisen uuman lujuuslaskennan tärkein tehtävä. Rusennuksen aiheuttamankalvojännityksen ja jännityskeskittymien määrittämiseen erilaisissa konstruktioissa etsittiin sopivimmat menetelmät kirjallisuudesta ja standardeista. Kalvojännitys voidaan määrittää luotettavasti käyttäen joko 45 asteen sääntöä tai standardin mukaista menetelmää ja jännityskonsentraatioiden suuruus saadaan kertomallakalvojännitys jännityskonsentraatiokertoimilla. Menetelmien toimivuus verifioitiin tekemällä kymmeniä uuman elementtimalleja erilaisin dimensioin ja reunaehdoin ja vertaamalla elementtimallien tuloksia käsin laskettuihin. Käsin lasketut jännitykset saatiin vastaamaan tarkasti elementtimallien tuloksia. Kaksiosaisen uuman lommahdus- ja väsymislaskentaa tutkittiin alustavasti. Kahdeksanpyöräisiä päätykannattajia käytetään suurissa siltanostureissa pienentämään pyöräkuormia ja radan rusennusjännityksiä. Kahdeksanpyöräiselle siltanosturin päätykannattajalle suunniteltiin elementtimallit molempiin rakenteesta käytettyihin konstruktioihin: nivelöityyn ja jäykkäkehäiseen malliin. Elementtimallien rakentamisessa hyödynnettiin jo olemassa olevia malleja, jolloin niiden lisääminen ohjelmakoodiin nopeutuu ja ne ovat varmasti yhteensopivia muiden laskentamoduuleiden kanssa. Elementtimallien värähtelyanalyysin reunaehtoja tarkasteltiin. Värähtelyanalyysin reunaehtoihin ei tutkimuksen perusteella tarvitse tehdä muutoksia, mutta staattisen analyysin reunaehdot kaipaavat vielä lisätutkimusta.

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.

Aortic blood pressure measured via EIT: investigation of different measurement settings.

Electrical impedance tomography (EIT) allows the measurement of intra-thoracic impedance changes related to cardiovascular activity. As a safe and low-cost imaging modality, EIT is an appealing candidate for non-invasive and continuous haemodynamic monitoring. EIT has recently been shown to allow the assessment of aortic blood pressure via the estimation of the aortic pulse arrival time (PAT). However, finding the aortic signal within EIT image sequences is a challenging task: the signal has a small amplitude and is difficult to locate due to the small size of the aorta and the inherent low spatial resolution of EIT. In order to most reliably detect the aortic signal, our objective was to understand the effect of EIT measurement settings (electrode belt placement, reconstruction algorithm). This paper investigates the influence of three transversal belt placements and two commonly-used difference reconstruction algorithms (Gauss-Newton and GREIT) on the measurement of aortic signals in view of aortic blood pressure estimation via EIT. A magnetic resonance imaging based three-dimensional finite element model of the haemodynamic bio-impedance properties of the human thorax was created. Two simulation experiments were performed with the aim to (1) evaluate the timing error in aortic PAT estimation and (2) quantify the strength of the aortic signal in each pixel of the EIT image sequences. Both experiments reveal better performance for images reconstructed with Gauss-Newton (with a noise figure of 0.5 or above) and a belt placement at the height of the heart or higher. According to the noise-free scenarios simulated, the uncertainty in the analysis of the aortic EIT signal is expected to induce blood pressure errors of at least ± 1.4 mmHg.

Virtuaaliprototyypin käyttö roottorin dynamiikan analysoinnissa

Työn tavoitteena oli selvittää kaupallisen dynamiikansimulointiohjelmiston soveltuvuus roottoridynamiikan analysointiin. Työssä keskityttiin erityisesti roottorin dynamiikkaan vaikuttavien epäideaalisuuksien mallintamiseen. Simulointitulosten tarkkuutta selvitettiin mittauksilla. Lisäksi vertailtiin yleiskäyttöisen dynamiikan simulointiohjelmiston ja roottoridynamiikan erikoisohjelmiston teoriaa. Tutkittava roottori oli paperikoneen putkitela. Telan joustavuus kuvattiin elementtimenetelmällä ratkaistujen moodien avulla. Elementtimallissa huomioitiin telan vaipan seinämänpaksuusvaihtelu, joka vaikuttaa telan massa- ja jäykkyysjakaumaan. 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 analysointiin ja roottorin epäideaalisuuksien mallintamiseen. Simulointimallilla saatiin esille murtolukukriittiset pyörimisnopeudet ja telan kriittinen pyörimisnopeus vastasi hyvin mittaustuloksia.

Henkilöauton korin mallintaminen virtuaalitestejä varten

Nykyään auton korien tuotekehitysprosessissa jää hyvin vähän aikaa korin rakenteiden optimointiin. Rakenteiden optimointi tehdään virtuaalitestauksen avulla siten, että yhdellä elementtimallilla tehdään erilaisia testejä ja optimoidaan useita omaisuuksia samanaikaisesti. Tässä tutkimuksessa tavoitteena on kehittää auton korin FEM-malli, jolla voidaan suorittaa mahdollisimman monta erilaista virtuaalitestiä. Rakennetun mallin toimivuutta testattiin tekemällä sille ympäriajoturvallisuustesti ja saatuja tuloksia verrattiin oikean auton testiin.

Kaksoislaatuisesta rakenneputkesta hitsatun X-liitoksen kestävyys ja muodonmuutoskyky

Rautaruukki Oyj:n uusi Ruukki double grade S420MH/S355J2H -rakenneputki yhdistää putkilajien S420MH ja S355J2H ominaisuudet. Eurocode 3:n suunnitteluohjeiden mukaan Ruukki double graden käyttöön on sovellettava lujuusluokan S420 mukaista hitsien mitoitusta ja rakenneputkiliitosten staattista mitoituskestävyyttä alentavaa varmuuskerrointa, kun halutaan hyödyntää lujuusluokan S420 ominaisuudet Ruukki double gradea käytettäessä. Tässä työssä tutkittiin S420-lujuusluokkaa vastaavien suunnittelumääräysten soveltamistarvetta Ruukki double grade -rakenneputken käyttöön. Työn tavoitteena oli arvioida Ruukki double grade -rakenneputkesta hitsaamalla valmistetun X-liitoksen kestävyyden ja muodonmuutoskyvyn riittävyys, kun hitsit mitoitettiin S355-lujuusluokkaa vastaavien Eurocode 3:n suunnitteluohjeiden mukaisesti. X-liitosten kestävyyttä ja muodonmuutoskykyä tutkittiin kokeellisesti -40 C lämpötilassa suoritettujen laboratoriokokeiden avulla sekä analyyttisen laskennan ja elementtimenetelmän keinoin. Kestävyyden riittävyyttä arvioitiin vertaamalla laboratoriokokeista saatuja liitosten kestävyyksiä Eurocode 3:n ja myötöviivateorian mukaisiin kestävyyksiin. Elementtimenetelmän ja laboratoriokokeen pohjalta piirrettyjä liitoksen voima-siirtymäkuvaajia vertailtiin keskenään, kun elementtimalli analysoitiin eri materiaalimalleilla. Lisäksi verrattiin elementtimallin voima-venymäkuvaajia liitoksen venymäliuskan arvoihin. Kaikki koesarjan S355-lujuusluokan mukaisilla hitsin a-mitoilla valmistetut liitokset täyttivät kestävyydelle ja muodonmuutoskyvylle asetetut vaatimukset. Täten tämän koesarjan perusteella ei ollut tarpeellista soveltaa S420-lujuusluokan mukaista hitsien mitoitusta ja staattista mitoituskestävyyttä alentavaa varmuuslukua Ruukki double grade -rakenneputkeen. Elementtimenetelmästä ja laboratoriokokeesta saadut liitoksen voima-siirtymäkuvaajat vastasivat hyvin toisiaan. Sen sijaan elementtimallin voima-venymä-kuvaajat eivät vastanneet venymäliuskojen mittausdataa kovin hyvin.

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.

Flexible multibody dynamics and intelligent control of a hydraulically driven hybrid redundant robot machine

The assembly and maintenance of the International Thermonuclear Experimental Reactor (ITER) vacuum vessel (VV) is highly challenging since the tasks performed by the robot involve welding, material handling, and machine cutting from inside the VV. The VV is made of stainless steel, which has poor machinability and tends to work harden very rapidly, and all the machining operations need to be carried out from inside of the ITER VV. A general industrial robot cannot be used due to its poor stiffness in the heavy duty machining process, and this will cause many problems, such as poor surface quality, tool damage, low accuracy. Therefore, one of the most suitable options should be a light weight mobile robot which is able to move around inside of the VV and perform different machining tasks by replacing different cutting tools. Reducing the mass of the robot manipulators offers many advantages: reduced material costs, reduced power consumption, the possibility of using smaller actuators, and a higher payload-to-robot weight ratio. Offsetting these advantages, the lighter weight robot is more flexible, which makes it more difficult to control. To achieve good machining surface quality, the tracking of the end effector must be accurate, and an accurate model for a more flexible robot must be constructed. This thesis studies the dynamics and control of a 10 degree-of-freedom (DOF) redundant hybrid robot (4-DOF serial mechanism and 6-DOF 6-UPS hexapod parallel mechanisms) hydraulically driven with flexible rods under the influence of machining forces. Firstly, the flexibility of the bodies is described using the floating frame of reference method (FFRF). A finite element model (FEM) provided the Craig-Bampton (CB) modes needed for the FFRF. A dynamic model of the system of six closed loop mechanisms was assembled using the constrained Lagrange equations and the Lagrange multiplier method. Subsequently, the reaction forces between the parallel and serial parts were used to study the dynamics of the serial robot. A PID control based on position predictions was implemented independently to control the hydraulic cylinders of the robot. Secondly, in machining, to achieve greater end effector trajectory tracking accuracy for surface quality, a robust control of the actuators for the flexible link has to be deduced. This thesis investigates the intelligent control of a hydraulically driven parallel robot part based on the dynamic model and two schemes of intelligent control for a hydraulically driven parallel mechanism based on the dynamic model: (1) a fuzzy-PID self-tuning controller composed of the conventional PID control and with fuzzy logic, and (2) adaptive neuro-fuzzy inference system-PID (ANFIS-PID) self-tuning of the gains of the PID controller, which are implemented independently to control each hydraulic cylinder of the parallel mechanism based on rod length predictions. The serial component of the hybrid robot can be analyzed using the equilibrium of reaction forces at the universal joint connections of the hexa-element. To achieve precise positional control of the end effector for maximum precision machining, the hydraulic cylinder should be controlled to hold the hexa-element. Thirdly, a finite element approach of multibody systems using the Special Euclidean group SE(3) framework is presented for a parallel mechanism with flexible piston rods under the influence of machining forces. The flexibility of the bodies is described using the nonlinear interpolation method with an exponential map. The equations of motion take the form of a differential algebraic equation on a Lie group, which is solved using a Lie group time integration scheme. The method relies on the local description of motions, so that it provides a singularity-free formulation, and no parameterization of the nodal variables needs to be introduced. The flexible slider constraint is formulated using a Lie group and used for modeling a flexible rod sliding inside a cylinder. The dynamic model of the system of six closed loop mechanisms was assembled using Hamilton’s principle and the Lagrange multiplier method. A linearized hydraulic control system based on rod length predictions was implemented independently to control the hydraulic cylinders. Consequently, the results of the simulations demonstrating the behavior of the robot machine are presented for each case study. In conclusion, this thesis studies the dynamic analysis of a special hybrid (serialparallel) robot for the above-mentioned special task involving the ITER and investigates different control algorithms that can significantly improve machining performance. These analyses and results provide valuable insight into the design and control of the parallel robot with flexible rods.

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.

Real-time simulation of multibody systems with applications for working mobile vehicles

This dissertation describes an approach for developing a real-time simulation for working mobile vehicles based on multibody modeling. The use of multibody modeling allows comprehensive description of the constrained motion of the mechanical systems involved and permits real-time solving of the equations of motion. By carefully selecting the multibody formulation method to be used, it is possible to increase the accuracy of the multibody model while at the same time solving equations of motion in real-time. In this study, a multibody procedure based on semi-recursive and augmented Lagrangian methods for real-time dynamic simulation application is studied in detail. In the semirecursive approach, a velocity transformation matrix is introduced to describe the dependent coordinates into relative (joint) coordinates, which reduces the size of the generalized coordinates. The augmented Lagrangian method is based on usage of global coordinates and, in that method, constraints are accounted using an iterative process. A multibody system can be modelled as either rigid or flexible bodies. When using flexible bodies, the system can be described using a floating frame of reference formulation. In this method, the deformation mode needed can be obtained from the finite element model. As the finite element model typically involves large number of degrees of freedom, reduced number of deformation modes can be obtained by employing model order reduction method such as Guyan reduction, Craig-Bampton method and Krylov subspace as shown in this study The constrained motion of the working mobile vehicles is actuated by the force from the hydraulic actuator. In this study, the hydraulic system is modeled using lumped fluid theory, in which the hydraulic circuit is divided into volumes. In this approach, the pressure wave propagation in the hoses and pipes is neglected. The contact modeling is divided into two stages: contact detection and contact response. Contact detection determines when and where the contact occurs, and contact response provides the force acting at the collision point. The friction between tire and ground is modelled using the LuGre friction model, which describes the frictional force between two surfaces. Typically, the equations of motion are solved in the full matrices format, where the sparsity of the matrices is not considered. Increasing the number of bodies and constraint equations leads to the system matrices becoming large and sparse in structure. To increase the computational efficiency, a technique for solution of sparse matrices is proposed in this dissertation and its implementation demonstrated. To assess the computing efficiency, augmented Lagrangian and semi-recursive methods are implemented employing a sparse matrix technique. From the numerical example, the results show that the proposed approach is applicable and produced appropriate results within the real-time period.