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.

On the Application of Beam on Elastic Foundation Theory to the Analysis of Stiffened Plate Strips

The main objective of this thesis is to show that plate strips subjected to transverse line loads can be analysed by using the beam on elastic foundation (BEF) approach. It is shown that the elastic behaviour of both the centre line section of a semi infinite plate supported along two edges, and the free edge of a cantilever plate strip can be accurately predicted by calculations based on the two parameter BEF theory. The transverse bending stiffness of the plate strip forms the foundation. The foundation modulus is shown, mathematically and physically, to be the zero order term of the fourth order differential equation governing the behaviour of BEF, whereas the torsion rigidity of the plate acts like pre tension in the second order term. Direct equivalence is obtained for harmonic line loading by comparing the differential equations of Levy's method (a simply supported plate) with the BEF method. By equating the second and zero order terms of the semi infinite BEF model for each harmonic component, two parameters are obtained for a simply supported plate of width B: the characteristic length, 1/ λ, and the normalized sum, n, being the effect of axial loading and stiffening resulting from the torsion stiffness, nlin. This procedure gives the following result for the first mode when a uniaxial stress field was assumed (ν = 0): 1/λ = √2B/π and nlin = 1. For constant line loading, which is the superimposition of harmonic components, slightly differing foundation parameters are obtained when the maximum deflection and bending moment values of the theoretical plate, with v = 0, and BEF analysis solutions are equated: 1 /λ= 1.47B/π and nlin. = 0.59 for a simply supported plate; and 1/λ = 0.99B/π and nlin = 0.25 for a fixed plate. The BEF parameters of the plate strip with a free edge are determined based solely on finite element analysis (FEA) results: 1/λ = 1.29B/π and nlin. = 0.65, where B is the double width of the cantilever plate strip. The stress biaxial, v > 0, is shown not to affect the values of the BEF parameters significantly the result of the geometric nonlinearity caused by in plane, axial and biaxial loading is studied theoretically by comparing the differential equations of Levy's method with the BEF approach. The BEF model is generalised to take into account the elastic rotation stiffness of the longitudinal edges. Finally, formulae are presented that take into account the effect of Poisson's ratio, and geometric non linearity, on bending behaviour resulting from axial and transverse inplane loading. It is also shown that the BEF parameters of the semi infinite model are valid for linear elastic analysis of a plate strip of finite length. The BEF model was verified by applying it to the analysis of bending stresses caused by misalignments in a laboratory test panel. In summary, it can be concluded that the advantages of the BEF theory are that it is a simple tool, and that it is accurate enough for specific stress analysis of semi infinite and finite plate bending problems.

Harustetun 110 kV kannatuspylvään konseptin kehittäminen

Diplomityössä kehitettiin harustetun 110 kV kannatuspylvään konsepti tuotteeksi. Pylväs on säänkestävästä teräksestä valmistettu putkipalkkirakenteinen I-pylväs. Tavoitteena oli suunnitella rakenteesta kokonaistaloudellisesti edullinen. Rakenteen suunnittelussa otettiin huomioon valmistus-, kuljetus- ja varastointi- sekä rakentamisnäkökohtia. Työssä perehdyttiin pylväsrakenteiden yksityiskohtiin, putkipalkkien liitosmenetelmiin ja pylvään jalan nivelöintiratkaisuihin. Säänkestävä rakennemateriaali otettiin huomioon rakennesuunnittelussa. Rakenteen lujuusteknisen suunnittelun apuna käytettiin epälineaarista elementtimenetelmää. Pylväsrakenteen käyttäytyminen mallinnettiin geometrisesti epälineaariseksi, ja liitosdetaljien analysointia varten kehitettiin epälineaarisia materiaalimalleja. Rakenteen värähtelykäyttäytyminen analysoitiin myös elementtimenetelmällä. Lopputuloksena saatiin aikaan pylväs, joka täyttää sille asetetut vaatimukset. Pylväs on helposti valmistettava, kuljetettava ja pystytettävä.

