A layered finite element for reinforced concrete beams with bond-slip effects

The behaviour of reinforced concrete members is affected by the slipping of steel bars inserted in the concrete matrix. A tension-stiffening effect and crack evolution occur from the beginning of slipping; thus, the assessment of those phenomena requires the introduction of a bond-slip interaction model. This work presents a beam-layered model, including the constitutive relationships of materials and their interaction, according to the CEB-FIP Model Code 1990. To eliminate the finite element sub-division procedure, a continuous slip function is imposed into the element domain. The results are continuous descriptions of bond stress in the steel-concrete interface, as well as concrete and steel stresses along the element. (C) 2007 Elsevier Ltd. All rights reserved.

Unconstrained Finite Element for Geometrical Nonlinear Dynamics of Shells

This paper presents a positional FEM formulation to deal with geometrical nonlinear dynamics of shells. The main objective is to develop a new FEM methodology based on the minimum potential energy theorem written regarding nodal positions and generalized unconstrained vectors not displacements and rotations. These characteristics are the novelty of the present work and avoid the use of large rotation approximations. A nondimensional auxiliary coordinate system is created, and the change of configuration function is written following two independent mappings from which the strain energy function is derived. This methodology is called positional and, as far as the authors' knowledge goes, is a new procedure to approximated geometrical nonlinear structures. In this paper a proof for the linear and angular momentum conservation property of the Newmark beta algorithm is provided for total Lagrangian description. The proposed shell element is locking free for elastic stress-strain relations due to the presence of linear strain variation along the shell thickness. The curved, high-order element together with an implicit procedure to solve nonlinear equations guarantees precision in calculations. The momentum conserving, the locking free behavior, and the frame invariance of the adopted mapping are numerically confirmed by examples. Copyright (C) 2009 H. B. Coda and R. R. Paccola.

EFFECT – Efficient finite element code

The theme of this dissertation is the finite element method applied to mechanical structures. A new finite element program is developed that, besides executing different types of structural analysis, also allows the calculation of the derivatives of structural performances using the continuum method of design sensitivities analysis, with the purpose of allowing, in combination with the mathematical programming algorithms found in the commercial software MATLAB, to solve structural optimization problems. The program is called EFFECT – Efficient Finite Element Code. The object-oriented programming paradigm and specifically the C ++ programming language are used for program development. The main objective of this dissertation is to design EFFECT so that it can constitute, in this stage of development, the foundation for a program with analysis capacities similar to other open source finite element programs. In this first stage, 6 elements are implemented for linear analysis: 2-dimensional truss (Truss2D), 3-dimensional truss (Truss3D), 2-dimensional beam (Beam2D), 3-dimensional beam (Beam3D), triangular shell element (Shell3Node) and quadrilateral shell element (Shell4Node). The shell elements combine two distinct elements, one for simulating the membrane behavior and the other to simulate the plate bending behavior. The non-linear analysis capability is also developed, combining the corotational formulation with the Newton-Raphson iterative method, but at this stage is only avaiable to solve problems modeled with Beam2D elements subject to large displacements and rotations, called nonlinear geometric problems. The design sensitivity analysis capability is implemented in two elements, Truss2D and Beam2D, where are included the procedures and the analytic expressions for calculating derivatives of displacements, stress and volume performances with respect to 5 different design variables types. Finally, a set of test examples were created to validate the accuracy and consistency of the result obtained from EFFECT, by comparing them with results published in the literature or obtained with the ANSYS commercial finite element code.

Predicted properties of the superheavy elements. III. Element 115, Eka-Bismuth

Element 115 is expected to be in group V-a of the periodic table and have most stable oxidation states of I and III. The oxidation state of I, which plays a minor role in bismuth chemistry, should be a major factor in 115 chemistry. This change will arise because of the large relativistic splitting of the spherically symmetric 7p_l/2 shell from the 7P_3/2 shell. Element 115 will therefore have a single 7p_3/2 electron outside a 7p^2_1/2 closed shell. The magnitude of the first ionization energy and ionic radius suggest a chemistry similar to Tl^+. Similar considerations suggest that 115^3+ will have a chemistry similar to Bi^3+. Hydrolysis will therefore be easy and relatively strongly complexing anions of strong acids will be needed in general to effect studies of complexation chemistry. Some other properties of 115 predicted are as follows: ionization potentials I 5.2 eV, II 18.1 eV, III 27.4 eV, IV 48.5 eV, 0 \rightarrow 5^+ 159 eV; heat of sublimation, 34 kcal (g-atom)^-1; atomic radius, 2.0 A; ionic radius, 115^+ 1.5 A, 115^3+ 1.0 A; entropy, 16 cal deg^-1 (g-atom)^-l (25°); standard electrode potential 115^+ |115, -1.5 V; melting and boiling points are similar to element 113.

