A computationally efficient shear deformable beam element for large deformation multibody applications

In this paper, a new two-dimensional shear deformable beam element based on the absolute nodal coordinate formulation is proposed. The nonlinear elastic forces of the beam element are obtained using a continuum mechanics approach without employing a local element coordinate system. In this study, linear polynomials are used to interpolate both the transverse and longitudinal components of the displacement. This is different from other absolute nodal-coordinate-based beam elements where cubic polynomials are used in the longitudinal direction. The accompanying defects of the phenomenon known as shear locking are avoided through the adoption of selective integration within the numerical integration method. The proposed element is verified using several numerical examples, and the results are compared to analytical solutions and the results for an existing shear deformable beam element. It is shown that by using the proposed element, accurate linear and nonlinear static deformations, as well as realistic dynamic behavior, can be achieved with a smaller computational effort than by using existing shear deformable two-dimensional beam elements.

Experimental insights into deformation dynamics and intermittency in rapid granular shear flows

This work concerns the experimental study of rapid granular shear flows in annular Couette geometry. The flow is induced by continuous driving of the horizontal plate at the top of the granular bed in an annulus. The compressive pressure, driving torque, instantaneous bed height and rotational speed of the shearing plate are measured. Moreover, local stress fluctuations are measured in a medium made of steel spheres 2 and 3 mm in diameter. Both monodisperse packing and bidisperse packing are investigated to reveal the influence of size diversity in intermittent features of granular materials. Experiments are conducted in an annulus that can contain up to 15 kg of spherical steel balls. The shearing granular medium takes place via the rotation of the upper plate which compresses the material loaded inside the annulus. Fluctuations of compressive force are locally measured at the bottom of the annulus using a piezoelectric sensor. Rapid shear flow experiments are pursued at different compressive forces and shear rates and the sensitivity of fluctuations are then investigated by different means through monodisperse and bidisperse packings. Another important feature of rapid granular shear flows is the formation of ordered structures upon shearing. It requires a certain range for the amount of granular material (uniform size distribution) loaded in the system in order to obtain stable flows. This is studied more deeply in this thesis. The results of the current work bring some new insights into deformation dynamics and intermittency in rapid granular shear flows. The experimental apparatus is modified in comparison to earlier investigations. The measurements produce data for various quantities continuously sampled from the start of shearing to the end. Static failure and dynamic shearing ofa granular medium is investigated. The results of this work revealed some important features of failure dynamics and structure formation in the system. Furthermore, some computer simulations are performed in a 2D annulus to examine the nature of kinetic energy dissipation. It is found that turbulent flow models can statistically represent rapid granular flows with high accuracy. In addition to academic outcomes and scientific publications our results have a number of technological applications associated with grinding, mining and massive grain storages.

The Sottunga-Jurmo shear zone : structure and deformation history of a crustal-scale ductile shear zone in SW Finland