Product costing model for laser welded hollow core steel panels

The aim of the study was to create an easily upgradable product costing model for laser welded hollow core steel panels to help in pricing decisions. The theory section includes a literature review to identify traditional and modern cost accounting methodologies, which are used by manufacturing companies. The theory section also presents the basics of steel panel structures and their manufacturing methods and manufacturing costs based on previous research. Activity-Based costing turned out to be the most appropriate methodology for the costing model because of wide product variations. Activity analysis and the determination of cost drivers based on observations and interviews were the key steps in the creation of the model. The created model was used to test how panel parameters affect the costs caused by the main manufacturing stages and materials. By comparing cost structures, it was possible to find the panel types that are the most economic and uneconomic to manufacture. A sensitivity analysis proved that the model gives sufficiently reliable cost information to support pricing decisions. More reliable cost information could be achieved by determining the cost drivers more accurately. Alternative methods for manufacturing the cores were compared with the model. The comparison proved that roll forming can be more advantageous and flexible than press brake bending. However, more extensive research showed that roll forming is possible only when the cores are designed to be manufactured by roll forming. Due to that fact, when new panels are designed consideration should be given to the possibility of using roll forming.

The energy efficiency of chilled beam systems in an office building at different climates

Besides the sustaining of healthy and comfortable indoor climate, the air conditioning system should also achieve for energy efficiency. The target indoor climate can be ob-tained with different systems; this study focuses on comparing the energy efficiency of different air conditioning room unit systems in different climates. The calculations are made with dynamic energy simulation software IDA ICE by comparing the indoor cli-mate and energy consumption of an office building with different systems in different climates. The aim of the study is to compare the energy efficiency of chilled beam systems to other common systems: variable air volume, fan coil and radiant ceiling systems. Besides the annual energy consumption also the sustainability of target indoor climate is compared between the simulations. Another aim is to provide conclusions to be used in the product development of the chilled beam systems’ energy efficiency. The adaptable chilled beam system and the radiant ceiling system prove to be energy efficient independent of the climate. The challenge of reliable comparison is that other systems are not able to reach the target indoor climate as well as the others. The complex calculation environment of the simulation software, made assumptions and excluding of the financial aspects complicate comparing the big picture. The results show that the development of the chilled beam systems should concentrate on energy efficient night heating, flexible demand based ventilation and capacity control and possibilities on integrating the best practices with other systems. 

Characterisation of Laser Beam and Paper Material Interaction