Numerical simulation of the punching shear behaviour of self-compacting fibre reinforced flat slabs

This paper presents the numerical simulations of the punching behaviour of centrally loaded steel fibre reinforced self-compacting concrete (SFRSCC) flat slabs. Eight half scaled slabs reinforced with different content of hooked-end steel fibres (0, 60, 75 and 90 kg/m3) and concrete strengths of 50 and 70 MPa were tested and numerically modelled. Moreover, a total of 54 three-point bending tests were carried out to assess the post-cracking flexural tensile strength. All the slabs had a relatively high conventional flexural reinforcement in order to promote the occurrence of punching failure mode. Neither of the slabs had any type of specific shear reinforcement rather than the contribution of the steel fibres. The numerical simulations were performed according to the Reissner-Mindlin theory under the finite element method framework. Regarding the classic formulation of the Reissner-Mindlin theory, in order to simulate the progressive damage induced by cracking, the shell element is discretized into layers, being assumed a plane stress state in each layer. The numerical results are, then, compared with the experimental ones and it is possible to notice that they accurately predict the experimental force-deflection relationship. The type of failure observed experimentally was also predicted in the numerical simulations.

Absoluuttisten solmukoordinaattien menetelmä kuorimaisten kappaleiden mallinnuksessa

Työn tavoitteena oli kehittää nopeasti konvergoiva kuorielementti epälineaarisesti joustavien kappaleiden analysointiin. Kuorielementti perustuu absoluuttisten solmukoordinaattien menetelmään ja se hyödyntää kaarevuuden kuvausta elastisten voimien määrityksessä. Kehitettyä elementtiä verrattiin kontinuumimekaniikalla kehitettyyn kuorielementtiin ja kaupallisen elementtimenetelmän kuorielementtiin. Yksinkertaisimman kuormitustapauksen tuloksia verrattiin teknisen taivutusteorian mukaiseen analyyttiseen ratkaisuun. Staattisten testien tulokset tässä työssä kehitetyllä kuorielementillä vastasivat hyvin kaupallisella elementtimenetelmällä saatuja tuloksia. Deformaatioiden ollessa geometrisesti lineaarisella alueella, kehitetyllä kuorielementillä saadut tulokset vastasivat paremmin sekä analyyttistä ratkaisua että kaupallisella elementtimenetelmällä saatuja tuloksia kuin aiemman kontinuumimekaniikkaan perustuvan kuorielementin tulokset. Kehitetyn kuorielementin ongelmana verrattuna kontinuumimekaniikkaan perustuvaan elementtiin on monimutkaisempi kinematiikan kuvaus. Tästä on seurauksena laskenta-ajan huomattava kasvaminen. Jatkossa kannattaisi keskittyä numeeristen ratkaisumenetelmien kehittämiseen.

Fatigue analysis of brackets with and without scallop for bridge and other deck structure using effective notch stress.

Due to functional requirement of a structural detail brackets with and without scallop are frequently used in bridges, decks, ships and offshore structure. Scallops are designed to serve as passage way for fluids, to reduce weld length and plate distortions. Moreover, scallops are used to avoid intersection of two or more welds for the fact that there is the presence of inventible inherent initial crack except for full penetrated weld and the formation of multi-axial stress state at the weld intersection. Welding all around the scallop corner increase the possibility of brittle fracture even for the case the bracket is not loaded by primary load. Avoiding of scallop will establish an initial crack in the corner if bracket is welded by fillet welds. If the two weld run pass had crossed, this would have given a 3D residual stress situation. Therefore the presences and absence of scallop necessitates the 3D FEA fatigue resistance of both types of brackets using effective notch stress approach ( ). FEMAP 10.1 with NX NASTRAN was used for the 3D FEA. The first and main objective of this research was to investigate and compare the fatigue resistance of brackets with and without scallop. The secondary goal was the fatigue design of scallops in case they cannot be avoided for some reason. The fatigue resistance for both types of brackets was determined based on approach using 1 mm fictitiously rounded radius based on IIW recommendation. Identical geometrical, boundary and loading conditions were used for the determination and comparison of fatigue resistance of both types of brackets using linear 3D FEA. Moreover the size effect of bracket length was also studied using 2D SHELL element FEA. In the case of brackets with scallop the flange plate weld toe at the corner of the scallop was found to exhibit the highest and made the flange plate weld toe critical for fatigue failure. Whereas weld root and weld toe at the weld intersections were the highly stressed location for brackets without scallop. Thus weld toe for brackets with scallop, and weld root and weld toe for brackets without scallop were found to be the critical area for fatigue failure. Employing identical parameters on both types of brackets, brackets without scallop had the highest except for full penetrated weld. Furthermore the fatigue resistance of brackets without scallop was highly affected by the lack of weld penetration length and it was found out that decreased as the weld penetration was increased. Despite the fact that the very presence of scallop reduces the stiffness and also same time induce stress concentration, based on the 3D FEA it is worth concluding that using scallop provided better fatigue resistance when both types of brackets were fillet welded. However brackets without scallop had the highest fatigue resistance when full penetration weld was used. This thesis also showed that weld toe for brackets with scallop was the only highly stressed area unlike brackets without scallop in which both weld toe and weld root were the critical locations for fatigue failure when different types of boundary conditions were used. Weld throat thickness, plate thickness, scallop radius, lack of weld penetration length, boundary condition and weld quality affected the fatigue resistance of both types of brackets. And as a result, bracket design procedure, especially welding quality and post weld treatment techniques significantly affect the fatigue resistance of both type of brackets.

