FLASH: A finite element computer code for variably saturated flow

International Nuclear Information System (INIS)

Baca, R.G.; Magnuson, S.O.

1992-05-01

A numerical model was developed for use in performance assessment studies at the INEL. The numerical model, referred to as the FLASH computer code, is designed to simulate two-dimensional fluid flow in fractured-porous media. The code is specifically designed to model variably saturated flow in an arid site vadose zone and saturated flow in an unconfined aquifer. In addition, the code also has the capability to simulate heat conduction in the vadose zone. This report presents the following: description of the conceptual frame-work and mathematical theory; derivations of the finite element techniques and algorithms; computational examples that illustrate the capability of the code; and input instructions for the general use of the code. The FLASH computer code is aimed at providing environmental scientists at the INEL with a predictive tool for the subsurface water pathway. This numerical model is expected to be widely used in performance assessments for: (1) the Remedial Investigation/Feasibility Study process and (2) compliance studies required by the US Department of Energy Order 5820.2A

A first course in finite elements

CERN Document Server

Fish, Jacob

2007-01-01

Developed from the authors, combined total of 50 years undergraduate and graduate teaching experience, this book presents the finite element method formulated as a general-purpose numerical procedure for solving engineering problems governed by partial differential equations.  Focusing on the formulation and application of the finite element method through the integration of finite element theory, code development, and software application, the book is both introductory and self-contained, as well as being a hands-on experience for any student. This authoritative text on Finite Elements:Adopts

Industrial applications of N3S finite element code

International Nuclear Information System (INIS)

Chabard, J.P.; Pot, G.; Martin, A.

1993-12-01

The Research and Development Division of EDF (French utilities) has been working since 1982 on N3S, a 3D finite element code for simulating turbulent incompressible flows (Chabard et al., 1992) which has many applications nowadays dealing with internal flows, thermal hydraulics (Delenne and Pot, 1993), turbomachinery (Combes and Rieutord, 1992). The size of these applications is larger and larger: calculations until 350 000 nodes are in progress (around 2 000 000 unknowns). To achieve so large applications, an important work has been done on the choice of efficient algorithms and on their implementation in order to reduce CPU time and memory allocation. The paper presents the central algorithm of the code, focusing on time and memory optimization. As an illustration, validation test cases and a recent industrial application are discussed. (authors). 11 figs., 2 tabs., 11 refs

Massively Parallel Finite Element Programming

KAUST Repository

Heister, Timo; Kronbichler, Martin; Bangerth, Wolfgang

2010-01-01

Today's large finite element simulations require parallel algorithms to scale on clusters with thousands or tens of thousands of processor cores. We present data structures and algorithms to take advantage of the power of high performance computers in generic finite element codes. Existing generic finite element libraries often restrict the parallelization to parallel linear algebra routines. This is a limiting factor when solving on more than a few hundreds of cores. We describe routines for distributed storage of all major components coupled with efficient, scalable algorithms. We give an overview of our effort to enable the modern and generic finite element library deal.II to take advantage of the power of large clusters. In particular, we describe the construction of a distributed mesh and develop algorithms to fully parallelize the finite element calculation. Numerical results demonstrate good scalability. © 2010 Springer-Verlag.

Massively Parallel Finite Element Programming

KAUST Repository

Heister, Timo

2010-01-01

Today\\'s large finite element simulations require parallel algorithms to scale on clusters with thousands or tens of thousands of processor cores. We present data structures and algorithms to take advantage of the power of high performance computers in generic finite element codes. Existing generic finite element libraries often restrict the parallelization to parallel linear algebra routines. This is a limiting factor when solving on more than a few hundreds of cores. We describe routines for distributed storage of all major components coupled with efficient, scalable algorithms. We give an overview of our effort to enable the modern and generic finite element library deal.II to take advantage of the power of large clusters. In particular, we describe the construction of a distributed mesh and develop algorithms to fully parallelize the finite element calculation. Numerical results demonstrate good scalability. © 2010 Springer-Verlag.

Evaluation of finite element codes for demonstrating the performance of radioactive material packages in hypothetical accident drop scenarios

International Nuclear Information System (INIS)

Tso, C.F.; Hueggenberg, R.

2004-01-01

Drop testing and analysis are the two methods for demonstrating the performance of packages in hypothetical drop accident scenarios. The exact purpose of the tests and the analyses, and the relative prominence of the two in the license application, may depend on the Competent Authority and will vary between countries. The Finite Element Method (FEM) is a powerful analysis tool. A reliable finite element (FE) code when used correctly and appropriately, will allow a package's behaviour to be simulated reliably. With improvements in computing power, and in sophistication and reliability of FE codes, it is likely that FEM calculations will increasingly be used as evidence of drop test performance when seeking Competent Authority approval. What is lacking at the moment, however, is a standardised method of assessing a FE code in order to determine whether it is sufficiently reliable or pessimistic. To this end, the project Evaluation of Codes for Analysing the Drop Test Performance of Radioactive Material Transport Containers, funded by the European Commission Directorate-General XVII (now Directorate-General for Energy and Transport) and jointly performed by Arup and Gesellschaft fuer Nuklear-Behaelter mbH, was carried out in 1998. The work consisted of three components: Survey of existing finite element software, with a view to finding codes that may be capable of analysing drop test performance of radioactive material packages, and to produce an inventory of them. Develop a set of benchmark problems to evaluate software used for analysing the drop test performance of packages. Evaluate the finite element codes by testing them against the benchmarks This paper presents a summary of this work

Evaluation of finite element codes for demonstrating the performance of radioactive material packages in hypothetical accident drop scenarios

Energy Technology Data Exchange (ETDEWEB)

Tso, C.F. [Arup (United Kingdom); Hueggenberg, R. [Gesellschaft fuer Nuklear-Behaelter mbH (Germany)

2004-07-01

Drop testing and analysis are the two methods for demonstrating the performance of packages in hypothetical drop accident scenarios. The exact purpose of the tests and the analyses, and the relative prominence of the two in the license application, may depend on the Competent Authority and will vary between countries. The Finite Element Method (FEM) is a powerful analysis tool. A reliable finite element (FE) code when used correctly and appropriately, will allow a package's behaviour to be simulated reliably. With improvements in computing power, and in sophistication and reliability of FE codes, it is likely that FEM calculations will increasingly be used as evidence of drop test performance when seeking Competent Authority approval. What is lacking at the moment, however, is a standardised method of assessing a FE code in order to determine whether it is sufficiently reliable or pessimistic. To this end, the project Evaluation of Codes for Analysing the Drop Test Performance of Radioactive Material Transport Containers, funded by the European Commission Directorate-General XVII (now Directorate-General for Energy and Transport) and jointly performed by Arup and Gesellschaft fuer Nuklear-Behaelter mbH, was carried out in 1998. The work consisted of three components: Survey of existing finite element software, with a view to finding codes that may be capable of analysing drop test performance of radioactive material packages, and to produce an inventory of them. Develop a set of benchmark problems to evaluate software used for analysing the drop test performance of packages. Evaluate the finite element codes by testing them against the benchmarks This paper presents a summary of this work.

Finite element analysis of piezoelectric materials

International Nuclear Information System (INIS)

Lowrie, F.; Stewart, M.; Cain, M.; Gee, M.

1999-01-01

This guide is intended to help people wanting to do finite element analysis of piezoelectric materials by answering some of the questions that are peculiar to piezoelectric materials. The document is not intended as a complete beginners guide for finite element analysis in general as this is better dealt with by the individual software producers. The guide is based around the commercial package ANSYS as this is a popular package amongst piezoelectric material users, however much of the information will still be useful to users of other finite element codes. (author)

FEHMN 1.0: Finite element heat and mass transfer code

International Nuclear Information System (INIS)

Zyvoloski, G.; Dash, Z.; Kelkar, S.

1991-04-01

A computer code is described which can simulate non-isothermal multiphase multicomponent flow in porous media. It is applicable to natural-state studies of geothermal systems and ground-water flow. The equations of heat and mass transfer for multiphase flow in porous and permeable media are solved using the finite element method. The permeability and porosity of the medium are allowed to depend on pressure and temperature. The code also has provisions for movable air and water phases and noncoupled tracers; that is, tracer solutions that do not affect the heat and mass transfer solutions. The tracers can be passive or reactive. The code can simulate two-dimensional, two-dimensional radial, or three-dimensional geometries. A summary of the equations in the model and the numerical solution procedure are provided in this report. A user's guide and sample problems are also included. The main use of FEHMN will be to assist in the understanding of flow fields in the saturated zone below the proposed Yucca Mountain Repository. 33 refs., 27 figs., 12 tabs

Modeling turbine-missile impacts using the HONDO finite-element code

International Nuclear Information System (INIS)

Schuler, K.W.

1981-11-01

Calculations have been performed using the dynamic finite element code HONDO to simulate a full scale rocket sled test. In the test a rocket sled was used to launch at a velocity of 150 m/s (490 ft/s), a 1527 kg (3366 lb) fragment of a steam turbine rotor disk into a structure which was a simplified model of a steam turbine casing. In the calculations the material behavior of and boundary conditions on the target structure were varied to assess its energy absorbing characteristics. Comparisons are made between the calculations and observations of missile velocity and strain histories of various points of the target structure

Axisym finite element code: modifications for pellet-cladding mechanical interaction analysis

International Nuclear Information System (INIS)

Engelman, G.P.

1978-10-01

Local strain concentrations in nuclear fuel rods are known to be potential sites for failure initiation. Assessment of such strain concentrations requires a two-dimensional analysis of stress and strain in both the fuel and the cladding during pellet-cladding mechanical interaction. To provide such a capability in the FRAP (Fuel Rod Analysis Program) codes, the AXISYM code (a small finite element program developed at the Idaho National Engineering Laboratory) was modified to perform a detailed fuel rod deformation analysis. This report describes the modifications which were made to the AXISYM code to adapt it for fuel rod analysis and presents comparisons made between the two-dimensional AXISYM code and the FRACAS-II code. FRACAS-II is the one-dimensional (generalized plane strain) fuel rod mechanical deformation subcode used in the FRAP codes. Predictions of these two codes should be comparable away from the fuel pellet free ends if the state of deformation at the pellet midplane is near that of generalized plane strain. The excellent agreement obtained in these comparisons checks both the correctness of the AXISYM code modifications as well as the validity of the assumption of generalized plane strain upon which the FRACAS-II subcode is based

ABAQUS/EPGEN - a general purpose finite element code with emphasis on nonlinear applications

International Nuclear Information System (INIS)

Hibbitt, H.D.

1984-01-01

The article contains a summary description of ABAQUS, a finite element program designed for general use in nonlinear as well as linear structural problems, in the context of its application to nuclear structural integrity analysis. The article begins with a discussion of the design criteria and methods upon which the code development has been based. The engineering modelling capabilities, currently implemented in the program - elements, constitutive models and analysis procedures - are then described. Finally, a few demonstration examples are presented, to illustrate some of the program's features that are of interest in structural integrity analysis associated with nuclear power plants. (orig.)

Finite Element Modelling of Seismic Liquefaction in Soils

NARCIS (Netherlands)

Galavi, V.; Petalas, A.; Brinkgreve, R.B.J.

2013-01-01

Numerical aspects of seismic liquefaction in soils as implemented in the finite element code, PLAXIS, is described in this paper. After description of finite element equations of dynamic problems, three practical dynamic boundary conditions, namely viscous boundary tractions, tied degrees of freedom

Modelling 3-D mechanical phenomena in a 1-D industrial finite element code: results and perspectives

International Nuclear Information System (INIS)

Guicheret-Retel, V.; Trivaudey, F.; Boubakar, M.L.; Masson, R.; Thevenin, Ph.

2005-01-01

Assessing fuel rod integrity in PWR reactors must enjoin two opposite goals: a one-dimensional finite element code (axial revolution symmetry) is needed to provide industrial results at the scale of the reactor core, while the main risk of cladding failure [e.g. pellet-cladding interaction (PCI)] is based on fully three-dimensional phenomena. First, parametric three-dimensional elastic calculations were performed to identify the relevant parameters (fragment number, contact pellet-cladding conditions, etc.) as regards PCI. Axial fragment number as well as friction coefficient are shown to play a major role in PCI as opposed to other parameters. Next, the main limitations of the one-dimensional hypothesis of the finite element code CYRANO3 are identified. To overcome these limitations, both two- and three-dimensional emulations of CYRANO3 were developed. These developments are shown to significantly improve the results provided by CYRANO3. (authors)

The finite element response matrix method

International Nuclear Information System (INIS)

Nakata, H.; Martin, W.R.

1983-02-01

A new technique is developed with an alternative formulation of the response matrix method implemented with the finite element scheme. Two types of response matrices are generated from the Galerkin solution to the weak form of the diffusion equation subject to an arbitrary current and source. The piecewise polynomials are defined in two levels, the first for the local (assembly) calculations and the second for the global (core) response matrix calculations. This finite element response matrix technique was tested in two 2-dimensional test problems, 2D-IAEA benchmark problem and Biblis benchmark problem, with satisfatory results. The computational time, whereas the current code is not extensively optimized, is of the same order of the well estabilished coarse mesh codes. Furthermore, the application of the finite element technique in an alternative formulation of response matrix method permits the method to easily incorporate additional capabilities such as treatment of spatially dependent cross-sections, arbitrary geometrical configurations, and high heterogeneous assemblies. (Author) [pt

Mechanical modelization of PCI with Fragema and CEA finite element codes

International Nuclear Information System (INIS)

Bernard, P.; Joseph, J.; Atabek, R.; Chantant, M.

1982-03-01

In order to modelize the PCI phenomenon during a power ramp test two finite element codes have been used by FRAGEMA and CEA, TITUS and VERDON. The results given, by the 3D equivalent method developed with TITUS, and VERDON are equivalent, in particular the strains and the equivalent Von Mises stresses at the pellet to pellet interface are quite similar. An evaluation was made to explain experimental ramp tests results. These results come from FRISCA 04bis and FRISCA 104 rods which were ramp tested in SILOE. The choice of the equivalent Von Mises stress seems to be quite a good criterion to explain the failure threshold

SQA of finite element method (FEM) codes used for analyses of pit storage/transport packages

Energy Technology Data Exchange (ETDEWEB)

Russel, E. [Lawrence Livermore National Lab., CA (United States)

1997-11-01

This report contains viewgraphs on the software quality assurance of finite element method codes used for analyses of pit storage and transport projects. This methodology utilizes the ISO 9000-3: Guideline for application of 9001 to the development, supply, and maintenance of software, for establishing well-defined software engineering processes to consistently maintain high quality management approaches.

A finite element conjugate gradient FFT method for scattering

Science.gov (United States)

Collins, Jeffery D.; Ross, Dan; Jin, J.-M.; Chatterjee, A.; Volakis, John L.

1991-01-01

Validated results are presented for the new 3D body of revolution finite element boundary integral code. A Fourier series expansion of the vector electric and mangnetic fields is employed to reduce the dimensionality of the system, and the exact boundary condition is employed to terminate the finite element mesh. The mesh termination boundary is chosen such that is leads to convolutional boundary operatores of low O(n) memory demand. Improvements of this code are discussed along with the proposed formulation for a full 3D implementation of the finite element boundary integral method in conjunction with a conjugate gradiant fast Fourier transformation (CGFFT) solution.

FLAME: A finite element computer code for contaminant transport n variably-saturated media

International Nuclear Information System (INIS)

Baca, R.G.; Magnuson, S.O.

1992-06-01

A numerical model was developed for use in performance assessment studies at the INEL. The numerical model referred to as the FLAME computer code, is designed to simulate subsurface contaminant transport in a variably-saturated media. The code can be applied to model two-dimensional contaminant transport in an and site vadose zone or in an unconfined aquifer. In addition, the code has the capability to describe transport processes in a porous media with discrete fractures. This report presents the following: description of the conceptual framework and mathematical theory, derivations of the finite element techniques and algorithms, computational examples that illustrate the capability of the code, and input instructions for the general use of the code. The development of the FLAME computer code is aimed at providing environmental scientists at the INEL with a predictive tool for the subsurface water pathway. This numerical model is expected to be widely used in performance assessments for: (1) the Remedial Investigation/Feasibility Study process and (2) compliance studies required by the US Department of energy Order 5820.2A

FLAME: A finite element computer code for contaminant transport n variably-saturated media

Energy Technology Data Exchange (ETDEWEB)

Baca, R.G.; Magnuson, S.O.

1992-06-01

A numerical model was developed for use in performance assessment studies at the INEL. The numerical model referred to as the FLAME computer code, is designed to simulate subsurface contaminant transport in a variably-saturated media. The code can be applied to model two-dimensional contaminant transport in an and site vadose zone or in an unconfined aquifer. In addition, the code has the capability to describe transport processes in a porous media with discrete fractures. This report presents the following: description of the conceptual framework and mathematical theory, derivations of the finite element techniques and algorithms, computational examples that illustrate the capability of the code, and input instructions for the general use of the code. The development of the FLAME computer code is aimed at providing environmental scientists at the INEL with a predictive tool for the subsurface water pathway. This numerical model is expected to be widely used in performance assessments for: (1) the Remedial Investigation/Feasibility Study process and (2) compliance studies required by the US Department of energy Order 5820.2A.

PLASTEF: a code for the numerical simulation of thermoelastoplastic behaviour of materials using the finite element method

International Nuclear Information System (INIS)

Basombrio, F.G.; Sanchez Sarmiento, G.

1978-01-01

A general code for solving two-dimensional thermo-elastoplastic problems in geometries of arbitrary shape using the finite element method, is presented. The initial stress incremental procedure was adopted, for given histories of load and temperature. Some classical applications are included. (Auth.)

FINELM: a multigroup finite element diffusion code

International Nuclear Information System (INIS)

Higgs, C.E.; Davierwalla, D.M.

1981-06-01

FINELM is a FORTRAN IV program to solve the Neutron Diffusion Equation in X-Y, R-Z, R-theta, X-Y-Z and R-theta-Z geometries using the method of Finite Elements. Lagrangian elements of linear or higher degree to approximate the spacial flux distribution have been provided. The method of dissections, coarse mesh rebalancing and Chebyshev acceleration techniques are available. Simple user defined input is achieved through extensive input subroutines. The input preparation is described followed by a program structure description. Sample test cases are provided. (Auth.)

Large Eddy Simulation of turbulent flows in compound channels with a finite element code

International Nuclear Information System (INIS)

Xavier, C.M.; Petry, A.P.; Moeller, S.V.

2011-01-01

This paper presents the numerical investigation of the developing flow in a compound channel formed by a rectangular main channel and a gap in one of the sidewalls. A three dimensional Large Eddy Simulation computational code with the classic Smagorinsky model is introduced, where the transient flow is modeled through the conservation equations of mass and momentum of a quasi-incompressible, isothermal continuous medium. Finite Element Method, Taylor-Galerkin scheme and linear hexahedrical elements are applied. Numerical results of velocity profile show the development of a shear layer in agreement with experimental results obtained with Pitot tube and hot wires. (author)

An introduction to the UNCLE finite element scheme

International Nuclear Information System (INIS)

Enderby, J.A.

1983-01-01

UNCLE is a completely general finite element scheme which provides common input, output, equation-solving and other facilities for a family of finite element codes for linear and non-linear stress analysis, heat transfer etc. This report describes the concepts on which UNCLE is based and gives a general account of the facilities provided. (author)

An introduction to the UNCLE finite element scheme

Energy Technology Data Exchange (ETDEWEB)

Enderby, J A [UK Atomic Energy Authority, Northern Division, Risley Nuclear Power Development Establishment, Risley, Warrington (United Kingdom)

1983-05-01

UNCLE is a completely general finite element scheme which provides common input, output, equation-solving and other facilities for a family of finite element codes for linear and non-linear stress analysis, heat transfer etc. This report describes the concepts on which UNCLE is based and gives a general account of the facilities provided. (author)

Implementation of advanced finite element technology in structural analysis computer codes

International Nuclear Information System (INIS)

Kohli, T.D.; Wiley, J.W.; Koss, P.W.

1975-01-01

Advances in finite element technology over the last several years have been rapid and have largely outstripped the ability of general purpose programs in the public domain to assimilate them. As a result, it has become the burden of the structural analyst to incorporate these advances himself. This paper discusses the implementation and extension of specific technological advances in Bechtel structural analysis programs. In general these advances belong in two categories: (1) the finite elements themselves and (2) equation solution algorithms. Improvements in the finite elements involve increased accuracy of the elements and extension of their applicability to various specialized modelling situations. Improvements in solution algorithms have been almost exclusively aimed at expanding problem solving capacity. (Auth.)

Numerical experiments in finite element analysis of thermoelastoplastic behaviour of materials. Further developments of the PLASTEF code

International Nuclear Information System (INIS)

Basombrio, F.G.; Sarmiento, G.S.

1980-01-01

In a previous paper the finite element code PLASTEF for the numerical simulation of thermoelastoplastic behaviour of materials was presented in its general outline. This code employs an initial stress incremental procedure for given histories of loads and temperature. It has been formulated for medium sized computers. The present work is an extension of the previous paper to consider additional aspects of the variable temperature case. Non-trivial tests of this type of situation are described. Finally, details are given of some concrete applications to the prediction of thermoelastoplastic collapse of nuclear fuel element cladding. (author)

Finite element solution of two dimensional time dependent heat equation

International Nuclear Information System (INIS)

Maaz

1999-01-01

A Microsoft Windows based computer code, named FHEAT, has been developed for solving two dimensional heat problems in Cartesian and Cylindrical geometries. The programming language is Microsoft Visual Basic 3.0. The code makes use of Finite element formulation for spatial domain and Finite difference formulation for time domain. Presently the code is capable of solving two dimensional steady state and transient problems in xy- and rz-geometries. The code is capable excepting both triangular and rectangular elements. Validation and benchmarking was done against hand calculations and published results. (author)

Analysis of lower head failure with simplified models and a finite element code

Energy Technology Data Exchange (ETDEWEB)

Koundy, V. [CEA-IPSN-DPEA-SEAC, Service d' Etudes des Accidents, Fontenay-aux-Roses (France); Nicolas, L. [CEA-DEN-DM2S-SEMT, Service d' Etudes Mecaniques et Thermiques, Gif-sur-Yvette (France); Combescure, A. [INSA-Lyon, Lab. Mecanique des Solides, Villeurbanne (France)

2001-07-01

The objective of the OLHF (OECD lower head failure) experiments is to characterize the timing, mode and size of lower head failure under high temperature loading and reactor coolant system pressure due to a postulated core melt scenario. Four tests have been performed at Sandia National Laboratories (USA), in the frame of an OECD project. The experimental results have been used to develop and validate predictive analysis models. Within the framework of this project, several finite element calculations were performed. In parallel, two simplified semi-analytical methods were developed in order to get a better understanding of the role of various parameters on the creep phenomenon, e.g. the behaviour of the lower head material and its geometrical characteristics on the timing, mode and location of failure. Three-dimensional modelling of crack opening and crack propagation has also been carried out using the finite element code Castem 2000. The aim of this paper is to present the two simplified semi-analytical approaches and to report the status of the 3D crack propagation calculations. (authors)

A non-linear, finite element, heat conduction code to calculate temperatures in solids of arbitrary geometry

International Nuclear Information System (INIS)

Tayal, M.

1987-01-01

Structures often operate at elevated temperatures. Temperature calculations are needed so that the design can accommodate thermally induced stresses and material changes. A finite element computer called FEAT has been developed to calculate temperatures in solids of arbitrary shapes. FEAT solves the classical equation for steady state conduction of heat. The solution is obtained for two-dimensional (plane or axisymmetric) or for three-dimensional problems. Gap elements are use to simulate interfaces between neighbouring surfaces. The code can model: conduction; internal generation of heat; prescribed convection to a heat sink; prescribed temperatures at boundaries; prescribed heat fluxes on some surfaces; and temperature-dependence of material properties like thermal conductivity. The user has a option of specifying the detailed variation of thermal conductivity with temperature. For convenience to the nuclear fuel industry, the user can also opt for pre-coded values of thermal conductivity, which are obtained from the MATPRO data base (sponsored by the U.S. Nuclear Regulatory Commission). The finite element method makes FEAT versatile, and enables it to accurately accommodate complex geometries. The optional link to MATPRO makes it convenient for the nuclear fuel industry to use FEAT, without loss of generality. Special numerical techniques make the code inexpensive to run, for the type of material non-linearities often encounter in the analysis of nuclear fuel. The code, however, is general, and can be used for other components of the reactor, or even for non-nuclear systems. The predictions of FEAT have been compared against several analytical solutions. The agreement is usually better than 5%. Thermocouple measurements show that the FEAT predictions are consistent with measured changes in temperatures in simulated pressure tubes. FEAT was also found to predict well, the axial variations in temperatures in the end-pellets(UO 2 ) of two fuel elements irradiated

Tests of a 3D Self Magnetic Field Solver in the Finite Element Gun Code MICHELLE

CERN Document Server

Nelson, Eric M

2005-01-01

We have recently implemented a prototype 3d self magnetic field solver in the finite-element gun code MICHELLE. The new solver computes the magnetic vector potential on unstructured grids. The solver employs edge basis functions in the curl-curl formulation of the finite-element method. A novel current accumulation algorithm takes advantage of the unstructured grid particle tracker to produce a compatible source vector, for which the singular matrix equation is easily solved by the conjugate gradient method. We will present some test cases demonstrating the capabilities of the prototype 3d self magnetic field solver. One test case is self magnetic field in a square drift tube. Another is a relativistic axisymmetric beam freely expanding in a round pipe.

Guided waves dispersion equations for orthotropic multilayered pipes solved using standard finite elements code.

Science.gov (United States)

Predoi, Mihai Valentin

2014-09-01

The dispersion curves for hollow multilayered cylinders are prerequisites in any practical guided waves application on such structures. The equations for homogeneous isotropic materials have been established more than 120 years ago. The difficulties in finding numerical solutions to analytic expressions remain considerable, especially if the materials are orthotropic visco-elastic as in the composites used for pipes in the last decades. Among other numerical techniques, the semi-analytical finite elements method has proven its capability of solving this problem. Two possibilities exist to model a finite elements eigenvalue problem: a two-dimensional cross-section model of the pipe or a radial segment model, intersecting the layers between the inner and the outer radius of the pipe. The last possibility is here adopted and distinct differential problems are deduced for longitudinal L(0,n), torsional T(0,n) and flexural F(m,n) modes. Eigenvalue problems are deduced for the three modes classes, offering explicit forms of each coefficient for the matrices used in an available general purpose finite elements code. Comparisons with existing solutions for pipes filled with non-linear viscoelastic fluid or visco-elastic coatings as well as for a fully orthotropic hollow cylinder are all proving the reliability and ease of use of this method. Copyright © 2014 Elsevier B.V. All rights reserved.

Recent uses of the finite element method in design/analysis of CANDU fuel

International Nuclear Information System (INIS)

Tayal, M.; Lim, D.

1985-06-01

Finite element codes FEAST and ELESTRES have been used to show: that initial pellet density can have a significant effect on the probability of fuel defect near end cap welds; that sheath stresses/strains are highly multiaxial near circumferential ridges; and that the multiaxiality affects sheath integrity significantly. The finite element thermal code FEAT was used to redesign bearing pads to obtain lower temperture; this eliminated crevice corrosion. FEAT was also used to assess the influences of braze voids and of end flux peaking. These analyses involved complex geometries. By using finite elements, we could obtain accurate assessments economically and rapidly. Finite element codes are also being developed for bowing, diffusion, flow patterns, and stress corrosion cracking

Implementation of second moment closure turbulence model for incompressible flows in the industrial finite element code N3S

International Nuclear Information System (INIS)

Pot, G.; Laurence, D.; Rharif, N.E.; Leal de Sousa, L.; Compe, C.

1995-12-01

This paper deals with the introduction of a second moment closure turbulence model (Reynolds Stress Model) in an industrial finite element code, N3S, developed at Electricite de France.The numerical implementation of the model in N3S will be detailed in 2D and 3D. Some details are given concerning finite element computations and solvers. Then, some results will be given, including a comparison between standard k-ε model, R.S.M. model and experimental data for some test case. (authors). 22 refs., 3 figs

Mixed finite element simulations in two-dimensional groundwater flow problems

International Nuclear Information System (INIS)

Kimura, Hideo

1989-01-01

A computer code of groundwater flow in two-dimensional porous media based on the mixed finite element method was developed for accurate approximations of Darcy velocities in safety evaluation of radioactive waste disposal. The mixed finite element procedure solves for both the Darcy velocities and pressure heads simultaneously in the Darcy equation and continuity equation. Numerical results of a single well pumping at a constant rate in a uniform flow field showed that the mixed finite element method gives more accurate Darcy velocities nearly 50 % on average error than standard finite element method. (author)

Implementation of thermo-viscoplastic constitutive equations into the finite element code ABAQUS

International Nuclear Information System (INIS)

Youn, Sam Son; Lee, Soon Bok; Kim, Jong Bum; Lee, Hyeong Yeon; Yoo, Bong

1998-01-01

Sophisticated viscoplatic constitutive laws describing material behavior at high temperature have been implemented in the general-purpose finite element code ABAQUS to predict the viscoplastic response of structures to cyclic loading. Because of the complexity of viscoplastic constitutive equation, the general implementation methods are developed. The solution of the non-linear system of algebraic equations arising from time discretization is determined using line-search and back-tracking in combination with Newton method. The time integration method of the constitutive equations is based on semi-implicit method with efficient time step control. For numerical examples, the viscoplastic model proposed by Chaboche is implemented and several applications are illustrated

A General Finite Element Scheme for Limit State Analysis and Optimization

DEFF Research Database (Denmark)

Damkilde, Lars

1999-01-01

Limit State analysis which is based on a perfect material behaviour is used in many different applications primarily within Structural Engineering and Geotechnics. The calculation methods have not reached the same level of automation such as Finite Element Analysis for elastic structures....... The computer based systems are more ad hoc based and are typically not well-integrated with pre- and postprocessors well-known from commercial Finite Element codes.A finite element based formulation of limit state analysis is presented which allows an easy integration with standard Finite Element codes...... for elastic analysis. In this way the user is able to perform a limit state analysis on the same model used for elastic analysis only adding data for the yield surface.The method is based on the lower-bound theorem and uses stress-based elements with a linearized yield surface. The mathematical problem...

Validation of the 3D finite element transport theory code EVENT for shielding applications

International Nuclear Information System (INIS)

Warner, Paul; Oliveira, R.E. de

2000-01-01

This paper is concerned with the validation of the 3D deterministic neutral-particle transport theory code EVENT for shielding applications. The code is based on the finite element-spherical harmonics (FE-P N ) method which has been extensively developed over the last decade. A general multi-group, anisotropic scattering formalism enables the code to address realistic steady state and time dependent, multi-dimensional coupled neutron/gamma radiation transport problems involving high scattering and deep penetration alike. The powerful geometrical flexibility and competitive computational effort makes the code an attractive tool for shielding applications. In recognition of this, EVENT is currently in the process of being adopted by the UK nuclear industry. The theory behind EVENT is described and its numerical implementation is outlined. Numerical results obtained by the code are compared with predictions of the Monte Carlo code MCBEND and also with the results from benchmark shielding experiments. In particular, results are presented for the ASPIS experimental configuration for both neutron and gamma ray calculations using the BUGLE 96 nuclear data library. (author)

Full scale seismic simulation of a nuclear reactor with parallel finite element analysis code for assembled structure

International Nuclear Information System (INIS)

Yamada, Tomonori

2010-01-01

The safety requirement of nuclear power plant attracts much attention nowadays. With the growing computing power, numerical simulation is one of key technologies to meet this safety requirement. Center for Computational Science and e-Systems of Japan Atomic Energy Agency has been developing a finite element analysis code for assembled structure to accurately evaluate the structural integrity of nuclear power plant in its entirety under seismic events. Because nuclear power plant is very huge assembled structure with tens of millions of mechanical components, the finite element model of each component is assembled into one structure and non-conforming meshes of mechanical components are bonded together inside the code. The main technique to bond these mechanical components is triple sparse matrix multiplication with multiple point constrains and global stiffness matrix. In our code, this procedure is conducted in a component by component manner, so that the working memory size and computing time for this multiplication are available on the current computing environment. As an illustrative example, seismic simulation of a real nuclear reactor of High Temperature engineering Test Reactor, which is located at the O-arai research and development center of JAEA, with 80 major mechanical components was conducted. Consequently, our code successfully simulated detailed elasto-plastic deformation of nuclear reactor and its computational performance was investigated. (author)

Numerical solution of multi group-Two dimensional- Adjoint equation with finite element method

International Nuclear Information System (INIS)

Poursalehi, N.; Khalafi, H.; Shahriari, M.; Minoochehr

2008-01-01

Adjoint equation is used for perturbation theory in nuclear reactor design. For numerical solution of adjoint equation, usually two methods are applied. These are Finite Element and Finite Difference procedures. Usually Finite Element Procedure is chosen for solving of adjoint equation, because it is more use able in variety of geometries. In this article, Galerkin Finite Element method is discussed. This method is applied for numerical solving multi group, multi region and two dimensional (X, Y) adjoint equation. Typical reactor geometry is partitioned with triangular meshes and boundary condition for adjoint flux is considered zero. Finally, for a case of defined parameters, Finite Element Code was applied and results were compared with Citation Code

Non-Linear Three Dimensional Finite Elements for Composite Concrete Structures

Directory of Open Access Journals (Sweden)

O. Kohnehpooshi

Full Text Available Abstract The current investigation focused on the development of effective and suitable modelling of reinforced concrete component with and without strengthening. The modelling includes physical and constitutive models. New interface elements have been developed, while modified constitutive law have been applied and new computational algorithm is utilised. The new elements are the Truss-link element to model the interaction between concrete and reinforcement bars, the interface element between two plate bending elements and the interface element to represent the interfacial behaviour between FRP, steel plates and concrete. Nonlinear finite-element (FE codes were developed with pre-processing. The programme was written using FORTRAN language. The accuracy and efficiency of the finite element programme were achieved by analyzing several examples from the literature. The application of the 3D FE code was further enhanced by carrying out the numerical analysis of the three dimensional finite element analysis of FRP strengthened RC beams, as well as the 3D non-linear finite element analysis of girder bridge. Acceptable distributions of slip, deflection, stresses in the concrete and FRP plate have also been found. These results show that the new elements are effective and appropriate to be used for structural component modelling.

Magnetic materials and 3D finite element modeling

CERN Document Server

Bastos, Joao Pedro A

2014-01-01

Magnetic Materials and 3D Finite Element Modeling explores material characterization and finite element modeling (FEM) applications. This book relates to electromagnetic analysis based on Maxwell’s equations and application of the finite element (FE) method to low frequency devices. A great source for senior undergraduate and graduate students in electromagnetics, it also supports industry professionals working in magnetics, electromagnetics, ferromagnetic materials science and electrical engineering. The authors present current concepts on ferromagnetic material characterizations and losses. They provide introductory material; highlight basic electromagnetics, present experimental and numerical modeling related to losses and focus on FEM applied to 3D applications. They also explain various formulations, and discuss numerical codes.

Elements of algebraic coding systems

CERN Document Server

Cardoso da Rocha, Jr, Valdemar

2014-01-01

Elements of Algebraic Coding Systems is an introductory text to algebraic coding theory. In the first chapter, you'll gain inside knowledge of coding fundamentals, which is essential for a deeper understanding of state-of-the-art coding systems. This book is a quick reference for those who are unfamiliar with this topic, as well as for use with specific applications such as cryptography and communication. Linear error-correcting block codes through elementary principles span eleven chapters of the text. Cyclic codes, some finite field algebra, Goppa codes, algebraic decoding algorithms, and applications in public-key cryptography and secret-key cryptography are discussed, including problems and solutions at the end of each chapter. Three appendices cover the Gilbert bound and some related derivations, a derivation of the Mac- Williams' identities based on the probability of undetected error, and two important tools for algebraic decoding-namely, the finite field Fourier transform and the Euclidean algorithm f...

Software requirements, design, and verification and validation for the FEHM application - a finite-element heat- and mass-transfer code

International Nuclear Information System (INIS)

Dash, Z.V.; Robinson, B.A.; Zyvoloski, G.A.

1997-07-01

The requirements, design, and verification and validation of the software used in the FEHM application, a finite-element heat- and mass-transfer computer code that can simulate nonisothermal multiphase multicomponent flow in porous media, are described. The test of the DOE Code Comparison Project, Problem Five, Case A, which verifies that FEHM has correctly implemented heat and mass transfer and phase partitioning, is also covered

Algorithms and data structures for massively parallel generic adaptive finite element codes

KAUST Repository

Bangerth, Wolfgang

2011-12-01

Today\\'s largest supercomputers have 100,000s of processor cores and offer the potential to solve partial differential equations discretized by billions of unknowns. However, the complexity of scaling to such large machines and problem sizes has so far prevented the emergence of generic software libraries that support such computations, although these would lower the threshold of entry and enable many more applications to benefit from large-scale computing. We are concerned with providing this functionality for mesh-adaptive finite element computations. We assume the existence of an "oracle" that implements the generation and modification of an adaptive mesh distributed across many processors, and that responds to queries about its structure. Based on querying the oracle, we develop scalable algorithms and data structures for generic finite element methods. Specifically, we consider the parallel distribution of mesh data, global enumeration of degrees of freedom, constraints, and postprocessing. Our algorithms remove the bottlenecks that typically limit large-scale adaptive finite element analyses. We demonstrate scalability of complete finite element workflows on up to 16,384 processors. An implementation of the proposed algorithms, based on the open source software p4est as mesh oracle, is provided under an open source license through the widely used deal.II finite element software library. © 2011 ACM 0098-3500/2011/12-ART10 $10.00.

European column buckling curves and finite element modelling including high strength steels

DEFF Research Database (Denmark)

Jönsson, Jeppe; Stan, Tudor-Cristian

2017-01-01

Eurocode allows for finite element modelling of plated steel structures, however the information in the code on how to perform the analysis or what assumptions to make is quite sparse. The present paper investigates the deterministic modelling of flexural column buckling using plane shell elements...... imperfections may be very conservative if considered by finite element analysis as described in the current Eurocode code. A suggestion is given for a slightly modified imperfection formula within the Ayrton-Perry formulation leading to adequate inclusion of modern high grade steels within the original four...... bucking curves. It is also suggested that finite element or frame analysis may be performed with equivalent column bow imperfections extracted directly from the Ayrton-Perry formulation....

Validation of High Displacement Piezoelectric Actuator Finite Element Models

Science.gov (United States)

Taleghani, B. K.

2000-01-01

The paper presents the results obtained by using NASTRAN(Registered Trademark) and ANSYS(Regitered Trademark) finite element codes to predict doming of the THUNDER piezoelectric actuators during the manufacturing process and subsequent straining due to an applied input voltage. To effectively use such devices in engineering applications, modeling and characterization are essential. Length, width, dome height, and thickness are important parameters for users of such devices. Therefore, finite element models were used to assess the effects of these parameters. NASTRAN(Registered Trademark) and ANSYS(Registered Trademark) used different methods for modeling piezoelectric effects. In NASTRAN(Registered Trademark), a thermal analogy was used to represent voltage at nodes as equivalent temperatures, while ANSYS(Registered Trademark) processed the voltage directly using piezoelectric finite elements. The results of finite element models were validated by using the experimental results.

Development of a computer code 'CRACK' for elastic and elastoplastic fracture mechanics analysis of 2-D structures by finite element technique

International Nuclear Information System (INIS)

Dutta, B.K.; Kakodkar, A.; Maiti, S.K.

1986-01-01

The fracture mechanics analysis of nuclear components is required to ensure prevention of sudden failure due to dynamic loadings. The linear elastic analysis near to a crack tip shows presence of stress singularity at the crack tip. The simulation of this singularity in numerical methods enhance covergence capability. In finite element technique this can be achieved by placing mid nodes of 8 noded or 6 noded isoparametric elements, at one fourth ditance from crack tip. Present report details this characteristic of finite element, implementation of this element in a code 'CRACK', implementation of J-integral to compute stress intensity factor and solution of number of cases for elastic and elastoplastic fracture mechanics analysis. 6 refs., 6 figures. (author)

Rn3D: A finite element code for simulating gas flow and radon transport in variably saturated, nonisothermal porous media

International Nuclear Information System (INIS)

Holford, D.J.

1994-01-01

This document is a user's manual for the Rn3D finite element code. Rn3D was developed to simulate gas flow and radon transport in variably saturated, nonisothermal porous media. The Rn3D model is applicable to a wide range of problems involving radon transport in soil because it can simulate either steady-state or transient flow and transport in one-, two- or three-dimensions (including radially symmetric two-dimensional problems). The porous materials may be heterogeneous and anisotropic. This manual describes all pertinent mathematics related to the governing, boundary, and constitutive equations of the model, as well as the development of the finite element equations used in the code. Instructions are given for constructing Rn3D input files and executing the code, as well as a description of all output files generated by the code. Five verification problems are given that test various aspects of code operation, complete with example input files, FORTRAN programs for the respective analytical solutions, and plots of model results. An example simulation is presented to illustrate the type of problem Rn3D is designed to solve. Finally, instructions are given on how to convert Rn3D to simulate systems other than radon, air, and water

Non-linear heat transfer computer code by finite element method

International Nuclear Information System (INIS)

Nagato, Kotaro; Takikawa, Noboru

1977-01-01

The computer code THETA-2D for the calculation of temperature distribution by the two-dimensional finite element method was made for the analysis of heat transfer in a high temperature structure. Numerical experiment was performed for the numerical integration of the differential equation of heat conduction. The Runge-Kutta method of the numerical experiment produced an unstable solution. A stable solution was obtained by the β method with the β value of 0.35. In high temperature structures, the radiative heat transfer can not be neglected. To introduce a term of the radiative heat transfer, a functional neglecting the radiative heat transfer was derived at first. Then, the radiative term was added after the discretion by variation method. Five model calculations were carried out by the computer code. Calculation of steady heat conduction was performed. When estimated initial temperature is 1,000 degree C, reasonable heat blance was obtained. In case of steady-unsteady temperature calculation, the time integral by THETA-2D turned out to be under-estimation for enthalpy change. With a one-dimensional model, the temperature distribution in a structure, in which heat conductivity is dependent on temperature, was calculated. Calculation with a model which has a void inside was performed. Finally, model calculation for a complex system was carried out. (Kato, T.)

ELLIPT2D: A Flexible Finite Element Code Written Python

International Nuclear Information System (INIS)

Pletzer, A.; Mollis, J.C.

2001-01-01

The use of the Python scripting language for scientific applications and in particular to solve partial differential equations is explored. It is shown that Python's rich data structure and object-oriented features can be exploited to write programs that are not only significantly more concise than their counter parts written in Fortran, C or C++, but are also numerically efficient. To illustrate this, a two-dimensional finite element code (ELLIPT2D) has been written. ELLIPT2D provides a flexible and easy-to-use framework for solving a large class of second-order elliptic problems. The program allows for structured or unstructured meshes. All functions defining the elliptic operator are user supplied and so are the boundary conditions, which can be of Dirichlet, Neumann or Robbins type. ELLIPT2D makes extensive use of dictionaries (hash tables) as a way to represent sparse matrices.Other key features of the Python language that have been widely used include: operator over loading, error handling, array slicing, and the Tkinter module for building graphical use interfaces. As an example of the utility of ELLIPT2D, a nonlinear solution of the Grad-Shafranov equation is computed using a Newton iterative scheme. A second application focuses on a solution of the toroidal Laplace equation coupled to a magnetohydrodynamic stability code, a problem arising in the context of magnetic fusion research

Synthesis of hydrocode and finite element technology for large deformation Lagrangian computation

International Nuclear Information System (INIS)

Goudreau, G.L.; Hallquist, J.O.

1979-08-01

Large deformation engineering analysis at Lawrence Livermore Laboratory has benefited from a synthesis of computational technology from the finite difference hydrocodes of the scientific weapons community and the structural finite element methodology of engineering. Two- and three-dimensional explicit and implicit Lagrangian continuum codes have been developed exploiting the strengths of each. The explicit methodology primarily exploits the primitive constant stress (or one point integration) brick element. Similarity and differences with the integral finite difference method are discussed. Choice of stress and finite strain measures, and selection of hour glass viscosity are also considered. The implicit codes also employ a Cauchy formulation, with Newton iteration and a symmetric tangent matrix. A library of finite strain material routines includes hypoelastic/plastic, hyperelastic, viscoelastic, as well as hydrodynamic behavior. Arbitrary finite element topology and a general slide-line treatment significantly extends Lagrangian hydrocode application. Computational experience spans weapons and non-weapons applications

Development of dynamic explicit crystallographic homogenization finite element analysis code to assess sheet metal formability

International Nuclear Information System (INIS)

Nakamura, Yasunori; Tam, Nguyen Ngoc; Ohata, Tomiso; Morita, Kiminori; Nakamachi, Eiji

2004-01-01

The crystallographic texture evolution induced by plastic deformation in the sheet metal forming process has a great influence on its formability. In the present study, a dynamic explicit finite element (FE) analysis code is newly developed by introducing a crystallographic homogenization method to estimate the polycrystalline sheet metal formability, such as the extreme thinning and 'earing'. This code can predict the plastic deformation induced texture evolution at the micro scale and the plastic anisotropy at the macro scale, simultaneously. This multi-scale analysis can couple the microscopic crystal plasticity inhomogeneous deformation with the macroscopic continuum deformation. In this homogenization process, the stress at the macro scale is defined by the volume average of those of the corresponding microscopic crystal aggregations in satisfying the equation of motion and compatibility condition in the micro scale 'unit cell', where the periodicity of deformation is satisfied. This homogenization algorithm is implemented in the conventional dynamic explicit finite element code by employing the updated Lagrangian formulation and the rate type elastic/viscoplastic constitutive equation.At first, it has been confirmed through a texture evolution analyses in cases of typical deformation modes that Taylor's 'constant strain homogenization algorithm' yields extreme concentration toward the preferred crystal orientations compared with our homogenization one. Second, we study the plastic anisotropy effects on 'earing' in the hemispherical cup deep drawing process of pure ferrite phase sheet metal. By the comparison of analytical results with those of Taylor's assumption, conclusions are drawn that the present newly developed dynamic explicit crystallographic homogenization FEM shows more reasonable prediction of plastic deformation induced texture evolution and plastic anisotropy at the macro scale

Automation of finite element methods

CERN Document Server

Korelc, Jože

2016-01-01

New finite elements are needed as well in research as in industry environments for the development of virtual prediction techniques. The design and implementation of novel finite elements for specific purposes is a tedious and time consuming task, especially for nonlinear formulations. The automation of this process can help to speed up this process considerably since the generation of the final computer code can be accelerated by order of several magnitudes. This book provides the reader with the required knowledge needed to employ modern automatic tools like AceGen within solid mechanics in a successful way. It covers the range from the theoretical background, algorithmic treatments to many different applications. The book is written for advanced students in the engineering field and for researchers in educational and industrial environments.

The computer code EURDYN - 1 M (release 1) for transient dynamic fluid-structure interaction. Pt.1: governing equations and finite element modelling

International Nuclear Information System (INIS)

Donea, J.; Fasoli-Stella, P.; Giuliani, S.; Halleux, J.P.; Jones, A.V.

1980-01-01

This report describes the governing equations and the finite element modelling used in the computer code EURDYN - 1 M. The code is a non-linear transient dynamic program for the analysis of coupled fluid-structure systems; It is designed for safety studies on LMFBR components (primary containment and fuel subassemblies)

Introduction of polycrystal constitutive laws in a finite element code with applications to zirconium forming

International Nuclear Information System (INIS)

Maudlin, P.J.; Tome, C.N.; Kaschner, G.C.; Gray, G.T. III

1998-01-01

In this work the authors simulate the compressive deformation of heavily textured zirconium sheet using a finite element code with the constitutive response given by a polycrystal self-consistent model. They show that the strong anisotropy of the response can be explained in terms of the texture and the relative activity of prismatic (easy) and pyramidal (hard) slip modes. The simulations capture the yield anisotropy observed for so-called through-thickness and in-plane compression tests in terMs of the loading curves and final specimen geometries

Solution of 2D and 3D hexagonal geometry benchmark problems by using the finite element diffusion code DIFGEN

International Nuclear Information System (INIS)

Gado, J.

1986-02-01

The four group, 2D and 3D hexagonal geometry HTGR benchmark problems and a 2D hexagonal geometry PWR (WWER) benchmark problem have been solved by using the finite element diffusion code DIFGEN. The hexagons (or hexagonal prisms) were subdivided into first order or second order triangles or quadrilaterals (or triangular or quadrilateral prisms). In the 2D HTGR case of the number of the inserted absorber rods was also varied (7, 6, 0 or 37 rods). The calculational results are in a good agreement with the results of other calculations. The larger systematic series of DIFGEN calculations have given a quantitative picture on the convergence properties of various finite element modellings of hexagonal grids in DIFGEN. (orig.)

3D unstructured mesh discontinuous finite element hydro

International Nuclear Information System (INIS)

Prasad, M.K.; Kershaw, D.S.; Shaw, M.J.

1995-01-01

The authors present detailed features of the ICF3D hydrodynamics code used for inertial fusion simulations. This code is intended to be a state-of-the-art upgrade of the well-known fluid code, LASNEX. ICF3D employs discontinuous finite elements on a discrete unstructured mesh consisting of a variety of 3D polyhedra including tetrahedra, prisms, and hexahedra. The authors discussed details of how the ROE-averaged second-order convection was applied on the discrete elements, and how the C++ coding interface has helped to simplify implementing the many physics and numerics modules within the code package. The author emphasized the virtues of object-oriented design in large scale projects such as ICF3D

Some aspects of the quality assurance of personnel carrying out finite element analysis

International Nuclear Information System (INIS)

Dickenson, P.W.

1990-01-01

In this paper, the need to assess the competence of personnel carrying out finite element analysis is emphasised. In carrying out its regulatory role on behalf of the Health and Safety Executive, the Nuclear Installations Inspectorate (NII) must be satisfied that appropriate standards are developed and maintained by the licensee. Since finite element methods have an important bearing on the acceptance of a safety case, it follows that relevant codes should be adequately validated and the personnel applying the code should be competent. Attention is drawn to the work of the Quality Assurance Working Group of the National Agency for Finite Element Methods and Standards (NAFEMS) who are active in this area. The paper also considers the methods that are available to assess the competence of personnel engaged in finite element methods. (author)

Aspects of Finite Element Simulation of Axi-Symmetric Hydromechanical Deep Drawing

DEFF Research Database (Denmark)

Jensen, Morten Rikard; Olovsson, Lars; Danckert, Joachim

1999-01-01

A new approach for the Finite Element modelling of the hydromechanical deep drawing process is evaluated. In the model a Finite Difference approximation of Reynold’s equation is solved for the fluid flow between the blank and the draw die in the flange region. The approach is implemented...... as a contact algorithm in an explicit Finite Element code, Exhale2D. The developed model is verified against experiments and good agreement is obtained. It is concluded that the developed model is a promising approach for simulating the hydromechanical deep drawing process using the Finite Element Method....

Validation of finite element code DELFIN by means of the zero power experiences at the nuclear power plant of Atucha I

International Nuclear Information System (INIS)

Grant, C.R.

1996-01-01

Code DELFIN, developed in CNEA, treats the spatial discretization using heterogeneous finite elements, allowing a correct treatment of the continuity of fluxes and currents among elements and a more realistic representation of the hexagonal lattice of the reactor. It can be used for fuel management calculation, Xenon oscillation and spatial kinetics. Using the HUEMUL code for cell calculation (which uses a generalized two dimensional collision probability theory and has the WIMS library incorporated in a data base), the zero power experiences performed in 1974 were calculated. (author). 8 refs., 9 figs., 3 tabs

Variational approach to probabilistic finite elements

Science.gov (United States)

Belytschko, T.; Liu, W. K.; Mani, A.; Besterfield, G.

1991-08-01

Probabilistic finite element methods (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Material model for non-linear finite element analyses of large concrete structures

NARCIS (Netherlands)

Engen, Morten; Hendriks, M.A.N.; Øverli, Jan Arve; Åldstedt, Erik; Beushausen, H.

2016-01-01

A fully triaxial material model for concrete was implemented in a commercial finite element code. The only required input parameter was the cylinder compressive strength. The material model was suitable for non-linear finite element analyses of large concrete structures. The importance of including

ABAQUS2MATLAB: A Novel Tool for Finite Element Post-Processing

DEFF Research Database (Denmark)

Martínez Pañeda, Emilio; Papazafeiropoulos, George; Muniz-Calvente, Miguel

2017-01-01

A novel piece of software is presented to connect Abaqus, a sophisticated finite element package, with Matlab, the most comprehensive program for mathematical analysis. This interface between these well-known codes not only benefits from the image processing and the integrated graph-plotting feat......A novel piece of software is presented to connect Abaqus, a sophisticated finite element package, with Matlab, the most comprehensive program for mathematical analysis. This interface between these well-known codes not only benefits from the image processing and the integrated graph...... to demonstrate its capabilities. The source code, detailed documentation and a large number of tutorials can be freely downloaded from www.abaqus2matlab.com....

User's Manual for the FEHM Application-A Finite-Element Heat- and Mass-Transfer Code

Energy Technology Data Exchange (ETDEWEB)

George A. Zyvoloski; Bruce A. Robinson; Zora V. Dash; Lynn L. Trease

1997-07-07

This document is a manual for the use of the FEHM application, a finite-element heat- and mass-transfer computer code that can simulate nonisothermal multiphase multicomponent flow in porous media. The use of this code is applicable to natural-state studies of geothermal systems and groundwater flow. A primary use of the FEHM application will be to assist in the understanding of flow fields and mass transport in the saturated and unsaturated zones below the proposed Yucca Mountain nuclear waste repository in Nevada. The equations of heat and mass transfer for multiphase flow in porous and permeable media are solved in the FEHM application by using the finite-element method. The permeability and porosity of the medium are allowed to depend on pressure and temperature. The code also has provisions for movable air and water phases and noncoupled tracers; that is, tracer solutions that do not affect the heat- and mass-transfer solutions. The tracers can be passive or reactive. The code can simulate two-dimensional, two-dimensional radial, or three-dimensional geometries. In fact, FEHM is capable of describing flow that is dominated in many areas by fracture and fault flow, including the inherently three-dimensional flow that results from permeation to and from faults and fractures. The code can handle coupled heat and mass-transfer effects, such as boiling, dryout, and condensation that can occur in the near-field region surrounding the potential repository and the natural convection that occurs through Yucca Mountain due to seasonal temperature changes. The code is also capable of incorporating the various adsorption mechanisms, ranging from simple linear relations to nonlinear isotherms, needed to describe the very complex transport processes at Yucca Mountain. This report outlines the uses and capabilities of the FEHM application, initialization of code variables, restart procedures, and error processing. The report describes all the data files, the input data

Coupling of smooth particle hydrodynamics with the finite element method

International Nuclear Information System (INIS)

Attaway, S.W.; Heinstein, M.W.; Swegle, J.W.

1994-01-01

A gridless technique called smooth particle hydrodynamics (SPH) has been coupled with the transient dynamics finite element code ppercase[pronto]. In this paper, a new weighted residual derivation for the SPH method will be presented, and the methods used to embed SPH within ppercase[pronto] will be outlined. Example SPH ppercase[pronto] calculations will also be presented. One major difficulty associated with the Lagrangian finite element method is modeling materials with no shear strength; for example, gases, fluids and explosive biproducts. Typically, these materials can be modeled for only a short time with a Lagrangian finite element code. Large distortions cause tangling of the mesh, which will eventually lead to numerical difficulties, such as negative element area or ''bow tie'' elements. Remeshing will allow the problem to continue for a short while, but the large distortions can prevent a complete analysis. SPH is a gridless Lagrangian technique. Requiring no mesh, SPH has the potential to model material fracture, large shear flows and penetration. SPH computes the strain rate and the stress divergence based on the nearest neighbors of a particle, which are determined using an efficient particle-sorting technique. Embedding the SPH method within ppercase[pronto] allows part of the problem to be modeled with quadrilateral finite elements, while other parts are modeled with the gridless SPH method. SPH elements are coupled to the quadrilateral elements through a contact-like algorithm. ((orig.))

Advanced finite element simulation with MSC Marc application of user subroutines

CERN Document Server

Javanbakht, Zia

2017-01-01

This book offers an in-depth insight into the general-purpose finite element program MSC Marc, which is distributed by MSC Software Corporation. It is a specialized program for nonlinear problems (implicit solver) which is common in academia and industry. The primary goal of this book is to provide a comprehensive introduction to a special feature of this software: the user can write user-subroutines in the programming language Fortran, which is the language of all classical finite element packages. This subroutine feature allows the user to replace certain modules of the core code and to implement new features such as constitutive laws or new elements. Thus, the functionality of commercial codes (‘black box’) can easily be extended by linking user written code to the main core of the program. This feature allows to take advantage of a commercial software package with the flexibility of a ‘semi-open’ code. .

Advances in dynamic relaxation techniques for nonlinear finite element analysis

International Nuclear Information System (INIS)

Sauve, R.G.; Metzger, D.R.

1995-01-01

Traditionally, the finite element technique has been applied to static and steady-state problems using implicit methods. When nonlinearities exist, equilibrium iterations must be performed using Newton-Raphson or quasi-Newton techniques at each load level. In the presence of complex geometry, nonlinear material behavior, and large relative sliding of material interfaces, solutions using implicit methods often become intractable. A dynamic relaxation algorithm is developed for inclusion in finite element codes. The explicit nature of the method avoids large computer memory requirements and makes possible the solution of large-scale problems. The method described approaches the steady-state solution with no overshoot, a problem which has plagued researchers in the past. The method is included in a general nonlinear finite element code. A description of the method along with a number of new applications involving geometric and material nonlinearities are presented. They include: (1) nonlinear geometric cantilever plate; (2) moment-loaded nonlinear beam; and (3) creep of nuclear fuel channel assemblies

TRIP: a finite element computer program for the solution of convection heat transfer problems

International Nuclear Information System (INIS)

Slagter, W.; Roodbergen, H.A.

1976-01-01

The theory and use of the finite element code TRIP are described. The code calculates temperature distributions in three-dimensional continua subjected to convection heat transfer. A variational principle for transport phenomena is applied to solve the convection heat transfer problem with temperature and heat flux boundary conditions. The finite element discretization technique is used to reduce the continuous spatial solution into a finite number of unknowns. The method is developed in detail to determine temperature distributions in coolant passages of fuel rod bundles which are idealized by hexahedral elements. The development of the TRIP code is discussed and the listing of the program is given in FORTRAN IV. An example is given to illustrate the validity and practicality of the method

Finite elements methods in mechanics

CERN Document Server

Eslami, M Reza

2014-01-01

This book covers all basic areas of mechanical engineering, such as fluid mechanics, heat conduction, beams, and elasticity with detailed derivations for the mass, stiffness, and force matrices. It is especially designed to give physical feeling to the reader for finite element approximation by the introduction of finite elements to the elevation of elastic membrane. A detailed treatment of computer methods with numerical examples are provided. In the fluid mechanics chapter, the conventional and vorticity transport formulations for viscous incompressible fluid flow with discussion on the method of solution are presented. The variational and Galerkin formulations of the heat conduction, beams, and elasticity problems are also discussed in detail. Three computer codes are provided to solve the elastic membrane problem. One of them solves the Poisson’s equation. The second computer program handles the two dimensional elasticity problems, and the third one presents the three dimensional transient heat conducti...

Abaqus2Matlab: A suitable tool for finite element post-processing

DEFF Research Database (Denmark)

Papazafeiropoulos, George; Muñiz-Calvente, Miguel; Martínez Pañeda, Emilio

2017-01-01

A suitable piece of software is presented to connect Abaqus, a sophisticated finite element package, with Matlab, the most comprehensive program for mathematical analysis. This interface between these well- known codes not only benefits from the image processing and the integrated graph-plotting ......A suitable piece of software is presented to connect Abaqus, a sophisticated finite element package, with Matlab, the most comprehensive program for mathematical analysis. This interface between these well- known codes not only benefits from the image processing and the integrated graph...... crack propagation in structural materials by means of a cohesive zone approach. The source code, detailed documentation and a large number of tutorials can be freely downloaded from www.abaqus2matlab.com ....

NCEL: two dimensional finite element code for steady-state temperature distribution in seven rod-bundle

International Nuclear Information System (INIS)

Hrehor, M.

1979-01-01

The paper deals with an application of the finite element method to the heat transfer study in seven-pin models of LMFBR fuel subassembly. The developed code NCEL solves two-dimensional steady state heat conduction equation in the whole subassembly model cross-section and enebles to perform the analysis of thermal behaviour in both normal and accidental operational conditions as eccentricity of the central rod or full or partial (porous) blockage of some part of the cross-flow area. The heat removal is simulated by heat sinks in coolant under conditions of subchannels slug flow approximation

Mesh-size errors in diffusion-theory calculations using finite-difference and finite-element methods

International Nuclear Information System (INIS)

Baker, A.R.

1982-07-01

A study has been performed of mesh-size errors in diffusion-theory calculations using finite-difference and finite-element methods. As the objective was to illuminate the issues, the study was performed for a 1D slab model of a reactor with one neutron-energy group for which analytical solutions were possible. A computer code SLAB was specially written to perform the finite-difference and finite-element calculations and also to obtain the analytical solutions. The standard finite-difference equations were obtained by starting with an expansion of the neutron current in powers of the mesh size, h, and keeping terms as far as h 2 . It was confirmed that these equations led to the well-known result that the criticality parameter varied with the square of the mesh size. An improved form of the finite-difference equations was obtained by continuing the expansion for the neutron current as far as the term in h 4 . In this case, the critical parameter varied as the fourth power of the mesh size. The finite-element solutions for 2 and 3 nodes per element revealed that the criticality parameter varied as the square and fourth power of the mesh size, respectively. Numerical results are presented for a bare reactive core of uniform composition with 2 zones of different uniform mesh and for a reactive core with an absorptive reflector. (author)

Parallel Finite Element Particle-In-Cell Code for Simulations of Space-charge Dominated Beam-Cavity Interactions

International Nuclear Information System (INIS)

Candel, A.; Kabel, A.; Ko, K.; Lee, L.; Li, Z.; Limborg, C.; Ng, C.; Prudencio, E.; Schussman, G.; Uplenchwar, R.

2007-01-01

Over the past years, SLAC's Advanced Computations Department (ACD) has developed the parallel finite element (FE) particle-in-cell code Pic3P (Pic2P) for simulations of beam-cavity interactions dominated by space-charge effects. As opposed to standard space-charge dominated beam transport codes, which are based on the electrostatic approximation, Pic3P (Pic2P) includes space-charge, retardation and boundary effects as it self-consistently solves the complete set of Maxwell-Lorentz equations using higher-order FE methods on conformal meshes. Use of efficient, large-scale parallel processing allows for the modeling of photoinjectors with unprecedented accuracy, aiding the design and operation of the next-generation of accelerator facilities. Applications to the Linac Coherent Light Source (LCLS) RF gun are presented

Application of the finite element method to the neutron transport equation

International Nuclear Information System (INIS)

Martin, W.R.

1976-01-01

This paper examines the theoretical and practical application of the finite element method to the neutron transport equation. It is shown that in principle the system of equations obtained by application of the finite element method can be solved with certain physical restrictions concerning the criticality of the medium. The convergence of this approximate solution to the exact solution with mesh refinement is examined, and a non-optical estimate of the convergence rate is obtained analytically. It is noted that the numerical results indicate a faster convergence rate and several approaches to obtain this result analytically are outlined. The practical application of the finite element method involved the development of a computer code capable of solving the neutron transport equation in 1-D plane geometry. Vacuum, reflecting, or specified incoming boundary conditions may be analyzed, and all are treated as natural boundary conditions. The time-dependent transport equation is also examined and it is shown that the application of the finite element method in conjunction with the Crank-Nicholson time discretization method results in a system of algebraic equations which is readily solved. Numerical results are given for several critical slab eigenvalue problems, including anisotropic scattering, and the results compare extremely well with benchmark results. It is seen that the finite element code is more efficient than a standard discrete ordinates code for certain problems. A problem with severe heterogeneities is considered and it is shown that the use of discontinuous spatial and angular elements results in a marked improvement in the results. Finally, time-dependent problems are examined and it is seen that the phenomenon of angular mode separation makes the numerical treatment of the transport equation in slab geometry a considerable challenge, with the result that the angular mesh has a dominant effect on obtaining acceptable solutions

Development of Multi-Scale Finite Element Analysis Codes for High Formability Sheet Metal Generation

International Nuclear Information System (INIS)

Nnakamachi, Eiji; Kuramae, Hiroyuki; Ngoc Tam, Nguyen; Nakamura, Yasunori; Sakamoto, Hidetoshi; Morimoto, Hideo

2007-01-01

In this study, the dynamic- and static-explicit multi-scale finite element (F.E.) codes are developed by employing the homogenization method, the crystalplasticity constitutive equation and SEM-EBSD measurement based polycrystal model. These can predict the crystal morphological change and the hardening evolution at the micro level, and the macroscopic plastic anisotropy evolution. These codes are applied to analyze the asymmetrical rolling process, which is introduced to control the crystal texture of the sheet metal for generating a high formability sheet metal. These codes can predict the yield surface and the sheet formability by analyzing the strain path dependent yield, the simple sheet forming process, such as the limit dome height test and the cylindrical deep drawing problems. It shows that the shear dominant rolling process, such as the asymmetric rolling, generates ''high formability'' textures and eventually the high formability sheet. The texture evolution and the high formability of the newly generated sheet metal experimentally were confirmed by the SEM-EBSD measurement and LDH test. It is concluded that these explicit type crystallographic homogenized multi-scale F.E. code could be a comprehensive tool to predict the plastic induced texture evolution, anisotropy and formability by the rolling process and the limit dome height test analyses

Basic Finite Element Method

International Nuclear Information System (INIS)

Lee, Byeong Hae

1992-02-01

This book gives descriptions of basic finite element method, which includes basic finite element method and data, black box, writing of data, definition of VECTOR, definition of matrix, matrix and multiplication of matrix, addition of matrix, and unit matrix, conception of hardness matrix like spring power and displacement, governed equation of an elastic body, finite element method, Fortran method and programming such as composition of computer, order of programming and data card and Fortran card, finite element program and application of nonelastic problem.

Maxwell's equations in axisymmetrical geometry: coupling H(curl) finite element in volume and H(div) finite element in surface. The numerical code FuMel

International Nuclear Information System (INIS)

Cambon, S.; Lacoste, P.

2011-01-01

We propose a finite element method to solve the axisymmetric scattering problem posed on a regular bounded domain. Here we shall show how to reduce the initial 3D problem into a truncated sum of 2D independent problems posed into a meridian plane of the object. Each of these problem results in the coupling of a partial differential equation into the interior domain and an integral equation on the surface simulating the free space. Then variational volume and boundary integral formulations of Maxwell's equation on regular surfaces are derived. We introduce some general finite element adapted to cylindrical coordinates and constructed from nodal and mixed finite element both for the interior (volume) and for the integral equation (surface). (authors)

Finite element simulation of ironing process under warm conditions

Directory of Open Access Journals (Sweden)

Swadesh Kumar Singh

2014-01-01

Full Text Available Metal forming is one of the most important steps in manufacturing of a large variety of products. Ironing in deep drawing is done by adjusting the clearance between the punch and the die and allow the material flow over the punch. In the present investigation effect of extent of ironing behavior on the characteristics of the product like thickness distribution with respect to temperature was studied. With the help of finite element simulation using explicit finite element code LS-DYNA the stress in the drawn cup were predicted in the drawn cup. To increase the accuracy in the simulation process, numbers of integration points were increased in the thickness direction and it was found that there is very close prediction of finite element results to that of experimental ones.

Development of a finite element code to solve thermo-hydro-mechanical coupling and simulate induced seismicity.

Science.gov (United States)

María Gómez Castro, Berta; De Simone, Silvia; Rossi, Riccardo; Larese De Tetto, Antonia; Carrera Ramírez, Jesús

2015-04-01

Coupled thermo-hydro-mechanical modeling is essential for CO2 storage because of (1) large amounts of CO2 will be injected, which will cause large pressure buildups and might compromise the mechanical stability of the caprock seal, (2) the most efficient technique to inject CO2 is the cold injection, which induces thermal stress changes in the reservoir and seal. These stress variations can cause mechanical failure in the caprock and can also trigger induced earthquakes. To properly assess these effects, numerical models that take into account the short and long-term thermo-hydro-mechanical coupling are an important tool. For this purpose, there is a growing need of codes that couple these processes efficiently and accurately. This work involves the development of an open-source, finite element code written in C ++ for correctly modeling the effects of thermo-hydro-mechanical coupling in the field of CO2 storage and in others fields related to these processes (geothermal energy systems, fracking, nuclear waste disposal, etc.), and capable to simulate induced seismicity. In order to be able to simulate earthquakes, a new lower dimensional interface element will be implemented in the code to represent preexisting fractures, where pressure continuity will be imposed across the fractures.

An object-oriented decomposition of the adaptive-hp finite element method

Energy Technology Data Exchange (ETDEWEB)

Wiley, J.C.

1994-12-13

Adaptive-hp methods are those which use a refinement control strategy driven by a local error estimate to locally modify the element size, h, and polynomial order, p. The result is an unstructured mesh in which each node may be associated with a different polynomial order and which generally require complex data structures to implement. Object-oriented design strategies and languages which support them, e.g., C++, help control the complexity of these methods. Here an overview of the major classes and class structure of an adaptive-hp finite element code is described. The essential finite element structure is described in terms of four areas of computation each with its own dynamic characteristics. Implications of converting the code for a distributed-memory parallel environment are also discussed.

FINEDAN - an explicit finite-element calculation code for two-dimensional analyses of fast dynamic transients in nuclear reactor technology

International Nuclear Information System (INIS)

Adamik, V.; Matejovic, P.

1989-01-01

The problems are discussed of nonstationary, nonlinear dynamics of the continuum. A survey is presented of calculation methods in the given area with emphasis on the area of impact problems. A description is presented of the explicit finite elements method and its application to two-dimensional Cartesian and cylindrical configurations. Using the method the explicit calculation code FINEDAN was written which was tested in a series of verification calculations for different configurations and different types of continuum. The main characteristics are presented of the code and of some, of its practical applications. Envisaged trends of the development of the code and its possible applications in the technology of nuclear reactors are given. (author). 9 figs., 4 tabs., 10 refs

Essentials of the finite element method for mechanical and structural engineers

CERN Document Server

Pavlou, Dimitrios G

2015-01-01

Fundamental coverage, analytic mathematics, and up-to-date software applications are hard to find in a single text on the finite element method (FEM). Dimitrios Pavlou's Essentials of the Finite Element Method: For Structural and Mechanical Engineers makes the search easier by providing a comprehensive but concise text for those new to FEM, or just in need of a refresher on the essentials. Essentials of the Finite Element Method explains the basics of FEM, then relates these basics to a number of practical engineering applications. Specific topics covered include linear spring elements, bar elements, trusses, beams and frames, heat transfer, and structural dynamics. Throughout the text, readers are shown step-by-step detailed analyses for finite element equations development. The text also demonstrates how FEM is programmed, with examples in MATLAB, CALFEM, and ANSYS allowing readers to learn how to develop their own computer code. Suitable for everyone from first-time BSc/MSc students to practicing mechanic...

Propel: A Discontinuous-Galerkin Finite Element Code for Solving the Reacting Navier-Stokes Equations

Science.gov (United States)

Johnson, Ryan; Kercher, Andrew; Schwer, Douglas; Corrigan, Andrew; Kailasanath, Kazhikathra

2017-11-01

This presentation focuses on the development of a Discontinuous Galerkin (DG) method for application to chemically reacting flows. The in-house code, called Propel, was developed by the Laboratory of Computational Physics and Fluid Dynamics at the Naval Research Laboratory. It was designed specifically for developing advanced multi-dimensional algorithms to run efficiently on new and innovative architectures such as GPUs. For these results, Propel solves for convection and diffusion simultaneously with detailed transport and thermodynamics. Chemistry is currently solved in a time-split approach using Strang-splitting with finite element DG time integration of chemical source terms. Results presented here show canonical unsteady reacting flow cases, such as co-flow and splitter plate, and we report performance for higher order DG on CPU and GPUs.

High-Fidelity RF Gun Simulations with the Parallel 3D Finite Element Particle-In-Cell Code Pic3P

Energy Technology Data Exchange (ETDEWEB)

Candel, A; Kabel, A.; Lee, L.; Li, Z.; Limborg, C.; Ng, C.; Schussman, G.; Ko, K.; /SLAC

2009-06-19

SLAC's Advanced Computations Department (ACD) has developed the first parallel Finite Element 3D Particle-In-Cell (PIC) code, Pic3P, for simulations of RF guns and other space-charge dominated beam-cavity interactions. Pic3P solves the complete set of Maxwell-Lorentz equations and thus includes space charge, retardation and wakefield effects from first principles. Pic3P uses higher-order Finite Elementmethods on unstructured conformal meshes. A novel scheme for causal adaptive refinement and dynamic load balancing enable unprecedented simulation accuracy, aiding the design and operation of the next generation of accelerator facilities. Application to the Linac Coherent Light Source (LCLS) RF gun is presented.

Finite element analysis of the neutron transport equation in spherical geometry

International Nuclear Information System (INIS)

Kim, Yong Ill; Kim, Jong Kyung; Suk, Soo Dong

1992-01-01

The Galerkin formulation of the finite element method is applied to the integral law of the first-order form of the one-group neutron transport equation in one-dimensional spherical geometry. Piecewise linear or quadratic Lagrange polynomials are utilized in the integral law for the angular flux to establish a set of linear algebraic equations. Numerical analyses are performed for the scalar flux distribution in a heterogeneous sphere as well as for the criticality problem in a uniform sphere. For the criticality problems in the uniform sphere, the results of the finite element method, with the use of continuous finite elements in space and angle, are compared with the exact solutions. In the heterogeneous problem, the scalar flux distribution obtained by using discontinuous angular and spatical finite elements is in good agreement with that from the ANISN code calculation. (Author)

Generalized finite elements

International Nuclear Information System (INIS)

Wachspress, E.

2009-01-01

Triangles and rectangles are the ubiquitous elements in finite element studies. Only these elements admit polynomial basis functions. Rational functions provide a basis for elements having any number of straight and curved sides. Numerical complexities initially associated with rational bases precluded extensive use. Recent analysis has reduced these difficulties and programs have been written to illustrate effectiveness. Although incorporation in major finite element software requires considerable effort, there are advantages in some applications which warrant implementation. An outline of the basic theory and of recent innovations is presented here. (authors)

Finite element simulation of thermal, elastic and plastic phenomena in fuel elements

International Nuclear Information System (INIS)

Soba, Alejandro; Denis, Alicia C.

1999-01-01

Taking as starting point an irradiation experiment of the first Argentine MOX fuel prototype, performed at the HFR reactor of Petten, Holland, the deformation suffered by the fuel element materials during burning has been numerically studied. Analysis of the pellet-cladding interaction is made by the finite element method. The code determines the temperature distribution and analyzes elastic and creep deformations, taking into account the dependency of the physical parameters of the problem on temperature. (author)

NACHOS: a finite element computer program for incompressible flow problems. Part I. Theoretical background

International Nuclear Information System (INIS)

Gartling, D.K.

1978-04-01

The theoretical background for the finite element computer program, NACHOS, is presented in detail. The NACHOS code is designed for the two-dimensional analysis of viscous incompressible fluid flows, including the effects of heat transfer. A general description of the fluid/thermal boundary value problems treated by the program is described. The finite element method and the associated numerical methods used in the NACHOS code are also presented. Instructions for use of the program are documented in SAND77-1334

Nonlinear Legendre Spectral Finite Elements for Wind Turbine Blade Dynamics: Preprint

Energy Technology Data Exchange (ETDEWEB)

Wang, Q.; Sprague, M. A.; Jonkman, J.; Johnson, N.

2014-01-01

This paper presents a numerical implementation and examination of new wind turbine blade finite element model based on Geometrically Exact Beam Theory (GEBT) and a high-order spectral finite element method. The displacement-based GEBT is presented, which includes the coupling effects that exist in composite structures and geometric nonlinearity. Legendre spectral finite elements (LSFEs) are high-order finite elements with nodes located at the Gauss-Legendre-Lobatto points. LSFEs can be an order of magnitude more efficient that low-order finite elements for a given accuracy level. Interpolation of the three-dimensional rotation, a major technical barrier in large-deformation simulation, is discussed in the context of LSFEs. It is shown, by numerical example, that the high-order LSFEs, where weak forms are evaluated with nodal quadrature, do not suffer from a drawback that exists in low-order finite elements where the tangent-stiffness matrix is calculated at the Gauss points. Finally, the new LSFE code is implemented in the new FAST Modularization Framework for dynamic simulation of highly flexible composite-material wind turbine blades. The framework allows for fully interactive simulations of turbine blades in operating conditions. Numerical examples showing validation and LSFE performance will be provided in the final paper.

Peridynamic Multiscale Finite Element Methods

Energy Technology Data Exchange (ETDEWEB)

Costa, Timothy [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Bond, Stephen D. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Littlewood, David John [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Moore, Stan Gerald [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2015-12-01

The problem of computing quantum-accurate design-scale solutions to mechanics problems is rich with applications and serves as the background to modern multiscale science research. The prob- lem can be broken into component problems comprised of communicating across adjacent scales, which when strung together create a pipeline for information to travel from quantum scales to design scales. Traditionally, this involves connections between a) quantum electronic structure calculations and molecular dynamics and between b) molecular dynamics and local partial differ- ential equation models at the design scale. The second step, b), is particularly challenging since the appropriate scales of molecular dynamic and local partial differential equation models do not overlap. The peridynamic model for continuum mechanics provides an advantage in this endeavor, as the basic equations of peridynamics are valid at a wide range of scales limiting from the classical partial differential equation models valid at the design scale to the scale of molecular dynamics. In this work we focus on the development of multiscale finite element methods for the peridynamic model, in an effort to create a mathematically consistent channel for microscale information to travel from the upper limits of the molecular dynamics scale to the design scale. In particular, we first develop a Nonlocal Multiscale Finite Element Method which solves the peridynamic model at multiple scales to include microscale information at the coarse-scale. We then consider a method that solves a fine-scale peridynamic model to build element-support basis functions for a coarse- scale local partial differential equation model, called the Mixed Locality Multiscale Finite Element Method. Given decades of research and development into finite element codes for the local partial differential equation models of continuum mechanics there is a strong desire to couple local and nonlocal models to leverage the speed and state of the

Finite element transport methods for criticality calculations - current status and potential applications

International Nuclear Information System (INIS)

Oliveira, C.R.E. de; Goddard, A.

1991-01-01

In this paper we review the current status of the finite element method applied to the solution of the neutron transport equation and we discuss its potential role in the field of criticality safety. We show that the method's ability in handling complex, irregular geometry in two- and three-dimensions coupled with its accurate solutions potentially renders it an attractive alternative to the longer-established Monte Carlo method. Details of the most favoured form of the method - that which combines finite elements in space and spherical harmonics in angle - are presented. This form of the method, which has been extensively investigated over the last decade by research groups at the University of London, has been numerically implemented in the finite element code EVENT. The code has among its main features the capability of solving fixed source eigenvalue and time-dependent complex geometry problems in two- and three-dimensions. Other features of the code include anisotropic up- and down-scatter, direct and/or adjoint solutions and access to standard data libraries. Numerical examples, ranging from simple criticality benchmark studies to the analysis of idealised three-dimensional reactor cores, are presented to demonstrate the potential of the method. (author)

Solution of the diffusion equations for several groups by the finite elements method

International Nuclear Information System (INIS)

Arredondo S, C.

1975-01-01

The code DELFIN has been implemented for the solution of the neutrons diffusion equations in two dimensions obtained by applying the approximation of several groups of energy. The code works with any number of groups and regions, and can be applied to thermal reactors as well as fast reactor. Providing it with the diffusion coefficients, the effective sections and the fission spectrum we obtain the results for the systems multiplying constant and the flows of each groups. The code was established using the method of finite elements, which is a form of resolution of the variational formulation of the equations applying the Ritz-Galerkin method with continuous polynomial functions by parts, in one case of the Lagrange type with rectangular geometry and up to the third grade. The obtained results and the comparison with the results in the literature, permit to reach the conclusion that it is convenient, to use the rectangular elements in all the cases where the geometry permits it, and demonstrate also that the finite elements method is better than the finite differences method. (author)

CHILES, Singularity Strength of Linear Elastic Bodies by Finite Elements Method

International Nuclear Information System (INIS)

Benzley, S.E.; Beisinger, Z.E.

1981-01-01

1 - Description of problem or function: CHILES is a finite element computer program that calculates the strength of singularities in linear elastic bodies. Plane stress, plane strain, and axisymmetric conditions are treated. Crack tip singularity problems are solved by this version of the code, but any type of integrable singularity may be properly modeled by modifying selected subroutines in the program. 2 - Method of solution: A generalized, quadrilateral finite element that includes a singular point at a corner node is incorporated in the code. The displacement formulation is used and inter-element compatibility is maintained so that monotone convergence is preserved. 3 - Restrictions on the complexity of the problem: CHILES allows three singular points to be modeled in the body being analyzed and each singular point may have coupled Mode I and II deformations. 1000 nodal points may be used

Introduction to nonlinear finite element analysis

CERN Document Server

Kim, Nam-Ho

2015-01-01

This book introduces the key concepts of nonlinear finite element analysis procedures. The book explains the fundamental theories of the field and provides instructions on how to apply the concepts to solving practical engineering problems. Instead of covering many nonlinear problems, the book focuses on three representative problems: nonlinear elasticity, elastoplasticity, and contact problems. The book is written independent of any particular software, but tutorials and examples using four commercial programs are included as appendices: ANSYS, NASTRAN, ABAQUS, and MATLAB. In particular, the MATLAB program includes all source codes so that students can develop their own material models, or different algorithms. This book also: ·         Presents clear explanations of nonlinear finite element analysis for elasticity, elastoplasticity, and contact problems ·         Includes many informative examples of nonlinear analyses so that students can clearly understand the nonlinear theory ·    ...

Quality-assurance study of the special-purpose finite-element program SPECTROM: II. Plasticity problems

International Nuclear Information System (INIS)

Callahan, G.D.; Fossum, A.F.

1982-11-01

General plasticity theory and solution techniques as are currently employed in RE/SPEC's finite element plasticity code SPECTROM-II are presented. Various yield functions are discussed and their differences are illustrated using example problems. Comparison of the results of SPECTROM-II with analytical solutions, numerical solutions, and the general purpose finite element program MARC-CDC show excellent agreement

User's manual for the FEHM application - A finite-element heat- and mass-transfer code

International Nuclear Information System (INIS)

Zyvoloski, G.A.; Robinson, B.A.; Dash, Z.V.; Trease, L.L.

1997-07-01

The use of this code is applicable to natural-state studies of geothermal systems and groundwater flow. A primary use of the FEHM application will be to assist in the understanding of flow fields and mass transport in the saturated and unsaturated zones below the proposed Yucca Mountain nuclear waste repository in Nevada. The equations of heat and mass transfer for multiphase flow in porous and permeable media are solved in the FEHM application by using the finite-element method. The permeability and porosity of the medium are allowed to depend on pressure and temperature. The code also has provisions for movable air and water phases and noncoupled tracers; that is, tracer solutions that do not affect the heat- and mass-transfer solutions. The tracers can be passive or reactive. The code can simulate two-dimensional, two-dimensional radial, or three-dimensional geometries. In fact, FEHM is capable of describing flow that is dominated in many areas by fracture and fault flow, including the inherently three-dimensional flow that results from permeation to and from faults and fractures. The code can handle coupled heat and mass-transfer effects, such as boiling, dryout, and condensation that can occur in the near-field region surrounding the potential repository and the natural convection that occurs through Yucca Mountain due to seasonal temperature changes. This report outlines the uses and capabilities of the FEHM application, initialization of code variables, restart procedures, and error processing. The report describes all the data files, the input data, including individual input records or parameters, and the various output files. The system interface is described, including the software environment and installation instructions

Finite element based electric motor design optimization

Science.gov (United States)

Campbell, C. Warren

1993-01-01

The purpose of this effort was to develop a finite element code for the analysis and design of permanent magnet electric motors. These motors would drive electromechanical actuators in advanced rocket engines. The actuators would control fuel valves and thrust vector control systems. Refurbishing the hydraulic systems of the Space Shuttle after each flight is costly and time consuming. Electromechanical actuators could replace hydraulics, improve system reliability, and reduce down time.

Finite element computational fluid mechanics

International Nuclear Information System (INIS)

Baker, A.J.

1983-01-01

This book analyzes finite element theory as applied to computational fluid mechanics. It includes a chapter on using the heat conduction equation to expose the essence of finite element theory, including higher-order accuracy and convergence in a common knowledge framework. Another chapter generalizes the algorithm to extend application to the nonlinearity of the Navier-Stokes equations. Other chapters are concerned with the analysis of a specific fluids mechanics problem class, including theory and applications. Some of the topics covered include finite element theory for linear mechanics; potential flow; weighted residuals/galerkin finite element theory; inviscid and convection dominated flows; boundary layers; parabolic three-dimensional flows; and viscous and rotational flows

Structural evaluation method for class 1 vessels by using elastic-plastic finite element analysis in code case of JSME rules on design and construction

International Nuclear Information System (INIS)

Asada, Seiji; Hirano, Takashi; Nagata, Tetsuya; Kasahara, Naoto

2008-01-01

A structural evaluation method by using elastic-plastic finite element analysis has been developed and published as a code case of Rules on Design and Construction for Nuclear Power Plants (The First Part: Light Water Reactor Structural Design Standard) in the JSME Codes for Nuclear Power Generation Facilities. Its title is 'Alternative Structural Evaluation Criteria for Class 1 Vessels Based on Elastic-Plastic Finite Element Analysis' (NC-CC-005). This code case applies elastic-plastic analysis to evaluation of such failure modes as plastic collapse, thermal ratchet, fatigue and so on. Advantage of this evaluation method is free from stress classification, consistently use of Mises stress and applicability to complex 3-dimensional structures which are hard to be treated by the conventional stress classification method. The evaluation method for plastic collapse has such variation as the Lower Bound Approach Method, Twice-Elastic-Slope Method and Elastic Compensation Method. Cyclic Yield Area (CYA) based on elastic analysis is applied to screening evaluation of thermal ratchet instead of secondary stress evaluation, and elastic-plastic analysis is performed when the CYA screening criteria is not satisfied. Strain concentration factors can be directly calculated based on elastic-plastic analysis. (author)

Monte Carlo method implemented in a finite element code with application to dynamic vacuum in particle accelerators

CERN Document Server

Garion, C

2009-01-01

Modern particle accelerators require UHV conditions during their operation. In the accelerating cavities, breakdowns can occur, releasing large amount of gas into the vacuum chamber. To determine the pressure profile along the cavity as a function of time, the time-dependent behaviour of the gas has to be simulated. To do that, it is useful to apply accurate three-dimensional method, such as Test Particles Monte Carlo. In this paper, a time-dependent Test Particles Monte Carlo is used. It has been implemented in a Finite Element code, CASTEM. The principle is to track a sample of molecules during time. The complex geometry of the cavities can be created either in the FE code or in a CAD software (CATIA in our case). The interface between the two softwares to export the geometry from CATIA to CASTEM is given. The algorithm of particle tracking for collisionless flow in the FE code is shown. Thermal outgassing, pumping surfaces and electron and/or ion stimulated desorption can all be generated as well as differ...

Critical state and magnetization loss in multifilamentary superconducting wire solved through the commercial finite element code ANSYS

Science.gov (United States)

Farinon, S.; Fabbricatore, P.; Gömöry, F.

2010-11-01

The commercially available finite element code ANSYS has been adapted to solve the critical state of single strips and multifilamentary tapes. We studied a special algorithm which approaches the critical state by an iterative adjustment of the material resistivity. Then, we proved its validity by comparing the results obtained for a thin strip to the Brand theory for the transport current and magnetization cases. Also, the challenging calculation of the magnetization loss of a real multifilamentary BSCCO tape showed the usefulness of our method. Finally, we developed several methods to enhance the speed of convergence, making the proposed process quite competitive in the existing survey of ac losses simulations.

Finite element and node point generation computer programs used for the design of toroidal field coils in tokamak fusion devices

International Nuclear Information System (INIS)

Smith, R.A.

1975-06-01

The structural analysis of toroidal field coils in Tokamak fusion machines can be performed with the finite element method. This technique has been employed for design evaluations of toroidal field coils on the Princeton Large Torus (PLT), the Poloidal Diverter Experiment (PDX), and the Tokamak Fusion Test Reactor (TFTR). The application of the finite element method can be simplified with computer programs that are used to generate the input data for the finite element code. There are three areas of data input where significant automation can be provided by supplementary computer codes. These concern the definition of geometry by a node point mesh, the definition of the finite elements from the geometric node points, and the definition of the node point force/displacement boundary conditions. The node point forces in a model of a toroidal field coil are computed from the vector cross product of the coil current and the magnetic field. The computer programs named PDXNODE and ELEMENT are described. The program PDXNODE generates the geometric node points of a finite element model for a toroidal field coil. The program ELEMENT defines the finite elements of the model from the node points and from material property considerations. The program descriptions include input requirements, the output, the program logic, the methods of generating complex geometries with multiple runs, computational time and computer compatibility. The output format of PDXNODE and ELEMENT make them compatible with PDXFORC and two general purpose finite element computer codes: (ANSYS) the Engineering Analysis System written by the Swanson Analysis Systems, Inc., and (WECAN) the Westinghouse Electric Computer Analysis general purpose finite element program. The Fortran listings of PDXNODE and ELEMENT are provided

Finite element simulations with ANSYS workbench 16

CERN Document Server

Lee , Huei-Huang

2015-01-01

Finite Element Simulations with ANSYS Workbench 16 is a comprehensive and easy to understand workbook. It utilizes step-by-step instructions to help guide readers to learn finite element simulations. Twenty seven real world case studies are used throughout the book. Many of these cases are industrial or research projects the reader builds from scratch. All the files readers may need if they have trouble are available for download on the publishers website. Companion videos that demonstrate exactly how to preform each tutorial are available to readers by redeeming the access code that comes in the book. Relevant background knowledge is reviewed whenever necessary. To be efficient, the review is conceptual rather than mathematical. Key concepts are inserted whenever appropriate and summarized at the end of each chapter. Additional exercises or extension research problems are provided as homework at the end of each chapter. A learning approach emphasizing hands-on experiences spreads through this entire book. A...

Finite elements and approximation

CERN Document Server

Zienkiewicz, O C

2006-01-01

A powerful tool for the approximate solution of differential equations, the finite element is extensively used in industry and research. This book offers students of engineering and physics a comprehensive view of the principles involved, with numerous illustrative examples and exercises.Starting with continuum boundary value problems and the need for numerical discretization, the text examines finite difference methods, weighted residual methods in the context of continuous trial functions, and piecewise defined trial functions and the finite element method. Additional topics include higher o

Finite element simulation of piezoelectric transformers.

Science.gov (United States)

Tsuchiya, T; Kagawa, Y; Wakatsuki, N; Okamura, H

2001-07-01

Piezoelectric transformers are nothing but ultrasonic resonators with two pairs of electrodes provided on the surface of a piezoelectric substrate in which electrical energy is carried in the mechanical form. The input and output electrodes are arranged to provide the impedance transformation, which results in the voltage transformation. As they are operated at a resonance, the electrical equivalent circuit approach has traditionally been developed in a rather empirical way and has been used for analysis and design. The present paper deals with the analysis of the piezoelectric transformers based on the three-dimensional finite element modelling. The PIEZO3D code that we have developed is modified to include the external loading conditions. The finite element approach is now available for a wide variety of the electrical boundary conditions. The equivalent circuit of lumped parameters can also be derived from the finite element method (FEM) solution if required. The simulation of the present transformers is made for the low intensity operation and compared with the experimental results. Demonstration is made for basic Rosen-type transformers in which the longitudinal mode of a plate plays an important role; in which the equivalent circuit of lumped constants has been used. However, there are many modes of vibration associated with the plate, the effect of which cannot always be ignored. In the experiment, the double resonances are sometimes observed in the vicinity of the operating frequency. The simulation demonstrates that this is due to the coupling of the longitudinal mode with the flexural mode. Thus, the simulation provides an invaluable guideline to the transformer design.

Modification of the finite element heat and mass transfer code (FEHMN) to model multicomponent reactive transport

International Nuclear Information System (INIS)

Viswanathan, H.S.

1995-01-01

The finite element code FEHMN is a three-dimensional finite element heat and mass transport simulator that can handle complex stratigraphy and nonlinear processes such as vadose zone flow, heat flow and solute transport. Scientists at LANL have been developed hydrologic flow and transport models of the Yucca Mountain site using FEHMN. Previous FEHMN simulations have used an equivalent K d model to model solute transport. In this thesis, FEHMN is modified making it possible to simulate the transport of a species with a rigorous chemical model. Including the rigorous chemical equations into FEHMN simulations should provide for more representative transport models for highly reactive chemical species. A fully kinetic formulation is chosen for the FEHMN reactive transport model. Several methods are available to computationally implement a fully kinetic formulation. Different numerical algorithms are investigated in order to optimize computational efficiency and memory requirements of the reactive transport model. The best algorithm of those investigated is then incorporated into FEHMN. The algorithm chosen requires for the user to place strongly coupled species into groups which are then solved for simultaneously using FEHMN. The complete reactive transport model is verified over a wide variety of problems and is shown to be working properly. The simulations demonstrate that gas flow and carbonate chemistry can significantly affect 14 C transport at Yucca Mountain. The simulations also provide that the new capabilities of FEHMN can be used to refine and buttress already existing Yucca Mountain radionuclide transport studies

Code conforming determination of cumulative usage factors for general elastic-plastic finite element analyses

International Nuclear Information System (INIS)

Rudolph, Juergen; Goetz, Andreas; Hilpert, Roland

2012-01-01

The procedures of fatigue analyses of several relevant nuclear and conventional design codes (ASME, KTA, EN, AD) for power plant components differentiate between an elastic, simplified elastic-plastic and elastic-plastic fatigue check. As a rule, operational load levels will exclude the purely elastic fatigue check. The application of the code procedure of the simplified elastic-plastic fatigue check is common practice. Nevertheless, resulting cumulative usage factors may be overly conservative mainly due to high code based plastification penalty factors Ke. As a consequence, the more complex and still code conforming general elastic-plastic fatigue analysis methodology based on non-linear finite element analysis (FEA) is applied for fatigue design as an alternative. The requirements of the FEA and the material law to be applied have to be clarified in a first step. Current design codes only give rough guidelines on these relevant items. While the procedure for the simplified elastic-plastic fatigue analysis and the associated code passages are based on stress related cycle counting and the determination of pseudo elastic equivalent stress ranges, an adaptation to elastic-plastic strains and strain ranges is required for the elastic-plastic fatigue check. The associated requirements are explained in detail in the paper. If the established and implemented evaluation mechanism (cycle counting according to the peak and valley respectively the rainflow method, calculation of stress ranges from arbitrary load-time histories and determination of cumulative usage factors based on all load events) is to be retained, a conversion of elastic-plastic strains and strain ranges into pseudo elastic stress ranges is required. The algorithm to be applied is described in the paper. It has to be implemented in the sense of an extended post processing operation of FEA e.g. by APDL scripts in ANSYS registered . Variations of principal stress (strain) directions during the loading

User`s manual for the FEHM application -- A finite-element heat- and mass-transfer code

Energy Technology Data Exchange (ETDEWEB)

Zyvoloski, G.A.; Robinson, B.A.; Dash, Z.V.; Trease, L.L.

1997-07-01

The use of this code is applicable to natural-state studies of geothermal systems and groundwater flow. A primary use of the FEHM application will be to assist in the understanding of flow fields and mass transport in the saturated and unsaturated zones below the proposed Yucca Mountain nuclear waste repository in Nevada. The equations of heat and mass transfer for multiphase flow in porous and permeable media are solved in the FEHM application by using the finite-element method. The permeability and porosity of the medium are allowed to depend on pressure and temperature. The code also has provisions for movable air and water phases and noncoupled tracers; that is, tracer solutions that do not affect the heat- and mass-transfer solutions. The tracers can be passive or reactive. The code can simulate two-dimensional, two-dimensional radial, or three-dimensional geometries. In fact, FEHM is capable of describing flow that is dominated in many areas by fracture and fault flow, including the inherently three-dimensional flow that results from permeation to and from faults and fractures. The code can handle coupled heat and mass-transfer effects, such as boiling, dryout, and condensation that can occur in the near-field region surrounding the potential repository and the natural convection that occurs through Yucca Mountain due to seasonal temperature changes. This report outlines the uses and capabilities of the FEHM application, initialization of code variables, restart procedures, and error processing. The report describes all the data files, the input data, including individual input records or parameters, and the various output files. The system interface is described, including the software environment and installation instructions.

Finite element modeling of TFTR poloidal field coils

International Nuclear Information System (INIS)

Baumgartner, J.A.; O'Toole, J.A.

1986-01-01

The Tokamak Fusion Test Reactor (TFTR) Poloidal Field (PF) coils were originally analyzed to TFTR design conditions. The coils have been reanalyzed by PPPL and Grumman to determine operating limits under as-built conditions. Critical stress levels, based upon data obtained from the reanalysis of each PF coil, are needed for input to the TFTR simulation code algorithms. The primary objective regarding structural integrity has been to ascertain the magnitude and location of critical internal stresses in each PF coil due to various combinations of electromagnetic and thermally induced loads. For each PF coil, a global finite element model (FEM) of a coil sector is being analyzed to obtain the basic coil internal loads and displacements. Subsequent fine mesh local models of the coil lead stem and lead spur regions produce the magnitudes and locations of peak stresses. Each copper turn and its surrounding insulation are modeled using solid finite elements. The corresponding electromagnetic and thermal analyses are similarly modeled. A series of test beams were developed to determine the best combination of MSC/NASTRAN-type finite elements for use in PF coil analysis. The results of this analysis compare favorably with those obtained by the earlier analysis which was limited in scope

Asymmetric Rolling Process Simulations by Dynamic Explicit Crystallographic Homogenized Finite Element Method

International Nuclear Information System (INIS)

Ngoc Tam, Nguyen; Nakamura, Yasunori; Terao, Toshihiro; Kuramae, Hiroyuki; Nakamachi, Eiji; Sakamoto, Hidetoshi; Morimoto, Hideo

2007-01-01

Recently, the asymmetric rolling (ASR) has been applied to the material processing of aluminum alloy sheet to control micro-crystal structure and texture in order to improve the mechanical properties. Previously, several studies aimed at high formability sheet generation have been carried out experimentally, but finite element simulations to predict the deformation induced texture evolution of the asymmetrically rolled sheet metals have not been investigated rigorously. In this study, crystallographic homogenized finite element (FE) codes are developed and applied to analyze the asymmetrical rolling processes. The textures of sheet metals were measured by electron back scattering diffraction (EBSD), and compared with FE simulations. The results from the dynamic explicit type Crystallographic homogenization FEM code shows that this type of simulation is a comprehensive tool to predict the plastic induced texture evolution

Modelling optimization involving different types of elements in finite element analysis

International Nuclear Information System (INIS)

Wai, C M; Rivai, Ahmad; Bapokutty, Omar

2013-01-01

Finite elements are used to express the mechanical behaviour of a structure in finite element analysis. Therefore, the selection of the elements determines the quality of the analysis. The aim of this paper is to compare and contrast 1D element, 2D element, and 3D element used in finite element analysis. A simple case study was carried out on a standard W460x74 I-beam. The I-beam was modelled and analyzed statically with 1D elements, 2D elements and 3D elements. The results for the three separate finite element models were compared in terms of stresses, deformation and displacement of the I-beam. All three finite element models yield satisfactory results with acceptable errors. The advantages and limitations of these elements are discussed. 1D elements offer simplicity although lacking in their ability to model complicated geometry. 2D elements and 3D elements provide more detail yet sophisticated results which require more time and computer memory in the modelling process. It is also found that the choice of element in finite element analysis is influence by a few factors such as the geometry of the structure, desired analysis results, and the capability of the computer

Modelling the attenuation in the ATHENA finite elements code for the ultrasonic testing of austenitic stainless steel welds.

Science.gov (United States)

Chassignole, B; Duwig, V; Ploix, M-A; Guy, P; El Guerjouma, R

2009-12-01

Multipass welds made in austenitic stainless steel, in the primary circuit of nuclear power plants with pressurized water reactors, are characterized by an anisotropic and heterogeneous structure that disturbs the ultrasonic propagation and makes ultrasonic non-destructive testing difficult. The ATHENA 2D finite element simulation code was developed to help understand the various physical phenomena at play. In this paper, we shall describe the attenuation model implemented in this code to give an account of wave scattering phenomenon through polycrystalline materials. This model is in particular based on the optimization of two tensors that characterize this material on the basis of experimental values of ultrasonic velocities attenuation coefficients. Three experimental configurations, two of which are representative of the industrial welds assessment case, are studied in view of validating the model through comparison with the simulation results. We shall thus provide a quantitative proof that taking into account the attenuation in the ATHENA code dramatically improves the results in terms of the amplitude of the echoes. The association of the code and detailed characterization of a weld's structure constitutes a remarkable breakthrough in the interpretation of the ultrasonic testing on this type of component.

A mixed finite element method for particle simulation in lasertron

International Nuclear Information System (INIS)

Le Meur, G.

1987-03-01

A particle simulation code is being developed with the aim to treat the motion of charged particles in electromagnetic devices, such as Lasertron. The paper describes the use of mixed finite element methods in computing the field components, without derivating them from scalar or vector potentials. Graphical results are shown

A mixed finite element method for particle simulation in Lasertron

International Nuclear Information System (INIS)

Le Meur, G.

1987-01-01

A particle simulation code is being developed with the aim to treat the motion of charged particles in electromagnetic devices, such as Lasertron. The paper describes the use of mixed finite element methods in computing the field components, without derivating them from scalar or vector potentials. Graphical results are shown

SANTOS - a two-dimensional finite element program for the quasistatic, large deformation, inelastic response of solids

Energy Technology Data Exchange (ETDEWEB)

Stone, C.M.

1997-07-01

SANTOS is a finite element program designed to compute the quasistatic, large deformation, inelastic response of two-dimensional planar or axisymmetric solids. The code is derived from the transient dynamic code PRONTO 2D. The solution strategy used to compute the equilibrium states is based on a self-adaptive dynamic relaxation solution scheme, which is based on explicit central difference pseudo-time integration and artificial mass proportional damping. The element used in SANTOS is a uniform strain 4-node quadrilateral element with an hourglass control scheme to control the spurious deformation modes. Finite strain constitutive models for many common engineering materials are included. A robust master-slave contact algorithm for modeling sliding contact is implemented. An interface for coupling to an external code is also provided. 43 refs., 22 figs.

MP Salsa: a finite element computer program for reacting flow problems. Part 1--theoretical development

Energy Technology Data Exchange (ETDEWEB)

Shadid, J.N.; Moffat, H.K.; Hutchinson, S.A.; Hennigan, G.L.; Devine, K.D.; Salinger, A.G.

1996-05-01

The theoretical background for the finite element computer program, MPSalsa, is presented in detail. MPSalsa is designed to solve laminar, low Mach number, two- or three-dimensional incompressible and variable density reacting fluid flows on massively parallel computers, using a Petrov-Galerkin finite element formulation. The code has the capability to solve coupled fluid flow, heat transport, multicomponent species transport, and finite-rate chemical reactions, and to solver coupled multiple Poisson or advection-diffusion- reaction equations. The program employs the CHEMKIN library to provide a rigorous treatment of multicomponent ideal gas kinetics and transport. Chemical reactions occurring in the gas phase and on surfaces are treated by calls to CHEMKIN and SURFACE CHEMKIN, respectively. The code employs unstructured meshes, using the EXODUS II finite element data base suite of programs for its input and output files. MPSalsa solves both transient and steady flows by using fully implicit time integration, an inexact Newton method and iterative solvers based on preconditioned Krylov methods as implemented in the Aztec solver library.

Automating the generation of finite element dynamical cores with Firedrake

Science.gov (United States)

Ham, David; Mitchell, Lawrence; Homolya, Miklós; Luporini, Fabio; Gibson, Thomas; Kelly, Paul; Cotter, Colin; Lange, Michael; Kramer, Stephan; Shipton, Jemma; Yamazaki, Hiroe; Paganini, Alberto; Kärnä, Tuomas

2017-04-01

The development of a dynamical core is an increasingly complex software engineering undertaking. As the equations become more complete, the discretisations more sophisticated and the hardware acquires ever more fine-grained parallelism and deeper memory hierarchies, the problem of building, testing and modifying dynamical cores becomes increasingly complex. Here we present Firedrake, a code generation system for the finite element method with specialist features designed to support the creation of geoscientific models. Using Firedrake, the dynamical core developer writes the partial differential equations in weak form in a high level mathematical notation. Appropriate function spaces are chosen and time stepping loops written at the same high level. When the programme is run, Firedrake generates high performance C code for the resulting numerics which are executed in parallel. Models in Firedrake typically take a tiny fraction of the lines of code required by traditional hand-coding techniques. They support more sophisticated numerics than are easily achieved by hand, and the resulting code is frequently higher performance. Critically, debugging, modifying and extending a model written in Firedrake is vastly easier than by traditional methods due to the small, highly mathematical code base. Firedrake supports a wide range of key features for dynamical core creation: A vast range of discretisations, including both continuous and discontinuous spaces and mimetic (C-grid-like) elements which optimally represent force balances in geophysical flows. High aspect ratio layered meshes suitable for ocean and atmosphere domains. Curved elements for high accuracy representations of the sphere. Support for non-finite element operators, such as parametrisations. Access to PETSc, a world-leading library of programmable linear and nonlinear solvers. High performance adjoint models generated automatically by symbolically reasoning about the forward model. This poster will present

Finite element modeling of piezoelectric elements with complex electrode configuration

International Nuclear Information System (INIS)

Paradies, R; Schläpfer, B

2009-01-01

It is well known that the material properties of piezoelectric materials strongly depend on the state of polarization of the individual element. While an unpolarized material exhibits mechanically isotropic material properties in the absence of global piezoelectric capabilities, the piezoelectric material properties become transversally isotropic with respect to the polarization direction after polarization. Therefore, for evaluating piezoelectric elements the material properties, including the coupling between the mechanical and the electromechanical behavior, should be addressed correctly. This is of special importance for the micromechanical description of piezoelectric elements with interdigitated electrodes (IDEs). The best known representatives of this group are active fiber composites (AFCs), macro fiber composites (MFCs) and the radial field diaphragm (RFD), respectively. While the material properties are available for a piezoelectric wafer with a homogeneous polarization perpendicular to its plane as postulated in the so-called uniform field model (UFM), the same information is missing for piezoelectric elements with more complex electrode configurations like the above-mentioned ones with IDEs. This is due to the inhomogeneous field distribution which does not automatically allow for the correct assignment of the material, i.e. orientation and property. A variation of the material orientation as well as the material properties can be accomplished by including the polarization process of the piezoelectric transducer in the finite element (FE) simulation prior to the actual load case to be investigated. A corresponding procedure is presented which automatically assigns the piezoelectric material properties, e.g. elasticity matrix, permittivity, and charge vector, for finite element models (FEMs) describing piezoelectric transducers according to the electric field distribution (field orientation and strength) in the structure. A corresponding code has been

Numerical modelling of the long-term evolution of EDZ. Development of material models, implementation in finite-element codes, and validation

International Nuclear Information System (INIS)

Pudewills, A.

2005-11-01

Construction of deep underground structures disturbs the initial stress field in the surrounding rock. This effect can generate microcracks and alter the hydromechanical properties of the rock salt around the excavations. For the long-term performance of an underground repository in rock salt, the evolution of the 'Excavation Disturbed Zone' (EDZ) and the hydromechanical behaviour of this zone represent important issues with respect to the integrity of the geological and technical barriers. Within the framework of the NF-PRO project, WP 4.4, attention focuses on the mathematical modelling of the development and evolution of the EDZ in the rock near a disposal drift due to its relevance on the integrity of the geological and technical barriers. To perform this task, finite-element codes containing a set of time- and temperature-dependent constitutive models have been improved. A new viscoplastic constitutive model for rock salt that can describe the damage of the rock has been implemented in the finite-element codes available. The model parameters were evaluated based on experimental results. Additionally, the long-term evolution of the EDZ around a gallery in a salt mine at about 700 m below the surface was analysed and the numerical results were compared with in-situ measurements. The calculated room closure, stress distribution and the increase of rock permeability in the EDZ were compared with in-situ data, thus providing confidence in the model used. (orig.)

Non-linear finite element modeling

DEFF Research Database (Denmark)

Mikkelsen, Lars Pilgaard

The note is written for courses in "Non-linear finite element method". The note has been used by the author teaching non-linear finite element modeling at Civil Engineering at Aalborg University, Computational Mechanics at Aalborg University Esbjerg, Structural Engineering at the University...

User's manual for DYNA2D: an explicit two-dimensional hydrodynamic finite-element code with interactive rezoning

Energy Technology Data Exchange (ETDEWEB)

Hallquist, J.O.

1982-02-01

This revised report provides an updated user's manual for DYNA2D, an explicit two-dimensional axisymmetric and plane strain finite element code for analyzing the large deformation dynamic and hydrodynamic response of inelastic solids. A contact-impact algorithm permits gaps and sliding along material interfaces. By a specialization of this algorithm, such interfaces can be rigidly tied to admit variable zoning without the need of transition regions. Spatial discretization is achieved by the use of 4-node solid elements, and the equations-of motion are integrated by the central difference method. An interactive rezoner eliminates the need to terminate the calculation when the mesh becomes too distorted. Rather, the mesh can be rezoned and the calculation continued. The command structure for the rezoner is described and illustrated by an example.

Structural modeling techniques by finite element method

International Nuclear Information System (INIS)

Kang, Yeong Jin; Kim, Geung Hwan; Ju, Gwan Jeong

1991-01-01

This book includes introduction table of contents chapter 1 finite element idealization introduction summary of the finite element method equilibrium and compatibility in the finite element solution degrees of freedom symmetry and anti symmetry modeling guidelines local analysis example references chapter 2 static analysis structural geometry finite element models analysis procedure modeling guidelines references chapter 3 dynamic analysis models for dynamic analysis dynamic analysis procedures modeling guidelines and modeling guidelines.

An evaluation of a translator for finite element data to resistor/capacitor data for the heat diffusion equation

International Nuclear Information System (INIS)

Manteufel, R.D.; Klein, D.E.; Yoshimura, H.R.

1988-01-01

This paper evaluates a translator for finite element data to resistor/capacitor data (FEM/RC) for the numerical solution of heat diffusion problems. The translator involves the derivation of thermal resistors and capacitors, implicit in the heat balance formulation of the finite difference method. It uses a finite element mesh, which consists of nodes and elements and is implicit in the Galerkin finite element method (GFEM). This hybrid translation method, FEM/RC, has been incorporated in Q/TRAN, a new thermal analysis computer code. This evaluation compares Q/TRAN, HEATING-6, and a research code employing GFEM on a purely mathematical, highly nonlinear steady-state conduction benchmark problem. The evaluation concludes that the FEM/RC technique has numerical characteristics that are consistent with comparable schemes for the benchmark problem. FEM/RC also accurately translates skewed meshes. Because FEM/RC generates resistors and capacitors, it appears to offer a more efficient method than the classical GFEM

Element-topology-independent preconditioners for parallel finite element computations

Science.gov (United States)

Park, K. C.; Alexander, Scott

1992-01-01

A family of preconditioners for the solution of finite element equations are presented, which are element-topology independent and thus can be applicable to element order-free parallel computations. A key feature of the present preconditioners is the repeated use of element connectivity matrices and their left and right inverses. The properties and performance of the present preconditioners are demonstrated via beam and two-dimensional finite element matrices for implicit time integration computations.

Discrete-ordinates finite-element method for atmospheric radiative transfer and remote sensing

International Nuclear Information System (INIS)

Gerstl, S.A.W.; Zardecki, A.

1985-01-01

Advantages and disadvantages of modern discrete-ordinates finite-element methods for the solution of radiative transfer problems in meteorology, climatology, and remote sensing applications are evaluated. After the common basis of the formulation of radiative transfer problems in the fields of neutron transport and atmospheric optics is established, the essential features of the discrete-ordinates finite-element method are described including the limitations of the method and their remedies. Numerical results are presented for 1-D and 2-D atmospheric radiative transfer problems where integral as well as angular dependent quantities are compared with published results from other calculations and with measured data. These comparisons provide a verification of the discrete-ordinates results for a wide spectrum of cases with varying degrees of absorption, scattering, and anisotropic phase functions. Accuracy and computational speed are also discussed. Since practically all discrete-ordinates codes offer a builtin adjoint capability, the general concept of the adjoint method is described and illustrated by sample problems. Our general conclusion is that the strengths of the discrete-ordinates finite-element method outweight its weaknesses. We demonstrate that existing general-purpose discrete-ordinates codes can provide a powerful tool to analyze radiative transfer problems through the atmosphere, especially when 2-D geometries must be considered

Finite element simulation of dynamic wetting flows as an interface formation process

KAUST Repository

Sprittles, J.E.; Shikhmurzaev, Y.D.

2013-01-01

A mathematically challenging model of dynamic wetting as a process of interface formation has been, for the first time, fully incorporated into a numerical code based on the finite element method and applied, as a test case, to the problem

Co-simulation coupling spectral/finite elements for 3D soil/structure interaction problems

Science.gov (United States)

Zuchowski, Loïc; Brun, Michael; De Martin, Florent

2018-05-01

The coupling between an implicit finite elements (FE) code and an explicit spectral elements (SE) code has been explored for solving the elastic wave propagation in the case of soil/structure interaction problem. The coupling approach is based on domain decomposition methods in transient dynamics. The spatial coupling at the interface is managed by a standard coupling mortar approach, whereas the time integration is dealt with an hybrid asynchronous time integrator. An external coupling software, handling the interface problem, has been set up in order to couple the FE software Code_Aster with the SE software EFISPEC3D.

Probabilistic finite elements

Science.gov (United States)

Belytschko, Ted; Wing, Kam Liu

1987-01-01

In the Probabilistic Finite Element Method (PFEM), finite element methods have been efficiently combined with second-order perturbation techniques to provide an effective method for informing the designer of the range of response which is likely in a given problem. The designer must provide as input the statistical character of the input variables, such as yield strength, load magnitude, and Young's modulus, by specifying their mean values and their variances. The output then consists of the mean response and the variance in the response. Thus the designer is given a much broader picture of the predicted performance than with simply a single response curve. These methods are applicable to a wide class of problems, provided that the scale of randomness is not too large and the probabilistic density functions possess decaying tails. By incorporating the computational techniques we have developed in the past 3 years for efficiency, the probabilistic finite element methods are capable of handling large systems with many sources of uncertainties. Sample results for an elastic-plastic ten-bar structure and an elastic-plastic plane continuum with a circular hole subject to cyclic loadings with the yield stress on the random field are given.

Finite element and finite difference methods in electromagnetic scattering

CERN Document Server

Morgan, MA

2013-01-01

This second volume in the Progress in Electromagnetic Research series examines recent advances in computational electromagnetics, with emphasis on scattering, as brought about by new formulations and algorithms which use finite element or finite difference techniques. Containing contributions by some of the world's leading experts, the papers thoroughly review and analyze this rapidly evolving area of computational electromagnetics. Covering topics ranging from the new finite-element based formulation for representing time-harmonic vector fields in 3-D inhomogeneous media using two coupled sca

Finite-Element 2D and 3D PIC Modeling of RF Devices with Applications to Multipacting

CERN Document Server

De Ford, John F; Petillo, John

2005-01-01

Multipacting currently limits the performance of many high power radio-frequency (RF) devices, particularly couplers and windows. Models have helped researchers understand and mitigate this problem in 2D structures, but useful multipacting models for complicated 3D structures are still a challenge. A combination of three recent technologies that have been developed in the Analyst and MICHELLE codes begin to address this challenge: high-order adaptive finite-element RF field calculations, advanced particle tracking on unstructured grids, and comprehensive secondary emission models. Analyst employs high-order adaptive finite-element methods to accurately compute driven RF fields and eigenmodes in complex geometries, particularly near edges, corners, and curved surfaces. To perform a multipacting analysis, we use the mesh and fields from Analyst in a modified version of the self-consistent, finite-element gun code MICHELLE. MICHELLE has both a fast, accurate, and reliable particle tracker for unstructured grids ...

Electrical machine analysis using finite elements

CERN Document Server

Bianchi, Nicola

2005-01-01

OUTLINE OF ELECTROMAGNETIC FIELDSVector AnalysisElectromagnetic FieldsFundamental Equations SummaryReferencesBASIC PRINCIPLES OF FINITE ELEMENT METHODSIntroductionField Problems with Boundary ConditionsClassical Method for the Field Problem SolutionThe Classical Residual Method (Galerkin's Method)The Classical Variational Method (Rayleigh-Ritz's Method)The Finite Element MethodReferencesAPPLICATIONS OF THE FINITE ELEMENT METHOD TO TWO-DIMENSIONAL FIELDSIntroductionLinear Interpolation of the Function fApplication of the Variational MethodSimple Descriptions of Electromagnetic FieldsAppendix: I

Analysis of piping systems by finite element method using code SAP-IV

International Nuclear Information System (INIS)

Cizelj, L.; Ogrizek, D.

1987-01-01

Due to extensive and multiple use of the computer code SAP-IV we have decided to install it on VAX 11/750 machine. Installation required a large quantity of programming due to great discrepancies between the CDC (the original program version) and the VAX. Testing was performed basically in the field of pipe elements, based on a comparison between results obtained with the codes PSAFE2, DOCIJEV, PIPESD and SAP -V. Besides, the model of reactor pressure vessel with 3-D thick shell elements was done. The capabilities show good agreement with the results of other programs mentioned above. Along with the package installation, the graphical postprocessors being developed for mesh plotting. (author)

Finite Element Analysis of Dam-Reservoir Interaction Using High-Order Doubly Asymptotic Open Boundary

Directory of Open Access Journals (Sweden)

Yichao Gao

2011-01-01

Full Text Available The dam-reservoir system is divided into the near field modeled by the finite element method, and the far field modeled by the excellent high-order doubly asymptotic open boundary (DAOB. Direct and partitioned coupled methods are developed for the analysis of dam-reservoir system. In the direct coupled method, a symmetric monolithic governing equation is formulated by incorporating the DAOB with the finite element equation and solved using the standard time-integration methods. In contrast, the near-field finite element equation and the far-field DAOB condition are separately solved in the partitioned coupled methodm, and coupling is achieved by applying the interaction force on the truncated boundary. To improve its numerical stability and accuracy, an iteration strategy is employed to obtain the solution of each step. Both coupled methods are implemented on the open-source finite element code OpenSees. Numerical examples are employed to demonstrate the performance of these two proposed methods.

Application of finite-element method to three-dimensional nuclear reactor analysis

International Nuclear Information System (INIS)

Cheung, K.Y.

1985-01-01

The application of the finite element method to solve a realistic one-or-two energy group, multiregion, three-dimensional static neutron diffusion problem is studied. Linear, quadratic, and cubic serendipity box-shape elements are used. The resulting sets of simultaneous algebraic equations with thousands of unknowns are solved by the conjugate gradient method, without forming the large coefficient matrix explicitly. This avoids the complicated data management schemes to store such a large coefficient matrix. Three finite-element computer programs: FEM-LINEAR, FEM-QUADRATIC and FEM-CUBIC were developed, using the linear, quadratic, and cubic box-shape elements respectively. They are self-contained, using simple nodal labeling schemes, without the need for separate finite element mesh generating routines. The efficiency and accuracy of these computer programs are then compared among themselves, and with other computer codes. The cubic element model is not recommended for practical usage because it gives almost identical results as the quadratic model, but it requires considerably longer computation time. The linear model is less accurate than the quadratic model, but it requires much shorter computation time. For a large 3-D problem, the linear model is to be preferred since it gives acceptable accuracy. The quadratic model may be used if improved accuracy is desired

Evaluation of the finite element software ABAQUS for biomechanical modelling of biphasic tissues.

Science.gov (United States)

Wu, J Z; Herzog, W; Epstein, M

1998-02-01

The biphasic cartilage model proposed by Mow et al. (1980) has proven successful to capture the essential mechanical features of articular cartilage. In order to analyse the joint contact mechanics in real, anatomical joints, the cartilage model needs to be implemented into a suitable finite element code to approximate the irregular surface geometries of such joints. However, systematic and extensive evaluation of the capacity of commercial software for modelling the contact mechanics with biphasic cartilage layers has not been made. This research was aimed at evaluating the commercial finite element software ABAQUS for analysing biphasic soft tissues. The solutions obtained using ABAQUS were compared with those obtained using other finite element models and analytical solutions for three numerical tests: an unconfined indentation test, a test with the contact of a spherical cartilage surface with a rigid plate, and an axi-symmetric joint contact test. It was concluded that the biphasic cartilage model can be implemented into the commercial finite element software ABAQUS to analyse practical joint contact problems with biphasic articular cartilage layers.

Development of polygon elements based on the scaled boundary finite element method

International Nuclear Information System (INIS)

Chiong, Irene; Song Chongmin

2010-01-01

We aim to extend the scaled boundary finite element method to construct conforming polygon elements. The development of the polygonal finite element is highly anticipated in computational mechanics as greater flexibility and accuracy can be achieved using these elements. The scaled boundary polygonal finite element will enable new developments in mesh generation, better accuracy from a higher order approximation and better transition elements in finite element meshes. Polygon elements of arbitrary number of edges and order have been developed successfully. The edges of an element are discretised with line elements. The displacement solution of the scaled boundary finite element method is used in the development of shape functions. They are shown to be smooth and continuous within the element, and satisfy compatibility and completeness requirements. Furthermore, eigenvalue decomposition has been used to depict element modes and outcomes indicate the ability of the scaled boundary polygonal element to express rigid body and constant strain modes. Numerical tests are presented; the patch test is passed and constant strain modes verified. Accuracy and convergence of the method are also presented and the performance of the scaled boundary polygonal finite element is verified on Cook's swept panel problem. Results show that the scaled boundary polygonal finite element method outperforms a traditional mesh and accuracy and convergence are achieved from fewer nodes. The proposed method is also shown to be truly flexible, and applies to arbitrary n-gons formed of irregular and non-convex polygons.

An object-oriented class design for the generalized finite element method programming

Directory of Open Access Journals (Sweden)

Dorival Piedade Neto

Full Text Available The Generalized Finite Element Method (GFEM is a numerical method based on the Finite Element Method (FEM, presenting as its main feature the possibility of improving the solution by means of local enrichment functions. In spite of its advantages, the method demands a complex data structure, which can be especially benefited by the Object-Oriented Programming (OOP. Even though the OOP for the traditional FEM has been extensively described in the technical literature, specific design issues related to the GFEM are yet little discussed and not clearly defined. In the present article it is described an Object-Oriented (OO class design for the GFEM, aiming to achieve a computational code that presents a flexible class structure, circumventing the difficulties associated to the method characteristics. The proposed design is evaluated by means of some numerical examples, computed using a code implemented in Python programming language.

Static Analysis of Steel Fiber Concrete Beam With Heterosis Finite Elements

Directory of Open Access Journals (Sweden)

James H. Haido

2014-08-01

Full Text Available Steel fiber is considered as the most commonly used constructional fibers in concrete structures. The formulation of new nonlinearities to predict the static performance of steel fiber concrete composite structures is considered essential. Present study is devoted to investigate the efficiency of utilizing heterosis finite elements analysis in static analysis of steel fibrous beams. New and simple material nonlinearities are proposed and used in the formulation of these elements. A computer program coded in FORTRAN was developed to perform current finite element static analysis with considering four cases of elements stiffness matrix determination. The results are compared with the experimental data available in literature in terms of central deflections, strains, and failure form, good agreement was found. Suitable outcomes have been observed in present static analysis with using of tangential stiffness matrix and stiffness matrix in second iteration of the load increment.

Finite element circuit theory of the numerical code EDDYMULT for solving eddy current problems in a multi-torus system

International Nuclear Information System (INIS)

Nakamura, Yukiharu; Ozeki, Takahisa

1986-07-01

The finite element circuit theory is extended to the general eddy current problem in a multi-torus system, which consists of various torus conductors and axisymmetric coil systems. The numerical procedures are devised to avoid practical restrictions of computer storage and computing time, that is, the reduction technique of eddy current eigen modes to save storage and the introduction of shape function into the double area integral of mode coupling to save time. The numerical code EDDYMULT based on the theory is developed to use in designing tokamak device from the viewpoints of the evaluation of electromagnetic loading on the device components and the control analysis of tokamak equilibrium. (author)

Development of finite element code for the analysis of coupled thermo-hydro-mechanical behaviors of a saturated-unsaturated medium

International Nuclear Information System (INIS)

Ohnishi, Y.; Shibata, H.; Kobsayashi, A.

1987-01-01

A model is presented which describes fully coupled thermo-hydro-mechanical behavior of a porous geologic medium. The mathematical formulation for the model utilizes the Biot theory for the consolidation and the energy balance equation. If the medium is in the condition of saturated-unsaturated flow, then the free surfaces are taken into consideration in the model. The model, incorporated in a finite element numerical procedure, was implemented in a two-dimensional computer code. The code was developed under the assumptions that the medium is poro-elastic and in the plane strain condition; that water in the ground does not change its phase; and that heat is transferred by conductive and convective flow. Analytical solutions pertaining to consolidation theory for soils and rocks, thermoelasticity for solids and hydrothermal convection theory provided verification of stress and fluid flow couplings, respectively, in the coupled model. Several types of problems are analyzed

COYOTE: a finite element computer program for nonlinear heat conduction problems

International Nuclear Information System (INIS)

Gartling, D.K.

1978-06-01

COYOTE is a finite element computer program designed for the solution of two-dimensional, nonlinear heat conduction problems. The theoretical and mathematical basis used to develop the code is described. Program capabilities and complete user instructions are presented. Several example problems are described in detail to demonstrate the use of the program

TAURUS, Post-processor of 3-D Finite Elements Plots

International Nuclear Information System (INIS)

Brown, B.E.; Hallquist, J.O.; Kennedy, T.

2002-01-01

Description of program or function: TAURUS reads the binary plot files generated by the LLNL three-dimensional finite element analysis codes, NIKE3D (NESC 9725), DYNA3D (NESC 9909), TACO3D (NESC 9838), TOPAZ3D (NESC9599) and GEMINI and plots contours, time histories, and deformed shapes. Contours of a large number of quantities may be plotted on meshes consisting of plate, shell, and solid type elements. TAURUS can compute a variety of strain measures, reaction forces along constrained boundaries, and momentum. TAURUS has three phases: initialization, geometry display with contouring, and time history processing

FINITE ELEMENT MODEL FOR PREDICTING RESIDUAL ...

African Journals Online (AJOL)

FINITE ELEMENT MODEL FOR PREDICTING RESIDUAL STRESSES IN ... the transverse residual stress in the x-direction (σx) had a maximum value of 375MPa ... the finite element method are in fair agreement with the experimental results.

A non-conformal finite element/finite volume scheme for the non-structured grid-based approximation of low Mach number flows

International Nuclear Information System (INIS)

Ansanay-Alex, G.

2009-01-01

The development of simulation codes aimed at a precise simulation of fires requires a precise approach of flame front phenomena by using very fine grids. The need to take different spatial scale into consideration leads to a local grid refinement and to a discretization with homogeneous grid for computing time and memory purposes. The author reports the approximation of the non-linear convection term, the scalar advection-diffusion in finite volumes, numerical simulations of a flow in a bent tube, of a three-dimensional laminar flame and of a low Mach number an-isotherm flow. Non conformal finite elements are also presented (Rannacher-Turek and Crouzeix-Raviart elements)

A 3D analysis of reinforced concrete structures by the finite element method

International Nuclear Information System (INIS)

Claure, J.D.; Campos Filho, A.

1995-01-01

Fundamental features of a computational model, based on the finite element methods, for the analysis of concrete structure are presented. The study comprehends short and long-term loading situations, where creep and shrinkage in concrete are considered. The reinforcement is inserted in the finite element model using an embedded model. A smeared crack model is used for the concrete cracking, which considers the contribution of concrete between cracks and allows the closing the cracks closing. The computational code MPGS (Multi-Purpose Graphic System) is used, to make easy the analysis and interpretation of the numeric results. (author). 8 refs., 4 figs

Finite element methods a practical guide

CERN Document Server

Whiteley, Jonathan

2017-01-01

This book presents practical applications of the finite element method to general differential equations. The underlying strategy of deriving the finite element solution is introduced using linear ordinary differential equations, thus allowing the basic concepts of the finite element solution to be introduced without being obscured by the additional mathematical detail required when applying this technique to partial differential equations. The author generalizes the presented approach to partial differential equations which include nonlinearities. The book also includes variations of the finite element method such as different classes of meshes and basic functions. Practical application of the theory is emphasised, with development of all concepts leading ultimately to a description of their computational implementation illustrated using Matlab functions. The target audience primarily comprises applied researchers and practitioners in engineering, but the book may also be beneficial for graduate students.

Finite Element Analysis of the Propagation of Acoustic Waves Along Waveguides Immersed in Water

Science.gov (United States)

Hladky-Hennion, A.-C.; Langlet, P.; de Billy, M.

1997-03-01

The finite element approach has previously been used, with the help of the ATILA code, to model the propagation of acoustic waves in waveguides [A.-C. Hladky-Hennion, Journal of Sound and Vibration, 194,119-136 (1996)]. In this paper an extension of the technique to the analysis of the propagation of acoustic waves in immersed waveguides is presented. In the proposed approach, the problem is reduced to a bidimensional problem, in which only the cross-section of the guide and the surrounding fluid domain are meshed by using finite elements. Then, wedges the top angles of which vary, are studied and the finite element results of the wedge wave speed are compared with experimental results. Finally, the conclusion indicates a way to extend this approach to waveguides of any cross-section.

Fracture Capabilities in Grizzly with the extended Finite Element Method (X-FEM)

Energy Technology Data Exchange (ETDEWEB)

Dolbow, John [Idaho National Lab. (INL), Idaho Falls, ID (United States); Zhang, Ziyu [Idaho National Lab. (INL), Idaho Falls, ID (United States); Spencer, Benjamin [Idaho National Lab. (INL), Idaho Falls, ID (United States); Jiang, Wen [Idaho National Lab. (INL), Idaho Falls, ID (United States)

2015-09-01

Efforts are underway to develop fracture mechanics capabilities in the Grizzly code to enable it to be used to perform deterministic fracture assessments of degraded reactor pressure vessels (RPVs). A capability was previously developed to calculate three-dimensional interaction- integrals to extract mixed-mode stress-intensity factors. This capability requires the use of a finite element mesh that conforms to the crack geometry. The eXtended Finite Element Method (X-FEM) provides a means to represent a crack geometry without explicitly fitting the finite element mesh to it. This is effected by enhancing the element kinematics to represent jump discontinuities at arbitrary locations inside of the element, as well as the incorporation of asymptotic near-tip fields to better capture crack singularities. In this work, use of only the discontinuous enrichment functions was examined to see how accurate stress intensity factors could still be calculated. This report documents the following work to enhance Grizzly’s engineering fracture capabilities by introducing arbitrary jump discontinuities for prescribed crack geometries; X-FEM Mesh Cutting in 3D: to enhance the kinematics of elements that are intersected by arbitrary crack geometries, a mesh cutting algorithm was implemented in Grizzly. The algorithm introduces new virtual nodes and creates partial elements, and then creates a new mesh connectivity; Interaction Integral Modifications: the existing code for evaluating the interaction integral in Grizzly was based on the assumption of a mesh that was fitted to the crack geometry. Modifications were made to allow for the possibility of a crack front that passes arbitrarily through the mesh; and Benchmarking for 3D Fracture: the new capabilities were benchmarked against mixed-mode three-dimensional fracture problems with known analytical solutions.

Finite element formulation for fluid-structure interaction in three-dimensional space

International Nuclear Information System (INIS)

Kulak, R.F.

1979-01-01

A development is presented for a three-dimension hexahedral hydrodynamic finite-element. Using trilinear shape functions and assuming a constant pressure field in each element, simple relations were obtained for internal nodal forces. Because the formulation was based upon a rate approach it was applicable to problems involving large displacements. This element was incorporated into an existing plate-shell finite element code. Diagonal mass matrices were used and the resulting discrete equations of motion were solved using explicit temporal integrator. Results for several problems were presented which compare numerical predictions to closed form analytical solutions. In addition, the fluid-structure interaction problem of a fluid-filled, cylindrical vessel containing internal cylinders was studied. The internal cylinders were cantilever supported from the top cover of the vessel and were periodically located circumferentially at a fixed radius. A pressurized cylindrical cavity located at the bottom of the vessel at its centerline provided the loading

A finite element method for calculating the 3-dimensional magnetic fields of cyclotron

International Nuclear Information System (INIS)

Zhao Xiaofeng

1986-01-01

A series of formula of the finite element method (scalar potential) for calculating the three-dimensional magnetic field of the main magnet of a sector focused cyclotron, and the realization method of the periodic boundary conditions in the code are given

An assessment of finite element synthesis model using 3D FEM

International Nuclear Information System (INIS)

Jagannathan, V.; Rastogi, B.P

1983-01-01

For the solution of multigroup diffusion theory equations a code system called FEMSYN has been developed. The code FEMSYN incorporates a finite difference (FD) module, a finite element (FE) module and one module based on single channel flux synthesis (SCFS) method, which uses trial functions generated by either FD or FE method. The 3D problems can be tackled by all the three modules. The 3D FE solutions are cheaper in comparison with fine mesh FD computations and are quite accurate. The accuracy is improved by proper choice of rectangular and triangular right prismatic elements and polynomial order in case of 3D FE analyses. Synthesis method using FE trial functions provides the fastest means of solution for 3D problems. In this paper, the FE synthesis method has been compared with 3D FE technique for the case of IAEA PWR benchmark problem. It is concluded from the above comparison that FE synthesis method can give at negligible computational efforts, accurate eigenvalue estimates and satisfactory prediction of reactor core power profiles

Finite element analysis of prestressed concrete reactor vessels

International Nuclear Information System (INIS)

Smith, P.D.; Cook, W.A.; Anderson, C.A.

1977-01-01

This paper discusses the development of a finite element code suitable for the safety analysis of prestressed concrete reactor vessels. The project has involved modification of a general purpose computer code to handle reinforced concrete structures as well as comparison of results obtained with the code against published experimental data. The NONSAP nonlinear structural analysis program was selected for the ease with which it can be modified to encompass problems peculiar to nuclear reactors. Pre- and post-processors have been developed for mesh generation and for graphical display of response variables. An out-of-core assembler and solver have been developed for the analysis of large three dimensional problems. The constitutive model for short term loads forms an orthotropic stress-strain relationship in which the concrete and the reinforcing steel are treated as a composite. The variation of stiffness and strength of concrete under multiaxial stress states is accounted for. Cracks are allowed to form at element integration points based on a three dimensional failure envelope in stress space. Composite tensile and shear properties across a crack are modified to account for bond degradation and for dowel action of the reinforcement. The constitutive law for creep is base on the expansion of the usual creep compliance function in the form of a Dirichlet exponential series. Empirical creep data are then fit to the Dirichlet series approximation by means of a least squares procedure. The incremental deformation process is subsequently reduced to a series of variable stiffness elasticity problems in which the past stress history is represented by a finite number of hidden material variables

Finite rotation shells basic equations and finite elements for Reissner kinematics

CERN Document Server

Wisniewski, K

2010-01-01

This book covers theoretical and computational aspects of non-linear shells. Several advanced topics of shell equations and finite elements - not included in standard textbooks on finite elements - are addressed, and the book includes an extensive bibliography.

A multigrid solution method for mixed hybrid finite elements

Energy Technology Data Exchange (ETDEWEB)

Schmid, W. [Universitaet Augsburg (Germany)

1996-12-31

We consider the multigrid solution of linear equations arising within the discretization of elliptic second order boundary value problems of the form by mixed hybrid finite elements. Using the equivalence of mixed hybrid finite elements and non-conforming nodal finite elements, we construct a multigrid scheme for the corresponding non-conforming finite elements, and, by this equivalence, for the mixed hybrid finite elements, following guidelines from Arbogast/Chen. For a rectangular triangulation of the computational domain, this non-conforming schemes are the so-called nodal finite elements. We explicitly construct prolongation and restriction operators for this type of non-conforming finite elements. We discuss the use of plain multigrid and the multilevel-preconditioned cg-method and compare their efficiency in numerical tests.

A benchmark study of 2D and 3D finite element calculations simulating dynamic pulse buckling tests of cylindrical shells under axial impact

International Nuclear Information System (INIS)

Hoffman, E.L.; Ammerman, D.J.

1993-01-01

A series of tests investigating dynamic pulse buckling of a cylindrical shell under axial impact is compared to several finite element simulations of the event. The purpose of the study is to compare the performance of the various analysis codes and element types with respect to a problem which is applicable to radioactive material transport packages, and ultimately to develop a benchmark problem to qualify finite element analysis codes for the transport package design industry

Domain decomposition methods for mortar finite elements

Energy Technology Data Exchange (ETDEWEB)

Widlund, O.

1996-12-31

In the last few years, domain decomposition methods, previously developed and tested for standard finite element methods and elliptic problems, have been extended and modified to work for mortar and other nonconforming finite element methods. A survey will be given of work carried out jointly with Yves Achdou, Mario Casarin, Maksymilian Dryja and Yvon Maday. Results on the p- and h-p-version finite elements will also be discussed.

Why do probabilistic finite element analysis ?

CERN Document Server

Thacker, Ben H

2008-01-01

The intention of this book is to provide an introduction to performing probabilistic finite element analysis. As a short guideline, the objective is to inform the reader of the use, benefits and issues associated with performing probabilistic finite element analysis without excessive theory or mathematical detail.

An efficient finite element solution for gear dynamics

International Nuclear Information System (INIS)

Cooley, C G; Parker, R G; Vijayakar, S M

2010-01-01

A finite element formulation for the dynamic response of gear pairs is proposed. Following an established approach in lumped parameter gear dynamic models, the static solution is used as the excitation in a frequency domain solution of the finite element vibration model. The nonlinear finite element/contact mechanics formulation provides accurate calculation of the static solution and average mesh stiffness that are used in the dynamic simulation. The frequency domain finite element calculation of dynamic response compares well with numerically integrated (time domain) finite element dynamic results and previously published experimental results. Simulation time with the proposed formulation is two orders of magnitude lower than numerically integrated dynamic results. This formulation admits system level dynamic gearbox response, which may include multiple gear meshes, flexible shafts, rolling element bearings, housing structures, and other deformable components.

Modification of the finite element heat and mass transfer code (FEHM) to model multicomponent reactive transport

International Nuclear Information System (INIS)

Viswanathan, H.S.

1996-08-01

The finite element code FEHMN, developed by scientists at Los Alamos National Laboratory (LANL), is a three-dimensional finite element heat and mass transport simulator that can handle complex stratigraphy and nonlinear processes such as vadose zone flow, heat flow and solute transport. Scientists at LANL have been developing hydrologic flow and transport models of the Yucca Mountain site using FEHMN. Previous FEHMN simulations have used an equivalent Kd model to model solute transport. In this thesis, FEHMN is modified making it possible to simulate the transport of a species with a rigorous chemical model. Including the rigorous chemical equations into FEHMN simulations should provide for more representative transport models for highly reactive chemical species. A fully kinetic formulation is chosen for the FEHMN reactive transport model. Several methods are available to computationally implement a fully kinetic formulation. Different numerical algorithms are investigated in order to optimize computational efficiency and memory requirements of the reactive transport model. The best algorithm of those investigated is then incorporated into FEHMN. The algorithm chosen requires for the user to place strongly coupled species into groups which are then solved for simultaneously using FEHMN. The complete reactive transport model is verified over a wide variety of problems and is shown to be working properly. The new chemical capabilities of FEHMN are illustrated by using Los Alamos National Laboratory's site scale model of Yucca Mountain to model two-dimensional, vadose zone 14 C transport. The simulations demonstrate that gas flow and carbonate chemistry can significantly affect 14 C transport at Yucca Mountain. The simulations also prove that the new capabilities of FEHMN can be used to refine and buttress already existing Yucca Mountain radionuclide transport studies

Finite element analysis theory and application with ANSYS

CERN Document Server

Moaveni, Saeed

2015-01-01

For courses in Finite Element Analysis, offered in departments of Mechanical or Civil and Environmental Engineering. While many good textbooks cover the theory of finite element modeling, Finite Element Analysis: Theory and Application with ANSYS is the only text available that incorporates ANSYS as an integral part of its content. Moaveni presents the theory of finite element analysis, explores its application as a design/modeling tool, and explains in detail how to use ANSYS intelligently and effectively. Teaching and Learning Experience This program will provide a better teaching and learning experience-for you and your students. It will help: *Present the Theory of Finite Element Analysis: The presentation of theoretical aspects of finite element analysis is carefully designed not to overwhelm students. *Explain How to Use ANSYS Effectively: ANSYS is incorporated as an integral part of the content throughout the book. *Explore How to Use FEA as a Design/Modeling Tool: Open-ended design problems help stude...

Finite element application to global reactor analysis

International Nuclear Information System (INIS)

Schmidt, F.A.R.

1981-01-01

The Finite Element Method is described as a Coarse Mesh Method with general basis and trial functions. Various consequences concerning programming and application of Finite Element Methods in reactor physics are drawn. One of the conclusions is that the Finite Element Method is a valuable tool in solving global reactor analysis problems. However, problems which can be described by rectangular boxes still can be solved with special coarse mesh programs more efficiently. (orig.) [de

Finite-Element Software for Conceptual Design

DEFF Research Database (Denmark)

Lindemann, J.; Sandberg, G.; Damkilde, Lars

2010-01-01

and research. Forcepad is an effort to provide a conceptual design and teaching tool in a finite-element software package. Forcepad is a two-dimensional finite-element application based on the same conceptual model as image editing applications such as Adobe Photoshop or Microsoft Paint. Instead of using...

Books and monographs on finite element technology

Science.gov (United States)

Noor, A. K.

1985-01-01

The present paper proviees a listing of all of the English books and some of the foreign books on finite element technology, taking into account also a list of the conference proceedings devoted solely to finite elements. The references are divided into categories. Attention is given to fundamentals, mathematical foundations, structural and solid mechanics applications, fluid mechanics applications, other applied science and engineering applications, computer implementation and software systems, computational and modeling aspects, special topics, boundary element methods, proceedings of symmposia and conferences on finite element technology, bibliographies, handbooks, and historical accounts.

Solving nonlinear nonstationary problem of heat-conductivity by finite element method

Directory of Open Access Journals (Sweden)

Антон Янович Карвацький

2016-11-01

Full Text Available Methodology and effective solving algorithm of non-linear dynamic problems of thermal and electric conductivity with significant temperature dependence of thermal and physical properties are given on the basis of finite element method (FEM and Newton linearization method. Discrete equations system FEM was obtained with the use of Galerkin method, where the main function is the finite element form function. The methodology based on successive solving problems of thermal and electrical conductivity has been examined in the work in order to minimize the requirements for calculating resources (RAM. in particular. Having used Mathcad software original programming code was developed to solve the given problem. After investigation of the received results, comparative analyses of accurate solution data and results of numerical solutions, obtained with the use of Matlab programming products, was held. The geometry of one fourth part of the finite sized cylinder was used to test the given numerical model. The discretization of the calculation part was fulfilled using the open programming software for automated Gmsh nets with tetrahedral units, while ParaView, which is an open programming code as well, was used to visualize the calculation results. It was found out that the maximum value violation of potential and temperature determination doesn`t exceed 0,2-0,83% in the given work according to the problem conditions

Probabilistic finite elements for fracture mechanics

Science.gov (United States)

Besterfield, Glen

1988-01-01

The probabilistic finite element method (PFEM) is developed for probabilistic fracture mechanics (PFM). A finite element which has the near crack-tip singular strain embedded in the element is used. Probabilistic distributions, such as expectation, covariance and correlation stress intensity factors, are calculated for random load, random material and random crack length. The method is computationally quite efficient and can be expected to determine the probability of fracture or reliability.

Coupling finite elements and reliability methods - application to safety evaluation of pressurized water reactor vessels

International Nuclear Information System (INIS)

Pitner, P.; Venturini, V.

1995-02-01

When reliability studies are extended form deterministic calculations in mechanics, it is necessary to take into account input parameters variabilities which are linked to the different sources of uncertainty. Integrals must then be calculated to evaluate the failure risk. This can be performed either by simulation methods, or by approximations ones (FORM/SORM). Model in mechanics often require to perform calculation codes. These ones must then be coupled with the reliability calculations. Theses codes can involve large calculation times when they are invoked numerous times during simulations sequences or in complex iterative procedures. Response surface method gives an approximation of the real response from a reduced number of points for which the finite element code is run. Thus, when it is combined with FORM/SORM methods, a coupling can be carried out which gives results in a reasonable calculation time. An application of response surface method to mechanics reliability coupling for a mechanical model which calls for a finite element code is presented. It corresponds to a probabilistic fracture mechanics study of a pressurized water reactor vessel. (authors). 5 refs., 3 figs

Use of heterogeneous finite elements generated by collision probability solutions to calculate a pool reactor core

International Nuclear Information System (INIS)

Calabrese, C.R.; Grant, C.R.

1990-01-01

This work presents comparisons between measured fluxes obtained by activation of Manganese foils in the light water, enriched uranium research pool reactor RA-2 MTR (Materials Testing Reactors) fuel element) and fluxes calculated by the finite element method FEM using DELFIN code, and describes the heterogeneus finite elements by a set of solutions of the transport equations for several different configurations obtained using the collision probability code HUEMUL. The agreement between calculated and measured fluxes is good, and the advantage of using FEM is showed because to obtain the flux distribution with same detail using an usual diffusion calculation it would be necessary 12000 mesh points against the 2000 points that FEM uses, hence the processing time is reduced in a factor ten. An interesting alternative to use in MTR fuel management is presented. (Author) [es

Demonstration of finite element simulations in MOOSE using crystallographic models of irradiation hardening and plastic deformation

Energy Technology Data Exchange (ETDEWEB)

Patra, Anirban [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Wen, Wei [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Martinez Saez, Enrique [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Tome, Carlos [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

2016-05-31

This report describes the implementation of a crystal plasticity framework (VPSC) for irradiation hardening and plastic deformation in the finite element code, MOOSE. Constitutive models for irradiation hardening and the crystal plasticity framework are described in a previous report [1]. Here we describe these models briefly and then describe an algorithm for interfacing VPSC with finite elements. Example applications of tensile deformation of a dog bone specimen and a 3D pre-irradiated bar specimen performed using MOOSE are demonstrated.

Introduction to finite and spectral element methods using Matlab

CERN Document Server

Pozrikidis, Constantine

2014-01-01

The Finite Element Method in One Dimension. Further Applications in One Dimension. High-Order and Spectral Elements in One Dimension. The Finite Element Method in Two Dimensions. Quadratic and Spectral Elements in Two Dimensions. Applications in Mechanics. Viscous Flow. Finite and Spectral Element Methods in Three Dimensions. Appendices. References. Index.

Development of a finite element model for ultrasonic NDT phenomena

International Nuclear Information System (INIS)

Lord, W.

1988-01-01

Ultrasonic NDT techniques are used extensively in the nuclear industry for the detection and characterization of defects in critical structural components such as pressure vessels and piping. The feasibility of applying finite element analysis methods to the problem of modeling ultrasound/defect interactions has been shown. Considerable work remains to be done before a full three-dimensional model is available for the prediction of realistic ultrasonic transducer signals from sound wave interaction with arbitrarily shaped defects in highly attenuative and anisotropic materials. However, a two-dimensional code has been developed that is capable of predicting finite aperture ultrasonic transducer signals associated with wave propagations in isotropic materials and that shows good qualitative agreement with corresponding experimental observations. This 2-D code has now been extended to include anisotropic materials such as centrifugally cast stainless steel (CCSS), a necessary step in the development of the full 3-D code. Results are given showing the capability of the 2-D code to predict the anomalous wave behavior normally associated with ultrasonic wave propagation in anisotropic materials. In addition, a new signal processing technique is discussed, based on the Wigner transformation, that shows promise for application to centrifugally cast stainless steel NDT problems

Complete permutation Gray code implemented by finite state machine

Directory of Open Access Journals (Sweden)

Li Peng

2014-09-01

Full Text Available An enumerating method of complete permutation array is proposed. The list of n! permutations based on Gray code defined over finite symbol set Z(n = {1, 2, …, n} is implemented by finite state machine, named as n-RPGCF. An RPGCF can be used to search permutation code and provide improved lower bounds on the maximum cardinality of a permutation code in some cases.

Comparison between a finite difference model (PUMA) and a finite element model (DELFIN) for simulation of the reactor of the atomic power plant of Atucha I

International Nuclear Information System (INIS)

Grant, C.R.

1996-01-01

The reactor code PUMA, developed in CNEA, simulates nuclear reactors discretizing space in finite difference elements. Core representation is performed by means a cylindrical mesh, but the reactor channels are arranged in an hexagonal lattice. That is why a mapping using volume intersections must be used. This spatial treatment is the reason of an overestimation of the control rod reactivity values, which must be adjusted modifying the incremental cross sections. Also, a not very good treatment of the continuity conditions between core and reflector leads to an overestimation of channel power of the peripherical fuel elements between 5 to 8 per cent. Another code, DELFIN, developed also in CNEA, treats the spatial discretization using heterogeneous finite elements, allowing a correct treatment of the continuity of fluxes and current among elements and a more realistic representation of the hexagonal lattice of the reactor. A comparison between results obtained using both methods in done in this paper. (author). 4 refs., 3 figs

Finite Element Methods and Their Applications

CERN Document Server

Chen, Zhangxin

2005-01-01

This book serves as a text for one- or two-semester courses for upper-level undergraduates and beginning graduate students and as a professional reference for people who want to solve partial differential equations (PDEs) using finite element methods. The author has attempted to introduce every concept in the simplest possible setting and maintain a level of treatment that is as rigorous as possible without being unnecessarily abstract. Quite a lot of attention is given to discontinuous finite elements, characteristic finite elements, and to the applications in fluid and solid mechanics including applications to porous media flow, and applications to semiconductor modeling. An extensive set of exercises and references in each chapter are provided.

Programming the finite element method

CERN Document Server

Smith, I M; Margetts, L

2013-01-01

Many students, engineers, scientists and researchers have benefited from the practical, programming-oriented style of the previous editions of Programming the Finite Element Method, learning how to develop computer programs to solve specific engineering problems using the finite element method. This new fifth edition offers timely revisions that include programs and subroutine libraries fully updated to Fortran 2003, which are freely available online, and provides updated material on advances in parallel computing, thermal stress analysis, plasticity return algorithms, convection boundary c

ABAQUS-EPGEN: a general-purpose finite-element code. Volume 1. User's manual

International Nuclear Information System (INIS)

Hibbitt, H.D.; Karlsson, B.I.; Sorensen, E.P.

1982-10-01

This document is the User's Manual for ABAQUS/EPGEN, a general purpose finite element computer program, designed specifically to serve advanced structural analysis needs. The program contains very general libraries of elements, materials and analysis procedures, and is highly modular, so that complex combinations of features can be put together to model physical problems. The program is aimed at production analysis needs, and for this purpose aspects such as ease-of-use, reliability, flexibility and efficiency have received maximum attention. The input language is designed to make it straightforward to describe complicated models; the analysis procedures are highly automated with the program choosing time or load increments based on user supplied tolerances and controls; and the program offers a wide range of post-processing options for display of the analysis results

Finite element code FENIA verification and application for 3D modelling of thermal state of radioactive waste deep geological repository

Science.gov (United States)

Butov, R. A.; Drobyshevsky, N. I.; Moiseenko, E. V.; Tokarev, U. N.

2017-11-01

The verification of the FENIA finite element code on some problems and an example of its application are presented in the paper. The code is being developing for 3D modelling of thermal, mechanical and hydrodynamical (THM) problems related to the functioning of deep geological repositories. Verification of the code for two analytical problems has been performed. The first one is point heat source with exponential heat decrease, the second one - linear heat source with similar behavior. Analytical solutions have been obtained by the authors. The problems have been chosen because they reflect the processes influencing the thermal state of deep geological repository of radioactive waste. Verification was performed for several meshes with different resolution. Good convergence between analytical and numerical solutions was achieved. The application of the FENIA code is illustrated by 3D modelling of thermal state of a prototypic deep geological repository of radioactive waste. The repository is designed for disposal of radioactive waste in a rock at depth of several hundred meters with no intention of later retrieval. Vitrified radioactive waste is placed in the containers, which are placed in vertical boreholes. The residual decay heat of radioactive waste leads to containers, engineered safety barriers and host rock heating. Maximum temperatures and corresponding times of their establishment have been determined.

Stress categorization in nozzle to pressure vessel connections finite elements models

International Nuclear Information System (INIS)

Albuquerque, Levi Barcelos de

1999-01-01

The ASME Boiler and Pressure Vessel Code, Section III , is the most important code for nuclear pressure vessels design. Its design criteria were developed to preclude the various pressure vessel failure modes throughout the so-called 'Design by Analysis', some of them by imposing stress limits. Thus, failure modes such as plastic collapse, excessive plastic deformation and incremental plastic deformation under cyclic loading (ratchetting) may be avoided by limiting the so-called primary and secondary stresses. At the time 'Design by Analysis' was developed (early 60's) the main tool for pressure vessel design was the shell discontinuity analysis, in which the results were given in membrane and bending stress distributions along shell sections. From that time, the Finite Element Method (FEM) has had a growing use in pressure vessels design. In this case, the stress results are neither normally separated in membrane and bending stress nor classified in primary and secondary stresses. This process of stress separation and classification in Finite Element (FE) results is what is called stress categorization. In order to perform the stress categorization to check results from FE models against the ASME Code stress limits, mainly from 3D solid FE models, several research works have been conducted. This work is included in this effort. First, a description of the ASME Code design criteria is presented. After that, a brief description of how the FEM can be used in pressure vessel design is showed. Several studies found in the literature on stress categorization for pressure vessel FE models are reviewed and commented. Then, the analyses done in this work are presented in which some typical nozzle to pressure vessel connections subjected to internal pressure and concentrated loads were modeled with solid finite elements. The results from linear elastic and limit load analyses are compared to each other and also with the results obtained by formulae for simple shell

Coupled Finite Volume and Finite Element Method Analysis of a Complex Large-Span Roof Structure

Science.gov (United States)

Szafran, J.; Juszczyk, K.; Kamiński, M.

2017-12-01

The main goal of this paper is to present coupled Computational Fluid Dynamics and structural analysis for the precise determination of wind impact on internal forces and deformations of structural elements of a longspan roof structure. The Finite Volume Method (FVM) serves for a solution of the fluid flow problem to model the air flow around the structure, whose results are applied in turn as the boundary tractions in the Finite Element Method problem structural solution for the linear elastostatics with small deformations. The first part is carried out with the use of ANSYS 15.0 computer system, whereas the FEM system Robot supports stress analysis in particular roof members. A comparison of the wind pressure distribution throughout the roof surface shows some differences with respect to that available in the engineering designing codes like Eurocode, which deserves separate further numerical studies. Coupling of these two separate numerical techniques appears to be promising in view of future computational models of stochastic nature in large scale structural systems due to the stochastic perturbation method.

3-D spherical harmonics code FFT3 by the finite Fourier transformation method

International Nuclear Information System (INIS)

Kobayashi, K.

1997-01-01

In the odd order spherical harmonics method, the rigorous boundary condition at the material interfaces is that the even moments of the angular flux and the normal components of the even order moments of current vectors must be continuous. However, it is difficult to derive spatial discretized equations by the finite difference or finite element methods, which satisfy this material interface condition. It is shown that using the finite Fourier transformation method, space discretized equations which satisfy this interface condition can be easily derived. The discrepancies of the flux distribution near void region between spherical harmonics method codes may be due to the difference of application of the material interface condition. (author)

Mixed Element Formulation for the Finite Element-Boundary Integral Method

National Research Council Canada - National Science Library

Meese, J; Kempel, L. C; Schneider, S. W

2006-01-01

A mixed element approach using right hexahedral elements and right prism elements for the finite element-boundary integral method is presented and discussed for the study of planar cavity-backed antennas...

Finite-element formulations for the thermal stress analysis of two- and three-dimensional thin ractor structures

International Nuclear Information System (INIS)

Kulak, R.F.; Kennedy, J.M.; Belytschko, T.B.; Schoeberle, D.F.

1977-01-01

This paper describes finite-element formulations for the thermal stress analysis of LMFBR structures. The first formulation is applicable to large displacement rotation problems in which the strains are small. For this formulation, a general temperature-dependent constituent relationship is derived from a Gibbs potential function and a temperature dependent yield surface. The temperature dependency of the yield surface is based upon a temperature-dependent, material-hardening model. The model uses a temperature-equivalent stress-plastic strain diagram which is generated from isothermal uniaxial stress-strain data. A second formulation is presented for problems characterized by both large displacement-rotations and large strains. Here a set of large strain hypoelastic-plastic relationships are developed to linearly relate the rate of stress to the rate of deformation. The temperature field is described through time-dependent values at mesh node points; the temperature fields in each element are then obtained by interpolation formulas. Hence, problems with both spatial and temporal dependent temperature fields can easily be treated. The above developments were incorporated into two ANL developed finite-element computer codes: the implicit version of STRAW and the 3D Implicit Structural Analysis Code. STRAW is a two-dimensional code with a plane stress/plane strain beam element. The 3D Implicit code has a triangular flat plate element which is capable of sustaining both membrane and bending loads. To insure numerical stability both codes are based on an iterative-incremental solution procedure with equilibrium checks based on an error in energy

Algebraic coding theory over finite commutative rings

CERN Document Server

Dougherty, Steven T

2017-01-01

This book provides a self-contained introduction to algebraic coding theory over finite Frobenius rings. It is the first to offer a comprehensive account on the subject. Coding theory has its origins in the engineering problem of effective electronic communication where the alphabet is generally the binary field. Since its inception, it has grown as a branch of mathematics, and has since been expanded to consider any finite field, and later also Frobenius rings, as its alphabet. This book presents a broad view of the subject as a branch of pure mathematics and relates major results to other fields, including combinatorics, number theory and ring theory. Suitable for graduate students, the book will be of interest to anyone working in the field of coding theory, as well as algebraists and number theorists looking to apply coding theory to their own work.

Development of finite element code for the analysis of coupled thermo-hydro-mechanical behaviors of saturated-unsaturated medium

International Nuclear Information System (INIS)

Ohnishi, Y.; Shibata, H.; Kobayashi, A.

1985-01-01

A model is presented which describes fully coupled thermo-hydro-mechanical behavior of porous geologic medium. The mathematical formulation for the model utilizes the Biot theory for the consolidation and the energy balance equation. The medium is in the condition of saturated-unsaturated flow, then the free surfaces are taken into consideration in the model. The model, incorporated in a finite element numerical procedure, was implemented in a two-dimensional computer code. The code was developed under the assumptions that the medium is poro-elastic and in plane strain condition; water in the ground does not change its phase; heat is transferred by conductive and convective flow. Analytical solutions pertaining to consolidation theory for soils and rocks, thermoelasticity for solids and hydrothermal convection theory provided verification of stress and fluid flow couplings, respectively in the coupled model. Several types of problems are analyzed. The one is a study of some of the effects of completely coupled thermo-hydro-mechanical behavior on the response of a saturated-unsaturated porous rock containing a buried heat source. Excavation of an underground opening which has radioactive wastes at elevated temperatures is modeled and analyzed. The results shows that the coupling phenomena can be estimated at some degree by the numerical procedure. The computer code has a powerful ability to analyze of the repository the complex nature of the repository

FINITE ELEMENT ANALYSIS OF STRUCTURES

Directory of Open Access Journals (Sweden)

PECINGINA OLIMPIA-MIOARA

2015-05-01

Full Text Available The application of finite element method is analytical when solutions can not be applied for deeper study analyzes static, dynamic or other types of requirements in different points of the structures .In practice it is necessary to know the behavior of the structure or certain parts components of the machine under the influence of certain factors static and dynamic . The application of finite element in the optimization of components leads to economic growth , to increase reliability and durability organs studied, thus the machine itself.

Finite Element Simulation of Sheet Metal Forming Process Using Local Interpolation for Tool Surfaces

International Nuclear Information System (INIS)

Hama, Takayuki; Takuda, Hirohiko; Takamura, Masato; Makinouchi, Akitake; Teodosiu, Cristian

2005-01-01

Treatment of contact between a sheet and tools is one of the most difficult problems to deal with in finite-element simulations of sheet forming processes. In order to obtain more accurate tool models without increasing the number of elements, this paper describes a new formulation for contact problems using interpolation proposed by Nagata for tool surfaces. A contact search algorithm between sheet nodes and the interpolated tool surfaces was developed and was introduced into the static-explicit elastoplastic finite-element method code STAMP3D. Simulations of a square cup deep drawing process with a very coarsely discretized punch model were carried out. The simulated results showed that the proposed algorithm gave the proper drawn shape, demonstrating the validity of the proposed algorithm

Surgery simulation using fast finite elements

DEFF Research Database (Denmark)

Bro-Nielsen, Morten

1996-01-01

This paper describes our recent work on real-time surgery simulation using fast finite element models of linear elasticity. In addition, we discuss various improvements in terms of speed and realism......This paper describes our recent work on real-time surgery simulation using fast finite element models of linear elasticity. In addition, we discuss various improvements in terms of speed and realism...

Piezoelectric Analysis of Saw Sensor Using Finite Element Method

Directory of Open Access Journals (Sweden)

Vladimír KUTIŠ

2013-06-01

Full Text Available In this contribution modeling and simulation of surface acoustic waves (SAW sensor using finite element method will be presented. SAW sensor is made from piezoelectric GaN layer and SiC substrate. Two different analysis types are investigated - modal and transient. Both analyses are only 2D. The goal of modal analysis, is to determine the eigenfrequency of SAW, which is used in following transient analysis. In transient analysis, wave propagation in SAW sensor is investigated. Both analyses were performed using FEM code ANSYS.

Using Finite Element Method

Directory of Open Access Journals (Sweden)

M.H.R. Ghoreishy

2008-02-01

Full Text Available This research work is devoted to the footprint analysis of a steel-belted radial tyre (185/65R14 under vertical static load using finite element method. Two models have been developed in which in the first model the tread patterns were replaced by simple ribs while the second model was consisted of details of the tread blocks. Linear elastic and hyper elastic (Arruda-Boyce material models were selected to describe the mechanical behavior of the reinforcing and rubbery parts, respectively. The above two finite element models of the tyre were analyzed under inflation pressure and vertical static loads. The second model (with detailed tread patterns was analyzed with and without friction effect between tread and contact surfaces. In every stage of the analysis, the results were compared with the experimental data to confirm the accuracy and applicability of the model. Results showed that neglecting the tread pattern design not only reduces the computational cost and effort but also the differences between computed deformations do not show significant changes. However, more complicated variables such as shape and area of the footprint zone and contact pressure are affected considerably by the finite element model selected for the tread blocks. In addition, inclusion of friction even in static state changes these variables significantly.

Modelling cohesive laws in finite element simulations via an adapted contact procedure in ABAQUS

DEFF Research Database (Denmark)

Feih, S.

2004-01-01

is not straightforward, and most existing publications consider theoretical and therefore simpler softening shapes. Two possible methods of bridging law approximation areexplained and compared in this report. The bridging laws were implemented in a numerical user subroutine in the finite element code ABAQUS. The main...

Dynamic pulse buckling of cylindrical shells under axial impact: A benchmark study of 2D and 3D finite element calculations

International Nuclear Information System (INIS)

Hoffman, E.L.; Ammerman, D.J.

1995-01-01

A series of tests investigating dynamic pulse buckling of a cylindrical shell under axial impact is compared to several 2D and 3D finite element simulations of the event. The purpose of the work is to investigate the performance of various analysis codes and element types on a problem which is applicable to radioactive material transport packages, and ultimately to develop a benchmark problem to qualify finite element analysis codes for the transport package design industry. During the pulse buckling tests, a buckle formed at each end of the cylinder, and one of the two buckles became unstable and collapsed. Numerical simulations of the test were performed using PRONTO, a Sandia developed transient dynamics analysis code, and ABAQUS/Explicit with both shell and continuum elements. The calculations are compared to the tests with respect to deformed shape and impact load history

Precise magnetostatic field using the finite element method

International Nuclear Information System (INIS)

Nascimento, Francisco Rogerio Teixeira do

2013-01-01

The main objective of this work is to simulate electromagnetic fields using the Finite Element Method. Even in the easiest case of electrostatic and magnetostatic numerical simulation some problems appear when the nodal finite element is used. It is difficult to model vector fields with scalar functions mainly in non-homogeneous materials. With the aim to solve these problems two types of techniques are tried: the adaptive remeshing using nodal elements and the edge finite element that ensure the continuity of tangential components. Some numerical analysis of simple electromagnetic problems with homogeneous and non-homogeneous materials are performed using first, the adaptive remeshing based in various error indicators and second, the numerical solution of waveguides using edge finite element. (author)

ANSYS mechanical APDL for finite element analysis

CERN Document Server

Thompson, Mary Kathryn

2017-01-01

ANSYS Mechanical APDL for Finite Element Analysis provides a hands-on introduction to engineering analysis using one of the most powerful commercial general purposes finite element programs on the market. Students will find a practical and integrated approach that combines finite element theory with best practices for developing, verifying, validating and interpreting the results of finite element models, while engineering professionals will appreciate the deep insight presented on the program's structure and behavior. Additional topics covered include an introduction to commands, input files, batch processing, and other advanced features in ANSYS. The book is written in a lecture/lab style, and each topic is supported by examples, exercises and suggestions for additional readings in the program documentation. Exercises gradually increase in difficulty and complexity, helping readers quickly gain confidence to independently use the program. This provides a solid foundation on which to build, preparing readers...

Finite elements of nonlinear continua

CERN Document Server

Oden, John Tinsley

1972-01-01

Geared toward undergraduate and graduate students, this text extends applications of the finite element method from linear problems in elastic structures to a broad class of practical, nonlinear problems in continuum mechanics. It treats both theory and applications from a general and unifying point of view.The text reviews the thermomechanical principles of continuous media and the properties of the finite element method, and then brings them together to produce discrete physical models of nonlinear continua. The mathematical properties of these models are analyzed, along with the numerical s

A finite element method for flow problems in blast loading

International Nuclear Information System (INIS)

Forestier, A.; Lepareux, M.

1984-06-01

This paper presents a numerical method which describes fast dynamic problems in flow transient situations as in nuclear plants. A finite element formulation has been chosen; it is described by a preprocessor in CASTEM system: GIBI code. For these typical flow problems, an A.L.E. formulation for physical equations is used. So, some applications are presented: the well known problem of shock tube, the same one in 2D case and a last application to hydrogen detonation

Finite element simulation of moisture movement and solute transport in a large caisson

International Nuclear Information System (INIS)

Huyakorn, P.S.; Jones, B.G.; Parker, J.C.; Wadsworth, T.D.; White, H.O. Jr.

1987-01-01

The results of the solute transport experiments performed on compacted, crushed Bandelier Tuff in caisson B of the experimental cluster described by DePoorter (1981) are simulated. Both one- and three-dimensional simulations of solute transport have been performed using two selected finite element codes. Results of bromide and iodide tracer experiments conducted during near-steady flow conditions have been analyzed for pulse additions made on December 6, 1984, and followed over a period of up to 60 days. In addition, a pulse addition of nonconservative strontium tracer on September 28, 1984, during questionably steady flow conditions has been analyzed over a period of 240 days. One-dimensional finite element flow and transport simulations were carried out assuming the porous medium to be homogeneous and the injection source uniformly distributed. To evaluate effects of the nonuniform source distribution and also to investigate effects of inhomogeneous porous medium properties, three dimensional finite element analyses of transport were carried out. Implications of the three-dimensional effects for the design and analysis of future tracer studies are discussed

A comparison study on the performance of lower order solid finite element for elastic analysis of plate and shell structures

International Nuclear Information System (INIS)

Lee, Young Jung; Lee, Sang Jin; Choun, Young Sun; Seo, Jeong Moon

2003-05-01

The objective of this research is to assess the performance of lower order solid finite elements which will be ultimately applied into the safety analysis of nuclear containment building. For the safety analysis of large structures such as nuclear containment building, efficient lower order finite element is necessarily required to calculate the structural response of containment building with low computational cost. In this study, the state of the art formulations of lower order solid finite element are throughly reviewed and the best possible solid finite element is adopted into the development of nuclear containment analysis system. Three 8-node solid finite elements based on standard strain-displacement relationship, B-bar method and EAS method are implemented as computer modules and completely tested with various plate and shell structures. The present results can be directly applied into the analysis code development for general reinforced concrete structures

The finite element method in engineering, 2nd edition

International Nuclear Information System (INIS)

Rao, S.S.

1986-01-01

This work provides a systematic introduction to the various aspects of the finite element method as applied to engineering problems. Contents include: introduction to finite element method; solution of finite element equations; solid and structural mechanics; static analysis; dynamic analysis; heat transfer; fluid mechanics and additional applications

A study on the nonlinear finite element analysis of reinforced concrete structures: shell finite element formulation

Energy Technology Data Exchange (ETDEWEB)

Lee, Sang Jin; Seo, Jeong Moon

2000-08-01

The main goal of this research is to establish a methodology of finite element analysis of containment building predicting not only global behaviour but also local failure mode. In this report, we summerize some existing numerical analysis techniques to be improved for containment building. In other words, a complete description of the standard degenerated shell finite element formulation is provided for nonlinear stress analysis of nuclear containment structure. A shell finite element is derived using the degenerated solid concept which does not rely on a specific shell theory. Reissner-Mindlin assumptions are adopted to consider the transverse shear deformation effect. In order to minimize the sensitivity of the constitutive equation to structural types, microscopic material model is adopted. The four solution algorithms based on the standard Newton-Raphson method are discussed. Finally, two numerical examples are carried out to test the performance of the adopted shell medel.

A study on the nonlinear finite element analysis of reinforced concrete structures: shell finite element formulation

International Nuclear Information System (INIS)

Lee, Sang Jin; Seo, Jeong Moon

2000-08-01

The main goal of this research is to establish a methodology of finite element analysis of containment building predicting not only global behaviour but also local failure mode. In this report, we summerize some existing numerical analysis techniques to be improved for containment building. In other words, a complete description of the standard degenerated shell finite element formulation is provided for nonlinear stress analysis of nuclear containment structure. A shell finite element is derived using the degenerated solid concept which does not rely on a specific shell theory. Reissner-Mindlin assumptions are adopted to consider the transverse shear deformation effect. In order to minimize the sensitivity of the constitutive equation to structural types, microscopic material model is adopted. The four solution algorithms based on the standard Newton-Raphson method are discussed. Finally, two numerical examples are carried out to test the performance of the adopted shell medel

Adaptive Multilevel Methods with Local Smoothing for $H^1$- and $H^{\\mathrm{curl}}$-Conforming High Order Finite Element Methods

KAUST Repository

Janssen, Bä rbel; Kanschat, Guido

2011-01-01

A multilevel method on adaptive meshes with hanging nodes is presented, and the additional matrices appearing in the implementation are derived. Smoothers of overlapping Schwarz type are discussed; smoothing is restricted to the interior of the subdomains refined to the current level; thus it has optimal computational complexity. When applied to conforming finite element discretizations of elliptic problems and Maxwell equations, the method's convergence rates are very close to those for the nonadaptive version. Furthermore, the smoothers remain efficient for high order finite elements. We discuss the implementation in a general finite element code using the example of the deal.II library. © 2011 Societ y for Industrial and Applied Mathematics.

Finite Element Analysis of Pipe T-Joint

OpenAIRE

P.M.Gedkar; Dr. D.V. Bhope

2012-01-01

This paper reports stress analysis of two pressurized cylindrical intersection using finite element method. The different combinations of dimensions of run pipe and the branch pipe are used to investigate thestresses in pipe at the intersection. In this study the stress analysis is accomplished by finite element package ANSYS.

Higher Order Lagrange Finite Elements In M3D

International Nuclear Information System (INIS)

Chen, J.; Strauss, H.R.; Jardin, S.C.; Park, W.; Sugiyama, L.E.; Fu, G.; Breslau, J.

2004-01-01

The M3D code has been using linear finite elements to represent multilevel MHD on 2-D poloidal planes. Triangular higher order elements, up to third order, are constructed here in order to provide M3D the capability to solve highly anisotropic transport problems. It is found that higher order elements are essential to resolve the thin transition layer characteristic of the anisotropic transport equation, particularly when the strong anisotropic direction is not aligned with one of the Cartesian coordinates. The transition layer is measured by the profile width, which is zero for infinite anisotropy. It is shown that only higher order schemes have the ability to make this layer converge towards zero when the anisotropy gets stronger and stronger. Two cases are considered. One has the strong transport direction partially aligned with one of the element edges, the other doesn't have any alignment. Both cases have the strong transport direction misaligned with the grid line by some angles

Introduction to finite element analysis using MATLAB and Abaqus

CERN Document Server

Khennane, Amar

2013-01-01

There are some books that target the theory of the finite element, while others focus on the programming side of things. Introduction to Finite Element Analysis Using MATLAB(R) and Abaqus accomplishes both. This book teaches the first principles of the finite element method. It presents the theory of the finite element method while maintaining a balance between its mathematical formulation, programming implementation, and application using commercial software. The computer implementation is carried out using MATLAB, while the practical applications are carried out in both MATLAB and Abaqus. MA

Hydro-mechanical coupled simulation of hydraulic fracturing using the eXtended Finite Element Method (XFEM)

Science.gov (United States)

Youn, Dong Joon

combining XFEM and random field is named as eXtended Random Finite Element Method (XRFEM). All the numerical analysis codes in this thesis are written in Fortran 2003, and these codes are applicable as a series of sub-modules within a suite of finite element codes developed by Smith and Griffiths (2004).

Nonlinear finite element modeling of corrugated board

Science.gov (United States)

A. C. Gilchrist; J. C. Suhling; T. J. Urbanik

1999-01-01

In this research, an investigation on the mechanical behavior of corrugated board has been performed using finite element analysis. Numerical finite element models for corrugated board geometries have been created and executed. Both geometric (large deformation) and material nonlinearities were included in the models. The analyses were performed using the commercial...

Bibliography for finite elements. [2200 references

Energy Technology Data Exchange (ETDEWEB)

Whiteman, J R [comp.

1975-01-01

This bibliography cites almost all of the significant papers on advances in the mathematical theory of finite elements. Reported are applications in aeronautical, civil, mechanical, nautical and nuclear engineering. Such topics as classical analysis, functional analysis, approximation theory, fluids, and diffusion are covered. Over 2200 references to publications up to the end of 1974 are included. Publications are listed alphabetically by author and also by keywords. In addition, finite element packages are listed.

Complex finite element sensitivity method for creep analysis

International Nuclear Information System (INIS)

Gomez-Farias, Armando; Montoya, Arturo; Millwater, Harry

2015-01-01

The complex finite element method (ZFEM) has been extended to perform sensitivity analysis for mechanical and structural systems undergoing creep deformation. ZFEM uses a complex finite element formulation to provide shape, material, and loading derivatives of the system response, providing an insight into the essential factors which control the behavior of the system as a function of time. A complex variable-based quadrilateral user element (UEL) subroutine implementing the power law creep constitutive formulation was incorporated within the Abaqus commercial finite element software. The results of the complex finite element computations were verified by comparing them to the reference solution for the steady-state creep problem of a thick-walled cylinder in the power law creep range. A practical application of the ZFEM implementation to creep deformation analysis is the calculation of the skeletal point of a notched bar test from a single ZFEM run. In contrast, the standard finite element procedure requires multiple runs. The value of the skeletal point is that it identifies the location where the stress state is accurate, regardless of the certainty of the creep material properties. - Highlights: • A novel finite element sensitivity method (ZFEM) for creep was introduced. • ZFEM has the capability to calculate accurate partial derivatives. • ZFEM can be used for identification of the skeletal point of creep structures. • ZFEM can be easily implemented in a commercial software, e.g. Abaqus. • ZFEM results were shown to be in excellent agreement with analytical solutions

Summary of the models and methods for the FEHM application - a finite-element heat- and mass-transfer code

International Nuclear Information System (INIS)

Zyvoloski, G.A.; Robinson, B.A.; Dash, Z.V.; Trease, L.L.

1997-07-01

The mathematical models and numerical methods employed by the FEHM application, a finite-element heat- and mass-transfer computer code that can simulate nonisothermal multiphase multi-component flow in porous media, are described. The use of this code is applicable to natural-state studies of geothermal systems and groundwater flow. A primary use of the FEHM application will be to assist in the understanding of flow fields and mass transport in the saturated and unsaturated zones below the proposed Yucca Mountain nuclear waste repository in Nevada. The component models of FEHM are discussed. The first major component, Flow- and Energy-Transport Equations, deals with heat conduction; heat and mass transfer with pressure- and temperature-dependent properties, relative permeabilities and capillary pressures; isothermal air-water transport; and heat and mass transfer with noncondensible gas. The second component, Dual-Porosity and Double-Porosity/Double-Permeability Formulation, is designed for problems dominated by fracture flow. Another component, The Solute-Transport Models, includes both a reactive-transport model that simulates transport of multiple solutes with chemical reaction and a particle-tracking model. Finally, the component, Constitutive Relationships, deals with pressure- and temperature-dependent fluid/air/gas properties, relative permeabilities and capillary pressures, stress dependencies, and reactive and sorbing solutes. Each of these components is discussed in detail, including purpose, assumptions and limitations, derivation, applications, numerical method type, derivation of numerical model, location in the FEHM code flow, numerical stability and accuracy, and alternative approaches to modeling the component

Finite-element formulations for the thermal stress analysis of two- and three-dimensional thin reactor structures

International Nuclear Information System (INIS)

Kulak, R.F.; Kennedy, J.M.; Belytschko, T.B.; Schoeberle, D.F.

1977-01-01

This paper describes finite-element formulations for the thermal stress analysis of LMFBR structures. The first formulation is applicable to large displacement rotation problems in which the strains are small. For this formulation, a general temperature-dependent constituent relationship is derived from a Gibbs potential and a temperature dependent surface. A second formulation is presented for problems characterized by both large displacement-rotations and large strains. Here a set of large strain hypoelastic-plastic relationships are developed to linearly relate the rate of stress to the rate of deformation. These developments were incorporated into two ANL developed finite-element computer codes: the implicit version of STRAW and the 3D Implicit Structural Analaysis code. A set of problems is presented to validate both the 3D and 2D programs and to illustrate their applicability to a variety of problems. (Auth.)

A suitable low-order, eight-node tetrahedral finite element for solids

Energy Technology Data Exchange (ETDEWEB)

Key, S.W.; Heinstein, M.S.; Stone, C.M.; Mello, F.J.; Blanford, M.L.; Budge, K.G.

1998-03-01

To use the all-tetrahedral mesh generation existing today, the authors have explored the creation of a computationally efficient eight-node tetrahedral finite element (a four-node tetrahedral finite element enriched with four mid-face nodal points). The derivation of the element`s gradient operator, studies in obtaining a suitable mass lumping, and the element`s performance in applications are presented. In particular they examine the eight-node tetrahedral finite element`s behavior in longitudinal plane wave propagation, in transverse cylindrical wave propagation, and in simulating Taylor bar impacts. The element samples only constant strain states and, therefore, has 12 hour-glass modes. In this regard it bears similarities to the eight-node, mean-quadrature hexahedral finite element. Comparisons with the results obtained from the mean-quadrature eight-node hexahedral finite element and the four-node tetrahedral finite element are included. Given automatic all-tetrahedral meshing, the eight-node, constant-strain tetrahedral finite element is a suitable replacement for the eight-node hexahedral finite element in those cases where mesh generation requires an inordinate amount of user intervention and direction to obtain acceptable mesh properties.

Finite element investigation of explosively formed projectiles (EFP)

International Nuclear Information System (INIS)

Ahmad, I.

1999-01-01

This thesis report represents the numerical simulation of explosively formed projectiles (EFP), a type of linear self-forging fragment device. The simulation is performed using a finite element code DYNA2D. It also explicates that how the shape, velocity and kinetic energy of an explosively formed projectile is effected by various parameters. Different parameters investigated are mesh density, material, thickness, contour and types of liner. Effect of shape of casing and material model is also analyzed. The shapes of projectiles at different times after detonation are shown. The maximum velocity and kinetic energy of the projectile have been used to ascertain the effect of above mentioned parameters. (author)

Finite element modeling of pellet-clad mechanical interaction with ABAQUS

International Nuclear Information System (INIS)

Cheon, C. S.; Lee, B. H.; Koo, Y. H.; Oh, J. Y.; Son, D. S.

2002-01-01

Pellet-clad mechanical interaction (PCMI) was modelled by an axisymmetric finite element method. Thermomechanical models of pellet and clad materials and a contact model for their interaction have been implemented in addition to the application of appropriate boundary conditions so that the FE model was configured. Temperature and displacement were evaluated through a coupled analysis using a general purposed FE code, ABAQUS. Also, a batch program has been developed to efficiently deal with a series of jobs such as making an interface with a fuel performance code, the generation of an input deck for ABAQUS code and its execution, and an interpretation of the output. Under various conditions, results from the present FE model were analyzed. Preliminary verification was conducted by comparing the clad elongation measured during an in-pile PCMI experiment with that calculated by means of the developed FE model

Application of Mass Lumped Higher Order Finite Elements

International Nuclear Information System (INIS)

J. Chen, H.R. Strauss, S.C. Jardin, W. Park, L.E. Sugiyama, G. Fu, J. Breslau

2005-01-01

There are many interesting phenomena in extended-MHD such as anisotropic transport, mhd, 2-fluid effects stellarator and hot particles. Any one of them challenges numerical analysts, and researchers are seeking for higher order methods, such as higher order finite difference, higher order finite elements and hp/spectral elements. It is true that these methods give more accurate solution than their linear counterparts. However, numerically they are prohibitively expensive. Here we give a successful solution of this conflict by applying mass lumped higher order finite elements. This type of elements not only keep second/third order accuracy but also scale closely to linear elements by doing mass lumping. This is especially true for second order lump elements. Full M3D and anisotropic transport models are studied

The finite element method its basis and fundamentals

CERN Document Server

Zienkiewicz, Olek C; Zhu, JZ

2013-01-01

The Finite Element Method: Its Basis and Fundamentals offers a complete introduction to the basis of the finite element method, covering fundamental theory and worked examples in the detail required for readers to apply the knowledge to their own engineering problems and understand more advanced applications. This edition sees a significant rearrangement of the book's content to enable clearer development of the finite element method, with major new chapters and sections added to cover: Weak forms Variational forms Multi-dimensional field prob

A suitable low-order, eight-node tetrahedral finite element for solids

International Nuclear Information System (INIS)

Key, S.W.; Heinstein, M.S.; Stone, C.M.; Mello, F.J.; Blanford, M.L.; Budge, K.G.

1998-03-01

To use the all-tetrahedral mesh generation existing today, the authors have explored the creation of a computationally efficient eight-node tetrahedral finite element (a four-node tetrahedral finite element enriched with four mid-face nodal points). The derivation of the element's gradient operator, studies in obtaining a suitable mass lumping, and the element's performance in applications are presented. In particular they examine the eight-node tetrahedral finite element's behavior in longitudinal plane wave propagation, in transverse cylindrical wave propagation, and in simulating Taylor bar impacts. The element samples only constant strain states and, therefore, has 12 hour-glass modes. In this regard it bears similarities to the eight-node, mean-quadrature hexahedral finite element. Comparisons with the results obtained from the mean-quadrature eight-node hexahedral finite element and the four-node tetrahedral finite element are included. Given automatic all-tetrahedral meshing, the eight-node, constant-strain tetrahedral finite element is a suitable replacement for the eight-node hexahedral finite element in those cases where mesh generation requires an inordinate amount of user intervention and direction to obtain acceptable mesh properties

Finite element formulation for a digital image correlation method

International Nuclear Information System (INIS)

Sun Yaofeng; Pang, John H. L.; Wong, Chee Khuen; Su Fei

2005-01-01

A finite element formulation for a digital image correlation method is presented that will determine directly the complete, two-dimensional displacement field during the image correlation process on digital images. The entire interested image area is discretized into finite elements that are involved in the common image correlation process by use of our algorithms. This image correlation method with finite element formulation has an advantage over subset-based image correlation methods because it satisfies the requirements of displacement continuity and derivative continuity among elements on images. Numerical studies and a real experiment are used to verify the proposed formulation. Results have shown that the image correlation with the finite element formulation is computationally efficient, accurate, and robust

Application of finite-element-methods in food processing

DEFF Research Database (Denmark)

Risum, Jørgen

2004-01-01

Presentation of the possible use of finite-element-methods in food processing. Examples from diffusion studies are given.......Presentation of the possible use of finite-element-methods in food processing. Examples from diffusion studies are given....

Accuracy of finite-element models for the stress analysis of multiple-holed moderator blocks

International Nuclear Information System (INIS)

Smith, P.D.; Sullivan, R.M.; Lewis, A.C.; Yu, H.J.

1981-01-01

Two steps have been taken to quantify and improve the accuracy in the analysis. First, the limitations of various approximation techniques have been studied with the aid of smaller benchmark problems containing fewer holes. Second, a new family of computer programs has been developed for handling such large problems. This paper describes the accuracy studies and the benchmark problems. A review is given of some proposed modeling techniques including local mesh refinement, homogenization, a special-purpose finite element, and substructuring. Some limitations of these approaches are discussed. The new finite element programs and the features that contribute to their efficiency are discussed. These include a standard architecture for out-of-core data processing and an equation solver that operates on a peripheral array processor. The central conclusions of the paper are: (1) modeling approximation methods such as local mesh refinement and homogenization tend to be unreliable, and they should be justified by a fine mesh benchmark analysis; and (2) finite element codes are now available that can achieve accurate solutions at a reasonable cost, and there is no longer a need to employ modeling approximations in the two-dimensional analysis of HTGR fuel elements. 10 figures

Adaptive Multilevel Methods with Local Smoothing for $H^1$- and $H^{\\mathrm{curl}}$-Conforming High Order Finite Element Methods

KAUST Repository

Janssen, Bärbel

2011-01-01

A multilevel method on adaptive meshes with hanging nodes is presented, and the additional matrices appearing in the implementation are derived. Smoothers of overlapping Schwarz type are discussed; smoothing is restricted to the interior of the subdomains refined to the current level; thus it has optimal computational complexity. When applied to conforming finite element discretizations of elliptic problems and Maxwell equations, the method\\'s convergence rates are very close to those for the nonadaptive version. Furthermore, the smoothers remain efficient for high order finite elements. We discuss the implementation in a general finite element code using the example of the deal.II library. © 2011 Societ y for Industrial and Applied Mathematics.

Two-dimensional calculation by finite element method of velocity field and temperature field development in fast reactor fuel assembly. II

International Nuclear Information System (INIS)

Schmid, J.

1985-11-01

A package of updated computer codes for velocity and temperature field calculations for a fast reactor fuel subassembly (or its part) by the finite element method is described. Isoparametric triangular elements of the second degree are used. (author)

Development of an anisotropic beam finite element for composite wind turbine blades in multibody system

DEFF Research Database (Denmark)

Kim, Taeseong; Hansen, Anders Melchior; Branner, Kim

2013-01-01

In this paper a new anisotropic beam finite element for composite wind turbine blades is developed and implemented into the aeroelastic nonlinear multibody code, HAWC2, intended to be used to investigate if use of anisotropic material layups in wind turbine blades can be tailored for improved...

Advances in 3D electromagnetic finite element modeling

International Nuclear Information System (INIS)

Nelson, E.M.

1997-01-01

Numerous advances in electromagnetic finite element analysis (FEA) have been made in recent years. The maturity of frequency domain and eigenmode calculations, and the growth of time domain applications is briefly reviewed. A high accuracy 3D electromagnetic finite element field solver employing quadratic hexahedral elements and quadratic mixed-order one-form basis functions will also be described. The solver is based on an object-oriented C++ class library. Test cases demonstrate that frequency errors less than 10 ppm can be achieved using modest workstations, and that the solutions have no contamination from spurious modes. The role of differential geometry and geometrical physics in finite element analysis is also discussed

On symmetric pyramidal finite elements

Czech Academy of Sciences Publication Activity Database

Liu, L.; Davies, K. B.; Yuan, K.; Křížek, Michal

2004-01-01

Roč. 11, 1-2 (2004), s. 213-227 ISSN 1492-8760 R&D Projects: GA AV ČR IAA1019201 Institutional research plan: CEZ:AV0Z1019905 Keywords : mesh generation * finite element method * composite elements Subject RIV: BA - General Mathematics Impact factor: 0.108, year: 2004

Comparison of finite-element stress analysis with experimental copper sphere impacts

International Nuclear Information System (INIS)

Frantz, C.E.; Hecker, S.S.; Stout, M.G.; Browning, R.V.

1980-07-01

Three copper spheres were impacted on targets of varying surface finishes at 100 m/s. Impact face friction was varied for each test and the impact was photographed with a high-speed camera. Postimpact strains and deformation were measured. A finite-element computer code, NONSAP, was used to model the impact. The best agreement between computer prediction and experiment was obtained using isoparametric elements, a graded mesh, and actual high-strain-rate copper stress-strain data. Frictional conditions at the impact face were also modeled by altering the standard NONSAP code. The most critical test of NONSAP was accurate prediction of experimental impact strains. The best agreement we could obtain had a maximum point-to-point error of 20%, although in general, the comparison was much better. Results of this research indicate that we must know more about material and impact interface friction in order to obtain reliable numerical predictions

Finite-element analysis of dynamic fracture

Science.gov (United States)

Aberson, J. A.; Anderson, J. M.; King, W. W.

1976-01-01

Applications of the finite element method to the two dimensional elastodynamics of cracked structures are presented. Stress intensity factors are computed for two problems involving stationary cracks. The first serves as a vehicle for discussing lumped-mass and consistent-mass characterizations of inertia. In the second problem, the behavior of a photoelastic dynamic tear test specimen is determined for the time prior to crack propagation. Some results of a finite element simulation of rapid crack propagation in an infinite body are discussed.

The finite element method in electromagnetics

CERN Document Server

Jin, Jianming

2014-01-01

A new edition of the leading textbook on the finite element method, incorporating major advancements and further applications in the field of electromagnetics The finite element method (FEM) is a powerful simulation technique used to solve boundary-value problems in a variety of engineering circumstances. It has been widely used for analysis of electromagnetic fields in antennas, radar scattering, RF and microwave engineering, high-speed/high-frequency circuits, wireless communication, electromagnetic compatibility, photonics, remote sensing, biomedical engineering, and space exploration. The

Finite elements for analysis and design

CERN Document Server

Akin, J E; Davenport, J H

1994-01-01

The finite element method (FEM) is an analysis tool for problem-solving used throughout applied mathematics, engineering, and scientific computing. Finite Elements for Analysis and Design provides a thoroughlyrevised and up-to-date account of this important tool and its numerous applications, with added emphasis on basic theory. Numerous worked examples are included to illustrate the material.Key Features* Akin clearly explains the FEM, a numerical analysis tool for problem-solving throughout applied mathematics, engineering and scientific computing* Basic theory has bee

Finite element implementation and numerical issues of strain gradient plasticity with application to metal matrix composites

DEFF Research Database (Denmark)

Frederiksson, Per; Gudmundson, Peter; Mikkelsen, Lars Pilgaard

2009-01-01

A framework of finite element equations for strain gradient plasticity is presented. The theoretical framework requires plastic strain degrees of freedom in addition to displacements and a plane strain version is implemented into a commercial finite element code. A couple of different elements...... of quadrilateral type are examined and a few numerical issues are addressed related to these elements as well as to strain gradient plasticity theories in general. Numerical results are presented for an idealized cell model of a metal matrix composite under shear loading. It is shown that strengthening due...... to fiber size is captured but strengthening due to fiber shape is not. A few modelling aspects of this problem are discussed as well. An analytic solution is also presented which illustrates similarities to other theories....

Finite Element Method in Machining Processes

CERN Document Server

Markopoulos, Angelos P

2013-01-01

Finite Element Method in Machining Processes provides a concise study on the way the Finite Element Method (FEM) is used in the case of manufacturing processes, primarily in machining. The basics of this kind of modeling are detailed to create a reference that will provide guidelines for those who start to study this method now, but also for scientists already involved in FEM and want to expand their research. A discussion on FEM, formulations and techniques currently in use is followed up by machining case studies. Orthogonal cutting, oblique cutting, 3D simulations for turning and milling, grinding, and state-of-the-art topics such as high speed machining and micromachining are explained with relevant examples. This is all supported by a literature review and a reference list for further study. As FEM is a key method for researchers in the manufacturing and especially in the machining sector, Finite Element Method in Machining Processes is a key reference for students studying manufacturing processes but al...

Inside finite elements

CERN Document Server

Weiser, Martin

2016-01-01

All relevant implementation aspects of finite element methods are discussed in this book. The focus is on algorithms and data structures as well as on their concrete implementation. Theory is covered as far as it gives insight into the construction of algorithms. Throughout the exercises a complete FE-solver for scalar 2D problems will be implemented in Matlab/Octave.

Multi Scale Finite Element Analyses By Using SEM-EBSD Crystallographic Modeling and Parallel Computing

International Nuclear Information System (INIS)

Nakamachi, Eiji

2005-01-01

A crystallographic homogenization procedure is introduced to the conventional static-explicit and dynamic-explicit finite element formulation to develop a multi scale - double scale - analysis code to predict the plastic strain induced texture evolution, yield loci and formability of sheet metal. The double-scale structure consists of a crystal aggregation - micro-structure - and a macroscopic elastic plastic continuum. At first, we measure crystal morphologies by using SEM-EBSD apparatus, and define a unit cell of micro structure, which satisfy the periodicity condition in the real scale of polycrystal. Next, this crystallographic homogenization FE code is applied to 3N pure-iron and 'Benchmark' aluminum A6022 polycrystal sheets. It reveals that the initial crystal orientation distribution - the texture - affects very much to a plastic strain induced texture and anisotropic hardening evolutions and sheet deformation. Since, the multi-scale finite element analysis requires a large computation time, a parallel computing technique by using PC cluster is developed for a quick calculation. In this parallelization scheme, a dynamic workload balancing technique is introduced for quick and efficient calculations

On the application of finite element method in the solution of steady state diffusion equation

International Nuclear Information System (INIS)

Ono, S.

1982-01-01

The solution of the steady state neutron diffusion equation is obtained by using the finite element method. Specifically the variational approach is used for one dimensional problems and the weighted residual method (Galerkin) for one and two dimensional problems. The spatial domain is divided into retangular elements and the neutron flux is approximated by linear (one dimensional case), and bilinear (two-dimensional case) functions. Numerical results are obtained with a FORTRAN IV computer program and compared with those obtained by the finite difference CITATION code. The results show that linear or bilinear functions, do not satisfactorily describe the differential parameters in highly heterogeneous reactor cases, but provide good results for integral parameters such as multiplication factor. (Author) [pt

A finite-element for the analysis of shell intersections

International Nuclear Information System (INIS)

Koves, W.J.; Nair, S.

1994-01-01

The analysis of discontinuity stresses at shell intersections is a problem of great importance in several major industries. Some of the major areas of interest are pressure-containing equipment, such as reactors and piping, in the nuclear and fossil power industry; pressure vessels and heat exchangers in the petrochemical industry; bracing in offshore oil platforms; and aerospace structures. A specialized shell-intersection finite element, which is compatible with adjoining shell elements, has been developed that has the capability of physically representing the complex three-dimensional geometry and stress state at shell intersections. The element geometry is a contoured shape that matches a wide variety of practical nozzle configurations used in ASME Code pressure vessel construction, and allows computational rigor. A closed-form theory of elasticity solution was used to compute the stress state and strain energy in the element. The concept of an energy-equivalent nodal displacement and force vector set was then developed to allow complete compatibility with adjoining shell elements and retain the analytical rigor within the element. This methodology provides a powerful and robust computation scheme that maintains the computational efficiency of shell element solutions. The shell-intersection element was then applied to the cylinder-sphere and cylinder-cylinder intersection problems

A geometric toolbox for tetrahedral finite element partitions

NARCIS (Netherlands)

Brandts, J.; Korotov, S.; Křížek, M.; Axelsson, O.; Karátson, J.

2011-01-01

In this work we present a survey of some geometric results on tetrahedral partitions and their refinements in a unified manner. They can be used for mesh generation and adaptivity in practical calculations by the finite element method (FEM), and also in theoretical finite element (FE) analysis.

Finite element analysis of rotating beams physics based interpolation

CERN Document Server

Ganguli, Ranjan

2017-01-01

This book addresses the solution of rotating beam free-vibration problems using the finite element method. It provides an introduction to the governing equation of a rotating beam, before outlining the solution procedures using Rayleigh-Ritz, Galerkin and finite element methods. The possibility of improving the convergence of finite element methods through a judicious selection of interpolation functions, which are closer to the problem physics, is also addressed. The book offers a valuable guide for students and researchers working on rotating beam problems – important engineering structures used in helicopter rotors, wind turbines, gas turbines, steam turbines and propellers – and their applications. It can also be used as a textbook for specialized graduate and professional courses on advanced applications of finite element analysis.

Numerical study of effects of the beam tube on laser fields with a three-dimensional simulation code using the finite element method

CERN Document Server

Sobajima, M; Yamazaki, T; Yoshikawa, K; Ohnishi, M; Toku, H; Masuda, K; Kitagaki, J; Nakamura, T

1999-01-01

In January 1997, the Beijing FEL observed large laser amplification at 8-18 mu m. However, through the collaborative work, it was found from both experiments and numerical simulations that the laser loss on the beam tube wall was not negligible, and that the saturation was not seen in the relatively long wavelength range because of this loss. This calls for further investigation on the effects of the beam tube of finite size. In order to include such effects self-consistently, we have developed a new three-dimensional code that can solve equations with the boundary conditions of the beam tube by using the Finite Element Method. Results show that the beam tube effects are dominant in deriving higher laser modes in the tube, compared with the optical guiding effects, and consequently reduced gain especially in the longer wavelength range, where the beam tube effects are greatly emphasized. It is also found that TEM sub 0 sub 2 mode is the most dominant higher mode in the beam tube, and is also the main cause of...

Impact of new computing systems on finite element computations

International Nuclear Information System (INIS)

Noor, A.K.; Fulton, R.E.; Storaasi, O.O.

1983-01-01

Recent advances in computer technology that are likely to impact finite element computations are reviewed. The characteristics of supersystems, highly parallel systems, and small systems (mini and microcomputers) are summarized. The interrelations of numerical algorithms and software with parallel architectures are discussed. A scenario is presented for future hardware/software environment and finite element systems. A number of research areas which have high potential for improving the effectiveness of finite element analysis in the new environment are identified

GOMESH, Finite Elements Structure Plot with Triangular Mesh

International Nuclear Information System (INIS)

Draper, J.

1977-01-01

1 - Nature of the physical problem solved: Graphical representation of calculations on structures with finite subdivision. 2 - Method of solution: GOMESH treats meshes with triangular basic elements. The program uses the same punched cards as those required for the input to the 'STAG' series of stress analysis codes and can prepare three basic mesh diagrams which differ in their mode of numbering. One objective of using these diagrams is to show up errors in the card deck by making them visually recognisable. Furthermore, digital tests are made within the program to check that certain requirements have been observed in the production of the lattice. The program 'GOMESH', can provide, superimposed in the graphical representation, stress and temperature values in numerical form, can represent the displacement of the mesh before and after a specified irradiation time, and give the directions and sense of the principal stresses occurring in the individual elements, in the form of arrows of varying length

3D visualization and finite element mesh formation from wood anatomy samples, Part II – Algorithm approach

Directory of Open Access Journals (Sweden)

Petr Koňas

2009-01-01

Full Text Available Paper presents new original application WOOD3D in form of program code assembling. The work extends the previous article “Part I – Theoretical approach” in detail description of implemented C++ classes of utilized projects Visualization Toolkit (VTK, Insight Toolkit (ITK and MIMX. Code is written in CMake style and it is available as multiplatform application. Currently GNU Linux (32/64b and MS Windows (32/64b platforms were released. Article discusses various filter classes for image filtering. Mainly Otsu and Binary threshold filters are classified for anatomy wood samples thresholding. Registration of images series is emphasized for difference of colour spaces compensation is included. Resulted work flow of image analysis is new methodological approach for images processing through the composition, visualization, filtering, registration and finite element mesh formation. Application generates script in ANSYS parametric design language (APDL which is fully compatible with ANSYS finite element solver and designer environment. The script includes the whole definition of unstructured finite element mesh formed by individual elements and nodes. Due to simple notation, the same script can be used for generation of geometrical entities in element positions. Such formed volumetric entities are prepared for further geometry approximation (e.g. by boolean or more advanced methods. Hexahedral and tetrahedral types of mesh elements are formed on user request with specified mesh options. Hexahedral meshes are formed both with uniform element size and with anisotropic character. Modified octree method for hexahedral mesh with anisotropic character was declared in application. Multicore CPUs in the application are supported for fast image analysis realization. Visualization of image series and consequent 3D image are realized in VTK format sufficiently known and public format, visualized in GPL application Paraview. Future work based on mesh

Hybrid finite element/waveguide mode analysis of passive RF devices

Science.gov (United States)

McGrath, Daniel T.

1993-07-01

A numerical solution for time-harmonic electromagnetic fields in two-port passive radio frequency (RF) devices has been developed, implemented in a computer code, and validated. Vector finite elements are used to represent the fields in the device interior, and field continuity across waveguide apertures is enforced by matching the interior solution to a sum of waveguide modes. Consequently, the mesh may end at the aperture instead of extending into the waveguide. The report discusses the variational formulation and its reduction to a linear system using Galerkin's method. It describes the computer code, including its interface to commercial CAD software used for geometry generation. It presents validation results for waveguide discontinuities, coaxial transitions, and microstrip circuits. They demonstrate that the method is an effective and versatile tool for predicting the performance of passive RF devices.

Finite size effects of a pion matrix element

International Nuclear Information System (INIS)

Guagnelli, M.; Jansen, K.; Palombi, F.; Petronzio, R.; Shindler, A.; Wetzorke, I.

2004-01-01

We investigate finite size effects of the pion matrix element of the non-singlet, twist-2 operator corresponding to the average momentum of non-singlet quark densities. Using the quenched approximation, they come out to be surprisingly large when compared to the finite size effects of the pion mass. As a consequence, simulations of corresponding nucleon matrix elements could be affected by finite size effects even stronger which could lead to serious systematic uncertainties in their evaluation

Effect of analysis parameters on non-linear implicit finite element analysis of marine corroded steel plate

Science.gov (United States)

Islam, Muhammad Rabiul; Sakib-Ul-Alam, Md.; Nazat, Kazi Kaarima; Hassan, M. Munir

2017-12-01

FEA results greatly depend on analysis parameters. MSC NASTRAN nonlinear implicit analysis code has been used in large deformation finite element analysis of pitted marine SM490A steel rectangular plate. The effect of two types actual pit shape on parameters of integrity of structure has been analyzed. For 3-D modeling, a proposed method for simulation of pitted surface by probabilistic corrosion model has been used. The result has been verified with the empirical formula proposed by finite element analysis of steel surface generated with different pitted data where analyses have been carried out by the code of LS-DYNA 971. In the both solver, an elasto-plastic material has been used where an arbitrary stress versus strain curve can be defined. In the later one, the material model is based on the J2 flow theory with isotropic hardening where a radial return algorithm is used. The comparison shows good agreement between the two results which ensures successful simulation with comparatively less energy and time.

Verification of Orthogrid Finite Element Modeling Techniques

Science.gov (United States)

Steeve, B. E.

1996-01-01

The stress analysis of orthogrid structures, specifically with I-beam sections, is regularly performed using finite elements. Various modeling techniques are often used to simplify the modeling process but still adequately capture the actual hardware behavior. The accuracy of such 'Oshort cutso' is sometimes in question. This report compares three modeling techniques to actual test results from a loaded orthogrid panel. The finite element models include a beam, shell, and mixed beam and shell element model. Results show that the shell element model performs the best, but that the simpler beam and beam and shell element models provide reasonable to conservative results for a stress analysis. When deflection and stiffness is critical, it is important to capture the effect of the orthogrid nodes in the model.

Review of Tomographic Imaging using Finite Element Method

Directory of Open Access Journals (Sweden)

Mohd Fua’ad RAHMAT

2011-12-01

Full Text Available Many types of techniques for process tomography were proposed and developed during the past 20 years. This paper review the techniques and the current state of knowledge and experience on the subject, aimed at highlighting the problems associated with the non finite element methods, such as the ill posed, ill conditioned which relates to the accuracy and sensitivity of measurements. In this paper, considerations for choice of sensors and its applications were outlined and descriptions of non finite element tomography systems were presented. The finite element method tomography system as obtained from recent works, suitable for process control and measurement were also presented.

New mixed finite-element methods

International Nuclear Information System (INIS)

Franca, L.P.

1987-01-01

New finite-element methods are proposed for mixed variational formulations. The methods are constructed by adding to the classical Galerkin method various least-squares like terms. The additional terms involve integrals over element interiors, and include mesh-parameter dependent coefficients. The methods are designed to enhance stability. Consistency is achieved in the sense that exact solutions identically satisfy the variational equations.Applied to several problems, simple finite-element interpolations are rendered convergent, including convenient equal-order interpolations generally unstable within the Galerkin approach. The methods are subdivided into two classes according to the manner in which stability is attained: (1) circumventing Babuska-Brezzi condition methods; (2) satisfying Babuska-Brezzi condition methods. Convergence is established for each class of methods. Applications of the first class of methods to Stokes flow and compressible linear elasticity are presented. The second class of methods is applied to the Poisson, Timoshenko beam and incompressible elasticity problems. Numerical results demonstrate the good stability and accuracy of the methods, and confirm the error estimates

Parallel Object-Oriented Computation Applied to a Finite Element Problem

Directory of Open Access Journals (Sweden)

Jon B. Weissman

1993-01-01

Full Text Available The conventional wisdom in the scientific computing community is that the best way to solve large-scale numerically intensive scientific problems on today's parallel MIMD computers is to use Fortran or C programmed in a data-parallel style using low-level message-passing primitives. This approach inevitably leads to nonportable codes and extensive development time, and restricts parallel programming to the domain of the expert programmer. We believe that these problems are not inherent to parallel computing but are the result of the programming tools used. We will show that comparable performance can be achieved with little effort if better tools that present higher level abstractions are used. The vehicle for our demonstration is a 2D electromagnetic finite element scattering code we have implemented in Mentat, an object-oriented parallel processing system. We briefly describe the application. Mentat, the implementation, and present performance results for both a Mentat and a hand-coded parallel Fortran version.

On higher order pyramidal finite elements

Czech Academy of Sciences Publication Activity Database

Liu, L.; Davies, K.B.; Křížek, Michal; Guan, L.

2011-01-01

Roč. 3, č. 2 (2011), s. 131-140 ISSN 2070-0733 R&D Projects: GA AV ČR(CZ) IAA100190803 Institutional research plan: CEZ:AV0Z10190503 Keywords : pyramidal polynomial basis functions * finite element method * composite elements * three-dimensional mortar elements Subject RIV: BA - General Mathematics Impact factor: 0.750, year: 2011

Finite element analysis for the initiation of lamellar tearing in welded joints

International Nuclear Information System (INIS)

Krieg, R.D.; Thomas, R.K.

1980-01-01

A numerical procedure using the finite element method is presented for predicting the initiation of lamellar tearing in fillet welded T-joints commonly employed in large structures. Starting with a prescribed geometry, the welding process is approximated by a known time-dependent volumetric heat source which simulates the arc heating and deposition of liquid metal. The transient nonlinear thermal and stress problems are then solved using finite element computer codes. Results of the elastic-plastic stress analysis are presented showing a predicted region of incipient lamellar tearing based on a ductile fracture theory which qualitatively agrees with the size and location of tears typically observed in photomicrographs. Additional insight into post failure crack length and stability is presented based on a simplified linear elastic fracture mechanics approach. Although the analysis procedure shows signs of promise, several weak points in the model are pointed out which must be improved before lamellar tearing can be quantified in an approach of this general type

Three-dimensional finite element impact analysis of a nuclear waste truck cask

International Nuclear Information System (INIS)

Miller, J.D.

1985-01-01

This paper presents a three-dimensional finite element impact analysis of a hypothetical accident event for the preliminary design of a shipping cask which is used to transport radioactive waste by standard tractor-semitrailer truck. The nonlinear dynamic structural analysis code DYNA3D run on Sandia's Cray-1 computer was used to calculate the effects of the cask's closure-end impacting a rigid frictionless surface on an edge of its external impact limiter after a 30-foot fall. The center of gravity of the cask (made of 304 stainless steel and depleted uranium) was assumed to be directly above the impact point. An elastic-plastic material constitutive model was used to calculate the nonlinear response of the cask components to the transient loading. Interactive color graphics (PATRAN and MOVIE BYU) were used throughout the analysis, proving to be extremely helpful for generation and verification of the geometry and boundary conditions of the finite element model and for interpretation of the analysis results. Results from the calculations show the cask sustained large localized deformations. However, these were almost entirely confined to the impact limiters built into the cask. The closure sections were determined to remain intact, and leakage would not be expected after the event. As an example of a large three-dimensional finite element dynamic impact calculation, this analysis can serve as an excellent benchmark for computer aided design procedures

Finite element approach to global gyrokinetic particle-in-cell simulations using magnetic coordinate

International Nuclear Information System (INIS)

Fivaz, M.; Brunner, S.; Ridder, G. de; Sauter, O.; Tran, T.M.; Vaclavik, J.; Villard, L.; Appert, K.

1997-08-01

We present a fully-global linear gyrokinetic simulation code (GYGLES) aimed at describing the instable spectrum of the ion-temperature-gradient modes in toroidal geometry. We formulate the Particle-In-Cell method with finite elements defined in magnetic coordinates, which provides excellent numerical convergence properties. The poloidal mode structure corresponding to k // =0 is extracted without approximation from the equations, which reduces drastically the numerical resolution needed. The code can simulate routinely modes with both very long and very short toroidal wavelengths, can treat realistic (MHD) equilibria of any size and runs on a massively parallel computer. (author) 10 figs., 28 refs

Validating the Galerkin least-squares finite element methods in predicting mixing flows in stirred tank reactors

International Nuclear Information System (INIS)

Johnson, K.; Bittorf, K.J.

2002-01-01

A novel approach for computer aided modeling and optimizing mixing process has been developed using Galerkin least-squares finite element technology. Computer aided mixing modeling and analysis involves Lagrangian and Eulerian analysis for relative fluid stretching, and energy dissipation concepts for laminar and turbulent flows. High quality, conservative, accurate, fluid velocity, and continuity solutions are required for determining mixing quality. The ORCA Computational Fluid Dynamics (CFD) package, based on a finite element formulation, solves the incompressible Reynolds Averaged Navier Stokes (RANS) equations. Although finite element technology has been well used in areas of heat transfer, solid mechanics, and aerodynamics for years, it has only recently been applied to the area of fluid mixing. ORCA, developed using the Galerkin Least-Squares (GLS) finite element technology, provides another formulation for numerically solving the RANS based and LES based fluid mechanics equations. The ORCA CFD package is validated against two case studies. The first, a free round jet, demonstrates that the CFD code predicts the theoretical velocity decay rate, linear expansion rate, and similarity profile. From proper prediction of fundamental free jet characteristics, confidence can be derived when predicting flows in a stirred tank, as a stirred tank reactor can be considered a series of free jets and wall jets. (author)

The finite-difference and finite-element modeling of seismic wave propagation and earthquake motion

International Nuclear Information System (INIS)

Moszo, P.; Kristek, J.; Galis, M.; Pazak, P.; Balazovijech, M.

2006-01-01

Numerical modeling of seismic wave propagation and earthquake motion is an irreplaceable tool in investigation of the Earth's structure, processes in the Earth, and particularly earthquake phenomena. Among various numerical methods, the finite-difference method is the dominant method in the modeling of earthquake motion. Moreover, it is becoming more important in the seismic exploration and structural modeling. At the same time we are convinced that the best time of the finite-difference method in seismology is in the future. This monograph provides tutorial and detailed introduction to the application of the finite-difference, finite-element, and hybrid finite-difference-finite-element methods to the modeling of seismic wave propagation and earthquake motion. The text does not cover all topics and aspects of the methods. We focus on those to which we have contributed. (Author)

A set of pathological tests to validate new finite elements

Indian Academy of Sciences (India)

M. Senthilkumar (Newgen Imaging) 1461 1996 Oct 15 13:05:22

The finite element method entails several approximations. Hence it ... researchers have designed several pathological tests to validate any new finite element. The .... Three dimensional thick shell elements using a hybrid/mixed formu- lation.

Experimental and finite element prediction of bursting pressure in compound cylinders

International Nuclear Information System (INIS)

Majzoobi, G.H.; Farrahi, G.H.; Pipelzadeh, M.K.; Akbari, A.

2004-01-01

Aluminium cylinders with a constant ratio of outer to inner radii, k=2.2, with different diametral interferences and various shrinkage radii were subjected to bursting and autofrettage pressures. Numerical simulations of the compound cylinders were also performed using the finite element code, NISA. The results can predict the optimum shrinkage radius to a reasonable accuracy with the use of finite element analysis. This radius corresponds to the situation when the maximum von-Mises stress at the internal radii of both the inner and outer cylinders become equal. It was shown that the maximum von-Mises stress across the wall of the cylinder is at the minimum at this shrinkage radius. The optimum diametral interference was found to be that which sufficiently brought the contact surface of the inner and outer cylinders to the point of yielding. Should the shrinkage pressure exceed the elastic limit, the pressure capacity of the cylinder will not be improved. The numerical and experimental results show that autofrettage had no effect on the bursting pressure of the thick-walled compound cylinder for the material tested

Crack Propagation by Finite Element Method

Directory of Open Access Journals (Sweden)

Luiz Carlos H. Ricardo

2018-01-01

Full Text Available Crack propagation simulation began with the development of the finite element method; the analyses were conducted to obtain a basic understanding of the crack growth. Today structural and materials engineers develop structures and materials properties using this technique. The aim of this paper is to verify the effect of different crack propagation rates in determination of crack opening and closing stress of an ASTM specimen under a standard suspension spectrum loading from FDandE SAE Keyhole Specimen Test Load Histories by finite element analysis. To understand the crack propagation processes under variable amplitude loading, retardation effects are observed

Finite element analysis of volumetrically heated fluids in an axisymmetric enclosure

International Nuclear Information System (INIS)

Gartling, D.K.

1979-01-01

A general purpose finite element computer code has been used to analyze the steady state and transient response of a confined fluid that is heated volumetrically. The numerical procedure is demonstrated to be capable of resolving flow fields of considerable complexity without undue computational expense. Results are discussed for a Grashof number range (4.0 x 10 4 to 4.0 x 10 6 ) in which the flow varies from a steady, single cell configuration to a multiple cell configuration that includes a periodic interaction

Evaluation on the structural soundness of the package for subsurface disposal by finite element method

International Nuclear Information System (INIS)

Itoh, Chihiro

2009-01-01

The structural analysis of the disposal package for low-level radioactive wastes with relatively high activities (called L1 waste in Japan) were performed against normal and hypothetical conditions. As a normal condition the external load due to lifting, stacking of the package and filling the space of disposal pit with mortar or something were considered. On the other hand, drop incident during handling and pressure due to some external force were taken up as hypothetical conditions. Using finite element code ABAQUS and three dimensional finite element model, structural analyses were carried out for the normal conditions. The results show that the maximum stresses occurred at the package due to the loads above mentioned were far less than the yield strength for all conditions. Therefore, it is confirmed that the disposal package keeps its integrity under the normal conditions. Analyses for load cases of 9 m drop onto the reinforced concrete slab and 5.9 m drop onto the embedded disposal package were performed by using finite element code LS-DYNA. Both results show that the strains at the impact zone of the package exceeded the fracture strain of the material but the damaged area was limited in the vicinity of impact zone. As a maximum external pressure, 4MPa was applied to the surface of the packages which were piled up in four layered in the disposal tunnel. According to the results of analyses by ABAQUS code the maximum strain occurred at the contact surfaces close to the welding zone between lid and body of the top package. However, the package stays in sound because the value of the maximum strain was less than the fracture strain of the materials. (author)

NONLINEAR FINITE ELEMENT ANALYSIS OF NONSEISMICALLY DETAILED INTERIOR RC BEAM-COLUMN CONNECTION UNDER REVERSED CYCLIC LOAD

Directory of Open Access Journals (Sweden)

Teeraphot Supaviriyakit

2017-11-01

Full Text Available This paper presents a nonlinear finite element analysis of non-seismically detailed RC beam column connections under reversed cyclic load. The test of half-scale nonductile reinforced concrete beam-column joints was conducted. The tested specimens represented those of the actual mid-rise reinforced concrete frame buildings designed according to the non-seismic provisions of the ACI building code.  The test results show that specimens representing small and medium column tributary area failed in brittle joint shear while specimen representing large column tributary area failed by ductile flexure though no ductile reinforcement details were provided. The nonlinear finite element analysis was applied to simulate the behavior of the specimens. The finite element analysis employs the smeared crack approach for modeling beam, column and joint, and employs the discrete crack approach for modeling the interface between beam and joint face. The nonlinear constitutive models of reinforced concrete elements consist of coupled tension-compression model to model normal force orthogonal and parallel to the crack and shear transfer model to capture the shear sliding mechanism. The FEM shows good comparison with test results in terms of load-displacement relations, hysteretic loops, cracking process and the failure mode of the tested specimens. The finite element analysis clarifies that the joint shear failure was caused by the collapse of principal diagonal concrete strut.

Finite element analysis of plastic recycling machine designed for ...

African Journals Online (AJOL)

... design was evaluated using finite element analysis (FEA) tool in Solid Works Computer ... Also, a minimum factor of safety value of 5.3 was obtained for shredder shaft ... Machine; Design; Recycling; Sustainability; Finite Element; Simulation ...

Charged particle tracking through electrostatic wire meshes using the finite element method

Energy Technology Data Exchange (ETDEWEB)

Devlin, L. J.; Karamyshev, O.; Welsch, C. P., E-mail: carsten.welsch@cockcroft.ac.uk [The Cockcroft Institute, Daresbury Laboratory, Warrington (United Kingdom); Department of Physics, University of Liverpool, Liverpool (United Kingdom)

2016-06-15

Wire meshes are used across many disciplines to accelerate and focus charged particles, however, analytical solutions are non-exact and few codes exist which simulate the exact fields around a mesh with physical sizes. A tracking code based in Matlab-Simulink using field maps generated using finite element software has been developed which tracks electrons or ions through electrostatic wire meshes. The fields around such a geometry are presented as an analytical expression using several basic assumptions, however, it is apparent that computational calculations are required to obtain realistic values of electric potential and fields, particularly when multiple wire meshes are deployed. The tracking code is flexible in that any quantitatively describable particle distribution can be used for both electrons and ions as well as other benefits such as ease of export to other programs for analysis. The code is made freely available and physical examples are highlighted where this code could be beneficial for different applications.

A particle finite element method for machining simulations

Science.gov (United States)

Sabel, Matthias; Sator, Christian; Müller, Ralf

2014-07-01

The particle finite element method (PFEM) appears to be a convenient technique for machining simulations, since the geometry and topology of the problem can undergo severe changes. In this work, a short outline of the PFEM-algorithm is given, which is followed by a detailed description of the involved operations. The -shape method, which is used to track the topology, is explained and tested by a simple example. Also the kinematics and a suitable finite element formulation are introduced. To validate the method simple settings without topological changes are considered and compared to the standard finite element method for large deformations. To examine the performance of the method, when dealing with separating material, a tensile loading is applied to a notched plate. This investigation includes a numerical analysis of the different meshing parameters, and the numerical convergence is studied. With regard to the cutting simulation it is found that only a sufficiently large number of particles (and thus a rather fine finite element discretisation) leads to converged results of process parameters, such as the cutting force.

A finite element primer for beginners the basics

CERN Document Server

Zohdi, Tarek I

2014-01-01

The purpose of this primer is to provide the basics of the Finite Element Method, primarily illustrated through a classical model problem, linearized elasticity. The topics covered are:(1) Weighted residual methods and Galerkin approximations,(2) A model problem for one-dimensional?linear elastostatics,(3) Weak formulations in one dimension,(4) Minimum principles in one dimension,(5) Error estimation in one dimension,(5) Construction of Finite Element basis functions in one dimension,(6) Gaussian Quadrature,(7) Iterative solvers and element by element data structures,(8) A model problem for th

On the reliability of finite element solutions

International Nuclear Information System (INIS)

Prasad, K.S.R.K.

1975-01-01

The extent of reliability of the finite element method for analysis of nuclear reactor structures, and that of reactor vessels in particular and the need for the engineer to guard against the pitfalls that may arise out of both physical and mathematical models have been high-lighted. A systematic way of checking the model to obtain reasonably accurate solutions is presented. Quite often sophisticated elements are suggested for specific design and stress concentration problems. The desirability or otherwise of these elements, their scope and utility vis-a-vis the use of large stack of conventional elements are discussed from the view point of stress analysts. The methods of obtaining a check on the reliability of the finite element solutions either through modelling changes or an extrapolation technique are discussed. (author)

Adaptive Smoothed Finite Elements (ASFEM) for history dependent material models

International Nuclear Information System (INIS)

Quak, W.; Boogaard, A. H. van den

2011-01-01

A successful simulation of a bulk forming process with finite elements can be difficult due to distortion of the finite elements. Nodal smoothed Finite Elements (NSFEM) are an interesting option for such a process since they show good distortion insensitivity and moreover have locking-free behavior and good computational efficiency. In this paper a method is proposed which takes advantage of the nodally smoothed field. This method, named adaptive smoothed finite elements (ASFEM), revises the mesh for every step of a simulation without mapping the history dependent material parameters. In this paper an updated-Lagrangian implementation is presented. Several examples are given to illustrate the method and to show its properties.

Finite Element Based Design and Optimization for Piezoelectric Accelerometers

DEFF Research Database (Denmark)

Liu, Bin; Kriegbaum, B.; Yao, Q.

1998-01-01

A systematic Finite Element design and optimisation procedure is implemented for the development of piezoelectric accelerometers. Most of the specifications of accelerometers can be obtained using the Finite Element simulations. The deviations between the simulated and calibrated sensitivities...

THE PRACTICAL ANALYSIS OF FINITE ELEMENTS METHOD ERRORS

Directory of Open Access Journals (Sweden)

Natalia Bakhova

2011-03-01

Full Text Available Abstract. The most important in the practical plan questions of reliable estimations of finite elementsmethod errors are considered. Definition rules of necessary calculations accuracy are developed. Methodsand ways of the calculations allowing receiving at economical expenditures of computing work the best finalresults are offered.Keywords: error, given the accuracy, finite element method, lagrangian and hermitian elements.

High accuracy 3D electromagnetic finite element analysis

International Nuclear Information System (INIS)

Nelson, Eric M.

1997-01-01

A high accuracy 3D electromagnetic finite element field solver employing quadratic hexahedral elements and quadratic mixed-order one-form basis functions will be described. The solver is based on an object-oriented C++ class library. Test cases demonstrate that frequency errors less than 10 ppm can be achieved using modest workstations, and that the solutions have no contamination from spurious modes. The role of differential geometry and geometrical physics in finite element analysis will also be discussed

High accuracy 3D electromagnetic finite element analysis

International Nuclear Information System (INIS)

Nelson, E.M.

1996-01-01

A high accuracy 3D electromagnetic finite element field solver employing quadratic hexahedral elements and quadratic mixed-order one-form basis functions will be described. The solver is based on an object-oriented C++ class library. Test cases demonstrate that frequency errors less than 10 ppm can be achieved using modest workstations, and that the solutions have no contamination from spurious modes. The role of differential geometry and geometrical physics in finite element analysis will also be discussed

ABOUT SOLUTION OF MULTIPOINT BOUNDARY PROBLEMS OF TWO-DIMENSIONAL STRUCTURAL ANALYSIS WITH THE USE OF COMBINED APPLICATION OF FINITE ELEMENT METHOD AND DISCRETE-CONTINUAL FINITE ELEMENT METHOD PART 2: SPECIAL ASPECTS OF FINITE ELEMENT APPROXIMATION

Directory of Open Access Journals (Sweden)

Pavel A. Akimov

2017-12-01

Full Text Available As is well known, the formulation of a multipoint boundary problem involves three main components: a description of the domain occupied by the structure and the corresponding subdomains; description of the conditions inside the domain and inside the corresponding subdomains, the description of the conditions on the boundary of the domain, conditions on the boundaries between subdomains. This paper is a continuation of another work published earlier, in which the formulation and general principles of the approximation of the multipoint boundary problem of a static analysis of deep beam on the basis of the joint application of the finite element method and the discrete-continual finite element method were considered. It should be noted that the approximation within the fragments of a domain that have regular physical-geometric parameters along one of the directions is expedient to be carried out on the basis of the discrete-continual finite element method (DCFEM, and for the approximation of all other fragments it is necessary to use the standard finite element method (FEM. In the present publication, the formulas for the computing of displacements partial derivatives of displacements, strains and stresses within the finite element model (both within the finite element and the corresponding nodal values (with the use of averaging are presented. Boundary conditions between subdomains (respectively, discrete models and discrete-continual models and typical conditions such as “hinged support”, “free edge”, “perfect contact” (twelve basic (basic variants are available are under consideration as well. Governing formulas for computing of elements of the corresponding matrices of coefficients and vectors of the right-hand sides are given for each variant. All formulas are fully adapted for algorithmic implementation.

Finite element simulation and testing of ISW CFRP anchorage

DEFF Research Database (Denmark)

Schmidt, Jacob Wittrup; Goltermann, Per; Hertz, Kristian Dahl

2013-01-01

is modelled in the 3D finite Element program ABAQUS, just as digital image correlation (DIC) testing was performed to verify the finite element simulation. Also a new optimized design was produced to ensure that the finite element simulation and anchorage behaviour correlated well. It is seen....... This paper presents a novel mechanical integrated sleeve wedge anchorage which seem very promising when perusing the scope of ultimate utilization of CFRP 8mm rods (with a tension capacity of approximately 140kN). Compression transverse to the CFRP is evaluated to prevent premature failure. The anchorage...

Finite element analysis of a finite-strain plasticity problem

International Nuclear Information System (INIS)

Crose, J.G.; Fong, H.H.

1984-01-01

A finite-strain plasticity analysis was performed of an engraving process in a plastic rotating band during the firing of a gun projectile. The aim was to verify a nonlinear feature of the NIFDI/RB code: plastic large deformation analysis of nearly incompressible materials using a deformation theory of plasticity approach and a total Lagrangian scheme. (orig.)

Hydrothermal analysis in engineering using control volume finite element method

CERN Document Server

Sheikholeslami, Mohsen

2015-01-01

Control volume finite element methods (CVFEM) bridge the gap between finite difference and finite element methods, using the advantages of both methods for simulation of multi-physics problems in complex geometries. In Hydrothermal Analysis in Engineering Using Control Volume Finite Element Method, CVFEM is covered in detail and applied to key areas of thermal engineering. Examples, exercises, and extensive references are used to show the use of the technique to model key engineering problems such as heat transfer in nanofluids (to enhance performance and compactness of energy systems),

A multiscale mortar multipoint flux mixed finite element method

KAUST Repository

Wheeler, Mary Fanett; Xue, Guangri; Yotov, Ivan

2012-01-01

In this paper, we develop a multiscale mortar multipoint flux mixed finite element method for second order elliptic problems. The equations in the coarse elements (or subdomains) are discretized on a fine grid scale by a multipoint flux mixed finite

Parallel iterative procedures for approximate solutions of wave propagation by finite element and finite difference methods

Energy Technology Data Exchange (ETDEWEB)

Kim, S. [Purdue Univ., West Lafayette, IN (United States)

1994-12-31

Parallel iterative procedures based on domain decomposition techniques are defined and analyzed for the numerical solution of wave propagation by finite element and finite difference methods. For finite element methods, in a Lagrangian framework, an efficient way for choosing the algorithm parameter as well as the algorithm convergence are indicated. Some heuristic arguments for finding the algorithm parameter for finite difference schemes are addressed. Numerical results are presented to indicate the effectiveness of the methods.

Global-Local Finite Element Analysis of Bonded Single-Lap Joints

Science.gov (United States)

Kilic, Bahattin; Madenci, Erdogan; Ambur, Damodar R.

2004-01-01

Adhesively bonded lap joints involve dissimilar material junctions and sharp changes in geometry, possibly leading to premature failure. Although the finite element method is well suited to model the bonded lap joints, traditional finite elements are incapable of correctly resolving the stress state at junctions of dissimilar materials because of the unbounded nature of the stresses. In order to facilitate the use of bonded lap joints in future structures, this study presents a finite element technique utilizing a global (special) element coupled with traditional elements. The global element includes the singular behavior at the junction of dissimilar materials with or without traction-free surfaces.

Solving the transport equation with quadratic finite elements: Theory and applications

International Nuclear Information System (INIS)

Ferguson, J.M.

1997-01-01

At the 4th Joint Conference on Computational Mathematics, the author presented a paper introducing a new quadratic finite element scheme (QFEM) for solving the transport equation. In the ensuing year the author has obtained considerable experience in the application of this method, including solution of eigenvalue problems, transmission problems, and solution of the adjoint form of the equation as well as the usual forward solution. He will present detailed results, and will also discuss other refinements of his transport codes, particularly for 3-dimensional problems on rectilinear and non-rectilinear grids

Accelerated finite element elastodynamic simulations using the GPU

Energy Technology Data Exchange (ETDEWEB)

Huthwaite, Peter, E-mail: p.huthwaite@imperial.ac.uk

2014-01-15

An approach is developed to perform explicit time domain finite element simulations of elastodynamic problems on the graphical processing unit, using Nvidia's CUDA. Of critical importance for this problem is the arrangement of nodes in memory, allowing data to be loaded efficiently and minimising communication between the independently executed blocks of threads. The initial stage of memory arrangement is partitioning the mesh; both a well established ‘greedy’ partitioner and a new, more efficient ‘aligned’ partitioner are investigated. A method is then developed to efficiently arrange the memory within each partition. The software is applied to three models from the fields of non-destructive testing, vibrations and geophysics, demonstrating a memory bandwidth of very close to the card's maximum, reflecting the bandwidth-limited nature of the algorithm. Comparison with Abaqus, a widely used commercial CPU equivalent, validated the accuracy of the results and demonstrated a speed improvement of around two orders of magnitude. A software package, Pogo, incorporating these developments, is released open source, downloadable from (http://www.pogo-fea.com/) to benefit the community. -- Highlights: •A novel memory arrangement approach is discussed for finite elements on the GPU. •The mesh is partitioned then nodes are arranged efficiently within each partition. •Models from ultrasonics, vibrations and geophysics are run. •The code is significantly faster than an equivalent commercial CPU package. •Pogo, the new software package, is released open source.

Accelerated finite element elastodynamic simulations using the GPU

International Nuclear Information System (INIS)

Huthwaite, Peter

2014-01-01

An approach is developed to perform explicit time domain finite element simulations of elastodynamic problems on the graphical processing unit, using Nvidia's CUDA. Of critical importance for this problem is the arrangement of nodes in memory, allowing data to be loaded efficiently and minimising communication between the independently executed blocks of threads. The initial stage of memory arrangement is partitioning the mesh; both a well established ‘greedy’ partitioner and a new, more efficient ‘aligned’ partitioner are investigated. A method is then developed to efficiently arrange the memory within each partition. The software is applied to three models from the fields of non-destructive testing, vibrations and geophysics, demonstrating a memory bandwidth of very close to the card's maximum, reflecting the bandwidth-limited nature of the algorithm. Comparison with Abaqus, a widely used commercial CPU equivalent, validated the accuracy of the results and demonstrated a speed improvement of around two orders of magnitude. A software package, Pogo, incorporating these developments, is released open source, downloadable from (http://www.pogo-fea.com/) to benefit the community. -- Highlights: •A novel memory arrangement approach is discussed for finite elements on the GPU. •The mesh is partitioned then nodes are arranged efficiently within each partition. •Models from ultrasonics, vibrations and geophysics are run. •The code is significantly faster than an equivalent commercial CPU package. •Pogo, the new software package, is released open source

The finite element structural analysis code SAP IV conversion from CDC to IBM

International Nuclear Information System (INIS)

Harrop, L.P.

1977-02-01

SAP IV is a general three dimensional, linear, static and dynamic finite element structural analysis program. The program which was obtained from the Earthquake Engineering Research Center, University of California, Berkeley, was written in FORTRAM for a CDC 6400. Its main use was anticipated to be the seismic analysis of reactor structures. SAP IV may also prove useful for fracture mechanics studies as well as the usual elastic stress analysis of structures. A brief description of SAP IV and a more detailed account of the FORTRAN conversion required to make SAP IV run successfully on the UKAEA Harwell IBM 370/168 are given. (author)

Finite element simulation of laser transmission welding of dissimilar ...

African Journals Online (AJOL)

user

materials between polyvinylidene fluoride and titanium ... finite element (FE) thermal model is developed to simulate the laser ... Keywords: Laser transmission welding, Temperature field, Weld dimension, Finite element analysis, Thermal modeling. 1. .... 4) The heating phenomena due to the phase changes are neglected.

Finite element analysis of tibial fractures

DEFF Research Database (Denmark)

Wong, Christian Nai En; Mikkelsen, Mikkel Peter W; Hansen, Leif Berner

2010-01-01

Project. The data consisted of 21,219 3D elements with a cortical shell and a trabecular core. Three types of load of torsion, a direct lateral load and axial compression were applied. RESULTS: The finite element linear static analysis resulted in relevant fracture localizations and indicated relevant...

Generalized multiscale finite element method. Symmetric interior penalty coupling

KAUST Repository

Efendiev, Yalchin R.; Galvis, Juan; Lazarov, Raytcho D.; Moon, M.; Sarkis, Marcus V.

2013-01-01

Motivated by applications to numerical simulations of flows in highly heterogeneous porous media, we develop multiscale finite element methods for second order elliptic equations. We discuss a multiscale model reduction technique in the framework of the discontinuous Galerkin finite element method. We propose two different finite element spaces on the coarse mesh. The first space is based on a local eigenvalue problem that uses an interior weighted L2-norm and a boundary weighted L2-norm for computing the "mass" matrix. The second choice is based on generation of a snapshot space and subsequent selection of a subspace of a reduced dimension. The approximation with these multiscale spaces is based on the discontinuous Galerkin finite element method framework. We investigate the stability and derive error estimates for the methods and further experimentally study their performance on a representative number of numerical examples. © 2013 Elsevier Inc.

Generalized multiscale finite element method. Symmetric interior penalty coupling

KAUST Repository

Efendiev, Yalchin R.

2013-12-01

Motivated by applications to numerical simulations of flows in highly heterogeneous porous media, we develop multiscale finite element methods for second order elliptic equations. We discuss a multiscale model reduction technique in the framework of the discontinuous Galerkin finite element method. We propose two different finite element spaces on the coarse mesh. The first space is based on a local eigenvalue problem that uses an interior weighted L2-norm and a boundary weighted L2-norm for computing the "mass" matrix. The second choice is based on generation of a snapshot space and subsequent selection of a subspace of a reduced dimension. The approximation with these multiscale spaces is based on the discontinuous Galerkin finite element method framework. We investigate the stability and derive error estimates for the methods and further experimentally study their performance on a representative number of numerical examples. © 2013 Elsevier Inc.

Implementation of the full viscoresistive magnetohydrodynamic equations in a nonlinear finite element code

Energy Technology Data Exchange (ETDEWEB)

Haverkort, J.W. [Centrum Wiskunde & Informatica, P.O. Box 94079, 1090 GB Amsterdam (Netherlands); Dutch Institute for Fundamental Energy Research, P.O. Box 6336, 5600 HH Eindhoven (Netherlands); Blank, H.J. de [Dutch Institute for Fundamental Energy Research, P.O. Box 6336, 5600 HH Eindhoven (Netherlands); Huysmans, G.T.A. [ITER Organization, Route de Vinon sur Verdon, 13115 Saint Paul Lez Durance (France); Pratt, J. [Dutch Institute for Fundamental Energy Research, P.O. Box 6336, 5600 HH Eindhoven (Netherlands); Koren, B., E-mail: b.koren@tue.nl [Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven (Netherlands)

2016-07-01

Numerical simulations form an indispensable tool to understand the behavior of a hot plasma that is created inside a tokamak for providing nuclear fusion energy. Various aspects of tokamak plasmas have been successfully studied through the reduced magnetohydrodynamic (MHD) model. The need for more complete modeling through the full MHD equations is addressed here. Our computational method is presented along with measures against possible problems regarding pollution, stability, and regularity. The problem of ensuring continuity of solutions in the center of a polar grid is addressed in the context of a finite element discretization of the full MHD equations. A rigorous and generally applicable solution is proposed here. Useful analytical test cases are devised to verify the correct implementation of the momentum and induction equation, the hyperdiffusive terms, and the accuracy with which highly anisotropic diffusion can be simulated. A striking observation is that highly anisotropic diffusion can be treated with the same order of accuracy as isotropic diffusion, even on non-aligned grids, as long as these grids are generated with sufficient care. This property is shown to be associated with our use of a magnetic vector potential to describe the magnetic field. Several well-known instabilities are simulated to demonstrate the capabilities of the new method. The linear growth rate of an internal kink mode and a tearing mode are benchmarked against the results of a linear MHD code. The evolution of a tearing mode and the resulting magnetic islands are simulated well into the nonlinear regime. The results are compared with predictions from the reduced MHD model. Finally, a simulation of a ballooning mode illustrates the possibility to use our method as an ideal MHD method without the need to add any physical dissipation.

A practical method for computation of ductile crack growth by means of finite elements and parametric 3D-modelling

International Nuclear Information System (INIS)

Baumjohann, F.; Kroening, J.

1999-01-01

The present paper originates from a contribution to the safety assessment of a reactor pressure vessel (RPV). Investigations evaluating the safety against brittle fracture (exclosure of crack initiation and arrest assessments) are completed by calculations concerning ductile crack extension. Crack geometries including the expected crack extension are generated parametrically by a computer code and are used for further calculations with finite element programs. J-integrals of ductile growing cracks located between two comparative contours are determined by interpolation. The comparative contours are loaded by instationary temperature and pressure fields and are evaluated in advance. Taking the stability condition into consideration, the ductile crack extension is determined by pursuing the equilibrium between loading and crack resistance. The automatic modelling and a mathematical program processing the finite element results evaluate the crack growth of the finite element results very effectively. (orig.)

Finite Element Analysis of Circular Plate using SolidWorks

International Nuclear Information System (INIS)

Kang, Yeo Jin; Jhung, Myung Jo

2011-01-01

Circular plates are used extensively in mechanical engineering for nuclear reactor internal components. The examples in the reactor vessel internals are upper guide structure support plate, fuel alignment plate, lower support plate etc. To verify the structural integrity of these plates, the finite element analyses are performed, which require the development of the finite element model. Sometimes it is very costly and time consuming to make the model especially for the beginners who start their engineering job for the structural analysis, necessitating a simple method to develop the finite element model for the pursuing structural analysis. Therefore in this study, the input decks are generated for the finite element analysis of a circular plate as shown in Fig. 1, which can be used for the structural analysis such as modal analysis, response spectrum analysis, stress analysis, etc using the commercial program Solid Works. The example problems are solved and the results are included for analysts to perform easily the finite element analysis of the mechanical plate components due to various loadings. The various results presented in this study would be helpful not only for the benchmark calculations and results comparisons but also as a part of the knowledge management for the future generation of young designers, scientists and computer analysts

A multiscale mortar multipoint flux mixed finite element method

KAUST Repository

Wheeler, Mary Fanett

2012-02-03

In this paper, we develop a multiscale mortar multipoint flux mixed finite element method for second order elliptic problems. The equations in the coarse elements (or subdomains) are discretized on a fine grid scale by a multipoint flux mixed finite element method that reduces to cell-centered finite differences on irregular grids. The subdomain grids do not have to match across the interfaces. Continuity of flux between coarse elements is imposed via a mortar finite element space on a coarse grid scale. With an appropriate choice of polynomial degree of the mortar space, we derive optimal order convergence on the fine scale for both the multiscale pressure and velocity, as well as the coarse scale mortar pressure. Some superconvergence results are also derived. The algebraic system is reduced via a non-overlapping domain decomposition to a coarse scale mortar interface problem that is solved using a multiscale flux basis. Numerical experiments are presented to confirm the theory and illustrate the efficiency and flexibility of the method. © EDP Sciences, SMAI, 2012.

High accuracy 3D electromagnetic finite element analysis

International Nuclear Information System (INIS)

Nelson, E.M.

1997-01-01

A high accuracy 3D electromagnetic finite element field solver employing quadratic hexahedral elements and quadratic mixed-order one-form basis functions will be described. The solver is based on an object-oriented C++ class library. Test cases demonstrate that frequency errors less than 10 ppm can be achieved using modest workstations, and that the solutions have no contamination from spurious modes. The role of differential geometry and geometrical physics in finite element analysis will also be discussed. copyright 1997 American Institute of Physics

Physical analysis and modelling of aerosols transport. implementation in a finite elements code. Experimental validation in laminar and turbulent flows

International Nuclear Information System (INIS)

Armand, Patrick

1995-01-01

The aim of this work consists in the Fluid Mechanics and aerosol Physics coupling. In the first part, the order of magnitude analysis of the particle dynamics is done. This particle is embedded in a non-uniform unsteady flow. Flow approximations around the inclusion are described. Corresponding aerodynamic drag formulae are expressed. Possible situations related to the problem data are extensively listed. In the second part, one studies the turbulent particles transport. Eulerian approach which is particularly well adapted to industrial codes is preferred in comparison with the Lagrangian methods. One chooses the two-fluid formalism in which career gas-particles slip is taken into account. Turbulence modelling gets through a k-epsilon modulated by the inclusions action on the flow. The model is implemented In a finite elements code. Finally, In the third part, one validates the modelling in laminar and turbulent cases. We compare simulations to various experiments (settling battery, inertial impaction in a bend, jets loaded with glass beads particles) which are taken in the literature or done by ourselves at the laboratory. The results are very close. It is a good point when it is thought of the particles transport model and associated software future use. (author) [fr

finite element model for predicting residual stresses in shielded

African Journals Online (AJOL)

eobe

This paper investigates the prediction of residual stresses developed ... steel plates through Finite Element Model simulation and experiments. ... The experimental values as measured by the X-Ray diffractometer were of ... Based on this, it can be concluded that Finite Element .... Comparison of Residual Stresses from X.

Thermo-elastic stress analysis of containment wall penetrations using improved finite element formulation

International Nuclear Information System (INIS)

Ramani, D.T.; Dimopoulos, A.; Heglin, B.M.

1979-01-01

An increased application of finite element techniques, particularly in evaluating structural integrity of nuclear containment walls around penetration points, has aroused considerable interest. Due to extreme thermal effects in the vicinity of penetrations, the concrete containment wall is subject to unwarranted cracking effects, which must be controlled in accordance with ASME-III Code. This paper essentially deals with a unique finite element method of analysis in which nonlinear heat transfer problem across the penetration assembly in the nuclear containment drywell wall, is formulated. Using this technique, thermal analysis, dealing with an evaluation of temperature distribution around axisymmetric penetration assembly accomodating main steam lines or other vital piping at 600 0 F, is carried out. The method of analysis considers steady-state heat transfer energy balance across the process-pipe, insulation layer, guard-pipe sleeve, two intermediate air layers and an axisymmetric opening in the concrete containment wall, the outer faces of which are maintained at ambient temperature of 120 0 F. (orig.)

Finite-Element Modeling of Timber Joints with Punched Metal Plate Fasteners

DEFF Research Database (Denmark)

Ellegaard, Peter

2006-01-01

The focus of this paper is to describe the idea and the theory behind a finite-element model developed for analysis of timber trusses with punched metal plate fasteners (nail plates). The finite-element model includes the semirigid and nonlinear behavior of the joints (nonlinear nail and plate...... elements) and contact between timber beams, if any (bilinear contact elements). The timber beams have linear-elastic properties. The section forces needed for design of the joints are given directly by the finite-element model, since special elements are used to model the nail groups and the nail plate...... the behavior of the joints very well at lower load levels. At higher load levels the stiffness is overestimated due to development of cracks in the timber and the linear-elastic timber properties in the finite-element model....

Finite element analysis of inelastic structural behavior

International Nuclear Information System (INIS)

Argyris, J.H.; Szimmat, J.; Willam, K.J.

1977-01-01

The paper describes recent achievements in the finite element analysis of inelastic material behavior. The main purpose is to examine the interaction of three disciplines; (i) the finite element formulation of large deformation problems in the light of a systematic linearization, (ii) the constitutive modelling of inelastic processes in the rate-dependent and rate-independent response regime and (iii) the numerical solution of nonlinear rate problems via incremental iteration techniques. In the first part, alternative finite element models are developed for the idealization of large deformation problems. A systematic approach is presented to linearize the field equations locally by an incremental procedure. The finite element formulation is then examined for the description of inelastic material processes. In the second part, nonlinear and inelastic material phenomena are classified and illustrated with representative examples of concrete and metal components. In particular, rate-dependent and rate-independent material behavior is examined and representative constitutive models are assessed for their mathematical characterization. Hypoelastic, elastoplastic and endochronic models are compared for the description rate-independent material phenomena. In the third part, the numerial solution of inelastic structural behavior is discussed. In this context, several incremental techniques are developed and compared for tracing the evolution of the inelastic process. The numerical procedures are examined with regard to stability and accuracy to assess the overall efficiency. The 'optimal' incremental technique is then contrasted with the computer storage requirements to retain the data for the 'memory-characteristics' of the constitutive model

Numerical experiment on finite element method for matching data

International Nuclear Information System (INIS)

Tokuda, Shinji; Kumakura, Toshimasa; Yoshimura, Koichi.

1993-03-01

Numerical experiments are presented on the finite element method by Pletzer-Dewar for matching data of an ordinary differential equation with regular singular points by using model equation. Matching data play an important role in nonideal MHD stability analysis of a magnetically confined plasma. In the Pletzer-Dewar method, the Frobenius series for the 'big solution', the fundamental solution which is not square-integrable at the regular singular point, is prescribed. The experiments include studies of the convergence rate of the matching data obtained by the finite element method and of the effect on the results of computation by truncating the Frobenius series at finite terms. It is shown from the present study that the finite element method is an effective method for obtaining the matching data with high accuracy. (author)

Analytical and finite element modeling of grounding systems

Energy Technology Data Exchange (ETDEWEB)

Luz, Mauricio Valencia Ferreira da [University of Santa Catarina (UFSC), Florianopolis, SC (Brazil)], E-mail: mauricio@grucad.ufsc.br; Dular, Patrick [University of Liege (Belgium). Institut Montefiore], E-mail: Patrick.Dular@ulg.ac.be

2007-07-01

Grounding is the art of making an electrical connection to the earth. This paper deals with the analytical and finite element modeling of grounding systems. An electrokinetic formulation using a scalar potential can benefit from floating potentials to define global quantities such as electric voltages and currents. The application concerns a single vertical grounding with one, two and three-layer soil, where the superior extremity stays in the surface of the soil. This problem has been modeled using a 2D axi-symmetric electrokinetic formulation. The grounding resistance obtained by finite element method is compared with the analytical one for one-layer soil. With the results of this paper it is possible to show that finite element method is a powerful tool in the analysis of the grounding systems in low frequencies. (author)

Determination of Nonlinear Stiffness Coefficients for Finite Element Models with Application to the Random Vibration Problem

Science.gov (United States)

Muravyov, Alexander A.

1999-01-01

In this paper, a method for obtaining nonlinear stiffness coefficients in modal coordinates for geometrically nonlinear finite-element models is developed. The method requires application of a finite-element program with a geometrically non- linear static capability. The MSC/NASTRAN code is employed for this purpose. The equations of motion of a MDOF system are formulated in modal coordinates. A set of linear eigenvectors is used to approximate the solution of the nonlinear problem. The random vibration problem of the MDOF nonlinear system is then considered. The solutions obtained by application of two different versions of a stochastic linearization technique are compared with linear and exact (analytical) solutions in terms of root-mean-square (RMS) displacements and strains for a beam structure.

Robust mixed finite element methods to deal with incompressibility in finite strain in an industrial framework

International Nuclear Information System (INIS)

Al-Akhrass, Dina

2014-01-01

Simulations in solid mechanics exhibit several difficulties, as dealing with incompressibility, with nonlinearities due to finite strains, contact laws, or constitutive laws. The basic motivation of our work is to propose efficient finite element methods capable of dealing with incompressibility in finite strain context, and using elements of low order. During the three last decades, many approaches have been proposed in the literature to overcome the incompressibility problem. Among them, mixed formulations offer an interesting theoretical framework. In this work, a three-field mixed formulation (displacement, pressure, volumetric strain) is investigated. In some cases, this formulation can be condensed in a two-field (displacement - pressure) mixed formulation. However, it is well-known that the discrete problem given by the Galerkin finite element technique, does not inherit the 'inf-sup' stability condition from the continuous problem. Hence, the interpolation orders in displacement and pressure have to be chosen in a way to satisfy the Brezzi-Babuska stability conditions when using Galerkin approaches. Interpolation orders must be chosen so as to satisfy this condition. Two possibilities are considered: to use stable finite element satisfying this requirement, or to use finite element that does not satisfy this condition, and to add terms stabilizing the FE Galerkin formulation. The latter approach allows the use of equal order interpolation. In this work, stable finite element P2/P1 and P2/P1/P1 are used as reference, and compared to P1/P1 and P1/P1/P1 formulations stabilized with a bubble function or with a VMS method (Variational Multi-Scale) based on a sub-grid-space orthogonal to the FE space. A finite strain model based on logarithmic strain is selected. This approach is extended to three and two field mixed formulations with stable or stabilized elements. These approaches are validated on academic cases and used on industrial cases. (author)

ZONE: a finite element mesh generator

International Nuclear Information System (INIS)

Burger, M.J.

1976-05-01

The ZONE computer program is a finite-element mesh generator which produces the nodes and element description of any two-dimensional geometry. The geometry is subdivided into a mesh of quadrilateral and triangular zones arranged sequentially in an ordered march through the geometry. The order of march can be chosen so that the minimum bandwidth is obtained. The node points are defined in terms of the x and y coordinates in a global rectangular coordinate system. The zones generated are quadrilaterals or triangles defined by four node points in a counterclockwise sequence. Node points defining the outside boundary are generated to describe pressure boundary conditions. The mesh that is generated can be used as input to any two-dimensional as well as any axisymmetrical structure program. The output from ZONE is essentially the input file to NAOS, HONDO, and other axisymmetric finite element programs. 14 figures

FREESURF: A three-dimensional finite-element model for simulating groundwater flow into and around an excavation

Energy Technology Data Exchange (ETDEWEB)

Weitzman, Morley

1992-07-15

A three-dimensional finite-element code was developed and used to simulate the flow of groundwater towards an excavation in a saturated porous medium, allowing for seepage faces. An iterative procedure was used to predict the movement of the water table and the seepage flux. The numerical solution agreed well with experimental results from a sandbox experiment. (auth)

SEACAS Theory Manuals: Part III. Finite Element Analysis in Nonlinear Solid Mechanics

Energy Technology Data Exchange (ETDEWEB)

Laursen, T.A.; Attaway, S.W.; Zadoks, R.I.

1999-03-01

This report outlines the application of finite element methodology to large deformation solid mechanics problems, detailing also some of the key technological issues that effective finite element formulations must address. The presentation is organized into three major portions: first, a discussion of finite element discretization from the global point of view, emphasizing the relationship between a virtual work principle and the associated fully discrete system, second, a discussion of finite element technology, emphasizing the important theoretical and practical features associated with an individual finite element; and third, detailed description of specific elements that enjoy widespread use, providing some examples of the theoretical ideas already described. Descriptions of problem formulation in nonlinear solid mechanics, nonlinear continuum mechanics, and constitutive modeling are given in three companion reports.

Application of hexagonal element scheme in finite element method to three-dimensional diffusion problem of fast reactors

International Nuclear Information System (INIS)

Ishiguro, Misako; Higuchi, Kenji

1983-01-01

The finite element method is applied in Galerkin-type approximation to three-dimensional neutron diffusion equations of fast reactors. A hexagonal element scheme is adopted for treating the hexagonal lattice which is typical for fast reactors. The validity of the scheme is verified by applying the scheme as well as alternative schemes to the neutron diffusion calculation of a gas-cooled fast reactor of actual scale. The computed results are compared with corresponding values obtained using the currently applied triangular-element and also with conventional finite difference schemes. The hexagonal finite element scheme is found to yield a reasonable solution to the problem taken up here, with some merit in terms of saving in computing time, but the resulting multiplication factor differs by 1% and the flux by 9% compared with the triangular mesh finite difference scheme. The finite element method, even in triangular element scheme, would appear to incur error in inadmissible amount and which could not be easily eliminated by refining the nodes. (author)

Validation of finite element code DELFIN by means of the zero power experiences at the nuclear power plant of Atucha I; Convalidacion del codigo DELFIN por medio de las experiencias a potencia cero de la central nuclear Atucha I

Energy Technology Data Exchange (ETDEWEB)

Grant, C R [Comision Nacional de Energia Atomica, San Martin (Argentina). Unidad de Actividad Reactores y Centrales Nucleares

1997-12-31

Code DELFIN, developed in CNEA, treats the spatial discretization using heterogeneous finite elements, allowing a correct treatment of the continuity of fluxes and currents among elements and a more realistic representation of the hexagonal lattice of the reactor. It can be used for fuel management calculation, Xenon oscillation and spatial kinetics. Using the HUEMUL code for cell calculation (which uses a generalized two dimensional collision probability theory and has the WIMS library incorporated in a data base), the zero power experiences performed in 1974 were calculated. (author). 8 refs., 9 figs., 3 tabs.

Field load and displacement boundary condition computer program used for the finite element analysis and design of toroidal field coils in a tokamak

International Nuclear Information System (INIS)

Smith, R.A.

1975-06-01

The design evaluation of toroidal field coils on the Princeton Large Torus (PLT), the Poloidal Diverter Experiment (PDX) and the Tokamak Fusion Test Reactor (TFTR) has been performed by structural analysis with the finite element method. The technique employed has been simplified with supplementary computer programs that are used to generate the input data for the finite element computer program. Significant automation has been provided by computer codes in three areas of data input. These are the definition of coil geometry by a mesh of node points, the definition of finite elements via the node points and the definition of the node point force/displacement boundary conditions. The computer programs by name that have been used to perform the above functions are PDXNODE, ELEMENT and PDXFORC. The geometric finite element modeling options for toroidal field coils provided by PDXNODE include one-fourth or one-half symmetric sections of circular coils, oval shaped coils or dee-shaped coils with or without a beveled wedging surface. The program ELEMENT which defines the finite elements for input to the finite element computer code can provide considerable time and labor savings when defining the model of coils of non-uniform cross-section or when defining the model of coils whose material properties are different in the R and THETA directions due to the laminations of alternate epoxy and copper windings. The modeling features provided by the program ELEMENT have been used to analyze the PLT and the TFTR toroidal field coils with integral support structures. The computer program named PDXFORC is described. It computes the node point forces in a model of a toroidal field coil from the vector crossproduct of the coil current and the magnetic field. The model can be of one-half or one-fourth symmetry to be consistent with the node model defined by PDXNODE, and the magnetic field is computed from toroidal or poloidal coils

Finite element approximation to the even-parity transport equation

International Nuclear Information System (INIS)

Lewis, E.E.

1981-01-01

This paper studies the finite element method, a procedure for reducing partial differential equations to sets of algebraic equations suitable for solution on a digital computer. The differential equation is cast into the form of a variational principle, the resulting domain then subdivided into finite elements. The dependent variable is then approximated by a simple polynomial, and these are linked across inter-element boundaries by continuity conditions. The finite element method is tailored to a variety of transport problems. Angular approximations are formulated, and the extent of ray effect mitigation is examined. Complex trial functions are introduced to enable the inclusion of buckling approximations. The ubiquitous curved interfaces of cell calculations, and coarse mesh methods are also treated. A concluding section discusses limitations of the work to date and suggests possible future directions

Linear finite element method for one-dimensional diffusion problems

Energy Technology Data Exchange (ETDEWEB)

Brandao, Michele A.; Dominguez, Dany S.; Iglesias, Susana M., E-mail: micheleabrandao@gmail.com, E-mail: dany@labbi.uesc.br, E-mail: smiglesias@uesc.br [Universidade Estadual de Santa Cruz (LCC/DCET/UESC), Ilheus, BA (Brazil). Departamento de Ciencias Exatas e Tecnologicas. Laboratorio de Computacao Cientifica

2011-07-01

We describe in this paper the fundamentals of Linear Finite Element Method (LFEM) applied to one-speed diffusion problems in slab geometry. We present the mathematical formulation to solve eigenvalue and fixed source problems. First, we discretized a calculus domain using a finite set of elements. At this point, we obtain the spatial balance equations for zero order and first order spatial moments inside each element. Then, we introduce the linear auxiliary equations to approximate neutron flux and current inside the element and architect a numerical scheme to obtain the solution. We offer numerical results for fixed source typical model problems to illustrate the method's accuracy for coarse-mesh calculations in homogeneous and heterogeneous domains. Also, we compare the accuracy and computational performance of LFEM formulation with conventional Finite Difference Method (FDM). (author)

Quadrature representation of finite element variational forms

DEFF Research Database (Denmark)

Ølgaard, Kristian Breum; Wells, Garth N.

2012-01-01

This chapter addresses the conventional run-time quadrature approach for the numerical integration of local element tensors associated with finite element variational forms, and in particular automated optimizations that can be performed to reduce the number of floating point operations...

Explicit dynamics for numerical simulation of crack propagation by the extended finite element method

International Nuclear Information System (INIS)

Menouillard, T.

2007-09-01

Computerized simulation is nowadays an integrating part of design and validation processes of mechanical structures. Simulation tools are more and more performing allowing a very acute description of the phenomena. Moreover, these tools are not limited to linear mechanics but are developed to describe more difficult behaviours as for instance structures damage which interests the safety domain. A dynamic or static load can thus lead to a damage, a crack and then a rupture of the structure. The fast dynamics allows to simulate 'fast' phenomena such as explosions, shocks and impacts on structure. The application domain is various. It concerns for instance the study of the lifetime and the accidents scenario of the nuclear reactor vessel. It is then very interesting, for fast dynamics codes, to be able to anticipate in a robust and stable way such phenomena: the assessment of damage in the structure and the simulation of crack propagation form an essential stake. The extended finite element method has the advantage to break away from mesh generation and from fields projection during the crack propagation. Effectively, crack is described kinematically by an appropriate strategy of enrichment of supplementary freedom degrees. Difficulties connecting the spatial discretization of this method with the temporal discretization of an explicit calculation scheme has then been revealed; these difficulties are the diagonal writing of the mass matrix and the associated stability time step. Here are presented two methods of mass matrix diagonalization based on the kinetic energy conservation, and studies of critical time steps for various enriched finite elements. The interest revealed here is that the time step is not more penalizing than those of the standard finite elements problem. Comparisons with numerical simulations on another code allow to validate the theoretical works. A crack propagation test in mixed mode has been exploited in order to verify the simulation

Development of a partitioned finite volume-finite element fluid-structure interaction scheme for strongly-coupled problems

CSIR Research Space (South Africa)

Suliman, Ridhwaan

2012-07-01

Full Text Available -linear deformations are accounted for. As will be demonstrated, the finite volume approach exhibits similar disad- vantages to the linear Q4 finite element formulation when undergoing bending. An enhanced finite volume approach is discussed and compared with finite...

The finite element analysis program MSC Marc/Mentat a first introduction

CERN Document Server

Öchsner, Andreas

2016-01-01

Based on simple examples, this book offers a short introduction to the general-purpose finite element program MSC Marc, a specialized program for non-linear problems (implicit solver) distributed by the MSC Software Corporation, which is commonly used in academia and industry. Today the documentation of all finite element programs includes a variety of step-by-step examples of differing complexity, and in addition, all software companies offer professional workshops on different topics. As such, rather than competing with these, the book focuses on providing simple examples, often single-element problems, which can easily be related to the theory that is discussed in finite element lectures. This makes it an ideal companion book to classical introductory courses on the finite element method.

A finite element thermohydrodynamic analyis of profile bore bearing

International Nuclear Information System (INIS)

Shah Nor bin Basri

1994-01-01

A finite element-based method is presented for analysing the thermohydrodynamic (THD) behaviour of profile bore bearing. A variational statement for the governing equation is derived and used to formulate a non-linear quadrilateral finite element of serendipity family. The predicted behaviour is compared with experimental evidence where possible and favorable correlation is obtained

Hualien forced vibration calculation with a finite element model

International Nuclear Information System (INIS)

Wang, F.; Gantenbein, F.; Nedelec, M.; Duretz, Ch.

1995-01-01

The forced vibration tests of the Hualien mock-up were useful to validate finite element models developed for soil-structure interaction. In this paper the two sets of tests with and without backfill were analysed. the methods used are based on finite element modeling for the soil. Two approaches were considered: calculation of soil impedance followed by the calculation of the transfer functions with a model taking into account the superstructure and the impedance; direct calculation of the soil-structure transfer functions, with the soil and the structure being represented in the same model by finite elements. Blind predictions and post-test calculations are presented and compared with the test results. (author). 4 refs., 8 figs., 2 tabs

The Finite Element Numerical Modelling of 3D Magnetotelluric

Directory of Open Access Journals (Sweden)

Ligang Cao

2014-01-01

Full Text Available The ideal numerical simulation of 3D magnetotelluric was restricted by the methodology complexity and the time-consuming calculation. Boundary values, the variation of weighted residual equation, and the hexahedral mesh generation method of finite element are three major causes. A finite element method for 3D magnetotelluric numerical modeling is presented in this paper as a solution for the problem mentioned above. In this algorithm, a hexahedral element coefficient matrix for magnetoelluric finite method is developed, which solves large-scale equations using preconditioned conjugate gradient of the first-type boundary conditions. This algorithm is verified using the homogeneous model, and the positive landform model, as well as the low resistance anomaly model.

A non-conformal finite element/finite volume scheme for the non-structured grid-based approximation of low Mach number flows; Un schema elements finis non-conformes/volumes finis pour l'approximation en maillages non-structures des ecoulements a faible nombre de Mach

Energy Technology Data Exchange (ETDEWEB)

Ansanay-Alex, G.

2009-06-17

The development of simulation codes aimed at a precise simulation of fires requires a precise approach of flame front phenomena by using very fine grids. The need to take different spatial scale into consideration leads to a local grid refinement and to a discretization with homogeneous grid for computing time and memory purposes. The author reports the approximation of the non-linear convection term, the scalar advection-diffusion in finite volumes, numerical simulations of a flow in a bent tube, of a three-dimensional laminar flame and of a low Mach number an-isotherm flow. Non conformal finite elements are also presented (Rannacher-Turek and Crouzeix-Raviart elements)

Bubble-Enriched Least-Squares Finite Element Method for Transient Advective Transport

Directory of Open Access Journals (Sweden)

Rajeev Kumar

2008-01-01

Full Text Available The least-squares finite element method (LSFEM has received increasing attention in recent years due to advantages over the Galerkin finite element method (GFEM. The method leads to a minimization problem in the L2-norm and thus results in a symmetric and positive definite matrix, even for first-order differential equations. In addition, the method contains an implicit streamline upwinding mechanism that prevents the appearance of oscillations that are characteristic of the Galerkin method. Thus, the least-squares approach does not require explicit stabilization and the associated stabilization parameters required by the Galerkin method. A new approach, the bubble enriched least-squares finite element method (BELSFEM, is presented and compared with the classical LSFEM. The BELSFEM requires a space-time element formulation and employs bubble functions in space and time to increase the accuracy of the finite element solution without degrading computational performance. We apply the BELSFEM and classical least-squares finite element methods to benchmark problems for 1D and 2D linear transport. The accuracy and performance are compared.

Finite element model updating of concrete structures based on imprecise probability

Science.gov (United States)

Biswal, S.; Ramaswamy, A.

2017-09-01

Imprecise probability based methods are developed in this study for the parameter estimation, in finite element model updating for concrete structures, when the measurements are imprecisely defined. Bayesian analysis using Metropolis Hastings algorithm for parameter estimation is generalized to incorporate the imprecision present in the prior distribution, in the likelihood function, and in the measured responses. Three different cases are considered (i) imprecision is present in the prior distribution and in the measurements only, (ii) imprecision is present in the parameters of the finite element model and in the measurement only, and (iii) imprecision is present in the prior distribution, in the parameters of the finite element model, and in the measurements. Procedures are also developed for integrating the imprecision in the parameters of the finite element model, in the finite element software Abaqus. The proposed methods are then verified against reinforced concrete beams and prestressed concrete beams tested in our laboratory as part of this study.

Shakedown analysis by finite element incremental procedures

International Nuclear Information System (INIS)

Borkowski, A.; Kleiber, M.

1979-01-01

It is a common occurence in many practical problems that external loads are variable and the exact time-dependent history of loading is unknown. Instead of it load is characterized by a given loading domain: a convex polyhedron in the n-dimensional space of load parameters. The problem is then to check whether a structure shakes down, i.e. responds elastically after a few elasto-plastic cycles, or not to a variable loading as defined above. Such check can be performed by an incremental procedure. One should reproduce incrementally a simple cyclic process which consists of proportional load paths that connect the origin of the load space with the corners of the loading domain. It was proved that if a structure shakes down to such loading history then it is able to adopt itself to an arbitrary load path contained in the loading domain. The main advantage of such approach is the possibility to use existing incremental finite-element computer codes. (orig.)

A three-dimensional cell-based smoothed finite element method for elasto-plasticity

International Nuclear Information System (INIS)

Lee, Kye Hyung; Im, Se Yong; Lim, Jae Hyuk; Sohn, Dong Woo

2015-01-01

This work is concerned with a three-dimensional cell-based smoothed finite element method for application to elastic-plastic analysis. The formulation of smoothed finite elements is extended to cover elastic-plastic deformations beyond the classical linear theory of elasticity, which has been the major application domain of smoothed finite elements. The finite strain deformations are treated with the aid of the formulation based on the hyperelastic constitutive equation. The volumetric locking originating from the nearly incompressible behavior of elastic-plastic deformations is remedied by relaxing the volumetric strain through the mean value. The comparison with the conventional finite elements demonstrates the effectiveness and accuracy of the present approach.

A three-dimensional cell-based smoothed finite element method for elasto-plasticity

Energy Technology Data Exchange (ETDEWEB)

Lee, Kye Hyung; Im, Se Yong [KAIST, Daejeon (Korea, Republic of); Lim, Jae Hyuk [KARI, Daejeon (Korea, Republic of); Sohn, Dong Woo [Korea Maritime and Ocean University, Busan (Korea, Republic of)

2015-02-15

This work is concerned with a three-dimensional cell-based smoothed finite element method for application to elastic-plastic analysis. The formulation of smoothed finite elements is extended to cover elastic-plastic deformations beyond the classical linear theory of elasticity, which has been the major application domain of smoothed finite elements. The finite strain deformations are treated with the aid of the formulation based on the hyperelastic constitutive equation. The volumetric locking originating from the nearly incompressible behavior of elastic-plastic deformations is remedied by relaxing the volumetric strain through the mean value. The comparison with the conventional finite elements demonstrates the effectiveness and accuracy of the present approach.

Some efficient Lagrangian mesh finite elements encoded in ZEPHYR for two dimensional transport calculations

International Nuclear Information System (INIS)

Mordant, Maurice.

1981-04-01

To solve a multigroup stationary neutron transport equation in two-dimensional geometries (X-Y), (R-O) or (R-Z) generally on uses discrete ordinates and rectangular meshes. The way to do it is then well known, well documented and somewhat obvious. If one needs to treat awkward geometries or distorted meshes, things are not so easy and the way to do it is no longer straightforward. We have studied this problem at Limeil Nuclear Center and as an alternative to Monte Carlo methods and code we have implemented in ZEPHYR code at least two efficient finite element solutions for Lagrangian meshes involving any kind of triangles and quadrilaterals

On Using Particle Finite Element for Hydrodynamics Problems Solving

Directory of Open Access Journals (Sweden)

E. V. Davidova

2015-01-01

Full Text Available The aim of the present research is to develop software for the Particle Finite Element Method (PFEM and its verification on the model problem of viscous incompressible flow simulation in a square cavity. The Lagrangian description of the medium motion is used: the nodes of the finite element mesh move together with the fluid that allows to consider them as particles of the medium. Mesh cells deform when in time-stepping procedure, so it is necessary to reconstruct the mesh to provide stability of the finite element numerical procedure.Meshing algorithm allows us to obtain the mesh, which satisfies the Delaunay criteria: it is called \\the possible triangles method". This algorithm is based on the well-known Fortune method of Voronoi diagram constructing for a certain set of points in the plane. The graphical representation of the possible triangles method is shown. It is suitable to use generalization of Delaunay triangulation in order to construct meshes with polygonal cells in case of multiple nodes close to be lying on the same circle.The viscous incompressible fluid flow is described by the Navier | Stokes equations and the mass conservation equation with certain initial and boundary conditions. A fractional steps method, which allows us to avoid non-physical oscillations of the pressure, provides the timestepping procedure. Using the finite element discretization and the Bubnov | Galerkin method allows us to carry out spatial discretization.For form functions calculation of finite element mesh with polygonal cells, \

Analysis of thin composite structures using an efficient hex-shell finite element

Energy Technology Data Exchange (ETDEWEB)

Shiri, Seddik [Universite Bordeaux, Pessac (France); Naceur, Hakim [Universite de valenciennes, Valenciennes (France)

2013-12-15

In this paper a general methodology for the modeling of material composite multilayered shell structures is proposed using a Hex-shell finite element modeling. The first part of the paper is devoted to the general FE formulation of the present composite 8-node Hex-shell element called SCH8, based only on displacement degrees of freedom. A particular attention is given to alleviate shear, trapezoidal and thickness locking, without resorting to the classical plane-stress assumption. The anisotropic material behavior of layered shells is modeled using a fully three dimensional elastic orthotropic material law in each layer, including the thickness stress component. Applications to laminate thick shell structures are studied to validate the methodology, and good results have been obtained in comparison with ABAQUS commercial code.

FINITE ELEMENT ANALYSIS OF ELEMENT ANALYSIS OF A FREE ...

African Journals Online (AJOL)

eobe

the stairs and to compare the finite element ana ... tual three dimensional behavior of the stair slab system. ..... due to its close relation of output with the propo .... flights. It is best not to consider any open well when .... thermodynamics of solids.

Finite element modelling of different CANDU fuel bundle types in various refuelling conditions

International Nuclear Information System (INIS)

Roman, M. R.; Ionescu, D. V.; Olteanu, G.; Florea, S.; Radut, A. C.

2016-01-01

The objective of this paper is to develop a finite element model for static strength analysis of the CANDU standard with 37 elements fuel bundle and the SEU43 with 43 elements fuel bundle design for various refuelling conditions. The computer code, ANSYS7.1, is used to simulate the axial compression in CANDU type fuel bundles subject to hydraulic drag loads, deflection of fuel elements, stresses and displacements in the end plates. Two possible situations for the fuelling machine side stops are considered in our analyses, as follows: the last fuel bundle is supported by the two side stops and a side stop can be blocked therefore, the last fuel bundle is supported by only one side stop. The results of the analyses performed are briefly presented and also illustrated in a graphical form. The finite element model developed in present study is verified against test results for endplate displacement and element bowing obtained from strength tests with fuel bundle string and fuelling machine side-stop simulators. Comparison of ANSYS model predictions with these experimental results led to a very good agreement. Despite the difference in hydraulic load between SEU43 and CANDU standard fuel bundles strings, the maximum stress in the SEU43 endplate is about the same with the maximum stress in the CANDU standard endplate. The comparative assessment reveals that SEU43 fuel bundle is able to withstand high flow rate without showing a significant geometric instability. (authors)

Finite element analysis for the impact behaviour of a cask interacting with a rigid pin

International Nuclear Information System (INIS)

Altes, J.; Geiser, H.; Voelzer, W.; Frenk, A.; Deeken, G.

1993-01-01

Full scale drop tests of casks to be licensed as type B packages according to the IAEA regulations for the safe transport of radioactive materials are expensive. Therefore efforts are being made to use computer codes for calculating the impact behaviour. But these codes have to be verified by experiments. Codes available for these calculations are for example DYNA3D and ABAQUS. In the paper results of both codes are compared. A 11 t ductile cast iron cask (type MOSAIK) without impact limiters was analysed dropping from a height of 1 m with its top onto a cylindrical steel pin. The results of the finite element calculations with both codes show good agreement. The ABAQUS results using the implicit method are in accordance with the explicit method, for which considerably shorter CPU times are noted. (author)

Finite element simulations of two rock mechanics tests

International Nuclear Information System (INIS)

Dahlke, H.J.; Lott, S.A.

1986-04-01

Rock mechanics tests are performed to determine in situ stress conditions and material properties of an underground rock mass. To design stable underground facilities for the permanent storage of high-level nuclear waste, determination of these properties and conditions is a necessary first step. However, before a test and its associated equipment can be designed, the engineer needs to know the range of expected values to be measured by the instruments. Sensitivity studies by means of finite element simulations are employed in this preliminary design phase to evaluate the pertinent parameters and their effects on the proposed measurements. The simulations, of two typical rock mechanics tests, the plate bearing test and the flat-jack test, by means of the finite element analysis, are described. The plate bearing test is used to determine the rock mass deformation modulus. The flat-jack test is used to determine the in situ stress conditions of the host rock. For the plate bearing test, two finite element models are used to simulate the classic problem of a load on an elastic half space and the actual problem of a plate bearing test in an underground tunnel of circular cross section. For the flat-jack simulation, a single finite element model is used to simulate both horizontal and vertical slots. Results will be compared to closed-form solutions available in the literature

Rolling Process Modeling Report: Finite-Element Prediction of Roll Separating Force and Rolling Defects

Energy Technology Data Exchange (ETDEWEB)

Soulami, Ayoub [Pacific Northwest National Lab. (PNNL), Richland, WA (United States); Lavender, Curt A. [Pacific Northwest National Lab. (PNNL), Richland, WA (United States); Paxton, Dean M. [Pacific Northwest National Lab. (PNNL), Richland, WA (United States); Burkes, Douglas [Pacific Northwest National Lab. (PNNL), Richland, WA (United States)

2014-04-23

Pacific Northwest National Laboratory (PNNL) has been investigating manufacturing processes for the uranium-10% molybdenum (U-10Mo) alloy plate-type fuel for the U.S. high-performance research reactors. This work supports the Convert Program of the U.S. Department of Energy’s National Nuclear Security Administration (DOE/NNSA) Global Threat Reduction Initiative. This report documents modeling results of PNNL’s efforts to perform finite-element simulations to predict roll separating forces and rolling defects. Simulations were performed using a finite-element model developed using the commercial code LS-Dyna. Simulations of the hot rolling of U-10Mo coupons encapsulated in low-carbon steel have been conducted following two different schedules. Model predictions of the roll-separation force and roll-pack thicknesses at different stages of the rolling process were compared with experimental measurements. This report discusses various attributes of the rolled coupons revealed by the model (e.g., dog-boning and thickness non-uniformity).

A simple finite element method for linear hyperbolic problems

International Nuclear Information System (INIS)

Mu, Lin; Ye, Xiu

2017-01-01

Here, we introduce a simple finite element method for solving first order hyperbolic equations with easy implementation and analysis. Our new method, with a symmetric, positive definite system, is designed to use discontinuous approximations on finite element partitions consisting of arbitrary shape of polygons/polyhedra. Error estimate is established. Extensive numerical examples are tested that demonstrate the robustness and flexibility of the method.

Finite element computation of plasma equilibria

International Nuclear Information System (INIS)

Rivier, M.

1977-01-01

The applicability of the finite element method is investigated for the numerical solution of the nonlinear Grad-Shafranov equation with free boundary for the flux function of a plasma at equilibrium. This method is based on the case of variational principles and finite dimensional subspaces whose elements are piecewise polynomial functions obtained by a Lagrange type interpolation procedure over a triangulation of the domain. Two cases of plasma pressure (exponential and quadratic including a vacuum region) were examined. In both cases the nonuniqueness of the solutions was shown in exhibiting a deeper solution in the case of exponential pressure function, and a non-constant solution for a quadratic pressure function. In order to get this ''other'' solution, two linearization methods were tested with two different constraints. Different cross sections are investigated

Flow Applications of the Least Squares Finite Element Method

Science.gov (United States)

Jiang, Bo-Nan

1998-01-01

The main thrust of the effort has been towards the development, analysis and implementation of the least-squares finite element method (LSFEM) for fluid dynamics and electromagnetics applications. In the past year, there were four major accomplishments: 1) special treatments in computational fluid dynamics and computational electromagnetics, such as upwinding, numerical dissipation, staggered grid, non-equal order elements, operator splitting and preconditioning, edge elements, and vector potential are unnecessary; 2) the analysis of the LSFEM for most partial differential equations can be based on the bounded inverse theorem; 3) the finite difference and finite volume algorithms solve only two Maxwell equations and ignore the divergence equations; and 4) the first numerical simulation of three-dimensional Marangoni-Benard convection was performed using the LSFEM.

Investigation of Shear Stud Performance in Flat Plate Using Finite Element Analysis

Directory of Open Access Journals (Sweden)

T.S. Viswanathan

2014-09-01

Full Text Available Three types of shear stud arrangement, respectively featuring an orthogonal, a radial and a critical perimeter pattern, were evaluated numerically. A numerical investigation was conducted using the finite element software ABAQUS to evaluate their ability to resist punching shear in a flat plate. The finite element analysis here is an application of the nonlinear analysis of reinforced concrete structures using three-dimensional solid finite elements. The nonlinear characteristics of concrete were achieved by employing the concrete damaged plasticity model in the finite element program. Transverse shear stress was evaluated using finite element analysis in terms of shear stress distribution for flat plate with and without shear stud reinforcement. The model predicted that shear studs placed along the critical perimeter are more effective compared to orthogonal and radial patterns.

Parallel direct solver for finite element modeling of manufacturing processes

DEFF Research Database (Denmark)

Nielsen, Chris Valentin; Martins, P.A.F.

2017-01-01

The central processing unit (CPU) time is of paramount importance in finite element modeling of manufacturing processes. Because the most significant part of the CPU time is consumed in solving the main system of equations resulting from finite element assemblies, different approaches have been...

A finite range coupled channel Born approximation code

International Nuclear Information System (INIS)

Nagel, P.; Koshel, R.D.

1978-01-01

The computer code OUKID calculates differential cross sections for direct transfer nuclear reactions in which multistep processes, arising from strongly coupled inelastic states in both the target and residual nuclei, are possible. The code is designed for heavy ion reactions where full finite range and recoil effects are important. Distorted wave functions for the elastic and inelastic scattering are calculated by solving sets of coupled differential equations using a Matrix Numerov integration procedure. These wave functions are then expanded into bases of spherical Bessel functions by the plane-wave expansion method. This approach allows the six-dimensional integrals for the transition amplitude to be reduced to products of two one-dimensional integrals. Thus, the inelastic scattering is treated in a coupled channel formalism while the transfer process is treated in a finite range born approximation formalism. (Auth.)

High-Order Curvilinear Finite Element Methods for Lagrangian Hydrodynamics [High Order Curvilinear Finite Elements for Lagrangian Hydrodynamics

Energy Technology Data Exchange (ETDEWEB)

Dobrev, Veselin A. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Kolev, Tzanio V. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Rieben, Robert N. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

2012-09-20

The numerical approximation of the Euler equations of gas dynamics in a movingLagrangian frame is at the heart of many multiphysics simulation algorithms. Here, we present a general framework for high-order Lagrangian discretization of these compressible shock hydrodynamics equations using curvilinear finite elements. This method is an extension of the approach outlined in [Dobrev et al., Internat. J. Numer. Methods Fluids, 65 (2010), pp. 1295--1310] and can be formulated for any finite dimensional approximation of the kinematic and thermodynamic fields, including generic finite elements on two- and three-dimensional meshes with triangular, quadrilateral, tetrahedral, or hexahedral zones. We discretize the kinematic variables of position and velocity using a continuous high-order basis function expansion of arbitrary polynomial degree which is obtained via a corresponding high-order parametric mapping from a standard reference element. This enables the use of curvilinear zone geometry, higher-order approximations for fields within a zone, and a pointwise definition of mass conservation which we refer to as strong mass conservation. Moreover, we discretize the internal energy using a piecewise discontinuous high-order basis function expansion which is also of arbitrary polynomial degree. This facilitates multimaterial hydrodynamics by treating material properties, such as equations of state and constitutive models, as piecewise discontinuous functions which vary within a zone. To satisfy the Rankine--Hugoniot jump conditions at a shock boundary and generate the appropriate entropy, we introduce a general tensor artificial viscosity which takes advantage of the high-order kinematic and thermodynamic information available in each zone. Finally, we apply a generic high-order time discretization process to the semidiscrete equations to develop the fully discrete numerical algorithm. Our method can be viewed as the high-order generalization of the so-called staggered

Analysis of Dynamic Fracture Parameters in Functionally Graded Material Plates with Cracks by Graded Finite Element Method and Virtual Crack Closure Technique

Directory of Open Access Journals (Sweden)

Li Ming Zhou

2016-01-01

Full Text Available Based on the finite element software ABAQUS and graded element method, we developed a dummy node fracture element, wrote the user subroutines UMAT and UEL, and solved the energy release rate component of functionally graded material (FGM plates with cracks. An interface element tailored for the virtual crack closure technique (VCCT was applied. Fixed cracks and moving cracks under dynamic loads were simulated. The results were compared to other VCCT-based analyses. With the implementation of a crack speed function within the element, it can be easily expanded to the cases of varying crack velocities, without convergence difficulty for all cases. Neither singular element nor collapsed element was required. Therefore, due to its simplicity, the VCCT interface element is a potential tool for engineers to conduct dynamic fracture analysis in conjunction with commercial finite element analysis codes.

Engineering computation of structures the finite element method

CERN Document Server

Neto, Maria Augusta; Roseiro, Luis; Cirne, José; Leal, Rogério

2015-01-01

This book presents theories and the main useful techniques of the Finite Element Method (FEM), with an introduction to FEM and many case studies of its use in engineering practice. It supports engineers and students to solve primarily linear problems in mechanical engineering, with a main focus on static and dynamic structural problems. Readers of this text are encouraged to discover the proper relationship between theory and practice, within the finite element method: Practice without theory is blind, but theory without practice is sterile. Beginning with elasticity basic concepts and the classical theories of stressed materials, the work goes on to apply the relationship between forces, displacements, stresses and strains on the process of modeling, simulating and designing engineered technical systems. Chapters discuss the finite element equations for static, eigenvalue analysis, as well as transient analyses. Students and practitioners using commercial FEM software will find this book very helpful. It us...

Modelling drawbeads with finite elements and verification

NARCIS (Netherlands)

Carleer, B.D.; Carleer, B.D.; Vreede, P.T.; Vreede, P.T.; Louwes, M.F.M.; Louwes, M.F.M.; Huetink, Han

1994-01-01

Drawbeads are commonly used in deep drawing processes to control the flow of the blank during the forming operation. In finite element simulations of deep drawing the drawbead geometries are seldom included because of the small radii; because of these small radii a very large number of elements is

Finite element prediction of the swift effect based on Taylor-type polycrystal plasticity models

OpenAIRE

Duchene, Laurent; Delannay, L.; Habraken, Anne

2004-01-01

This paper describes the main concepts of the stress-strain interpolation model that has been implemented in the non-linear finite element code Lagamine. This model consists in a local description of the yield locus based on the texture of the material through the full constraints Taylor’s model. The prediction of the Swift effect is investigated: the influence of the texture evolution is shown up. The LAMEL model is also investigated for the Swift effect prediction. Peer reviewed

SAFE: A computer code for the steady-state and transient thermal analysis of LMR fuel elements

International Nuclear Information System (INIS)

Hayes, S.L.

1993-12-01

SAFE is a computer code developed for both the steady-state and transient thermal analysis of single LMR fuel elements. The code employs a two-dimensional control-volume based finite difference methodology with fully implicit time marching to calculate the temperatures throughout a fuel element and its associated coolant channel for both the steady-state and transient events. The code makes no structural calculations or predictions whatsoever. It does, however, accept as input structural parameters within the fuel such as the distributions of porosity and fuel composition, as well as heat generation, to allow a thermal analysis to be performed on a user-specified fuel structure. The code was developed with ease of use in mind. An interactive input file generator and material property correlations internal to the code are available to expedite analyses using SAFE. This report serves as a complete design description of the code as well as a user's manual. A sample calculation made with SAFE is included to highlight some of the code's features. Complete input and output files for the sample problem are provided

Documentation of SPECTROM-55: A finite element thermohydrogeological analysis program

International Nuclear Information System (INIS)

Osnes, J.D.; Ratigan, J.L.; Loken, M.C.; Parrish, D.K.

1985-12-01

SPECTROM-55 is a finite element computer program developed by RE/SPEC Inc. for analyses of coupled heat and fluid transfer through fully saturated porous media. The theoretical basis for the mathematical model, the implementation of the mathematical model into the computer code, the verification and validation efforts with the computer code, and the code support and continuing documentation are described in this document. The program is especially suited for analyses of the regional hydrogeology in the vicinity of a heat-generating nuclear waste repository. These applications typically involve forced and free convection in a ground-water flow system. The program provides transient or steady-state temperatures, pressures, and fluid velocities resulting from the application of a variety of initial and boundary conditions to bodies with complicated shapes. The boundary conditions include constant heat and fluid fluxes, convective heat transfer, constant temperature, and constant pressure. Initial temperatures and pressures can be specified. Composite systems of anisotropic materials, such as geologic strata, can be defined in either planar or axisymmetric configurations. Constant or varying volumetric heat generation, such as decaying heat generation from radioactive waste, can be specified

Finite rate chemistry for USA-series codes - rormulation and applications

International Nuclear Information System (INIS)

Palaniswamy, S.; Chakravarthy, S.R.; Ota, D.K.

1989-01-01

The USA-series of CFD codes are based on unified solution algorithms including explicit and implicit formulations, factorization and relaxation approaches, time marching and space marching methodologies, etc., in order to be able to solve a very wide class of CFD problems using a single framework. Euler or Navier-Stokes equations are solved using a finite-volume treatment with upwind Total Variation Diminishing discretization for the inviscid terms. Recently, these codes have been enlarged to also unify different aerothermodynamic options (perfect gas, real gas including equilibrium and nonequlibrium chemistry). This paper describes aspects of the finite-rate-chemistry capability. 27 references

A finite element solution method for quadrics parallel computer

International Nuclear Information System (INIS)

Zucchini, A.

1996-08-01

A distributed preconditioned conjugate gradient method for finite element analysis has been developed and implemented on a parallel SIMD Quadrics computer. The main characteristic of the method is that it does not require any actual assembling of all element equations in a global system. The physical domain of the problem is partitioned in cells of n p finite elements and each cell element is assigned to a different node of an n p -processors machine. Element stiffness matrices are stored in the data memory of the assigned processing node and the solution process is completely executed in parallel at element level. Inter-element and therefore inter-processor communications are required once per iteration to perform local sums of vector quantities between neighbouring elements. A prototype implementation has been tested on an 8-nodes Quadrics machine in a simple 2D benchmark problem

ABAQUS-EPGEN: a general-purpose finite element code. Volume 3. Example problems manual

International Nuclear Information System (INIS)

Hibbitt, H.D.; Karlsson, B.I.; Sorensen, E.P.

1983-03-01

This volume is the Example and Verification Problems Manual for ABAQUS/EPGEN. Companion volumes are the User's, Theory and Systems Manuals. This volume contains two major parts. The bulk of the manual (Sections 1-8) contains worked examples that are discussed in detail, while Appendix A documents a large set of basic verification cases that provide the fundamental check of the elements in the code. The examples in Sections 1-8 illustrate and verify significant aspects of the program's capability. Most of these problems provide verification, but they have also been chosen to allow discussion of modeling and analysis techniques. Appendix A contains basic verification cases. Each of these cases verifies one element in the program's library. The verification consists of applying all possible load or flux types (including thermal loading of stress elements), and all possible foundation or film/radiation conditions, and checking the resulting force and stress solutions or flux and temperature results. This manual provides program verification. All of the problems described in the manual are run and the results checked, for each release of the program, and these verification results are made available

Evaluation of Concrete Cylinder Tests Using Finite Elements

DEFF Research Database (Denmark)

Saabye Ottosen, Niels

1984-01-01

Nonlinear axisymmetric finite element analyses are performed on the uniaxial compressive test of concrete cylinders. The models include thick steel loading plates, and cylinders with height‐to‐diameter ratios (h/d) ranging from 1‐3 are treated. A simple constitutive model of the concrete is emplo......Nonlinear axisymmetric finite element analyses are performed on the uniaxial compressive test of concrete cylinders. The models include thick steel loading plates, and cylinders with height‐to‐diameter ratios (h/d) ranging from 1‐3 are treated. A simple constitutive model of the concrete...... uniaxial strength the use of geometrically matched loading plates seems to be advantageous. Finally, it is observed that for variations of the element size within limits otherwise required to obtain a realistic analysis, the results are insensitive to the element size....

Finite element analysis of structures through unified formulation

CERN Document Server

Carrera, Erasmo; Petrolo, Marco; Zappino, Enrico

2014-01-01

The finite element method (FEM) is a computational tool widely used to design and analyse  complex structures. Currently, there are a number of different approaches to analysis using the FEM that vary according to the type of structure being analysed: beams and plates may use 1D or 2D approaches, shells and solids 2D or 3D approaches, and methods that work for one structure are typically not optimized to work for another. Finite Element Analysis of Structures Through Unified Formulation deals with the FEM used for the analysis of the mechanics of structures in the case of linear elasticity. The novelty of this book is that the finite elements (FEs) are formulated on the basis of a class of theories of structures known as the Carrera Unified Formulation (CUF). It formulates 1D, 2D and 3D FEs on the basis of the same ''fundamental nucleus'' that comes from geometrical relations and Hooke''s law, and presents both 1D and 2D refined FEs that only have displacement variables as in 3D elements. It also covers 1D...

Exponentially-convergent Monte Carlo via finite-element trial spaces

International Nuclear Information System (INIS)

Morel, Jim E.; Tooley, Jared P.; Blamer, Brandon J.

2011-01-01

Exponentially-Convergent Monte Carlo (ECMC) methods, also known as adaptive Monte Carlo and residual Monte Carlo methods, were the subject of intense research over a decade ago, but they never became practical for solving the realistic problems. We believe that the failure of previous efforts may be related to the choice of trial spaces that were global and thus highly oscillatory. As an alternative, we consider finite-element trial spaces, which have the ability to treat fully realistic problems. As a first step towards more general methods, we apply piecewise-linear trial spaces to the spatially-continuous two-stream transport equation. Using this approach, we achieve exponential convergence and computationally demonstrate several fundamental properties of finite-element based ECMC methods. Finally, our results indicate that the finite-element approach clearly deserves further investigation. (author)

Finite element formulation for thermal stress analysis of thin reactor structures

International Nuclear Information System (INIS)

Kulak, R.F.; Kennedy, J.M.; Belytschko, T.B.

1978-01-01

This paper describes the formulation of a finite-element procedure for the thermal stress analysis of thin wall reactor components. A general temperature-dependent constituent relationship is derived from a Gibbs potential function and a temperature-dependent yield surface. This form of constitutive relationship is applicable to problems of small strain. A similar form of a hypoelastic-plastic type is also developed for large strains. The variation of the yield surface with temperature is based upon a temperature-dependent, work-hardening model. The model uses a temperature-equivalent stress-plastic strain diagram which is generated from isothermal unaxial stress-strain data. The above constitutive relationships are incorporated into two computer codes and a previously developed numerical algorithm is used with minor modifications. A set of problems is presented validating the thermal analysis capability of the computer codes to a variety of problem types. (Auth.)

Some studies on weld bead geometries for laser spot welding process using finite element analysis

International Nuclear Information System (INIS)

Siva Shanmugam, N.; Buvanashekaran, G.; Sankaranarayanasamy, K.

2012-01-01

Highlights: → In this study, a 2 kW Nd:YAG laser welding system is used to conduct laser spot welding trials. → The size and shape of the laser spot weld is predicted using finite element simulation. → The heat input is assumed to be a three-dimensional conical Gaussian heat source. → The result highlights the effect of beam incident angle on laser spot welds. → The achieved results of numerical simulation are almost identical with a real weldment. -- Abstract: Nd:YAG laser beam welding is a high power density welding process which has the capability to focus the beam to a very small spot diameter of about 0.4 mm. It has favorable characteristics namely, low heat input, narrow heat affected zone and lower distortions, as compared to conventional welding processes. In this study, finite element method (FEM) is applied for predicting the weld bead geometry i.e. bead length (BL), bead width (BW) and depth of penetration (DP) in laser spot welding of AISI 304 stainless steel sheet of thickness 2.5 mm. The input parameters of laser spot welding such as beam power, incident angle of the beam and beam exposure time are varied for conducting experimental trials and numerical simulations. Temperature-dependent thermal properties of AISI 304 stainless steel, the effect of latent heat of fusion, and the convective and radiative aspects of boundary conditions are considered while developing the finite element model. The heat input to the developed model is assumed to be a three-dimensional conical Gaussian heat source. Finite-element simulations of laser spot welding were carried out by using Ansys Parametric Design Language (APDL) available in finite-element code, ANSYS. The results of the numerical analysis provide the shape of the weld beads for different ranges of laser input parameters that are subsequently compared with the results obtained through experimentation and it is found that they are in good agreement.

Automatic mesh generation for finite element calculations in the case of thermal loads

International Nuclear Information System (INIS)

Cords, H.; Zimmermann, R.

1975-01-01

The presentation describes a method to generate finite element nodal point networks on the basis of isothermals and flux lines. Such a mesh provides a relatively fine partitioning at regions where pronounced temperature variations exist. In case of entirely thermal loads a net of this kind is advantageous since the refinement is provided at exactly those locations where high stress levels are expected. In the present contribution the method was employed to analyze the structural behavior of a nuclear fuel element under operating conditions. The graphite block fuel elements for high temperature reactors are of prismatic shape with a large number of parallel bores in the axial direction. Some of these bores are open at both ends and cooling is effected by helium flowing through. Blind holes contain the fuel as compacts or cartridges. The basic temperature distribution in a horizontal section of the block was obtained by the boundary point least squares method which yields analytical expressions for both temperature and thermal flux. The corresponding computer code was presented at an earlier SMiRT conference. The method is particularly useful for regular arrays of heat sources and sinks as encountered in heat exchanger problems. The generated mesh matches the requirements of a subsequent structural analysis with finite elements provided there are no other than thermal loads

Analysis of Tube Drawing Process – A Finite Element Approach ...

African Journals Online (AJOL)

In this paper the effect of die semi angle on drawing load in cold tube drawing has been investigated numerically using the finite element method. The equation governing the stress distribution was derived and solved using Galerkin finite element method. An isoparametric formulation for the governing equation was utilized ...

A finite element for plates and shells

International Nuclear Information System (INIS)

Muller, A.; Feijoo, R.A.; Bevilacqua, L.

1981-08-01

A simple triangular finite element for plates and shells, is presented. Since the rotation fields are assumed independent of the displacement fields, the element allows one to solve thick shells problems. In the limit for thin shell, the Kirchoff-Love hypothesis is automatically satisfied, thus enlarging its range of application. (Author) [pt

COMPUTER EXPERIMENTS WITH FINITE ELEMENTS OF HIGHER ORDER

Directory of Open Access Journals (Sweden)

Khomchenko A.

2017-12-01

Full Text Available The paper deals with the problem of constructing the basic functions of a quadrilateral finite element of the fifth order by the means of the computer algebra system Maple. The Lagrangian approximation of such a finite element contains 36 nodes: 20 nodes perimeter and 16 internal nodes. Alternative models with reduced number of internal nodes are considered. Graphs of basic functions and cognitive portraits of lines of zero level are presented. The work is aimed at studying the possibilities of using modern information technologies in the teaching of individual mathematical disciplines.

Fourier analysis of finite element preconditioned collocation schemes

Science.gov (United States)

Deville, Michel O.; Mund, Ernest H.

1990-01-01

The spectrum of the iteration operator of some finite element preconditioned Fourier collocation schemes is investigated. The first part of the paper analyses one-dimensional elliptic and hyperbolic model problems and the advection-diffusion equation. Analytical expressions of the eigenvalues are obtained with use of symbolic computation. The second part of the paper considers the set of one-dimensional differential equations resulting from Fourier analysis (in the tranverse direction) of the 2-D Stokes problem. All results agree with previous conclusions on the numerical efficiency of finite element preconditioning schemes.

Stochastic Finite Elements in Reliability-Based Structural Optimization

DEFF Research Database (Denmark)

Sørensen, John Dalsgaard; Engelund, S.

Application of stochastic finite elements in structural optimization is considered. It is shown how stochastic fields modelling e.g. the modulus of elasticity can be discretized in stochastic variables and how a sensitivity analysis of the reliability of a structural system with respect to optimi......Application of stochastic finite elements in structural optimization is considered. It is shown how stochastic fields modelling e.g. the modulus of elasticity can be discretized in stochastic variables and how a sensitivity analysis of the reliability of a structural system with respect...

Mechanical strength calculation of the disk type windings with elastic couplings by the finite element method

International Nuclear Information System (INIS)

Sivkova, G.N.; Spirchenko, Yu.V.; Chvartatskij, P.V.

1981-01-01

Stressed-deformed state of toroidal field coils of the disc type with elastic couplings of the tokamaks has been investigated with provision for the effect of the central core pliability by means of the two-dimensional version of the finite element method. Numerical solution of the finite element method is performed by means of the ES 1040 computer according to the computer code permitting taking account of boundary conditions of elastic support. The calculation has been performed using as the example the project of T-20 facility coil of the disc type. Consideration of pliability of the central core of the facility inductor is accomplished by the introduction of additional rigidities to the complete matrix of rigidity. Scheme of the structure distretization includes 141 units, 211 elements. The accuracy of solution depends on the reduction accuracy of the volume load to unit forces and on the number of finite elements. Analysis of the solution convergence is performed by the comparison of solutions obtained for three different schemes of the disk discretization without regard for the inductor pliability. The comparative analysis of the results shows that transfer epures for all the three discretization versions practically coincide and stresses differ not more than by 10%. On the whole the above investigation has demonstrated good convergence of the problem solution [ru

Comparison of finite element J-integral evaluations for the blunt crack model and the sharp crack model

International Nuclear Information System (INIS)

Pan, Y.C.; Kennedy, J.M.

1983-01-01

In assessing the safety of a liquid metal fast breeder reactor (LMFBR), a major concern is that of hot sodium coming into contact with either unprotected concrete or steel-lined concrete equipment cells and containment structures. An aspect of this is the potential of concrete cracking which would significantly influence the safety assessment. Concrete cracking in finite element analysis can be modeled as a blunt crack in which the crack is assumed to be uniformly distributed throughout the area of the element. A blunt crack model based on the energy release rate and the effective strength concepts which was insensitive to the element size was presented by Bazant and Cedolin. Some difficulties were encountered in incorporating their approach into a general purpose finite element code. An approach based on the J-integral to circumvent some of the difficulties was proposed by Pan, Marchertas, and Kennedy. Alternatively, cracking can also be modeled as a sharp crack where the crack surface is treated as the boundary of the finite element mesh. The sharp crack model is adopted by most researchers and its J-integral has been well established. It is desirable to establish the correlation between the J-integrals, or the energy release rates, for the blunt crack model and the sharp crack model so that data obtained from one model can be used on the other

Finite-element solutions of the AER-2 rod ejection benchmark by CRONOS

International Nuclear Information System (INIS)

Kolev, N.P.; Lenain, R.; Fedon-Magnaud, C.

2001-01-01

The finite-element option in CRONOS was used to analyse the AER-2 rod-ejection benchmark for WWER-440. The objective is to obtain spatially converged solutions by means of node subdivision and approximation refinement. This paper presents the first phase of analysis dealing with the initial and just-ejected states used for calculation of the initial reactivity. Fine-mesh and extrapolated to zero mesh size solutions were obtained and verified by comparison to MAG code solutions. These differences provide potential for large deviations in the transient results and deserve further attention in reactor safety analysis (Authors)

Modelling bucket excavation by finite element

Science.gov (United States)

Pecingina, O. M.

2015-11-01

Changes in geological components of the layers from lignite pits have an impact on the sustainability of the cup path elements and under the action of excavation force appear efforts leading to deformation of the entire assembly. Application of finite element method in the optimization of components leads to economic growth, to increase the reliability and durability of the studied machine parts thus the machine. It is obvious usefulness of knowledge the state of mechanical tensions that the designed piece or the assembly not to break under the action of tensions that must cope during operation. In the course of excavation work on all bucket cutting force components, the first coming into contact with the material being excavated cutting edge. Therefore in the study with finite element analysis is retained only cutting edge. To study the field of stress and strain on the cutting edge will be created geometric patterns for each type of cup this will be subject to static analysis. The geometric design retains the cutting edge shape and on this on the tooth cassette location will apply an areal force on the abutment tooth. The cutting edge real pattern is subjected to finite element study for the worst case of rock cutting by symmetrical and asymmetrical cups whose profile is different. The purpose of this paper is to determine the displacement and tensions field for both profiles considering the maximum force applied on the cutting edge and the depth of the cutting is equal with the width of the cutting edge of the tooth. It will consider the worst case when on the structure will act both the tangential force and radial force on the bucket profile. For determination of stress and strain field on the form design of cutting edge profile will apply maximum force assuming uniform distribution and on the edge surface force will apply a radial force. After geometric patterns discretization on the cutting knives and determining stress field, can be seen that at the

Introduction to the Explicit Finite Element Method for Nonlinear Transient Dynamics

CERN Document Server

Wu, Shen R

2012-01-01

A systematic introduction to the theories and formulations of the explicit finite element method As numerical technology continues to grow and evolve with industrial applications, understanding the explicit finite element method has become increasingly important, particularly in the areas of crashworthiness, metal forming, and impact engineering. Introduction to the Explicit FiniteElement Method for Nonlinear Transient Dynamics is the first book to address specifically what is now accepted as the most successful numerical tool for nonlinear transient dynamics. The book aids readers in master

A study on the improvement of shape optimization associated with the modification of a finite element

International Nuclear Information System (INIS)

Sung, Jin Il; Yoo, Jeong Hoon

2002-01-01

In this paper, we investigate the effect and the importance of the accuracy of finite element analysis in the shape optimization based on the finite element method and improve the existing finite element which has inaccuracy in some cases. And then, the shape optimization is performed by using the improved finite element. One of the main stream to improve finite element is the prevention of locking phenomenon. In case of bending dominant problems, finite element solutions cannot be reliable because of shear locking phenomenon. In the process of shape optimization, the mesh distortion is large due to the change of the structure outline. So, we have to raise the accuracy of finite element analysis for the large mesh distortion. We cannot guarantee the accurate result unless the finite element itself is accurate or the finite elements are remeshed. So, we approach to more accurate shape optimization to diminish these inaccuracies by improving the existing finite element. The shape optimization using the modified finite element is applied to a two and three dimensional simple beam. Results show that the modified finite element has improved the optimization results

Finite element analysis of car hood for impact test by using ...

African Journals Online (AJOL)

Finite element analysis of car hood for impact test by using solidworks software ... high safety and at the same time can be built according to market demands. ... Keywords: finite element analysis; impact test; Solidworks; automation, car hood.

A Full-wave Model for Wave Propagation and Dissipation in the Inner Magnetosphere Using the Finite Element Method

International Nuclear Information System (INIS)

Valeo, Ernest; Johnson, Jay R.; Kim, Eun-Hwa; Phillips, Cynthia

2012-01-01

A wide variety of plasma waves play an important role in the energization and loss of particles in the inner magnetosphere. Our ability to understand and model wave-particle interactions in this region requires improved knowledge of the spatial distribution and properties of these waves as well as improved understanding of how the waves depend on changes in solar wind forcing and/or geomagnetic activity. To this end, we have developed a two-dimensional, finite element code that solves the full wave equations in global magnetospheric geometry. The code describes three-dimensional wave structure including mode conversion when ULF, EMIC, and whistler waves are launched in a two-dimensional axisymmetric background plasma with general magnetic field topology. We illustrate the capabilities of the code by examining the role of plasmaspheric plumes on magnetosonic wave propagation; mode conversion at the ion-ion and Alfven resonances resulting from external, solar wind compressions; and wave structure and mode conversion of electromagnetic ion cyclotron waves launched in the equatorial magnetosphere, which propagate along the magnetic field lines toward the ionosphere. We also discuss advantages of the finite element method for resolving resonant structures, and how the model may be adapted to include nonlocal kinetic effects.

Finite element formulation for dynamics of planar flexible multi-beam system

International Nuclear Information System (INIS)

Liu Zhuyong; Hong Jiazhen; Liu Jinyang

2009-01-01

In some previous geometric nonlinear finite element formulations, due to the use of axial displacement, the contribution of all the elements lying between the reference node of zero axial displacement and the element to the foreshortening effect should be taken into account. In this paper, a finite element formulation is proposed based on geometric nonlinear elastic theory and finite element technique. The coupling deformation terms of an arbitrary point only relate to the nodal coordinates of the element at which the point is located. Based on Hamilton principle, dynamic equations of elastic beams undergoing large overall motions are derived. To investigate the effect of coupling deformation terms on system dynamic characters and reduce the dynamic equations, a complete dynamic model and three reduced models of hub-beam are prospected. When the Cartesian deformation coordinates are adopted, the results indicate that the terms related to the coupling deformation in the inertia forces of dynamic equations have small effect on system dynamic behavior and may be neglected, whereas the terms related to coupling deformation in the elastic forces are important for system dynamic behavior and should be considered in dynamic equation. Numerical examples of the rotating beam and flexible beam system are carried out to demonstrate the accuracy and validity of this dynamic model. Furthermore, it is shown that a small number of finite elements are needed to obtain a stable solution using the present coupling finite element formulation

Isogeometric finite element analysis of poroelasticity

NARCIS (Netherlands)

Irzal, F.; Remmers, J.J.C.; Verhoosel, C.V.; Borst, de R.

2013-01-01

We present an alternative numerical approach for predicting the behaviour of a deformable fluid-saturated porous medium. The conventional finite element technology is replaced by isogeometric analysis that uses non-uniform rational B-splines. The ability of these functions to provide higher-order

Fast finite elements for surgery simulation

DEFF Research Database (Denmark)

Bro-Nielsen, Morten

1997-01-01

This paper discusses volumetric deformable models for modeling human body parts and organs in surgery simulation systems. These models are built using finite element models for linear elastic materials. To achieve real-time response condensation has been applied to the system stiffness matrix...

Simplicial Finite Elements in Higher Dimensions

Czech Academy of Sciences Publication Activity Database

Brandts, J.; Korotov, S.; Křížek, Michal

2007-01-01

Roč. 52, č. 3 (2007), s. 251-265 ISSN 0862-7940 R&D Projects: GA ČR GA201/04/1503 Institutional research plan: CEZ:AV0Z10190503 Keywords : n-simplex * finite element method * superconvergence Subject RIV: BA - General Mathematics

Investigating ASME allowable loads with finite element analyses

International Nuclear Information System (INIS)

Mattar Neto, Miguel; Bezerra, Luciano M.; Miranda, Carlos A. de J.; Cruz, Julio R.B.

1995-01-01

The evaluation of nuclear components using finite element analysis (FEA) does not generally fall into the shell type verification adopted by the ASME Code. Consequently, the demonstration that the modes of failure are avoided sometimes is not straightforward. Allowable limits, developed by limit load theory, require the computation of shell membrane and bending stresses. How to calculate these stresses from FEA is not necessarily self-evident. One approach to be considered is to develop recommendations in a case-by-case basis for the most common pressure vessel geometries and loads based on comparison between the results of elastic and plastic FEA. In this paper, FE analyses of common 2D and complex 3D geometries are examined and discussed. It will be clear that in the cases studied, stress separation and categorization are not self-evident and simple tasks to undertake. Certain unclear recommendations of ASME Code can lead the stress analyst to non conservative designs as will be demonstrated in this paper. At the endo of this paper, taking into account comparison between elastic and elastic-plastic FE results from ANSYS some observations, suggestions and conclusions about the degree of conservatism of the ASME recommendations will be addressed. (author)

Finite element analysis of human joints

Energy Technology Data Exchange (ETDEWEB)

Bossart, P.L.; Hollerbach, K.

1996-09-01

Our work focuses on the development of finite element models (FEMs) that describe the biomechanics of human joints. Finite element modeling is becoming a standard tool in industrial applications. In highly complex problems such as those found in biomechanics research, however, the full potential of FEMs is just beginning to be explored, due to the absence of precise, high resolution medical data and the difficulties encountered in converting these enormous datasets into a form that is usable in FEMs. With increasing computing speed and memory available, it is now feasible to address these challenges. We address the first by acquiring data with a high resolution C-ray CT scanner and the latter by developing semi-automated method for generating the volumetric meshes used in the FEM. Issues related to tomographic reconstruction, volume segmentation, the use of extracted surfaces to generate volumetric hexahedral meshes, and applications of the FEM are described.

Finite element analysis of human joints

International Nuclear Information System (INIS)

Bossart, P.L.; Hollerbach, K.

1996-09-01

Our work focuses on the development of finite element models (FEMs) that describe the biomechanics of human joints. Finite element modeling is becoming a standard tool in industrial applications. In highly complex problems such as those found in biomechanics research, however, the full potential of FEMs is just beginning to be explored, due to the absence of precise, high resolution medical data and the difficulties encountered in converting these enormous datasets into a form that is usable in FEMs. With increasing computing speed and memory available, it is now feasible to address these challenges. We address the first by acquiring data with a high resolution C-ray CT scanner and the latter by developing semi-automated method for generating the volumetric meshes used in the FEM. Issues related to tomographic reconstruction, volume segmentation, the use of extracted surfaces to generate volumetric hexahedral meshes, and applications of the FEM are described

Finite-element blunt-crack propagation: a modified J-integral approach

International Nuclear Information System (INIS)

Pan, Y.C.; Marchertas, A.H.; Kennedy, J.M.

1983-01-01

In assessing the safety of a liquid metal fast breeder reactor (LMFBR), a major concern is the behavior of concrete structures subjected to high temperatures. The potential of concrete cracking is an important parameter which could significantly influence the safety assessment of thermally attacked concrete. A new modified J-integral approach for the blunt crack model has been derived to provide a general procedure to accurately predict the direction of crack growth. This formulation has been incorporated into the coupled heat transfer-stress analysis finite element code TEMP-STRESS. A description of the formulation is presented in this paper. Results for the problems of a Mode I and mixed mode crack in a plate using regular and slanted meshes subjected to uniaxial and shear loading are presented

Finite element method for computational fluid dynamics with any type of elements; Finite Element Methode zur numerischen Stroemungsberechnung mit beliebigen Elementen

Energy Technology Data Exchange (ETDEWEB)

Steibler, P.

2000-07-01

The unsteady, turbulent flow is to be calculated in a complex geometry. For this purpose a stabilized finite element formulation in which the same functions for velocity and pressure are used is developed. Thus the process remains independent of the type of elements. This simplifies the application. Above all, it is easier to deal with the boundary conditions. The independency from the elements is also achieved by the extended uzawa-algorithm which uses quadratic functions for velocity and an element-constant pressure. This method is also programmed. In order to produce the unstructured grids, an algorithm is implemented which produces meshes consisting of triangular and tetrahedral elements with flow-dependent adaptation. With standard geometries both calculation methods are compared with results. Finally the flow in a draft tube of a Kaplan turbine is calculated and compared with results from model tests. (orig.) [German] Die instationaere, turbulente Stroemung in einer komplexen Geometrie soll berechnet werden. Dazu wird eine Stabilisierte Finite Element Formulierung entwickelt, bei der die gleichen Ansatzfunktionen fuer Geschwindigkeiten und Druck verwendet werden. Das Verfahren wird damit unabhaengig von der Form der Elemente. Dies vereinfacht die Anwendung. Vor allem wird der Umgang mit den Randbedingungen erleichert. Die Elementunabhaengigkeit erreicht man auch mit dem erweiterten Uzawa-Algorithmus, welcher quadratische Ansatzfunktionen fuer die Geschwindigkeiten und elementweisen konstanten Druck verwendet. Dieses Verfahren wird ebenso implementiert. Zur Erstellung der unstrukturierten Gitter wird ein Algorithmus erzeugt, der Netze aus Dreiecks- und Tetraederelementen erstellt, welche stroemungsabhaengige Groessen besitzen koennen. Anhand einiger Standardgeometrien werden die beiden Berechnungsmethoden mit Ergebnissen aus der Literatur verglichen. Als praxisrelevantes Beispiel wird abschliessend die Stroemung in einem Saugrohr einer Kaplanturbine berechnet

A Note on Symplectic, Multisymplectic Scheme in Finite Element Method

Institute of Scientific and Technical Information of China (English)

GUO Han-Ying; JI Xiao-Mei; LI Yu-Qi; WU Ke

2001-01-01

We find that with uniform mesh, the numerical schemes derived from finite element method can keep a preserved symplectic structure in one-dimensional case and a preserved multisymplectic structure in two-dimensional case respectively. These results are in fact the intrinsic reason why the numerical experiments show that such finite element algorithms are accurate in practice.``

A Novel Polygonal Finite Element Method: Virtual Node Method

Science.gov (United States)

Tang, X. H.; Zheng, C.; Zhang, J. H.

2010-05-01

Polygonal finite element method (PFEM), which can construct shape functions on polygonal elements, provides greater flexibility in mesh generation. However, the non-polynomial form of traditional PFEM, such as Wachspress method and Mean Value method, leads to inexact numerical integration. Since the integration technique for non-polynomial functions is immature. To overcome this shortcoming, a great number of integration points have to be used to obtain sufficiently exact results, which increases computational cost. In this paper, a novel polygonal finite element method is proposed and called as virtual node method (VNM). The features of present method can be list as: (1) It is a PFEM with polynomial form. Thereby, Hammer integral and Gauss integral can be naturally used to obtain exact numerical integration; (2) Shape functions of VNM satisfy all the requirements of finite element method. To test the performance of VNM, intensive numerical tests are carried out. It found that, in standard patch test, VNM can achieve significantly better results than Wachspress method and Mean Value method. Moreover, it is observed that VNM can achieve better results than triangular 3-node elements in the accuracy test.

Use of simple finite elements for mechanical systems impact analysis based on stereomechanics, stress wave propagation, and energy method approaches

International Nuclear Information System (INIS)

McCoy, Michael L.; Moradi, Rasoul; Lankarani, Hamid M.

2011-01-01

This paper examines the effectiveness of analyzing impact events in mechanical systems for design purposes using simple or low ordered finite elements. Traditional impact dynamics analyses of mechanical systems namely stereomechanics, energy method, stress-wave propagation and contact mechanics approaches are limited to very simplified geometries and provide basic analyses in making predictions and understanding the dominant features of the impact in a mechanical system. In engineering practice, impacted systems present a complexity of geometry, stiffness, mass distributions, contact areas and impact angles that are impossible to analyze and design with the traditional impact dynamics methods. In real cases, the effective tool is the finite element (FE) method. The high-end FEA codes though may be not available for typical engineer/designer. This paper provides information on whether impact events of mechanical systems can be successfully modeled using simple or low-order finite elements. FEA models using simple elements are benchmarked against theoretical impact problems and published experimental impact results. As a case study, an FE model using simple plastic beam elements is further tested to predict stresses and deflections in an experimental structural impact

Preliminary Formulation of Finite Element Solution for the 1-D, 1-G Time Dependent Neutron Diffusion Equation without Consideration about Delay Neutron

Energy Technology Data Exchange (ETDEWEB)

Ryu, Eun Hyun; Song, Yong Mann; Park, Joo Hwan [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of)

2013-05-15

If time-dependent equation is solved with the FEM, the limitation of the input geometry will disappear. It has often been pointed out that the numerical methods implemented in the RFSP code are not state-of-the-art. Although an acceleration method such as the Coarse Mesh Finite Difference (CMFD) for Finite Difference Method (FDM) does not exist for the FEM, one should keep in mind that the number of time steps for the transient simulation is not large. The rigorous formulation in this study will richen the theoretical basis of the FEM and lead to an extension of the dynamics code to deal with a more complicated problem. In this study, the formulation for the 1-D, 1-G Time Dependent Neutron Diffusion Equation (TDNDE) without consideration of the delay neutron will first be done. A problem including one multiplying medium will be solved. Also several conclusions from a comparison between the numerical and analytic solutions, a comparison between solutions with various element orders, and a comparison between solutions with different time differencing will be made to be certain about the formulation and FEM solution. By investigating various cases with different values of albedo, theta, and the order of elements, it can be concluded that the finite element solution is agree well with the analytic solution. The higher the element order used, the higher the accuracy improvements are obtained.

Comparison of 3-D finite elements for incompressible fluid flow

International Nuclear Information System (INIS)

Robichaud, M.; Tanguy, P.A.

1985-01-01

In recent years, the finite element method applied to the solution of incompressible fluid flow has been in constant evolution. In the present state-of-the-art, 2-D problems are solved routinely and reliable results are obtained at a reasonable cost. In 3-D the finite element method is still undergoing active research and many methods have been proposed to solve the Navier-Stokes equations at 'low cost'. These methods have in common the choice of the element which has a trilinear velocity and a discontinuous constant pressure (Q1-PO). The prohibitive cost of 3-D finite element method in fluid flow is the reason for this choice: the Q1-PO is the simplest and the cheapest 3-D element. However, as mentioned in (5) and (6), it generates 'spurious' pressure modes phenomenon called checkerboarding. On regular mesh these spurious modes can be filtered but on distorted mesh the pressure solution is meaningless. (author)

Finite element analysis of an inflatable torus considering air mass structural element

Science.gov (United States)

Gajbhiye, S. C.; Upadhyay, S. H.; Harsha, S. P.

2014-01-01

Inflatable structures, also known as gossamer structures, are at high boom in the current space technology due to their low mass and compact size comparing to the traditional spacecraft designing. Internal pressure becomes the major source of strength and rigidity, essentially stiffen the structure. However, inflatable space based membrane structure are at high risk to the vibration disturbance due to their low structural stiffness and material damping. Hence, the vibration modes of the structure should be known to a high degree of accuracy in order to provide better control authority. In the past, most of the studies conducted on the vibration analysis of gossamer structures used inaccurate or approximate theories in modeling the internal pressure. The toroidal shaped structure is one of the important key element in space application, helps to support the reflector in space application. This paper discusses the finite-element analysis of an inflated torus. The eigen-frequencies are obtained via three-dimensional small-strain elasticity theory, based on extremum energy principle. The two finite-element model (model-1 and model-2) have cases have been generated using a commercial finite-element package. The structure model-1 with shell element and model-2 with the combination of the mass of enclosed fluid (air) added to the shell elements have been taken for the study. The model-1 is computed with present analytical approach to understand the convergence rate and the accuracy. The convergence study is made available for the symmetric modes and anti-symmetric modes about the centroidal-axis plane, meeting the eigen-frequencies of an inflatable torus with the circular cross section. The structural model-2 is introduced with air mass element and analyzed its eigen-frequency with different aspect ratio and mode shape response using in-plane and out-plane loading condition are studied.

Error-controlled adaptive finite elements in solid mechanics

National Research Council Canada - National Science Library

Stein, Erwin; Ramm, E

2003-01-01

... . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Error-controlled Adaptive Finite-element-methods . . . . . . . . . . . . Missing Features and Properties of Today's General Purpose FE Programs for Structural...

Modelling of single walled carbon nanotube cylindrical structures with finite element method simulations

Science.gov (United States)

Günay, E.

2016-04-01

In this study, the modulus of elasticity and shear modulus values of single-walled carbon nanotubes SWCNTs were modelled by using both finite element method and the Matlab code. Initially, cylindrical armchair and zigzag single walled 3D space frames were demonstrated as carbon nanostructures. Thereafter, macro programs were written by the Matlab code producing the space truss for zigzag and armchair models. 3D space frames were introduced to the ANSYS software and then tension, compression and additionally torsion tests were performed on zigzag and armchair carbon nanotubes with BEAM4 element in obtaining the exact values of elastic and shear modulus values. In this study, two different boundary conditions were tested and especially used in torsion loading. The equivalent shear modulus data was found by averaging the corresponding values obtained from ten different nodal points on the nanotube path. Finally, in this study it was determined that the elastic constant values showed proportional changes by increasing the carbon nanotube diameters up to a certain level but beyond this level these values remained stable.

Modelling of single walled carbon nanotube cylindrical structures with finite element method simulations

International Nuclear Information System (INIS)

Günay, E.

2016-01-01

In this study, the modulus of elasticity and shear modulus values of single-walled carbon nanotubes SWCNTs were modelled by using both finite element method and the Matlab code. Initially, cylindrical armchair and zigzag single walled 3D space frames were demonstrated as carbon nanostructures. Thereafter, macro programs were written by the Matlab code producing the space truss for zigzag and armchair models. 3D space frames were introduced to the ANSYS software and then tension, compression and additionally torsion tests were performed on zigzag and armchair carbon nanotubes with BEAM4 element in obtaining the exact values of elastic and shear modulus values. In this study, two different boundary conditions were tested and especially used in torsion loading. The equivalent shear modulus data was found by averaging the corresponding values obtained from ten different nodal points on the nanotube path. Finally, in this study it was determined that the elastic constant values showed proportional changes by increasing the carbon nanotube diameters up to a certain level but beyond this level these values remained stable.

Modelling of single walled carbon nanotube cylindrical structures with finite element method simulations

Energy Technology Data Exchange (ETDEWEB)

Günay, E. [Gazi University, Mechanical Engineering Department, 06570, Ankara (Turkey)

2016-04-21

In this study, the modulus of elasticity and shear modulus values of single-walled carbon nanotubes SWCNTs were modelled by using both finite element method and the Matlab code. Initially, cylindrical armchair and zigzag single walled 3D space frames were demonstrated as carbon nanostructures. Thereafter, macro programs were written by the Matlab code producing the space truss for zigzag and armchair models. 3D space frames were introduced to the ANSYS software and then tension, compression and additionally torsion tests were performed on zigzag and armchair carbon nanotubes with BEAM4 element in obtaining the exact values of elastic and shear modulus values. In this study, two different boundary conditions were tested and especially used in torsion loading. The equivalent shear modulus data was found by averaging the corresponding values obtained from ten different nodal points on the nanotube path. Finally, in this study it was determined that the elastic constant values showed proportional changes by increasing the carbon nanotube diameters up to a certain level but beyond this level these values remained stable.

Stress analysis of heated concrete using finite elements

International Nuclear Information System (INIS)

Majumdar, P.; Gupta, A.; Marchertas, A.

1994-01-01

Described is a finite element analysis of concrete, which is subjected to rapid heating. Using thermal mass transport calculation, the moisture content, temperature and pore pressure distribution over space and time is obtained first. From these effects, stress at various points of the concrete are computed using the finite element method. Contribution to the stress formulation comes from three components, namely the thermal expansion, pore pressure, and the shrinkage of concrete due to moisture loss (from dehydration). The material properties of concrete are assumed to be homogeneous, elastic, and cracking is not taken into consideration. (orig.)

Finite element simulation of internal flows with heat transfer using a ...

Indian Academy of Sciences (India)

Unknown

Velocity correction method; finite element simulation; turbulent .... CFD, developments in turbulence modeling have been only evolutionary and ...... variables are made dimensionless using appropriate combinations of Uav, H, ...... Srinivas M 1994 Finite element analysis of internal flows with heat transfer Ph D thesis, Indian.

Integral finite element analysis of turntable bearing with flexible rings

Science.gov (United States)

Deng, Biao; Liu, Yunfei; Guo, Yuan; Tang, Shengjin; Su, Wenbin; Lei, Zhufeng; Wang, Pengcheng

2018-03-01

This paper suggests a method to calculate the internal load distribution and contact stress of the thrust angular contact ball turntable bearing by FEA. The influence of the stiffness of the bearing structure and the plastic deformation of contact area on the internal load distribution and contact stress of the bearing is considered. In this method, the load-deformation relationship of the rolling elements is determined by the finite element contact analysis of a single rolling element and the raceway. Based on this, the nonlinear contact between the rolling elements and the inner and outer ring raceways is same as a nonlinear compression spring and bearing integral finite element analysis model including support structure was established. The effects of structural deformation and plastic deformation on the built-in stress distribution of slewing bearing are investigated on basis of comparing the consequences of load distribution, inner and outer ring stress, contact stress and other finite element analysis results with the traditional bearing theory, which has guiding function for improving the design of slewing bearing.

Finite Element Simulation of Fracture Toughness Test

International Nuclear Information System (INIS)

Chu, Seok Jae; Liu, Cong Hao

2013-01-01

Finite element simulations of tensile tests were performed to determine the equivalent stress - equivalent plastic strain curves, critical equivalent stresses, and critical equivalent plastic strains. Then, the curves were used as inputs to finite element simulations of fracture toughness tests to determine the plane strain fracture toughness. The critical COD was taken as the COD when the equivalent plastic strain at the crack tip reached a critical value, and it was used as a crack growth criterion. The relationship between the critical COD and the critical equivalent plastic strain or the reduction of area was found. The relationship between the plane strain fracture toughness and the product of the critical equivalent stress and the critical equivalent plastic strain was also found

Semianalytic Design Sensitivity Analysis of Nonlinear Structures With a Commercial Finite Element Package

International Nuclear Information System (INIS)

Lee, Tae Hee; Yoo, Jung Hun; Choi, Hyeong Cheol

2002-01-01

A finite element package is often used as a daily design tool for engineering designers in order to analyze and improve the design. The finite element analysis can provide the responses of a system for given design variables. Although finite element analysis can quite well provide the structural behaviors for given design variables, it cannot provide enough information to improve the design such as design sensitivity coefficients. Design sensitivity analysis is an essential step to predict the change in responses due to a change in design variables and to optimize a system with the aid of the gradient-based optimization techniques. To develop a numerical method of design sensitivity analysis, analytical derivatives that are based on analytical differentiation of the continuous or discrete finite element equations are effective but analytical derivatives are difficult because of the lack of internal information of the commercial finite element package such as shape functions. Therefore, design sensitivity analysis outside of the finite element package is necessary for practical application in an industrial setting. In this paper, the semi-analytic method for design sensitivity analysis is used for the development of the design sensitivity module outside of a commercial finite element package of ANSYS. The direct differentiation method is employed to compute the design derivatives of the response and the pseudo-load for design sensitivity analysis is effectively evaluated by using the design variation of the related internal nodal forces. Especially, we suggest an effective method for stress and nonlinear design sensitivity analyses that is independent of the commercial finite element package is also discussed. Numerical examples are illustrated to show the accuracy and efficiency of the developed method and to provide insights for implementation of the suggested method into other commercial finite element packages