It is known already from 1970´s that laser beam is suitable for processing paper materials. In this thesis, term paper materials mean all wood-fibre based materials, like dried pulp, copy paper, newspaper, cardboard, corrugated board, tissue paper etc. Accordingly, laser processing in this thesis means all laser treatments resulting material removal, like cutting, partial cutting, marking, creasing, perforation etc. that can be used to process paper materials. Laser technology provides many advantages for processing of paper materials: non-contact method, freedom of processing geometry, reliable technology for non-stop production etc. Especially packaging industry is very promising area for laser processing applications. However, there are only few industrial laser processing applications worldwide even in beginning of 2010´s. One reason for small-scale use of lasers in paper material manufacturing is that there is a shortage of published research and scientific articles. Another problem, restraining the use of laser for processing of paper materials, is colouration of paper material i.e. the yellowish and/or greyish colour of cut edge appearing during cutting or after cutting. These are the main reasons for selecting the topic of this thesis to concern characterization of interaction of laser beam and paper materials. This study was carried out in Laboratory of Laser Processing at Lappeenranta University of Technology (Finland). Laser equipment used in this study was TRUMPF TLF 2700 carbon dioxide laser that produces a beam with wavelength of 10.6 μm with power range of 190-2500 W (laser power on work piece). Study of laser beam and paper material interaction was carried out by treating dried kraft pulp (grammage of 67 g m-2) with different laser power levels, focal plane postion settings and interaction times. Interaction between laser beam and dried kraft pulp was detected with different monitoring devices, i.e. spectrometer, pyrometer and active illumination imaging system. This way it was possible to create an input and output parameter diagram and to study the effects of input and output parameters in this thesis. When interaction phenomena are understood also process development can be carried out and even new innovations developed. Fulfilling the lack of information on interaction phenomena can assist in the way of lasers for wider use of technology in paper making and converting industry. It was concluded in this thesis that interaction of laser beam and paper material has two mechanisms that are dependent on focal plane position range. Assumed interaction mechanism B appears in range of average focal plane position of 3.4 mm and 2.4 mm and assumed interaction mechanism A in range of average focal plane position of 0.4 mm and -0.6 mm both in used experimental set up. Focal plane position 1.4 mm represents midzone of these two mechanisms. Holes during laser beam and paper material interaction are formed gradually: first small hole is formed to interaction area in the centre of laser beam cross-section and after that, as function of interaction time, hole expands, until interaction between laser beam and dried kraft pulp is ended. By the image analysis it can be seen that in beginning of laser beam and dried kraft pulp material interaction small holes off very good quality are formed. It is obvious that black colour and heat affected zone appear as function of interaction time. This reveals that there still are different interaction phases within interaction mechanisms A and B. These interaction phases appear as function of time and also as function of peak intensity of laser beam. Limit peak intensity is the value that divides interaction mechanism A and B from one-phase interaction into dual-phase interaction. So all peak intensity values under limit peak intensity belong to MAOM (interaction mechanism A one-phase mode) or to MBOM (interaction mechanism B onephase mode) and values over that belong to MADM (interaction mechanism A dual-phase mode) or to MBDM (interaction mechanism B dual-phase mode). Decomposition process of cellulose is evolution of hydrocarbons when temperature is between 380- 500°C. This means that long cellulose molecule is split into smaller volatile hydrocarbons in this temperature range. As temperature increases, decomposition process of cellulose molecule changes. In range of 700-900°C, cellulose molecule is mainly decomposed into H2 gas; this is why this range is called evolution of hydrogen. Interaction in this range starts (as in range of MAOM and MBOM), when a small good quality hole is formed. This is due to “direct evaporation” of pulp via decomposition process of evolution of hydrogen. And this can be seen can be seen in spectrometer as high intensity peak of yellow light (in range of 588-589 nm) which refers to temperature of ~1750ºC. Pyrometer does not detect this high intensity peak since it is not able to detect physical phase change from solid kraft pulp to gaseous compounds. As interaction time between laser beam and dried kraft pulp continues, hypothesis is that three auto ignition processes occurs. Auto ignition of substance is the lowest temperature in which it will spontaneously ignite in a normal atmosphere without an external source of ignition, such as a flame or spark. Three auto ignition processes appears in range of MADM and MBDM, namely: 1. temperature of auto ignition of hydrogen atom (H2) is 500ºC, 2. temperature of auto ignition of carbon monoxide molecule (CO) is 609ºC and 3. temperature of auto ignition of carbon atom (C) is 700ºC. These three auto ignition processes leads to formation of plasma plume which has strong emission of radiation in range of visible light. Formation of this plasma plume can be seen as increase of intensity in wavelength range of ~475-652 nm. Pyrometer shows maximum temperature just after this ignition. This plasma plume is assumed to scatter laser beam so that it interacts with larger area of dried kraft pulp than what is actual area of beam cross-section. This assumed scattering reduces also peak intensity. So result shows that assumably scattered light with low peak intensity is interacting with large area of hole edges and due to low peak intensity this interaction happens in low temperature. So interaction between laser beam and dried kraft pulp turns from evolution of hydrogen to evolution of hydrocarbons. This leads to black colour of hole edges.

Capacity of an I-beam made of S960 QC based on Eurocode 3