Performance Assessment of Sandwich Structures with Debonds and Dents

A sandwich construction is a special form of the laminated composite consisting of light weight core, sandwiched between two stiff thin face sheets. Due to high stiffness to weight ratio, sandwich construction is widely adopted in aerospace industries. As a process dependent bonded structure, the most severe defects associated with sandwich construction are debond (skin core bond failure) and dent (locally deformed skin associated with core crushing). Reasons for debond may be attributed to initial manufacturing flaws or in service loads and dent can be caused by tool drops or impacts by foreign objects. This paper presents an evaluation on the performance of honeycomb sandwich cantilever beam with the presence of debond or dent, using layered finite element models. Dent is idealized by accounting core crushing in the core thickness along with the eccentricity of the skin. Debond is idealized using multilaminate modeling at debond location with contact element between the laminates. Vibration and buckling behavior of metallic honeycomb sandwich beam with and without damage are carried out. Buckling load factor, natural frequency, mode shape and modal strain energy are evaluated using finite element package ANSYS 13.0. Study shows that debond affect the performance of the structure more severely than dent. Reduction in the fundamental frequencies due to the presence of dent or debond is not significant for the case considered. But the debond reduces the buckling load factor significantly. Dent of size 8-20% of core thickness shows 13% reduction in buckling load capacity of the sandwich column. But debond of the same size reduced the buckling load capacity by about 90%. This underscores the importance of detecting these damages in the initiation level itself to avoid catastrophic failures. Influence of the damages on fundamental frequencies, mode shape and modal strain energy are examined. Effectiveness of these parameters as a damage detection tool for sandwich structure is also assessed

Ontologies and Methods for Interoperability of Engineering Analysis Models (EAMS) in an E-Design Environment

ABSTRACT ONTOLOGIES AND METHODS FOR INTEROPERABILITY OF ENGINEERING ANALYSIS MODELS (EAMS) IN AN E-DESIGN ENVIRONMENT SEPTEMBER 2007 NEELIMA KANURI, B.S., BIRLA INSTITUTE OF TECHNOLOGY AND SCIENCES PILANI INDIA M.S., UNIVERSITY OF MASSACHUSETTS AMHERST Directed by: Professor Ian Grosse Interoperability is the ability of two or more systems to exchange and reuse information efficiently. This thesis presents new techniques for interoperating engineering tools using ontologies as the basis for representing, visualizing, reasoning about, and securely exchanging abstract engineering knowledge between software systems. The specific engineering domain that is the primary focus of this report is the modeling knowledge associated with the development of engineering analysis models (EAMs). This abstract modeling knowledge has been used to support integration of analysis and optimization tools in iSIGHT FD , a commercial engineering environment. ANSYS , a commercial FEA tool, has been wrapped as an analysis service available inside of iSIGHT-FD. Engineering analysis modeling (EAM) ontology has been developed and instantiated to form a knowledge base for representing analysis modeling knowledge. The instances of the knowledge base are the analysis models of real world applications. To illustrate how abstract modeling knowledge can be exploited for useful purposes, a cantilever I-Beam design optimization problem has been used as a test bed proof-of-concept application. Two distinct finite element models of the I-beam are available to analyze a given beam design- a beam-element finite element model with potentially lower accuracy but significantly reduced computational costs and a high fidelity, high cost, shell-element finite element model. The goal is to obtain an optimized I-beam design at minimum computational expense. An intelligent KB tool was developed and implemented in FiPER . This tool reasons about the modeling knowledge to intelligently shift between the beam and the shell element models during an optimization process to select the best analysis model for a given optimization design state. In addition to improved interoperability and design optimization, methods are developed and presented that demonstrate the ability to operate on ontological knowledge bases to perform important engineering tasks. One such method is the automatic technical report generation method which converts the modeling knowledge associated with an analysis model to a flat technical report. The second method is a secure knowledge sharing method which allocates permissions to portions of knowledge to control knowledge access and sharing. Both the methods acting together enable recipient specific fine grain controlled knowledge viewing and sharing in an engineering workflow integration environment, such as iSIGHT-FD. These methods together play a very efficient role in reducing the large scale inefficiencies existing in current product design and development cycles due to poor knowledge sharing and reuse between people and software engineering tools. This work is a significant advance in both understanding and application of integration of knowledge in a distributed engineering design framework.

