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.

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...

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 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

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),

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

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

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

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

Performance prediction of a permanent magnet servomotor using finite element analysis

International Nuclear Information System (INIS)

O'Dwyer, J.; Brady, G.; Byrne, J.V.

1998-01-01

Finite element analysis package ''MAGSOLVE'' is employed to predict the performance of a brushless permanent magnet servomotor. Equations are developed for a scaling method to predict motor rated stall torque. The torque v. current characteristic is determined using the Maxwell stress method and compared to experimental results. Less than 2% difference is observed over the current range. The solution is substantially independent of the location of the stress contour in the airgap. Valuable insight is also obtained regarding demagnetisation and torque ripple attributes of the motor. (orig.)

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.

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

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.

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....

Assessing performance and validating finite element simulations using probabilistic knowledge

Energy Technology Data Exchange (ETDEWEB)

Dolin, Ronald M.; Rodriguez, E. A. (Edward A.)

2002-01-01

Two probabilistic approaches for assessing performance are presented. The first approach assesses probability of failure by simultaneously modeling all likely events. The probability each event causes failure along with the event's likelihood of occurrence contribute to the overall probability of failure. The second assessment method is based on stochastic sampling using an influence diagram. Latin-hypercube sampling is used to stochastically assess events. The overall probability of failure is taken as the maximum probability of failure of all the events. The Likelihood of Occurrence simulation suggests failure does not occur while the Stochastic Sampling approach predicts failure. The Likelihood of Occurrence results are used to validate finite element predictions.

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...

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.

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.

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.

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

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

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)

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.

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.

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....

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 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)

Sensitivity analysis of the Galerkin finite element method neutron diffusion solver to the shape of the elements

Energy Technology Data Exchange (ETDEWEB)