The capacity of beams is a very important factor in the study of durability of structures and structural members. The capacity of a high-strength steel I-beam made of S960 QC was investigated in this study. The investigation included assessment of the service limits and ultimate limits of the steel beam. The thesis was done according to European standards for steel structures, Eurocode 3. An analytical method was used to determine the throat thickness, deformation, elastic and plastic moment capacities as well as the fatigue life of the beam. The results of the analytical method were compared with those obtained by Finite Element Analysis (FEA). Elastic moment capacity obtained by the analytical method was 172 kNm. FEA and the analytical method predicted the maximum lateral-torsional buckling (LTB) capacity in the range of 90-93 kNm and the probability of failure as a result of LTB is estimated to be 50%. The lateral buckling capacity meant that the I-beam can carry a safe load of 300 kN instead of the initial load of 600 kN. The beam is liable to fail shortly after exceeding the elastic moment capacity. Based on results in of the different approaches, it was noted that FEA predicted higher deformation values on the load-deformation curve than the analytical results. However, both FEA and the analytical methods predicted identical results for nominal stress range and moment capacities. Fatigue life was estimated to be in the range of 53000-64000 cycles for bending stress range using crack propagation equation and strength-life approach. As Eurocode 3 is limited to steel grades up to S690, results for S960 must be verified with experimental data and appropriate design rules.

Chaos and Unpredictability in the Vibration of an Elasto-Plastic Beam

This contribution discusses the nonlinear dynamics of a pin-ended elasto-plastic beam with both kinematic and isotropic hardening. An iterative numerical procedure based on the operator split technique is developed in order to deal with the nonlinearities in the equations of motion. Free and forced responses for harmonic sinusoidal and square wave excitations are investigated. Numerical simulations present many interesting behaviors such as jump phenomena, sensitivity to initial conditions, chaos and transient chaos. These results indicate that there are practical problems in predicting the response of the beam even when periodic steady state response is expected.

Experimental Vibration Characteristics of a Beam with Nonlinear Support Using Receptance Coupling Analysis

This paper applies the Multi-Harmonic Nonlinear Receptance Coupling Approach (MUHANORCA) (Ferreira 1998) to evaluate the frequency response characteristics of a beam which is clamped at one end and supported at the other end by a nonlinear cubic stiffness joint. In order to apply the substructure coupling technique, the problem was characterised by coupling a clamped linear beam with a nonlinear cubic stiffness joint. The experimental results were obtained by a sinusoidal excitation with a special force control algorithm where the level of the fundamental force is kept constant and the level of the harmonics is kept zero for all the frequencies measured.

Ultimate capacity of stainless steel rectangular hollow section X- joints at low temperatures

Growing demand for stainless steel construction materials has increased the popularity of substitutive materials for austenitic stainless steels. The lean duplex grades have taken their place in building of structures exposed to corrosive environments. Since the duplex grades are relatively new materials, the current codes and norms do not fully cover the newest duplex grades. The joints tested in this thesis were designed and studied according to Eurocode 3, even though all the materials are not yet accepted to the standards. The main objective in this thesis was to determine the differences of the used materials in behaviour under loading at low temperatures. Tests in which the deformation and strength properties of the joints were determined were done at the temperature of -46°C, which is the requirement of temperature for structures designed according to Norsok standards. Results show that replacing the austenitic grade with the lean duplex grade is acceptable.

SolidWorks Simulations ja Femap FE-analyysiohjelmien vertailu RHS-putkipalkkien mallinnuksen avulla

Tässä kandidaatintyössä vertailtiin SolidWorks Simulations ja Femap FE-analyysiohjelmien tuottamia tuloksia. Vertailu suoritettiin yksinkertaisten RHS-putkipalkkien mallinnuksen avulla. Mallinnuksessa käytettiin viittä erilaista kuormitustapausta, jotka sisälsivät mm. taivutus-, vääntö- ja jännityskeskittymätarkastelut. Ohjelmien tuottamia tuloksia verrattiin mahdollisuuksien mukaan analyyttisin menetelmin tehtyihin vertailulaskelmiin. Työn tarkoituksena oli selvittää tuottivatko eri ohjelmat erilaisia tuloksia samoista kuormitustapauksista ja kuinka paljon tulokset erosivat.