Analysis of Combined Bending and Punching Tests of Reinforced Concrete Slabs within IMPACT III Project

The paper reports on a collaborative effort between the Swiss Federal Nuclear Safety Inspectorate (ENSI) and their consultants Principia and Stangenberg. As part of the IMPACT III project, reduced scale impact tests of reinforced concrete structures were carried out. The simulation of test X3 is presented here and the numerical results are compared with those obtained in the test, carried out in August 2013. The general object is to improve the safety of nuclear facilities and, more specifically, to demonstrate the capabilities of current simulation techniques to reproduce the behaviour of a reinforced concrete structure impacted by a soft missile. The missile is a steel tube with a mass of 50 kg and travelling at 140 m/s. The target is a 250 mm thick, 2,1 m by 2,1 m reinforced concrete wall, held in a stiff supporting frame. The reinforcement includes both longitudinal and transverse rebars. Calculations were carried out before and after the test with Abaqus (Principia) and SOFiSTiK (Stangenberg). In the Abaqus simulation the concrete is modelled using solid elements and a damaged plasticity formulation, the rebars with embedded beam elements, and the missile with shell elements. In SOFiSTiK the target is modelled with non-linear, layered shell elements for the reinforcement on both sides; non-linear shear deformations of shell/plate elements are approximately included. The results generally indicate a good agreement between calculations and measurements.

Suppression of spurious intermediate frequency modes in under-integrated elements by combined stiffness/viscous stabilization

Solutions employing perturbation stiffness or viscous hourglass control with one-point quadrature finite elements often exhibit spurious modes in the intermediate frequency range. These spurious frequencies are demonstrated in several examples and their origin is explained. Then it is shown that by critically damping the hourglass modes, these spurious mid-range frequency modes can be suppressed. Estimates of the hourglass frequency and damping coefficients are provided for the plane 4-node quadrilateral and a 4-node shell element. Results are presented that show almost complete annihilation of spurious intermediate frequency modes for both linear and non-linear problems. Copyright (c) 2005 John Wiley & Sons, Ltd.

Stress Concentration Factor and Fatigue Life Evaluation in offshore tubular KT-Joints

In the last few decades, offshore field has grown fast especially after the notable development of technologies, explorations of oil and gas in deep water and the high concern of offshore companies in renewable energy mainly Wind Energy. Fatigue damage was noticed as one of the main problems causing failure of offshore structures. The purpose of this research is to focus on the evaluation of Stress Concentration Factor and its influence on Fatigue Life for 2 tubular KT-Joints in offshore Jacket structure using different calculation methods. The work is done by using analytical calculations, mainly Efthymiou’s formulations, and numerical solutions, FEM analysis, using ABAQUS software. As for the analytical formulations, the calculations were done according to the geometrical parameters of each method using excel sheets. As for the numerical model, 2 different types of tubular KT-Joints are present where for each model 5 shell element type, 3 solid element type and 3 solid-with-weld element type models were built on ABAQUS. Meshing was assigned according to International Institute of Welding (IIW) recommendations, 5 types of mesh element, to evaluate the Hot-spot stresses. 23 different types of unitary loading conditions were assigned, 9 axial, 7 in-plane bending moment and 7 out-plane bending moment loads. The extraction of Hot-spot stresses and the evaluation of the Stress Concentration Factor were done using PYTHON scripting and MATLAB. Then, the fatigue damage evaluation for a critical KT tubular joint based on Simplified Fatigue Damage Rule and Local Approaches (Strain Damage Parameter and Stress Damage Parameter) methods were calculated according to the maximum Stress Concentration Factor conducted from DNV and FEA methods. In conclusion, this research helped us to compare different results of Stress Concentration Factor and Fatigue Life using different methods and provided us with a general overview about what to study next in the future.