The aim of this study is to gain a better understanding of the structure and the deformation history of a NW-SE trending regional, crustal-scale shear structure in the Åland archipelago, SW Finland, called the Sottunga-Jurmo shear zone (SJSZ). Approaches involving e.g. structural geology, geochronology, geochemistry and metamorphic petrology were utilised in order to reconstruct the overall deformation history of the study area. The study therefore describes several features of the shear zone including structures, kinematics and lithologies within the study area, the ages of the different deformation phases (ductile to brittle) within the shear zone, as well as some geothermobarometric results. The results indicate that the SJSZ outlines a major crustal discontinuity between the extensively migmatized rocks NE of the shear zone and the unmigmatised, amphibolite facies rocks SW of the zone. The main SJSZ shows overall dextral lateral kinematics with a SW-side up vertical component and deformation partitioning into pure shear and simple shear dominated deformation styles that was intensified toward later stages of the deformation history. The deformation partitioning resulted in complex folding and refolding against the SW margin of the SJSZ, including conical and sheath folds, and in a formation of several minor strike-slip shear zones both parallel and conjugate to the main SJSZ in order to accommodate the regional transpressive stresses. Different deformation phases within the study area were dated by SIMS (zircon U-Pb), ID-TIMS (titanite U-Pb) and 40Ar/39Ar (pseudotachylyte wholerock) methods. The first deformation phase within the ca. 1.88 Ga rocks of the study area is dated at ca. 1.85 Ga, and the shear zone was reactivated twice within the ductile regime (at ca. 1.83 Ga and 1.79 Ga), during which the strain was successively increasingly partitioned into the main SJSZ and the minor shear zones. The age determinations suggest that the orogenic processes within the study area did not occur in a temporal continuum; instead, the metamorphic zircon rims and titanites show distinct, 10-20 Ma long breaks in deformation between phases of active deformation. The results of this study further imply slow cooling of the rocks through 600-700ºC so that at 1.79 Ga, 2 the temperature was still at least 600ºC. The highest recorded metamorphic pressures are 6.4-7.1 kbar. At the late stages or soon after the last ductile phase (ca. 1.79 Ga), relatively high-T mylonites and ultramylonites were formed, witnessing extreme deformation partitioning and high strain rates. After the rocks reached lower amphibolite facies to amphibolite-greenschist facies transitional conditions (ca. 500-550ºC), they cooled rapidly, probably due to crustal uplift and exhumation. The shear zone was reactivated at least once within the semi-brittle to brittle regime between ca. 1.79 Ga and 1.58 Ga, as evidenced by cataclasites and pseudotachylytes. In summary, the results of this study suggest that the Sottunga-Jurmo shear zone (and the South Finland shear zone) defines a major crustal discontinuity, and played a central role in accommodating the regional stresses during and after the Svecofennian orogeny.

The Sottunga-Jurmo shear zone : structure and deformation history of a crustal-scale ductile shear zone in SW Finland

The aim of this study is to gain a better understanding of the structure and the deformation history of a NW-SE trending regional, crustal-scale shear structure in the Åland archipelago, SW Finland, called the Sottunga-Jurmo shear zone (SJSZ). Approaches involving e.g. structural geology, geochronology, geochemistry and metamorphic petrology were utilised in order to reconstruct the overall deformation history of the study area. The study therefore describes several features of the shear zone including structures, kinematics and lithologies within the study area, the ages of the different deformation phases (ductile to brittle) within the shear zone, as well as some geothermobarometric results. The results indicate that the SJSZ outlines a major crustal discontinuity between the extensively migmatized rocks NE of the shear zone and the unmigmatised, amphibolite facies rocks SW of the zone. The main SJSZ shows overall dextral lateral kinematics with a SW-side up vertical component and deformation partitioning into pure shear and simple shear dominated deformation styles that was intensified toward later stages of the deformation history. The deformation partitioning resulted in complex folding and refolding against the SW margin of the SJSZ, including conical and sheath folds, and in a formation of several minor strike-slip shear zones both parallel and conjugate to the main SJSZ in order to accommodate the regional transpressive stresses. Different deformation phases within the study area were dated by SIMS (zircon U-Pb), ID-TIMS (titanite U-Pb) and 40Ar/39Ar (pseudotachylyte wholerock) methods. The first deformation phase within the ca. 1.88 Ga rocks of the study area is dated at ca. 1.85 Ga, and the shear zone was reactivated twice within the ductile regime (at ca. 1.83 Ga and 1.79 Ga), during which the strain was successively increasingly partitioned into the main SJSZ and the minor shear zones. The age determinations suggest that the orogenic processes within the study area did not occur in a temporal continuum; instead, the metamorphic zircon rims and titanites show distinct, 10-20 Ma long breaks in deformation between phases of active deformation. The results of this study further imply slow cooling of the rocks through 600-700ºC so that at 1.79 Ga, 2 the temperature was still at least 600ºC. The highest recorded metamorphic pressures are 6.4-7.1 kbar. At the late stages or soon after the last ductile phase (ca. 1.79 Ga), relatively high-T mylonites and ultramylonites were formed, witnessing extreme deformation partitioning and high strain rates. After the rocks reached lower amphibolite facies to amphibolite-greenschist facies transitional conditions (ca. 500-550ºC), they cooled rapidly, probably due to crustal uplift and exhumation. The shear zone was reactivated at least once within the semi-brittle to brittle regime between ca. 1.79 Ga and 1.58 Ga, as evidenced by cataclasites and pseudotachylytes. In summary, the results of this study suggest that the Sottunga-Jurmo shear zone (and the South Finland shear zone) defines a major crustal discontinuity, and played a central role in accommodating the regional stresses during and after the Svecofennian orogeny.

