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)

Elastic-plastic and creep analyses by assumed stress finite elements

International Nuclear Information System (INIS)

Pian, T.H.H.; Spilker, R.L.; Lee, S.W.

1975-01-01

A formulation is presented of incremental finite element solutions for both initial stress and initial strain problems based on modified complementary energy principle with relaxed inter-element continuity requirement. The corresponding finite element model is the assumed stress hybrid model which has stress parameters in the interior of each element and displacements at the individual nodes as unknowns. The formulation includes an important consideration that the states of stress and strain and the beginning of each increment may not satisfy the equilibrium and compatibility equations. These imbalance and mismatch conditions all lead to correction terms for the equivalent nodal forces of the matrix equations. The initial stress method is applied to elastic-plastic analysis of structures. In this case the stress parameters for the individual elements can be eliminated resulting to a system of equations with only nodal displacements as unknowns. Two different complementary energy principles can be formulated, in one of which the equilibrium of the final state of stress is maintained while in the other the equilibrium of the stress increments is maintained. Each of these two different formulations can be combined with different iterative schemes to be used at each incremental steps of the elastic-plastic analysis. It is also indicated clearly that for the initial stress method the state of stress at the beginning of each increments is in general, not in equilibrium and an imbalance correction is needed. Results of a comprehensive evaluation of various solution procedures by the initial stress method using the assumed stress hybrid elements are presented. The example used is the static response of a thick wall cylinder of elastic-perfectly plastic material under internal pressure. Solid of revolution elements with rectangular cross sections are used

Finite element analysis of prestressed concrete reactor vessels

International Nuclear Information System (INIS)

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

1977-01-01

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

Creep analysis of adhesively bonded single lap joint using finite element method

International Nuclear Information System (INIS)

Zehsaz, Mohammad; Vakilitahami, Tahami Farid; Saeimisadigh, Mohammad Ali

2014-01-01

Adhesive joints are being used widely in engineering industries due to the increasing demand for designing lightweight structures. Because of the physical properties of the most adhesives, they creep even at room temperature. Therefore, the creep behavior of a single lap adhesive joint is studied in this paper. For this purpose, using the experimental data, creep constitutive equations for the adhesive has been obtained. Then, these equations have been employed to investigate the creep behavior of the joint. The results show that due to the creep straining, the stresses in the joint corners, decrease. However, creep strain accumulates in these areas which this in turn may lead to separation of adhesive from adherent. In order to eliminate the effect of strain accumulation, two modifying methods have been proposed in this paper: increasing the layer thickness and using filleted joints.

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 finite element method with contact for tensile analysis in fuel rods

International Nuclear Information System (INIS)

Tanajura, C.A.S.; Galeao, A.C.N.R.

1987-01-01