Shell curvature as an initial deformation: A geometrically exact finite element approach

An alternative approach for the analysis of arbitrarily curved shells is developed in this paper based on the idea of initial deformations. By `alternative` we mean that neither differential geometry nor the concept of degeneration is invoked here to describe the shell surface. We begin with a flat reference configuration for the shell mid-surface, after which the initial (curved) geometry is mapped as a stress-free deformation from the plane position. The actual motion of the shell takes place only after this initial mapping. In contrast to classical works in the literature, this strategy enables the use of only orthogonal frames within the theory and therefore objects such as Christoffel symbols, the second fundamental form or three-dimensional degenerated solids do not enter the formulation. Furthermore, the issue of physical components of tensors does not appear. Another important aspect (but not exclusive of our scheme) is the possibility to describe exactly the initial geometry. The model is kinematically exact, encompasses finite strains in a totally consistent manner and is here discretized under the light of the finite element method (although implementation via mesh-free techniques is also possible). Assessment is made by means of several numerical simulations. Copyright (C) 2009 John Wiley & Sons, Ltd.

Finite element analysis of heat transfer and mineralization in layered hydrothermal systems with upward throughflow

We use the finite element method to solve the coupled problem between convective pore-fluid flow, heat transfer and mineralization in layered hydrothermal systems with upward throughflow. In particular, we present the improved rock alteration index (IRAI) concept for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in the systems. To validate the numerical method used in the computation, analytical solutions to a benchmark problem have been derived. After the numerical method is validated, it is used to investigate the pattern of pore-fluid Aom, the distribution of temperature and the mineralization pattern of gold minerals in a layered hydrothermal system with upward throughflow. The related numerical results have demonstrated that the present concept of IRAI is useful and applicable for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in hydrothermal systems. (C) 2000 Elsevier Science S.A. All rights reserved.

U-Pb-Hf-trace element systematics and geochronology of zircon from a granulite-facies metamorphosed mafic-ultramafic layered complex in Central Brazil

Zircon recrystallization is a common process during high-grade metamorphism and promotes partial or complete resetting of the original isotopic and chemical characteristics of the mineral and thus complicates U-Pb geochronological interpretation. In Central Brazil, this complexity may be illustrated by three composite mafic-ultramafic intrusions metamorphosed under amphibolite-to-granulite conditions. Their ages of emplacement and metamorphic ages have been a matter of controversy for the last thirty years. The Serra da Malacacheta and Barro Alto complexes make up the southernmost of these layered bodies and four samples from distinct rock types were investigated in order to verify the consequences of metamorphic alteration of zircon for U-Pb dating. Cathodoluminescent imaging reveals internal features which are typical of concomitant dissolution-reprecipitation processes, such as convolute zoning and inward-moving recrystallization fronts, even in samples in which partially preserved igneous textures are observed. Due to this extensive alteration, LA-ICPMS U-Pb isotopic analysis yielded inconclusive data. However, in situ Hf isotopic and trace-element analyses help to clarify the real meaning of the geochronological data. Low Lu/Hf (<0.004) and homogeneous (176)Hf/(177)Hf(t) values imply that the zircon populations within individual samples have crystallized in a single episode, despite the observed variations in age values. Trace element signatures of zircon grains from garnet-bearing samples reveal that they were unreactive to the development of the peak metamorphism mineral assemblage and, thus, the main chemical feature in such grains is attributed to a coupled dissolution-reprecipitation process. However, in the Cafelandia amphibolite an additional alteration process is identified, probably related to the influx of late-stage fluids. Combined isotopic and geochemical investigation on zircon grains allowed the distinction of two magmatic events. The first corresponds to the crystallization of the Serra da Malacacheta Complex and characterizes a juvenile magmatism at similar to 1.3 Ga. The younger episode, recognized in the Barro Alto Complex, is dated at ca. 800 Ma and is represented by mafic and ultramafic rocks showing intense contamination with continental crust, implying that the emplacement took place, most likely, in a continental back-arc setting. Altered zircon domains as well as titanite grains date the metamorphic event at ca. 760-750 Ma. (C) 2011 Elsevier B.V. All rights reserved.