Pituussuunnassa hitsatun rakenneputken väsymislujuus HF –hitsausta käytettäessä

Tässä diplomityössä määritettiin pituussuunnassa hitsatun rakenneputken väsymislujuus HF –hitsausta käytettäessä. Työn tavoitteena oli vertailla eri jälkikäsittelyjen ja kahden eri materiaalin vaikutusta väsymislujuuteen. Materiaaleina kokeissa oli Ruukki double grade ja Optim 700 Plus MH. Eurokoodi 3 ja IIW määrittelee kiinnityshitseille väsymiskestävyysluokaksi 125-140, kun SN -käyrän kaltevuus, m = 3. Koetuloksia vertailtiin näihin valmiiksi määriteltyihin väsymiskestävyysluokkiin. Ruukki Metals Oy toimitti koemateriaalin kahdenlaisia väsytyskokeita varten, joiden avulla väsymislujuus voitiin määrittää. Ensimmäinen koe suoritettiin värähtelyllä hyödyntäen rakenneputken ominaistaajuutta, jolloin kuormituksen rajajännityssuhde on -1. Toinen väsytyskoe suoritettiin rakenneputkista leikatuille väsytyssauvoille, missä pituussuuntainen HF –hitsi on keskellä sauvaa. Tämä väsytyskoe suoritettiin vetotykytyksellä rajajännityssuhteella 0,1. Väsytyskokeissa havaittiin, että kiinnityshitsi täyttää sille asetetut vaatimukset, koska murtumat tapahtuivat perusaineesta. Tällöin koetuloksia vertaillaan Eurokoodi 3:n perusaineen väsymisluokkaan 160, SN -käyrän kaltevuudella m = 5. Kaikkien koetulosten väsymiskestävyysluokan keskiarvoksi saatiin 185 95 %:n todennäköisyydellä. Tutkimuksessa saatiin arvokasta tietoa kahden eri materiaalin väsymislujuudesta. Materiaalilla ei kuitenkaan havaittu vaikutusta lujuuteen. Eri jälkikäsittelyt eli juuren puolen metallipurseen höyläys tai pituussuuntaisen hitsisauman hehkutus eivät vaikuttaneet väsymislujuuteen. Kokeissa väsymismurtuma ydintyi rakenneputken pinnalta olevista alkuvioista, joten putkien käsittelyä valmistuksen jälkeen voitaisiin parantaa.

Distortion and internal warping torsion of a double symmetric hollow section

This thesis is made in cooperation with Laboratory of Steel Structures and the steel company SSAB. Maximization of the benefits of high-strength steel usually requires the usage of thin wall thicknesses. This means the failures related to buckling, distortion and warping stand out. One must be aware of these phenomena to design thin-walled structures stressed with forces such as torsional loading. It is also important to take into account small stress ranges when evaluating the accurate fatigue strength of structures. The objective of this thesis is to clarify the theory of the uniform and non-uniform torsion. This paper focuses on warping due to the non-uniform torsion in double symmetric box girder and structural hollow section. The arisen stress states are explained and researched using the finite element method. Another research target is the distortion in double symmetric box girder due to torsion, and the restraining effect of transverse diaphragms at the load end. Multiple transverse diaphragms are used to study more efficient restraining against warping and distortion than a common one end plate structure.

Thick Section Laser Beam Welding of Structural Steels: Methods for Improving Welding Efficiency