Computational studies of deformation in stainless steel rings due to torch heating

A support ring of AISI 304L stainless steel that holds vertical, parallel wires arranged in a circle forming a cylinder is studied. The wires are attached to the ring with heat-induced shrinkage. When the ring is heated with a torch the heat affected zone tries to expand while the adjacent cool structure obstructs the expansion causing upsetting. During cooling, the ring shrinks smaller than its original size clamping the wires. The most important requirement for the ring is that it should be as round as possible and the deformations should occur as overall shrinkage in the ring diameter. A three-dimensional nonlinear transient sequential thermo-structural Abaqus model is used together with a Fortran code that enters the heat flux to each affected element. The local and overall deformations in one ring inflicted by the heating are studied with a small amount of inspection on residual stresses. A variety of different cases are chosen to be studied with the model constructed to provide directional knowledge; torch flux with the means of speed, location of the wires, heating location and structural factors. The decrease of heating speed increases heat flux that rises the temperature increasing shrinkage. In a single progressive heating uneven distribution of shrinkage appears to the start/end region that can be partially fixed with using speeded heating’s to strengthen the heating of that region. Location of the wires affect greatly to the caused shrinkage unlike heating location. The ring structure affects also greatly to the shrinkage; smaller diameter, bigger ring height, thinner thickness and greater number of wires increase shrinkage.

Cardiac Imaging in the Diagnosis of Coronary Artery Disease: The Value of Quantitative Imaging of Myocardial Perfusion and Deformation

Atherosclerosis is a chronic and progressive disease of the vasculature. Increasing coronary atherosclerosis can lead to obstructive coronary artery disease (CAD) or myocardial infarction. Computed tomography angiography (CTA) allows noninvasive assessment of coronary anatomy and quantitation of atherosclerotic burden. Myocardial blood flow (MBF) can be accurately measured in absolute terms (mL/g/min) by positron emission tomography (PET) with [15O] H O as a radiotracer. We studied the coronary microvascular dysfunction as a risk factor for future coronary calcification in healthy young men by measuring the coronary flow reserve (CFR) which is the ratio between resting and hyperemic MBF. Impaired vasodilator function was not linked with accelerated atherosclerosis 11 years later. Currently, there is a global interest in quantitative PET perfusion imaging. We established optimal thresholds of [15O] H O PET perfusion for diagnosis of CAD (hyperemic MBF of 2.3 mL/g/min and CFR of 2.5) in the first multicenter study of this type (Turku, Amsterdam and Uppsala). In myocardial bridging a segment of the coronary artery travels inside the myocardium and can be seen as intramural course (CTA) or systolic compression (invasive coronary angiography). Myocardial bridging is frequently linked with proximal atherosclerotic plaques. We used quantitative [15O] H O PET perfusion to evaluate the hemodynamic effects of myocardial bridging. Myocardial bridging was not associated with decreased absolute MBF or increased atherosclerotic burden. Speckle tracking allows quantitative echocardiographic imaging of myocardial deformation. Speckle tracking during dobutamine stress echocardiography was feasible and comparable to subjective wall motion analysis in the diagnosis of CAD. In addition, it correctly risk stratified patients with multivessel disease and extensive ischemia.

Laser-MAG-hybridihitsauksen käyttöönotto laivan rungon valmistuksessa