Hosseini, Seyed Abolfaz [Dept. of Energy Engineering, Sharif University of Technology, Tehran (Iran, Islamic Republic of)

2017-02-15

The purpose of the present study is the presentation of the appropriate element and shape function in the solution of the neutron diffusion equation in two-dimensional (2D) geometries. To this end, the multigroup neutron diffusion equation is solved using the Galerkin finite element method in both rectangular and hexagonal reactor cores. The spatial discretization of the equation is performed using unstructured triangular and quadrilateral finite elements. Calculations are performed using both linear and quadratic approximations of shape function in the Galerkin finite element method, based on which results are compared. Using the power iteration method, the neutron flux distributions with the corresponding eigenvalue are obtained. The results are then validated against the valid results for IAEA-2D and BIBLIS-2D benchmark problems. To investigate the dependency of the results to the type and number of the elements, and shape function order, a sensitivity analysis of the calculations to the mentioned parameters is performed. It is shown that the triangular elements and second order of the shape function in each element give the best results in comparison to the other states.

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...

Node-based finite element method for large-scale adaptive fluid analysis in parallel environments

International Nuclear Information System (INIS)

Toshimitsu, Fujisawa; Genki, Yagawa

2003-01-01

In this paper, a FEM-based (finite element method) mesh free method with a probabilistic node generation technique is presented. In the proposed method, all computational procedures, from the mesh generation to the solution of a system of equations, can be performed fluently in parallel in terms of nodes. Local finite element mesh is generated robustly around each node, even for harsh boundary shapes such as cracks. The algorithm and the data structure of finite element calculation are based on nodes, and parallel computing is realized by dividing a system of equations by the row of the global coefficient matrix. In addition, the node-based finite element method is accompanied by a probabilistic node generation technique, which generates good-natured points for nodes of finite element mesh. Furthermore, the probabilistic node generation technique can be performed in parallel environments. As a numerical example of the proposed method, we perform a compressible flow simulation containing strong shocks. Numerical simulations with frequent mesh refinement, which are required for such kind of analysis, can effectively be performed on parallel processors by using the proposed method. (authors)

Node-based finite element method for large-scale adaptive fluid analysis in parallel environments

Energy Technology Data Exchange (ETDEWEB)

Toshimitsu, Fujisawa [Tokyo Univ., Collaborative Research Center of Frontier Simulation Software for Industrial Science, Institute of Industrial Science (Japan); Genki, Yagawa [Tokyo Univ., Department of Quantum Engineering and Systems Science (Japan)

2003-07-01

In this paper, a FEM-based (finite element method) mesh free method with a probabilistic node generation technique is presented. In the proposed method, all computational procedures, from the mesh generation to the solution of a system of equations, can be performed fluently in parallel in terms of nodes. Local finite element mesh is generated robustly around each node, even for harsh boundary shapes such as cracks. The algorithm and the data structure of finite element calculation are based on nodes, and parallel computing is realized by dividing a system of equations by the row of the global coefficient matrix. In addition, the node-based finite element method is accompanied by a probabilistic node generation technique, which generates good-natured points for nodes of finite element mesh. Furthermore, the probabilistic node generation technique can be performed in parallel environments. As a numerical example of the proposed method, we perform a compressible flow simulation containing strong shocks. Numerical simulations with frequent mesh refinement, which are required for such kind of analysis, can effectively be performed on parallel processors by using the proposed method. (authors)

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 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

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.

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.

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.

Stochastic Finite Element Analysis of Non-Linear Structures Modelled by Plasticity Theory

DEFF Research Database (Denmark)

Frier, Christian; Sørensen, John Dalsgaard

2003-01-01

A Finite Element Reliability Method (FERM) is introduced to perform reliability analyses on two-dimensional structures in plane stress, modeled by non-linear plasticity theory. FERM is a coupling between the First Order Reliability Method (FORM) and the Finite Element Method (FEM). FERM can be us...

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...

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.

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.

Optical Finite Element Processor

Science.gov (United States)

Casasent, David; Taylor, Bradley K.

1986-01-01

A new high-accuracy optical linear algebra processor (OLAP) with many advantageous features is described. It achieves floating point accuracy, handles bipolar data by sign-magnitude representation, performs LU decomposition using only one channel, easily partitions and considers data flow. A new application (finite element (FE) structural analysis) for OLAPs is introduced and the results of a case study presented. Error sources in encoded OLAPs are addressed for the first time. Their modeling and simulation are discussed and quantitative data are presented. Dominant error sources and the effects of composite error sources are analyzed.

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)

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

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.

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

Stochastic Finite Elements in Reliability-Based Structural Optimization

DEFF Research Database (Denmark)

Sørensen, John Dalsgaard; Engelund, S.

1995-01-01

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...... to optimization variables can be performed. A computer implementation is described and an illustrative example is given....

Finite element modelling

International Nuclear Information System (INIS)

Tonks, M.R.; Williamson, R.; Masson, R.

2015-01-01

The Finite Element Method (FEM) is a numerical technique for finding approximate solutions to boundary value problems. While FEM is commonly used to solve solid mechanics equations, it can be applied to a large range of BVPs from many different fields. FEM has been used for reactor fuels modelling for many years. It is most often used for fuel performance modelling at the pellet and pin scale, however, it has also been used to investigate properties of the fuel material, such as thermal conductivity and fission gas release. Recently, the United Stated Department Nuclear Energy Advanced Modelling and Simulation Program has begun using FEM as the basis of the MOOSE-BISON-MARMOT Project that is developing a multi-dimensional, multi-physics fuel performance capability that is massively parallel and will use multi-scale material models to provide a truly predictive modelling capability. (authors)

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.

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.

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

Evaluating the performance of the particle finite element method in parallel architectures

Science.gov (United States)

Gimenez, Juan M.; Nigro, Norberto M.; Idelsohn, Sergio R.

2014-05-01

This paper presents a high performance implementation for the particle-mesh based method called particle finite element method two (PFEM-2). It consists of a material derivative based formulation of the equations with a hybrid spatial discretization which uses an Eulerian mesh and Lagrangian particles. The main aim of PFEM-2 is to solve transport equations as fast as possible keeping some level of accuracy. The method was found to be competitive with classical Eulerian alternatives for these targets, even in their range of optimal application. To evaluate the goodness of the method with large simulations, it is imperative to use of parallel environments. Parallel strategies for Finite Element Method have been widely studied and many libraries can be used to solve Eulerian stages of PFEM-2. However, Lagrangian stages, such as streamline integration, must be developed considering the parallel strategy selected. The main drawback of PFEM-2 is the large amount of memory needed, which limits its application to large problems with only one computer. Therefore, a distributed-memory implementation is urgently needed. Unlike a shared-memory approach, using domain decomposition the memory is automatically isolated, thus avoiding race conditions; however new issues appear due to data distribution over the processes. Thus, a domain decomposition strategy for both particle and mesh is adopted, which minimizes the communication between processes. Finally, performance analysis running over multicore and multinode architectures are presented. The Courant-Friedrichs-Lewy number used influences the efficiency of the parallelization and, in some cases, a weighted partitioning can be used to improve the speed-up. However the total cputime for cases presented is lower than that obtained when using classical Eulerian strategies.

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

Finite element stress analysis of brick-mortar masonry under ...

African Journals Online (AJOL)

Stress analysis of a brick-mortar couplet as a substitute for brick wall structure has been performed by finite element method, and algorithm for determining the element stiffness matrix for a plane stress problem using the displacement approach was developed. The nodal displacements were derived for the stress in each ...

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...

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.

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.

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...

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.

Finite-element pre-analysis for pressurized thermoshock tests

International Nuclear Information System (INIS)

Keinaenen, H.; Talja, H.; Lehtonen, M.; Rintamaa, R.; Bljumin, A.; Timofeev, B.

1992-05-01

The behaviour of a model pressure vessel is studied in a pressurized thermal shock loading. The tests were performed at the Prometey Institute in St. Petersburg. The calculations were performed at the Technical Research Centre of Finland. The report describes the preliminary finite-element analyses for the fourth, fifth and sixth thermoshock tests with the first model pressure vessel. Seven pressurized thermoshock tests were made with the same model using five different flaw geometries. In the first three tests the flaw was actually a blunt notch. In the two following tests (tests 4 and 5) a sharp pre-crack was produced before the test. In the last two test (tests 6 and 7) the old crack was used. According to the measurements and post-test ultrasonic examination of the crack front, the sixth test led to significant crack extension. Both temperatures and stresses were calculated using the finite-element method. The calculations were made using the idealized initial flaw geometry and preliminary material data. Both two-and three-dimensional models were used in the calculations. J-integral values were calculated from the elastic-plastic finite-element results. The stress intensity factor values were evaluated on the basis of the calculated J-integrals and compared with the preliminary material fracture toughness data obtained from the Prometey Institute

Finite-element simulation of blood perfusion in muscle tissue during compression and sustained contraction.

Science.gov (United States)

Vankan, W J; Huyghe, J M; Slaaf, D W; van Donkelaar, C C; Drost, M R; Janssen, J D; Huson, A

1997-09-01

Mechanical interaction between tissue stress and blood perfusion in skeletal muscles plays an important role in blood flow impediment during sustained contraction. The exact mechanism of this interaction is not clear, and experimental investigation of this mechanism is difficult. We developed a finite-element model of the mechanical behavior of blood-perfused muscle tissue, which accounts for mechanical blood-tissue interaction in maximally vasodilated vasculature. Verification of the model was performed by comparing finite-element results of blood pressure and flow with experimental measurements in a muscle that is subject to well-controlled mechanical loading conditions. In addition, we performed simulations of blood perfusion during tetanic, isometric contraction and maximal vasodilation in a simplified, two-dimensional finite-element model of a rat calf muscle. A vascular waterfall in the venous compartment was identified as the main cause for blood flow impediment both in the experiment and in the finite-element simulations. The validated finite-element model offers possibilities for detailed analysis of blood perfusion in three-dimensional muscle models under complicated loading conditions.

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

Extended finite element method and its application in heterogeneous materials with inclusions

International Nuclear Information System (INIS)

Du Chengbin; Jiang Shouyan; Ying Zongquan

2010-01-01

To simplify the technology of finite element mesh generation for particle reinforced material, enrichment techniques is used to account for the material interfaces in the framework of extended finite element method (XFEM). The geometry of material distribution is described by level set function, which allows one to model the internal boundaries of the microstructure without the adaptation of the mesh. The enrichment function is used to improve the shape function of classical finite element method (FEM) for the nodes supporting the elements cut by the interface. The key issue of XFEM including constructing displacement pattern, establishment of the governing equation and scheme of numerical integration is also presented. It is not necessarily matching the internal features of the inclusions using XFEM, so the generation of finite element mesh can be performed easily. Finally, a plate with multi-circular inclusions under uniaxial tension is simulated by XFEM and FEM, respectively. The results show that XFEM is highly effective and efficient.

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

Finite element model updating of natural fibre reinforced composite structure in structural dynamics

Directory of Open Access Journals (Sweden)

Sani M.S.M.

2016-01-01

Full Text Available Model updating is a process of making adjustment of certain parameters of finite element model in order to reduce discrepancy between analytical predictions of finite element (FE and experimental results. Finite element model updating is considered as an important field of study as practical application of finite element method often shows discrepancy to the test result. The aim of this research is to perform model updating procedure on a composite structure as well as trying improving the presumed geometrical and material properties of tested composite structure in finite element prediction. The composite structure concerned in this study is a plate of reinforced kenaf fiber with epoxy. Modal properties (natural frequency, mode shapes, and damping ratio of the kenaf fiber structure will be determined using both experimental modal analysis (EMA and finite element analysis (FEA. In EMA, modal testing will be carried out using impact hammer test while normal mode analysis using FEA will be carried out using MSC. Nastran/Patran software. Correlation of the data will be carried out before optimizing the data from FEA. Several parameters will be considered and selected for the model updating procedure.

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.

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)

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

Improving the performance of finite element simulations on the wheel–rail interaction by using a coupling strategy

NARCIS (Netherlands)

Ma, Y.; Markine, V.L.; Ahad Mashal, Abdul; Ren, Mingfa

2018-01-01

Over the past few years, a number of implicit/explicit finite element models have been introduced for the purpose of tackling the problems of wheel–rail interaction. Yet, most of those finite element models encounter common numerical difficulties. For instance, initial gaps/penetrations between two

A holistic 3D finite element simulation model for thermoelectric power generator element

International Nuclear Information System (INIS)

Wu, Guangxi; Yu, Xiong

2014-01-01

Highlights: • Development of a holistic simulation model for the thermoelectric energy harvester. • Account for delta Seebeck coefficient and carrier charge densities variations. • Solution of thermo-electric coupling problem with finite element method. • Model capable of predicting phenomena not captured by traditional models. • A simulation tool for design of innovative TEM materials and structures. - Abstract: Harvesting the thermal energy stored in the ambient environment provides a potential sustainable energy source. Thermoelectric power generators have advantages of having no moving parts, being durable, and light-weighted. These unique features are advantageous for many applications (i.e., carry-on medical devices, embedded infrastructure sensors, aerospace, transportation, etc.). To ensure the efficient applications of thermoelectric energy harvesting system, the behaviors of such systems need to be fully understood. Finite element simulations provide important tools for such purpose. Although modeling the performance of thermoelectric modules has been conducted by many researchers, due to the complexity in solving the coupled problem, the influences of the effective Seebeck coefficient and carrier density variations on the performance of thermoelectric system are generally neglected. This results in an overestimation of the power generator performance under strong-ionization temperature region. This paper presents an advanced simulation model for thermoelectric elements that considers the effects of both factors. The mathematical basis of this model is firstly presented. Finite element simulations are then implemented on a thermoelectric power generator unit. The characteristics of the thermoelectric power generator and their relationship to its performance are discussed under different working temperature regions. The internal physics processes of the TEM harvester are analyzed from the results of computational simulations. The new model

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

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.

Mechanical Performance of Natural / Natural Fiber Reinforced Hybrid Composite Materials Using Finite Element Method Based Micromechanics and Experiments

OpenAIRE

Rahman, Muhammad Ziaur

2017-01-01

A micromechanical analysis of the representative volume element (RVE) of a unidirectional flax/jute fiber reinforced epoxy composite is performed using finite element analysis (FEA). To do so, first effective mechanical properties of flax fiber and jute fiber are evaluated numerically and then used in evaluating the effective properties of ax/jute/epoxy hybrid composite. Mechanics of Structure Genome (MSG), a new homogenization tool developed in Purdue University, is used to calculate the hom...

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 multilevel correction adaptive finite element method for Kohn-Sham equation

Science.gov (United States)

Hu, Guanghui; Xie, Hehu; Xu, Fei

2018-02-01

In this paper, an adaptive finite element method is proposed for solving Kohn-Sham equation with the multilevel correction technique. In the method, the Kohn-Sham equation is solved on a fixed and appropriately coarse mesh with the finite element method in which the finite element space is kept improving by solving the derived boundary value problems on a series of adaptively and successively refined meshes. A main feature of the method is that solving large scale Kohn-Sham system is avoided effectively, and solving the derived boundary value problems can be handled efficiently by classical methods such as the multigrid method. Hence, the significant acceleration can be obtained on solving Kohn-Sham equation with the proposed multilevel correction technique. The performance of the method is examined by a variety of numerical experiments.

Probabilistic finite elements for fracture and fatigue analysis

Science.gov (United States)

Liu, W. K.; Belytschko, T.; Lawrence, M.; Besterfield, G. H.

1989-01-01

The fusion of the probabilistic finite element method (PFEM) and reliability analysis for probabilistic fracture mechanics (PFM) is presented. A comprehensive method for determining the probability of fatigue failure for curved crack growth was developed. The criterion for failure or performance function is stated as: the fatigue life of a component must exceed the service life of the component; otherwise failure will occur. An enriched element that has the near-crack-tip singular strain field embedded in the element is used to formulate the equilibrium equation and solve for the stress intensity factors at the crack-tip. Performance and accuracy of the method is demonstrated on a classical mode 1 fatigue problem.

Finite element investigation of the prestressed jointed concrete ...

African Journals Online (AJOL)

Precast prestressed concrete pavement (PCP) technology is of recent origin, and the information on PCP performance is not available in literature. This research presents a finite-element analysis of the potential benefits of prestressing on the jointed concrete pavements (JCP). With using a 3-dimensional (3D) ...

OpenSeesPy: Python library for the OpenSees finite element framework

Science.gov (United States)

Zhu, Minjie; McKenna, Frank; Scott, Michael H.

2018-01-01

OpenSees, an open source finite element software framework, has been used broadly in the earthquake engineering community for simulating the seismic response of structural and geotechnical systems. The framework allows users to perform finite element analysis with a scripting language and for developers to create both serial and parallel finite element computer applications as interpreters. For the last 15 years, Tcl has been the primary scripting language to which the model building and analysis modules of OpenSees are linked. To provide users with different scripting language options, particularly Python, the OpenSees interpreter interface was refactored to provide multi-interpreter capabilities. This refactoring, resulting in the creation of OpenSeesPy as a Python module, is accomplished through an abstract interface for interpreter calls with concrete implementations for different scripting languages. Through this approach, users are able to develop applications that utilize the unique features of several scripting languages while taking advantage of advanced finite element analysis models and algorithms.

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....

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

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

A Timoshenko Piezoelectric Beam Finite Element with Consistent Performance Irrespective of Geometric and Material Configurations

Directory of Open Access Journals (Sweden)

Litesh N. Sulbhewar

Full Text Available Abstract The conventional Timoshenko piezoelectric beam finite elements based on First-order Shear Deformation Theory (FSDT do not maintain the accuracy and convergence consistently over the applicable range of material and geometric properties. In these elements, the inaccuracy arises due to the induced potential effects in the transverse direction and inefficiency arises due to the use of independently assumed linear polynomial interpolation of the field variables in the longitudinal direction. In this work, a novel FSDT-based piezoelectric beam finite element is proposed which is devoid of these deficiencies. A variational formulation with consistent through-thickness potential is developed. The governing equilibrium equations are used to derive the coupled field relations. These relations are used to develop a polynomial interpolation scheme which properly accommodates the bending-extension, bending-shear and induced potential couplings to produce accurate results in an efficient manner. It is noteworthy that this consistently accurate and efficient beam finite element uses the same nodal variables as of conventional FSDT formulations available in the literature. Comparison of numerical results proves the consistent accuracy and efficiency of the proposed formulation irrespective of geometric and material configurations, unlike the conventional formulations.

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

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.

Solid Modeling and Finite Element Analysis of an Overhead Crane Bridge

Directory of Open Access Journals (Sweden)

C. Alkin

2005-01-01

Full Text Available The design of an overhead crane bridge with a double box girder has been investigated and a case study of a crane with 35 ton capacity and 13 m span length has been conducted. In the initial phase of the case study, conventional design calculations proposed by F. E. M. Rules and DIN standards were performed to verify the stress and deflection levels. The crane design was modeled using both solids and surfaces. Finite element meshes with 4-node tetrahedral and 4-node quadrilateral shell elements were generated from the solid and shell models, respectively. After a comparison of the finite element analyses, the conventional calculations and performance of the existing crane, the analysis with quadratic shell elements was found to give the most realistic results. As a result of this study, a design optimization method for an overhead crane is proposed. 

A preliminary investigation of finite-element modeling for composite rotor blades

Science.gov (United States)

Lake, Renee C.; Nixon, Mark W.

1988-01-01

The results from an initial phase of an in-house study aimed at improving the dynamic and aerodynamic characteristics of composite rotor blades through the use of elastic couplings are presented. Large degree of freedom shell finite element models of an extension twist coupled composite tube were developed and analyzed using MSC/NASTRAN. An analysis employing a simplified beam finite element representation of the specimen with the equivalent engineering stiffness was additionally performed. Results from the shell finite element normal modes and frequency analysis were compared to those obtained experimentally, showing an agreement within 13 percent. There was appreciable degradation in the frequency prediction for the torsional mode, which is elastically coupled. This was due to the absence of off-diagonal coupling terms in the formulation of the equivalent engineering stiffness. Parametric studies of frequency variation due to small changes in ply orientation angle and ply thickness were also performed. Results showed linear frequency variations less than 2 percent per 1 degree variation in the ply orientation angle, and 1 percent per 0.0001 inch variation in the ply thickness.

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

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...

Summary compilation of shell element performance versus formulation.

Energy Technology Data Exchange (ETDEWEB)

Heinstein, Martin Wilhelm; Hales, Jason Dean (Idaho National Laboratory, Idaho Falls, ID); Breivik, Nicole L.; Key, Samuel W. (FMA Development, LLC, Great Falls, MT)

2011-07-01

This document compares the finite element shell formulations in the Sierra Solid Mechanics code. These are finite elements either currently in the Sierra simulation codes Presto and Adagio, or expected to be added to them in time. The list of elements are divided into traditional two-dimensional, plane stress shell finite elements, and three-dimensional solid finite elements that contain either modifications or additional terms designed to represent the bending stiffness expected to be found in shell formulations. These particular finite elements are formulated for finite deformation and inelastic material response, and, as such, are not based on some of the elegant formulations that can be found in an elastic, infinitesimal finite element setting. Each shell element is subjected to a series of 12 verification and validation test problems. The underlying purpose of the tests here is to identify the quality of both the spatially discrete finite element gradient operator and the spatially discrete finite element divergence operator. If the derivation of the finite element is proper, the discrete divergence operator is the transpose of the discrete gradient operator. An overall summary is provided from which one can rank, at least in an average sense, how well the individual formulations can be expected to perform in applications encountered year in and year out. A letter grade has been assigned albeit sometimes subjectively for each shell element and each test problem result. The number of A's, B's, C's, et cetera assigned have been totaled, and a grade point average (GPA) has been computed, based on a 4.0-system. These grades, combined with a comparison between the test problems and the application problem, can be used to guide an analyst to select the element with the best shell formulation.

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.

Finite element analysis of nonlinear creeping flows

International Nuclear Information System (INIS)

Loula, A.F.D.; Guerreiro, J.N.C.

1988-12-01

Steady-state creep problems with monotone constitutive laws are studied. Finite element approximations are constructed based on mixed Petrov-Galerkin formulations for constrained problems. Stability, convergence and a priori error estimates are proved for equal-order discontinuous stress and continuous velocity interpolations. Numerical results are presented confirming the rates of convergence predicted in the analysis and the good performance of this formulation. (author) [pt

Finite-element analysis of elastic sound-proof coupling thermal state

Science.gov (United States)

Tsyss, V. G.; Strokov, I. M.; Sergaeva, M. Yu

2018-01-01

The aim is in calculated determining of the elastic rubber-metal element thermal state of soundproof coupling ship shafting under variable influence during loads in time. Thermal coupling calculation is performed with finite element method using NX Simens software with Nastran solver. As a result of studies, the following results were obtained: - a volumetric picture of the temperature distribution over the array of the deformed coupling body is obtained; - time to reach steady-state thermal coupling mode has been determined; - dependences of maximum temperature and time to reach state on the established operation mode on rotation frequency and ambient temperature are determined. The findings prove the conclusion that usage of finite element analysis modern software can significantly speed up problem solving.

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

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.

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)

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.

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

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

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.

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.

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.

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

SPLAI: Computational Finite Element Model for Sensor Networks

Directory of Open Access Journals (Sweden)

Ruzana Ishak

2006-01-01

Full Text Available Wireless sensor network refers to a group of sensors, linked by a wireless medium to perform distributed sensing task. The primary interest is their capability in monitoring the physical environment through the deployment of numerous tiny, intelligent, wireless networked sensor nodes. Our interest consists of a sensor network, which includes a few specialized nodes called processing elements that can perform some limited computational capabilities. In this paper, we propose a model called SPLAI that allows the network to compute a finite element problem where the processing elements are modeled as the nodes in the linear triangular approximation problem. Our model also considers the case of some failures of the sensors. A simulation model to visualize this network has been developed using C++ on the Windows environment.

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

The reliability of finite element analysis results of the low impact test in predicting the energy absorption performance of thin-walled structures

Energy Technology Data Exchange (ETDEWEB)

Alipour, R.; Nejadx, Farokhi A.; Izman, S. [Universiti Teknologi Malaysia, Johor Bahru (Malaysia)

2015-05-15

The application of dual phase steels (DPS) such as DP600 in the form of thin-walled structure in automotive components is being continuously increased as vehicle designers utilize modern steel grades and low weight structures to improve structural performance, make automotive light and reinforce crash performance. Preventing cost enhancement of broad investigations in this area can be gained by using computers in structural analysis in order to substitute lots of experiments with finite element analysis (FEA). Nevertheless, it necessitates to be certified that selected method including element type and solution methodology is capable of predicting real condition. In this paper, numerical and experimental studies are done to specify the effect of element type selection and solution methodology on the results of finite element analysis in order to investigate the energy absorption behavior of a DP600 thin-walled structure with three different geometries under a low impact loading. The outcomes indicated the combination of implicit method and solid elements is in better agreement with the experiments. In addition, using a combination of shell element types with implicit method reduces the time of simulation remarkably, although the error of results compared to the experiments increased to some extent.

Multilayer Finite-Element Model Application to Define the Bearing Structure Element Stress State of Launch Complexes

Directory of Open Access Journals (Sweden)

V. A. Zverev

2016-01-01

Full Text Available The article objective is to justify the rationale for selecting the multilayer finite element model parameters of the bearing structure of a general-purpose launch complex unit.A typical design element of the launch complex unit, i.e. a mount of the hydraulic or pneumatic cylinder, block, etc. is under consideration. The mount represents a set of the cantilevered axis and external structural cage. The most loaded element of the cage is disk to which a moment is transferred from the cantilevered axis due to actuator effort acting on it.To calculate the stress-strain state of disk was used a finite element method. Five models of disk mount were created. The only difference in models was the number of layers of the finite elements through the thickness of disk. There were models, which had one, three, five, eight, and fourteen layers of finite elements through the thickness of disk. For each model, we calculated the equivalent stresses arising from the action of the test load. Disk models were formed and calculated using the MSC Nastran complex software.The article presents results in the table to show data of equivalent stresses in each of the multi-layered models and graphically to illustrate the changing equivalent stresses through the thickness of disk.Based on these results we have given advice on selecting the proper number of layers in the model allowing a desirable accuracy of results with the lowest run time. In addition, it is concluded that there is a need to use the multi-layer models in assessing the performance of structural elements in case the stress exceeds the allowable one in their surface layers.

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.

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)

Finite element predictions of active buckling control of stiffened panels

Science.gov (United States)

Thompson, Danniella M.; Griffin, O. H., Jr.

1993-04-01

Materials systems and structures that can respond 'intelligently' to their environment are currently being proposed and investigated. A series of finite element analyses was performed to investigate the potential for active buckling control of two different stiffened panels by embedded shape memory alloy (SMA) rods. Changes in the predicted buckling load increased with the magnitude of the actuation level for a given structural concept. Increasing the number of actuators for a given concept yielded greater predicted increases in buckling load. Considerable control authority was generated with a small number of actuators, with greater authority demonstrated for those structural concepts where the activated SMA rods could develop greater forces and moments on the structure. Relatively simple and inexpensive analyses were performed with standard finite elements to determine such information, indicating the viability of these types of models for design purposes.

Stress analysis in pressure vessels by mixed finite element methods taking into account shear deformation

International Nuclear Information System (INIS)

Franca, L.P.; Toledo, E.M.; Loula, A.F.D.; Garcia, E.L.M.

1988-12-01

A new finite element method is employed to approximate axisymmetric shell problems. This formulation enhances stability and accuracy, from thin to moderately thick shells, compared to the correspondent Galerkin finite element approximations. Numerical results illustrate the good performance of the present method on some typical pressure vessels aplications. (author) [pt

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.

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 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 ...

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)

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)

Dual and primal mixed Petrov-Galerkin finite element methods in heat transfer problems

International Nuclear Information System (INIS)

Loula, A.F.D.; Toledo, E.M.

1988-12-01

New mixed finite element formulations for the steady state heat transfer problem are presented with no limitation in the choice of conforming finite element spaces. Adding least square residual forms of the governing equations of the classical Galerkin formulation the original saddle point problem is transformed into a minimization problem. Stability analysis, error estimates and numerical results are presented, confirming the error estimates and the good performance of this new formulation. (author) [pt

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.

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.

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...

Parallel eigenanalysis of finite element models in a completely connected architecture

Science.gov (United States)

Akl, F. A.; Morel, M. R.

1989-01-01

A parallel algorithm is presented for the solution of the generalized eigenproblem in linear elastic finite element analysis, (K)(phi) = (M)(phi)(omega), where (K) and (M) are of order N, and (omega) is order of q. The concurrent solution of the eigenproblem is based on the multifrontal/modified subspace method and is achieved in a completely connected parallel architecture in which each processor is allowed to communicate with all other processors. The algorithm was successfully implemented on a tightly coupled multiple-instruction multiple-data parallel processing machine, Cray X-MP. A finite element model is divided into m domains each of which is assumed to process n elements. Each domain is then assigned to a processor or to a logical processor (task) if the number of domains exceeds the number of physical processors. The macrotasking library routines are used in mapping each domain to a user task. Computational speed-up and efficiency are used to determine the effectiveness of the algorithm. The effect of the number of domains, the number of degrees-of-freedom located along the global fronts and the dimension of the subspace on the performance of the algorithm are investigated. A parallel finite element dynamic analysis program, p-feda, is documented and the performance of its subroutines in parallel environment is analyzed.

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 stress and deflection analysis of CDF yike and end plug

International Nuclear Information System (INIS)

Wands, R.; Grimson, J.; Kephart, R.; Theriot, D.

1982-01-01

A large detector is being designed to study anti pp collisions at center-of-mass energies of up to 2000 GeV as part of the Fermilab Collider Detector Facility (CDF). The central detector of this facility consists of a solenoid, calorimeter yoke, and a variety of particle measurement devices. The yoke will be a large steel structure that will provide the magnetic flux return path as well as support structure for calorimetry and other instrumentation. It must resist both electromagnetic and gravitational loads while exhibiting only small elastic deformations. The instrumented endplugs of the yoke are subjected to large electromagnetic loads. Moreover, due to the presence of wire chambers within these plugs, they must also be particularly stiff. The purpose of this paper is to present the results of a finite element stress and deflection analysis of these structures under various anticipated load conditions. The PATRAN-G finite element modeling program, installed on a CDF-VAX 11/780 and operating from a Ramtek 6212 colorgraphics terminal, was used to generate the analysis models. The actual finite element analysis was performed by the ANSYS general purpose finite element program, installed on the Fermilab Cyber 175's

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.

Finite element modeling of micromachined MEMS photon devices

Science.gov (United States)

Evans, Boyd M., III; Schonberger, D. W.; Datskos, Panos G.

1999-09-01

The technology of microelectronics that has evolved over the past half century is one of great power and sophistication and can now be extended to many applications (MEMS and MOEMS) other than electronics. An interesting application of MEMS quantum devices is the detection of electromagnetic radiation. The operation principle of MEMS quantum devices is based on the photoinduced stress in semiconductors, and the photon detection results from the measurement of the photoinduced bending. These devices can be described as micromechanical photon detectors. In this work, we have developed a technique for simulating electronic stresses using finite element analysis. We have used our technique to model the response of micromechanical photon devices to external stimuli and compared these results with experimental data. Material properties, geometry, and bimaterial design play an important role in the performance of micromechanical photon detectors. We have modeled these effects using finite element analysis and included the effects of bimaterial thickness coating, effective length of the device, width, and thickness.

Finite Element Modeling of Micromachined MEMS Photon Devices

International Nuclear Information System (INIS)

Datskos, P.G.; Evans, B.M.; Schonberger, D.

1999-01-01

The technology of microelectronics that has evolved over the past half century is one of great power and sophistication and can now be extended to many applications (MEMS and MOEMS) other than electronics. An interesting application of MEMS quantum devices is the detection of electromagnetic radiation. The operation principle of MEMS quantum devices is based on the photoinduced stress in semiconductors, and the photon detection results from the measurement of the photoinduced bending. These devices can be described as micromechanical photon detectors. In this work, we have developed a technique for simulating electronic stresses using finite element analysis. We have used our technique to model the response of micromechanical photon devices to external stimuli and compared these results with experimental data. Material properties, geometry, and bimaterial design play an important role in the performance of micromechanical photon detectors. We have modeled these effects using finite element analysis and included the effects of bimaterial thickness coating, effective length of the device, width, and thickness

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.

Finite element analysis of cutting tools prior to fracture in hard turning operations

International Nuclear Information System (INIS)

Cakir, M. Cemal; I Sik, Yahya

2005-01-01

In this work cutting FEA of cutting tools prior to fracture is investigated. Fracture is the catastrophic end of the cutting edge that should be avoided for the cutting tool in order to have a longer tool life. This paper presents finite element modelling of a cutting tool just before its fracture. The data used in FEA are gathered from a tool breakage system that detects the fracture according to the variations of the cutting forces measured by a three-dimensional force dynamometer. The workpiece material used in the experiments is cold work tool steel, AISI O1 (60 HRC) and the cutting tool material is uncoated tungsten carbide (DNMG 150608). In order to investigate the cutting tool conditions in longitudinal external turning operations prior to fracture, static and dynamic finite element analyses are conducted. After the static finite element analysis, the modal and harmonic response analyses are carried on and the dynamic behaviours of the cutting tool structure are investigated. All FE analyses were performed using a commercial finite element package ANSYS

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...

Understanding the forming and performance of aluminum beverage cans through finite element modelling

International Nuclear Information System (INIS)

MacEwen, S.R.; Langille, A.; Hamstra, P.; Wu, P-D.; Savoie, J.; Nardini, D.

2000-01-01

The past decade has witnessed major changes in the aluminum beverage can business, for both the canmakers and the suppliers of aluminum sheet. The major driving force for change was 'light-weighting', with the gauge of can body stock decreasing from 0.0124inches to 0.0102inches or less at present. The resulting reduction in weight of a can had to be achieved without compromise to the performance of the can, as defined by dome reversal pressure, axial buckle load and drop resistance. Furthermore, the fundamentals of the draw-and- iron process used to manufacture cans had to remain essentially unchanged, despite the new problems arising from downgauging. This presentation will demonstrate how finite element modeling of can forming processes and can performance tests has assisted in understanding the interaction of the tooling with the sheet and how the workhardening and residual stresses produced by the forming operations affect subsequent performance. The roles of sheet characterization and model validation will be discussed for several of the major canmaking steps, including cup draw and redraw and the forming of the bottom profile. (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.

Finite cover method with mortar elements for elastoplasticity problems

Science.gov (United States)

Kurumatani, M.; Terada, K.

2005-06-01

Finite cover method (FCM) is extended to elastoplasticity problems. The FCM, which was originally developed under the name of manifold method, has recently been recognized as one of the generalized versions of finite element methods (FEM). Since the mesh for the FCM can be regular and squared regardless of the geometry of structures to be analyzed, structural analysts are released from a burdensome task of generating meshes conforming to physical boundaries. Numerical experiments are carried out to assess the performance of the FCM with such discretization in elastoplasticity problems. Particularly to achieve this accurately, the so-called mortar elements are introduced to impose displacement boundary conditions on the essential boundaries, and displacement compatibility conditions on material interfaces of two-phase materials or on joint surfaces between mutually incompatible meshes. The validity of the mortar approximation is also demonstrated in the elastic-plastic FCM.

Finite element modelling and updating of friction stir welding (FSW joint for vibration analysis

Directory of Open Access Journals (Sweden)

Zahari Siti Norazila

2017-01-01

Full Text Available Friction stir welding of aluminium alloys widely used in automotive and aerospace application due to its advanced and lightweight properties. The behaviour of FSW joints plays a significant role in the dynamic characteristic of the structure due to its complexities and uncertainties therefore the representation of an accurate finite element model of these joints become a research issue. In this paper, various finite elements (FE modelling technique for prediction of dynamic properties of sheet metal jointed by friction stir welding will be presented. Firstly, nine set of flat plate with different series of aluminium alloy; AA7075 and AA6061 joined by FSW are used. Nine set of specimen was fabricated using various types of welding parameters. In order to find the most optimum set of FSW plate, the finite element model using equivalence technique was developed and the model validated using experimental modal analysis (EMA on nine set of specimen and finite element analysis (FEA. Three types of modelling were engaged in this study; rigid body element Type 2 (RBE2, bar element (CBAR and spot weld element connector (CWELD. CBAR element was chosen to represent weld model for FSW joints due to its accurate prediction of mode shapes and contains an updating parameter for weld modelling compare to other weld modelling. Model updating was performed to improve correlation between EMA and FEA and before proceeds to updating, sensitivity analysis was done to select the most sensitive updating parameter. After perform model updating, total error of the natural frequencies for CBAR model is improved significantly. Therefore, CBAR element was selected as the most reliable element in FE to represent FSW weld joint.

Free material stiffness design of laminated composite structures using commercial finite element analysis codes

DEFF Research Database (Denmark)

Henrichsen, Søren Randrup; Lindgaard, Esben; Lund, Erik

2015-01-01

In this work optimum stiffness design of laminated composite structures is performed using the commercially available programs ANSYS and MATLAB. Within these programs a Free Material Optimization algorithm is implemented based on an optimality condition and a heuristic update scheme. The heuristic...... update scheme is needed because commercially available finite element analysis software is used. When using a commercial finite element analysis code it is not straight forward to implement a computationally efficient gradient based optimization algorithm. Examples considered in this work are a clamped......, where full access to the finite element analysis core is granted. This comparison displays the possibility of using commercially available programs for stiffness design of laminated composite structures....

Analysis of Tire Tractive Performance on Deformable Terrain by Finite Element-Discrete Element Method

Science.gov (United States)

Nakashima, Hiroshi; Takatsu, Yuzuru

The goal of this study is to develop a practical and fast simulation tool for soil-tire interaction analysis, where finite element method (FEM) and discrete element method (DEM) are coupled together, and which can be realized on a desktop PC. We have extended our formerly proposed dynamic FE-DE method (FE-DEM) to include practical soil-tire system interaction, where not only the vertical sinkage of a tire, but also the travel of a driven tire was considered. Numerical simulation by FE-DEM is stable, and the relationships between variables, such as load-sinkage and sinkage-travel distance, and the gross tractive effort and running resistance characteristics, are obtained. Moreover, the simulation result is accurate enough to predict the maximum drawbar pull for a given tire, once the appropriate parameter values are provided. Therefore, the developed FE-DEM program can be applied with sufficient accuracy to interaction problems in soil-tire systems.

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

An embedded crack in a constant strain triangle utilizing extended finite element concepts

DEFF Research Database (Denmark)

Olesen, J.F.; Poulsen, P.N.

2013-01-01

This paper revisits the formulation of the CST element with an embedded discrete crack taking advantage of the direct formulations developed within the framework of the extended finite element method, XFEM. The result is a simple element for modeling cohesive fracture processes in quasi-brittle m......This paper revisits the formulation of the CST element with an embedded discrete crack taking advantage of the direct formulations developed within the framework of the extended finite element method, XFEM. The result is a simple element for modeling cohesive fracture processes in quasi......-element discontinuity of displacements. The formulation is based on a variational principle of virtual work involving only the interpolation of displacements. The good performance of the element is demonstrated through the comparison with three benchmark tests in which a single crack is propagated: The center cracked...

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.

Assessment of finite element and smoothed particles hydrodynamics methods for modeling serrated chip formation in hardened steel

Directory of Open Access Journals (Sweden)

Usama Umer

2016-05-01

Full Text Available This study aims to perform comparative analyses in modeling serrated chip morphologies using traditional finite element and smoothed particles hydrodynamics methods. Although finite element models are being employed in predicting machining performance variables for the last two decades, many drawbacks and limitations exist with the current finite element models. The problems like excessive mesh distortions, high numerical cost of adaptive meshing techniques, and need of geometric chip separation criteria hinder its practical implementation in metal cutting industries. In this study, a mesh free method, namely, smoothed particles hydrodynamics, is implemented for modeling serrated chip morphology while machining AISI H13 hardened tool steel. The smoothed particles hydrodynamics models are compared with the traditional finite element models, and it has been found that the smoothed particles hydrodynamics models have good capabilities in handling large distortions and do not need any geometric or mesh-based chip separation criterion.

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)

TITUS: a general finite element system

International Nuclear Information System (INIS)

Bougrelle, P.

1983-01-01

TITUS is a general finite element structural analysis system which performs linear/non-linear, static/dynamic analyses of heat-transfer/thermo-mechanical problems. One of the major features of TITUS is that it was designed by engineers, to address engineers in an industrial environment. This has resulted in an easy to use system, with a high-level free-formatted problem oriented language, a large selection of pre- and post processors and sophisticated graphic capabilities. TITUS has many references in civil, mechanical and nuclear engineering applications. The TITUS system is available on various types of machines, from large mainframes to mini computers

3D finite element simulation of optical modes in VCSELs

OpenAIRE

Rozova, M.; Pomplun, J.; Zschiedrich, L.; Schmidt, F.; Burger, S.

2011-01-01

We present a finite element method (FEM) solver for computation of optical resonance modes in VCSELs. We perform a convergence study and demonstrate that high accuracies for 3D setups can be attained on standard computers. We also demonstrate simulations of thermo-optical effects in VCSELs.

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)

contact stress analysis of involute spur gear by finite element method

African Journals Online (AJOL)

USER

section using C- programming. Bending stress analysis has been performed using finite element analysis with ANSYS software. Comparison of bending stress analysis has been performed for symmetric and asymmetric spur gear tooth at critical section. Mushin J. Jweeg, et.al. [7] used 2D contact stress FEA model to ...

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.

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

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...

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 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.

Coupled thermomechanical behavior of graphene using the spring-based finite element approach

Energy Technology Data Exchange (ETDEWEB)

Georgantzinos, S. K., E-mail: sgeor@mech.upatras.gr; Anifantis, N. K., E-mail: nanif@mech.upatras.gr [Machine Design Laboratory, Department of Mechanical Engineering and Aeronautics, University of Patras, Rio, 26500 Patras (Greece); Giannopoulos, G. I., E-mail: ggiannopoulos@teiwest.gr [Materials Science Laboratory, Department of Mechanical Engineering, Technological Educational Institute of Western Greece, 1 Megalou Alexandrou Street, 26334 Patras (Greece)

2016-07-07

The prediction of the thermomechanical behavior of graphene using a new coupled thermomechanical spring-based finite element approach is the aim of this work. Graphene sheets are modeled in nanoscale according to their atomistic structure. Based on molecular theory, the potential energy is defined as a function of temperature, describing the interatomic interactions in different temperature environments. The force field is approached by suitable straight spring finite elements. Springs simulate the interatomic interactions and interconnect nodes located at the atomic positions. Their stiffness matrix is expressed as a function of temperature. By using appropriate boundary conditions, various different graphene configurations are analyzed and their thermo-mechanical response is approached using conventional finite element procedures. A complete parametric study with respect to the geometric characteristics of graphene is performed, and the temperature dependency of the elastic material properties is finally predicted. Comparisons with available published works found in the literature demonstrate the accuracy of the proposed method.

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.

Transient thermal finite element analysis of CFC–Cu ITER monoblock using X-ray tomography data

Energy Technology Data Exchange (ETDEWEB)

Evans, Ll.M., E-mail: llion.evans@ccfe.ac.uk [CCFE, Culham Science Centre, Abingdon, Oxon OX14 3DB (United Kingdom); School of Materials, University of Manchester, Grosvenor Street, Manchester M1 7HS (United Kingdom); Margetts, L. [School of Earth, Atmospheric and Environmental Sciences, University of Manchester, Williamson Building, Manchester M13 9PL (United Kingdom); Casalegno, V. [Department of Applied Science and Technology, Politecnico di Torino, Corso Duca degli Abruzzi 24, I-10129 Torino (Italy); Lever, L.M. [IT Services for Research, University of Manchester, Devonshire House, Oxford Road, Manchester M13 9PL (United Kingdom); Bushell, J.; Lowe, T.; Wallwork, A. [School of Materials, University of Manchester, Grosvenor Street, Manchester M1 7HS (United Kingdom); Young, P. [Simpleware Ltd., Bradninch Hall, Castle Street, Exeter EX4 3PL (United Kingdom); Lindemann, A. [NETZSCH-Gerätebau GmbH, Wittelsbacherstraße 42, D-95100 Selb, Bayern (Germany); Schmidt, M.; Mummery, P.M. [School of Mechanical, Aerospace and Civil Engineering (MACE), University of Manchester, Manchester M13 9PL (United Kingdom)

2015-11-15

Highlights: • Thermal performance of a fusion power heat exchange component was investigated. • Microstructures effecting performance were determined using X-ray tomography. • This data was used to perform a microstructurally faithful finite element analysis. • FEA demonstrated that manufacturing defects had an appreciable effect on performance. • This image-based modelling showed which regions could be targeted for improvements. - Abstract: The thermal performance of a carbon fibre composite-copper monoblock, a sub-component of a fusion reactor divertor, was investigated by finite element analysis. High-accuracy simulations were created using an emerging technique, image-based finite element modelling, which converts X-ray tomography data into micro-structurally faithful models, capturing details such as manufacturing defects. For validation, a case study was performed where the thermal analysis by laser flash of a carbon fibre composite-copper disc was simulated such that computational and experimental results could be compared directly. Results showed that a high resolution image-based simulation (102 million elements of 32 μm width) provided increased accuracy over a low resolution image-based simulation (0.6 million elements of 194 μm width) and idealised computer aided design simulations. Using this technique to analyse a monoblock mock-up, it was possible to detect and quantify the effects of debonding regions at the carbon fibre composite-copper interface likely to impact both component performance and expected lifetime. These features would not have been accounted for in idealised computer aided design simulations.

Transient thermal finite element analysis of CFC–Cu ITER monoblock using X-ray tomography data

International Nuclear Information System (INIS)

Evans, Ll.M.; Margetts, L.; Casalegno, V.; Lever, L.M.; Bushell, J.; Lowe, T.; Wallwork, A.; Young, P.; Lindemann, A.; Schmidt, M.; Mummery, P.M.

2015-01-01

Highlights: • Thermal performance of a fusion power heat exchange component was investigated. • Microstructures effecting performance were determined using X-ray tomography. • This data was used to perform a microstructurally faithful finite element analysis. • FEA demonstrated that manufacturing defects had an appreciable effect on performance. • This image-based modelling showed which regions could be targeted for improvements. - Abstract: The thermal performance of a carbon fibre composite-copper monoblock, a sub-component of a fusion reactor divertor, was investigated by finite element analysis. High-accuracy simulations were created using an emerging technique, image-based finite element modelling, which converts X-ray tomography data into micro-structurally faithful models, capturing details such as manufacturing defects. For validation, a case study was performed where the thermal analysis by laser flash of a carbon fibre composite-copper disc was simulated such that computational and experimental results could be compared directly. Results showed that a high resolution image-based simulation (102 million elements of 32 μm width) provided increased accuracy over a low resolution image-based simulation (0.6 million elements of 194 μm width) and idealised computer aided design simulations. Using this technique to analyse a monoblock mock-up, it was possible to detect and quantify the effects of debonding regions at the carbon fibre composite-copper interface likely to impact both component performance and expected lifetime. These features would not have been accounted for in idealised computer aided design simulations.

Finite element study of a lumbar intervertebral disc nucleus replacement device

Directory of Open Access Journals (Sweden)

Jessica S Coogan

2016-12-01

Full Text Available Nucleus replacement technologies are a minimally invasive alternative to spinal fusion and total disc replacement that have the potential to reduce pain and restore motion for patients with degenerative disc disease. Finite element modeling can be used to determine the biomechanics associated with nucleus replacement technologies. The current study focuses on a new nucleus replacement device designed as a conforming silicone implant with an internal void. A validated finite element model of the human lumbar L3-L4 motion segment was developed and used to investigate the influence of the nucleus replacement device on spine biomechanics. In addition, the effect of device design changes on biomechanics was determined. A 3D, L3-L4 finite element model was constructed from medical imaging data. Models were created with the normal intact nucleus, the nucleus replacement device, and a solid silicone implant. Probabilistic analysis was performed on the normal model to provide quantitative validation metrics. Sensitivity analysis was performed on the silicone Shore A durometer of the device. Models were loaded under axial compression followed by flexion/extension, lateral bending, or axial rotation. Compressive displacement, endplate stresses, reaction moment, and annulus stresses were determined and compared between the different models. The novel nucleus replacement device resulted in similar compressive displacement, endplate stress, and annulus stress and slightly higher reaction moment compared with the normal nucleus. The solid implant resulted in decreased displacement, increased endplate stress, decreased annulus stress, and decreased reaction moment compared with the novel device. With increasing silicone durometer, compressive displacement decreased, endplate stress increased, reaction moment increased, and annulus stress decreased. Finite element analysis was used to show that the novel nucleus replacement device results in similar biomechanics

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)

The impact of process variables and wear characteristics on the cutting tool performance using Finite Element Analysis

OpenAIRE

Jiang, Xiaoheng

2016-01-01

The frequent failure of cutting tool in the cutting process may cause a huge loss of money and time especially for hard to machine materials such as titanium alloys. Thus this study is mainly focused on the impact of wear characteristics and process variables on the cutting tool which is ignored by most of researchers. A thermo-mechanical finite element model of orthogonal metal cutting with segment chip formation is presented. This model can be used to predict the process performance in the ...

Integrating a logarithmic-strain based hyper-elastic formulation into a three-field mixed finite element formulation to deal with incompressibility in finite-strain elasto-plasticity

International Nuclear Information System (INIS)

Dina Al Akhrass; Bruchon, Julien; Drapier, Sylvain; Fayolle, Sebastien

2014-01-01

This paper deals with the treatment of incompressibility in solid mechanics in finite-strain elasto-plasticity. A finite-strain model proposed by Miehe, Apel and Lambrecht, which is based on a logarithmic strain measure and its work-conjugate stress tensor is chosen. Its main interest is that it allows for the adoption of standard constitutive models established in a small-strain framework. This model is extended to take into account the plastic incompressibility constraint intrinsically. In that purpose, an extension of this model to a three-field mixed finite element formulation is proposed, involving displacements, a strain variable and pressure as nodal variables with respect to standard finite element. Numerical examples of finite-strain problems are presented to assess the performance of the formulation. To conclude, an industrial case for which the classical under-integrated elements fail is considered. (authors)

A piecewise linear finite element discretization of the diffusion equation for arbitrary polyhedral grids

Energy Technology Data Exchange (ETDEWEB)

Bailey, Teresa S. [Texas A and M University, Department of Nuclear Engineering, College Station, TX 77843-3133 (United States)], E-mail: baileyte@tamu.edu; Adams, Marvin L. [Texas A and M University, Department of Nuclear Engineering, College Station, TX 77843-3133 (United States)], E-mail: mladams@tamu.edu; Yang, Brian [Lawrence Livermore National Laboratory, Livermore, CA 94551 (United States); Zika, Michael R. [Lawrence Livermore National Laboratory, Livermore, CA 94551 (United States)], E-mail: zika@llnl.gov

2008-04-01

We develop a piecewise linear (PWL) Galerkin finite element spatial discretization for the multi-dimensional radiation diffusion equation. It uses recently introduced piecewise linear weight and basis functions in the finite element approximation and it can be applied on arbitrary polygonal (2D) or polyhedral (3D) grids. We first demonstrate some analytical properties of the PWL method and perform a simple mode analysis to compare the PWL method with Palmer's vertex-centered finite-volume method and with a bilinear continuous finite element method. We then show that this new PWL method gives solutions comparable to those from Palmer's. However, since the PWL method produces a symmetric positive-definite coefficient matrix, it should be substantially more computationally efficient than Palmer's method, which produces an asymmetric matrix. We conclude that the Galerkin PWL method is an attractive option for solving diffusion equations on unstructured grids.

A piecewise linear finite element discretization of the diffusion equation for arbitrary polyhedral grids

International Nuclear Information System (INIS)

Bailey, Teresa S.; Adams, Marvin L.; Yang, Brian; Zika, Michael R.

2008-01-01

We develop a piecewise linear (PWL) Galerkin finite element spatial discretization for the multi-dimensional radiation diffusion equation. It uses recently introduced piecewise linear weight and basis functions in the finite element approximation and it can be applied on arbitrary polygonal (2D) or polyhedral (3D) grids. We first demonstrate some analytical properties of the PWL method and perform a simple mode analysis to compare the PWL method with Palmer's vertex-centered finite-volume method and with a bilinear continuous finite element method. We then show that this new PWL method gives solutions comparable to those from Palmer's. However, since the PWL method produces a symmetric positive-definite coefficient matrix, it should be substantially more computationally efficient than Palmer's method, which produces an asymmetric matrix. We conclude that the Galerkin PWL method is an attractive option for solving diffusion equations on unstructured grids

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

Distributed Finite Element Analysis Using a Transputer Network

Science.gov (United States)

Watson, James; Favenesi, James; Danial, Albert; Tombrello, Joseph; Yang, Dabby; Reynolds, Brian; Turrentine, Ronald; Shephard, Mark; Baehmann, Peggy

1989-01-01

The principal objective of this research effort was to demonstrate the extraordinarily cost effective acceleration of finite element structural analysis problems using a transputer-based parallel processing network. This objective was accomplished in the form of a commercially viable parallel processing workstation. The workstation is a desktop size, low-maintenance computing unit capable of supercomputer performance yet costs two orders of magnitude less. To achieve the principal research objective, a transputer based structural analysis workstation termed XPFEM was implemented with linear static structural analysis capabilities resembling commercially available NASTRAN. Finite element model files, generated using the on-line preprocessing module or external preprocessing packages, are downloaded to a network of 32 transputers for accelerated solution. The system currently executes at about one third Cray X-MP24 speed but additional acceleration appears likely. For the NASA selected demonstration problem of a Space Shuttle main engine turbine blade model with about 1500 nodes and 4500 independent degrees of freedom, the Cray X-MP24 required 23.9 seconds to obtain a solution while the transputer network, operated from an IBM PC-AT compatible host computer, required 71.7 seconds. Consequently, the $80,000 transputer network demonstrated a cost-performance ratio about 60 times better than the $15,000,000 Cray X-MP24 system.

Finite element simulations and experiments of ballistic impacts on high performance PE composite material

NARCIS (Netherlands)

Herlaar, K.; Jagt-Deutekom, M.J. van der; Jacobs, M.J.N.

2005-01-01

The use of lightweight composite armour concepts is essential for the protection of future combat systems, both vehicles and personal. The design of such armour systems is challenging due to the complex material behaviour. Finite element simulations can be used to help understand the important

Nonlinear finite element formulation for analyzing shape memory alloy cylindrical panels

International Nuclear Information System (INIS)

Mirzaeifar, R; Shakeri, M; Sadighi, M

2009-01-01

In this paper, a general incremental displacement based finite element formulation capable of modeling material nonlinearities based on first-order shear deformation theory (FSDT) is developed for cylindrical shape memory alloy (SMA) shells. The Boyd–Lagoudas phenomenological model with polynomial hardening in conjunction with 3D incremental convex cutting plane explicit algorithm is implemented for preparing the SMA constitutive model in the finite element formulation. Several numerical examples are presented for demonstrating the performance of the proposed formulation in stress, deflection and phase transformation analysis of pseudoelastic behavior of shape memory cylindrical panels with various boundary conditions. Also, it is shown that the presented formulation can be implemented for studying plates and beams with rectangular cross section

Computation of 2D compressible flows with a finite element method

International Nuclear Information System (INIS)

Montagne, J.L.

1981-04-01

When the homogeneous modelisation of the two phase flow is used the set of equations describing the flow is similar to an Euler system. Mixed finite elements are appropriate to discretize the equations. First, main properties of this kind of elements are reminded. Then, some properties of semi-implicite schemes on stability and entropy are given. Numerical tests have been performed, and the scheme gave satisfactory results

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

Compatible-strain mixed finite element methods for incompressible nonlinear elasticity

Science.gov (United States)

Faghih Shojaei, Mostafa; Yavari, Arash

2018-05-01

We introduce a new family of mixed finite elements for incompressible nonlinear elasticity - compatible-strain mixed finite element methods (CSFEMs). Based on a Hu-Washizu-type functional, we write a four-field mixed formulation with the displacement, the displacement gradient, the first Piola-Kirchhoff stress, and a pressure-like field as the four independent unknowns. Using the Hilbert complexes of nonlinear elasticity, which describe the kinematics and the kinetics of motion, we identify the solution spaces of the independent unknown fields. In particular, we define the displacement in H1, the displacement gradient in H (curl), the stress in H (div), and the pressure field in L2. The test spaces of the mixed formulations are chosen to be the same as the corresponding solution spaces. Next, in a conforming setting, we approximate the solution and the test spaces with some piecewise polynomial subspaces of them. Among these approximation spaces are the tensorial analogues of the Nédélec and Raviart-Thomas finite element spaces of vector fields. This approach results in compatible-strain mixed finite element methods that satisfy both the Hadamard compatibility condition and the continuity of traction at the discrete level independently of the refinement level of the mesh. By considering several numerical examples, we demonstrate that CSFEMs have a good performance for bending problems and for bodies with complex geometries. CSFEMs are capable of capturing very large strains and accurately approximating stress and pressure fields. Using CSFEMs, we do not observe any numerical artifacts, e.g., checkerboarding of pressure, hourglass instability, or locking in our numerical examples. Moreover, CSFEMs provide an efficient framework for modeling heterogeneous solids.

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.

Improvement of implicit finite element code performance in deep drawing simulations by dynamics contributions

NARCIS (Netherlands)

Meinders, Vincent T.; van den Boogaard, Antonius H.; Huetink, Han

2003-01-01

To intensify the use of implicit finite element codes for solving large scale problems, the computation time of these codes has to be decreased drastically. A method is developed which decreases the computational time of implicit codes by factors. The method is based on introducing inertia effects

Development of efficient finite elements for structural integrity analysis of solid rocket motor propellant grains

International Nuclear Information System (INIS)

Marimuthu, R.; Nageswara Rao, B.

2013-01-01

Solid propellant rocket motors (SRM) are regularly used in the satellite launch vehicles which consist of mainly three different structural materials viz., solid propellant, liner, and casing materials. It is essential to assess the structural integrity of solid propellant grains under the specified gravity, thermal and pressure loading conditions. For this purpose finite elements developed following the Herrmann formulation are: twenty node brick element (BH20), eight node quadrilateral plane strain element (PH8) and, eight node axi-symmetric solid of revolution element (AH8). The time-dependent nature of the solid propellant grains is taken into account utilizing the direct inverse method of Schepary to specify the effective Young's modulus and Poisson's ratio. The developed elements are tested considering various problems prior to implementation in the in-house software package (viz., Finite Element Analysis of STructures, FEAST). Several SRM configurations are analyzed to assess the structural integrity under different loading conditions. Finite element analysis results are found to be in good agreement with those obtained earlier from MARC software. -- Highlights: • Developed efficient Herrmann elements. • Accuracy of finite elements demonstrated solving several known solution problems. • Time dependent structural response obtained using the direct inverse method of Schepary. • Performed structural analysis of grains under gravity, thermal and pressure loads

Thermo-mechanical finite element analyses of bolted cask lid structures

International Nuclear Information System (INIS)

Wieser, G.; Qiao Linan; Eberle, A.; Voelzke, H.

2004-01-01

The analysis of complex bolted cask lid structures under mechanical or thermal accident conditions is important for the evaluation of cask integrity and leak-tightness in package design assessment according to the Transport Regulations or in aircraft crash scenarios. In this context BAM is developing methods based on Finite Elements to calculate the effects of mechanical impacts onto the bolted lid structures as well as effects caused by severe fire scenarios. I n case of fire it might be not enough to perform only a thermal heat transfer analysis. The complex cask design in connection with a severe hypothetical time-temperature-curve representing an accident fire scenario will create a strong transient heating up of the cask body and its lid system. This causes relative displacements between the seals and its counterparts that can be analyzed by a so-called thermo-mechanical calculation. Although it is currently not possible to correlate leakage rates with results from deformation analyses directly an appropriate Finite Element model of the considered type of metallic lid seal has been developed. For the present it is possible to estimate the behaviour of the seal based on the calculated relative displacements at its seating and the behaviour of the lid bolts under the impact load or the temperature field respectively. Except of the lid bolts the geometry of the cask and the mechanical loading is axial-symmetric which simplifies the analysis considerably and a two-dimensional Finite Element model with substitute lid bolts may be used. The substitute bolts are modelled as one-dimensional truss or beam elements. An advanced two-dimensional bolt submodel represents the bolts with plane stress continuum elements. This paper discusses the influence of different bolt modelling on the relative displacements at the seating of the seals. Besides this, the influence of bolt modelling, thermal properties and detail in geometry of the two-dimensional Finite Element models on

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.

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.

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, \

The Cost of Continuity: Performance of Iterative Solvers on Isogeometric Finite Elements

KAUST Repository

Collier, Nathan; Dalcin, Lisandro; Pardo, David; Calo, Victor M.

2013-01-01

In this paper we study how the use of a more continuous set of basis functions affects the cost of solving systems of linear equations resulting from a discretized Galerkin weak form. Specifically, we compare performance of linear solvers when discretizing using Co B-splines, which span traditional finite element spaces, and Cp-1 B-splines, which represent maximum continuity We provide theoretical estimates for the increase in cost of the matrix-vector product as well as for the construction and application of black-box preconditioners. We accompany these estimates with numerical results and study their sensitivity to various grid parameters such as element size h and polynomial order of approximation p in addition to the aforementioned continuity of the basis. Finally, we present timing results for a range of preconditioning options for the Laplace problem. We conclude that the matrix-vector product operation is at most 33p2/8 times more expensive for the more continuous space, although for moderately low p, this number is significantly reduced. Moreover, if static condensation is not employed, this number further reduces to at most a value of 8, even for high p. Preconditioning options can be up to p3 times more expensive to set up, although this difference significantly decreases for some popular preconditioners such as incomplete LU factorization. © 2013 Society for Industrial and Applied Mathematics.

The Cost of Continuity: Performance of Iterative Solvers on Isogeometric Finite Elements

KAUST Repository

Collier, Nathan

2013-03-19

In this paper we study how the use of a more continuous set of basis functions affects the cost of solving systems of linear equations resulting from a discretized Galerkin weak form. Specifically, we compare performance of linear solvers when discretizing using Co B-splines, which span traditional finite element spaces, and Cp-1 B-splines, which represent maximum continuity We provide theoretical estimates for the increase in cost of the matrix-vector product as well as for the construction and application of black-box preconditioners. We accompany these estimates with numerical results and study their sensitivity to various grid parameters such as element size h and polynomial order of approximation p in addition to the aforementioned continuity of the basis. Finally, we present timing results for a range of preconditioning options for the Laplace problem. We conclude that the matrix-vector product operation is at most 33p2/8 times more expensive for the more continuous space, although for moderately low p, this number is significantly reduced. Moreover, if static condensation is not employed, this number further reduces to at most a value of 8, even for high p. Preconditioning options can be up to p3 times more expensive to set up, although this difference significantly decreases for some popular preconditioners such as incomplete LU factorization. © 2013 Society for Industrial and Applied Mathematics.

Finite element analysis of permanent magnet motors

International Nuclear Information System (INIS)

Boglietti, A.; Chiampi, M.; Tartaglia, M.; Chiarabaglio, D.

1989-01-01

The analysis of permanent magnet D.C. brushless motors, supplied by current control inverters, is developed employing a finite element package tailored for such devices. The study is devoted to predicting the performance of a set of four poles machines, under different operating conditions (no-load, rated load). The over-load conditions are also considered including the saturation effect. Moreover the influence of such design parameters, as the tooth shape and the number of magnet segments, is investigated. Computed results are found in satisfactory agreement with experimental ones

The Mixed Finite Element Multigrid Method for Stokes Equations

Science.gov (United States)

Muzhinji, K.; Shateyi, S.; Motsa, S. S.

2015-01-01

The stable finite element discretization of the Stokes problem produces a symmetric indefinite system of linear algebraic equations. A variety of iterative solvers have been proposed for such systems in an attempt to construct efficient, fast, and robust solution techniques. This paper investigates one of such iterative solvers, the geometric multigrid solver, to find the approximate solution of the indefinite systems. The main ingredient of the multigrid method is the choice of an appropriate smoothing strategy. This study considers the application of different smoothers and compares their effects in the overall performance of the multigrid solver. We study the multigrid method with the following smoothers: distributed Gauss Seidel, inexact Uzawa, preconditioned MINRES, and Braess-Sarazin type smoothers. A comparative study of the smoothers shows that the Braess-Sarazin smoothers enhance good performance of the multigrid method. We study the problem in a two-dimensional domain using stable Hood-Taylor Q 2-Q 1 pair of finite rectangular elements. We also give the main theoretical convergence results. We present the numerical results to demonstrate the efficiency and robustness of the multigrid method and confirm the theoretical results. PMID:25945361

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.

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

A finite element model for estimating the techno-economic performance of the building-integrated photovoltaic blind

International Nuclear Information System (INIS)

Park, Hyo Seon; Koo, Choongwan; Hong, Taehoon; Oh, Jeongyoon; Jeong, Kwangbok

2016-01-01

Highlights: • A FEM was developed to estimate the techno-economic performance of the BIPB. • The mean absolute percentage error of the FEM_4_-_n_o_d_eBIPB was determined to be 4.54%. • In implementing the BIPB with the GC_i_n_c_l_._S_R_E_C plan, it was superior to the others. • Users can understand the operating mechanism of the proposed model (FEM_4_-_n_o_d_eBIPB). • The proposed model can be extended to any other country in the global environment. - Abstract: This study aims to develop the four-node-based finite element model for estimating the techno-economic performance of the building-integrated photovoltaic blind (FEM_4_-_n_o_d_eBIPB), which can be used by decision-maker in the early design phase. In developing the proposed model, this study uses various research methodologies such as energy simulation, finite element method, life cycle cost analysis, policy analysis, and visual basic application. Compared to the simulation results, the mean absolute percentage error of the proposed model was determined to be 4.54%, showing that the prediction accuracy of the proposed model was found to be excellent. Furthermore, the practical application was conducted for the ‘S’ elementary school facility in South Korea, which allows potential readers to easily and clearly understand the operating mechanism of the proposed model as well as its usability and extendability. The proposed model can be used to conduct the detailed analysis of the techno-economic performance of the BIPB by the type of utilization plan and to determine the optimal strategy for maximizing the value of the investment. Furthermore, the proposed research framework can be extended to any other technology, industry, and country in the global environment.

Assessment of Structural Behavior of Non-corroded and Corroded RCC Beams Using Finite Element Method

Directory of Open Access Journals (Sweden)

Anand Parande

2008-09-01

Full Text Available A three dimensional finite element model is developed to examine the structural behaviour of corroded reinforced concrete beam and non corroded reinforced concrete beam. Non linear finite element analysis is performed using the ANSYS program. SOLID 65, LINK 8 element represent concrete and discrete reinforcing steel bars, based on each component actual characteristics, non linear material properties are defined for both elements. The effect of corrosion in reinforced concrete is studied by finite element analysis; an approach is developed to model the corrosion product expansion causing concrete cover cracking for this, beam has been modeled using ANSYS and using this data the beam has been casted with M20 concrete after 28 days the beam will be tested for flexural strength. The comparison between ANSYS prediction and field data are made in terms of deflection, stress, strain, bond strength and crack pattern of concrete beam.

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....

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...

Elastoplastic finite element analysis for wet multidisc brake during lasting braking

Directory of Open Access Journals (Sweden)

Ji Zhanling

2015-01-01

Full Text Available Addressed to serious heat degradation problem of the braking continuously performed in the drag brake application for a long time, finite element analysis for bidirectional thermal-structure coupling is adopted to investigate temperature and stress when material properties are temperature-dependent. Based on the constitutive relations of heat transfer and strain-stress, three-dimensional transient finite element equilibrium equations with many kinds of boundary conditions for bidirectional thermal-structure coupling were derived. And it was originally presented that start time, location, severity and evolution laws of plastic deformation were depicted using dimensionless stress distribution contour with the yield limit related to temperature. The change laws of plastic element number and contact area versus braking time were expressed by plasticity ratio and contact ratio curves, respectively. The laws revealed by the numerical calculation results are in accordance with the objective perception and reasoning.

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

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

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...

Finite-element analysis and comparison of the AC loss performance of BSCCO and YBCO conductors

International Nuclear Information System (INIS)

Stavrev, Svetlomir; Grilli, Francesco; Dutoit, Bertrand; Ashworth, Stephen

2006-01-01

The AC loss performance of two BSCCO and two YBCO conductors of different geometry, characterized by the same self-field critical current of 150 A, is analysed and compared quantitatively. The comparison is made using the finite-element method with a nonlinear B-dependent E-J relation. A new shell-region model is utilised for the simulations of thin YBCO strips. Different AC working conditions are simulated: self-field, applied external field, and combined transport current and external perpendicular field application. Magnetic field and current density profiles are investigated in order to illustrate the reasons for the loss difference in the conductors. Depending on the application, the advantages of using BSCCO or YBCO conductors with specific geometry are outlined

Time-independent hybrid enrichment for finite element solution of transient conduction–radiation in diffusive grey media

Energy Technology Data Exchange (ETDEWEB)

Mohamed, M. Shadi, E-mail: m.s.mohamed@durham.ac.uk [School of Engineering and Computing Sciences, University of Durham, South Road, Durham DH1 3LE (United Kingdom); Seaid, Mohammed; Trevelyan, Jon [School of Engineering and Computing Sciences, University of Durham, South Road, Durham DH1 3LE (United Kingdom); Laghrouche, Omar [Institute for Infrastructure and Environment, Heriot-Watt University, Edinburgh EH14 4AS (United Kingdom)

2013-10-15

We investigate the effectiveness of the partition-of-unity finite element method for transient conduction–radiation problems in diffusive grey media. The governing equations consist of a semi-linear transient heat equation for the temperature field and a stationary diffusion approximation to the radiation in grey media. The coupled equations are integrated in time using a semi-implicit method in the finite element framework. We show that for the considered problems, a combination of hyperbolic and exponential enrichment functions based on an approximation of the boundary layer leads to improved accuracy compared to the conventional finite element method. It is illustrated that this approach can be more efficient than using h adaptivity to increase the accuracy of the finite element method near the boundary walls. The performance of the proposed partition-of-unity method is analyzed on several test examples for transient conduction–radiation problems in two space dimensions.

Studies of biaxial mechanical properties and nonlinear finite element modeling of skin.

Science.gov (United States)

Shang, Xituan; Yen, Michael R T; Gaber, M Waleed

2010-06-01

The objective of this research is to conduct mechanical property studies of skin from two individual but potentially connected aspects. One is to determine the mechanical properties of the skin experimentally by biaxial tests, and the other is to use the finite element method to model the skin properties. Dynamic biaxial tests were performed on 16 pieces of abdominal skin specimen from rats. Typical biaxial stress-strain responses show that skin possesses anisotropy, nonlinearity and hysteresis. To describe the stress-strain relationship in forms of strain energy function, the material constants of each specimen were obtained and the results show a high correlation between theory and experiments. Based on the experimental results, a finite element model of skin was built to model the skin's special properties including anisotropy and nonlinearity. This model was based on Arruda and Boyce's eight-chain model and Bischoff et al.'s finite element model of skin. The simulation results show that the isotropic, nonlinear eight-chain model could predict the skin's anisotropic and nonlinear responses to biaxial loading by the presence of an anisotropic prestress state.

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

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...

The Dynamic Response of an Euler-Bernoulli Beam on an Elastic Foundation by Finite Element Analysis using the Exact Stiffness Matrix

International Nuclear Information System (INIS)

Kim, Jeong Soo; Kim, Moon Kyum

2012-01-01

In this study, finite element analysis of beam on elastic foundation, which received great attention of researchers due to its wide applications in engineering, is performed for estimating dynamic responses of shallow foundation using exact stiffness matrix. First, element stiffness matrix based on the closed solution of beam on elastic foundation is derived. Then, we performed static finite element analysis included exact stiffness matrix numerically, comparing results from the analysis with some exact analysis solutions well known for verification. Finally, dynamic finite element analysis is performed for a shallow foundation structure under rectangular pulse loading using trapezoidal method. The dynamic analysis results exist in the reasonable range comparing solution of single degree of freedom problem under a similar condition. The results show that finite element analysis using exact stiffness matrix is evaluated as a good tool of estimating the dynamic response of structures on elastic foundation.

Finite Element Analysis of Reinforced Concrete Beam-Column Connections with Governing Joint Shear Failure Mode

Directory of Open Access Journals (Sweden)

M.A. Najafgholipour

Full Text Available Abstract Reinforced concrete (RC beam-column connections especially those without transverse reinforcement in joint region can exhibit brittle behavior when intensive damage is concentrated in the joint region during an earthquake event. Brittle behavior in the joint region can compromise the ductile design philosophy and the expected overall performance of structure when subjected to seismic loading. Considering the importance of joint shear failure influences on strength, ductility and stability of RC moment resisting frames, a finite element modeling which focuses on joint shear behavior is presented in this article. Nonlinear finite element analysis (FEA of RC beam-column connections is performed in order to investigate the joint shear failure mode in terms of joint shear capacity, deformations and cracking pattern. A 3D finite element model capable of appropriately modeling the concrete stress-strain behavior, tensile cracking and compressive damage of concrete and indirect modeling of steel-concrete bond is used. In order to define nonlinear behavior of concrete material, the concrete damage plasticity is applied to the numerical model as a distributed plasticity over the whole geometry. Finite element model is then verified against experimental results of two non-ductile beam-column connections (one exterior and one interior which are vulnerable to joint shear failure. The comparison between experimental and numerical results indicates that the FE model is able to simulate the performance of the beam-column connections and is able to capture the joint shear failure in RC beam-column connections.

Finite element simulations with ANSYS workbench 17 theory, applications, case studies

CERN Document Server

Lee, Huei-Huang

2017-01-01

Finite Element Simulations with ANSYS Workbench 17 is a comprehensive and easy to understand workbook. Printed in full color, it utilizes rich graphics and step-by-step instructions to guide you through learning how to perform finite element simulations using ANSYS Workbench. Twenty seven real world case studies are used throughout the book. Many of these case studies are industrial or research projects that you build from scratch. Prebuilt project files are available for download should you run into any problems. Companion videos, that demonstrate exactly how to perform each tutorial, are also available. 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 though this entire boo...

Finite element simulations with ANSYS Workbench 18 theory, applications, case studies

CERN Document Server

Lee,\tHuei-huang

2018-01-01

Finite Element Simulations with ANSYS Workbench 18 is a comprehensive and easy to understand workbook. Printed in full color, it utilizes rich graphics and step-by-step instructions to guide you through learning how to perform finite element simulations using ANSYS Workbench. Twenty seven real world case studies are used throughout the book. Many of these case studies are industrial or research projects that you build from scratch. Prebuilt project files are available for download should you run into any problems. Companion videos, that demonstrate exactly how to perform each tutorial, are also available. 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 is utilized though this entire...

Simulation of three-dimensional, time-dependent, incompressible flows by a finite element method

International Nuclear Information System (INIS)

Chan, S.T.; Gresho, P.M.; Lee, R.L.; Upson, C.D.

1981-01-01

A finite element model has been developed for simulating the dynamics of problems encountered in atmospheric pollution and safety assessment studies. The model is based on solving the set of three-dimensional, time-dependent, conservation equations governing incompressible flows. Spatial discretization is performed via a modified Galerkin finite element method, and time integration is carried out via the forward Euler method (pressure is computed implicitly, however). Several cost-effective techniques (including subcycling, mass lumping, and reduced Gauss-Legendre quadrature) which have been implemented are discussed. Numerical results are presented to demonstrate the applicability of the model

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 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.

Thermal buckling comparative analysis using Different FE (Finite Element) tools

Energy Technology Data Exchange (ETDEWEB)

Banasiak, Waldemar; Labouriau, Pedro [INTECSEA do Brasil, Rio de Janeiro, RJ (Brazil); Burnett, Christopher [INTECSEA UK, Surrey (United Kingdom); Falepin, Hendrik [Fugro Engineers SA/NV, Brussels (Belgium)

2009-12-19

High operational temperature and pressure in offshore pipelines may lead to unexpected lateral movements, sometimes call lateral buckling, which can have serious consequences for the integrity of the pipeline. The phenomenon of lateral buckling in offshore pipelines needs to be analysed in the design phase using FEM. The analysis should take into account many parameters, including operational temperature and pressure, fluid characteristic, seabed profile, soil parameters, coatings of the pipe, free spans etc. The buckling initiation force is sensitive to small changes of any initial geometric out-of-straightness, thus the modeling of the as-laid state of the pipeline is an important part of the design process. Recently some dedicated finite elements programs have been created making modeling of the offshore environment more convenient that has been the case with the use of general purpose finite element software. The present paper aims to compare thermal buckling analysis of sub sea pipeline performed using different finite elements tools, i.e. general purpose programs (ANSYS, ABAQUS) and dedicated software (SAGE Profile 3D) for a single pipeline resting on an the seabed. The analyses considered the pipeline resting on a flat seabed with a small levels of out-of straightness initiating the lateral buckling. The results show the quite good agreement of results of buckling in elastic range and in the conclusions next comparative analyses with sensitivity cases are recommended. (author)

Finite element analysis of the cyclic indentation of bilayer enamel

International Nuclear Information System (INIS)

Jia, Yunfei; Xuan, Fu-zhen; Chen, Xiaoping; Yang, Fuqian

2014-01-01

Tooth enamel is often subjected to repeated contact and often experiences contact deformation in daily life. The mechanical strength of the enamel determines the biofunctionality of the tooth. Considering the variation of the rod arrangement in outer and inner enamel, we approximate enamel as a bilayer structure and perform finite element analysis of the cyclic indentation of the bilayer structure, to mimic the repeated contact of enamel during mastication. The dynamic deformation behaviour of both the inner enamel and the bilayer enamel is examined. The material parameters of the inner and outer enamel used in the analysis are obtained by fitting the finite element results with the experimental nanoindentation results. The penetration depth per cycle at the quasi-steady state is used to describe the depth propagation speed, which exhibits a two-stage power-law dependence on the maximum indentation load and the amplitude of the cyclic load, respectively. The continuous penetration of the indenter reflects the propagation of the plastic zone during cyclic indentation, which is related to the energy dissipation. The outer enamel serves as a protective layer due to its great resistance to contact deformation in comparison to the inner enamel. The larger equivalent plastic strain and lower stresses in the inner enamel during cyclic indentation, as calculated from the finite element analysis, indicate better crack/fracture resistance of the inner enamel. (paper)

Finite element analysis of the cyclic indentation of bilayer enamel

Science.gov (United States)

Jia, Yunfei; Xuan, Fu-zhen; Chen, Xiaoping; Yang, Fuqian

2014-04-01

Tooth enamel is often subjected to repeated contact and often experiences contact deformation in daily life. The mechanical strength of the enamel determines the biofunctionality of the tooth. Considering the variation of the rod arrangement in outer and inner enamel, we approximate enamel as a bilayer structure and perform finite element analysis of the cyclic indentation of the bilayer structure, to mimic the repeated contact of enamel during mastication. The dynamic deformation behaviour of both the inner enamel and the bilayer enamel is examined. The material parameters of the inner and outer enamel used in the analysis are obtained by fitting the finite element results with the experimental nanoindentation results. The penetration depth per cycle at the quasi-steady state is used to describe the depth propagation speed, which exhibits a two-stage power-law dependence on the maximum indentation load and the amplitude of the cyclic load, respectively. The continuous penetration of the indenter reflects the propagation of the plastic zone during cyclic indentation, which is related to the energy dissipation. The outer enamel serves as a protective layer due to its great resistance to contact deformation in comparison to the inner enamel. The larger equivalent plastic strain and lower stresses in the inner enamel during cyclic indentation, as calculated from the finite element analysis, indicate better crack/fracture resistance of the inner enamel.

Generalized multiscale finite element methods (GMsFEM)

KAUST Repository

Efendiev, Yalchin R.; Galvis, Juan; Hou, Thomasyizhao

2013-01-01

In this paper, we propose a general approach called Generalized Multiscale Finite Element Method (GMsFEM) for performing multiscale simulations for problems without scale separation over a complex input space. As in multiscale finite element methods (MsFEMs), the main idea of the proposed approach is to construct a small dimensional local solution space that can be used to generate an efficient and accurate approximation to the multiscale solution with a potentially high dimensional input parameter space. In the proposed approach, we present a general procedure to construct the offline space that is used for a systematic enrichment of the coarse solution space in the online stage. The enrichment in the online stage is performed based on a spectral decomposition of the offline space. In the online stage, for any input parameter, a multiscale space is constructed to solve the global problem on a coarse grid. The online space is constructed via a spectral decomposition of the offline space and by choosing the eigenvectors corresponding to the largest eigenvalues. The computational saving is due to the fact that the construction of the online multiscale space for any input parameter is fast and this space can be re-used for solving the forward problem with any forcing and boundary condition. Compared with the other approaches where global snapshots are used, the local approach that we present in this paper allows us to eliminate unnecessary degrees of freedom on a coarse-grid level. We present various examples in the paper and some numerical results to demonstrate the effectiveness of our method. © 2013 Elsevier Inc.

Generalized multiscale finite element methods (GMsFEM)

KAUST Repository

Efendiev, Yalchin R.

2013-10-01

In this paper, we propose a general approach called Generalized Multiscale Finite Element Method (GMsFEM) for performing multiscale simulations for problems without scale separation over a complex input space. As in multiscale finite element methods (MsFEMs), the main idea of the proposed approach is to construct a small dimensional local solution space that can be used to generate an efficient and accurate approximation to the multiscale solution with a potentially high dimensional input parameter space. In the proposed approach, we present a general procedure to construct the offline space that is used for a systematic enrichment of the coarse solution space in the online stage. The enrichment in the online stage is performed based on a spectral decomposition of the offline space. In the online stage, for any input parameter, a multiscale space is constructed to solve the global problem on a coarse grid. The online space is constructed via a spectral decomposition of the offline space and by choosing the eigenvectors corresponding to the largest eigenvalues. The computational saving is due to the fact that the construction of the online multiscale space for any input parameter is fast and this space can be re-used for solving the forward problem with any forcing and boundary condition. Compared with the other approaches where global snapshots are used, the local approach that we present in this paper allows us to eliminate unnecessary degrees of freedom on a coarse-grid level. We present various examples in the paper and some numerical results to demonstrate the effectiveness of our method. © 2013 Elsevier Inc.

A finite strain Eulerian formulation for compressible and nearly incompressible hyperelasticity using high-order B-spline finite elements

KAUST Repository

Duddu, Ravindra

2011-10-05

We present a numerical formulation aimed at modeling the nonlinear response of elastic materials using large deformation continuum mechanics in three dimensions. This finite element formulation is based on the Eulerian description of motion and the transport of the deformation gradient. When modeling a nearly incompressible solid, the transport of the deformation gradient is decomposed into its isochoric part and the Jacobian determinant as independent fields. A homogeneous isotropic hyperelastic solid is assumed and B-splines-based finite elements are used for the spatial discretization. A variational multiscale residual-based approach is employed to stabilize the transport equations. The performance of the scheme is explored for both compressible and nearly incompressible applications. The numerical results are in good agreement with theory illustrating the viability of the computational scheme. © 2011 John Wiley & Sons, Ltd.

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

Failure analysis of natural gas buried X65 steel pipeline under deflection load using finite element method

International Nuclear Information System (INIS)

Liu, P.F.; Zheng, J.Y.; Zhang, B.J.; Shi, P.

2010-01-01

A 3D parametric finite element model of the pipeline and soil is established using finite element method to perform the failure analysis of natural gas buried X65 steel pipeline under deflection load. The pipeline is assumed to be loaded in a parabolic deflection displacement along the axial direction. Based on the true stress-strain constitutive relationship of X65 steel, the elastic-plastic finite element analysis employs the arc-length algorithm and non-linear stabilization algorithm respectively to simulate the strain softening properties of pipeline after plastic collapse. Besides, effects of the soil types and model sizes on the maximum deflection displacement of pipeline are investigated. The proposed finite element method serves as a base available for the safety design and evaluation as well as engineering acceptance criterion for the failure of pipeline due to deflection.

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...

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

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 ...

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.

Model Reduction in Dynamic Finite Element Analysis of Lightweight Structures

DEFF Research Database (Denmark)

Flodén, Ola; Persson, Kent; Sjöström, Anders

2012-01-01

models may be created by assembling models of floor and wall structures into large models of complete buildings. When assembling the floor and wall models, the number of degrees of freedom quickly increases to exceed the limits of computer capacity, at least in a reasonable amount of computational time...... Hz. Three different methods of model reduction were investigated; Guyan reduction, component mode synthesis and a third approach where a new finite element model was created with structural elements. Eigenvalue and steady-state analyses were performed in order to compare the errors...

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.))

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

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)

Aranha: a 2D mesh generator for triangular finite elements

International Nuclear Information System (INIS)

Fancello, E.A.; Salgado, A.C.; Feijoo, R.A.

1990-01-01

A method for generating unstructured meshes for linear and quadratic triangular finite elements is described in this paper. Some topics on the C language data structure used in the development of the program Aranha are also presented. The applicability for adaptive remeshing is shown and finally several examples are included to illustrate the performance of the method in irregular connected planar domains. (author)

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

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

Fiber-reinforced materials: finite elements for the treatment of the inextensibility constraint

Science.gov (United States)

Auricchio, Ferdinando; Scalet, Giulia; Wriggers, Peter

2017-12-01

The present paper proposes a numerical framework for the analysis of problems involving fiber-reinforced anisotropic materials. Specifically, isotropic linear elastic solids, reinforced by a single family of inextensible fibers, are considered. The kinematic constraint equation of inextensibility in the fiber direction leads to the presence of an undetermined fiber stress in the constitutive equations. To avoid locking-phenomena in the numerical solution due to the presence of the constraint, mixed finite elements based on the Lagrange multiplier, perturbed Lagrangian, and penalty method are proposed. Several boundary-value problems under plane strain conditions are solved and numerical results are compared to analytical solutions, whenever the derivation is possible. The performed simulations allow to assess the performance of the proposed finite elements and to discuss several features of the developed formulations concerning the effective approximation for the displacement and fiber stress fields, mesh convergence, and sensitivity to penalty parameters.

Discrete/Finite Element Modelling of Rock Cutting with a TBM Disc Cutter

Science.gov (United States)

Labra, Carlos; Rojek, Jerzy; Oñate, Eugenio

2017-03-01

This paper presents advanced computer simulation of rock cutting process typical for excavation works in civil engineering. Theoretical formulation of the hybrid discrete/finite element model has been presented. The discrete and finite element methods have been used in different subdomains of a rock sample according to expected material behaviour, the part which is fractured and damaged during cutting is discretized with the discrete elements while the other part is treated as a continuous body and it is modelled using the finite element method. In this way, an optimum model is created, enabling a proper representation of the physical phenomena during cutting and efficient numerical computation. The model has been applied to simulation of the laboratory test of rock cutting with a single TBM (tunnel boring machine) disc cutter. The micromechanical parameters have been determined using the dimensionless relationships between micro- and macroscopic parameters. A number of numerical simulations of the LCM test in the unrelieved and relieved cutting modes have been performed. Numerical results have been compared with available data from in-situ measurements in a real TBM as well as with the theoretical predictions showing quite a good agreement. The numerical model has provided a new insight into the cutting mechanism enabling us to investigate the stress and pressure distribution at the tool-rock interaction. Sensitivity analysis of rock cutting performed for different parameters including disc geometry, cutting velocity, disc penetration and spacing has shown that the presented numerical model is a suitable tool for the design and optimization of rock cutting process.

Domain decomposition solvers for nonlinear multiharmonic finite element equations

KAUST Repository

Copeland, D. M.

2010-01-01

In many practical applications, for instance, in computational electromagnetics, the excitation is time-harmonic. Switching from the time domain to the frequency domain allows us to replace the expensive time-integration procedure by the solution of a simple elliptic equation for the amplitude. This is true for linear problems, but not for nonlinear problems. However, due to the periodicity of the solution, we can expand the solution in a Fourier series. Truncating this Fourier series and approximating the Fourier coefficients by finite elements, we arrive at a large-scale coupled nonlinear system for determining the finite element approximation to the Fourier coefficients. The construction of fast solvers for such systems is very crucial for the efficiency of this multiharmonic approach. In this paper we look at nonlinear, time-harmonic potential problems as simple model problems. We construct and analyze almost optimal solvers for the Jacobi systems arising from the Newton linearization of the large-scale coupled nonlinear system that one has to solve instead of performing the expensive time-integration procedure. © 2010 de Gruyter.

Heat Conduction Analysis Using Semi Analytical Finite Element Method

International Nuclear Information System (INIS)

Wargadipura, A. H. S.

1997-01-01

Heat conduction problems are very often found in science and engineering fields. It is of accrual importance to determine quantitative descriptions of this important physical phenomena. This paper discusses the development and application of a numerical formulation and computation that can be used to analyze heat conduction problems. The mathematical equation which governs the physical behaviour of heat conduction is in the form of second order partial differential equations. The numerical resolution used in this paper is performed using the finite element method and Fourier series, which is known as semi-analytical finite element methods. The numerical solution results in simultaneous algebraic equations which is solved using the Gauss elimination methodology. The computer implementation is carried out using FORTRAN language. In the final part of the paper, a heat conduction problem in a rectangular plate domain with isothermal boundary conditions in its edge is solved to show the application of the computer program developed and also a comparison with analytical solution is discussed to assess the accuracy of the numerical solution obtained

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

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.

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.``

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...

Finite element modelling for mode-I fracture behaviour of CFRP

Science.gov (United States)

Chetan, H. C.; Kattimani, Subhaschandra; Murigendrappa, S. M.

2018-04-01

Debonding is a major failure mechanism in Carbon Fiber Reinforced Polymer (CFRP) due to presence of many adhesion joins, in between many layers. In the current study a finite element simulation is carried out using Virtual Crack Closure Technique (VCCT) and Cohesive Zone Modelling (CZM) using Abaqus as analysis tool. A comparative study is performed in to order analyze convergence of results from CZM and VCCT. It was noted that CZM results matched well with published literature. The results from VCCT were also in good comparison with experimental data of published literature, but were seen to be overestimated. Parametric study is performed to evaluate the variation of input parameters like initial stiffness, element size, peak stress and energy release rate `G'. From the numerical evaluation, it was noted that CZM simulation relies largely on element size and peak stress.

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.

Finite Macro-Element Mesh Deformation in a Structured Multi-Block Navier-Stokes Code

Science.gov (United States)

Bartels, Robert E.

2005-01-01

A mesh deformation scheme is developed for a structured multi-block Navier-Stokes code consisting of two steps. The first step is a finite element solution of either user defined or automatically generated macro-elements. Macro-elements are hexagonal finite elements created from a subset of points from the full mesh. When assembled, the finite element system spans the complete flow domain. Macro-element moduli vary according to the distance to the nearest surface, resulting in extremely stiff elements near a moving surface and very pliable elements away from boundaries. Solution of the finite element system for the imposed boundary deflections generally produces smoothly varying nodal deflections. The manner in which distance to the nearest surface has been found to critically influence the quality of the element deformation. The second step is a transfinite interpolation which distributes the macro-element nodal deflections to the remaining fluid mesh points. The scheme is demonstrated for several two-dimensional applications.

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.

The aetiology of pressure sores: combining animal experiments and finite element modelling

NARCIS (Netherlands)

Bosboom, E.M.H.; Oomens, C.W.J.; Bouten, C.V.C.; Janssen, J.D.; Kuipers, H.

1999-01-01

Animal experiments were performed to study the relationship between strictly controlled external mechanical loads and the onset of tissue damage. The finite element model was used to determine the local stresses and strains within the tissues during the experiments. By comparison of the results of

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

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

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.)

Probabilistic finite elements for transient analysis in nonlinear continua

Science.gov (United States)

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

1985-01-01

The probabilistic finite element method (PFEM), which is a combination of finite element methods and second-moment analysis, is formulated for linear and nonlinear continua with inhomogeneous random fields. Analogous to the discretization of the displacement field in finite element methods, the random field is also discretized. The formulation is simplified by transforming the correlated variables to a set of uncorrelated variables through an eigenvalue orthogonalization. Furthermore, it is shown that a reduced set of the uncorrelated variables is sufficient for the second-moment analysis. Based on the linear formulation of the PFEM, the method is then extended to transient analysis in nonlinear continua. The accuracy and efficiency of the method is demonstrated by application to a one-dimensional, elastic/plastic wave propagation problem. The moments calculated compare favorably with those obtained by Monte Carlo simulation. Also, the procedure is amenable to implementation in deterministic FEM based computer programs.

Adaptive finite-element ballooning analysis of bipolar ionized fields

International Nuclear Information System (INIS)

Al-Hamouz, Z.M.

1995-01-01

This paper presents an adaptive finite-element iterative method for the analysis of the ionized field around high-voltage bipolar direct-current (HVDC) transmission line conductors without resort to Deutsch's assumption. A new iterative finite-element ballooning technique is proposed to solve Poisson's equation wherein the commonly used artificial boundary around the transmission line conductors is simulated at infinity. Unlike all attempts reported in the literature for the solution of ionized field, the constancy of the conductors' surface field at the corona onset value is directly implemented in the finite-element formulation. In order to investigate the effectiveness of the proposed method, a laboratory model was built. It has been found that the calculated V-I characteristics and the ground-plane current density agreed well with those measured experimentally. The simplicity in computer programming in addition to the low number of iterations required to achieve convergence characterize this method of analysis

A simple nodal force distribution method in refined finite element meshes

Energy Technology Data Exchange (ETDEWEB)

Park, Jai Hak [Chungbuk National University, Chungju (Korea, Republic of); Shin, Kyu In [Gentec Co., Daejeon (Korea, Republic of); Lee, Dong Won [Korea Atomic Energy Research Institute, Daejeon (Korea, Republic of); Cho, Seungyon [National Fusion Research Institute, Daejeon (Korea, Republic of)

2017-05-15

In finite element analyses, mesh refinement is frequently performed to obtain accurate stress or strain values or to accurately define the geometry. After mesh refinement, equivalent nodal forces should be calculated at the nodes in the refined mesh. If field variables and material properties are available at the integration points in each element, then the accurate equivalent nodal forces can be calculated using an adequate numerical integration. However, in certain circumstances, equivalent nodal forces cannot be calculated because field variable data are not available. In this study, a very simple nodal force distribution method was proposed. Nodal forces of the original finite element mesh are distributed to the nodes of refined meshes to satisfy the equilibrium conditions. The effect of element size should also be considered in determining the magnitude of the distributing nodal forces. A program was developed based on the proposed method, and several example problems were solved to verify the accuracy and effectiveness of the proposed method. From the results, accurate stress field can be recognized to be obtained from refined meshes using the proposed nodal force distribution method. In example problems, the difference between the obtained maximum stress and target stress value was less than 6 % in models with 8-node hexahedral elements and less than 1 % in models with 20-node hexahedral elements or 10-node tetrahedral elements.

Straightened cervical lordosis causes stress concentration: a finite element model study

Energy Technology Data Exchange (ETDEWEB)

Wei, Wei; Shi, Shiyuan; Fei, Jun; Wang, Yifan; Chen, Chunyue [Hangzhou Red Cross Hospital, Hangzhou, Zhejiang, (China); Liao, Shenhui [School of Information Science and Engineering, Central South University, Changsha, Hunan (China)

2013-03-15

In this study, we propose a finite element analysis of the complete cervical spine with straightened and normal physiological curvature by using a specially designed modelling system. An accurate finite element model is established to recommend plausible approaches to treatment of cervical spondylosis through the finite element analysis results. There are few reports of biomechanics influence of the straightened cervical curve. It is difficult to measure internal responses of cervical spine directly. However, the finite element method has been reported to have the capability to quantify both external and internal responses to mechanical loading, such as the strain and stress distribution of spinal components. We choose a subject with a straightened cervical spine from whom to collect the CT scan data, which formed the basis of the finite element analysis. By using a specially designed modelling system, a high quality finite element model of the complete cervical spine with straightened curvature was generated, which was then mapped to reconstruct a normal physiological curvature model by a volumetric mesh deformation method based on discrete differential properties. Then, the same boundary conditions were applied to do a comparison. The result demonstrated that the active movement range of straightened cervical spine decreased by 24–33 %, but the stress increased by 5–95 %. The stress was concentrated at the facet joint cartilage, uncovertebral joint and the disk. The results suggest that cervical lordosis may have a direct impact on cervical spondylosis treatment. These results may be useful for clinical treatment of cervical spondylosis with straightened curvature.

Straightened cervical lordosis causes stress concentration: a finite element model study

International Nuclear Information System (INIS)

Wei, Wei; Shi, Shiyuan; Fei, Jun; Wang, Yifan; Chen, Chunyue; Liao, Shenhui

2013-01-01

In this study, we propose a finite element analysis of the complete cervical spine with straightened and normal physiological curvature by using a specially designed modelling system. An accurate finite element model is established to recommend plausible approaches to treatment of cervical spondylosis through the finite element analysis results. There are few reports of biomechanics influence of the straightened cervical curve. It is difficult to measure internal responses of cervical spine directly. However, the finite element method has been reported to have the capability to quantify both external and internal responses to mechanical loading, such as the strain and stress distribution of spinal components. We choose a subject with a straightened cervical spine from whom to collect the CT scan data, which formed the basis of the finite element analysis. By using a specially designed modelling system, a high quality finite element model of the complete cervical spine with straightened curvature was generated, which was then mapped to reconstruct a normal physiological curvature model by a volumetric mesh deformation method based on discrete differential properties. Then, the same boundary conditions were applied to do a comparison. The result demonstrated that the active movement range of straightened cervical spine decreased by 24–33 %, but the stress increased by 5–95 %. The stress was concentrated at the facet joint cartilage, uncovertebral joint and the disk. The results suggest that cervical lordosis may have a direct impact on cervical spondylosis treatment. These results may be useful for clinical treatment of cervical spondylosis with straightened curvature.

Finite-element analysis of NiTi wire deflection during orthodontic levelling treatment

International Nuclear Information System (INIS)

Razali, M F; Mahmud, A S; Abdullah, J; Mokhtar, N

2016-01-01

Finite-element analysis is an important product development tool in medical devices industry for design and failure analysis of devices. This tool helps device designers to quickly explore various design options, optimizing specific designs and providing a deeper insight how a device is actually performing. In this study, three-dimensional finite-element models of superelastic nickel-titanium arch wire engaged in a three brackets system were developed. The aim was to measure the effect of binding friction developed on wire-bracket interaction towards the remaining recovery force available for tooth movement. Uniaxial and three brackets bending test were modelled and validated against experimental works. The prediction made by the three brackets bending models shows good agreement with the experimental results. (paper)

Finite-element analysis of NiTi wire deflection during orthodontic levelling treatment

Science.gov (United States)

Razali, M. F.; Mahmud, A. S.; Mokhtar, N.; Abdullah, J.

2016-02-01

Finite-element analysis is an important product development tool in medical devices industry for design and failure analysis of devices. This tool helps device designers to quickly explore various design options, optimizing specific designs and providing a deeper insight how a device is actually performing. In this study, three-dimensional finite-element models of superelastic nickel-titanium arch wire engaged in a three brackets system were developed. The aim was to measure the effect of binding friction developed on wire-bracket interaction towards the remaining recovery force available for tooth movement. Uniaxial and three brackets bending test were modelled and validated against experimental works. The prediction made by the three brackets bending models shows good agreement with the experimental results.

Parallel algorithms for solving the diffusion equation by finite elements methods and by nodal methods

International Nuclear Information System (INIS)

Coulomb, F.

1989-06-01

The aim of this work is to study methods for solving the diffusion equation, based on a primal or mixed-dual finite elements discretization and well suited for use on multiprocessors computers; domain decomposition methods are the subject of the main part of this study, the linear systems being solved by the block-Jacobi method. The origin of the diffusion equation is explained in short, and various variational formulations are reminded. A survey of iterative methods is given. The elemination of the flux or current is treated in the case of a mixed method. Numerical tests are performed on two examples of reactors, in order to compare mixed elements and Lagrange elements. A theoretical study of domain decomposition is led in the case of Lagrange finite elements, and convergence conditions for the block-Jacobi method are derived; the dissection decomposition is previously the purpose of a particular numerical analysis. In the case of mixed-dual finite elements, a study is led on examples and is confirmed by numerical tests performed for the dissection decomposition; furthermore, after being justified, decompositions along axes of symmetry are numerically tested. In the case of a decomposition into two subdomains, the dissection decomposition and the decomposition with an integrated interface are compared. Alternative directions methods are defined; the convergence of those relative to Lagrange elements is shown; in the case of mixed elements, convergence conditions are found [fr

Finite element model for heat conduction in jointed rock masses

International Nuclear Information System (INIS)

Gartling, D.K.; Thomas, R.K.

1981-01-01

A computatonal procedure for simulating heat conduction in a fractured rock mass is proposed and illustrated in the present paper. The method makes use of a simple local model for conduction in the vicinity of a single open fracture. The distributions of fractures and fracture properties within the finite element model are based on a statistical representation of geologic field data. Fracture behavior is included in the finite element computation by locating local, discrete fractures at the element integration points

Mixed finite element - discontinuous finite volume element discretization of a general class of multicontinuum models

Science.gov (United States)

Ruiz-Baier, Ricardo; Lunati, Ivan

2016-10-01

We present a novel discretization scheme tailored to a class of multiphase models that regard the physical system as consisting of multiple interacting continua. In the framework of mixture theory, we consider a general mathematical model that entails solving a system of mass and momentum equations for both the mixture and one of the phases. The model results in a strongly coupled and nonlinear system of partial differential equations that are written in terms of phase and mixture (barycentric) velocities, phase pressure, and saturation. We construct an accurate, robust and reliable hybrid method that combines a mixed finite element discretization of the momentum equations with a primal discontinuous finite volume-element discretization of the mass (or transport) equations. The scheme is devised for unstructured meshes and relies on mixed Brezzi-Douglas-Marini approximations of phase and total velocities, on piecewise constant elements for the approximation of phase or total pressures, as well as on a primal formulation that employs discontinuous finite volume elements defined on a dual diamond mesh to approximate scalar fields of interest (such as volume fraction, total density, saturation, etc.). As the discretization scheme is derived for a general formulation of multicontinuum physical systems, it can be readily applied to a large class of simplified multiphase models; on the other, the approach can be seen as a generalization of these models that are commonly encountered in the literature and employed when the latter are not sufficiently accurate. An extensive set of numerical test cases involving two- and three-dimensional porous media are presented to demonstrate the accuracy of the method (displaying an optimal convergence rate), the physics-preserving properties of the mixed-primal scheme, as well as the robustness of the method (which is successfully used to simulate diverse physical phenomena such as density fingering, Terzaghi's consolidation

OPTIMIZATION OF I-SECTION PROFILE DESIGN BY THE FINITE ELEMENT METHOD

Directory of Open Access Journals (Sweden)

Patryk Różyło

2016-03-01

Full Text Available This paper discusses the problem of design optimization for an I-section profile. The optimization process was performed using the Abaqus program. The numerical analysis of a strictly static problem was based on the finite element method. The scope of the analysis involved both determination of stresses and displacements in the profile and structure topology optimization. The main focus of the numerical analysis was put on reducing profile volume while maintaining the same load and similar stresses prior to and after optimization. The solution of the optimization problem is just an example of the potential of using this method in combination with the finite element method in the Abaqus environment. Nowadays numerical analysis is the most effective cost-reducing alternative to experimental tests and it enables structure examination by means of a computer.

hpGEM -- A software framework for discontinuous Galerkin finite element methods

NARCIS (Netherlands)

Pesch, L.; Bell, A.; Sollie, W.E.H.; Ambati, V.R.; Bokhove, Onno; van der Vegt, Jacobus J.W.

2006-01-01

hpGEM, a novel framework for the implementation of discontinuous Galerkin finite element methods, is described. We present structures and methods that are common for many (discontinuous) finite element methods and show how we have implemented the components as an object-oriented framework. This

Orthodontic treatment: Introducing finite element analysis

NARCIS (Netherlands)

Driel, W.D. van; Leeuwen, E.J. van

1998-01-01

The aim of orthodontic treatment is the displacement of teeth by means ofspecial appliances, like braces and brackets. Through these appliances the orthodontist can apply a set of forces to the teeth which wilt result in its displacement through the jawbone. Finite Element analysis of this process

A finite element modeling method for predicting long term corrosion rates

International Nuclear Information System (INIS)

Fu, J.W.; Chan, S.

1984-01-01

For the analyses of galvanic corrosion, pitting and crevice corrosion, which have been identified as possible corrosion processes for nuclear waste isolation, a finite element method has been developed for the prediction of corrosion rates. The method uses a finite element mesh to model the corrosive environment and the polarization curves of metals are assigned as the boundary conditions to calculate the corrosion cell current distribution. A subroutine is used to calculate the chemical change with time in the crevice or the pit environments. In this paper, the finite element method is described along with experimental confirmation

Domain decomposition based iterative methods for nonlinear elliptic finite element problems

Energy Technology Data Exchange (ETDEWEB)

Cai, X.C. [Univ. of Colorado, Boulder, CO (United States)

1994-12-31

The class of overlapping Schwarz algorithms has been extensively studied for linear elliptic finite element problems. In this presentation, the author considers the solution of systems of nonlinear algebraic equations arising from the finite element discretization of some nonlinear elliptic equations. Several overlapping Schwarz algorithms, including the additive and multiplicative versions, with inexact Newton acceleration will be discussed. The author shows that the convergence rate of the Newton`s method is independent of the mesh size used in the finite element discretization, and also independent of the number of subdomains into which the original domain in decomposed. Numerical examples will be presented.

Precise magnetostatic field using the finite element method; Calculo de campos magnetostaticos com precisao utilizando o metodo dos elementos finitos

Energy Technology Data Exchange (ETDEWEB)

Nascimento, Francisco Rogerio Teixeira do

2013-07-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)

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 ·    ...

Finite element analysis of multilayer coextrusion.

Energy Technology Data Exchange (ETDEWEB)

Hopkins, Matthew Morgan; Schunk, Peter Randall; Baer, Thomas A. (Proctor & Gamble Company, West Chester, OH); Mrozek, Randy A. (Army Research Laboratory, Adelphi, MD); Lenhart, Joseph Ludlow (Army Research Laboratory, Adelphi, MD); Rao, Rekha Ranjana; Collins, Robert (Oak Ridge National Laboratory); Mondy, Lisa Ann

2011-09-01

Multilayer coextrusion has become a popular commercial process for producing complex polymeric products from soda bottles to reflective coatings. A numerical model of a multilayer coextrusion process is developed based on a finite element discretization and two different free-surface methods, an arbitrary-Lagrangian-Eulerian (ALE) moving mesh implementation and an Eulerian level set method, to understand the moving boundary problem associated with the polymer-polymer interface. The goal of this work is to have a numerical capability suitable for optimizing and troubleshooting the coextrusion process, circumventing flow instabilities such as ribbing and barring, and reducing variability in layer thickness. Though these instabilities can be both viscous and elastic in nature, for this work a generalized Newtonian description of the fluid is used. Models of varying degrees of complexity are investigated including stability analysis and direct three-dimensional finite element free surface approaches. The results of this work show how critical modeling can be to reduce build test cycles, improve material choices, and guide mold design.

Finite element analysis of ARPS structures

International Nuclear Information System (INIS)

Ruhkamp, J.D.; McDougal, J.R.; Kramer, D.P.

1998-01-01

Algor finite element software was used to determine the stresses and deflections in the metallic walls of Advanced Radioisotope Power Systems (ARPS) designs. The preliminary design review of these systems often neglects the structural integrity of the design which can effect fabrication and the end use of the design. Before finite element analysis (FEA) was run on the canister walls of the thermophotovoltaic (TPV) generator, hand calculations were used to approximate the stresses and deflections in a flat plate. These results compared favorably to the FEA results of a similar size flat plate. The AMTEC (Alkali Metal Thermal-to-Electric Conversion) cells were analyzed by FEA and the results compared to two cells that were mechanically tested. The mechanically tested cells buckled in the thin sections, one at the top and one in the lower section. The FEA predicted similar stress and shape results but the critical buckling load was found to be very shape dependent

Performance of mixed formulations for the particle finite element method in soil mechanics problems

Science.gov (United States)

Monforte, Lluís; Carbonell, Josep Maria; Arroyo, Marcos; Gens, Antonio

2017-07-01

This paper presents a computational framework for the numerical analysis of fluid-saturated porous media at large strains. The proposal relies, on one hand, on the particle finite element method (PFEM), known for its capability to tackle large deformations and rapid changing boundaries, and, on the other hand, on constitutive descriptions well established in current geotechnical analyses (Darcy's law; Modified Cam Clay; Houlsby hyperelasticity). An important feature of this kind of problem is that incompressibility may arise either from undrained conditions or as a consequence of material behaviour; incompressibility may lead to volumetric locking of the low-order elements that are typically used in PFEM. In this work, two different three-field mixed formulations for the coupled hydromechanical problem are presented, in which either the effective pressure or the Jacobian are considered as nodal variables, in addition to the solid skeleton displacement and water pressure. Additionally, several mixed formulations are described for the simplified single-phase problem due to its formal similitude to the poromechanical case and its relevance in geotechnics, since it may approximate the saturated soil behaviour under undrained conditions. In order to use equal-order interpolants in displacements and scalar fields, stabilization techniques are used in the mass conservation equation of the biphasic medium and in the rest of scalar equations. Finally, all mixed formulations are assessed in some benchmark problems and their performances are compared. It is found that mixed formulations that have the Jacobian as a nodal variable perform better.

Modal representation of geometrically nonlinear behavior by the finite element method

International Nuclear Information System (INIS)

Nagy, D.A.

1977-01-01

A method is presented for representing mild geometrically nonlinear static behavior of thin-type structures, within the finite element method, in terms of linear elastic and linear (bifurcation) buckling analysis results for structural loading or geometry situations which violate the idealized restrictive (perfect) interpretation of linear behavior up to bifurcation. Formulation of the finite element displacement method for material linearity but retaining the full, nonlinear strain-displacement relations (geometric nonlinearity) leads to highly nonlinear equations relating the unknown nodal generalized displacements r to the applied loading R. Restriction to small strains alone does not linearize these equations for thin-type structural configurations; only explicitly requiring that all products of displacement gadients be much smaller than the gadients themselves reduces the equations to the familiar linear form Ksub(e)r=R, where Ksub(e) is the elastic stiffness. Assuming then that the solutions r of the linear equations also satisfies the full nonlinear equations (i.e., that the above explicit requirement is satisfied), a second solution to the full equations can be sought for a one-parameter loading path lambdaR, leading to the well-known linear (bifurcation) buckling eigenvalue problem Ksub(e)X=-Ksub(g)XΛ where Ksub(g) is the geometric stiffness, X the matrix whose columns are the eigenvectors (so-called buckling mode shapes) and Λ is a diagonal matrix of eigenvalues lambda(i) (so-called load scale factors). From the viewpoint of the practising structural analyst using finite element software, the method presented here gives broader and deeper significance to an existing linear (bifurcation) buckling analysis capability, in that the additional computations are minimal beyond those already required for a linear static and buckling analysis, and should be easily performable within any well-designed general purpose finite element system

A short summary on finite element modelling of fatigue crack closure

Energy Technology Data Exchange (ETDEWEB)

Singh, Konjengbam Darunkumar [Indian Institute of Technology, Guwahati (India); Parry, Matthew Roger [Airbus Operations Ltd, Bristol(United Kingdom); Sinclair, Ian [University of Southampton, Southampton (United Kingdom)

2011-12-15

This paper presents a short summary pertaining to the finite element modelling of fatigue crack closure. Several key issues related to finite element modelling of fatigue crack closure are highlighted: element type, mesh refinement, stabilization of crack closure, crack-tip node release scheme, constitutive model, specimen geometry, stress-states (i.e., plane stress, plane strain), crack closure monitoring. Reviews are presented for both straight and deflected cracks.

Overlapping Schwarz for Nonlinear Problems. An Element Agglomeration Nonlinear Additive Schwarz Preconditioned Newton Method for Unstructured Finite Element Problems

Energy Technology Data Exchange (ETDEWEB)

Cai, X C; Marcinkowski, L; Vassilevski, P S

2005-02-10

This paper extends previous results on nonlinear Schwarz preconditioning ([4]) to unstructured finite element elliptic problems exploiting now nonlocal (but small) subspaces. The non-local finite element subspaces are associated with subdomains obtained from a non-overlapping element partitioning of the original set of elements and are coarse outside the prescribed element subdomain. The coarsening is based on a modification of the agglomeration based AMGe method proposed in [8]. Then, the algebraic construction from [9] of the corresponding non-linear finite element subproblems is applied to generate the subspace based nonlinear preconditioner. The overall nonlinearly preconditioned problem is solved by an inexact Newton method. Numerical illustration is also provided.

Finite element analysis of the Girkmann problem using the modern hp-version and the classical h-version

KAUST Repository

Niemi, Antti

2011-06-03

We perform finite element analysis of the so called Girkmann problem in structural mechanics. The problem involves an axially symmetric spherical shell stiffened with a foot ring and is approached (1) by using the axisymmetric formulation of linear elasticity theory and (2) by using a dimensionally reduced shell-ring model. In the first approach the problem is solved with a fully automatic hp-adaptive finite element solver whereas the classical h-version of the finite element method is used in the second approach. We study the convergence behaviour of the different numerical models and show that accurate stress resultants can be obtained with both models by using effective post-processing formulas. © Springer-Verlag London Limited 2011.

Investigation of Forming Performance of Laminated Steel Sheets Using Finite Element Analyses

International Nuclear Information System (INIS)

Liu Wenning; Sun Xin; Ruokolainen, Robert; Gayden Xiaohong

2007-01-01

Laminated steel sheets have been used in automotive structures for reducing in-cabin noise. However, due to the marked difference in material properties of the different laminated layers, integrating laminated steel parts into the manufacturing processes can be challenging. Especially, the behavior of laminated sheets during forming processes is very different from that of monolithic steel sheets. During the deep-draw forming process, large shear deformation and corresponding high interfacial stress may initiate and propagate interfacial cracks between the core polymer and the metal skin, hence degrading the performance of the laminated sheets. In this paper, the formability of the laminated steel sheets is investigated by means of numerical analysis. The goal of this work is to gain insight into the relationship between the individual properties of the laminated sheet layers and the corresponding formability of the laminated sheet as a whole, eventually leading to reliable design and successful forming process development of such materials. Finite element analyses of laminate sheet forming are presented. Effects of polymer core thickness and viscoelastic properties of the polymer core, as well as punching velocity, are also investigated

Generalized multiscale finite element methods: Oversampling strategies

KAUST Repository

Efendiev, Yalchin R.; Galvis, Juan; Li, Guanglian; Presho, Michael

2014-01-01

In this paper, we propose oversampling strategies in the generalized multiscale finite element method (GMsFEM) framework. The GMsFEM, which has been recently introduced in Efendiev et al. (2013b) [Generalized Multiscale Finite Element Methods, J. Comput. Phys., vol. 251, pp. 116-135, 2013], allows solving multiscale parameter-dependent problems at a reduced computational cost by constructing a reduced-order representation of the solution on a coarse grid. The main idea of the method consists of (1) the construction of snapshot space, (2) the construction of the offline space, and (3) construction of the online space (the latter for parameter-dependent problems). In Efendiev et al. (2013b) [Generalized Multiscale Finite Element Methods, J. Comput. Phys., vol. 251, pp. 116-135, 2013], it was shown that the GMsFEM provides a flexible tool to solve multiscale problems with a complex input space by generating appropriate snapshot, offline, and online spaces. In this paper, we develop oversampling techniques to be used in this context (see Hou and Wu (1997) where oversampling is introduced for multiscale finite element methods). It is known (see Hou and Wu (1997)) that the oversampling can improve the accuracy of multiscale methods. In particular, the oversampling technique uses larger regions (larger than the target coarse block) in constructing local basis functions. Our motivation stems from the analysis presented in this paper, which shows that when using oversampling techniques in the construction of the snapshot space and offline space, GMsFEM will converge independent of small scales and high contrast under certain assumptions. We consider the use of a multiple eigenvalue problems to improve the convergence and discuss their relation to single spectral problems that use oversampled regions. The oversampling procedures proposed in this paper differ from those in Hou and Wu (1997). In particular, the oversampling domains are partially used in constructing local

PROBABILISTIC FINITE ELEMENT ANALYSIS OF A HEAVY DUTY RADIATOR UNDER INTERNAL PRESSURE LOADING

Directory of Open Access Journals (Sweden)

ROBIN ROY P.

2017-09-01

Full Text Available Engine cooling is vital in keeping the engine at most efficient temperature for the different vehicle speed and operating road conditions. Radiator is one of the key components in the heavy duty engine cooling system. Heavy duty radiator is subjected to various kinds of loading such as pressure, thermal, vibration, internal erosion, external corrosion, creep. Pressure cycle durability is one of the most important characteristic in the design of heavy duty radiator. Current design methodologies involve design of heavy duty radiator using the nominal finite element approach which does not take into account of the variations occurring in the geometry, material and boundary condition, leading to over conservative and uneconomical designs of radiator system. A new approach is presented in the paper to integrate traditional linear finite element method and probabilistic approach to design a heavy duty radiator by including the uncertainty in the computational model. As a first step, nominal run is performed with input design variables and desired responses are extracted. A probabilistic finite elementanalysis is performed to identify the robust designs and validated for reliability. Probabilistic finite element includes the uncertainty of the material thickness, dimensional and geometrical variation. Gaussian distribution is employed to define the random variation and uncertainty. Monte Carlo method is used to generate the random design points.Output response distributions of the random design points are post-processed using different statistical and probability technique to find the robust design. The above approach of systematic virtual modelling and analysis of the data helps to find efficient and reliable robust design.

Groundwater flow analysis using mixed hybrid finite element method for radioactive waste disposal facilities

International Nuclear Information System (INIS)

Aoki, Hiroomi; Shimomura, Masanori; Kawakami, Hiroto; Suzuki, Shunichi

2011-01-01

In safety assessments of radioactive waste disposal facilities, ground water flow analysis are used for calculating the radionuclide transport pathway and the infiltration flow rate of groundwater into the disposal facilities. For this type of calculations, the mixed hybrid finite element method has been used and discussed about the accuracy of ones in Europe. This paper puts great emphasis on the infiltration flow rate of groundwater into the disposal facilities, and describes the accuracy of results obtained from mixed hybrid finite element method by comparing of local water mass conservation and the reliability of the element breakdown numbers among the mixed hybrid finite element method, finite volume method and nondegenerated finite element method. (author)

Finite element elastic-plastic analysis of LMFBR components

International Nuclear Information System (INIS)

Levy, A.; Pifko, A.; Armen, H. Jr.

1978-01-01

The present effort involves the development of computationally efficient finite element methods for accurately predicting the isothermal elastic-plastic three-dimensional response of thick and thin shell structures subjected to mechanical and thermal loads. This work will be used as the basis for further development of analytical tools to be used to verify the structural integrity of liquid metal fast breeder reactor (LMFBR) components. The methods presented here have been implemented into the three-dimensional solid element module (HEX) of the Grumman PLANS finite element program. These methods include the use of optimal stress points as well as a variable number of stress points within an element. This allows monitoring the stress history at many points within an element and hence provides an accurate representation of the elastic-plastic boundary using a minimum number of degrees of freedom. Also included is an improved thermal stress analysis capability in which the temperature variation and corresponding thermal strain variation are represented by the same functional form as the displacement variation. Various problems are used to demonstrate these improved capabilities. (Auth.)

Finite Element Residual Stress Analysis of Planetary Gear Tooth

Directory of Open Access Journals (Sweden)

Jungang Wang

2013-01-01

Full Text Available A method to simulate residual stress field of planetary gear is proposed. In this method, the finite element model of planetary gear is established and divided to tooth zone and profile zone, whose different temperature field is set. The gear's residual stress simulation is realized by the thermal compression stress generated by the temperature difference. Based on the simulation, the finite element model of planetary gear train is established, the dynamic meshing process is simulated, and influence of residual stress on equivalent stress of addendum, pitch circle, and dedendum of internal and external meshing planetary gear tooth profile is analyzed, according to non-linear contact theory, thermodynamic theory, and finite element theory. The results show that the equivalent stresses of planetary gear at both meshing and nonmeshing surface are significantly and differently reduced by residual stress. The study benefits fatigue cracking analysis and dynamic optimization design of planetary gear train.

2D Finite Element Model of a CIGS Module

Energy Technology Data Exchange (ETDEWEB)

Janssen, G.J.M.; Slooff, L.H.; Bende, E.E. [ECN Solar Energy, P.O.Box 1, NL-1755 ZG Petten (Netherlands)

2012-06-15

The performance of thin-film CIGS (Copper indium gallium selenide) modules is often limited due to inhomogeneities in CIGS layers. A 2-dimensional Finite Element Model for CIGS modules is presented that predicts the impact of such inhomogeneities on the module performance. Results are presented of a module with a region of poor diode characteristics. It is concluded that according to this model the effects of poor diodes depend strongly on their location in the module and on their dispersion over the module surface. Due to its generic character the model can also be applied to other series connections of photovoltaic cells.

2D - Finite element model of a CIGS module

Energy Technology Data Exchange (ETDEWEB)

Janssen, G.J.M.; Slooff, L.H.; Bende, E.E. [ECN Solar Energy, Petten (Netherlands)

2012-09-15

The performance of thin-film CIGS modules is often limited due to inhomogeneities in CIGS layers. A 2-dimensional Finite Element Model for CIGS modules is demonstrated that predicts the impact of such inhomogeneities on the module performance. Results are presented of a module with a region of poor diode characteristics. It is concluded that according to this model the effects of poor diodes depend strongly on their location in the module and on their dispersion over the module surface. Due to its generic character the model can also be applied to other series connections of photovoltaic cells.

Finite element modeling of the filament winding process using ABAQUS

OpenAIRE

Miltenberger, Louis C.

1992-01-01

A comprehensive stress model of the filament winding fabrication process, previously implemented in the finite element program, WACSAFE, was implemented using the ABAQUS finite element software package. This new implementation, referred to as the ABWACSAFE procedure, consists of the ABAQUS software and a pre/postprocessing routine that was developed to prepare necessary ABAQUS input files and process ABAQUS displacement results for stress and strain computation. The ABWACSAF...

A parallel finite element procedure for contact-impact problems using edge-based smooth triangular element and GPU

Science.gov (United States)

Cai, Yong; Cui, Xiangyang; Li, Guangyao; Liu, Wenyang

2018-04-01

The edge-smooth finite element method (ES-FEM) can improve the computational accuracy of triangular shell elements and the mesh partition efficiency of complex models. In this paper, an approach is developed to perform explicit finite element simulations of contact-impact problems with a graphical processing unit (GPU) using a special edge-smooth triangular shell element based on ES-FEM. Of critical importance for this problem is achieving finer-grained parallelism to enable efficient data loading and to minimize communication between the device and host. Four kinds of parallel strategies are then developed to efficiently solve these ES-FEM based shell element formulas, and various optimization methods are adopted to ensure aligned memory access. Special focus is dedicated to developing an approach for the parallel construction of edge systems. A parallel hierarchy-territory contact-searching algorithm (HITA) and a parallel penalty function calculation method are embedded in this parallel explicit algorithm. Finally, the program flow is well designed, and a GPU-based simulation system is developed, using Nvidia's CUDA. Several numerical examples are presented to illustrate the high quality of the results obtained with the proposed methods. In addition, the GPU-based parallel computation is shown to significantly reduce the computing time.

Finite element reliability analysis of fatigue life

International Nuclear Information System (INIS)

Harkness, H.H.; Belytschko, T.; Liu, W.K.

1992-01-01

Fatigue reliability is addressed by the first-order reliability method combined with a finite element method. Two-dimensional finite element models of components with cracks in mode I are considered with crack growth treated by the Paris law. Probability density functions of the variables affecting fatigue are proposed to reflect a setting where nondestructive evaluation is used, and the Rosenblatt transformation is employed to treat non-Gaussian random variables. Comparisons of the first-order reliability results and Monte Carlo simulations suggest that the accuracy of the first-order reliability method is quite good in this setting. Results show that the upper portion of the initial crack length probability density function is crucial to reliability, which suggests that if nondestructive evaluation is used, the probability of detection curve plays a key role in reliability. (orig.)

Finite element electromagnetic field computation on the Sequent Symmetry 81 parallel computer

International Nuclear Information System (INIS)

Ratnajeevan, S.; Hoole, H.

1990-01-01

Finite element field analysis algorithms lend themselves to parallelization and this fact is exploited in this paper to implement a finite element analysis program for electromagnetic field computation on the Sequent Symmetry 81 parallel computer with three processors. In terms of waiting time, the maximum gains are to be made in matrix solution and therefore this paper concentrates on the gains in parallelizing the solution part of finite element analysis. An outline of how parallelization could be exploited in most finite element operations is given in this paper although the actual implemention of parallelism on the Sequent Symmetry 81 parallel computer was in sparsity computation, matrix assembly and the matrix solution areas. In all cases, the algorithms were modified suit the parallel programming application rather than allowing the compiler to parallelize on existing algorithms

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

Advanced finite element method in structural engineering

CERN Document Server

Long, Yu-Qiu; Long, Zhi-Fei

2009-01-01

This book systematically introduces the research work on the Finite Element Method completed over the past 25 years. Original theoretical achievements and their applications in the fields of structural engineering and computational mechanics are discussed.

Investigation of the Behavior of Steel Shear Walls Using Finite Elements Analysis

Directory of Open Access Journals (Sweden)

K. Abubakri

2016-10-01

Full Text Available Currently, steel shear walls are considered by engineers as an economic method against lateral loads imposed by wind and earthquake in tall structures. Accordingly, there is a growing need to develop accurate methods alongside approximation methods to estimate the behavior of these structural elements. The finite element technique is one of the strongest numerical methods in analysis of solid mechanics problems. Finite element analysis however requires high technical knowledge of the behavioral models of materials. Therefore, it is less used by designers for certain structural elements such as steel shear walls. This study examines the failure mechanism of steel shear walls using finite elements analysis and validates this modeling by comparing the results with experimental studies.

Application of finite element techniques in predicting the acoustic properties of turbofan inlets

Science.gov (United States)

Majjigi, R. K.; Sigman, R. K.; Zinn, B. T.

1978-01-01

An analytical technique was developed for predicting the acoustic performance of turbofan inlets carrying a subsonic axisymmetric steady flow. The finite element method combined with the method of weighted residuals is used in predicting the acoustic properties of variable area, annular ducts with or without acoustic treatments along their walls. An approximate solution for the steady inviscid flow field is obtained using an integral method for calculating the incompressible potential flow field in the inlet with a correction to account for compressibility effects. The accuracy of the finite element technique was assessed by comparison with available analytical solutions for the problems of plane and spinning wave propagation through a hard walled annular cylinder with a constant mean flow.

Ring tests on high density polyethylene: Full investigation assisted by finite element modeling

International Nuclear Information System (INIS)

Laiarinandrasana, L.; Devilliers, C.; Oberti, S.; Gaudichet, E.; Fayolle, B.; Lucatelli, J.M.

2011-01-01

In order to characterize the mechanical behavior of HDPE pipes, the ASTM D 2290-04 standard recommends carrying out tensile tests on notched rings, cut out from the pipe. This very simple test is also utilized to investigate the aging effect of the pipe by determining the strain at failure. Comparison between full ring and notched ring mechanical responses are discussed. Constitutive modeling including strain rate effects was performed by finite element analysis. This allowed a better understanding of the stress state in the cross section perpendicular to the loading direction. Additionally, the influence of a thin layer of oxidized HDPE in the inner wall of the ring was studied in the light of the finite element results.

Matlab and C programming for Trefftz finite element methods

CERN Document Server

Qin, Qing-Hua

2008-01-01

Although the Trefftz finite element method (FEM) has become a powerful computational tool in the analysis of plane elasticity, thin and thick plate bending, Poisson's equation, heat conduction, and piezoelectric materials, there are few books that offer a comprehensive computer programming treatment of the subject. Collecting results scattered in the literature, MATLAB® and C Programming for Trefftz Finite Element Methods provides the detailed MATLAB® and C programming processes in applications of the Trefftz FEM to potential and elastic problems. The book begins with an introduction to th

Preconditioning for Mixed Finite Element Formulations of Elliptic Problems

KAUST Repository

Wildey, Tim; Xue, Guangri

2013-01-01

In this paper, we discuss a preconditioning technique for mixed finite element discretizations of elliptic equations. The technique is based on a block-diagonal approximation of the mass matrix which maintains the sparsity and positive definiteness of the corresponding Schur complement. This preconditioner arises from the multipoint flux mixed finite element method and is robust with respect to mesh size and is better conditioned for full permeability tensors than a preconditioner based on a diagonal approximation of the mass matrix. © Springer-Verlag Berlin Heidelberg 2013.

A Kohn–Sham equation solver based on hexahedral finite elements

International Nuclear Information System (INIS)

Fang Jun; Gao Xingyu; Zhou Aihui

2012-01-01

We design a Kohn–Sham equation solver based on hexahedral finite element discretizations. The solver integrates three schemes proposed in this paper. The first scheme arranges one a priori locally-refined hexahedral mesh with appropriate multiresolution. The second one is a modified mass-lumping procedure which accelerates the diagonalization in the self-consistent field iteration. The third one is a finite element recovery method which enhances the eigenpair approximations with small extra work. We carry out numerical tests on each scheme to investigate the validity and efficiency, and then apply them to calculate the ground state total energies of nanosystems C 60 , C 120 , and C 275 H 172 . It is shown that our solver appears to be computationally attractive for finite element applications in electronic structure study.

Investigation of faulted tunnel models by combined photoelasticity and finite element analysis

International Nuclear Information System (INIS)

Ladkany, S.G.; Huang, Yuping

1994-01-01

Models of square and circular tunnels with short faults cutting through their surfaces are investigated by photoelasticity. These models, when duplicated by finite element analysis can predict the stress states of square or circular faulted tunnels adequately. Finite element analysis, using gap elements, may be used to investigate full size faulted tunnel system

Simple one-dimensional finite element algorithm with multi-dimensional capabilities

International Nuclear Information System (INIS)

Pepper, D.W.; Baker, A.J.

1978-01-01

The application of the finite element procedure for the solution of partial differential equations is gaining widespread acceptance. The ability of the finite element procedure to solve problems which are arbitrarily shaped as well as the alleviation of boundary condition problems is well known. By using local interpolation functionals over each subdomain, or element, a set of linearized algebraic equations are obtained which can be solved using any direct, iterative, or inverse numerical technique. Subsequent use of an explicit or implicit integration procedure permits closure of the solution over the global domain

Microstructural analysis of alumina chromium composites by X-ray tomography and 3-D finite element simulation of thermal stresses

International Nuclear Information System (INIS)

Geandier, G.; Hazotte, A.; Denis, S.; Mocellin, A.; Maire, E.

2003-01-01

X-ray microtomography is used to measure volume fraction and connectivity of the metallic phase in an alumina-chromium composite. Reconstructed images are used as input data for a finite element calculation of the residual thermal stresses. Results confirm the main trends shown by similar calculations previously performed on less-realistic finite element models

Microstructural analysis of alumina chromium composites by X-ray tomography and 3-D finite element simulation of thermal stresses

Energy Technology Data Exchange (ETDEWEB)

Geandier, G.; Hazotte, A.; Denis, S.; Mocellin, A.; Maire, E

2003-04-14

X-ray microtomography is used to measure volume fraction and connectivity of the metallic phase in an alumina-chromium composite. Reconstructed images are used as input data for a finite element calculation of the residual thermal stresses. Results confirm the main trends shown by similar calculations previously performed on less-realistic finite element models.

Review on Finite Element Method * ERHUNMWUN, ID ...

African Journals Online (AJOL)

ADOWIE PERE

ABSTRACT: In this work, we have discussed what Finite Element Method (FEM) is, its historical development, advantages and ... residual procedures, are examples of the direct approach ... The paper centred on the "stiffness and deflection of ...

Finite element evaluation of erosion/corrosion affected reducing elbow

International Nuclear Information System (INIS)

Basavaraju, C.

1996-01-01

Erosion/corrosion is a primary source for wall thinning or degradation of carbon steel piping systems in service. A number of piping failures in the power industry have been attributed to erosion/corrosion. Piping elbow is one of such susceptible components for erosion/corrosion because of increased flow turbulence due to its geometry. In this paper, the acceptability of a 12 in. x 8 in. reducing elbow in RHR service water pump discharge piping, which experienced significant degradation due to wall thinning in localized areas, was evaluated using finite element analysis methodology. Since the simplified methods showed very small margin and recommended replacement of the elbow, a detailed 3-D finite element model was built using shell elements and analyzed for internal pressure and moment loadings. The finite element analysis incorporated the U.T. measured wall thickness data at various spots that experienced wall thinning. The results showed that the elbow is acceptable as-is until the next fuel cycle. FEA, though cumbersome, and time consuming is a valuable analytical tool in making critical decisions with regard to component replacement of border line situation cases, eliminating some conservatism while not compromising the safety

Micro-finite-element method to assess elastic properties of trabecular bone at micro- and macroscopic level.

Science.gov (United States)

Rieger, R; Auregan, J C; Hoc, T

2018-03-01

The objective of the present study is to assess the mechanical behavior of trabecular bone based on microCT imaging and micro-finite-element analysis. In this way two methods are detailed: (i) direct determination of macroscopic elastic property of trabecular bone; (ii) inverse approach to assess mechanical properties of trabecular bone tissue. Thirty-five females and seven males (forty-two subjects) mean aged (±SD) 80±11.7 years from hospitals of Assistance publique-Hôpitaux de Paris (AP-HP) diagnosed with osteoporosis following a femoral neck fracture due to a fall from standing were included in this study. Fractured heads were collected during hip replacement surgery. Standardized bone cores were removed from the femoral head's equator by a trephine in a water bath. MicroCT images acquisition and analysis were performed with CTan ® software and bone volume fraction was then determined. Micro-finite-element simulations were per-formed using Abaqus 6.9-2 ® software in order to determine the macroscopic mechanical behaviour of the trabecular bone. After microCT acquisition, a longitudinal compression test was performed and the experimental macroscopic Young's Modulus was extracted. An inverse approach based on the whole trabecular bone's mechanical response and micro-finite-element analysis was performed to determine microscopic mechanical properties of trabecular bone. In the present study, elasticity of the tissue was shown to be similar to that of healthy tissue but with a lower yield stress. Classical histomorphometric analysis form microCT imaging associated with an inverse micro-finite-element method allowed to assess microscopic mechanical trabecular bone parameters. Copyright © 2017 Elsevier Masson SAS. All rights reserved.

Application of the finite element method to neutronics problems with inhomogeneous boundray conditions

International Nuclear Information System (INIS)

Yoo, K.J.

1982-01-01

The albedo boundary conditions are incorporated into the finite element method using bicubic Hermite element functions in order to reduce the computer memory and computation time in two-group diffusion calculations by excluding the reflector regions in computation space. The basis functions at the core-reflector interfaces are newly established to satisfy the albedo boundary conditions, and then the ''weak'' form of two-group diffusion equations is discretized using the principle of the weighted residual method in combination with the Galerkin approximation. The discretized two-group diffusion equation is then solved by the Gaussian elimination method with the scaled column pivoting algorithm in one-dimensional problem and Gauss-Seidel method in two-dimensional problem. Prior to the application of the method to two-group diffusion problems, the same method is applied to the one-speed neutron transport equation in a bare slab reactor with the vacuum boundary condition to confirm its usefulness in the diffusion calculations. To investigate the applicability of our diffusion method, several numerical calculations are performed: two-dimensional IAEA benchmark problem and two-dimensional ZION problem. The results are compared with the available results from the conventional finite difference and other finite element methods. If the albedo values are appropriately adjusted, our results of the two-dimensional IAEA benchmark problem are agreed within 0.002% of ksub(eff) with the fine mesh PDQ results. Comparing with CITATION results, one-eighth of core memory and one-fifteenth of computing time are required to obtain the same accuracy even though no acceleration technique is used in the present case. Also, it is found that the results are comparable with the other finite element results. However, no significant saving is obtained in computation time comparing with the other finite element results, where the reflector regions are explicity included. This mainly comes from

Modeling and Finite Element Analysis of Load-Carrying Performance of a Wind Turbine Considering the Influence of Assembly Factors

Directory of Open Access Journals (Sweden)

Jianmei Wang

2017-03-01

Full Text Available In this work, a wind turbine shrink disk is used as the research object to investigate load-carrying performance of a multi-layer interference fit, and the theoretical model and finite element model are constructed. According to those models, a MW-level turbine shrink disk is designed, and a test device is developed to apply torque to this turbine shrink disk by hydraulic jack. Then, the circumferential slip between the contact surfaces is monitored and the slip of all contact surfaces is zero. This conclusion verifies the reasonability of the proposed models. The effect of the key influencing factors, such as machining deviation, assembly clearance and propel stroke, were analyzed. The contact pressure and load torque of the mating surfaces were obtained by building typical models with different parameters using finite element analysis (FEA. The results show that the minimum assembly clearance and the machining deviation within the machining range have little influence on load-carrying performance of multi-layer interference fit, while having a greater influence on the maximum assembly clearance and the propel stroke. The results also show that the load-carrying performance of a multiple-layer interference fit can be ensured only if the key factors are set within a reasonable design range. To avoid the abnormal operation of equipment caused by insufficient load torque, the propel stroke during practical assembly should be at least 0.95 times the designed propel stroke, which is significant in guiding the design and assembly of the multi-layer interference fit.

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.

Behaviour of Lagrangian triangular mixed fluid finite elements

Indian Academy of Sciences (India)

The behaviour of mixed fluid finite elements, formulated based on the Lagrangian frame of reference, is investigated to understand the effects of locking due to incompressibility and irrotational constraints. For this purpose, both linear and quadratic mixed triangular fluid elements are formulated. It is found that there exists a ...

A fully coupled finite element framework for thermal fracturing simulation in subsurface cold CO2 injection

Directory of Open Access Journals (Sweden)

Shunde Yin

2018-03-01

Simulation of thermal fracturing during cold CO2 injection involves the coupled processes of heat transfer, mass transport, rock deforming as well as fracture propagation. To model such a complex coupled system, a fully coupled finite element framework for thermal fracturing simulation is presented. This framework is based on the theory of non-isothermal multiphase flow in fracturing porous media. It takes advantage of recent advances in stabilized finite element and extended finite element methods. The stabilized finite element method overcomes the numerical instability encountered when the traditional finite element method is used to solve the convection dominated heat transfer equation, while the extended finite element method overcomes the limitation with traditional finite element method that a model has to be remeshed when a fracture is initiated or propagating and fracturing paths have to be aligned with element boundaries.

Finite element modelling of solidification phenomena

Indian Academy of Sciences (India)

Unknown

Abstract. The process of solidification process is complex in nature and the simulation of such process is required in industry before it is actually undertaken. Finite element method is used to simulate the heat transfer process accompanying the solidification process. The metal and the mould along with the air gap formation ...

Improving torque per kilogram magnet of permanent magnet couplings using finite element analysis

DEFF Research Database (Denmark)

Högberg, Stig; Jensen, Bogi Bech; Bendixen, Flemming Buus

2013-01-01

This paper presents the methodology and subsequent findings of a performance-improvement routine that employs automated finite element (FE) analysis to increase the torque-per-kilogram-magnet (TPKM) of a permanent magnet coupling (PMC). The routine is applied to a commercially available cylindrical...

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 modelling of fibre Bragg grating strain sensors and experimental validation

Science.gov (United States)

Malik, Shoaib A.; Mahendran, Ramani S.; Harris, Dee; Paget, Mark; Pandita, Surya D.; Machavaram, Venkata R.; Collins, David; Burns, Jonathan M.; Wang, Liwei; Fernando, Gerard F.

2009-03-01

Fibre Bragg grating (FBG) sensors continue to be used extensively for monitoring strain and temperature in and on engineering materials and structures. Previous researchers have also developed analytical models to predict the loadtransfer characteristics of FBG sensors as a function of applied strain. The general properties of the coating or adhesive that is used to surface-bond the FBG sensor to the substrate has also been modelled using finite element analysis. In this current paper, a technique was developed to surface-mount FBG sensors with a known volume and thickness of adhesive. The substrates used were aluminium dog-bone tensile test specimens. The FBG sensors were tensile tested in a series of ramp-hold sequences until failure. The reflected FBG spectra were recorded using a commercial instrument. Finite element analysis was performed to model the response of the surface-mounted FBG sensors. In the first instance, the effect of the mechanical properties of the adhesive and substrate were modelled. This was followed by modelling the volume of adhesive used to bond the FBG sensor to the substrate. Finally, the predicted values obtained via finite element modelling were correlated to the experimental results. In addition to the FBG sensors, the tensile test specimens were instrumented with surface-mounted electrical resistance strain gauges.

Energy flow in plate assembles by hierarchical version of finite element method

DEFF Research Database (Denmark)

Wachulec, Marcin; Kirkegaard, Poul Henning

method has been proposed. In this paper a modified hierarchical version of finite element method is used for modelling of energy flow in plate assembles. The formulation includes description of in-plane forces so that planes lying in different planes can be modelled. Two examples considered are: L......The dynamic analysis of structures in medium and high frequencies are usually focused on frequency and spatial averages of energy of components, and not the displacement/velocity fields. This is especially true for structure-borne noise modelling. For the analysis of complicated structures...... the finite element method has been used to study the energy flow. The finite element method proved its usefulness despite the computational expense. Therefore studies have been conducted in order to simplify and reduce the computations required. Among others, the use of hierarchical version of finite element...

Finite element bending behaviour of discretely delaminated ...

African Journals Online (AJOL)

user

due to their light weight, high specific strength and stiffness properties. ... cylindrical shell roofs respectively using finite element method with centrally located .... where { }ε and { }γ are the direct and shear strains in midplane and { }κ denotes ...

Multisymplectic Structure-Preserving in Simple Finite Element Method in High Dimensional Case

Institute of Scientific and Technical Information of China (English)

BAI Yong-Qiang; LIU Zhen; PEI Ming; ZHENG Zhu-Jun

2003-01-01

In this paper, we study a finite element scheme of some semi-linear elliptic boundary value problems inhigh-dimensional space. With uniform mesh, we find that, the numerical scheme derived from finite element method cankeep a preserved multisymplectic structure.

A Dual Super-Element Domain Decomposition Approach for Parallel Nonlinear Finite Element Analysis

Science.gov (United States)

Jokhio, G. A.; Izzuddin, B. A.

2015-05-01

This article presents a new domain decomposition method for nonlinear finite element analysis introducing the concept of dual partition super-elements. The method extends ideas from the displacement frame method and is ideally suited for parallel nonlinear static/dynamic analysis of structural systems. In the new method, domain decomposition is realized by replacing one or more subdomains in a "parent system," each with a placeholder super-element, where the subdomains are processed separately as "child partitions," each wrapped by a dual super-element along the partition boundary. The analysis of the overall system, including the satisfaction of equilibrium and compatibility at all partition boundaries, is realized through direct communication between all pairs of placeholder and dual super-elements. The proposed method has particular advantages for matrix solution methods based on the frontal scheme, and can be readily implemented for existing finite element analysis programs to achieve parallelization on distributed memory systems with minimal intervention, thus overcoming memory bottlenecks typically faced in the analysis of large-scale problems. Several examples are presented in this article which demonstrate the computational benefits of the proposed parallel domain decomposition approach and its applicability to the nonlinear structural analysis of realistic structural systems.

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

Numerical solution of recirculating flow by a simple finite element recursion relation

Energy Technology Data Exchange (ETDEWEB)

Pepper, D W; Cooper, R E

1980-01-01

A time-split finite element recursion relation, based on linear basis functions, is used to solve the two-dimensional equations of motion. Recirculating flow in a rectangular cavity and free convective flow in an enclosed container are analyzed. The relation has the advantage of finite element accuracy and finite difference speed and simplicity. Incorporating dissipation parameters in the functionals decreases numerical dispersion and improves phase lag.

Optimization of tissue physical parameters for accurate temperature estimation from finite-element simulation of radiofrequency ablation

International Nuclear Information System (INIS)

Subramanian, Swetha; Mast, T Douglas

2015-01-01

Computational finite element models are commonly used for the simulation of radiofrequency ablation (RFA) treatments. However, the accuracy of these simulations is limited by the lack of precise knowledge of tissue parameters. In this technical note, an inverse solver based on the unscented Kalman filter (UKF) is proposed to optimize values for specific heat, thermal conductivity, and electrical conductivity resulting in accurately simulated temperature elevations. A total of 15 RFA treatments were performed on ex vivo bovine liver tissue. For each RFA treatment, 15 finite-element simulations were performed using a set of deterministically chosen tissue parameters to estimate the mean and variance of the resulting tissue ablation. The UKF was implemented as an inverse solver to recover the specific heat, thermal conductivity, and electrical conductivity corresponding to the measured area of the ablated tissue region, as determined from gross tissue histology. These tissue parameters were then employed in the finite element model to simulate the position- and time-dependent tissue temperature. Results show good agreement between simulated and measured temperature. (note)

Optimization of tissue physical parameters for accurate temperature estimation from finite-element simulation of radiofrequency ablation.

Science.gov (United States)

Subramanian, Swetha; Mast, T Douglas

2015-10-07

Computational finite element models are commonly used for the simulation of radiofrequency ablation (RFA) treatments. However, the accuracy of these simulations is limited by the lack of precise knowledge of tissue parameters. In this technical note, an inverse solver based on the unscented Kalman filter (UKF) is proposed to optimize values for specific heat, thermal conductivity, and electrical conductivity resulting in accurately simulated temperature elevations. A total of 15 RFA treatments were performed on ex vivo bovine liver tissue. For each RFA treatment, 15 finite-element simulations were performed using a set of deterministically chosen tissue parameters to estimate the mean and variance of the resulting tissue ablation. The UKF was implemented as an inverse solver to recover the specific heat, thermal conductivity, and electrical conductivity corresponding to the measured area of the ablated tissue region, as determined from gross tissue histology. These tissue parameters were then employed in the finite element model to simulate the position- and time-dependent tissue temperature. Results show good agreement between simulated and measured temperature.

Fluid-structure finite-element vibrational analysis

Science.gov (United States)

Feng, G. C.; Kiefling, L.

1974-01-01

A fluid finite element has been developed for a quasi-compressible fluid. Both kinetic and potential energy are expressed as functions of nodal displacements. Thus, the formulation is similar to that used for structural elements, with the only differences being that the fluid can possess gravitational potential, and the constitutive equations for fluid contain no shear coefficients. Using this approach, structural and fluid elements can be used interchangeably in existing efficient sparse-matrix structural computer programs such as SPAR. The theoretical development of the element formulations and the relationships of the local and global coordinates are shown. Solutions of fluid slosh, liquid compressibility, and coupled fluid-shell oscillation problems which were completed using a temporary digital computer program are shown. The frequency correlation of the solutions with classical theory is excellent.

Finite element analysis of degraded concrete structures - Workshop proceedings

International Nuclear Information System (INIS)

1999-09-01

This workshop is related to the finite element analysis of degraded concrete structures. It is composed of three sessions. The first session (which title is: the use of finite element analysis in safety assessments) comprises six papers which titles are: Historical Development of Concrete Finite Element Modeling for Safety Evaluation of Accident-Challenged and Aging Concrete Structures; Experience with Finite Element Methods for Safety Assessments in Switzerland; Stress State Analysis of the Ignalina NPP Confinement System; Prestressed Containment: Behaviour when Concrete Cracking is Modelled; Application of FEA for Design and Support of NPP Containment in Russia; Verification Problems of Nuclear Installations Safety Software of Strength Analysis (NISS SA). The second session (title: concrete containment structures under accident loads) comprises seven papers which titles are: Two Application Examples of Concrete Containment Structures under Accident Load Conditions Using Finite Element Analysis; What Kind of Prediction for Leak rates for Nuclear Power Plant Containments in Accidental Conditions; Influence of Different Hypotheses Used in Numerical Models for Concrete At Elevated Temperatures on the Predicted Behaviour of NPP Core Catchers Under Severe Accident Conditions; Observations on the Constitutive Modeling of Concrete Under Multi-Axial States at Elevated Temperatures; Analyses of a Reinforced Concrete Containment with Liner Corrosion Damage; Program of Containment Concrete Control During Operation for the Temelin Nuclear Power Plant; Static Limit Load of a Deteriorated Hyperbolic Cooling Tower. The third session (concrete structures under extreme environmental load) comprised five papers which titles are: Shear Transfer Mechanism of RC Plates After Cracking; Seismic Back Calculation of an Auxiliary Building of the Nuclear Power Plant Muehleberg, Switzerland; Seismic Behaviour of Slightly Reinforced Shear Wall Structures; FE Analysis of Degraded Concrete

Finite element analysis of adanced composite structures containing mechanically fastened joints

International Nuclear Information System (INIS)

Baumann, E.

1982-01-01

Although the usual engineering practice is to ignore joint effects in finite element models of overall structures, there are times when the inclusion of fastener effects in a model is necessary for accurate analysis. This paper describes some simple but accurate methods for accommodating this modeling requirement. The approach involves correlation of test results from a few composite mechanically fastened joints with finite element analyses of joints. It is assumed that if the fastener actions in the test articles can be properly predicted by simple finite element techniques, then these same techniques can be applied to large overall structure models. During the course of this test-analysis effort it was determined that it is possible to obtain correct results for overall structure-joint analyses by using simple modeling concepts provided special care is employed. Also, some emphasis is given in this paper to the importance of properly reducing test data in order to obtain meaningful correlations with finite element analysis. Finally, for those interested, the appendix contains brief descriptions of the test results and failure modes explored in the test program. (orig.)

Multisymplectic Structure－Preserving in Simple Finite Element Method in High Dimensional Case

Institute of Scientific and Technical Information of China (English)

BAIYong-Qiang; LIUZhen; PEIMing; ZHENGZhu-Jun

2003-01-01

In this paper, we study a finite element scheme of some semi-linear elliptic boundary value problems in high-dhnensjonal space. With uniform mesh, we find that, the numerical scheme derived from finite element method can keep a preserved multisymplectic structure.

A Finite Element Analysis of Optimal Variable Thickness Sheets

DEFF Research Database (Denmark)

Petersson, Joakim S

1996-01-01

A quasimixed Finite Element (FE) method for maximum stiffness of variablethickness sheets is analysed. The displacement is approximated with ninenode Lagrange quadrilateral elements and the thickness is approximated aselementwise constant. One is guaranteed that the FE displacement solutionswill ...

Supercomputer implementation of finite element algorithms for high speed compressible flows. Progress report, period ending 30 June 1986

International Nuclear Information System (INIS)

Thornton, E.A.; Ramakrishnan, R.

1986-06-01

Prediction of compressible flow phenomena using the finite element method is of recent origin and considerable interest. Two shock capturing finite element formulations for high speed compressible flows are described. A Taylor-Galerkin formulation uses a Taylor series expansion in time coupled with a Galerkin weighted residual statement. The Taylor-Galerkin algorithms use explicit artificial dissipation, and the performance of three dissipation models are compared. A Petrov-Galerkin algorithm has as its basis the concepts of streamline upwinding. Vectorization strategies are developed to implement the finite element formulations on the NASA Langley VPS-32. The vectorization scheme results in finite element programs that use vectors of length of the order of the number of nodes or elements. The use of the vectorization procedure speeds up processing rates by over two orders of magnitude. The Taylor-Galerkin and Petrov-Galerkin algorithms are evaluated for 2D inviscid flows on criteria such as solution accuracy, shock resolution, computational speed and storage requirements. The convergence rates for both algorithms are enhanced by local time-stepping schemes. Extension of the vectorization procedure for predicting 2D viscous and 3D inviscid flows are demonstrated. Conclusions are drawn regarding the applicability of the finite element procedures for realistic problems that require hundreds of thousands of nodes

Thermal stresses in rectangular plates: variational and finite element solutions

International Nuclear Information System (INIS)

Laura, P.A.A.; Gutierrez, R.H.; Sanchez Sarmiento, G.; Basombrio, F.G.

1978-01-01

This paper deals with the development of an approximate method for the analysis of thermal stresses in rectangular plates (plane stress problem) and an evaluation of the relative accuracy of the finite element method. The stress function is expanded in terms of polynomial coordinate functions which identically satisfy the boundary conditions, and a variational approach is used to determine the expansion coefficients. The results are in good agreement with a finite element approach. (Auth.)

Coupled porohyperelastic mass transport (PHEXPT) finite element models for soft tissues using ABAQUS.

Science.gov (United States)

Vande Geest, Jonathan P; Simon, B R; Rigby, Paul H; Newberg, Tyler P

2011-04-01

Finite element models (FEMs) including characteristic large deformations in highly nonlinear materials (hyperelasticity and coupled diffusive/convective transport of neutral mobile species) will allow quantitative study of in vivo tissues. Such FEMs will provide basic understanding of normal and pathological tissue responses and lead to optimization of local drug delivery strategies. We present a coupled porohyperelastic mass transport (PHEXPT) finite element approach developed using a commercially available ABAQUS finite element software. The PHEXPT transient simulations are based on sequential solution of the porohyperelastic (PHE) and mass transport (XPT) problems where an Eulerian PHE FEM is coupled to a Lagrangian XPT FEM using a custom-written FORTRAN program. The PHEXPT theoretical background is derived in the context of porous media transport theory and extended to ABAQUS finite element formulations. The essential assumptions needed in order to use ABAQUS are clearly identified in the derivation. Representative benchmark finite element simulations are provided along with analytical solutions (when appropriate). These simulations demonstrate the differences in transient and steady state responses including finite deformations, total stress, fluid pressure, relative fluid, and mobile species flux. A detailed description of important model considerations (e.g., material property functions and jump discontinuities at material interfaces) is also presented in the context of finite deformations. The ABAQUS-based PHEXPT approach enables the use of the available ABAQUS capabilities (interactive FEM mesh generation, finite element libraries, nonlinear material laws, pre- and postprocessing, etc.). PHEXPT FEMs can be used to simulate the transport of a relatively large neutral species (negligible osmotic fluid flux) in highly deformable hydrated soft tissues and tissue-engineered materials.

A finite element modeling of a multifunctional hybrid composite beam with viscoelastic materials

Science.gov (United States)

Wang, Ya; Inman, Daniel J.

2013-04-01

The multifunctional hybrid composite structure studied here consists of a ceramic outer layer capable of withstanding high temperatures, a functionally graded ceramic layer combining shape memory alloy (SMA) properties of NiTi together with Ti2AlC (called Graded Ceramic/Metal Composite, or GCMeC), and a high temperature sensor patch, followed by a polymer matrix composite laced with vascular cooling channels all held together with various epoxies. Due to the recoverable nature of SMA and adhesive properties of Ti2AlC, the damping behavior of the GCMeC is largely viscoelastic. This paper presents a finite element formulation for this multifunctional hybrid structure with embedded viscoelastic material. In order to implement the viscoelastic model into the finite element formulation, a second order three parameter Golla-Hughes-McTavish (GHM) method is used to describe the viscoelastic behavior. Considering the parameter identification, a strategy to estimate the fractional order of the time derivative and the relaxation time is outlined. The curve-fitting aspects of both GHM and ADF show good agreement with experimental data obtained from dynamic mechanics analysis. The performance of the finite element of the layered multifunctional beam is verified through experimental model analysis.

Bolted Ribs Analysis for the ITER Vacuum Vessel using Finite Element Submodelling Techniques

Energy Technology Data Exchange (ETDEWEB)

Zarzalejos, José María, E-mail: jose.zarzalejos@ext.f4e.europa.eu [External at F4E, c/Josep Pla, n.2, Torres Diagonal Litoral, Edificio B3, E-08019, Barcelona (Spain); Fernández, Elena; Caixas, Joan; Bayón, Angel [F4E, c/Josep Pla, n.2, Torres Diagonal Litoral, Edificio B3, E-08019, Barcelona (Spain); Polo, Joaquín [Iberdrola Ingeniería y Construcción, Avenida de Manoteras 20, 28050 Madrid (Spain); Guirao, Julio [Numerical Analysis Technologies, S L., Marqués de San Esteban 52, Entlo, 33209 Gijon (Spain); García Cid, Javier [Iberdrola Ingeniería y Construcción, Avenida de Manoteras 20, 28050 Madrid (Spain); Rodríguez, Eduardo [Mechanical Engineering Department EPSIG, University of Oviedo, Gijon (Spain)

2014-10-15

Highlights: • The ITER Vacuum Vessel Bolted Ribs assemblies are modelled using Finite Elements. • Finite Element submodelling techniques are used. • Stress results are obtained for all the assemblies and a post-processing is performed. • All the elements of the assemblies are compliant with the regulatory provisions. • Submodelling is a time-efficient solution to verify the structural integrity of this type of structures. - Abstract: The ITER Vacuum Vessel (VV) primary function is to enclose the plasmas produced by the ITER Tokamak. Since it acts as the first radiological barrier of the plasma, it is classified as a class 2 welded box structure, according to RCC-MR 2007. The VV is made of an inner and an outer D-shape, 60 mm-thick double shell connected through thick massive bars (housings) and toroidal and poloidal structural stiffening ribs. In order to provide neutronic shielding to the ex-vessel components, the space between shells is filled with borated steel plates, called In-Wall Shielding (IWS) blocks, and water. In general, these blocks are connected to the IWS ribs which are connected to adjacent housings. The development of a Finite Element model of the ITER VV including all its components in detail is unaffordable from the computational point of view due to the large number of degrees of freedom it would require. This limitation can be overcome by using submodelling techniques to simulate the behaviour of the bolted ribs assemblies. Submodelling is a Finite Element technique which allows getting more accurate results in a given region of a coarse model by generating an independent, finer model of the region under study. In this paper, the methodology and several simulations of the VV bolted ribs assemblies using submodelling techniques are presented. A stress assessment has been performed for the elements involved in the assembly considering possible types of failure and including stress classification and categorization techniques to analyse

Finite element analysis of three dimensional crack growth by the use of a boundary element sub model

DEFF Research Database (Denmark)

Lucht, Tore

2009-01-01

A new automated method to model non-planar three dimensional crack growth is proposed which combines the advantages of both the boundary element method and the finite element method. The proposed method links the two methods by a submodelling strategy in which the solution of a global finite...... element model containing an approximation of the crack is interpolated to a much smaller boundary element model containing a fine discretization of the real crack. The method is validated through several numerical comparisons and by comparison to crack growth measured in a test specimen for an engineering...

Real-time volumetric deformable models for surgery simulation using finite elements and condensation

DEFF Research Database (Denmark)

Bro-Nielsen, Morten; Cotin, S.

1996-01-01

This paper discusses the application of SD solid volumetric Finite Element models to surgery simulation. In particular it introduces three new ideas for solving the problem of achieving real-time performance for these models. The simulation system we have developed is described and we demonstrate...

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

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

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.

A three-dimensional thermal finite element analysis of AISI 304 stainless steel and copper dissimilar weldment

Science.gov (United States)

Singh, Gurdeep; Saxena, Ravindra K.; Pandey, Sunil

2018-04-01

The aim of this study to developed a 3-D thermal finite element model for dissimilar material welding of AISI-304 stainless steel and copper. Welding of similar material is widely studied using experimental and numerical methods but the problem becomes trivial for the welding of dissimilar materials especially in ferrous and nonferrous materials. Finite element analysis of dissimilar material welding is a cost-effective method for the understanding and analysis of the process. The finite element analysis has been performed to predict the heat affected zone and temperature distribution in AISI-304 stainless steel and copper dissimilar weldment using MSC Marc 2017®. Due to the difference in physical properties of these materials the behavior of heat affected zone and temperature distribution are perceived to be different. To verify the accuracy of the thermal finite element model, the welding process was simulated with butt-welded joints having same dimensions and parameters from Attarha and Far [1]. It is found from the study that the heat affected zone is larger in copper weld pads than in AISI 304 stainless steel due to large difference in thermal conductivity of these two weld pads.

Finite element analysis of residual stress in plasma-sprayed ceramic

International Nuclear Information System (INIS)

Mullen, R.L.; Hendricks, R.C.; McDonald, G.

1985-01-01

Residual stress in a ZrO 2 -Y 2 O 3 ceramic coating resulting from the plasma spraying operation is calculated. The calculations were done using the finite element method. Both thermal and mechanical analysis were performed. The resulting residual stress field was compared to the measurements obtained by Hendricks and McDonald. Reasonable agreement between the predicted and measured moment occurred. However, the resulting stress field is not in pure bending

Coupled Analytical-Finite Element Methods for Linear Electromagnetic Actuator Analysis

Directory of Open Access Journals (Sweden)

K. Srairi

2005-09-01

Full Text Available In this paper, a linear electromagnetic actuator with moving parts is analyzed. The movement is considered through the modification of boundary conditions only using coupled analytical and finite element analysis. In order to evaluate the dynamic performance of the device, the coupling between electric, magnetic and mechanical phenomena is established. The displacement of the moving parts and the inductor current are determined when the device is supplied by capacitor discharge voltage.

Eigensolution of finite element problems in a completely connected parallel architecture

Science.gov (United States)

Akl, Fred A.; Morel, Michael R.

1989-01-01

A parallel algorithm for the solution of the generalized eigenproblem in linear elastic finite element analysis, (K)(phi)=(M)(phi)(omega), where (K) and (M) are of order N, and (omega) is of order q is presented. The parallel algorithm is based on a completely connected parallel architecture in which each processor is allowed to communicate with all other processors. The algorithm has been successfully implemented on a tightly coupled multiple-instruction-multiple-data (MIMD) parallel processing computer, Cray X-MP. A finite element model is divided into m domains each of which is assumed to process n elements. Each domain is then assigned to a processor, or to a logical processor (task) if the number of domains exceeds the number of physical processors. The macro-tasking library routines are used in mapping each domain to a user task. Computational speed-up and efficiency are used to determine the effectiveness of the algorithm. The effect of the number of domains, the number of degrees-of-freedom located along the global fronts and the dimension of the subspace on the performance of the algorithm are investigated. For a 64-element rectangular plate, speed-ups of 1.86, 3.13, 3.18 and 3.61 are achieved on two, four, six and eight processors, respectively.

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.

A non conforming finite element method for computing eigenmodes of resonant cavities

International Nuclear Information System (INIS)

Touze, F.; Le Meur, G.

1990-06-01

We present here a non conforming finite element in R 3 . This finite element, built on tetrahedrons, is particularly suited for computing eigenmodes. The main advantage of this element is that it preserves some structural properties of the space in which the solutions of the Maxwell's equations are to be found. Numerical results are presented for both two-dimensional and three-dimensional cases

Finite element method - theory and applications

International Nuclear Information System (INIS)

Baset, S.

1992-01-01

This paper summarizes the mathematical basis of the finite element method. Attention is drawn to the natural development of the method from an engineering analysis tool into a general numerical analysis tool. A particular application to the stress analysis of rubber materials is presented. Special advantages and issues associated with the method are mentioned. (author). 4 refs., 3 figs

Application of Finite Element Method to Determine the Performances of a Permanent Magnet Synchronous Motor for Driving a Bicycle

Directory of Open Access Journals (Sweden)

Nicolae Digă

2014-09-01

Full Text Available In this paper, the authors present a case study in which was analyzed by finite element method a permanent magnet synchronous motor for driving a bicycle using the analysis and simulation software ANSYS Electromagnetics Low Frequency of ANSYS Inc. Company. Modelling and simulation with ANSYS ® Maxwell 2D of electromagnetic field in the studied motor was conducted for different initial positions (internal angle rotor-stator configured (set δ1. It was identified the internal angle for which the performances of PMSM are very close to those obtained by computation.

Finite groups with three conjugacy class sizes of some elements

Indian Academy of Sciences (India)

Conjugacy class sizes; p-nilpotent groups; finite groups. 1. Introduction. All groups ... group G has exactly two conjugacy class sizes of elements of prime power order. .... [5] Huppert B, Character Theory of Finite Groups, de Gruyter Exp. Math.

A least squares principle unifying finite element, finite difference and nodal methods for diffusion theory

International Nuclear Information System (INIS)

Ackroyd, R.T.

1987-01-01

A least squares principle is described which uses a penalty function treatment of boundary and interface conditions. Appropriate choices of the trial functions and vectors employed in a dual representation of an approximate solution established complementary principles for the diffusion equation. A geometrical interpretation of the principles provides weighted residual methods for diffusion theory, thus establishing a unification of least squares, variational and weighted residual methods. The complementary principles are used with either a trial function for the flux or a trial vector for the current to establish for regular meshes a connection between finite element, finite difference and nodal methods, which can be exact if the mesh pitches are chosen appropriately. Whereas the coefficients in the usual nodal equations have to be determined iteratively, those derived via the complementary principles are given explicitly in terms of the data. For the further development of the connection between finite element, finite difference and nodal methods, some hybrid variational methods are described which employ both a trial function and a trial vector. (author)

Dislocation dynamics in non-convex domains using finite elements with embedded discontinuities

Science.gov (United States)

Romero, Ignacio; Segurado, Javier; LLorca, Javier

2008-04-01

The standard strategy developed by Van der Giessen and Needleman (1995 Modelling Simul. Mater. Sci. Eng. 3 689) to simulate dislocation dynamics in two-dimensional finite domains was modified to account for the effect of dislocations leaving the crystal through a free surface in the case of arbitrary non-convex domains. The new approach incorporates the displacement jumps across the slip segments of the dislocations that have exited the crystal within the finite element analysis carried out to compute the image stresses on the dislocations due to the finite boundaries. This is done in a simple computationally efficient way by embedding the discontinuities in the finite element solution, a strategy often used in the numerical simulation of crack propagation in solids. Two academic examples are presented to validate and demonstrate the extended model and its implementation within a finite element program is detailed in the appendix.

Dislocation dynamics in non-convex domains using finite elements with embedded discontinuities

International Nuclear Information System (INIS)

Romero, Ignacio; Segurado, Javier; LLorca, Javier

2008-01-01

The standard strategy developed by Van der Giessen and Needleman (1995 Modelling Simul. Mater. Sci. Eng. 3 689) to simulate dislocation dynamics in two-dimensional finite domains was modified to account for the effect of dislocations leaving the crystal through a free surface in the case of arbitrary non-convex domains. The new approach incorporates the displacement jumps across the slip segments of the dislocations that have exited the crystal within the finite element analysis carried out to compute the image stresses on the dislocations due to the finite boundaries. This is done in a simple computationally efficient way by embedding the discontinuities in the finite element solution, a strategy often used in the numerical simulation of crack propagation in solids. Two academic examples are presented to validate and demonstrate the extended model and its implementation within a finite element program is detailed in the appendix

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

Mathematical aspects of finite element methods for incompressible viscous flows

Science.gov (United States)

Gunzburger, M. D.

1986-01-01

Mathematical aspects of finite element methods are surveyed for incompressible viscous flows, concentrating on the steady primitive variable formulation. The discretization of a weak formulation of the Navier-Stokes equations are addressed, then the stability condition is considered, the satisfaction of which insures the stability of the approximation. Specific choices of finite element spaces for the velocity and pressure are then discussed. Finally, the connection between different weak formulations and a variety of boundary conditions is explored.

Extension to linear dynamics for hybrid stress finite element formulation based on additional displacements

Science.gov (United States)

Sumihara, K.

Based upon legitimate variational principles, one microscopic-macroscopic finite element formulation for linear dynamics is presented by Hybrid Stress Finite Element Method. The microscopic application of Geometric Perturbation introduced by Pian and the introduction of infinitesimal limit core element (Baby Element) have been consistently combined according to the flexible and inherent interpretation of the legitimate variational principles initially originated by Pian and Tong. The conceptual development based upon Hybrid Finite Element Method is extended to linear dynamics with the introduction of physically meaningful higher modes.

Evaluation of Round Window Stimulation Performance in Otosclerosis Using Finite Element Modeling

Directory of Open Access Journals (Sweden)

Shanguo Yang

2016-01-01

Full Text Available Round window (RW stimulation is a new type of middle ear implant’s application for treating patients with middle ear disease, such as otosclerosis. However, clinical outcomes show a substantial degree of variability. One source of variability is the variation in the material properties of the ear components caused by the disease. To investigate the influence of the otosclerosis on the performance of the RW stimulation, a human ear finite element model including middle ear and cochlea was established based on a set of microcomputerized tomography section images of a human temporal bone. Three characteristic changes of the otosclerosis in the auditory system were simulated in the FE model: stapedial annular ligament stiffness enlargement, stapedial abnormal bone growth, and partial fixation of the malleus. The FE model was verified by comparing the model-predicted results with published experimental measurements. The equivalent sound pressure (ESP of RW stimulation was calculated via comparing the differential intracochlear pressure produced by the RW stimulation and the normal eardrum sound stimulation. The results show that the increase of stapedial annular ligament and partial fixation of the malleus decreases RW stimulation’s ESP prominently at lower frequencies. In contrast, the stapedial abnormal bone growth deteriorates RW stimulation’s ESP severely at higher frequencies.

Finite element modelling of aluminum alloy 2024-T3 under transverse impact loading

Science.gov (United States)

Abdullah, Ahmad Sufian; Kuntjoro, Wahyu; Yamin, A. F. M.

2017-12-01

Fiber metal laminate named GLARE is a new aerospace material which has great potential to be widely used in future lightweight aircraft. It consists of aluminum alloy 2024-T3 and glass-fiber reinforced laminate. In order to produce reliable finite element model of impact response or crashworthiness of structure made of GLARE, one can initially model and validate the finite element model of the impact response of its constituents separately. The objective of this study was to develop a reliable finite element model of aluminum alloy 2024-T3 under low velocity transverse impact loading using commercial software ABAQUS. Johnson-Cook plasticity and damage models were used to predict the alloy's material properties and impact behavior. The results of the finite element analysis were compared to the experiment that has similar material and impact conditions. Results showed good correlations in terms of impact forces, deformation and failure progressions which concluded that the finite element model of 2024-T3 aluminum alloy under low velocity transverse impact condition using Johnson-Cook plastic and damage models was reliable.

An Eulerian-Lagrangian finite-element method for modeling crack growth in creeping materials

International Nuclear Information System (INIS)

Lee Hae Sung.

1991-01-01

This study is concerned with the development of finite-element-solution methods for analysis of quasi-static, ductile crack growth in history-dependent materials. The mixed Eulerian-Langrangian description (ELD) kinematic model is shown to have several desirable properties for modeling inelastic crack growth. Accordingly, a variational statement based on the ELD for history-dependent materials is developed, and a new moving-grid finite-element method based on the variational statement is presented. The moving-grid finite-element method based on the variational statement is presented. The moving-grid finite-element method is applied to the analysis of transient, quasi-static, mode-III crack growth in creeping materials. A generalized Petrov-Galerkin method (GPG) is developed that simultaneously stabilizes the statement to admit L 2 basis functions for the nonlinear strain field. Quasi-static, model-III crack growth in creeping materials under small-scale-yielding (SSY) conditions is considered. The GPG/ELD moving-grid finite-element formulation is used to model a transient crack-growth problem. The GPG/ELD results compare favorably with previously-published numerical results and the asymptotic solutions

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 finite element method for neutron transport

International Nuclear Information System (INIS)

Ackroyd, R.T.

1983-01-01

A completely boundary-free maximum principle for the first-order Boltzmann equation is derived from the completely boundary-free maximum principle for the mixed-parity Boltzmann equation. When continuity is imposed on the trial function for directions crossing interfaces the completely boundary-free principle for the first-order Boltzmann equation reduces to a maximum principle previously established directly from first principles and indirectly by the Euler-Lagrange method. Present finite element methods for the first-order Boltzmann equation are based on a weighted-residual method which permits the use of discontinuous trial functions. The new principle for the first-order equation can be used as a basis for finite-element methods with the same freedom from boundary conditions as those based on the weighted-residual method. The extremum principle as the parent of the variationally-derived weighted-residual equations ensures their good behaviour. (author)

Isogeometric finite element data structures based on Bézier extraction of T-splines

NARCIS (Netherlands)

Scott, M.A.; Borden, M.J.; Verhoosel, C.V.; Sederberg, T.W.; Hughes, T.J.R.

2011-01-01

We develop finite element data structures for T-splines based on Bézier extraction generalizing our previous work for NURBS. As in traditional finite element analysis, the extracted Bézier elements are defined in terms of a fixed set of polynomial basis functions, the so-called Bernstein basis. The

Rational bases and generalized barycentrics applications to finite elements and graphics

CERN Document Server

Wachspress, Eugene

2016-01-01

This three-part volume explores theory for construction of rational interpolation functions for continuous patchwork approximation.  Authored by the namesake of the Wachspress Coordinates, the book develops construction of basis functions for a broad class of elements which have widespread graphics and finite element application. Part one is the 1975 book A Rational Finite Element Basis (with minor updates and corrections) written by Dr. Wachspress.  Part two describes theoretical advances since 1975 and includes analysis of elements not considered previously.  Part three consists of annotated MATLAB programs implementing theory presented in parts one and two.

2D and 3D finite element analysis of buffer-backfill interaction

International Nuclear Information System (INIS)

Leoni, M.

2013-08-01

Methods for backfilling and sealing of disposal tunnels in an underground repository for spent nuclear fuel are studied in cooperation between Finland (Posiva Oy) and Sweden (Svensk Kaernbraenslehantering AB, SKB) in 'BAckfilling and CLOsure of the deep repository' (Baclo) programme. Baclo phase III included modelling task force SP1: Finite element modelling of deformation of the backfill due to swelling of the buffer. The objective of the finite element modelling of the backfill was to study the interaction between the buffer and backfilling. The calculations aimed to find out how large deformations can happen in the buffer-backfill interface causing loosening of the buffer bentonite above the canister. The criterion used was that the saturated density of the buffer right above the canister should be higher than 1990 kg/m 3 . This report presents the results of finite element numerical analyses carried out by Wesi Geotecnica Srl. The modelling calculations were conducted with the so-called OL1-2 deposition tunnel geometry (Juvankoski 2009). Several parameters have been considered, varying from geometry variations to different mechanical constitutive models for different components of the model. In all analyses it has been assumed that the buffer material is fully saturated, thus exerting the isotropic swelling pressure estimated in the range 7 MPa .. 15 MPa, against a fully-dry backfill, which is no doubt the 'worst case scenario' with the highest risk to lead in decrease in dry density of the buffer. Friedland clay has been considered for backfill blocks and 30/70 mixture for foundation bed on which backfill blocks are installed. Preliminarily, finite element analyses have been performed with newly released PLAXIS 2D 2010 within the assumption of axial symmetry, the purpose of this first set of calculations being the evaluation of most relevant parameters influencing the deformations of buffer material. Hence, full 3D calculations have been performed with PLAXIS

Transport and dispersion of pollutants in surface impoundments: a finite element model

International Nuclear Information System (INIS)

Yeh, G.T.

1980-07-01

A surface impoundment model in finite element (SIMFE) is presented to enable the simulation of flow circulations and pollutant transport and dispersion in natural or artificial lakes, reservoirs or ponds with any number of islands. This surface impoundment model consists of two sub-models: hydrodynamic and pollutant transport models. Both submodels are simulated by the finite element method. While the hydrodynamic model is solved by the standard Galerkin finite element scheme, the pollutant transport model can be solved by any of the twelve optional finite element schemes built in the program. Theoretical approximations and the numerical algorithm of SIMFE are described. Detail instruction of the application are given and listing of FORTRAN IV source program are provided. Two sample problems are given. One is for an idealized system with a known solution to show the accuracy and partial validation of the models. The other is applied to Prairie Island for a set of hypothetical input data, typifying a class of problems to which SIMFE may be applied

Numerical simulation of mechatronic sensors and actuators finite elements for computational multiphysics

CERN Document Server

Kaltenbacher, Manfred

2015-01-01

Like the previous editions also the third edition of this book combines the detailed physical modeling of mechatronic systems and their precise numerical simulation using the Finite Element (FE) method. Thereby, the basic chapter concerning the Finite Element (FE) method is enhanced, provides now also a description of higher order finite elements (both for nodal and edge finite elements) and a detailed discussion of non-conforming mesh techniques. The author enhances and improves many discussions on principles and methods. In particular, more emphasis is put on the description of single fields by adding the flow field. Corresponding to these field, the book is augmented with the new chapter about coupled flow-structural mechanical systems. Thereby, the discussion of computational aeroacoustics is extended towards perturbation approaches, which allows a decomposition of flow and acoustic quantities within the flow region. Last but not least, applications are updated and restructured so that the book meets mode...

Transport and dispersion of pollutants in surface impoundments: a finite element model

Energy Technology Data Exchange (ETDEWEB)

Yeh, G.T.

1980-07-01

A surface impoundment model in finite element (SIMFE) is presented to enable the simulation of flow circulations and pollutant transport and dispersion in natural or artificial lakes, reservoirs or ponds with any number of islands. This surface impoundment model consists of two sub-models: hydrodynamic and pollutant transport models. Both submodels are simulated by the finite element method. While the hydrodynamic model is solved by the standard Galerkin finite element scheme, the pollutant transport model can be solved by any of the twelve optional finite element schemes built in the program. Theoretical approximations and the numerical algorithm of SIMFE are described. Detail instruction of the application are given and listing of FORTRAN IV source program are provided. Two sample problems are given. One is for an idealized system with a known solution to show the accuracy and partial validation of the models. The other is applied to Prairie Island for a set of hypothetical input data, typifying a class of problems to which SIMFE may be applied.

Finite element analysis on badminton racket design parameters

CERN Document Server

Nasruddin, Fakhrizal Azmy; Syahrom, Ardiyansyah; Abdul Kadir, Mohammed Rafiq; Omar, Abdul Hafidz; Öchsner, Andreas

2016-01-01

This work identifies the characteristics of racket design parameters that influence racket performance.  It presents the finite element analysis of several designs of badminton rackets and compares them to experimental results for validation. Designing a racket requires a comprehensive understanding of racket performance characteristics. Essentially, racket performance is related to the sweet spot, which is the spot on the racket head that produces the most power and control when it strikes a shuttlecock. Determining a coefficient of restitution can help to identify the sweet spot on a racket. By analyzing several head shape designs, it becomes apparent that isometric head shape rackets produce better coefficients of restitution compared to oval and round ones. It is recommended that the racket design consist of low string tension, stiffer racket shafts and bigger head size in order to produce higher shuttlecock speed.

Development and validation of a weight-bearing finite element model for total knee replacement.

Science.gov (United States)

Woiczinski, M; Steinbrück, A; Weber, P; Müller, P E; Jansson, V; Schröder, Ch

2016-01-01

Total knee arthroplasty (TKA) is a successful procedure for osteoarthritis. However, some patients (19%) do have pain after surgery. A finite element model was developed based on boundary conditions of a knee rig. A 3D-model of an anatomical full leg was generated from magnetic resonance image data and a total knee prosthesis was implanted without patella resurfacing. In the finite element model, a restarting procedure was programmed in order to hold the ground reaction force constant with an adapted quadriceps muscle force during a squat from 20° to 105° of flexion. Knee rig experimental data were used to validate the numerical model in the patellofemoral and femorotibial joint. Furthermore, sensitivity analyses of Young's modulus of the patella cartilage, posterior cruciate ligament (PCL) stiffness, and patella tendon origin were performed. Pearson's correlations for retropatellar contact area, pressure, patella flexion, and femorotibial ap-movement were near to 1. Lowest root mean square error for retropatellar pressure, patella flexion, and femorotibial ap-movement were found for the baseline model setup with Young's modulus of 5 MPa for patella cartilage, a downscaled PCL stiffness of 25% compared to the literature given value and an anatomical origin of the patella tendon. The results of the conducted finite element model are comparable with the experimental results. Therefore, the finite element model developed in this study can be used for further clinical investigations and will help to better understand the clinical aspects after TKA with an unresurfaced patella.