Laser beam welding (LBW) is applicable for a wide range of industrial sectors and has a history of fifty years. However, it is considered an unusual method with applications typically limited to welding of thin sheet metal. With a new generation of high power lasers there has been a renewed interest in thick section LBW (also known as keyhole laser welding). There was a growing body of publications during 2001-2011 that indicates an increasing interest in laser welding for many industrial applications, and in last ten years, an increasing number of studies have examined the ways to increase the efficiency of the process. Expanding the thickness range and efficiency of LBW makes the process a possibility for industrial applications dealing with thick metal welding: shipbuilding, offshore structures, pipelines, power plants and other industries. The advantages provided by LBW, such as high process speed, high productivity, and low heat input, may revolutionize these industries and significantly reduce the process costs. The research to date has focused on either increasing the efficiency via optimizing process parameters, or on the process fundamentals, rather than on process and workpiece modifications. The argument of this thesis is that the efficiency of the laser beam process can be increased in a straightforward way in the workshop conditions. Throughout this dissertation, the term “efficiency” is used to refer to welding process efficiency, specifically, an increase in efficiency refers an increase in weld’s penetration depth without increasing laser power level or decreasing welding speed. These methods are: modifications of the workpiece – edge surface roughness and air gap between the joining plates; modification of the ambient conditions – local reduction of the pressure in the welding zone; modification of the welding process – preheating of the welding zone. Approaches to improve the efficiency are analyzed and compared both separately and combined. These experimentally proven methods confirm previous findings and contribute additional evidence which expand the opportunities for laser beam welding applications. The focus of this research was primarily on the effects of edge surface roughness preparation and pre-set air gap between the plates on weld quality and penetration depth. To date, there has been no reliable evidence that such modifications of the workpiece give a positive effect on the welding efficiency. Other methods were tested in combination with the two methods mentioned above. The most promising - combining with reduced pressure method - resulted in at least 100% increase in efficiency. The results of this thesis support the idea that joining those methods in one modified process will provide the modern engineering with a sufficient tool for many novel applications with potential benefits to a range of industries.

Anaerobic performance in ice hockey : the effect of skate blade radius of hollow

The purpose of the study was to investigate the effect of skate blade radius of hollow (ROH) on anaerobic performance, specifically during the acceleration and stopping phases of an on-ice skating test. Fifteen, male Junior B hockey players (mean age 19 y ± 1.46) were recruited to participate. On-icc testing required each participant to complete an on-ice anaerobic performance test [Reed Repeat Skate (RRS)) on three separate days. During each on-ice test, the participant's skate blades were sharpened to one of three, randomly assigned, ROH values (0.63 cm, 1.27 cm, 1.90 cm). Performance times were recorded during each RRS and used to calculate anaerobic variables [anaerobic power (W), anaerobic capacity (W), and fatigue index (s, %)). Each RRS was video recorded for the purpose of motion analysis. Video footage was imported into Peak Motus™ to measure kinematic variables of the acceleration and stopping phases. The specific variables calculated from the acceleration phase were: average velocity over 6 m (m/s), average stride length (m), and mean stride rate (strides/s). The specific variables calculated from the stopping phase were: velocity at initiation of stopping (rn/s), stopping distance (m), stopping time (s). A repeated measures ANOV A was used to assess differences in mean performance and kinematic variables across the three selected hollows. Further analysis was conducted to assess differences in trial by trial performance and kinematic variables for all hollows. The primary findings of the study suggested that skate blade ROH can have a significant effect on kinematic variables, namely stride length and stride rate during the acceleration phase and stopping distance and stopping time during the stopping phase of an on-ice anaerobic performance test. During the acceleration phase, no significant difdifferences were found in SR and SL across the three selected hollows. Mean SR on the 1.27 cm hollow was significantly slower than both the 0.63 cm and 1.90 cm hollows and SL was significantly longer when skating on the 1.27 cm hollow in comparison to the 1.90 cm hollow. During the stopping phase, stopping distance on the 0.63 cm hollow (4.12 m ± 0.14) was significantly shorter than both the 1.27 cm hollow (4.43 m ± 0.08) (p < 0.05) and the 1.90 cm ho])ow (4.35 m ± 0.12) (p < 0.05). Mean ST was also significantly shorter when stopping on the 0.63 cm hollow then both the 1.27 cm and 1.90 cm hollows. Trial by trial results clearly illustrated the affect of fatigue on kinematic variables; AV, SR, IV decreased from trial 1 to 6. There was no significant effect on anaerobic performance variables during the RRS. Altering the skate blade ROH has a significant and practical affect on accelerating and stopping performance will be discussed in this paper.