Perinteisten kaarihitsausmenetelmien suhteellisen suuri lämmöntuonti aiheuttaa huomattavia muodonmuutoksia laivan rungon valmistusprosessin alkuvaiheessa. Muodonmuutosten seurauksena rakenteiden mitta- ja muototarkkuus heikkenee, mikä lisää oikaisu- ja sovitustyötä myöhemmissä työvaiheissa. Hitsausmuodonmuutoksia voidaan vähentää siirtymällä käyttämään laser-MAG-hybridihitsausta, jossa lämmöntuonti on merkittävästi pienempi kuin kaarihitsauksessa. Näin kyetään oleellisesti leikkaamaan oikaisu- ja sovitustyöstä syntyviä kustannuksia. Tämän diplomityön tavoitteena oli kehittää tuotantovalmiiksi kuitulaser- ja MAG-hitsauksen yhdistelmäprosessi Aker Yards Oy:n Turun telakalla loppuvuoden 2006 aikana. Hitsauslaitteiston asennus oli valmistunut kesäkuussa 2006, minkä jälkeen aloitettiin luokituslaitoksen hyväksymän koeohjelman hitsaukset. Käyttöönotto suunnitelmaan sisältyvä koehitsausohjelma oli laadittu Det Norske Veritaksen julkaisemaa ohjetta (Guidelines no. 19) mukaillen. Ensimmäiseksi määritettiin hitsauskokeiden avulla prosessille laadun ja tehokkuuden suhteen optimaalinen railogeometria. Seuraavaksi optimoitiin prosessin hitsausparametrit 6 mm:n aineenpaksuudelle hyödyntäen Taguchi-koesuunnittelumenetelmää. Tämän jälkeen optimiparametreilla hitsattiin koekappale väsytyskokeisiin, jotka suoritettiin Teknillisen korkeakoulun laivalaboratoriossa. Väsytyskoetulokset täyttivät luokituslaitoksen vaatimukset. Myös hitsauksen menetelmäkoe suoritettiin hyväksytetysti. Viimeinen koeohjelman mukainen hitsauskoesarja tehtiin prosessiparametrien sallittujen vaihtelurajojen määrittämiseksi. Diplomityön tavoite täyttyi joulukuussa 2006, jolloin 'laivan kansipaneeli hitsattiin ensimmäistä kertaa uudella hitsausprosessilla. Hitsauksen laatu korreloi hyvin menetelmäkokeen tulosten kanssa ¿ hitsit olivat tasalaatuisia ja ne täyttivät B-luokan vaatimukset.

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.

Development of a Testing Method for Convertibility of Paperboard

Työssä tutkittiin muovattujen kartonkivuokien sekä muovattujen kartonkinäytteiden rinnastettavuutta. Puristusvaiheen prosessiolosuhteiden miellettiin vaikuttavan eniten multidimensionaliseen muodonmuutokseen. Multidimensionaalista muodonmuutosta simuloitiin uudella muovaamiseen soveltuvalla muovauslaitteella. Kirjallisuusosassa keskeisiä teemoja ovat kartongin muovaus sekä kuitupohjaisen materiaalin reologinen käyttäytyminen. Kirjallisuusosassa esitellään lisäksi yksi tekninen sovellus, jonka avulla kyetään ennustamaan kuitumateriaalin muovautuvuutta sekä mittaamaan tapahtunutta muodonmuutosta. Prosessiparametrien teoreettista vaikutustakuituihin tarkastellaan myös kirjallisuusosassa. Kokeellisessa osassa toteutettiin kartonkivuokien valmistus puristamalla. Vastaavilla prosessiparametreilla muovattiin myös pienemmät testinäytteet. Perinteiset yksidimensionaliset deformaatiomittaukset toteutettiin lujuusominaisuuksien laboratoriomäärityksinä. Myös kitka, joka toimii tärkeänä muuttujana prässäysprosessissa, mitattiin laboratorio-olosuhteissa. Tämän työn tulokset osoittavat uuden kehitetyn muovausmenetelmän toimivuuden. Asema-voima kuvaajat ovat selkeitä sekä helposti luettavia. Tuloksissa havaittiin materiaalin muovauspotentiaalin sekä asema-voima kuvaajan välillä vallitseva yhteys. Erittäin merkittävä huomio oli myös, että muovipäällystetyllä kartongilla oli yhteys päällystämättömän kartongin asema-voima kuvaajaan. Tämä tulos osoittaa, että muovipäällystetyn kartongin muovautuvuutta voi olla mahdollista ennustaa pohjakartongin muovautuvuustulosten perusteella. Perinteiset yksidimensionaliset laboratoriomittaukset eivät kykene antamaan riittävää informaatiota muovautuvuuden ennustamiseen. Tästä näkökulmasta on tärkeää että kartongin multidimensionalista muotoutuvuutta voidaankin tutkia kehitetyllä muovausmenetelmällä.