Elements for mechanical analysis of fuel rod of a PWR type reactor, are presented. The rod, consists basically in a cylindrical coating of zircalloy which contains pilling of UO 2 pellets, is submitted to strong internal and external pressures, intense temperature gradients and neutron flux. These conditions lead several phenomena in the pellet (swelling, fracture, densification, creep) and in the cladding (embrittlement, corrosion, creep) which undergo deformations leading them to contact the restriction for the interpenetration is included in the problem without restriction by Lagrange multipliers. Considering a non-linear problem, due to the surface of contact to be not known a priori, the numerical solutions were obtained using the finite element method. (M.C.K.) [pt

Inelastic finite element cyclic analysis of a nozzle-to-cylinder intersection

International Nuclear Information System (INIS)

Barsoum, R.S.; Loomis, R.W.

1976-01-01

A finite element elastic-plastic and creep analysis of a nozzle-to-cylinder intersection subject to cyclic thermal shock, internal pressure, and mechanical loads is presented. The nozzle configuration is that of the intermediate heat exchanger (IHX) for the Fast Flux Test Facility (FFTF). The analysis was performed using the general purpose program MARC. Both the elastic and inelastic results of the analysis are presented. The intention of this study to analytically investigate the applicability of simplified ratchetting and creep-fatigue rules for LMFBR components, as a part of a program covering various geometries and loadings

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

International Nuclear Information System (INIS)

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

1995-01-01

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

Viscoelastic creep of high-temperature concrete

International Nuclear Information System (INIS)

Pfeiffer, P.A.; Marchertas, A.H.; Bazant, Z.P.

1985-01-01

Presented in this report is the analytical model for analysis of high temperature creep response of concrete. The creep law used is linear (viscoelastic), the temperature and moisture effects on the creep rate and also aging are included. Both constant and transient temperature as well as constant and transient moisture conditions are considered. Examples are presented to correlate experimental data with parameters of the analytical model by the use of a finite element scheme

Finite element analysis of structures at high temperatures with special application to plane steel beams and frames

International Nuclear Information System (INIS)

Peterson, A.

1984-01-01

Nonlinear analysis of structures at high temperatures is studied. Both geometric and material nonlinearities are taken into account. Continuum mechanics relations are used to derive general finite element equations. An alternative formulation to Total Lagrangian (TL) and Updated Lagrangian (UL) formulations named Partially updated Lagrangian (PL) formulation is presented. An isotropic small strain constitutive model using the von Mises yield criterion is derived for high temperature conditions. The model developed can be characterized as combined elastic-plastic-viscoplastic. The strain components are treated separately but plastic strains and viscoplastic (creep)strains are allowed to interact. A new formulation of the creep behaviour is given. Both primary and secondary creep are considered. As an application of the derived finite element equations and the constitutive model steel beams and frames are studied. The theory is implemented in a computer program, CAMFEM. The program operates on a command language with possibilities to store user-defined matrices on files and to create macro commands. Comparison with experimental observation shows that the present theory well describes experimentally observed phenomena. (author)

Engineering C-integral estimates for generalised creep behaviour and finite element validation

International Nuclear Information System (INIS)

Kim, Yun-Jae; Kim, Jin-Su; Huh, Nam-Su; Kim, Young-Jin

2002-01-01

This paper proposes an engineering method to estimate the creep C-integral for realistic creep laws to assess defective components operating at elevated temperatures. The proposed estimation method is mainly for the steady-state C * -integral, but a suggestion is also given for estimating the transient C(t)-integral. The reference stress approach is the basis of the proposed equation, but an enhancement in terms of accuracy is made through the definition of the reference stress. The proposed estimation equations are compared with extensive elastic-creep FE results employing various creep-deformation constitutive laws for six different geometries, including two-dimensional, axi-symmetric and three-dimensional geometries. Overall good agreement between the proposed method and the FE results provides confidence in the use of the proposed method for defect assessment of components at elevated temperatures. Moreover, it is shown that for surface cracks the proposed method can be used to estimate C * at any location along the crack front

Elasto-viscoplastic finite element model for prestressed concrete structures

International Nuclear Information System (INIS)

Prates Junior, N.P.; Silva, C.S.B.; Campos Filho, A.; Gastal, F.P.S.L.

1995-01-01

This paper presents a computational model, based on the finite element method, for the study of reinforced and prestressed concrete structures under plane stress states. It comprehends short and long-term loading situations, where creep and shrinkage in concrete and steel relaxation are considered. Elasto-viscoplastic constitutive models are used to describe the behavior of the materials. The model includes prestressing and no prestressing reinforcement, on situation with pre- and post-tension with and without bond. A set of prestressed concrete slab elements were tested under instantaneous and long-term loading. The experimental data for deflections, deformations and ultimate strength are used to compare and validate the results obtained through the proposed model. (author). 11 refs., 5 figs

Boundary element method for modelling creep behaviour

International Nuclear Information System (INIS)

Zarina Masood; Shah Nor Basri; Abdel Majid Hamouda; Prithvi Raj Arora

2002-01-01

A two dimensional initial strain direct boundary element method is proposed to numerically model the creep behaviour. The boundary of the body is discretized into quadratic element and the domain into quadratic quadrilaterals. The variables are also assumed to have a quadratic variation over the elements. The boundary integral equation is solved for each boundary node and assembled into a matrix. This matrix is solved by Gauss elimination with partial pivoting to obtain the variables on the boundary and in the interior. Due to the time-dependent nature of creep, the solution has to be derived over increments of time. Automatic time incrementation technique and backward Euler method for updating the variables are implemented to assure stability and accuracy of results. A flowchart of the solution strategy is also presented. (Author)

Analysis of indentation creep

Science.gov (United States)

Don S. Stone; Joseph E. Jakes; Jonathan Puthoff; Abdelmageed A. Elmustafa

2010-01-01

Finite element analysis is used to simulate cone indentation creep in materials across a wide range of hardness, strain rate sensitivity, and work-hardening exponent. Modeling reveals that the commonly held assumption of the hardness strain rate sensitivity (mΗ) equaling the flow stress strain rate sensitivity (mσ...

A Numerical Procedure to Obtain the Creep Parameters of the Thermal Barrier Coating

Directory of Open Access Journals (Sweden)

Shifeng Wen

2014-05-01

Full Text Available Three-point bending creep test was used to understand the creep behavior of typical thin film/substrate systems—thermal barrier coating (TBC systems. Firstly, a simplified model, which does not consider the local effect, has been set up to get an analytical relationship. The important result is that creep stress exponent of materials is equal to the creep load exponent of the steady-state deflection rate of BC specimens. Secondly, in order to consider the local effect of bending, the finite element method (FEM has been carried out. FEM calculation shows that there is a steady stage of the creep deflection under a constant applied load. And the exponent of the steady-state creep deflection rate to the applied load is found to be equal to the creep stress exponent of materials. The creep constant of the materials can be obtained by a set of trials with assumed creep constants of materials and can be finally determined by the best fit method. Finally, the finite element results show that the influences of the friction, the thickness of TBCs, and the modulus ratio of TBC to the substrate on stress distribution are important.

Basic Finite Element Method

International Nuclear Information System (INIS)

Lee, Byeong Hae

1992-02-01

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

FEAST: a two-dimensional non-linear finite element code for calculating stresses

International Nuclear Information System (INIS)

Tayal, M.

1986-06-01

The computer code FEAST calculates stresses, strains, and displacements. The code is two-dimensional. That is, either plane or axisymmetric calculations can be done. The code models elastic, plastic, creep, and thermal strains and stresses. Cracking can also be simulated. The finite element method is used to solve equations describing the following fundamental laws of mechanics: equilibrium; compatibility; constitutive relations; yield criterion; and flow rule. FEAST combines several unique features that permit large time-steps in even severely non-linear situations. The features include a special formulation for permitting many finite elements to simultaneously cross the boundary from elastic to plastic behaviour; accomodation of large drops in yield-strength due to changes in local temperature and a three-step predictor-corrector method for plastic analyses. These features reduce computing costs. Comparisons against twenty analytical solutions and against experimental measurements show that predictions of FEAST are generally accurate to ± 5%

Implementation of creep-fatigue model into finite-element code to assess cooled turbine blade.

CSIR Research Space (South Africa)

Dedekind, MO

1994-01-01

Full Text Available Turbine blades which are designed with airfoil cooling are subject to thermo-mechanical fatigue as well as creep damage. These problems arise due to thermal cycling and high operating temperatures in service. An implementation of fatigue and creep...

Generalized finite elements

International Nuclear Information System (INIS)

Wachspress, E.

2009-01-01

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

Effects of segregation of primary alloying elements on the creep response in magnesium alloys

DEFF Research Database (Denmark)

Huang, Y.D.; Dieringa, H.; Hort, N.

2008-01-01

The segregation of primary alloying elements deteriorates the high temperature creep resistance of magnesium alloys. Annealing at high temperatures alleviating their segregations can improve the creep resistance. Present investigation on the effect of segregation of primary alloying elements...... on the creep response may provide some useful information about how to improve the creep resistance of magnesium alloys in the future. (c) 2008 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved....

Inelastic finite element analysis of a pipe-elbow assembly (benchmark problem 2)

Energy Technology Data Exchange (ETDEWEB)

Knapp, H P [Internationale Atomreaktorbau GmbH (INTERATOM) Bergisch Gladbach (Germany); Prij, J [Netherlands Energy Research Foundation (ECN) Petten (Netherlands)

1979-06-01

In the scope of the international benchmark problem effort on piping systems, benchmark problem 2 consisting of a pipe elbow assembly, subjected to a time dependent in-plane bending moment, was analysed using the finite element program MARC. Numerical results are presented and a comparison with experimental results is made. It is concluded that the main reason for the deviation between the calculated and measured values is due to the fact that creep-plasticity interaction is not taken into account in the analysis. (author)

Creep fracture mechanics analysis for through-wall cracked pipes under widespread creep condition

International Nuclear Information System (INIS)

Huh, Nam Su; Kim, Yun Jae; Kim, Young Jin

2003-01-01

This paper compares engineering estimation schemes of C * and creep COD for circumferential and axial through-wall cracked pipes at elevated temperatures with detailed 3-D elastic-creep finite element results. Engineering estimation schemes included the GE/EPRI method, the reference stress method where reference stress is defined based on the plastic limit load and the enhanced reference stress method where the reference stress is defined based on the optimized reference load. Systematic investigations are made not only on the effect of creep-deformation behaviour on C * and creep COD, but also on effects of the crack location, the pipe geometry, the crack length and the loading mode. Comparison of the FE results with engineering estimations provides that for idealized power law creep, estimated C * and COD rate results from the GE/EPRI method agree best with FE results. For general creep-deformation laws where either primary or tertiary creep is important and thus the GE/EPRI method is hard to apply, on the other hand, the enhanced reference stress method provides more accurate and robust estimations for C * and COD rate than the reference stress method

Microtomography and creep modeling of a short fiber reinforced aluminum piston alloy

Energy Technology Data Exchange (ETDEWEB)

Marks, Esteban; Requena, Guillermo; Degischer, Hans-Peter [Institute of Materials Science and Technology, Vienna University of Technology Karlsplatz 13/E308, A-1040 Vienna (Austria); Boller, Elodie [European Synchrotron Radiation Facility, Grenoble (France)

2011-03-15

Interconnectivity between eutectic silicon and short fibers in an AlSi12CuNiMg/Al2O3/15s composite increases with long-term creep exposure time due to diffusion. It is also observed that the stationary creep rate decreases significantly compared to the initial stationary creep rate. Interconnectivity between rigid phases is analyzed and three-dimensional statistical functions are applied to deduce the representative volume elements of the composite in different conditions. Unit cell geometrical models are generated to simulate the stationary creep behavior using the finite element method. (Copyright copyright 2011 WILEY-VCH Verlag GmbH and Co. KGaA, Weinheim)

Application of viscoplastic constitutive equations in finite element programs

International Nuclear Information System (INIS)

Hornberger, K.; Stamm, H.

1987-04-01

The general mathematical formulation of frequently used viscoplastic constitutive equations is explained and Robinson's model is discussed in more detail. The implementation of viscoplastic constitutive equations into Finite Element programs (such as ABAQUS) is described using Robinson's model as an example. For the numerical integration both an explicit (explicit Euler) and an implicit (generalized midpoint rule) integration scheme is utilized in combination with a time step control strategy. In the implicit integration scheme, convergence in solving a system of nonlinear algebraic equation is improved introducing a projection method. The efficiency of the implemented procedures is demonstrated for different homogeneous load cases as well as for creep loading and strain controlled cyclic loading of a perforated plate. (orig./HP) [de

Creep buckling: an experiment, an 'exact' solution and some simple thoughts

International Nuclear Information System (INIS)

Heller, P.; Anderson, R.G.

1986-01-01

The paper presents attempts to analyse and understand a carefully conducted creep buckling experiment. The analysis was conducted using the ABAQUS Finite Element Code coupled to a number of plausible creep laws. The results show good agreement between ABAQUS runs and experimental deflections but it is difficult to reproduce the early loads. A simple model of buckling analysis for n-power creep laws is derived as an aid to understanding the development of the deflections for non-linear creep laws. In particular, the model suggests why deflections develop so rapidly and how the creep deflection development relates to the elastic behaviour. (author)

Numerical treatment of creep crack growth

International Nuclear Information System (INIS)

Kienzler, R.; Hollstein, T.

1990-06-01

To accomplish the safety analysis and to predict the lifetime of high-termpature components with flaws, several concepts have been proposed to correlate creep-crack initiation and growth with fracture mechanics parameters. The concepts of stress-intensity factor K, reference stress σ ref , line integral C * , and others will be discussed. Among them, the C * -integral concept seems to have the widest range of applicability, if large creep zones develop and steady state creep conditions can be assumed. The numerical evaluation of C * by the virtual crack extension method is described. The methods are demonstrated by two- and three-dimensional finite element simulations including creep crack growth. As for ductile fracture experiments, plane stress and plane strain simulations are bounds to the three-dimensional simulations which agree well with corresponding experiments. (orig.)

Survey of creep properties of copper intended for nuclear waste disposal

Energy Technology Data Exchange (ETDEWEB)

Andersson-Oestling, Henrik C.M. (Swerea KIMAB AB, Stockholm (Sweden)); Sandstroem, Rolf (Materials Science and Engineering, School of Industrial Engineering and Management, Royal Inst. of Technology (KTH), Stockholm (Sweden))

2009-12-15

Creep in copper for application in canisters for nuclear waste disposal is surveyed. The importance of phosphorus doping to obtain adequate properties is demonstrated experimentally as well as explained theoretically. Creep tests results for electron beam and friction stir welds are compared. The latter type of welds has properties that are close to those of parent metal. The relation between slow strain rate tensile and creep is described. Fundamental constitutive equations are presented that are suitable for finite element modelling. These equations are used to simulate creep deformation in canisters

Creep relaxation of fuel pin bending and ovalling stresses

International Nuclear Information System (INIS)

Chan, D.P.; Jackson, R.J.

1979-06-01

Analytical methods for calculating fuel pin cladding bending and ovalling stresses due to pin bundle-duct mechanical interaction taking into account nonlinear creep are presented. Calculated results are in close agreement with finite element results by MARC-CDC program. The methods are used to investigate the effect of creep on the FTR fuel cladding bending and ovalling stresses. It is concluded that the cladding of 316 SS 20% CW and reference design has high creep rates in the FTR core region to keep the bending and ovalling stresses to low levels

Massively Parallel Finite Element Programming

KAUST Repository

Heister, Timo; Kronbichler, Martin; Bangerth, Wolfgang

2010-01-01

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

Massively Parallel Finite Element Programming

KAUST Repository

Heister, Timo

2010-01-01

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

FY17 Status Report on the Micromechanical Finite Element Modeling of Creep Fracture of Grade 91 Steel

Energy Technology Data Exchange (ETDEWEB)

Messner, M. C. [Argonne National Lab. (ANL), Argonne, IL (United States); Truster, T. J. [Univ. of Tennessee, Knoxville, TN (United States); Cochran, K. B. [DR& C Inc.; Parks, D. M. [DR& C Inc.; Sham, T. -L. [Argonne National Lab. (ANL), Argonne, IL (United States)

2017-09-01

Advanced reactors designed to operate at higher temperatures than current light water reactors require structural materials with high creep strength and creep-fatigue resistance to achieve long design lives. Grade 91 is a ferritic/martensitic steel designed for long creep life at elevated temperatures. It has been selected as a candidate material for sodium fast reactor intermediate heat exchangers and other advanced reactor structural components. This report focuses on the creep deformation and rupture life of Grade 91 steel. The time required to complete an experiment limits the availability of long-life creep data for Grade 91 and other structural materials. Design methods often extrapolate the available shorter-term experimental data to longer design lives. However, extrapolation methods tacitly assume the underlying material mechanisms causing creep for long-life/low-stress conditions are the same as the mechanisms controlling creep in the short-life/high-stress experiments. A change in mechanism for long-term creep could cause design methods based on extrapolation to be non-conservative. The goal for physically-based microstructural models is to accurately predict material response in experimentally-inaccessible regions of design space. An accurate physically-based model for creep represents all the material mechanisms that contribute to creep deformation and damage and predicts the relative influence of each mechanism, which changes with loading conditions. Ideally, the individual mechanism models adhere to the material physics and not an empirical calibration to experimental data and so the model remains predictive for a wider range of loading conditions. This report describes such a physically-based microstructural model for Grade 91 at 600° C. The model explicitly represents competing dislocation and diffusional mechanisms in both the grain bulk and grain boundaries. The model accurately recovers the available experimental creep curves at higher stresses

Finite strain transient creep of D16T alloy: identification and validation employing heterogeneous tests

Science.gov (United States)

Shutov, A. V.; Larichkin, A. Yu

2017-10-01

A cyclic creep damage model, previously proposed by the authors, is modified for a better description of the transient creep of D16T alloy observed in the finite strain range under rapidly changing stresses. The new model encompasses the concept of kinematic hardening, which allows us to account for the creep-induced anisotropy. The model kinematics is based on the nested multiplicative split of the deformation gradient, proposed by Lion. The damage evolution is accounted for by the classical Kachanov-Rabotnov approach. The material parameters are identified using experimental data on cyclic torsion of thick-walled samples with different holding times between load reversals. For the validation of the proposed material model, an additional experiment is analyzed. Although this additional test is not involved in the identification procedure, the proposed cyclic creep damage model describes it accurately.

Finite element methods in incompressible, adiabatic, and compressible flows from fundamental concepts to applications

CERN Document Server

Kawahara, Mutsuto

2016-01-01

This book focuses on the finite element method in fluid flows. It is targeted at researchers, from those just starting out up to practitioners with some experience. Part I is devoted to the beginners who are already familiar with elementary calculus. Precise concepts of the finite element method remitted in the field of analysis of fluid flow are stated, starting with spring structures, which are most suitable to show the concepts of superposition/assembling. Pipeline system and potential flow sections show the linear problem. The advection–diffusion section presents the time-dependent problem; mixed interpolation is explained using creeping flows, and elementary computer programs by FORTRAN are included. Part II provides information on recent computational methods and their applications to practical problems. Theories of Streamline-Upwind/Petrov–Galerkin (SUPG) formulation, characteristic formulation, and Arbitrary Lagrangian–Eulerian (ALE) formulation and others are presented with practical results so...

Finite element analysis of irradiation-induced dilation of the fuel subassembly duct in LMFBR

International Nuclear Information System (INIS)

Gao Fuhai; Fu Hao; Li Nan; Yang Kongli; Wang Mingzhen

2013-01-01

Background: The calculation of irradiation-induced dilation of the fuel subassembly duct in LMFBR is important for fast reactor core design.. Purpose: To investigate how to calculate the dilation by using finite element method (FEM). Methods: First, irradiation-induced creep and swelling material models are introduced. Then, a theoretical solution based on a simplified bending plate model is briefly given. Finally, a stress update scheme for the adopted material models is presented and furthermore embedded into ABAQUS user interface UMAT to conduct finite element analysis. Both solutions are compared and discussed. Results: FEM successfully predicts the duct dilation and its solution agrees well with theoretical one in small deformation. Conclusions: The proposed stress update scheme is effective, The accuracy of the theory solution declines when dilation becomes larger. The maximum stress occurs at the duct corner point, and the location has stress relaxation effect. (authors)

Finite element computational fluid mechanics

International Nuclear Information System (INIS)

Baker, A.J.

1983-01-01

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

Creep fatigue assessment for EUROFER components

Energy Technology Data Exchange (ETDEWEB)

Özkan, Furkan, E-mail: oezkan.furkan@partner.kit.edu; Aktaa, Jarir

2015-11-15

Highlights: • Design rules for creep fatigue assessment are developed to EUROFER components. • Creep fatigue assessment tool is developed in FORTRAN code with coupling MAPDL. • Durability of the HCPB-TBM design is discussed under typical fusion reactor loads. - Abstract: Creep-fatigue of test blanket module (TBM) components built from EUROFER is evaluated based on the elastic analysis approach in ASME Boiler Pressure Vessel Code (BPVC). The required allowable number of cycles design fatigue curve and stress-to-rupture curve to estimate the creep-fatigue damage are used from the literature. Local stress, strain and temperature inputs for the analysis of creep-fatigue damage are delivered by the finite element code ANSYS utilizing the Mechanical ANSYS Parametric Design Language (MAPDL). A developed external FORTRAN code used as a post processor is coupled with MAPDL. Influences of different pulse durations (hold-times) and irradiation on creep-fatigue damage for the preliminary design of the Helium Cooled Pebble Bed Test Blanket Module (HCPB-TBM) are discussed for the First Wall component of the TBM box.

Relationship between strain and central deflection in small punch creep specimens

International Nuclear Information System (INIS)

Yang Zhen; Wang Zhiwen

2003-01-01

Acquiring information about creep strain directly from small punch creep tests is difficult because the deformation behaviour of the small punch specimen is complicated. A routine is suggested in the present paper to treat this problem indirectly. Based on a finite element analysis, it is proposed that the relationship of central deflection δ to central creep strain ε c of a specimen subjected to creep can be represented approximately by the relationship of central deflection δ to central (elastic-plastic) strain ε of a specimen not subjected to creep. With this hypothesis, the δ∼ε c relation of the small punch creep specimen is obtained by resorting to a rigid-plastic membrane stretch forming model. Finally, small punch creep test results are used to evaluate creep strain and creep strain rate by taking advantage of this δ∼ε c relation

Hastelloy X fuel element creep relaxation and residual effects

International Nuclear Information System (INIS)

Castle, R.A.

1971-01-01

A worst case, seven element, asymmetric fuel, thermal environment was assumed and a creep relaxation analysis generated. The fuel element clad is .020 inch Hastelloy X. The contact load decreased from 11.6 pounds to 5.87 pounds in 100,000 hours. The residual stresses were then computed for various shutdown times. (U.S.)

A simple model for indentation creep

Science.gov (United States)

Ginder, Ryan S.; Nix, William D.; Pharr, George M.

2018-03-01

A simple model for indentation creep is developed that allows one to directly convert creep parameters measured in indentation tests to those observed in uniaxial tests through simple closed-form relationships. The model is based on the expansion of a spherical cavity in a power law creeping material modified to account for indentation loading in a manner similar to that developed by Johnson for elastic-plastic indentation (Johnson, 1970). Although only approximate in nature, the simple mathematical form of the new model makes it useful for general estimation purposes or in the development of other deformation models in which a simple closed-form expression for the indentation creep rate is desirable. Comparison to a more rigorous analysis which uses finite element simulation for numerical evaluation shows that the new model predicts uniaxial creep rates within a factor of 2.5, and usually much better than this, for materials creeping with stress exponents in the range 1 ≤ n ≤ 7. The predictive capabilities of the model are evaluated by comparing it to the more rigorous analysis and several sets of experimental data in which both the indentation and uniaxial creep behavior have been measured independently.

Homogenization of steady-state creep of porous metals using three-dimensional microstructural reconstructions

DEFF Research Database (Denmark)

Kwok, Kawai; Boccaccini, Dino; Persson, Åsa Helen

2016-01-01

The effective steady-state creep response of porous metals is studied by numerical homogenization and analytical modeling in this paper. The numerical homogenization is based on finite element models of three-dimensional microstructures directly reconstructed from tomographic images. The effects ...... model, and closely matched by the Gibson-Ashby compression and the Ramakrishnan-Arunchalam creep models. [All rights reserved Elsevier]....

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

Constitutive model of creep in polycrystalline halite based on workhardening and recovery

International Nuclear Information System (INIS)

Munson, D.E.

1993-01-01

A multimechanism constitutive model of creep has been developed which incorporates the workhardening and recovery transient creep behavior. This model has been applied to the creep of polycrystalline halite. The specific application of the model is in the calculation of the closure of underground rooms in layered salt deposits. Through the use of finite element calculations, this model, with appropriate laboratory material parameters and a Tresca flow potential, has predicted the measured closure of a number of large in situ experimental rooms

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

Computational simulation of the creep-rupture process in filamentary composite materials

Science.gov (United States)

Slattery, Kerry T.; Hackett, Robert M.

1991-01-01

A computational simulation of the internal damage accumulation which causes the creep-rupture phenomenon in filamentary composite materials is developed. The creep-rupture process involves complex interactions between several damage mechanisms. A statistically-based computational simulation using a time-differencing approach is employed to model these progressive interactions. The finite element method is used to calculate the internal stresses. The fibers are modeled as a series of bar elements which are connected transversely by matrix elements. Flaws are distributed randomly throughout the elements in the model. Load is applied, and the properties of the individual elements are updated at the end of each time step as a function of the stress history. The simulation is continued until failure occurs. Several cases, with different initial flaw dispersions, are run to establish a statistical distribution of the time-to-failure. The calculations are performed on a supercomputer. The simulation results compare favorably with the results of creep-rupture experiments conducted at the Lawrence Livermore National Laboratory.

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2001-07-01

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

Experimental verification of creep analyses for prestressed concrete reactor vessels

International Nuclear Information System (INIS)

Aoyagi, Y.; Abe, H.; Ohnuma, H.

1977-01-01

The authors proposed a new method of creep analysis based on the theory of strain hardening, which assumes that accumulated creep at a given time influences the creep after that. This method was applied to calculate step-by-step the behaviors of uniaxial creep of concrete under variable temperatures and stresses, creep in reinforced concrete specimens and the behaviors of prestressed concrete beams under themal gradients. The experimental and calculated results agreed fairly well. Further, this method was incorporated in the finite element creep analysis for the prestressed concrete hollow cylinder and the full scale model. The calculated strain changes with time pursued closely those obtained by experiments. The above led to the conclusion that from the viewpoint of both accuracy and computation time the strain hardening method proposed by the authors may be judged advantageous for practical usages

Influence of creep and cyclic oxidation in thermal barrier coatings

Energy Technology Data Exchange (ETDEWEB)

Seiler, Philipp; Baeker, Martin; Roesler, Joachim [Technische Univ. Braunschweig (Germany). Inst. fuer Werkstoffe

2012-01-15

The lifetime of thermal barrier coating systems is limited by cracks close to the interfaces, causing delamination. To study the failure mechanisms, a simplified model system is analysed which consists of a bond-coat bulk material, a thermally grown oxide, and an yttria-stabilised zirconia topcoat. The stresses in the model system are calculated using a finite element model which covers the simulation of full thermal cycles, creep in all layers, and the anisotropic oxidation during dwelling. Creep in the oxide and the thermal barrier coating is varied with the use of different creep parameter sets. The influence of creep in the bondcoat is analysed by using two different bond-coat materials: fast creeping Fecralloy and slow creeping oxide dispersion strengthened MA956. It is shown that creep in the bondcoat influences the lifetime of the coatings. Furthermore, a fast creeping thermally grown oxide benefits the lifetime of the coating system. (orig.)

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.

Design and analysis of multiaxial creep tests

International Nuclear Information System (INIS)

Mallett, R.H.; Dhalla, A.K.; Yocolano, J.T.

1974-01-01

A procedure is described for presenting the complete data as obtained from tests of thin-walled tubular creep test specimens. Thereafter, a procedure for processing the data is presented. The processed data is based in part upon results of detailed inelastic finite element analyses performed to determine uniform and constant stress quantities and effective gage lengths. (U.S.)

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

Creep and shrinkage analysis for concrete spent fuel dry storage module

Energy Technology Data Exchange (ETDEWEB)

Zhang, D. [Atomic Energy of Canada Limited, Chalk River, Ontario (Canada)], E-mail: zhangd@aecl.ca

2009-07-01

CANDU reactors are designed in Canada and are built and operated worldwide to produce electricity economically with no emission of green house gases. This paper presents creep and shrinkage analysis for a concrete spent fuel dry storage module of a CANDU nuclear power plant. Creep and shrinkage analysis was performed using a method outlined in American Concrete Institute (ACI) code, and then the creep and shrinkage strains were analyzed in a finite element model to obtain the structural behavior of the concrete module. This demonstrated that the creep and shrinkage analysis for concrete spent fuel dry storage is reasonable. AECL's spent fuel dry storage module is adequate to resist the time-dependent effects due to creep and shrinkage of concrete. (author)

Creep and shrinkage analysis for concrete spent fuel dry storage module

International Nuclear Information System (INIS)

Zhang, D.

2009-01-01

CANDU reactors are designed in Canada and are built and operated worldwide to produce electricity economically with no emission of green house gases. This paper presents creep and shrinkage analysis for a concrete spent fuel dry storage module of a CANDU nuclear power plant. Creep and shrinkage analysis was performed using a method outlined in American Concrete Institute (ACI) code, and then the creep and shrinkage strains were analyzed in a finite element model to obtain the structural behavior of the concrete module. This demonstrated that the creep and shrinkage analysis for concrete spent fuel dry storage is reasonable. AECL's spent fuel dry storage module is adequate to resist the time-dependent effects due to creep and shrinkage of concrete. (author)

Irradiation creep lifetime analysis on first wall structure materials for CFETR

Energy Technology Data Exchange (ETDEWEB)

Ma, Bing; Peng, Lei, E-mail: penglei@ustc.edu.cn; Zhang, Xiansheng; Shi, Jingyi; Zhan, Jie

2017-05-15

Fusion reactor first wall services on the conditions of high surface heat flux and intense neutron irradiation. For China Fusion Engineering Test Reactor (CFETR) with high duty time factor, it is important to analyze the irradiation effect on the creep lifetime of the main candidate structure materials for first wall, i.e. ferritic/martensitic steel, austenite steel and oxide dispersion strengthened steel. The allowable irradiation creep lifetime was evaluated with Larson-Miller Parameter (LMP) model and finite element method. The results show that the allowable irradiation creep lifetime decreases with increasing of surface heat flux, first wall thickness and inlet coolant temperature. For the current CFETR conceptual design, the lifetime is not limited by thermal creep or irradiation creep, which indicated the room for design parameters optimization.

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.

Creep characteristics in thick welded joints and their improvements. 2. Applicability of a simple model for creep analysis of thick welded joints

International Nuclear Information System (INIS)

Nakacho, Keiji; Ueda, Yukio; Kinugawa, Junichi; Yamazaki, Masayoshi

1997-01-01

Reliable predictions of the creep behavior of thick welded joints are very important to secure the safety of elevated temperature vessels like nuclear reactors. Creep behavior is very complex, thus it is difficult to perform the experiment and conduct the theoretical analysis. A simple accurate model for theoretical analysis was developed in the first report. The simple model is constructed of seven one-dimensional finite elements which can analyze not only one-dimensional stress creep behavior but also the three-dimensional situation. The simple model is verified by comparing the analyzed results with the experimental ones in this report. The model is easy to treat, and needs only a little labor and computation time to predict the creep curve and the local strain for a thick welded joint. (author)

Element-topology-independent preconditioners for parallel finite element computations

Science.gov (United States)

Park, K. C.; Alexander, Scott

1992-01-01

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

Creep buckling analysis of shells

International Nuclear Information System (INIS)

Stone, C.M.; Nickell, R.E.

1977-01-01

The current study was conducted in an effort to determine the degree of conservatism or lack of conservatism in current ASME design rules concerning time-dependent (creep) buckling. In the course of this investigation, certain observations were made concerning the numerical solution of creep buckling problems. It was demonstrated that a nonlinear finite element code could be used to solve the time-dependent buckling problem. A direct method of solution was presented which proved to be computationally efficient and provided answers which agreed very well with available analytical solutions. It was observed that the calculated buckling times could vary widely for small errors in computed displacements. The presence of high creep strain rates contributed to the prediction of early buckling times when calculated during the primary creep stage. The predicted time estimates were found to increase with time until the secondary stage was reached and the estimates approached the critical times predicted without primary creep. It can be concluded, therefore, that for most nuclear piping components, whose primary creep stage is small compared to the secondary stage, the effect of primary creep is negligible and can be omitted from the calculations. In an evaluation of the past and current ASME design rules for time-dependent, load controlled buckling, it was concluded that current use of design load safety factors is not equivalent to a safety factor of ten on service life for low creep exponents

Simplified inelastic (plastic and creep) analysis of pipe elbows subjected to inplane and out-of-plane bending

International Nuclear Information System (INIS)

Mello, R.M.; Scheller, J.D.

1975-01-01

The inelastic analysis of elbows typical of use in Liquid Metal Fast Reactor piping systems is described. Detailed information on stresses, plastic and creep strains, and deformations throughout the elbow bodies resulting from elastic/plastic, elastic/plastic/creep, and elastic/plastic-creep/relaxation material behavior was obtained for the cases of applied inplane and out-of-plane moment loading on the elbows. Some conclusions are made concerning the nature of the resulting stresses in the elbows. The simplified pipe-bend element in the MARC nonlinear finite element program is used as the analytical tool. This element permits nonlinear analysis of piping elbows and pipeline systems at realistic cost. 16 refs

Experimental and finite element analysis of bond-slip in reinforced concrete

Directory of Open Access Journals (Sweden)

A. R. V. WOLENSKI

Full Text Available Abstract The modeling of reinforced concrete structures has taken advantage of the increasing progress on Computational Mechanics, in such way that complex phenomena, such as cracking and crushing, creep, reinforcement yielding, steel-concrete bond loss, can be modeled in a reasonable realistic way, using the proper set of numerical and computational resources. Among several options, the ones based on the Finite Element Method (FEM allow complex analysis simulations of reinforced concrete structures, including the interaction of different nonlinear effects. This paper deals with the nonlinear finite element analysis of the bond-slip between reinforcing steel and concrete, taking into account an experimental study previously performed. The FEM analysis presented uses a combination of resources where the material behavior of concrete is described by the Microplane Constitutive Model, and an embedded reinforcement model is used to represent steel inside the concrete and take into account the effect of bond-slip. The FEM models were created using the INSANE (INteractive Structural ANalysis Environment computational system, open source software that has a set of FEM tools for nonlinear analysis of reinforced concrete structures. The correlations between numerical-experimentals results and several parameters validate the proposed combination of resources and identifies the significance of various effects on the response.

Probabilistic finite elements

Science.gov (United States)

Belytschko, Ted; Wing, Kam Liu

1987-01-01

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

Finite element and finite difference methods in electromagnetic scattering

CERN Document Server

Morgan, MA

2013-01-01

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

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

Creep simulation of adhesively bonded joints using modified generalized time hardening model

Energy Technology Data Exchange (ETDEWEB)

Sadigh, Mohammad Ali Saeimi [Azarbaijan Shahid Madani University, Tabriz (Iran, Islamic Republic of)

2016-04-15

Creep behavior of double lap adhesively bonded joints was investigated using experimental tests and numerical analysis. Firstly, uniaxial creep tests were carried out to obtain the creep characteristics and constitutive parameters of the adhesive at different stress and temperature levels. Generalized time hardening model was used to predict the creep behavior of the adhesive. This model was modified to simulate the creep behavior at different stress and temperature levels. Secondly, the developed model was used to simulate the creep behavior of bonded joints using finite element based numerical analysis. Creep deformations of the joints were measured experimentally and good agreement was observed in comparison with the results obtained using numerical simulation. Afterward, stress redistribution due to the creep along the adhesively bonded joint was obtained numerically. It was observed that temperature level had a significant effect on the stress redistribution along the adhesive thickness.

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

International Nuclear Information System (INIS)

Chiong, Irene; Song Chongmin

2010-01-01

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

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.

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.

Recent Advances in Creep Modelling of the Nickel Base Superalloy, Alloy 720Li.

Science.gov (United States)

Harrison, William; Whittaker, Mark; Williams, Steve

2013-03-20

Recent work in the creep field has indicated that the traditional methodologies involving power law equations are not sufficient to describe wide ranging creep behaviour. More recent approaches such as the Wilshire equations however, have shown promise in a wide range of materials, particularly in extrapolation of short term results to long term predictions. In the aerospace industry however, long term creep behaviour is not critical and more focus is required on the prediction of times to specific creep strains. The current paper illustrates the capability of the Wilshire equations to recreate full creep curves in a modern nickel superalloy. Furthermore, a finite-element model based on this method has been shown to accurately predict stress relaxation behaviour allowing more accurate component lifing.

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)

Prediction of long time creep rupture properties of welded joints using the results of short duration creep crack incubation tests

Energy Technology Data Exchange (ETDEWEB)

Hosseini, E.

2013-07-01

This dissertation submitted to the Swiss Federal Institute of Technology ETH in Zurich examines the applicability of the LICON methodology for the prediction of long-time creep rupture strength of a dissimilar metal weld. The LICON methodology is an approach for predicting the lifetime of materials under creep loading conditions. It predicts long-time uniaxial creep strength using the results from several short duration creep crack incubation tests in conjunction with the outcome of a mechanical analysis on the test-piece. This study has re-examined the previous application of the LICON methodology for 9%Cr and 1CrMoV steels. It has shown that application of the original Lion method (based on reference stress solutions) for certain materials is not appropriate. This study therefore proposes a new development for the Lion approach which uses finite-element analysis to account for the generated multiaxial stress states within welded uniaxial test-pieces.

Prediction of long time creep rupture properties of welded joints using the results of short duration creep crack incubation tests

International Nuclear Information System (INIS)

Hosseini, E.

2013-01-01

This dissertation submitted to the Swiss Federal Institute of Technology ETH in Zurich examines the applicability of the LICON methodology for the prediction of long-time creep rupture strength of a dissimilar metal weld. The LICON methodology is an approach for predicting the lifetime of materials under creep loading conditions. It predicts long-time uniaxial creep strength using the results from several short duration creep crack incubation tests in conjunction with the outcome of a mechanical analysis on the test-piece. This study has re-examined the previous application of the LICON methodology for 9%Cr and 1CrMoV steels. It has shown that application of the original Lion method (based on reference stress solutions) for certain materials is not appropriate. This study therefore proposes a new development for the Lion approach which uses finite-element analysis to account for the generated multiaxial stress states within welded uniaxial test-pieces

Creep-fatigue monitoring system for header ligaments of fossil power plants

International Nuclear Information System (INIS)

Chen, K.L.; Deardorf, A.F.; Copeland, J.F.; Pflasterer, R.; Beckerdite, G.

1993-01-01

The cracking of headers (primary and secondary superheater outlet, and reheater outlet headers) at ligament locations is an important issue for fossil power plants. A model for crack initiation and growth has been developed, based on creep-fatigue damage mechanisms. This cracking model is included in a creep-fatigue monitoring system to assess header structural integrity under high temperature operating conditions. The following principal activities are required to achieve this goal: (1) the development of transfer functions and (2) the development of a ligament cracking model. The first task is to develop stress transfer functions to convert measured (monitored) temperatures, pressures and flow rates into stresses to be used to compute damage. Elastic three-dimensional finite element analyses were performed to study transient thermal stress behavior. The sustained pressure stress redistribution due to high temperature creep was studied by nonlinear finite element analyses. The preceding results are used to derive Green's functions and pressure stress gradient transfer functions for monitoring at the juncture of the tube with the header inner surface, and for crack growth at the ligaments. The virtual crack closure method is applied to derive a stress intensity factor K solution for a corner crack at the tube/header juncture. Similarly, using the reference stress method, the steady state creep crack growth parameter C * is derived for a header corner crack. The C * solution for a small corner crack in a header can be inserted directed into the available C t solution, along with K to provide the complete transient creep solution

Finite rotation shells basic equations and finite elements for Reissner kinematics

CERN Document Server

Wisniewski, K

2010-01-01

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

A multigrid solution method for mixed hybrid finite elements

Energy Technology Data Exchange (ETDEWEB)

Schmid, W. [Universitaet Augsburg (Germany)

1996-12-31

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

Domain decomposition methods for mortar finite elements

Energy Technology Data Exchange (ETDEWEB)

Widlund, O.

1996-12-31

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

Why do probabilistic finite element analysis ?

CERN Document Server

Thacker, Ben H

2008-01-01

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

An efficient finite element solution for gear dynamics

International Nuclear Information System (INIS)

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

2010-01-01

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

Finite element analysis theory and application with ANSYS

CERN Document Server

Moaveni, Saeed

2015-01-01

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

Finite element application to global reactor analysis

International Nuclear Information System (INIS)

Schmidt, F.A.R.

1981-01-01

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

Finite-Element Software for Conceptual Design

DEFF Research Database (Denmark)

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

2010-01-01

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

Books and monographs on finite element technology

Science.gov (United States)

Noor, A. K.

1985-01-01

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

Impression creep properties of a semi-solid processed magnesium-aluminum alloy containing calcium and rare earth elements

International Nuclear Information System (INIS)

Nami, B.; Razavi, H.; Miresmaeili, S.M.; Mirdamadi, Sh.; Shabestari, S.G.

2011-01-01

The creep properties of a thixoformed magnesium-aluminum alloy containing calcium and rare earth elements were studied under shear modulus-normalized stresses ranging from 0.0225 to 0.035 at temperatures of 150-212 o C using the impression creep technique. Analysis of the creep mechanism based on a power-law equation indicated that pipe diffusion-controlled dislocation climb is the dominant mechanism during creep. The alloy has a better creep resistance than high-pressure die-cast magnesium-aluminum alloy.

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.

Creep cavitation in the neighborhood of stress concentrations

International Nuclear Information System (INIS)

Liu, T.S.; Delph, T.J.; Fields, R.J.

1983-01-01

The results of several experiments into the formation and distribution of creep cavitation in the neighborhood of stress concentrations is reported. Of particular interest is the use of an image analyzing computer to construct quantitative maps of cavity sizes and distributions. Comparisons are drawn in one case with the results of a finite element simulation, and some degree of overall agreement is noted. (orig.)

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.

Creep crack growth by grain boundary cavitation under monotonic and cyclic loading

Science.gov (United States)

Wen, Jian-Feng; Srivastava, Ankit; Benzerga, Amine; Tu, Shan-Tung; Needleman, Alan

2017-11-01

Plane strain finite deformation finite element calculations of mode I crack growth under small scale creep conditions are carried out. Attention is confined to isothermal conditions and two time histories of the applied stress intensity factor: (i) a monononic increase to a plateau value subsequently held fixed; and (ii) a cyclic time variation. The crack growth calculations are based on a micromechanics constitutive relation that couples creep deformation and damage due to grain boundary cavitation. Grain boundary cavitation, with cavity growth due to both creep and diffusion, is taken as the sole failure mechanism contributing to crack growth. The influence on the crack growth rate of loading history parameters, such as the magnitude of the applied stress intensity factor, the ratio of the applied minimum to maximum stress intensity factors, the loading rate, the hold time and the cyclic loading frequency, are explored. The crack growth rate under cyclic loading conditions is found to be greater than under monotonic creep loading with the plateau applied stress intensity factor equal to its maximum value under cyclic loading conditions. Several features of the crack growth behavior observed in creep-fatigue tests naturally emerge, for example, a Paris law type relation is obtained for cyclic loading.

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.

Ratchetting in the creep range

International Nuclear Information System (INIS)

Ponter, A.R.S.; Cocks, A.C.F.; Clement, G.; Roche, R.; Corradi, L.; Franchi, A.

1985-01-01

This report attempts to present a ''State of the Art'' of this problem from three contracting and complementary points of view which reflect separate traditions within the discipline of structural analysis. Part I gives a brief summary of the essential elements of the three constitutive parts and a set of conclusions and recommendations are then formulated. Part II is an attempt by a group at CEA Saclay, France, to distil from available experimental data a set of rules expressed in terms of the stress classifications of the ASME codes, which will ensure the prevention of excessive creep ratchetting. The resulting stresses to an effective (or reference) stress and the creep assessment is then made in terms of the creep produced by the effective stress. They aim at analytical procedures for LMFBR components that operate in the creep region and are subject to considerable thermal transients. Part III by Ponter and Cocks of the University of Leicester is a theoretical study of the problem using bounding and other approximate techniques. The problem is studied in a sequence of increasingly complex problems commencing with an isothermal structure subjected to constant load and terminating in a structure subjected to arbitrary cyclic thermal loading. The results are expressed in terms of a reference stress derived from a plastic shakedown solution, and a reference history of temperature. These techniques are capable of providing assessment of the creep deformation of a structure when the plastic shakedown properties of the structures are known. The particular circumstances which occur in a LMFBR are emphasized. Part IV by Corradi and Franchi discusses the methods by which finite element solution may be calculated. These are surveyed with particular reference to the numerical problems involved and the relationship between computational procedure and the form of the constitutive equation. 162 refs

Non-Contact Measurements of Creep Properties of Refractory Materials

Science.gov (United States)

Lee, Jonghyun; Bradshaw, Richard C.; Hyers, Robert W.; Rogers, Jan R.; Rathz, Thomas J.; Wall, James J.; Choo, Hahn; Liaw, Peter

2006-01-01

State-of-the-art technologies for hypersonic aircraft, nuclear electric/thermal propulsion for spacecraft, and more efficient jet engines are driving ever more demanding needs for high-temperature (>2000 C) materials. At such high temperatures, creep rises as one of the most important design factors to be considered. Since conventional measurement techniques for creep resistance are limited to about 17OO0C, a new technique is in demand for higher temperatures. This paper presents a non-contact method using electrostatic levitation (ESL) which is applicable to both metallic and non-metallic materials. The samples were rotated quickly enough to cause creep deformation by centrifugal acceleration. The deformation of the samples was captured with a high speed camera and then the images were analyzed to estimate creep resistance. Finite element analyses were performed and compared to the experiments to verify the new method. Results are presented for niobium and tungsten, representative refractory materials at 2300 C and 2700 C respectively.

Pretest 3D finite element analysis of the WIPP Intermediate Scale Borehole Test

International Nuclear Information System (INIS)

Arguello, J.G.

1991-11-01

A three dimensional pretest finite element analysis of the Intermediate Scale Borehole Test has been performed. In the analysis, the 7.7 years simulation period includes the mining of Rooms C1 and C2, and the N1420 cross drift, at time zero; drilling of the borehole between the two rooms at 5.7 years; and 2 years of post-drilling response. An all salt configuration was used in the calculation. The 1984 Waste Isolation Pilot Plant (WIPP) reference elastic-secondary creep law, with reduced elastic moduli, was used to model the creeping response of the salt. Results show that after mining of the rooms and cross drift a relatively high von Mises stress state exists around the perimeter of the pillar. However, by 5.7 years, or immediately prior to drilling of the borehole, the pillar has relaxed to an almost uniform von Mises stress of about 7--8 MPa. After the borehole is drilled, a relatively high von Mises stress field is once again set up in the immediate vicinity of the hole. This drives the creep closure of the borehole. The hole closes more in the vertical direction than in the horizontal direction, resulting in ovalling of the hole. At the end of the simulation, the von Mises stress around the borehole is still higher than that in the remained of the pillar. Thus, the closure rates are relatively high at the end of the simulation time

Estimation of creep life of thick welded joints using a simple model. Creep characteristics in thick welded joint and their improvements. 2

International Nuclear Information System (INIS)

Nakacho, Keiji; Yamazaki, Masayoshi

2001-01-01

The information of the creep behavior of the thick welded joint is very important to secure the safety of the elevated temperature vessels like the nuclear reactors. The creep behavior of the thick welded point is very complex, thence it is difficult to practice the experiment or the theoretical analysis. A simple accurate model for theoretical analysis was developed in the first study. The simple model is constructed of several one-dimensional finite elements which can analyze three-dimensional creep behavior under a assumption. The model is easy to treat, and needs only a little labor and computation time to simulate the creep curve and local strain of the thick welded joint. In this second study, the capability of the model is expanded to estimate the creep life of the thick welded joint. New model can easily estimate the time of the rupture of the thick welded joint. It is verified comparing the result with the experimental one that the model can accurately predict the creep life. The histories of the local strains to the rupture time may be observed in the simulation by using the model. The information will be useful to improve the creep characteristics of the joints. (author)

Creep and Recovery Behaviour of Polyolefin-Rubber Nanocomposites Developed for Additive Manufacturing

Directory of Open Access Journals (Sweden)

Fugen Daver

2016-12-01

Full Text Available Nanocomposite application in automotive engineering materials is subject to continual stress fields together with recovery periods, under extremes of temperature variations. The aim is to prepare and characterize polyolefin-rubber nanocomposites developed for additive manufacturing in terms of their time-dependent deformation behaviour as revealed in creep-recovery experiments. The composites consisted of linear low density polyethylene and functionalized rubber particles. Maleic anhydride compatibilizer grafted to polyethylene was used to enhance adhesion between the polyethylene and rubber; and multi-walled carbon nanotubes were introduced to impart electrical conductivity. Various compositions of nanocomposites were tested under constant stress in creep and recovery. A four-element mechanistic Burger model was employed to model the creep phase of the composites, while a Weibull distribution function was employed to model the recovery phase of the composites. Finite element analysis using Abaqus enabled numerical modelling of the creep phase of the composites. Both analytical and numerical solutions were found to be consistent with the experimental results. Creep and recovery were dependent on: (i composite composition; (ii compatibilizers content; (iii carbon nanotubes that formed a percolation network.

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

On the derivation of a creep law from isothermal bore hole convergence

International Nuclear Information System (INIS)

Prij, J.; Mengelers, J.H.J.

1981-01-01

Some analytical as well as numerical aspects relevant to the creep behaviour of cavity-like structures in salt domes are presented. Two finite element models are presented for the modelling of the bore hole configuration, both dealing with the problem of a correct choice of the amount of salts which must be taken into account. A numerical procedure is suggested to derive a material creep law from measured bore hole convergence. This procedure is applied on convergence measurement in the ASSE mine (Germany) leading to a secondary creep law (depsilon/dt)sup(c)=8.8x10 -11 sigmasup(5.5) (sigma in MPa, (depsilon/dt)sup(c) in days -1 ) which describes the transient convergence behaviour correctly. Some questions concerning the uniqueness of the derived creep law are discussed

Evaluation of conversion relationships for impression creep test at elevated temperatures

International Nuclear Information System (INIS)

Hyde, T.H.; Sun, W.

2009-01-01

This paper contains some results related to the evaluation of the conversion relationships between impression creep test data and conventional uniaxial creep test date, for determining the secondary creep properties at elevated temperature. Some important aspects, including conversion factors, specimen dimensions, typical test results and validity of the test technique etc are briefly reviewed. The method used to determine the conversion factors is based on a reference stress approach using the results of finite element (FE) analyses; this is described in the paper. The conversion factors (reference parameters) obtained from 2-dimensional (2D) and 3-dimensional (3D) FE analyses are compared and the effects of specimen geometry, on the conversion relationships, are assessed. The recommendations on the use of these conversion factors, in practical impression creep testing, are given. Proposals for future exploitation of the technique are addressed.

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

Multiaxial Creep-Fatigue and Creep-Ratcheting Failures of Grade 91 and Haynes 230 Alloys Toward Addressing Design Issues of Gen IV Nuclear Power Plants

Energy Technology Data Exchange (ETDEWEB)

Hassan, Tasnim [North Carolina State Univ., Raleigh, NC (United States); Lissenden, Cliff [Penn State Univ., University Park, PA (United States); Carroll, Laura [Idaho National Lab. (INL), Idaho Falls, ID (United States)

2015-04-01

The proposed research will develop systematic sets of uniaxial and multiaxial experimental data at a very high temperature (850-950°C) for Alloy 617. The loading histories to be prescribed in the experiments will induce creep-fatigue and creep-ratcheting failure mechanisms. These experimental responses will be scrutinized in order to quantify the influences of temperature and creep on fatigue and ratcheting failures. A unified constitutive model (UCM) will be developed and validated against these experimental responses. The improved UCM will be incorporated into the widely used finite element commercial software packages ANSYS. The modified ANSYS will be validated so that it can be used for evaluating the very high temperature ASME-NH design-by-analysis methodology for Alloy 617 and thereby addressing the ASME-NH design code issues.

Multiaxial Creep-Fatigue and Creep-Ratcheting Failures of Grade 91 and Haynes 230 Alloys Toward Addressing Design Issues of Gen IV Nuclear Power Plants

International Nuclear Information System (INIS)

Hassan, Tasnim; Lissenden, Cliff; Carroll, Laura

2015-01-01

The proposed research will develop systematic sets of uniaxial and multiaxial experimental data at a very high temperature (850-950°C) for Alloy 617. The loading histories to be prescribed in the experiments will induce creep-fatigue and creep-ratcheting failure mechanisms. These experimental responses will be scrutinized in order to quantify the influences of temperature and creep on fatigue and ratcheting failures. A unified constitutive model (UCM) will be developed and validated against these experimental responses. The improved UCM will be incorporated into the widely used finite element commercial software packages ANSYS. The modified ANSYS will be validated so that it can be used for evaluating the very high temperature ASME-NH design-by-analysis methodology for Alloy 617 and thereby addressing the ASME-NH design code issues.

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.

Precise magnetostatic field using the finite element method

International Nuclear Information System (INIS)

Nascimento, Francisco Rogerio Teixeira do

2013-01-01

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

Creep Crack Initiation and Growth Behavior for Ni-Base Superalloys

Science.gov (United States)

Nagumo, Yoshiko; Yokobori, A. Toshimitsu, Jr.; Sugiura, Ryuji; Ozeki, Go; Matsuzaki, Takashi

The structural components which are used in high temperature gas turbines have various shapes which may cause the notch effect. Moreover, the site of stress concentration might have the heterogeneous microstructural distribution. Therefore, it is necessary to clarify the creep fracture mechanism for these materials in order to predict the life of creep fracture with high degree of accuracy. In this study, the creep crack growth tests were performed using in-situ observational testing machine with microscope to observe the creep damage formation and creep crack growth behavior. The materials used are polycrystalline Ni-base superalloy IN100 and directionally solidified Ni-base superalloy CM247LC which were developed for jet engine turbine blades and gas turbine blades in electric power plants, respectively. The microstructural observation of the test specimens was also conducted using FE-SEM/EBSD. Additionally, the analyses of two-dimensional elastic-plastic creep finite element using designed methods were conducted to understand the effect of microstructural distribution on creep damage formation. The experimental and analytical results showed that it is important to determine the creep crack initiation and early crack growth to predict the life of creep fracture and it is indicated that the highly accurate prediction of creep fracture life could be realized by measuring notch opening displacement proposed as the RNOD characteristic.

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

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

Stress Calculation of a TRISO Coated Particle Fuel by Using a Poisson's Ratio in Creep Condition

International Nuclear Information System (INIS)

Cho, Moon-Sung; Kim, Y. M.; Lee, Y. W.; Jeong, K. C.; Kim, Y. K.; Oh, S. C.; Kim, W. K.

2007-01-01

KAERI, which has been carrying out the Korean VHTR (Very High Temperature modular gas cooled Reactor) project since 2004, has been developing a performance analysis code for the TRISO coated particle fuel named COPA (COated Particle fuel Analysis). COPA predicts temperatures, stresses, a fission gas release and failure probabilities of a coated particle fuel in normal operating conditions. KAERI, on the other hand, is developing an ABAQUS based finite element(FE) model to cover the non-linear behaviors of a coated particle fuel such as cracking or debonding of the TRISO coating layers. Using the ABAQUS based FE model, verification calculations were carried out for the IAEA CRP-6 benchmark problems involving creep, swelling, and pressure. However, in this model the Poisson's ratio for elastic solution was used for creep strain calculation. In this study, an improvement is made for the ABAQUS based finite element model by using the Poisson's ratio in creep condition for the calculation of the creep strain rate. As a direct input of the coefficient in a creep condition is impossible, a user subroutine for the ABAQUS solution is prepared in FORTRAN for use in the calculations of the creep strain of the coating layers in the radial and hoop directions of the spherical fuel. This paper shows the calculation results of a TRISO coated particle fuel subject to an irradiation condition assumed as in the Miller's publication in comparison with the results obtained from the old FE model used in the CRP-6 benchmark calculations

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

Energy Technology Data Exchange (ETDEWEB)

Lee, Sang Jin; Seo, Jeong Moon

2000-08-01

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

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

International Nuclear Information System (INIS)

Lee, Sang Jin; Seo, Jeong Moon

2000-08-01

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

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.

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

GOLIA-RK, Structure Stress for Isotropic Materials with Creep and Temperature Fields

International Nuclear Information System (INIS)

Donea, J.; Giuliani, S.

1976-01-01

1 - Nature of the physical problem solved: Stress analysis of complex structures in presence of creep, dimensional changes and thermal field. Plane stress, plane strain, generalized plane strain and axisymmetric problems can be solved. The material is assumed to be either isotropic or transversely isotropic. Any laws of material behaviour can easily be incorporated by the user (see subroutines WIGNER and CLAW). 2 - Method of solution: Finite element method using triangular elements with linear local fields. The equations for the displacements are solved by Choleski's method. An algorithm is incorporated to calculate automatically the successive time steps in a creep problem. 3 - Restrictions on the complexity of the problem: Maximum number of elements is 700. Maximum number of nodal points is 400. The indexes of two adjacent nodes are not permitted to differ by more than 19

Nonlinear finite element modeling of corrugated board

Science.gov (United States)

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

1999-01-01

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

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.

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.

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

Finite element formulation of fluctuating hydrodynamics for fluids filled with rigid particles using boundary fitted meshes

Energy Technology Data Exchange (ETDEWEB)

De Corato, M., E-mail: marco.decorato@unina.it [Dipartimento di Ingegneria Chimica, dei Materiali e della Produzione Industriale, Università di Napoli Federico II, Piazzale Tecchio 80, 80125 Napoli (Italy); Slot, J.J.M., E-mail: j.j.m.slot@tue.nl [Department of Mathematics and Computer Science, Eindhoven University of Technology, PO Box 513, 5600 MB Eindhoven (Netherlands); Hütter, M., E-mail: m.huetter@tue.nl [Department of Mechanical Engineering, Eindhoven University of Technology, PO Box 513, 5600 MB Eindhoven (Netherlands); D' Avino, G., E-mail: gadavino@unina.it [Dipartimento di Ingegneria Chimica, dei Materiali e della Produzione Industriale, Università di Napoli Federico II, Piazzale Tecchio 80, 80125 Napoli (Italy); Maffettone, P.L., E-mail: pierluca.maffettone@unina.it [Dipartimento di Ingegneria Chimica, dei Materiali e della Produzione Industriale, Università di Napoli Federico II, Piazzale Tecchio 80, 80125 Napoli (Italy); Hulsen, M.A., E-mail: m.a.hulsen@tue.nl [Department of Mechanical Engineering, Eindhoven University of Technology, PO Box 513, 5600 MB Eindhoven (Netherlands)

2016-07-01

In this paper, we present a finite element implementation of fluctuating hydrodynamics with a moving boundary fitted mesh for treating the suspended particles. The thermal fluctuations are incorporated into the continuum equations using the Landau and Lifshitz approach [1]. The proposed implementation fulfills the fluctuation–dissipation theorem exactly at the discrete level. Since we restrict the equations to the creeping flow case, this takes the form of a relation between the diffusion coefficient matrix and friction matrix both at the particle and nodal level of the finite elements. Brownian motion of arbitrarily shaped particles in complex confinements can be considered within the present formulation. A multi-step time integration scheme is developed to correctly capture the drift term required in the stochastic differential equation (SDE) describing the evolution of the positions of the particles. The proposed approach is validated by simulating the Brownian motion of a sphere between two parallel plates and the motion of a spherical particle in a cylindrical cavity. The time integration algorithm and the fluctuating hydrodynamics implementation are then applied to study the diffusion and the equilibrium probability distribution of a confined circle under an external harmonic potential.

The finite element method its basis and fundamentals

CERN Document Server

Zienkiewicz, Olek C; Zhu, JZ

2013-01-01

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

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

International Nuclear Information System (INIS)

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

1998-03-01

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

Finite element formulation for a digital image correlation method

International Nuclear Information System (INIS)

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

2005-01-01

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

Application of finite-element-methods in food processing

DEFF Research Database (Denmark)

Risum, Jørgen

2004-01-01

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

Analysis of structures based on a characteristic-strain model of creep

Energy Technology Data Exchange (ETDEWEB)

Bolton, J. [Alstom Power, Newbold Road, Rugby CV21 2NH (United Kingdom)], E-mail: janjohn.bolton@virgin.net

2008-01-15

A companion paper [Bolton J. In: A characteristic-strain model for creep, ECCC/I.Mech.E. conference on creep and fracture in high-temperature components, London, September 2005] describes a creep model based on a constant 'characteristic strain' at any temperature. The present paper discusses the application of such a model, first to simple structures and then to engineering components of general form under steady loading. A basis is proposed for identifying the stress within a structure, or within the critical part of a structure, which can be considered to govern both its overall and local deformations. The concept is similar to skeletal-point stress but is more readily applied to components of any shape. The implementation of the concept of 'structural stress' is discussed in the context of finite-element creep calculations. Consideration is given to the analysis of cracked structures, where very high strains at the crack tip must be accommodated.

Analysis of structures based on a characteristic-strain model of creep

Energy Technology Data Exchange (ETDEWEB)

Bolton, J. [Alstom Power, Newbold Road, Rugby CV21 2NH (United Kingdom)], E-mail: janjohn.bolton@virgin.net

2008-01-15

A companion paper [Bolton J. In: A characteristic-strain model for creep, ECCC/I.Mech.E. conference on creep and fracture in high-temperature components, London, September 2005] describes a creep model based on a constant 'characteristic strain' at any temperature. The present paper discusses the application of such a model, first to simple structures and then to engineering components of general form under steady loading. A basis is proposed for identifying the stress within a structure, or within the critical part of a structure, which can be considered to govern both its overall and local deformations. The concept is similar to skeletal-point stress but is more readily applied to components of any shape. The implementation of the concept of 'structural stress' is discussed in the context of finite-element creep calculations. Consideration is given to the analysis of cracked structures, where very high strains at the crack tip must be accommodated.

Analysis of structures based on a characteristic-strain model of creep

International Nuclear Information System (INIS)

Bolton, J.

2008-01-01

A companion paper [Bolton J. In: A characteristic-strain model for creep, ECCC/I.Mech.E. conference on creep and fracture in high-temperature components, London, September 2005] describes a creep model based on a constant 'characteristic strain' at any temperature. The present paper discusses the application of such a model, first to simple structures and then to engineering components of general form under steady loading. A basis is proposed for identifying the stress within a structure, or within the critical part of a structure, which can be considered to govern both its overall and local deformations. The concept is similar to skeletal-point stress but is more readily applied to components of any shape. The implementation of the concept of 'structural stress' is discussed in the context of finite-element creep calculations. Consideration is given to the analysis of cracked structures, where very high strains at the crack tip must be accommodated

Improved creep strength of nickel-base superalloys by optimized γ/γ′ partitioning behavior of solid solution strengthening elements

International Nuclear Information System (INIS)

Pröbstle, M.; Neumeier, S.; Feldner, P.; Rettig, R.; Helmer, H.E.; Singer, R.F.; Göken, M.

2016-01-01

Solid solution strengthening of the γ matrix is one key factor for improving the creep strength of single crystal nickel-base superalloys at high temperatures. Therefore a strong partitioning of solid solution hardening elements to the matrix is beneficial for high temperature creep strength. Different Rhenium-free alloys which are derived from CMSX-4 are investigated. The alloys have been characterized regarding microstructure, phase compositions as well as creep strength. It is found that increasing the Titanium (Ti) as well as the Tungsten (W) content causes a stronger partitioning of the solid solution strengtheners, in particular W, to the γ phase. As a result the creep resistance is significantly improved. Based on these ideas, a Rhenium-free alloy with an optimized chemistry regarding the partitioning behavior of W is developed and validated in the present study. It shows comparable creep strength to the Rhenium containing second generation alloy CMSX-4 in the high temperature / low stress creep regime and is less prone to the formation of deleterious topologically close packed (TCP) phases. This more effective usage of solid solution strengtheners can enhance the creep properties of nickel-base superalloys while reducing the content of strategic elements like Rhenium.

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

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

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

Effects of Surface Modification and Bulk Geometry on the Biotribological Behavior of Cross-Linked Polyethylene: Wear Testing and Finite Element Analysis.

Science.gov (United States)

Watanabe, Kenichi; Kyomoto, Masayuki; Saiga, Kenichi; Taketomi, Shuji; Inui, Hiroshi; Kadono, Yuho; Takatori, Yoshio; Tanaka, Sakae; Ishihara, Kazuhiko; Moro, Toru

2015-01-01

The wear and creep deformation resistances of polymeric orthopedic bearing materials are both important for extending their longevity. In this study, we evaluated the wear and creep deformation resistances, including backside damage, of different polyethylene (PE) materials, namely, conventional PE, cross-linked PE (CLPE), and poly(2-methacryloyloxyethyl phosphorylcholine)- (PMPC-) grafted CLPE, through wear tests and finite element analysis. The gravimetric and volumetric degrees of wear of disks (3 or 6 mm in thickness) of these materials against a cobalt-chromium-molybdenum alloy pin were examined using a multidirectional pin-on-disk tester. Cross-linking and PMPC grafting decreased the gravimetric wear of the PE disks significantly. The volumetric wear at the bearing surface and the volumetric penetration in the backside of the 3-mm thick PE disk were higher than those of the 6-mm thick PE disk, regardless of the bearing material. The geometrical changes induced in the PE disks consisted of creep, because the calculated internal von Mises stress at the bearing side of all disks and that at the backside of the 3-mm thick disks exceeded their actual yield strengths. A highly hydrated bearing surface layer, formed by PMPC grafting, and a cross-linking-strengthened substrate of adequate thickness are essential for increasing the wear and creep deformation resistances.

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.

Improved methods of creep-fatigue life assessment of components

Energy Technology Data Exchange (ETDEWEB)

Scholz, Alfred; Berger, Christina [Inst. fuer Werkstoffkunde (IfW), Technische Univ. Darmstadt (Germany)

2009-07-01

The improvement of life assessment methods contributes to a reduction of efforts at design and an effective long term operation of high temperature components, reduces technical risk and increases high economical advantages. Creep-fatigue at multi-stage loading, covering cold start, warm start and hot start cycles in typical loading sequences e.g. for medium loaded power plants, was investigated here. At hold times creep and stress relaxation, respectively, lead to an acceleration of crack initiation. Creep fatigue life time can be calculated by a modified damage accumulation rule, which considers the fatigue fraction rule for fatigue damage and the life fraction rule for creep damage. Mean stress effects, internal stress and interaction effects of creep and fatigue are considered. Along with the generation of advanced creep data, fatigue data and creep fatigue data as well scatter band analyses are necessary in order to generate design curves and lower bound properties inclusive. Besides, in order to improve lifing methods the enhancement of modelling activities for deformation and life time are important. For verification purposes, complex experiments at variable creep conditions as well as at creep fatigue interaction under multi-stage loading are of interest. Generally, the development of methods to transfer uniaxial material properties to multiaxial loading situations is a current challenge. For specific design purposes, a constitutive material model is introduced which is implemented as an user subroutine for Finite Element applications due to start-up and shut-down phases of components. Identification of material parameters have been performed by Neural Networks. (orig.)

Finite Element Modelling of Seismic Liquefaction in Soils

NARCIS (Netherlands)

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

2013-01-01

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

Earthquake induced rock shear through a deposition hole when creep is considered - first model. Effect on the canister and the buffer

Energy Technology Data Exchange (ETDEWEB)

Hernelind, Jan [5T Engineering AB, Vaesteraas (Sweden)

2006-08-15

March, 2000, a study regarding 'Earthquake induced rock shear through a deposition hole' was performed. Existing fractures crossing a deposition hole may be activated and sheared by an earthquake. The effect of such a rock shear has been investigated in a project that includes both laboratory tests and finite element calculations. The buffer material in a deposition hole acts as a cushion between the canister and the rock, which reduces the effect of a rock shear substantially. Lower density of the buffer yields softer material and reduced effect on the canister. However, at the high density that is suggested for a repository the stiffness of the buffer is rather high. The stiffness is also a function of the rate of shear, which means that there may be a substantial damage on the canister at very high shear rates. The rock shear has been modeled with finite element calculations with the code ABAQUS. A three-dimensional finite element mesh of the buffer and the canister has been created and simulation of a rock shear has been performed. The rock shear has been assumed to take place perpendicular to the canister at the quarter point. The shear calculations have been driven to a total shear of 20 cm. This report summarizes the effect of considering creep in the canister for one of the previous cases. Two different creep models have been used - the first one has been suggested by K Pettersson and the second one has been suggested by R Sandstroem. Both have been implemented in the FE-code ABAQUS as a user supplied subroutine CREEP. This report summarizes results obtained by using the first model suggested by K Pettersson. As can be seen from the obtained results using the first creep model (in the following named creep{sub k}p) the effect of creep in copper doesn't affect stresses and strains in the buffer and the steel part very much. However, especially the stresses in the canister are highly affected.

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

International Nuclear Information System (INIS)

Baker, A.R.

1982-07-01

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

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)

The effect of alloying elements on the creep and impact properties of high Cr steels

International Nuclear Information System (INIS)

Kim, S. H.; Song, B. J.; Ryu, W. S.

2000-01-01

The effect of minor alloying elements on the creep and impact properties in high Cr steels has been studied. The addition of W and N in creased the creep rupture strength without the decrease of the impact toughness. During deformation, growth of lath width and agglomeration of precipitates and precipitation of Laves phase occurred. These microstructural changes made the steels soften. The degree of softening was delayed by the addition of W and N. In W added steel, the Laves phase had a important role in increasing the creep rupture strength. But the impact toughness was rapidly degraded by the addition of W after aging at 600 .deg. C for 5000 hours. So it needs to evaluate more accurately the effect of Laves phase on creep and impact properties. In N added steel, V(C,N) was precipitated in lath boundary and interior of lath. The size of the precipitates was 20-50nm. The increase of creep rupture strength in N added steel may be due to the precipitate of the V(C,N). So it needs more test to clarify the effect of N on the creep and impact properties

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

Effect of Friction Coefficient on the Small Punch Creep Behavior of AISI 316L Stainless Steel

International Nuclear Information System (INIS)

Kim, Bum-Joon; Cho, Nam-Hyuck; Kim, Moon-K; Lim, Byeong-Soo

2011-01-01

Small punch creep testing has received attention due to the convenience of using smaller specimens than those of conventional uniaxial creep tests, which enables creep testing on developing or currently operational components. However, precedent studies have shown that it is necessary to consider friction between the punch and specimen when computing uniaxial equivalent stress from a finite element model. In this study, small punch creep behaviors of AISI 316L stainless steel, which is widely used in high temperature-high pressure machineries, have been compared for the two different ceramic balls such as Si 3 N 4 and Al 2 O 3 . The optimal range of the friction coefficient is 0.4⁓0.5 at 650°C for the best fit between experimental and simulation data of AISI 316 L stainless steel. The higher the friction coefficient, the longer the creep rupture time is. Therefore, the type of ceramic ball used must be specified for standardization of small punch creep testing.

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.

Verification of Orthogrid Finite Element Modeling Techniques

Science.gov (United States)

Steeve, B. E.

1996-01-01

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

Review of Tomographic Imaging using Finite Element Method

Directory of Open Access Journals (Sweden)

Mohd Fua’ad RAHMAT

2011-12-01

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

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

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

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)

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

Effect of welding on creep damage evolution in P91B steel

Energy Technology Data Exchange (ETDEWEB)

Baral, J., E-mail: jayshree2k4@gmail.com [Metallurgical and Materials Engineering, Indian Institute of Kharagpur, WB 721302 (India); Swaminathan, J. [CSIR–National Metallurgical Laboratory, Jamshedpur 831007 (India); Chakrabarti, D.; Ghosh, R.N. [Metallurgical and Materials Engineering, Indian Institute of Kharagpur, WB 721302 (India)

2017-07-15

Study of creep behavior of base metal (without weld) and welded specimens of P91B steel over a range of temperatures (600–650 °C) and stresses (50–180 MPa) showed similar values of minimum creep-rates for both specimens at higher stress regime (>100 MPa) whilst, significantly higher creep rates in the case of welded specimens at lower stress regime. Considering that welded specimen is comprised of two distinct structural regimes, i.e. weld affected zone and base metal, a method has been proposed for estimating the material parameters describing creep behavior of those regimes. Stress–strain distribution across welded specimen predicted from finite element analysis based on material parameters revealed preferential accumulation of stress and creep strain at the interface between weld zone and base metal. This is in-line with the experimental finding that creep rupture preferentially occurs at inter-critical heat affected zone in welded specimens owing to ferrite-martensite structure with coarse Cr{sub 23}C{sub 6} particles. - Highlights: •Comparison of creep properties of welded and virgin specimens of P91B steel. •At lower stresses (<100 MPa) welded samples show higher minimum creep-rate. •Creep rupture at inter-critical heat affected zone (IC-HAZ) in welded specimens. •FEA showing accumulation of creep strain in weld/base metal interface. •Precipitate free soft ferrite matrix accumulates strain and weakens IC-HAZ.

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 particle finite element method for machining simulations

Science.gov (United States)

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

2014-07-01

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

A finite element primer for beginners the basics

CERN Document Server

Zohdi, Tarek I

2014-01-01

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

On the reliability of finite element solutions

International Nuclear Information System (INIS)

Prasad, K.S.R.K.

1975-01-01

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

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

International Nuclear Information System (INIS)

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

2011-01-01

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

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

THE PRACTICAL ANALYSIS OF FINITE ELEMENTS METHOD ERRORS

Directory of Open Access Journals (Sweden)

Natalia Bakhova

2011-03-01

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

High accuracy 3D electromagnetic finite element analysis

International Nuclear Information System (INIS)

Nelson, Eric M.

1997-01-01

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

High accuracy 3D electromagnetic finite element analysis

International Nuclear Information System (INIS)

Nelson, E.M.

1996-01-01

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

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

Directory of Open Access Journals (Sweden)

Pavel A. Akimov

2017-12-01

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

Finite element simulation and testing of ISW CFRP anchorage

DEFF Research Database (Denmark)

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

2013-01-01

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

3D Finite Element Modeling for Possible Creeping Behavior of Gas Hydrate-related Slipstream Submarine Slide, offshore Vancouver Island, Canada

Science.gov (United States)

LONG, S.; He, T.; Lan, K.; Spence, G.; Yelisetti, S.

2017-12-01

Natural gas hydrate-related submarine landslides have been identified on worldwide continental slope. Being a potential risk for marine environment and engineering projects, it has been a hot topic of hydrate research in recent decades. The study target is Slipstream submarine landslide, one of the slope failures on the frontal ridges of the Northern Cascadia accretionary margin, off Vancouver Island, Canada. The previous studies of P- & S-wave velocity structure based on OBS (Ocean Bottom Seismometer) data of SeaJade (Seafloor Earthquake Array - Japan Canada Cascadia Experiment) project indicated that there are two high concentration gas-hydrate layers within the ridge, one is at a depth of 100 mbsf (meter beneath the seafloor) with anomalous high P-wave velocities and the other is just above the prominent BSR (bottom-simulating reflector) at a depth of 265-275 mbsf. In this study we investigated the possible creeping behavior of gas hydrate layer to examine the critical instability of the ridge slope using the finite element method for self weight and additional stress (e.g., mega earthquake) conditions. The elastic and elasticoplasticity moduli of gas hydrate layer were obtained from laboratory measurements for different uniaxial pressure tests, which indicated that the sediments behave elastically for uniaxial pressures below 6 MPa, but elasticoplastically between 6-6.77 MPa. The modeled shear stress distribution indicated that the current sliding surface is more likely connected with the shallow high-velocity gas hydrate layer and sliding process related with gas hydrate starts from the toe of the slope and then progressively retreats to the place of current headwall, in a series of triangular blocks or wedges. Since the study area is in the earthquake belt, the large seismic acceleration will greatly affect the stress field and pore pressure distribution within the ridge, and the landslide is going to happen and supposedly at the shallow high-velocity gas

Ferritic Alloys with Extreme Creep Resistance via Coherent Hierarchical Precipitates.

Science.gov (United States)

Song, Gian; Sun, Zhiqian; Li, Lin; Xu, Xiandong; Rawlings, Michael; Liebscher, Christian H; Clausen, Bjørn; Poplawsky, Jonathan; Leonard, Donovan N; Huang, Shenyan; Teng, Zhenke; Liu, Chain T; Asta, Mark D; Gao, Yanfei; Dunand, David C; Ghosh, Gautam; Chen, Mingwei; Fine, Morris E; Liaw, Peter K

2015-11-09

There have been numerous efforts to develop creep-resistant materials strengthened by incoherent particles at high temperatures and stresses in response to future energy needs for steam turbines in thermal-power plants. However, the microstructural instability of the incoherent-particle-strengthened ferritic steels limits their application to temperatures below 900 K. Here, we report a novel ferritic alloy with the excellent creep resistance enhanced by coherent hierarchical precipitates, using the integrated experimental (transmission-electron microscopy/scanning-transmission-electron microscopy, in-situ neutron diffraction, and atom-probe tomography) and theoretical (crystal-plasticity finite-element modeling) approaches. This alloy is strengthened by nano-scaled L21-Ni2TiAl (Heusler phase)-based precipitates, which themselves contain coherent nano-scaled B2 zones. These coherent hierarchical precipitates are uniformly distributed within the Fe matrix. Our hierarchical structure material exhibits the superior creep resistance at 973 K in terms of the minimal creep rate, which is four orders of magnitude lower than that of conventional ferritic steels. These results provide a new alloy-design strategy using the novel concept of hierarchical precipitates and the fundamental science for developing creep-resistant ferritic alloys. The present research will broaden the applications of ferritic alloys to higher temperatures.

Hydrothermal analysis in engineering using control volume finite element method

CERN Document Server

Sheikholeslami, Mohsen

2015-01-01

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

A multiscale mortar multipoint flux mixed finite element method

KAUST Repository

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

2012-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

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

1994-12-31

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

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

Science.gov (United States)

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

2004-01-01

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

Creep-fatigue assessment of a thermina test specimen using the R5 procedure

International Nuclear Information System (INIS)

Booth, P.; Budden, P.J.; Bate, S.K.

1997-01-01

A creep-fatigue life assessment of an axisymmetric 316 stainless steel test specimen under constant mechanical and cyclic thermal shock loading using the R5 Procedure is described in this paper. This test was carried out at CEA, France, and formed part of the 'Thermina' series. Stress analysis has been carried out using both full inelastic finite element analysis and also the simplified shakedown methods, based on elastic calculation, within R5. The estimates of strain range and the stress at the start of the creep dwell have then been used with R5 to estimate creep and fatigue damage per cycle and hence to make predictions of component life. The predicted lives are compared with the lives observed in the tests. The simplified R5 estimate of life, based on development of a crack of depth 200 microns, is 260 cycles using best-estimate material properties. Experimentally, cracks of depth at least 150 microns were observed in between 526 and 650 cycles, for two similar tests. The simplified R5 route therefore leads to an estimate of life which is conservative but not unduly so on this component. Detailed cyclic inelastic analysis using the ORNL constitutive model and the ABAQUS finite element code to estimate the strain range and dwell stress led to a best estimate of 618 cycles to crack initiation using R5. (author). 16 refs, 11 figs, 4 tabs

Ratchetting and creep-fatigue evaluation for nozzle-to-cylinder intersection

International Nuclear Information System (INIS)

Barsoum, R.S.; Loomis, R.W.; Stewart, B.D.

1976-01-01

The study is part of an analytical investigation on the applicability of the simplified ratchetting and creep-fatigue rules to LMFBR component geometry. Both the detailed inelastic rules and the simplified elastic rules are applied to the results obtained from a three-dimensional finite element analysis of the nozzle-to-cylinder intersection. The results of both evaluations are compared at several locations on the surface, and an assessment of the degree of conservatism of the simplified methods is discussed

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

Transformation-Induced Creep and Creep Recovery of Shape Memory Alloy.

Science.gov (United States)

Takeda, Kohei; Tobushi, Hisaaki; Pieczyska, Elzbieta A

2012-05-22

If the shape memory alloy is subjected to the subloop loading under the stress-controlled condition, creep and creep recovery can appear based on the martensitic transformation. In the design of shape memory alloy elements, these deformation properties are important since the deflection of shape memory alloy elements can change under constant stress. The conditions for the progress of the martensitic transformation are discussed based on the kinetics of the martensitic transformation for the shape memory alloy. During loading under constant stress rate, temperature increases due to the stress-induced martensitic transformation. If stress is held constant during the martensitic transformation stage in the loading process, temperature decreases and the condition for the progress of the martensitic transformation is satisfied, resulting in the transformation-induced creep deformation. If stress is held constant during the reverse transformation stage in the unloading process, creep recovery appears due to the reverse transformation. The details for these thermomechanical properties are investigated experimentally for TiNi shape memory alloy, which is most widely used in practical applications. The volume fraction of the martensitic phase increases in proportion to an increase in creep strain.

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.

An experimental and theoretical investigation of creep buckling

International Nuclear Information System (INIS)

Ohya, H.

1977-01-01

The purpose of the present paper is to investigate creep buckling phenomena and the methods of analysis. Creep buckling experiments were performed on aluminum alloy 2024-T4 cylindrical shells having radius to thickness ratios of 16, 25, 50 and 80, in single, double and triple step axial compression at 250 0 C. It was observed that buckling occurred at one of the edges and the buckling mode depended on the radius to thickness ratio and also on the applied stress level. Thicker cylinders buckled in axisymmetric mode. Thinner ones under higher applied stress levels buckled in the asymmetric mode, whereas they under lower applied stress levels buckled in the axisymmetric mode. Creep buckling times were obtained from end shortening record of the cylinders. Experimental results were compared with theoretical values obtained by the following two methods. One is a simplified method to estimate buckling times, proposed by Gerard et al., Papirno et al. and others. The method is based on the fact that the creep buckling solutions are analogous to those of plastic buckling under a certain assumption. It was found that the bukling times could be reasonably estimated by this simplified method. The other is a finite element computer program for axisymmetric thin shells. This program is based on the incremental theory and can treat thermoelastoplastic creep analysis of axisymmetric thin shells with large deflection. Creep deformation behavior of cylindrical shells under axial compression and buckling times were calculated by the program and the effects of plasticity on buckling times were also investigated

Effects of Surface Modification and Bulk Geometry on the Biotribological Behavior of Cross-Linked Polyethylene: Wear Testing and Finite Element Analysis

Directory of Open Access Journals (Sweden)

Kenichi Watanabe

2015-01-01

phosphorylcholine- (PMPC- grafted CLPE, through wear tests and finite element analysis. The gravimetric and volumetric degrees of wear of disks (3 or 6 mm in thickness of these materials against a cobalt-chromium-molybdenum alloy pin were examined using a multidirectional pin-on-disk tester. Cross-linking and PMPC grafting decreased the gravimetric wear of the PE disks significantly. The volumetric wear at the bearing surface and the volumetric penetration in the backside of the 3-mm thick PE disk were higher than those of the 6-mm thick PE disk, regardless of the bearing material. The geometrical changes induced in the PE disks consisted of creep, because the calculated internal von Mises stress at the bearing side of all disks and that at the backside of the 3-mm thick disks exceeded their actual yield strengths. A highly hydrated bearing surface layer, formed by PMPC grafting, and a cross-linking-strengthened substrate of adequate thickness are essential for increasing the wear and creep deformation resistances.

Study of the concrete tensile creep: application for the containment vessel of the nuclear power plants (PWR)

International Nuclear Information System (INIS)

Reviron, Nanthilde

2009-01-01

The aim of this work is to study experimentally and to conduct numerical simulations on the creep of concrete subjected to tensile stresses. The main purpose is to predict the behaviour of containment vessels of nuclear power plants (PWR) in the case of decennial test or accident. In order to satisfy to these industrial needs, it is necessary to characterize the behaviour of concrete under uniaxial tension. Thus, an important experimental study of tensile creep in concrete has been performed for different loading levels (50%, 70% and 90% of the tensile strength). In these tests, load was kept constant during 3 days. Several tests were performed: measurements of elastic properties and strength (in tension and in compression), monitoring of drying, shrinkage, basic creep and drying creep strains. Moreover, compressive creep tests were also performed and showed a difference with tensile creep. Furthermore, decrease of tensile strength and failure under tensile creep for large loading levels were observed. A numerical model has been proposed and developed in Cast3m finite element code. (author)

COVE-1: a finite difference creep collapse code for oval fuel pin cladding material

International Nuclear Information System (INIS)

Mohr, C.L.

1975-03-01

COVE-1 is a time-dependent incremental creep collapse code that estimates the change in ovality of a fuel pin cladding tube. It uses a finite difference method of solving the differential equations which describe the deflection of the tube walls as a function of time. The physical problem is nonlinear, both with respect to geometry and material properties, which requires the use of an incremental, analytical, path-dependent solution. The application of this code is intended primarily for tubes manufactured from Zircaloy. Therefore, provision has been made to include some of the effects of anisotropy in the flow equations for inelastic incremental deformations. 10 references. (U.S.)

Finite Element Analysis of Circular Plate using SolidWorks

International Nuclear Information System (INIS)

Kang, Yeo Jin; Jhung, Myung Jo

2011-01-01

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

A multiscale mortar multipoint flux mixed finite element method

KAUST Repository

Wheeler, Mary Fanett

2012-02-03

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

High accuracy 3D electromagnetic finite element analysis

International Nuclear Information System (INIS)

Nelson, E.M.

1997-01-01

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

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.

TSAAS: finite-element thermal and stress analysis of plane and axisymmetric solids with orthotropic temperature-dependent material properties

Energy Technology Data Exchange (ETDEWEB)

Browning, R.V.; Anderson, C.A.

1982-02-01

The finite element method is used to determine the temperatures, displacements, stresses, and strains in axisymmetric solids with orthotropic, temperature-dependent material properties under axisymmetric thermal and mechanical loads. The mechanical loads can be surface pressures, surface shears, and nodal point forces as well as an axial or centripetal acceleration. The continuous solid is replaced by a system of ring elements with triangular or quadrilateral cross sections. Accordingly, the method is valid for solids that are composed of many different materials and that have complex geometry. Nonlinear mechanical behavior as typified by plastic, locking, or creeping materials can be approximated. Two dimensional mesh generation, plotting, and editing features allow the computer program to be readily used. In addition to a stress analysis program that is based on a modified version of the SAAS code, TSAAS can carry out a transient thermal analysis with the finite element mesh used in stress analysis. An implicit time differencing scheme allows the use of arbitrary time steps with consequent fast running times. At specified times, the program will return to SAAS for thermal stress analysis. Nonlinear thermal properties and Arrhenius reaction kinetics are also incorporated into TSAAS. Several versions of TSAAS are in use at Los Alamos, running on CDC-7600, CRAY-1 and VAX 11/780 computers. This report describes the nominal TSAAS; other versions may have some unique features.

Transient effects in friction fractal asperity creep

CERN Document Server

Goedecke, Andreas

2013-01-01

Transient friction effects determine the behavior of a wide class of mechatronic systems. Classic examples are squealing brakes, stiction in robotic arms, or stick-slip in linear drives. To properly design and understand mechatronic systems of this type, good quantitative models of transient friction effects are of primary interest. The theory developed in this book approaches this problem bottom-up, by deriving the behavior of macroscopic friction surfaces from the microscopic surface physics. The model is based on two assumptions: First, rough surfaces are inherently fractal, exhibiting roughness on a wide range of scales. Second, transient friction effects are caused by creep enlargement of the real area of contact between two bodies. This work demonstrates the results of extensive Finite Element analyses of the creep behavior of surface asperities, and proposes a generalized multi-scale area iteration for calculating the time-dependent real contact between two bodies. The toolset is then demonstrated both...

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)

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)

Creep Properties of Walikukun (Schouthenia ovata Timber Beams

Directory of Open Access Journals (Sweden)

Ali Awaludin

2016-09-01

Full Text Available This study presents an evaluation of creep constants of Walikukun (Schoutheniaovata timber beams when rheological model of four solid elements, which is obtained byassembling Kelvin and Maxwell bodies in parallel configuration, was adopted. Creep behaviorobtained by this method was further discussed and compared with creep behavior developedusing phenomenological model of the previous study. Creep data of previous study was deformationmeasurement of Walikukun beams having cross-section of 15 mm by 20 mm with a clearspan of 550 mm loaded for three weeks period under two different room conditions: with andwithout Air Conditioner. Creep behavior given by both four solid elements model and phenomenological(in this case are power functions had good agreement during the period of creepmeasurement, but they give different prediction of creep factor beyond this period. The powerfunction of phenomenological model could give a reasonable creep prediction, while for the foursolid elements model a necessary modification is required to adjust its long-term creep behavior.

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

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.

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

Energy Technology Data Exchange (ETDEWEB)

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

1999-03-01

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

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

International Nuclear Information System (INIS)

Ishiguro, Misako; Higuchi, Kenji

1983-01-01

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

Finite element approximation to the even-parity transport equation

International Nuclear Information System (INIS)

Lewis, E.E.

1981-01-01

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

Linear finite element method for one-dimensional diffusion problems

Energy Technology Data Exchange (ETDEWEB)

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

2011-07-01

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

Quadrature representation of finite element variational forms

DEFF Research Database (Denmark)

Ølgaard, Kristian Breum; Wells, Garth N.

2012-01-01

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

A phenomenological creep model for nickel-base single crystal superalloys at intermediate temperatures

Science.gov (United States)

Gao, Siwen; Wollgramm, Philip; Eggeler, Gunther; Ma, Anxin; Schreuer, Jürgen; Hartmaier, Alexander

2018-07-01

For the purpose of good reproduction and prediction of creep deformation of nickel-base single crystal superalloys at intermediate temperatures, a phenomenological creep model is developed, which accounts for the typical γ/γ‧ microstructure and the individual thermally activated elementary deformation processes in different phases. The internal stresses from γ/γ‧ lattice mismatch and deformation heterogeneity are introduced through an efficient method. The strain hardening, the Orowan stress, the softening effect due to dislocation climb along γ/γ‧ interfaces and the formation of dislocation ribbons, and the Kear–Wilsdorf-lock effect as key factors in the main flow rules are formulated properly. By taking the cube slip in \\{100\\} slip systems and \\{111\\} twinning mechanisms into account, the creep behavior for [110] and [111] loading directions are well captured. Without specific interaction and evolution of dislocations, the simulations of this model achieve a good agreement with experimental creep results and reproduce temperature, stress and crystallographic orientation dependences. It can also be used as the constitutive relation at material points in finite element calculations with complex boundary conditions in various components of superalloys to predict creep behavior and local stress distributions.

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

Life prediction of simple structures subject to cyclic primary and secondary loading resulting in creep and platicity

International Nuclear Information System (INIS)

Otter, N.R.; Jones, R.T.

1979-01-01

High temperature reactors are subject to cyclic mechanical and thermal loadings resulting from start up and shut down operations. The design must therefore guard against structural failure resulting from excessive deformation and creep-fatigue damage. Before any simplified inelastic analysis techniques can be applied, their validity needs to be examined under situations representative of the reactor. For this to be carried out it is necessary to determine the behaviour of components, initially geometrically simple, subject to loadings, cyclic primary and secondary in nature, which result in creep and plasticity. Beam-like structures have been investigated on a finite element basis with the aim of determining how cyclic plasticity, creep enhancement and plastic ratchetting vary in relationship with modified shakedown criteria, magnitude of loading and hold time. (orig.)

Model for transient creep of southeastern New Mexico rock salt

International Nuclear Information System (INIS)

Herrmann, W.; Wawersik, W.R.; Lauson, H.S.

1980-11-01

In a previous analysis, existing experimental data pertaining to creep tests on rock salt from the Salado formation of S.E. New Mexico were fitted to an exponential transient creep law. While very early time portions of creep strain histories were not fitted very well for tests at low temperatures and stresses, initial creep rates in particular generally being underestimated, the exponential creep law has the property that the transient creep strain approaches a finite limit with time, and is therefore desirable from a creep modelling point of view. In this report, an analysis of transient creep is made. It is found that exponential transient creep can be related to steady-state creep through a universal creep curve. The resultant description is convenient for creep analyses where very early time behavior is not important

Micromechanical Modeling of Grain Boundaries Damage in a Copper Alloy Under Creep

International Nuclear Information System (INIS)

Voese, Markus

2015-01-01

In order to include the processes on the scale of the grain structure into the description of the creep behaviour of polycrystalline materials, the damage development of a single grain boundary has been initially investigated in the present work. For this purpose, a special simulationmethod has been used, whose resolution procedure based on holomorphic functions. The mechanisms taken into account for the simulations include nucleation, growth by grain boundary diffusion, coalescence and shrinkage until complete sintering of grain boundary cavities. These studies have then been used to develop a simplified cavitation model, which describes the grain boundary damage by two state variables and the time-dependent development by a mechanism-oriented rate formulation. To include the influence of grain boundaries within continuum mechanical considerations of polycrystals, an interface model has been developed, that incorporates both damage according to the simplified cavitation model and grain boundary sliding in dependence of a phenomenological grain boundary viscosity. Furthermore a micromechanical model of a polycrystal has been developed that allows to include a material's grain structure into the simulation of the creep behaviour by means of finite element simulations. Thereby, the deformations of individual grains are expressed by a viscoplastic single crystal model and the grain boundaries are described by the proposed interface model. The grain structure is represented by a finite element model, in which the grain boundaries are modelled by cohesive elements. From the evaluation of experimental creep data, the micromechanical model of a polycrystal has been calibrated for a copper-antimony alloy at a temperature of 823 K. Thereby, the adjustment of the single crystal model has been carried out on the basis of creep rates of pure copper single crystal specimens. The experimental determination of grain boundary sliding and grain boundary porosity for coarse

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.

Research on dynamic creep strain and settlement prediction under the subway vibration loading.

Science.gov (United States)

Luo, Junhui; Miao, Linchang

2016-01-01

This research aims to explore the dynamic characteristics and settlement prediction of soft soil. Accordingly, the dynamic shear modulus formula considering the vibration frequency was utilized and the dynamic triaxial test conducted to verify the validity of the formula. Subsequently, the formula was applied to the dynamic creep strain function, with the factors influencing the improved dynamic creep strain curve of soft soil being analyzed. Meanwhile, the variation law of dynamic stress with sampling depth was obtained through the finite element simulation of subway foundation. Furthermore, the improved dynamic creep strain curve of soil layer was determined based on the dynamic stress. Thereafter, it could to estimate the long-term settlement under subway vibration loading by norms. The results revealed that the dynamic shear modulus formula is straightforward and practical in terms of its application to the vibration frequency. The values predicted using the improved dynamic creep strain formula closed to the experimental values, whilst the estimating settlement closed to the measured values obtained in the field test.

Time dependent voiding mechanisms in polyamide 6 submitted to high stress triaxiality: experimental characterisation and finite element modelling

Science.gov (United States)

Selles, Nathan; King, Andrew; Proudhon, Henry; Saintier, Nicolas; Laiarinandrasana, Lucien

2017-08-01

Double notched round bars made of semi-crystalline polymer polyamide 6 (PA6) were submitted to monotonic tensile and creep tests. The two notches had a root radius of 0.45 mm, which imposes a multiaxial stress state and a state of high triaxiality in the net (minimal) section of the specimens. Tests were carried out until the failure occurred from one of the notches. The other one, unbroken but deformed under steady strain rate or steady load, was inspected using the Synchrotron Radiation Computed Tomography (SRCT) technique. These 3D through thickness inspections allowed the study of microstructural evolution at the peak stress for the monotonic tensile test and at the beginning of the tertiary creep for the creep tests. Cavitation features were assessed with a micrometre resolution within the notched region. Spatial distributions of void volume fraction ( Vf) and void morphology were studied. Voiding mechanisms were similar under steady strain rates and steady loads. The maximum values of Vf were located between the axis of revolution of the specimens and the notch surface and voids were considered as flat cylinders with a circular basis perpendicular to the loading direction. A model, based on porous plasticity, was used to simulate the mechanical response of this PA6 material under high stress triaxiality. Both macroscopic behaviour (loading curves) and voiding micro-mechanisms (radial distributions of void volume fraction) were accurately predicted using finite element simulations.

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

Directory of Open Access Journals (Sweden)

Rajeev Kumar

2008-01-01

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

Finite element model updating of concrete structures based on imprecise probability

Science.gov (United States)

Biswal, S.; Ramaswamy, A.

2017-09-01

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

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

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.

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.

Structural impact of creep in tungsten monoblock divertor target at 20 MW/m2

Directory of Open Access Journals (Sweden)

Muyuan Li

2018-01-01

Full Text Available In order to increase erosion lifetime of the divertor target, in the 2nd design phase of R&D work package ‘Divertor’ for European DEMO, armor thickness of tungsten monoblock divertor target is increased from 5 mm to 8 mm. By increasing armor thickness, surface temperature increases nearly linearly, which makes effect of creep no longer negligible at slow transients of 20 MW/m2. In this work, structural impact of creep in tungsten monoblock divertor target is for the first time quantitatively analyzed with the aid of finite element method. The numerical simulations have revealed that creep results in an increase of inelastic strain accumulation. With increasing armor thickness, tensile surface stress along x-axis (the longer edge at the plasma-facing surface of tungsten monoblock reduces, while surface stress along z-axis (axial direction of the cooling tube changes from tensile to compressive. Creep will accelerate this change. With increasing grain size, creep strain accumulation at loading surface increases due to higher creep rates, while plastic strain accumulation decreases. Creep can mitigate the risk of deep cracking by reducing the driving force for crack opening, and has a positive impact for preventing the contact between the upper parts of neighboring monoblocks in high heat flux tests.

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

A simple finite element method for linear hyperbolic problems

International Nuclear Information System (INIS)

Mu, Lin; Ye, Xiu

2017-01-01

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

Finite element computation of plasma equilibria

International Nuclear Information System (INIS)

Rivier, M.

1977-01-01

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

Flow Applications of the Least Squares Finite Element Method

Science.gov (United States)

Jiang, Bo-Nan

1998-01-01

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

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

Directory of Open Access Journals (Sweden)

T.S. Viswanathan

2014-09-01

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

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

Application of the cracked pipe element to creep crack growth prediction

Energy Technology Data Exchange (ETDEWEB)

Brochard, J.; Charras, T.

1997-04-01

The modification of a computer code for leak before break analysis is very briefly described. The CASTEM2000 code was developed for ductile fracture assessment of piping systems with postulated circumferential through-wall cracks under static or dynamic loading. The modification extends the capabilities of the cracked pipe element to the determination of fracture parameters under creep conditions (C*, {phi}c and {Delta}c). The model has the advantage of evaluating significant secondary effects, such as those from thermal loading.

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

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

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

Comparison study of inelastic analyses for high temperature structure subjected to cyclic creep loading

International Nuclear Information System (INIS)

Kim, J. B.; Lee, H. Y.; Lee, J. H.

2002-01-01

It is necessary to develop a reliable numerical analysis method to simulate the plasticity and creep behavior of LMR high temperature structures. Since general purpose finite element analysis codes such as ABAQUS and ANSYS provide various models for plastic hardening and creep equation of Norton's power law, it is possible to perform the separate iscoplasticity analysis. In this study, the high temperature structural analysis program(NONSTA-VP) implementing Chaboche's unified visco plasticity equation into ABAQUS has been developed and the viscoplastic response of the 316 SS plate having a circular hole subjected to a cyclic creep loading has been analyzed. The results among the separate visco plasticity analyses and the unified visco plasticity analysis using NONSTA-VP have been compared and the results from NONSTA-VP shows remarkable responses of stress relaxation and creep behavior during hold time compared to those from separate visco plasticity analyses. Also, it is anticipated to reduce the conservatism arising from using elastic approach for creep-fatigue damage analysis since the stress range and the strain range from the unified visco plasticity analysis has been greatly reduced compared to those from separate visco plasticity analyses and elastic analysis

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

Analysis of the forever-experiments using a finite element model

International Nuclear Information System (INIS)

Altstadt, E.; Moessner, T.; Sehgal, B.R.; Willschuetz, H.G.

2000-01-01

For future nuclear power plants it is demanded that there are no consequences for the environment and the population even in the closest vicinity of the plant during and after every possible accident scenario. This includes the hypothetical scenario of a severe accident with subsequent core meltdown and formation of a melt pool in the reactor pressure vessel (RPV) lower plenum of a Light Water Reactor (LWR). One accident management strategy is to stabilize the in-vessel debris configuration in the RPV as one major barrier against uncontrolled release of heat and radio nuclides. This strategy also applies to existing plants. Therefore, an axis-symmetric Finite Element (FE) model is developed based on the multipurpose code ANSYS/Multiphysics trademark . Using the Computational Fluid Dynamics (CFD) module the temperature field within the melt pool and within the vessel wall is evaluated. The transient structural mechanical calculations are then performed applying a creep model which takes into account large temperature, stress and strain variations. Taking into account both - experimental and numerical results - provides a good opportunity to improve the simulation and the understanding of the vessel failure mechanisms. Of particular interest are (i) the time to failure and (ii) the location and mode of failure. (orig.)

Influence of friction on stress and strain distributions in small punch creep test models

Czech Academy of Sciences Publication Activity Database

Dymáček, Petr; Seitl, Stanislav; Milička, Karel; Dobeš, Ferdinand

417-418, - (2010), s. 561-564 ISSN 1013-9826. [International Conference on Fracture and Damage Mechanics /8./. Malta, 08.09.2009-10.09.2009] R&D Projects: GA AV ČR(CZ) IAA200410801; GA AV ČR(CZ) 1QS200410502 Institutional research plan: CEZ:AV0Z20410507 Keywords : small punch test * creep * chromium steel * finite element method Subject RIV: JG - Metallurgy www.scientific.net/KEM.417-418.561

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

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

Analysis of Simple Creep Stress Calculation Methods for Creep Life Assessment

Energy Technology Data Exchange (ETDEWEB)

Seo, Jun Min; Lee, Han Sang; Kim, Yun Jae [Korea Univ., Daejeon (Korea, Republic of)

2017-08-15

Creep analysis takes much more time than elastic or elastic-plastic analysis. In this study, we conducted elastic and elastic-plastic analysis and compared the results with creep analysis results. In the elastic analysis, we used primary stress, which can be classified by the Mα-tangent method and stress intensities recommended in the ASME code. In the elastic-plastic analysis, we calculated the parameters recommended in the R5 code. For the FE models, a bending load, uniaxial load, and biaxial load were applied to the cross shaped welded plate, and a bending load and internal pressure were applied to the elbow pipe. To investigate the element size sensitivity, we conducted FE analysis for various element sizes for the cases where bending load was applied to the cross shaped welded plate. There was no significant difference between the creep.

Application of Linear Viscoelastic Properties in Semianalytical Finite Element Method with Recursive Time Integration to Analyze Asphalt Pavement Structure

Directory of Open Access Journals (Sweden)

Pengfei Liu

2018-01-01

Full Text Available Traditionally, asphalt pavements are considered as linear elastic materials in finite element (FE method to save computational time for engineering design. However, asphalt mixture exhibits linear viscoelasticity at small strain and low temperature. Therefore, the results derived from the elastic analysis will inevitably lead to discrepancies from reality. Currently, several FE programs have already adopted viscoelasticity, but the high hardware demands and long execution times render them suitable primarily for research purposes. Semianalytical finite element method (SAFEM was proposed to solve the abovementioned problem. The SAFEM is a three-dimensional FE algorithm that only requires a two-dimensional mesh by incorporating the Fourier series in the third dimension, which can significantly reduce the computational time. This paper describes the development of SAFEM to capture the viscoelastic property of asphalt pavements by using a recursive formulation. The formulation is verified by comparison with the commercial FE software ABAQUS. An application example is presented for simulations of creep deformation of the asphalt pavement. The investigation shows that the SAFEM is an efficient tool for pavement engineers to fast and reliably predict asphalt pavement responses; furthermore, the SAFEM provides a flexible, robust platform for the future development in the numerical simulation of asphalt pavements.

Evaluation of Concrete Cylinder Tests Using Finite Elements

DEFF Research Database (Denmark)

Saabye Ottosen, Niels

1984-01-01

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

Finite element analysis of structures through unified formulation

CERN Document Server

Carrera, Erasmo; Petrolo, Marco; Zappino, Enrico

2014-01-01

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

Exponentially-convergent Monte Carlo via finite-element trial spaces

International Nuclear Information System (INIS)

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

2011-01-01

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

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

Revision of Drucker-Prager cap creep modelling of pebble beds in fusion blankets

International Nuclear Information System (INIS)

Hofer, D.; Kamlah, M.; Hermsmeyer, S.

2004-01-01

A continuum model commonly used in soil mechanics analysis is compiled by use of a finite element software and has been used to simulate the thermomechanical behaviour of pebble beds. The Drucker-Prager Cap theory accounts for inelastic volume change, cap hardening, nonlinear elasticity and pressure dependent shear failure. The hardening mechanism allows for defining the hydrostatic pressure yield stress as a function of the volumetric inelastic strain. Volumetric creep is considered in order to simulate the pebble bed behaviour at high temperatures. Here, the strain hardening option has been used for the consolidation creep mechanism. The model has been calibrated using the fitting curves of the oedometric test given by Reimann et al. The fitted data has been used to calculate a pebble bed with simplified boundary conditions loaded by non-uniform volumetric heating. This calculation demonstrated that the model is capable of representing creep behaviour under volumetric heating conditions. (author)

Drucker-Prager-Cap creep modelling of pebble beds in fusion blankets

International Nuclear Information System (INIS)

Hofer, D.; Kamlah, M.

2005-01-01

Modelling of thermal and mechanical behaviour of pebble beds for fusion blankets is an important issue to understand the interaction of solid breeder and beryllium pebble beds with the surrounding structural material. Especially the differing coefficients of thermal expansion of these materials cause high stresses and strains during irradiation induced volumetric heating. To describe this process, the coupled thermomechanical behaviour of both pebble bed materials has to be modelled. Additionally, creep has to be considered contributing to bed deformations and stress relaxation. Motivated by experiments, we use a continuum mechanical approach called Drucker-Prager/Cap theory to model the macroscopic pebble bed behaviour. The model accounts for pressure dependent shear failure, inelastic hardening, and volumetric creep. The elastic part is described by a nonlinear elasticity law. The model has been implemented by user-defined routines in the commercial finite-element code ABAQUS. To check the numerics, the implementation is compared to an analytical solution. Furthermore, the Drucker-Prager/Cap tool is applied to a single ceramic breeder bed subject to creep under volumetric heating

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

Multi Resolution In-Situ Testing and Multiscale Simulation for Creep Fatigue Damage Analysis of Alloy 617

Energy Technology Data Exchange (ETDEWEB)

Liu, Yongming [Arizona State Univ., Tempe, AZ (United States). School for Engineering of Matter, Transport and Energy; Oskay, Caglar [Vanderbilt Univ., Nashville, TN (United States). Dept. of Civil and Environmental Engineering

2017-04-30

This report outlines the research activities that were carried out for the integrated experimental and simulation investigation of creep-fatigue damage mechanism and life prediction of Nickel-based alloy, Inconel 617 at high temperatures (950° and 850°). First, a novel experimental design using a hybrid control technique is proposed. The newly developed experimental technique can generate different combinations of creep and fatigue damage by changing the experimental design parameters. Next, detailed imaging analysis and statistical data analysis are performed to quantify the failure mechanisms of the creep fatigue of alloy 617 at high temperatures. It is observed that the creep damage is directly associated with the internal voids at the grain boundaries and the fatigue damage is directly related to the surface cracking. It is also observed that the classical time fraction approach does not has a good correlation with the experimental observed damage features. An effective time fraction parameter is seen to have an excellent correlation with the material microstructural damage. Thus, a new empirical damage interaction diagram is proposed based on the experimental observations. Following this, a macro level viscoplastic model coupled with damage is developed to simulate the stress/strain response under creep fatigue loadings. A damage rate function based on the hysteresis energy and creep energy is proposed to capture the softening behavior of the material and a good correlation with life prediction and material hysteresis behavior is observed. The simulation work is extended to include the microstructural heterogeneity. A crystal plasticity finite element model considering isothermal and large deformation conditions at the microstructural scale has been developed for fatigue, creep-fatigue as well as creep deformation and rupture at high temperature. The model considers collective dislocation glide and climb of the grains and progressive damage accumulation of

Creep analysis by the path function method

International Nuclear Information System (INIS)

Akin, J.E.; Pardue, R.M.

1977-01-01

The finite element method has become a common analysis procedure for the creep analysis of structures. The most recent programs are designed to handle a general class of material properties and are able to calculate elastic, plastic, and creep components of strain under general loading histories. The constant stress approach is too crude a model to accurately represent the actual behaviour of the stress for large time steps. The true path of a point in the effective stress-effective strain (sigmasup(e)-epsilonsup(c)) plane is often one in which the slope is rapidly changing. Thus the stress level quickly moves away from the initial stress level and then gradually approaches the final one. The result is that the assumed constant stress level quickly becomes inaccurate. What is required is a better method of approximation of the true path in the sigmasup(e)-epsilonsup(c) space. The method described here is called the path function approach because it employs an assumed function to estimate the motion of points in the sigmasup(e)-epsilonsup(c) space. (Auth.)

A finite element for plates and shells

International Nuclear Information System (INIS)

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

1981-08-01

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

COMPUTER EXPERIMENTS WITH FINITE ELEMENTS OF HIGHER ORDER

Directory of Open Access Journals (Sweden)

Khomchenko A.

2017-12-01

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

Fourier analysis of finite element preconditioned collocation schemes

Science.gov (United States)

Deville, Michel O.; Mund, Ernest H.

1990-01-01

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

Stochastic Finite Elements in Reliability-Based Structural Optimization

DEFF Research Database (Denmark)

Sørensen, John Dalsgaard; Engelund, S.

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

Modelling bucket excavation by finite element

Science.gov (United States)

Pecingina, O. M.

2015-11-01

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

Introduction to the Explicit Finite Element Method for Nonlinear Transient Dynamics

CERN Document Server

Wu, Shen R

2012-01-01

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

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

International Nuclear Information System (INIS)

Sung, Jin Il; Yoo, Jeong Hoon

2002-01-01

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

New evaluation method of crack growth in SiC/SiC composites using interface elements

International Nuclear Information System (INIS)

Serizawa, H.; Ando, M.; Lewinsohn, C.A.; Murakawa, H.

2000-01-01

Crack propagation behavior in SiC/SiC composites was analyzed using a new computer simulation method that included time-dependent interface elements. The simulation method was used to describe the time-dependent crack growth in SiC/SiC composites under four-point bending of single-edge-notched beam bend-bars. Two methods were used to simulate time-dependent crack growth in SiC/SiC composites due to fiber creep. In one method, the creep property was introduced into the interface elements by the general method of finite element method (FEM) analysis. In the second method, a new technique making the best use of the potential function was used to represent crack closure tractions due to creeping fibers. The stage-II slow crack growth of a general creep deformation was simulated by both methods. Additionally, stage-III crack growth and the transition from stage-II to stage-III could be simulated by the new method. The new method has the potential to completely simulate time-dependent crack growth behavior in SiC/SiC composites due to fiber creep

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.

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

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)

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)

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

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

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.

Creep behavior under internal pressure of zirconium alloy cladding oxidized in steam at high temperature

International Nuclear Information System (INIS)

Chosson, Raphael

2014-01-01

During hypothetical Loss-Of-Coolant-Accident (LOCA) scenarios, zirconium alloy fuel cladding tubes creep under internal pressure and are oxidized on their outer surface at high temperature (HT). Claddings become stratified materials: zirconia and oxygen-stabilized α phase, called α(O), are formed on the outer surface of the cladding whereas the inner part remains in the β domain. The strengthening effect of oxidation on the cladding creep behavior under internal pressure has been highlighted at HT. In order to model this effect, the creep behavior of each layer had to be determined. This study focused on the characterization of the creep behavior of the α(O) phase at HT, through axial creep tests performed under vacuum on model materials, containing from 2 to 7 wt.% of oxygen and representative of the α(O) phase. For the first time, two creep flow regimes have been observed in this phase. Underlying physical mechanisms and relevant microstructural parameters have been discussed for each regime. The strengthening effect due to oxygen on the α(O) phase creep behavior at HT has been quantified and creep flow equations have been identified. A ductile to brittle transition criterion has been also suggested as a function of temperature and oxygen content. Relevance of the creep flow equations for each layer, identified in this study or from the literature, has been discussed. Then, a finite element model, describing the oxidized cladding as a stratified material, has been built. Based on this model, a fraction of the experimental strengthening during creep is predicted. (author) [fr

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)

Three-dimensional finite element models of the human pubic symphysis with viscohyperelastic soft tissues.

Science.gov (United States)

Li, Zuoping; Alonso, Jorge E; Kim, Jong-Eun; Davidson, James S; Etheridge, Brandon S; Eberhardt, Alan W

2006-09-01

Three-dimensional finite element (FE) models of human pubic symphyses were constructed from computed tomography image data of one male and one female cadaver pelvis. The pubic bones, interpubic fibrocartilaginous disc and four pubic ligaments were segmented semi-automatically and meshed with hexahedral elements using automatic mesh generation schemes. A two-term viscoelastic Prony series, determined by curve fitting results of compressive creep experiments, was used to model the rate-dependent effects of the interpubic disc and the pubic ligaments. Three-parameter Mooney-Rivlin material coefficients were calculated for the discs using a heuristic FE approach based on average experimental joint compression data. Similarly, a transversely isotropic hyperelastic material model was applied to the ligaments to capture average tensile responses. Linear elastic isotropic properties were assigned to bone. The applicability of the resulting models was tested in bending simulations in four directions and in tensile tests of varying load rates. The model-predicted results correlated reasonably with the joint bending stiffnesses and rate-dependent tensile responses measured in experiments, supporting the validity of the estimated material coefficients and overall modeling approach. This study represents an important and necessary step in the eventual development of biofidelic pelvis models to investigate symphysis response under high-energy impact conditions, such as motor vehicle collisions.

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

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.

Finite Element Simulation of Fracture Toughness Test

International Nuclear Information System (INIS)

Chu, Seok Jae; Liu, Cong Hao

2013-01-01

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

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

International Nuclear Information System (INIS)

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

2002-01-01

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

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

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

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

Rupture of Al matrix in U-Mo/Al dispersion fuel by fission induced creep

Energy Technology Data Exchange (ETDEWEB)

Jeong, Gwan Yoon; Sohn, Dong Seong [UNIST, Daejeon (Korea, Republic of); Kim, Yeon Soo [Argonne National Laboratory, Argonnge (United States); Lee, Kyu Hong [KAERI, Daejeon (Korea, Republic of)

2016-05-15

This phenomenon was found specifically in the dispersion fuel plate with Si addition in the Al matrix to suppress interaction layer (IL) formation between UMo and Al. It is known that the stresses induced by fission induced swelling in U-Mo fuel particles are relieved by creep deformation of the IL, surrounding the fuel particles, that has a much higher creep rate than the Al matrix. Thus, when IL growth is suppressed, the stress is instead exerted on the Al matrix. The observed rupture in the Al matrix is believed to be caused when the stress exceeded the rupture strength of the Al matrix. In this study, the possibility of creep rupture of the Al matrix between the neighboring U-Mo fuel particles was examined using the ABAQUS finite element analysis (FEA) tool. The predicted rupture time for a plate was much shorter than its irradiation life indicating a rupture during the irradiation. The higher stress leads Al matrix to early creep rupture in this plate for which the Al matrix with lower creep strain rate does not effectively relieve the stress caused by the swelling of the U-Mo fuel particles. For the other plate, no rupture was predicted for the given irradiation condition. The effect of creeping of the continuous phase on the state of stress is significant.

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.

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

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.

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

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.

Improved Creep Measurements for Ultra-High Temperature Materials

Science.gov (United States)

Hyers, Robert W.; Ye, X.; Rogers, Jan R.

2010-01-01

Our team has developed a novel approach to measuring creep at extremely high temperatures using electrostatic levitation (ESL). This method has been demonstrated on niobium up to 2300 C, while ESL has melted tungsten (3400 C). This method has been extended to lower temperatures and higher stresses and applied to new materials, including a niobium-based superalloy, MASC. High-precision machined spheres of the sample are levitated in the NASA MSFC ESL, a national user facility and heated with a laser. The samples are rotated with an induction motor at up to 30,000 revolutions per second. The rapid rotation loads the sample through centripetal acceleration, producing a shear stress of about 60 MPa at the center, causing the sample to deform. The deformation of the sample is captured on high-speed video, which is analyzed by machine-vision software from the University of Massachusetts. The deformations are compared to finite element models to determine the constitutive constants in the creep relation. Furthermore, the non-contact method exploits stress gradients within the sample to determine the stress exponent in a single test.

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

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

Creep model of unsaturated sliding zone soils and long-term deformation analysis of landslides

Science.gov (United States)

Zou, Liangchao; Wang, Shimei; Zhang, Yeming

2015-04-01

Sliding zone soil is a special soil layer formed in the development of a landslide. Its creep behavior plays a significant role in long-term deformation of landslides. Due to rainfall infiltration and reservoir water level fluctuation, the soils in the slide zone are often in unsaturated state. Therefore, the investigation of creep behaviors of the unsaturated sliding zone soils is of great importance for understanding the mechanism of the long-term deformation of a landslide in reservoir areas. In this study, the full-process creep curves of the unsaturated soils in the sliding zone in different net confining pressure, matric suctions and stress levels were obtained from a large number of laboratory triaxial creep tests. A nonlinear creep model for unsaturated soils and its three-dimensional form was then deduced based on the component model theory and unsaturated soil mechanics. This creep model was validated with laboratory creep data. The results show that this creep model can effectively and accurately describe the nonlinear creep behaviors of the unsaturated sliding zone soils. In order to apply this creep model to predict the long-term deformation process of landslides, a numerical model for simulating the coupled seepage and creep deformation of unsaturated sliding zone soils was developed based on this creep model through the finite element method (FEM). By using this numerical model, we simulated the deformation process of the Shuping landslide located in the Three Gorges reservoir area, under the cycling reservoir water level fluctuation during one year. The simulation results of creep displacement were then compared with the field deformation monitoring data, showing a good agreement in trend. The results show that the creeping deformations of landslides have strong connections with the changes of reservoir water level. The creep model of unsaturated sliding zone soils and the findings obtained by numerical simulations in this study are conducive to

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

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.

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.

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.

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.

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

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

Stress Analysis of a TRISO Coated Particle Fuel by Using ABAQUS Finite Element Visco-Elastoplastic Solutions

International Nuclear Information System (INIS)

Cho, Moon Sung; Kim, Y. M.; Lee, Y. W.

2006-01-01

The fundamental design for a gas-cooled reactor relies on an understanding of the behavior of a coated particle fuel. KAERI, which has been carrying out the Korean VHTR (Very High Temperature modular gas cooled Reactor) Project since 2004, is developing a fuel performance analysis code for a VHTR named COPA (COated Particle fuel Analysis). A validation of COPA was attempted by comparing its benchmark results with the visco-elastic solutions obtained from the ABAQUS code calculations for the IAEA-CRP-6 TRISO coated particle benchmark problems involving a creep, swelling, and pressure. However, the ABAQUS finite element model used for the above-mentioned analysis did not consider the material nonlinearity of the SiC coating layer that showed stress levels higher than the assumed yield point of the material. In this study, a consideration of the material nonlinearity is included in the ABAQUS model to obtain the visco-elastoplastic solutions and the results are compared with the visco-elastic solutions obtained from the previous ABAQUS model

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.

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

Modelling of cladding creep collapse

International Nuclear Information System (INIS)

Koundy, V.; Forgeron, T.; Hivroz, J.

1993-01-01

The effects of the initial ovality and pressure level on the collapse time of Zircaloy-4 tubing subjected to uniform external pressure were examined experimentally and analytically. Experiments were performed on end closed tubes with two metallurgical states: stress relieved and recrystallized. Numerical simulations were accomplished with a specific computer program based on an analytical approach and the calculated results were compared with the experimental ones. As a comparison, the finite element method is also partially examined in this analysis. Numerical collapse times are in good agreement with regard to experimental results in the case of stress relieved structure. They seem to be too conservative in the case of a recrystallized metallurgical state and the use of the anisotropic option ameliorates numerical results. Sensibility of numerical solutions to the formulation of primary creep laws are presented

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.

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

Performance and scalability of finite-difference and finite-element wave-propagation modeling on Intel's Xeon Phi

NARCIS (Netherlands)

Zhebel, E.; Minisini, S.; Kononov, A.; Mulder, W.A.

2013-01-01

With the rapid developments in parallel compute architectures, algorithms for seismic modeling and imaging need to be reconsidered in terms of parallelization. The aim of this paper is to compare scalability of seismic modeling algorithms: finite differences, continuous mass-lumped finite elements

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.

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

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

Prediction and Monitoring Systems of Creep-Fracture Behavior of 9Cr-1Mo Steels for Reactor Pressure Vessels

International Nuclear Information System (INIS)

Potirniche, Gabriel; Barlow, Fred D.; Charit, Indrajit; Rink, Karl

2013-01-01

A recent workshop on next-generation nuclear plant (NGNP) topics underscored the need for research studies on the creep fracture behavior of two materials under consideration for reactor pressure vessel (RPV) applications: 9Cr-1Mo and SA-5XX steels. This research project will provide a fundamental understanding of creep fracture behavior of modified 9Cr-1Mo steel welds for through modeling and experimentation and will recommend a design for an RPV structural health monitoring system. Following are the specific objectives of this research project: Characterize metallurgical degradation in welded modified 9Cr-1Mo steel resulting from aging processes and creep service conditions; Perform creep tests and characterize the mechanisms of creep fracture process; Quantify how the microstructure degradation controls the creep strength of welded steel specimens; Perform finite element (FE) simulations using polycrystal plasticity to understand how grain texture affects the creep fracture properties of welds; Develop a microstructure-based creep fracture model to estimate RPVs service life; Manufacture small, prototypic, cylindrical pressure vessels, subject them to degradation by aging, and measure their leak rates; Simulate damage evolution in creep specimens by FE analyses; Develop a model that correlates gas leak rates from welded pressure vessels with the amount of microstructural damage; Perform large-scale FE simulations with a realistic microstructure to evaluate RPV performance at elevated temperatures and creep strength; Develop a fracture model for the structural integrity of RPVs subjected to creep loads; and Develop a plan for a non-destructive structural health monitoring technique and damage detection device for RPVs.

Prediction and Monitoring Systems of Creep-Fracture Behavior of 9Cr-1Mo Steels for Teactor Pressure Vessels

Energy Technology Data Exchange (ETDEWEB)

Potirniche, Gabriel [Univ. of Idaho, Moscow, ID (United States); Barlow, Fred D. [Univ. of Idaho, Moscow, ID (United States); Charit, Indrajit [Univ. of Idaho, Moscow, ID (United States); Rink, Karl [Univ. of Idaho, Moscow, ID (United States)

2013-11-26

A recent workshop on next-generation nuclear plant (NGNP) topics underscored the need for research studies on the creep fracture behavior of two materials under consideration for reactor pressure vessel (RPV) applications: 9Cr-1Mo and SA-5XX steels. This research project will provide a fundamental understanding of creep fracture behavior of modified 9Cr-1Mo steel welds for through modeling and experimentation and will recommend a design for an RPV structural health monitoring system. Following are the specific objectives of this research project: Characterize metallurgical degradation in welded modified 9Cr-1Mo steel resulting from aging processes and creep service conditions; Perform creep tests and characterize the mechanisms of creep fracture process; Quantify how the microstructure degradation controls the creep strength of welded steel specimens; Perform finite element (FE) simulations using polycrystal plasticity to understand how grain texture affects the creep fracture properties of welds; Develop a microstructure-based creep fracture model to estimate RPVs service life; Manufacture small, prototypic, cylindrical pressure vessels, subject them to degradation by aging, and measure their leak rates; Simulate damage evolution in creep specimens by FE analyses; Develop a model that correlates gas leak rates from welded pressure vessels with the amount of microstructural damage; Perform large-scale FE simulations with a realistic microstructure to evaluate RPV performance at elevated temperatures and creep strength; Develop a fracture model for the structural integrity of RPVs subjected to creep loads; and Develop a plan for a non-destructive structural health monitoring technique and damage detection device for RPVs.

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

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

Implications of Jeffreys-Lomnitz Transient Creep

Science.gov (United States)

Strick, Ellis

1984-01-01

In 1958 Jeffreys proposed a power law generalization of the logarithmic transient creep earlier attributed to Lomnitz. Although Jeffreys' power law form was admittedly defective in that it became unbounded at infinite time, he did apply it to the viscoelastic behavior of the earth-moon system. Since then it has been successfully applied by many investigators to mantle rehology and Chandler wobble. Experimental seismic studies indicate that most rock types exhibit the almost constant Q behavior which Lomnitz showed to be associated with his logarithmic creep. In this paper, we study not only the Q behavior related to Jeffreys' power law creep but also other mechanical properties such as a precise spring-dashpot ladder network realization are developed. In addition, a very simple physically realizable modification of this ladder network leads to a boundedness at long times of Jeffreys' creep in a manner which does not affect his successful application at finite times.

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 quantitative prediction model of SCC rate for nuclear structure materials in high temperature water based on crack tip creep strain rate

International Nuclear Information System (INIS)

Yang, F.Q.; Xue, H.; Zhao, L.Y.; Fang, X.R.

2014-01-01

Highlights: • Creep is considered to be the primary mechanical factor of crack tip film degradation. • The prediction model of SCC rate is based on crack tip creep strain rate. • The SCC rate calculated at the secondary stage of creep is recommended. • The effect of stress intensity factor on SCC growth rate is discussed. - Abstract: The quantitative prediction of stress corrosion cracking (SCC) of structure materials is essential in safety assessment of nuclear power plants. A new quantitative prediction model is proposed by combining the Ford–Andresen model, a crack tip creep model and an elastic–plastic finite element method. The creep at the crack tip is considered to be the primary mechanical factor of protective film degradation, and the creep strain rate at the crack tip is suggested as primary mechanical factor in predicting the SCC rate. The SCC rates at secondary stage of creep are recommended when using the approach introduced in this study to predict the SCC rates of materials in high temperature water. The proposed approach can be used to understand the SCC crack growth in structural materials of light water reactors

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

High temperature creep behavior in the (α + β) phase temperature range of M5 alloy

International Nuclear Information System (INIS)

Trego, G.

2011-01-01

The isothermal steady-state creep behavior of a M5 thin sheet alloy in a vacuum environment was investigated in the (α + β) temperature, low-stress (1-10 MPa) range. To this aim, the simplest approach consists in identifying α and β creep flow rules in their respective single-phase temperature ranges and extrapolating them in the two-phase domain. However, the (α + β) experimental behavior may fall outside any bounds calculated using such creep flow data. Here, the model was improved for each phase by considering two microstructural effects: (i) Grain size: Thermo-mechanical treatments applied on the material yielded various controlled grain size distributions. Creep tests in near-α and near-β ranges evidenced a strong grain-size effect, especially in the diffusional creep regime. (ii) Chemical contrast between the two phases in the (α + β) range: From thermodynamic calculations and microstructural investigations, the β phase is enriched in Nb and depleted in O (the reverse being true for the α phase). Thus, creep tests were performed on model Zr-Nb-O thin sheets with Nb and O concentrations representative of each phase in the considered temperature range. New α and β creep flow equations were developed from this extended experimental database and used to compute, via a finite element model, the creep rates of the two-phase material. The 3D morphology of phases (β grains nucleated at α grain boundaries) was explicitly introduced in the computations. The effect of phase morphology on the macroscopic creep flow was shown using this specific morphology, compared to other typical morphologies and to experimental data. (author) [fr

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.

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.

Creep lifing methodologies applied to a single crystal superalloy by use of small scale test techniques

International Nuclear Information System (INIS)

Jeffs, S.P.; Lancaster, R.J.; Garcia, T.E.

2015-01-01

In recent years, advances in creep data interpretation have been achieved either by modified Monkman–Grant relationships or through the more contemporary Wilshire equations, which offer the opportunity of predicting long term behaviour extrapolated from short term results. Long term lifing techniques prove extremely useful in creep dominated applications, such as in the power generation industry and in particular nuclear where large static loads are applied, equally a reduction in lead time for new alloy implementation within the industry is critical. The latter requirement brings about the utilisation of the small punch (SP) creep test, a widely recognised approach for obtaining useful mechanical property information from limited material volumes, as is typically the case with novel alloy development and for any in-situ mechanical testing that may be required. The ability to correlate SP creep results with uniaxial data is vital when considering the benefits of the technique. As such an equation has been developed, known as the k SP method, which has been proven to be an effective tool across several material systems. The current work now explores the application of the aforementioned empirical approaches to correlate small punch creep data obtained on a single crystal superalloy over a range of elevated temperatures. Finite element modelling through ABAQUS software based on the uniaxial creep data has also been implemented to characterise the SP deformation and help corroborate the experimental results

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.

Creep properties of phosphorus alloyed oxygen free copper under multiaxial stress state

International Nuclear Information System (INIS)

Rui Wu; Sandstroem, Rolf; Seitisleam, Facredin

2009-10-01

Phosphorus alloyed oxygen free copper (Cu-OFP) canisters are planned to be used for spent nuclear fuel in Sweden. The copper canisters will be subjected to creep under multiaxial stress states in the repository. Creep tests have therefore been carried out at 75 deg C using double notch specimens with notch acuities of 0.5, 2, 5, and 18.8, respectively. The creep lifetime for notched specimens is considerably longer than that for the smooth one at a given net section stress, indicating that the investigated Cu-OFP is notch insensitive (notch strengthening). The notch strengthening factor in time is, for instance, greater than 70 at 180 MPa for the bluntest notch (notch acuity = 0.5). The creep lifetime is notch acuity dependent. The sharper the notch, the longer the creep lifetime is. The creep deformation is to a significant extent concentrated to the region around the notches. Different deformation on the two notches is observed. Both axial and radial strains on the failed notch are several times larger than those on the unbroken one. Linear relation between the axial and the radial strains on the notches is found. Transgranular failure is predominant, independent of stress, rupture time, and notch acuity. Adjacent to fracture, elongated grains along the stress direction, separate pores and cavities are often visible. On the unbroken notch, fewer separate cavities and cracks are only seen intergranularly for the sharper notches (notch acuity > 2). To interpret the tests for the notched creep specimens, finite element computations have been performed. A fundamental model for primary and secondary creep without fitting parameters has been used as constitutive equation. The FEM-modelling could represent the creep strain versus time curves for the notched specimens in a satisfactory way. In these curves the strain on loading is included. From the FEM-computations a stationary creep stress could be assessed, which is close to the reference stress. For a given

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

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.

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.

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.

Examination of the creep behaviour of ceramic fuel elements under neutron irradiation

International Nuclear Information System (INIS)

Brucklacher, D.

1978-01-01

This paper examines the creeping of UO 2 , UO 2 -PuO 2 and UN under neutron irradiation. It starts with the experimental results about the relation between the thermal creep rate and the load, the temperature, as well as characteristic material values, stoichiometry, grain size and porosity. These correlation are first qualitatively discussed and then compared with the statements of actual quantitative equations. From the models and theories on which these equations are based a modified Nabarro-Heering-equation results for the correlation between the creep rate of ceramic fuels, stress, temperature and the fission rate. In the experimental part of the examination, length-changes of creep samples of UO 2 , (U,Pu)O 2 and UN were measured in specially developed irradiation creep casings in different reactors. The measuring data were corrected and evaluated considering the thermal expansion effects, irregular temperature distribution and swelling effects in such a way that the dependences of the creep rate of UO 2 , UO 2 -PuO 2 and UN under irradiation on stress, temperature, fission rate, burn-up and porosity is obtained. It shows that creeping of fuels under irradiation at high temperatures is equivalent to thermally activated creeping, while at low temperature the creep rate induced by irradiation is much higher than the condition without irradiation. The increment of oxidic nuclear fuels is greater than in UN, the stress dependence on low burn-up is proportional in both cases, and the influence of temperature is quite small. (orig.) [de

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.

[Research Progress and Prospect of Applications of Finite Element Method in Lumbar Spine Biomechanics].

Science.gov (United States)

Zhang, Zhenjun; Li, Yang; Liao, Zhenhua; Liu, Weiqiang

2016-12-01

Based on the application of finite element analysis in spine biomechanics,the research progress of finite element method applied in lumbar spine mechanics is reviewed and the prospect is forecasted.The related works,including lumbar ontology modeling,clinical application research,and occupational injury and protection,are summarized.The main research areas of finite element method are as follows:new accurate modeling process,the optimized simulation method,diversified clinical effect evaluation,and the clinical application of artificial lumbar disc.According to the recent research progress,the application prospects of finite element method,such as automation and individuation of modeling process,evaluation and analysis of new operation methods and simulation of mechanical damage and dynamic response,are discussed.The purpose of this paper is to provide the theoretical reference and practical guidance for the clinical lumbar problems by reviewing the application of finite element method in the field of the lumbar spine biomechanics.

Finite element analyses for RF photoinjector gun cavities

International Nuclear Information System (INIS)

Marhauser, F.

2006-01-01

This paper details electromagnetical, thermal and structural 3D Finite Element Analyses (FEA) for normal conducting RF photoinjector gun cavities. The simulation methods are described extensively. Achieved results are presented. (orig.)

Finite element analyses for RF photoinjector gun cavities

Energy Technology Data Exchange (ETDEWEB)

Marhauser, F. [Berliner Elektronenspeicherring-Gesellschaft fuer Synchrotronstrahlung mbH (BESSY), Berlin (Germany)

2006-07-01

This paper details electromagnetical, thermal and structural 3D Finite Element Analyses (FEA) for normal conducting RF photoinjector gun cavities. The simulation methods are described extensively. Achieved results are presented. (orig.)

Validation Assessment of a Glass-to-Metal Seal Finite-Element Model

Energy Technology Data Exchange (ETDEWEB)

Jamison, Ryan Dale [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Buchheit, Thomas E. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Emery, John M [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Romero, Vicente J. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Stavig, Mark E. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Newton, Clay S. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Brown, Arthur [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2017-10-01

Sealing glasses are ubiquitous in high pressure and temperature engineering applications, such as hermetic feed-through electrical connectors. A common connector technology are glass-to-metal seals where a metal shell compresses a sealing glass to create a hermetic seal. Though finite-element analysis has been used to understand and design glass-to-metal seals for many years, there has been little validation of these models. An indentation technique was employed to measure the residual stress on the surface of a simple glass-to-metal seal. Recently developed rate- dependent material models of both Schott 8061 and 304L VAR stainless steel have been applied to a finite-element model of the simple glass-to-metal seal. Model predictions of residual stress based on the evolution of material models are shown. These model predictions are compared to measured data. Validity of the finite- element predictions is discussed. It will be shown that the finite-element model of the glass-to-metal seal accurately predicts the mean residual stress in the glass near the glass-to-metal interface and is valid for this quantity of interest.

A cohesive finite element formulation for modelling fracture and ...

Indian Academy of Sciences (India)

cohesive elements experience material softening and lose their stress carrying capacity. A few simple ..... In the present work, a Lagrangian finite element procedure is employed. In this formu clation ...... o, is related to 'c o by,. 't o='c o ¼ 1 ہ. 1.

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

International Nuclear Information System (INIS)

Dickenson, P.W.

1990-01-01

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

A methodology to analyze the creep behaviour of nuclear fuel waste containers

Energy Technology Data Exchange (ETDEWEB)

Dutton, R [Atomic Energy of Canada Ltd., Pinawa, MB (Canada). Whiteshell Labs.