A finite element field solver for dipole modes

International Nuclear Information System (INIS)

Nelson, E.M.

1992-01-01

A finite element field solver for dipole modes in axisymmetric structures has been written. The second-order elements used in this formulation yield accurate mode frequencies with no spurious modes. Quasi-periodic boundaries are included to allow travelling waves in periodic structures. The solver is useful in applications requiring precise frequency calculations such as detuned accelerator structures for linear colliders. Comparisons are made with measurements and with the popular but less accurate field solver URMEL. (author). 7 refs., 4 figs

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

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.

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)

A simple nodal force distribution method in refined finite element meshes

Energy Technology Data Exchange (ETDEWEB)

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

2017-05-15

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

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

OpenAIRE

Abubakri, K.; Veladi, H.

2016-01-01

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

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.

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

Finite Element Modelling of a Field-Sensed Magnetic Suspended System for Accurate Proximity Measurement Based on a Sensor Fusion Algorithm with Unscented Kalman Filter.

Science.gov (United States)

Chowdhury, Amor; Sarjaš, Andrej

2016-09-15

The presented paper describes accurate distance measurement for a field-sensed magnetic suspension system. The proximity measurement is based on a Hall effect sensor. The proximity sensor is installed directly on the lower surface of the electro-magnet, which means that it is very sensitive to external magnetic influences and disturbances. External disturbances interfere with the information signal and reduce the usability and reliability of the proximity measurements and, consequently, the whole application operation. A sensor fusion algorithm is deployed for the aforementioned reasons. The sensor fusion algorithm is based on the Unscented Kalman Filter, where a nonlinear dynamic model was derived with the Finite Element Modelling approach. The advantage of such modelling is a more accurate dynamic model parameter estimation, especially in the case when the real structure, materials and dimensions of the real-time application are known. The novelty of the paper is the design of a compact electro-magnetic actuator with a built-in low cost proximity sensor for accurate proximity measurement of the magnetic object. The paper successively presents a modelling procedure with the finite element method, design and parameter settings of a sensor fusion algorithm with Unscented Kalman Filter and, finally, the implementation procedure and results of real-time operation.

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

Peridynamic Multiscale Finite Element Methods

Energy Technology Data Exchange (ETDEWEB)

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

2015-12-01

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

Finite-element time evolution operator for the anharmonic oscillator

Science.gov (United States)

Milton, Kimball A.

1995-01-01

The finite-element approach to lattice field theory is both highly accurate (relative errors approximately 1/N(exp 2), where N is the number of lattice points) and exactly unitary (in the sense that canonical commutation relations are exactly preserved at the lattice sites). In this talk I construct matrix elements for dynamical variables and for the time evolution operator for the anharmonic oscillator, for which the continuum Hamiltonian is H = p(exp 2)/2 + lambda q(exp 4)/4. Construction of such matrix elements does not require solving the implicit equations of motion. Low order approximations turn out to be extremely accurate. For example, the matrix element of the time evolution operator in the harmonic oscillator ground state gives a results for the anharmonic oscillator ground state energy accurate to better than 1 percent, while a two-state approximation reduces the error to less than 0.1 percent.

3-D finite element analysis of claw-poled stepping motor

International Nuclear Information System (INIS)

Kawase, Yoshihiro; Yamaguchi, Tadashi; Mizuno; Koike, Yoshikazu

2002-01-01

Stepping motors are widely used for various electric instruments. It is necessary for the optimum design to analyze the magnetic field accurately. The 3-D finite element method with edge elements taking into account the rotation of the rotor has been applied to analyze the magnetic field of a claw-poled stepping motor. (Author)

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.

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

Nonlinear magnetohydrodynamics simulation using high-order finite elements

International Nuclear Information System (INIS)

Plimpton, Steven James; Schnack, D.D.; Tarditi, A.; Chu, M.S.; Gianakon, T.A.; Kruger, S.E.; Nebel, R.A.; Barnes, D.C.; Sovinec, C.R.; Glasser, A.H.

2005-01-01

A conforming representation composed of 2D finite elements and finite Fourier series is applied to 3D nonlinear non-ideal magnetohydrodynamics using a semi-implicit time-advance. The self-adjoint semi-implicit operator and variational approach to spatial discretization are synergistic and enable simulation in the extremely stiff conditions found in high temperature plasmas without sacrificing the geometric flexibility needed for modeling laboratory experiments. Growth rates for resistive tearing modes with experimentally relevant Lundquist number are computed accurately with time-steps that are large with respect to the global Alfven time and moderate spatial resolution when the finite elements have basis functions of polynomial degree (p) two or larger. An error diffusion method controls the generation of magnetic divergence error. Convergence studies show that this approach is effective for continuous basis functions with p (ge) 2, where the number of test functions for the divergence control terms is less than the number of degrees of freedom in the expansion for vector fields. Anisotropic thermal conduction at realistic ratios of parallel to perpendicular conductivity (x(parallel)/x(perpendicular)) is computed accurately with p (ge) 3 without mesh alignment. A simulation of tearing-mode evolution for a shaped toroidal tokamak equilibrium demonstrates the effectiveness of the algorithm in nonlinear conditions, and its results are used to verify the accuracy of the numerical anisotropic thermal conduction in 3D magnetic topologies.

Finite Element Method for Analysis of Material Properties

DEFF Research Database (Denmark)

Rauhe, Jens Christian

and the finite element method. The material microstructure of the heterogeneous material is non-destructively determined using X-ray microtomography. A software program has been generated which uses the X-ray tomographic data as an input for the mesh generation of the material microstructure. To obtain a proper...... which are used for the determination of the effective properties of the heterogeneous material. Generally, the properties determined using the finite element method coupled with X-ray microtomography are in good agreement with both experimentally determined properties and properties determined using......The use of cellular and composite materials have in recent years become more and more common in all kinds of structural components and accurate knowledge of the effective properties is therefore essential. In this wok the effective properties are determined using the real material microstructure...

[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 elements and approximation

CERN Document Server

Zienkiewicz, O C

2006-01-01

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

A first course in finite elements

CERN Document Server

Fish, Jacob

2007-01-01

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

Mixed hybrid finite elements and streamline computation for the potential flow problem

NARCIS (Netherlands)

Kaasschieter, E.F.; Huijben, A.J.M.

1992-01-01

An important class of problems in mathematical physics involves equations of the form -¿ · (A¿¿) = f. In a variety of problems it is desirable to obtain an accurate approximation of the flow quantity u = -A¿¿. Such an accurate approximation can be determined by the mixed finite element method. In

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

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.

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.

Non-linear finite element modeling

DEFF Research Database (Denmark)

Mikkelsen, Lars Pilgaard

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

Structural modeling techniques by finite element method

International Nuclear Information System (INIS)

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

1991-01-01

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

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.

A Posteriori Finite Element Bounds for Sensitivity Derivatives of Partial-Differential-Equation Outputs. Revised

Science.gov (United States)

Lewis, Robert Michael; Patera, Anthony T.; Peraire, Jaume

1998-01-01

We present a Neumann-subproblem a posteriori finite element procedure for the efficient and accurate calculation of rigorous, 'constant-free' upper and lower bounds for sensitivity derivatives of functionals of the solutions of partial differential equations. The design motivation for sensitivity derivative error control is discussed; the a posteriori finite element procedure is described; the asymptotic bounding properties and computational complexity of the method are summarized; and illustrative numerical results are presented.

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

Directory of Open Access Journals (Sweden)

Zahari Siti Norazila

2017-01-01

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

Finite element calculation of the interaction energy of shape memory alloy

International Nuclear Information System (INIS)

Yang, Seung Yong

2004-01-01

Strain energy due to the mechanical interaction between self-accommodation groups of martensitic phase transformation is called interaction energy. Evaluation of the interaction energy should be accurate since the energy appears in constitutive models for predicting the mechanical behavior of shape memory alloy. In this paper, the interaction energy is evaluated in terms of theoretical formulation and explicit finite element calculation. A simple example with two habit plane variants was considered. It was shown that the theoretical formulation assuming elastic interaction between the self-accommodation group and matrix gives larger interaction energy than explicit finite element calculation in which transformation softening is accounted for

Accurate Modeling of a Transverse Flux Permanent Magnet Generator Using 3D Finite Element Analysis

DEFF Research Database (Denmark)

Hosseini, Seyedmohsen; Moghani, Javad Shokrollahi; Jensen, Bogi Bech

2011-01-01

method is then proposed that reveals the behavior of the generator under any load. Finally, torque calculations are carried out using three dimensional finite element analyses. It is shown that although in the single-phase generator the cogging torque is very high, this can be improved significantly...

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.

Hermitian Mindlin Plate Wavelet Finite Element Method for Load Identification

Directory of Open Access Journals (Sweden)

Xiaofeng Xue

2016-01-01

Full Text Available A new Hermitian Mindlin plate wavelet element is proposed. The two-dimensional Hermitian cubic spline interpolation wavelet is substituted into finite element functions to construct frequency response function (FRF. It uses a system’s FRF and response spectrums to calculate load spectrums and then derives loads in the time domain via the inverse fast Fourier transform. By simulating different excitation cases, Hermitian cubic spline wavelets on the interval (HCSWI finite elements are used to reverse load identification in the Mindlin plate. The singular value decomposition (SVD method is adopted to solve the ill-posed inverse problem. Compared with ANSYS results, HCSWI Mindlin plate element can accurately identify the applied load. Numerical results show that the algorithm of HCSWI Mindlin plate element is effective. The accuracy of HCSWI can be verified by comparing the FRF of HCSWI and ANSYS elements with the experiment data. The experiment proves that the load identification of HCSWI Mindlin plate is effective and precise by using the FRF and response spectrums to calculate the loads.

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

Incorporation of exact boundary conditions into a discontinuous galerkin finite element method for accurately solving 2d time-dependent maxwell equations

KAUST Repository

Sirenko, Kostyantyn

2013-01-01

A scheme that discretizes exact absorbing boundary conditions (EACs) to incorporate them into a time-domain discontinuous Galerkin finite element method (TD-DG-FEM) is described. The proposed TD-DG-FEM with EACs is used for accurately characterizing transient electromagnetic wave interactions on two-dimensional waveguides. Numerical results demonstrate the proposed method\\'s superiority over the TD-DG-FEM that employs approximate boundary conditions and perfectly matched layers. Additionally, it is shown that the proposed method can produce the solution with ten-eleven digit accuracy when high-order spatial basis functions are used to discretize the Maxwell equations as well as the EACs. © 1963-2012 IEEE.

Finite element model updating of a small steel frame using neural networks

International Nuclear Information System (INIS)

Zapico, J L; González, M P; Alonso, R; González-Buelga, A

2008-01-01

This paper presents an experimental and analytical dynamic study of a small-scale steel frame. The experimental model was physically built and dynamically tested on a shaking table in a series of different configurations obtained from the original one by changing the mass and by causing structural damage. Finite element modelling and parameterization with physical meaning is iteratively tried for the original undamaged configuration. The finite element model is updated through a neural network, the natural frequencies of the model being the net input. The updating process is made more accurate and robust by using a regressive procedure, which constitutes an original contribution of this work. A novel simplified analytical model has been developed to evaluate the reduction of bending stiffness of the elements due to damage. The experimental results of the rest of the configurations have been used to validate both the updated finite element model and the analytical one. The statistical properties of the identified modal data are evaluated. From these, the statistical properties and a confidence interval for the estimated model parameters are obtained by using the Latin Hypercube sampling technique. The results obtained are successful: the updated model accurately reproduces the low modes identified experimentally for all configurations, and the statistical study of the transmission of errors yields a narrow confidence interval for all the identified parameters

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

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 cover method with mortar elements for elastoplasticity problems

Science.gov (United States)

Kurumatani, M.; Terada, K.

2005-06-01

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

Electromagnetic Extended Finite Elements for High-Fidelity Multimaterial Problems LDRD Final Report

Energy Technology Data Exchange (ETDEWEB)

Siefert, Christopher [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Bochev, Pavel Blagoveston [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Kramer, Richard Michael Jack [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Voth, Thomas Eugene [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Cox, James [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2014-09-01

Surface effects are critical to the accurate simulation of electromagnetics (EM) as current tends to concentrate near material surfaces. Sandia EM applications, which include exploding bridge wires for detonator design, electromagnetic launch of flyer plates for material testing and gun design, lightning blast-through for weapon safety, electromagnetic armor, and magnetic flux compression generators, all require accurate resolution of surface effects. These applications operate in a large deformation regime, where body-fitted meshes are impractical and multimaterial elements are the only feasible option. State-of-the-art methods use various mixture models to approximate the multi-physics of these elements. The empirical nature of these models can significantly compromise the accuracy of the simulation in this very important surface region. We propose to substantially improve the predictive capability of electromagnetic simulations by removing the need for empirical mixture models at material surfaces. We do this by developing an eXtended Finite Element Method (XFEM) and an associated Conformal Decomposition Finite Element Method (CDFEM) which satisfy the physically required compatibility conditions at material interfaces. We demonstrate the effectiveness of these methods for diffusion and diffusion-like problems on node, edge and face elements in 2D and 3D. We also present preliminary work on h -hierarchical elements and remap algorithms.

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

CERN Document Server

De Ford, John F; Petillo, John

2005-01-01

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

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

International Nuclear Information System (INIS)

Cheung, K.Y.

1985-01-01

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

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.

Use of the finite element displacement method to solve solid-fluid interaction vibration problems

International Nuclear Information System (INIS)

Brown, S.J.; Hsu, K.H.

1978-01-01

It is shown through comparison to experimental, theoretical, and other finite element formulations that the finite element displacement method can solve accurately and economically a certain class of solid-fluid eigenvalue problems. The problems considered are small displacements in the absence of viscous damping and are 2-D and 3-D in nature. In this study the advantages of the finite element method (in particular the displacement formulation) is apparent in that a large structure consisting of the cylinders, support flanges, fluid, and other experimental boundaries could be modeled to yield good correlation to experimental data. The ability to handle large problems with standard structural programs is the key advantage of the displacement fluid method. The greatest obstacle is the inability of the analyst to inhibit those rotational degrees of freedom that are unnecessary to his fluid-structure vibration problem. With judicious use of element formulation, boundary conditions and modeling, the displacement finite element method can be successfully used to predict solid-fluid response to vibration and seismic loading

Finite element methods a practical guide

CERN Document Server

Whiteley, Jonathan

2017-01-01

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

Finite element analysis of piezoelectric materials

International Nuclear Information System (INIS)

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

1999-01-01

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

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 Timoshenko Piezoelectric Beam Finite Element with Consistent Performance Irrespective of Geometric and Material Configurations

Directory of Open Access Journals (Sweden)

Litesh N. Sulbhewar

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

A multigrid solution method for mixed hybrid finite elements

Energy Technology Data Exchange (ETDEWEB)

Schmid, W. [Universitaet Augsburg (Germany)

1996-12-31

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

A deep learning approach to estimate stress distribution: a fast and accurate surrogate of finite-element analysis.

Science.gov (United States)

Liang, Liang; Liu, Minliang; Martin, Caitlin; Sun, Wei

2018-01-01

Structural finite-element analysis (FEA) has been widely used to study the biomechanics of human tissues and organs, as well as tissue-medical device interactions, and treatment strategies. However, patient-specific FEA models usually require complex procedures to set up and long computing times to obtain final simulation results, preventing prompt feedback to clinicians in time-sensitive clinical applications. In this study, by using machine learning techniques, we developed a deep learning (DL) model to directly estimate the stress distributions of the aorta. The DL model was designed and trained to take the input of FEA and directly output the aortic wall stress distributions, bypassing the FEA calculation process. The trained DL model is capable of predicting the stress distributions with average errors of 0.492% and 0.891% in the Von Mises stress distribution and peak Von Mises stress, respectively. This study marks, to our knowledge, the first study that demonstrates the feasibility and great potential of using the DL technique as a fast and accurate surrogate of FEA for stress analysis. © 2018 The Author(s).

Domain decomposition methods for mortar finite elements

Energy Technology Data Exchange (ETDEWEB)

Widlund, O.

1996-12-31

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

An assessment of finite element synthesis model using 3D FEM

International Nuclear Information System (INIS)

Jagannathan, V.; Rastogi, B.P

1983-01-01

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

A Floating Node Method for the Modelling of Discontinuities Within a Finite Element

Science.gov (United States)

Pinho, Silvestre T.; Chen, B. Y.; DeCarvalho, Nelson V.; Baiz, P. M.; Tay, T. E.

2013-01-01

This paper focuses on the accurate numerical representation of complex networks of evolving discontinuities in solids, with particular emphasis on cracks. The limitation of the standard finite element method (FEM) in approximating discontinuous solutions has motivated the development of re-meshing, smeared crack models, the eXtended Finite Element Method (XFEM) and the Phantom Node Method (PNM). We propose a new method which has some similarities to the PNM, but crucially: (i) does not introduce an error on the crack geometry when mapping to natural coordinates; (ii) does not require numerical integration over only part of a domain; (iii) can incorporate weak discontinuities and cohesive cracks more readily; (iv) is ideally suited for the representation of multiple and complex networks of (weak, strong and cohesive) discontinuities; (v) leads to the same solution as a finite element mesh where the discontinuity is represented explicitly; and (vi) is conceptually simpler than the PNM.

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.

Automatic procedure for realistic 3D finite element modelling of human brain for bioelectromagnetic computations

International Nuclear Information System (INIS)

Aristovich, K Y; Khan, S H

2010-01-01

Realistic computer modelling of biological objects requires building of very accurate and realistic computer models based on geometric and material data, type, and accuracy of numerical analyses. This paper presents some of the automatic tools and algorithms that were used to build accurate and realistic 3D finite element (FE) model of whole-brain. These models were used to solve the forward problem in magnetic field tomography (MFT) based on Magnetoencephalography (MEG). The forward problem involves modelling and computation of magnetic fields produced by human brain during cognitive processing. The geometric parameters of the model were obtained from accurate Magnetic Resonance Imaging (MRI) data and the material properties - from those obtained from Diffusion Tensor MRI (DTMRI). The 3D FE models of the brain built using this approach has been shown to be very accurate in terms of both geometric and material properties. The model is stored on the computer in Computer-Aided Parametrical Design (CAD) format. This allows the model to be used in a wide a range of methods of analysis, such as finite element method (FEM), Boundary Element Method (BEM), Monte-Carlo Simulations, etc. The generic model building approach presented here could be used for accurate and realistic modelling of human brain and many other biological objects.

Finite-element semi-discretization of linearized compressible and resistive MHD

International Nuclear Information System (INIS)

Kerner, W.; Jakoby, A.; Lerbinger, K.

1985-08-01

The full resistive MHD equations are linearized around an equilibrium with cylindrical symmetry and solved numerically as an initial-value problem. The semi-discretization using cubic and quadratic finite elements for the spatial discretization and a fully implicit time advance yields very accurate results even for small values of the resistivity. In the application different phenomena such as waves, resistive instabilities and overstable modes are addressed. (orig.)

Scientific use of the finite element method in Orthodontics

Science.gov (United States)

Knop, Luegya; Gandini, Luiz Gonzaga; Shintcovsk, Ricardo Lima; Gandini, Marcia Regina Elisa Aparecida Schiavon

2015-01-01

INTRODUCTION: The finite element method (FEM) is an engineering resource applied to calculate the stress and deformation of complex structures, and has been widely used in orthodontic research. With the advantage of being a non-invasive and accurate method that provides quantitative and detailed data on the physiological reactions possible to occur in tissues, applying the FEM can anticipate the visualization of these tissue responses through the observation of areas of stress created from applied orthodontic mechanics. OBJECTIVE: This article aims at reviewing and discussing the stages of the finite element method application and its applicability in Orthodontics. RESULTS: FEM is able to evaluate the stress distribution at the interface between periodontal ligament and alveolar bone, and the shifting trend in various types of tooth movement when using different types of orthodontic devices. Therefore, it is necessary to know specific software for this purpose. CONCLUSIONS: FEM is an important experimental method to answer questions about tooth movement, overcoming the disadvantages of other experimental methods. PMID:25992996

Finite-element analysis of flawed and unflawed pipe tests

International Nuclear Information System (INIS)

James, R.J.; Nickell, R.E.; Sullaway, M.F.

1989-12-01

Contemporary versions of the general purpose, nonlinear finite element program ABAQUS have been used in structural response verification exercises on flawed and unflawed austenitic stainless steel and ferritic steel piping. Among the topics examined, through comparison between ABAQUS calculations and test results, were: (1) the effect of using variations in the stress-strain relationship from the test article material on the calculated response; (2) the convergence properties of various finite element representations of the pipe geometry, using shell, beam and continuum models; (3) the effect of test system compliance; and (4) the validity of ABAQUS J-integral routines for flawed pipe evaluations. The study was culminated by the development and demonstration of a ''macroelement'' representation for the flawed pipe section. The macroelement can be inserted into an existing piping system model, in order to accurately treat the crack-opening and crack-closing static and dynamic response. 11 refs., 20 figs., 1 tab

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.

A novel hybrid stress-function finite element method immune to severe mesh distortion

International Nuclear Information System (INIS)

Cen Song; Zhou Mingjue; Fu Xiangrong

2010-01-01

This paper introduces a hybrid stress-function finite element method proposed recently for developing 2D finite element models immune to element shapes. Deferent from the first version of the hybrid-stress element constructed by Pian, the stress function φ of 2D elastic or fracture problem is regarded as the functional variable of the complementary energy functional. Then, the basic analytical solutions of φ are taken as the trial functions for finite element models, and meanwhile, the corresponding unknown stress-function constants are introduced. By using the principle of minimum complementary energy, these unknown stress-function constants can be expressed in terms of the displacements along element edges. Finally, the complementary energy functional can be rewritten in terms of element nodal displacement vector, and thus, the element stiffness matrix of such hybrid-function element can be obtained. As examples, two (8- and 12-node) quadrilateral plane elements and an arbitrary polygonal crack element are constructed by employing different basic analytical solutions of different stress functions. Numerical results show that, the 8- and 12-node plane models can produce the exact solutions for pure bending and linear bending problems, respectively, even the element shape degenerates into triangle and concave quadrangle; and the crack element can also predict accurate results with very low computational cost in analysis of stress-singularity problems.

Band-limited Green's Functions for Quantitative Evaluation of Acoustic Emission Using the Finite Element Method

Science.gov (United States)

Leser, William P.; Yuan, Fuh-Gwo; Leser, William P.

2013-01-01

A method of numerically estimating dynamic Green's functions using the finite element method is proposed. These Green's functions are accurate in a limited frequency range dependent on the mesh size used to generate them. This range can often match or exceed the frequency sensitivity of the traditional acoustic emission sensors. An algorithm is also developed to characterize an acoustic emission source by obtaining information about its strength and temporal dependence. This information can then be used to reproduce the source in a finite element model for further analysis. Numerical examples are presented that demonstrate the ability of the band-limited Green's functions approach to determine the moment tensor coefficients of several reference signals to within seven percent, as well as accurately reproduce the source-time function.

Probabilistic finite elements for fracture mechanics

Science.gov (United States)

Besterfield, Glen

1988-01-01

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

Finite element analysis of beam-to-column joints in steel frames under cyclic loading

Directory of Open Access Journals (Sweden)

Elsayed Mashaly

2011-03-01

Full Text Available The aim of this paper is to present a simple and accurate three-dimensional (3D finite element model (FE capable of predicting the actual behavior of beam-to-column joints in steel frames subjected to lateral loads. The software package ANSYS is used to model the joint. The bolted extended-end-plate connection was chosen as an important type of beam–column joints. The extended-end-plate connection is chosen for its complexity in the analysis and behavior due to the number of connection components and their inheritable non-linear behavior. Two experimental tests in the literature were chosen to verify the finite element model. The results of both the experimental and the proposed finite element were compared. One of these tests was monotonically loaded, whereas the second was cyclically loaded. The finite element model is improved to enhance the defects of the finite element model used. These defects are; the long time need for the analysis and the inability of the contact element type to follow the behavior of moment–rotation curve under cyclic loading. As a contact element, the surface-to-surface element is used instead of node-to-node element to enhance the model. The FE results show good correlation with the experimental one. An attempt to improve a new technique for modeling bolts is conducted. The results show that this technique is supposed to avoid the defects above, give much less elements number and less solution time than the other modeling techniques.

Introduction to finite and spectral element methods using Matlab

CERN Document Server

Pozrikidis, Constantine

2014-01-01

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

Finite Element Methods and Their Applications

CERN Document Server

Chen, Zhangxin

2005-01-01

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

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

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

Compatible-strain mixed finite element methods for incompressible nonlinear elasticity

Science.gov (United States)

Faghih Shojaei, Mostafa; Yavari, Arash

2018-05-01

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

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.

A multipoint flux mixed finite element method on distorted quadrilaterals and hexahedra

KAUST Repository

Wheeler, Mary

2011-11-06

In this paper, we develop a new mixed finite element method for elliptic problems on general quadrilateral and hexahedral grids that reduces to a cell-centered finite difference scheme. A special non-symmetric quadrature rule is employed that yields a positive definite cell-centered system for the pressure by eliminating local velocities. The method is shown to be accurate on highly distorted rough quadrilateral and hexahedral grids, including hexahedra with non-planar faces. Theoretical and numerical results indicate first-order convergence for the pressure and face fluxes. © 2011 Springer-Verlag.

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

KAUST Repository

Niemi, Antti

2011-06-03

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

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)

Mesh-morphing algorithms for specimen-specific finite element modeling.

Science.gov (United States)

Sigal, Ian A; Hardisty, Michael R; Whyne, Cari M

2008-01-01

Despite recent advances in software for meshing specimen-specific geometries, considerable effort is still often required to produce and analyze specimen-specific models suitable for biomechanical analysis through finite element modeling. We hypothesize that it is possible to obtain accurate models by adapting a pre-existing geometry to represent a target specimen using morphing techniques. Here we present two algorithms for morphing, automated wrapping (AW) and manual landmarks (ML), and demonstrate their use to prepare specimen-specific models of caudal rat vertebrae. We evaluate the algorithms by measuring the distance between target and morphed geometries and by comparing response to axial loading simulated with finite element (FE) methods. First a traditional reconstruction process based on microCT was used to obtain two natural specimen-specific FE models. Next, the two morphing algorithms were used to compute mappings from the surface of one model, the source, to the other, the target, and to use this mapping to morph the source mesh to produce a target mesh. The microCT images were then used to assign element-specific material properties. In AW the mappings were obtained by wrapping the source and target surfaces with an auxiliary triangulated surface. In ML, landmarks were manually placed on corresponding locations on the surfaces of both source and target. Both morphing algorithms were successful in reproducing the shape of the target vertebra with a median distance between natural and morphed models of 18.8 and 32.2 microm, respectively, for AW and ML. Whereas AW-morphing produced a surface more closely resembling that of the target, ML guaranteed correspondence of the landmark locations between source and target. Morphing preserved the quality of the mesh producing models suitable for FE simulation. Moreover, there were only minor differences between natural and morphed models in predictions of deformation, strain and stress. We therefore conclude that

ANSYS mechanical APDL for finite element analysis

CERN Document Server

Thompson, Mary Kathryn

2017-01-01

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

Finite elements of nonlinear continua

CERN Document Server

Oden, John Tinsley

1972-01-01

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

A locally conservative non-negative finite element formulation for anisotropic advective-diffusive-reactive systems

Science.gov (United States)

Mudunuru, M. K.; Shabouei, M.; Nakshatrala, K.

2015-12-01

Advection-diffusion-reaction (ADR) equations appear in various areas of life sciences, hydrogeological systems, and contaminant transport. Obtaining stable and accurate numerical solutions can be challenging as the underlying equations are coupled, nonlinear, and non-self-adjoint. Currently, there is neither a robust computational framework available nor a reliable commercial package known that can handle various complex situations. Herein, the objective of this poster presentation is to present a novel locally conservative non-negative finite element formulation that preserves the underlying physical and mathematical properties of a general linear transient anisotropic ADR equation. In continuous setting, governing equations for ADR systems possess various important properties. In general, all these properties are not inherited during finite difference, finite volume, and finite element discretizations. The objective of this poster presentation is two fold: First, we analyze whether the existing numerical formulations (such as SUPG and GLS) and commercial packages provide physically meaningful values for the concentration of the chemical species for various realistic benchmark problems. Furthermore, we also quantify the errors incurred in satisfying the local and global species balance for two popular chemical kinetics schemes: CDIMA (chlorine dioxide-iodine-malonic acid) and BZ (Belousov--Zhabotinsky). Based on these numerical simulations, we show that SUPG and GLS produce unphysical values for concentration of chemical species due to the violation of the non-negative constraint, contain spurious node-to-node oscillations, and have large errors in local and global species balance. Second, we proposed a novel finite element formulation to overcome the above difficulties. The proposed locally conservative non-negative computational framework based on low-order least-squares finite elements is able to preserve these underlying physical and mathematical properties

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

International Nuclear Information System (INIS)

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

2005-01-01

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

The finite element method in engineering, 2nd edition

International Nuclear Information System (INIS)

Rao, S.S.

1986-01-01

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

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

Energy Technology Data Exchange (ETDEWEB)

Lee, Sang Jin; Seo, Jeong Moon

2000-08-01

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

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

International Nuclear Information System (INIS)

Lee, Sang Jin; Seo, Jeong Moon

2000-08-01

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

Extended Finite Element Method with Simplified Spherical Harmonics Approximation for the Forward Model of Optical Molecular Imaging

Directory of Open Access Journals (Sweden)

Wei Li

2012-01-01

Full Text Available An extended finite element method (XFEM for the forward model of 3D optical molecular imaging is developed with simplified spherical harmonics approximation (SPN. In XFEM scheme of SPN equations, the signed distance function is employed to accurately represent the internal tissue boundary, and then it is used to construct the enriched basis function of the finite element scheme. Therefore, the finite element calculation can be carried out without the time-consuming internal boundary mesh generation. Moreover, the required overly fine mesh conforming to the complex tissue boundary which leads to excess time cost can be avoided. XFEM conveniences its application to tissues with complex internal structure and improves the computational efficiency. Phantom and digital mouse experiments were carried out to validate the efficiency of the proposed method. Compared with standard finite element method and classical Monte Carlo (MC method, the validation results show the merits and potential of the XFEM for optical imaging.

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

International Nuclear Information System (INIS)

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

1991-01-01

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

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

Nonlinear finite element modeling of corrugated board

Science.gov (United States)

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

1999-01-01

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

Bibliography for finite elements. [2200 references

Energy Technology Data Exchange (ETDEWEB)

Whiteman, J R [comp.

1975-01-01

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

Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media

Energy Technology Data Exchange (ETDEWEB)

Gao, Kai, E-mail: kaigao87@gmail.com [Department of Geology and Geophysics, Texas A& M University, College Station, TX 77843 (United States); Fu, Shubin, E-mail: shubinfu89@gmail.com [Department of Mathematics, Texas A& M University, College Station, TX 77843 (United States); Gibson, Richard L., E-mail: gibson@tamu.edu [Department of Geology and Geophysics, Texas A& M University, College Station, TX 77843 (United States); Chung, Eric T., E-mail: tschung@math.cuhk.edu.hk [Department of Mathematics, The Chinese University of Hong Kong, Shatin, NT (Hong Kong); Efendiev, Yalchin, E-mail: efendiev@math.tamu.edu [Department of Mathematics, Texas A& M University, College Station, TX 77843 (United States); Numerical Porous Media SRI Center (NumPor), King Abdullah University of Science and Technology, Thuwal (Saudi Arabia)

2015-08-15

It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.

Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media

International Nuclear Information System (INIS)

Gao, Kai; Fu, Shubin; Gibson, Richard L.; Chung, Eric T.; Efendiev, Yalchin

2015-01-01

It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system

Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media

KAUST Repository

Gao, Kai

2015-04-14

It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both boundaries and the interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.

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.

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

Development and experimental verification of a finite element method for accurate analysis of a surface acoustic wave device

Science.gov (United States)

Mohibul Kabir, K. M.; Matthews, Glenn I.; Sabri, Ylias M.; Russo, Salvy P.; Ippolito, Samuel J.; Bhargava, Suresh K.

2016-03-01

Accurate analysis of surface acoustic wave (SAW) devices is highly important due to their use in ever-growing applications in electronics, telecommunication and chemical sensing. In this study, a novel approach for analyzing the SAW devices was developed based on a series of two-dimensional finite element method (FEM) simulations, which has been experimentally verified. It was found that the frequency response of the two SAW device structures, each having slightly different bandwidth and center lobe characteristics, can be successfully obtained utilizing the current density of the electrodes via FEM simulations. The two SAW structures were based on XY Lithium Niobate (LiNbO3) substrates and had two and four electrode finger pairs in both of their interdigital transducers, respectively. Later, SAW devices were fabricated in accordance with the simulated models and their measured frequency responses were found to correlate well with the obtained simulations results. The results indicated that better match between calculated and measured frequency response can be obtained when one of the input electrode finger pairs was set at zero volts and all the current density components were taken into account when calculating the frequency response of the simulated SAW device structures.

FEHM, Finite Element Heat and Mass Transfer Code

International Nuclear Information System (INIS)

Zyvoloski, G.A.

2002-01-01

1 - Description of program or function: FEHM is a numerical simulation code for subsurface transport processes. It models 3-D, time-dependent, multiphase, multicomponent, non-isothermal, reactive flow through porous and fractured media. It can accurately represent complex 3-D geologic media and structures and their effects on subsurface flow and transport. Its capabilities include flow of gas, water, and heat; flow of air, water, and heat; multiple chemically reactive and sorbing tracers; finite element/finite volume formulation; coupled stress module; saturated and unsaturated media; and double porosity and double porosity/double permeability capabilities. 2 - Methods: FEHM uses a preconditioned conjugate gradient solution of coupled linear equations and a fully implicit, fully coupled Newton Raphson solution of nonlinear equations. It has the capability of simulating transport using either a advection/diffusion solution or a particle tracking method. 3 - Restriction on the complexity of the problem: Disk space and machine memory are the only limitations

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

Finite element modeling of nonlinear piezoelectric energy harvesters with magnetic interaction

International Nuclear Information System (INIS)

Upadrashta, Deepesh; Yang, Yaowen

2015-01-01

Piezoelectric energy harvesting from ambient vibrations is a potential technology for powering wireless sensors and low power electronic devices. The conventional linear harvesters suffer from narrow operational bandwidth. Many attempts have been made especially using the magnetic interaction to broaden the bandwidth of harvesters. The finite element (FE) modeling has been used only for analyzing the linear harvesters in the literature. The main difficulties in extending the FE modeling to analyze the nonlinear harvesters involving magnetic interaction are developing the mesh needed for magnetic interaction in dynamic problems and the high demand on computational resource needed for solving the coupled electrical–mechanical–magnetic problem. In this paper, an innovative method is proposed to model the magnetic interaction without inclusion of the magnetic module. The magnetic force is modeled using the nonlinear spring element available in ANSYS finite element analysis (FEA) package, thus simplifying the simulation of nonlinear piezoelectric energy harvesters as an electromechanically coupled problem. Firstly, an FE model of a monostable nonlinear harvester with cantilever configuration is developed and the results are validated with predictions from the theoretical model. Later, the proposed technique of FE modeling is extended to a complex 2-degree of freedom nonlinear energy harvester for which an accurate analytical model is difficult to derive. The performance predictions from FEA are compared with the experimental results. It is concluded that the proposed modeling technique is able to accurately analyze the behavior of nonlinear harvesters with magnetic interaction. (paper)

Dispersion analysis of the Pn -Pn-1DG mixed finite element pair for atmospheric modelling

Science.gov (United States)

Melvin, Thomas

2018-02-01

Mixed finite element methods provide a generalisation of staggered grid finite difference methods with a framework to extend the method to high orders. The ability to generate a high order method is appealing for applications on the kind of quasi-uniform grids that are popular for atmospheric modelling, so that the method retains an acceptable level of accuracy even around special points in the grid. The dispersion properties of such schemes are important to study as they provide insight into the numerical adjustment to imbalance that is an important component in atmospheric modelling. This paper extends the recent analysis of the P2 - P1DG pair, that is a quadratic continuous and linear discontinuous finite element pair, to higher polynomial orders and also spectral element type pairs. In common with the previously studied element pair, and also with other schemes such as the spectral element and discontinuous Galerkin methods, increasing the polynomial order is found to provide a more accurate dispersion relation for the well resolved part of the spectrum but at the cost of a number of unphysical spectral gaps. The effects of these spectral gaps are investigated and shown to have a varying impact depending upon the width of the gap. Finally, the tensor product nature of the finite element spaces is exploited to extend the dispersion analysis into two-dimensions.

Computation of stress intensity factors for nozzle corner cracks by various finite element procedures

International Nuclear Information System (INIS)

Broekhoven, M.J.G.

1975-01-01

The present study aims at deriving accurate K-factors for a series of 5 elliptical nozzle corner cracks of increasing size by various finite element procedures, using a three-level recursive substructuring scheme to perform the computations in an economic way on an intermediate size computer (IBM 360/65 system). A nozzle on a flat plate has been selected for subsequent experimental verification, this configuration being considered an adequate simulation of a nozzle on a shallow shell. The computations have been performed with the ASKA finite element system using mainly HEXEC-27 (incomplete quartic) elements. The geometry has been subdivided into 5 subnets with a total of 3515 nodal points and 6250 unknowns, two main nets and one hyper net. Each crack front is described by 11 nodal points and all crack front nodes are inserted in the hyper net, which allows for the realization of the successive crack geometries by changing only a relatively small hyper net (615 to 725 unknowns). Output data have been interpreted in terms of K-factors by the global energy method, the displacement method and the stress method. Besides, a stiffness derivative procedure, recently developed at Brown University, which takes full advantage of the finite element formulation to calculate local K-factors, has been applied. Finally it has been investigated whether sufficiently accurate results can be obtained by analyzing a considerably smaller part than one half of the geometry (as strictly required by symmetry considerations), using fixed boundary conditions derived from a far cheaper analysis of the uncracked structure

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

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

International Nuclear Information System (INIS)

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

1988-01-01

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

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

A high-order multiscale finite-element method for time-domain acoustic-wave modeling

Science.gov (United States)

Gao, Kai; Fu, Shubin; Chung, Eric T.

2018-05-01

Accurate and efficient wave equation modeling is vital for many applications in such as acoustics, electromagnetics, and seismology. However, solving the wave equation in large-scale and highly heterogeneous models is usually computationally expensive because the computational cost is directly proportional to the number of grids in the model. We develop a novel high-order multiscale finite-element method to reduce the computational cost of time-domain acoustic-wave equation numerical modeling by solving the wave equation on a coarse mesh based on the multiscale finite-element theory. In contrast to existing multiscale finite-element methods that use only first-order multiscale basis functions, our new method constructs high-order multiscale basis functions from local elliptic problems which are closely related to the Gauss-Lobatto-Legendre quadrature points in a coarse element. Essentially, these basis functions are not only determined by the order of Legendre polynomials, but also by local medium properties, and therefore can effectively convey the fine-scale information to the coarse-scale solution with high-order accuracy. Numerical tests show that our method can significantly reduce the computation time while maintain high accuracy for wave equation modeling in highly heterogeneous media by solving the corresponding discrete system only on the coarse mesh with the new high-order multiscale basis functions.

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

International Nuclear Information System (INIS)

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

1981-01-01

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

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.

A Finite Element Method for Simulation of Compressible Cavitating Flows

Science.gov (United States)

Shams, Ehsan; Yang, Fan; Zhang, Yu; Sahni, Onkar; Shephard, Mark; Oberai, Assad

2016-11-01

This work focuses on a novel approach for finite element simulations of multi-phase flows which involve evolving interface with phase change. Modeling problems, such as cavitation, requires addressing multiple challenges, including compressibility of the vapor phase, interface physics caused by mass, momentum and energy fluxes. We have developed a mathematically consistent and robust computational approach to address these problems. We use stabilized finite element methods on unstructured meshes to solve for the compressible Navier-Stokes equations. Arbitrary Lagrangian-Eulerian formulation is used to handle the interface motions. Our method uses a mesh adaptation strategy to preserve the quality of the volumetric mesh, while the interface mesh moves along with the interface. The interface jump conditions are accurately represented using a discontinuous Galerkin method on the conservation laws. Condensation and evaporation rates at the interface are thermodynamically modeled to determine the interface velocity. We will present initial results on bubble cavitation the behavior of an attached cavitation zone in a separated boundary layer. We acknowledge the support from Army Research Office (ARO) under ARO Grant W911NF-14-1-0301.

Elastic-plastic finite element analyses for reducers with constant-depth internal circumferential surface cracks

International Nuclear Information System (INIS)

Wu, Szu-Ying; Tsai, Bor-Jiun; Chen, Jien-Jong

2015-01-01

In this study, a 3-D automatic elastic-plastic finite element mesh generator is established to accurately predict the J-integral value of an arbitrary reducer with a constant-depth internal circumferential surface crack under bending and axial force. The contact pairs are used on the crack surfaces to simulate the actual contact behaviors of the crack model under loadings. In order to verify the accuracy of the proposed elastic-plastic finite element model for a reducer with a surface crack, the cracked straight pipe models are generated according to a special modeling procedure for a flawed reducer. The J-integral values along the crack front of surface crack are calculated and compared with the straight pipe models which have been verified in the previous published studies. Based on the comparison of computed results, good agreements are obtained to show the accuracy of present numerical models. More confidence on using the 3-D elastic-plastic finite element analysis for reducers with internal circumferential surface cracks can be thus established in this work

Finite Element Modelling of Seismic Liquefaction in Soils

NARCIS (Netherlands)

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

2013-01-01

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

A globally well-posed finite element algorithm for aerodynamics applications

Science.gov (United States)

Iannelli, G. S.; Baker, A. J.

1991-01-01

A finite element CFD algorithm is developed for Euler and Navier-Stokes aerodynamic applications. For the linear basis, the resultant approximation is at least second-order-accurate in time and space for synergistic use of three procedures: (1) a Taylor weak statement, which provides for derivation of companion conservation law systems with embedded dispersion-error control mechanisms; (2) a stiffly stable second-order-accurate implicit Rosenbrock-Runge-Kutta temporal algorithm; and (3) a matrix tensor product factorization that permits efficient numerical linear algebra handling of the terminal large-matrix statement. Thorough analyses are presented regarding well-posed boundary conditions for inviscid and viscous flow specifications. Numerical solutions are generated and compared for critical evaluation of quasi-one- and two-dimensional Euler and Navier-Stokes benchmark test problems.

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.

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%

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

International Nuclear Information System (INIS)

Moczo, P.; Kristek, J.; Pazak, P.; Balazovjech, M.; Moczo, P.; Kristek, J.; Galis, M.

2007-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 (FD), finite-element (FE), and hybrid FD-FE 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. We present alternative formulations of equation of motion for a smooth elastic continuum. We then develop alternative formulations for a canonical problem with a welded material interface and free surface. We continue with a model of an earthquake source. We complete the general theoretical introduction by a chapter on the constitutive laws for elastic and viscoelastic media, and brief review of strong formulations of the equation of motion. What follows is a block of chapters on the finite-difference and finite-element methods. We develop FD targets for the free surface and welded material interface. We then present various FD schemes for a smooth continuum, free surface, and welded interface. We focus on the staggered-grid and mainly optimally-accurate FD schemes. We also present alternative formulations of the FE method. We include the FD and FE implementations of the traction-at-split-nodes method for simulation of dynamic rupture propagation. The FD modeling is applied to the model of the deep sedimentary Grenoble basin, France. The FD and FE methods are combined in the hybrid FD-FE method. The hybrid

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

Overview of adaptive finite element analysis in computational geodynamics

Science.gov (United States)

May, D. A.; Schellart, W. P.; Moresi, L.

2013-10-01

The use of numerical models to develop insight and intuition into the dynamics of the Earth over geological time scales is a firmly established practice in the geodynamics community. As our depth of understanding grows, and hand-in-hand with improvements in analytical techniques and higher resolution remote sensing of the physical structure and state of the Earth, there is a continual need to develop more efficient, accurate and reliable numerical techniques. This is necessary to ensure that we can meet the challenge of generating robust conclusions, interpretations and predictions from improved observations. In adaptive numerical methods, the desire is generally to maximise the quality of the numerical solution for a given amount of computational effort. Neither of these terms has a unique, universal definition, but typically there is a trade off between the number of unknowns we can calculate to obtain a more accurate representation of the Earth, and the resources (time and computational memory) required to compute them. In the engineering community, this topic has been extensively examined using the adaptive finite element (AFE) method. Recently, the applicability of this technique to geodynamic processes has started to be explored. In this review we report on the current status and usage of spatially adaptive finite element analysis in the field of geodynamics. The objective of this review is to provide a brief introduction to the area of spatially adaptive finite analysis, including a summary of different techniques to define spatial adaptation and of different approaches to guide the adaptive process in order to control the discretisation error inherent within the numerical solution. An overview of the current state of the art in adaptive modelling in geodynamics is provided, together with a discussion pertaining to the issues related to using adaptive analysis techniques and perspectives for future research in this area. Additionally, we also provide a

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.

Soft tissue deformation estimation by spatio-temporal Kalman filter finite element method.

Science.gov (United States)

Yarahmadian, Mehran; Zhong, Yongmin; Gu, Chengfan; Shin, Jaehyun

2018-01-01

Soft tissue modeling plays an important role in the development of surgical training simulators as well as in robot-assisted minimally invasive surgeries. It has been known that while the traditional Finite Element Method (FEM) promises the accurate modeling of soft tissue deformation, it still suffers from a slow computational process. This paper presents a Kalman filter finite element method to model soft tissue deformation in real time without sacrificing the traditional FEM accuracy. The proposed method employs the FEM equilibrium equation and formulates it as a filtering process to estimate soft tissue behavior using real-time measurement data. The model is temporally discretized using the Newmark method and further formulated as the system state equation. Simulation results demonstrate that the computational time of KF-FEM is approximately 10 times shorter than the traditional FEM and it is still as accurate as the traditional FEM. The normalized root-mean-square error of the proposed KF-FEM in reference to the traditional FEM is computed as 0.0116. It is concluded that the proposed method significantly improves the computational performance of the traditional FEM without sacrificing FEM accuracy. The proposed method also filters noises involved in system state and measurement data.

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

From Laser Scanning to Finite Element Analysis of Complex Buildings by Using a Semi-Automatic Procedure.

Science.gov (United States)

Castellazzi, Giovanni; D'Altri, Antonio Maria; Bitelli, Gabriele; Selvaggi, Ilenia; Lambertini, Alessandro

2015-07-28

In this paper, a new semi-automatic procedure to transform three-dimensional point clouds of complex objects to three-dimensional finite element models is presented and validated. The procedure conceives of the point cloud as a stacking of point sections. The complexity of the clouds is arbitrary, since the procedure is designed for terrestrial laser scanner surveys applied to buildings with irregular geometry, such as historical buildings. The procedure aims at solving the problems connected to the generation of finite element models of these complex structures by constructing a fine discretized geometry with a reduced amount of time and ready to be used with structural analysis. If the starting clouds represent the inner and outer surfaces of the structure, the resulting finite element model will accurately capture the whole three-dimensional structure, producing a complex solid made by voxel elements. A comparison analysis with a CAD-based model is carried out on a historical building damaged by a seismic event. The results indicate that the proposed procedure is effective and obtains comparable models in a shorter time, with an increased level of automation.

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

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2015-09-01

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

Electric field calculations in brain stimulation based on finite elements

DEFF Research Database (Denmark)

Windhoff, Mirko; Opitz, Alexander; Thielscher, Axel

2013-01-01

The need for realistic electric field calculations in human noninvasive brain stimulation is undisputed to more accurately determine the affected brain areas. However, using numerical techniques such as the finite element method (FEM) is methodologically complex, starting with the creation...... of accurate head models to the integration of the models in the numerical calculations. These problems substantially limit a more widespread application of numerical methods in brain stimulation up to now. We introduce an optimized processing pipeline allowing for the automatic generation of individualized...... the successful usage of the pipeline in six subjects, including field calculations for transcranial magnetic stimulation and transcranial direct current stimulation. The quality of the head volume meshes is validated both in terms of capturing the underlying anatomy and of the well-shapedness of the mesh...

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.

Three-dimensional parallel edge-based finite element modeling of electromagnetic data with field redatuming

DEFF Research Database (Denmark)

Cai, Hongzhu; Čuma, Martin; Zhdanov, Michael

2015-01-01

This paper presents a parallelized version of the edge-based finite element method with a novel post-processing approach for numerical modeling of an electromagnetic field in complex media. The method uses an unstructured tetrahedral mesh which can reduce the number of degrees of freedom signific......This paper presents a parallelized version of the edge-based finite element method with a novel post-processing approach for numerical modeling of an electromagnetic field in complex media. The method uses an unstructured tetrahedral mesh which can reduce the number of degrees of freedom...... significantly. The linear system of finite element equations is solved using parallel direct solvers which are robust for ill-conditioned systems and efficient for multiple source electromagnetic (EM) modeling. We also introduce a novel approach to compute the scalar components of the electric field from...... the tangential components along each edge based on field redatuming. The method can produce a more accurate result as compared to conventional approach. We have applied the developed algorithm to compute the EM response for a typical 3D anisotropic geoelectrical model of the off-shore HC reservoir with complex...

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.

A piecewise bi-linear discontinuous finite element spatial discretization of the Sn transport equation

International Nuclear Information System (INIS)

Bailey, Teresa S.; Warsa, James S.; Chang, Jae H.; Adams, Marvin L.

2011-01-01

We present a new spatial discretization of the discrete-ordinates transport equation in two dimensional Cartesian (X-Y) geometry for arbitrary polygonal meshes. The discretization is a discontinuous finite element method (DFEM) that utilizes piecewise bi-linear (PWBL) basis functions, which are formally introduced in this paper. We also present a series of numerical results on quadrilateral and polygonal grids and compare these results to a variety of other spatial discretization that have been shown to be successful on these grid types. Finally, we note that the properties of the PWBL basis functions are such that the leading-order piecewise bi-linear discontinuous finite element (PWBLD) solution will satisfy a reasonably accurate diffusion discretization in the thick diffusion limit, making the PWBLD method a viable candidate for many different classes of transport problems. (author)

A Piecewise Bi-Linear Discontinuous Finite Element Spatial Discretization of the Sn Transport Equation

International Nuclear Information System (INIS)

Bailey, T.S.; Chang, J.H.; Warsa, J.S.; Adams, M.L.

2010-01-01

We present a new spatial discretization of the discrete-ordinates transport equation in two-dimensional Cartesian (X-Y) geometry for arbitrary polygonal meshes. The discretization is a discontinuous finite element method (DFEM) that utilizes piecewise bi-linear (PWBL) basis functions, which are formally introduced in this paper. We also present a series of numerical results on quadrilateral and polygonal grids and compare these results to a variety of other spatial discretizations that have been shown to be successful on these grid types. Finally, we note that the properties of the PWBL basis functions are such that the leading-order piecewise bi-linear discontinuous finite element (PWBLD) solution will satisfy a reasonably accurate diffusion discretization in the thick diffusion limit, making the PWBLD method a viable candidate for many different classes of transport problems.

A Piecewise Bi-Linear Discontinuous Finite Element Spatial Discretization of the Sn Transport Equation

Energy Technology Data Exchange (ETDEWEB)

Bailey, T S; Chang, J H; Warsa, J S; Adams, M L

2010-12-22

We present a new spatial discretization of the discrete-ordinates transport equation in two-dimensional Cartesian (X-Y) geometry for arbitrary polygonal meshes. The discretization is a discontinuous finite element method (DFEM) that utilizes piecewise bi-linear (PWBL) basis functions, which are formally introduced in this paper. We also present a series of numerical results on quadrilateral and polygonal grids and compare these results to a variety of other spatial discretizations that have been shown to be successful on these grid types. Finally, we note that the properties of the PWBL basis functions are such that the leading-order piecewise bi-linear discontinuous finite element (PWBLD) solution will satisfy a reasonably accurate diffusion discretization in the thick diffusion limit, making the PWBLD method a viable candidate for many different classes of transport problems.

Finite element modeling and experimentation of bone drilling forces

International Nuclear Information System (INIS)

Lughmani, W A; Bouazza-Marouf, K; Ashcroft, I

2013-01-01

Bone drilling is an essential part of many orthopaedic surgery procedures, including those for internal fixation and for attaching prosthetics. Estimation and control of bone drilling forces are critical to prevent drill breakthrough, excessive heat generation, and mechanical damage to the bone. This paper presents a 3D finite element (FE) model for prediction of thrust forces experienced during bone drilling. The model incorporates the dynamic characteristics involved in the process along with the accurate geometrical considerations. The average critical thrust forces and torques obtained using FE analysis, for set of machining parameters are found to be in good agreement with the experimental results

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

Generalized multiscale finite element methods (GMsFEM)

KAUST Repository

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

2013-01-01

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

Generalized multiscale finite element methods (GMsFEM)

KAUST Repository

Efendiev, Yalchin R.

2013-10-01

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

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

International Nuclear Information System (INIS)

Johnson, K.; Bittorf, K.J.

2002-01-01

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

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

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.

Direct Calculation of Permeability by High-Accurate Finite Difference and Numerical Integration Methods

KAUST Repository

Wang, Yi

2016-07-21

Velocity of fluid flow in underground porous media is 6~12 orders of magnitudes lower than that in pipelines. If numerical errors are not carefully controlled in this kind of simulations, high distortion of the final results may occur [1-4]. To fit the high accuracy demands of fluid flow simulations in porous media, traditional finite difference methods and numerical integration methods are discussed and corresponding high-accurate methods are developed. When applied to the direct calculation of full-tensor permeability for underground flow, the high-accurate finite difference method is confirmed to have numerical error as low as 10-5% while the high-accurate numerical integration method has numerical error around 0%. Thus, the approach combining the high-accurate finite difference and numerical integration methods is a reliable way to efficiently determine the characteristics of general full-tensor permeability such as maximum and minimum permeability components, principal direction and anisotropic ratio. Copyright © Global-Science Press 2016.

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.

Neurosurgery simulation using non-linear finite element modeling and haptic interaction

Science.gov (United States)

Lee, Huai-Ping; Audette, Michel; Joldes, Grand R.; Enquobahrie, Andinet

2012-02-01

Real-time surgical simulation is becoming an important component of surgical training. To meet the realtime requirement, however, the accuracy of the biomechancial modeling of soft tissue is often compromised due to computing resource constraints. Furthermore, haptic integration presents an additional challenge with its requirement for a high update rate. As a result, most real-time surgical simulation systems employ a linear elasticity model, simplified numerical methods such as the boundary element method or spring-particle systems, and coarse volumetric meshes. However, these systems are not clinically realistic. We present here an ongoing work aimed at developing an efficient and physically realistic neurosurgery simulator using a non-linear finite element method (FEM) with haptic interaction. Real-time finite element analysis is achieved by utilizing the total Lagrangian explicit dynamic (TLED) formulation and GPU acceleration of per-node and per-element operations. We employ a virtual coupling method for separating deformable body simulation and collision detection from haptic rendering, which needs to be updated at a much higher rate than the visual simulation. The system provides accurate biomechancial modeling of soft tissue while retaining a real-time performance with haptic interaction. However, our experiments showed that the stability of the simulator depends heavily on the material property of the tissue and the speed of colliding objects. Hence, additional efforts including dynamic relaxation are required to improve the stability of the system.

Assessment of natural frequency of installed offshore wind turbines using nonlinear finite element model considering soil-monopile interaction

Directory of Open Access Journals (Sweden)

Djillali Amar Bouzid

2018-04-01

Full Text Available A nonlinear finite element model is developed to examine the lateral behaviors of monopiles, which support offshore wind turbines (OWTs chosen from five different offshore wind farms in Europe. The simulation is using this model to accurately estimate the natural frequency of these slender structures, as a function of the interaction of the foundations with the subsoil. After a brief introduction to the wind power energy as a reliable alternative in comparison to fossil fuel, the paper focuses on concept of natural frequency as a primary indicator in designing the foundations of OWTs. Then the range of natural frequencies is provided for a safe design purpose. Next, an analytical expression of an OWT natural frequency is presented as a function of soil-monopile interaction through monopile head springs characterized by lateral stiffness KL, rotational stiffness KR and cross-coupling stiffness KLR, of which the differences are discussed. The nonlinear pseudo three-dimensional finite element vertical slices model has been used to analyze the lateral behaviors of monopiles supporting the OWTs of different wind farm sites considered. Through the monopiles head movements (displacements and rotations, the values of KL, KR and KLR were obtained and substituted in the analytical expression of natural frequency for comparison. The comparison results between computed and measured natural frequencies showed an excellent agreement for most cases. This confirms the convenience of the finite element model used for the accurate estimation of the monopile head stiffness. Keywords: Nonlinear finite element analysis, Vertical slices model, Monopiles under horizontal loading, Natural frequency, Monopile head stiffness, Offshore wind turbines (OWTs

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

A Posteriori Error Estimation for Finite Element Methods and Iterative Linear Solvers

Energy Technology Data Exchange (ETDEWEB)

Melboe, Hallgeir

2001-10-01

This thesis addresses a posteriori error estimation for finite element methods and iterative linear solvers. Adaptive finite element methods have gained a lot of popularity over the last decades due to their ability to produce accurate results with limited computer power. In these methods a posteriori error estimates play an essential role. Not only do they give information about how large the total error is, they also indicate which parts of the computational domain should be given a more sophisticated treatment in order to reduce the error. A posteriori error estimates are traditionally aimed at estimating the global error, but more recently so called goal oriented error estimators have been shown a lot of interest. The name reflects the fact that they estimate the error in user-defined local quantities. In this thesis the main focus is on global error estimators for highly stretched grids and goal oriented error estimators for flow problems on regular grids. Numerical methods for partial differential equations, such as finite element methods and other similar techniques, typically result in a linear system of equations that needs to be solved. Usually such systems are solved using some iterative procedure which due to a finite number of iterations introduces an additional error. Most such algorithms apply the residual in the stopping criterion, whereas the control of the actual error may be rather poor. A secondary focus in this thesis is on estimating the errors that are introduced during this last part of the solution procedure. The thesis contains new theoretical results regarding the behaviour of some well known, and a few new, a posteriori error estimators for finite element methods on anisotropic grids. Further, a goal oriented strategy for the computation of forces in flow problems is devised and investigated. Finally, an approach for estimating the actual errors associated with the iterative solution of linear systems of equations is suggested. (author)

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.

Highly accurate symplectic element based on two variational principles

Science.gov (United States)

Qing, Guanghui; Tian, Jia

2018-02-01

For the stability requirement of numerical resultants, the mathematical theory of classical mixed methods are relatively complex. However, generalized mixed methods are automatically stable, and their building process is simple and straightforward. In this paper, based on the seminal idea of the generalized mixed methods, a simple, stable, and highly accurate 8-node noncompatible symplectic element (NCSE8) was developed by the combination of the modified Hellinger-Reissner mixed variational principle and the minimum energy principle. To ensure the accuracy of in-plane stress results, a simultaneous equation approach was also suggested. Numerical experimentation shows that the accuracy of stress results of NCSE8 are nearly the same as that of displacement methods, and they are in good agreement with the exact solutions when the mesh is relatively fine. NCSE8 has advantages of the clearing concept, easy calculation by a finite element computer program, higher accuracy and wide applicability for various linear elasticity compressible and nearly incompressible material problems. It is possible that NCSE8 becomes even more advantageous for the fracture problems due to its better accuracy of stresses.

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.

Finite element Fourier and Abbe transform methods for generalization of aperture function and geometry in Fraunhofer diffraction theory

International Nuclear Information System (INIS)

Kraus, H.G.

1991-01-01

This paper discusses methods for calculating Fraunhofer intensity fields resulting from diffraction through one- and two-dimensional apertures are presented. These methods are based on the geometric concept of finite elements and on Fourier and Abbe transforms. The geometry of the two-dimensional diffracting aperture(s) is based on biquadratic isoparametric elements, which are used to define aperture(s) of complex geometry. These elements are also used to build complex amplitude and phase functions across the aperture(s) which may be of continuous or discontinuous form. The transform integrals are accurately and efficiently integrated numerically using Gaussian quadrature. The power of these methods is most evident in two dimensions, where several examples are presented which include secondary obstructions, straight and curved secondary spider supports, multiple-mirror arrays, synthetic aperture arrays, segmented mirrors, apertures covered by screens, apodization, and phase plates. Typically, the finite element Abbe transform method results in significant gains in computational efficiency over the finite element Fourier transform method, but is also subject to some loss in generality

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.

A stabilized second-order time accurate finite element formulation for incompressible viscous flow with heat transfer; Uma formulacao de elementos finitos estabilizada de segunda ordem no tempo para escoamentos viscosos com transferencia de calor

Energy Technology Data Exchange (ETDEWEB)

Curi, Marcos Filardy

2011-07-01

In view of the problem of global warming and the search for clean energy sources, a worldwide expansion on the use of nuclear energy is foreseen. Thus, the development of science and technology regarding nuclear power plants is essential, in particular in the field of reactor engineering. Fluid mechanics and heat transfer play an important role in the development of nuclear reactors. Computational Fluid Mechanics (CFD) is becoming ever more important in the optimization of cost and safety of the designs. This work presents a stabilized second-order time accurate finite element formulation for incompressible flows with heat transfer. A second order time discretization precedes a spatial discretization using finite elements. The terms that stabilize the finite element method arise naturally from the discretization process, rather than being introduced a priori in the variational formulation. The method was implemented in the program 'ns{sub n}ew{sub s}olvec2d{sub av}2{sub M}PI' written in FORTRAN90, developed in the Parallel Computing Laboratory at the Institute of Nuclear Engineering (LCP/IEN). Numerical solutions of some representative examples, including free, mixed and forced convection, demonstrate that the proposed stabilized formulation attains very good agreement with experimental and computational results available in the literature. (author)

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)

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

Directory of Open Access Journals (Sweden)

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

2016-11-01

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

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)

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.

Development of Graphical User Interface for Finite Element Analysis of Static Loading of a Column using MATLAB

Directory of Open Access Journals (Sweden)

Moses Omolayo PETINRIN

2010-12-01

Full Text Available In this work, the capability of MATLAB software package to develop graphical user interface (GUI package was demonstrated. A GUI was successfully developed using MATLAB programming language to study the behaviour of a suspended column under uniaxial static loading by solving the numerical model created based on the finite element method (FEM. The comparison between the exact solution from previous researches and the numerical analysis showed good agreement. The column average strain, average stress and average load are equivalent but more accurate to the ones obtained when the whole column is taken as one element (two nodes for one dimensional linear finite element problem. It was established in this work that MATLAB is not only a software package for numerical computation but also for application development.

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.

Analysis of the jet pipe electro-hydraulic servo valve with finite element methods

Directory of Open Access Journals (Sweden)

Kaiyu Zhao

2018-01-01

Full Text Available The dynamic characteristics analysis about the jet pipe electro-hydraulic servo valve based on experience and mathematical derivation was difficult and not so precise. So we have analysed the armature feedback components, torque motor and jet pipe receiver in electrohydraulic servo valve by sophisticated finite element analysis tools respectively and have got physical meaning data on these parts. Then the data were fitted by Matlab and the mathematical relationships among them were calculated. We have done the dynamic multi-physical fields’ Simulink co-simulation using above mathematical relationship, and have got the input-output relationship of the overall valve, the frequency response and step response. This work can show the actual working condition accurately. At the same time, we have considered the materials and the impact of the critical design dimensions in the finite element analysis process. It provides some new ideas to the overall design of jet pipe electro-hydraulic servo valve.

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)

Hybrid determination of mixed-mode stress intensity factors on discontinuous finite-width plate by finite element and photoelasticity

International Nuclear Information System (INIS)

Baek, Tae Hyun; Chen, Lei; Hong, Dong Pyo

2011-01-01

For isotropic material structure, the stress in the vicinity of crack tip is generally much higher than the stress far away from it. This phenomenon usually leads to stress concentration and fracture of structure. Previous researches and studies show that the stress intensity factor is one of most important parameter for crack growth and propagation. This paper provides a convenient numerical method, which is called hybrid photoelasticity method, to accurately determine the stress field distribution in the vicinity of crack tip and mixed-mode stress intensity factors. The model was simulated by finite element method and isochromatic data along straight lines far away from the crack tip were calculated. By using the isochromatic data obtained from finite element method and a conformal mapping procedure, stress components and photoelastic fringes in the hybrid region were calculated. To easily compare calculated photoelastic fringes with experiment results, the fringe patterns were reconstructed, doubled and sharpened. Good agreement shows that the method presented in this paper is reliable and convenient. This method can then directly be applied to obtain mixed mode stress intensity factors from the experimentally measured isochromatic data along the straight lines

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)

Matrix equation decomposition and parallel solution of systems resulting from unstructured finite element problems in electromagnetics

Energy Technology Data Exchange (ETDEWEB)

Cwik, T. [California Institute of Technology, Pasadena, CA (United States); Katz, D.S. [Cray Research, El Segundo, CA (United States)

1996-12-31

Finite element modeling has proven useful for accurately simulating scattered or radiated electromagnetic fields from complex three-dimensional objects whose geometry varies on the scale of a fraction of an electrical wavelength. An unstructured finite element model of realistic objects leads to a large, sparse, system of equations that needs to be solved efficiently with regard to machine memory and execution time. Both factorization and iterative solvers can be used to produce solutions to these systems of equations. Factorization leads to high memory requirements that limit the electrical problem size of three-dimensional objects that can be modeled. An iterative solver can be used to efficiently solve the system without excessive memory use and in a minimal amount of time if the convergence rate is controlled.

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

Science.gov (United States)

Nakashima, Hiroshi; Takatsu, Yuzuru

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

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 three-dimensional finite element model for biomechanical analysis of the hip.

Science.gov (United States)

Chen, Guang-Xing; Yang, Liu; Li, Kai; He, Rui; Yang, Bin; Zhan, Yan; Wang, Zhi-Jun; Yu, Bing-Nin; Jian, Zhe

2013-11-01

The objective of this study was to construct a three-dimensional (3D) finite element model of the hip. The images of the hip were obtained from Chinese visible human dataset. The hip model includes acetabular bone, cartilage, labrum, and bone. The cartilage of femoral head was constructed using the AutoCAD and Solidworks software. The hip model was imported into ABAQUS analysis system. The contact surface of the hip joint was meshed. To verify the model, the single leg peak force was loaded, and contact area of the cartilage and labrum of the hip and pressure distribution in these structures were observed. The constructed 3D hip model reflected the real hip anatomy. Further, this model reflected biomechanical behavior similar to previous studies. In conclusion, this 3D finite element hip model avoids the disadvantages of other construction methods, such as imprecision of cartilage construction and the absence of labrum. Further, it provides basic data critical for accurately modeling normal and abnormal loads, and the effects of abnormal loads on the hip.

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

CSIR Research Space (South Africa)

Suliman, Ridhwaan

2012-07-01

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

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

CERN Document Server

Öchsner, Andreas

2016-01-01

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

A finite element thermohydrodynamic analyis of profile bore bearing

International Nuclear Information System (INIS)

Shah Nor bin Basri

1994-01-01

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

Hualien forced vibration calculation with a finite element model

International Nuclear Information System (INIS)

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

1995-01-01

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

The Finite Element Numerical Modelling of 3D Magnetotelluric

Directory of Open Access Journals (Sweden)

Ligang Cao

2014-01-01

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

A parallel finite-volume finite-element method for transient compressible turbulent flows with heat transfer

International Nuclear Information System (INIS)

Masoud Ziaei-Rad

2010-01-01

In this paper, a two-dimensional numerical scheme is presented for the simulation of turbulent, viscous, transient compressible flows in the simultaneously developing hydraulic and thermal boundary layer region. The numerical procedure is a finite-volume-based finite-element method applied to unstructured grids. This combination together with a new method applied for the boundary conditions allows for accurate computation of the variables in the entrance region and for a wide range of flow fields from subsonic to transonic. The Roe-Riemann solver is used for the convective terms, whereas the standard Galerkin technique is applied for the viscous terms. A modified κ-ε model with a two-layer equation for the near-wall region combined with a compressibility correction is used to predict the turbulent viscosity. Parallel processing is also employed to divide the computational domain among the different processors to reduce the computational time. The method is applied to some test cases in order to verify the numerical accuracy. The results show significant differences between incompressible and compressible flows in the friction coefficient, Nusselt number, shear stress and the ratio of the compressible turbulent viscosity to the molecular viscosity along the developing region. A transient flow generated after an accidental rupture in a pipeline was also studied as a test case. The results show that the present numerical scheme is stable, accurate and efficient enough to solve the problem of transient wall-bounded flow.

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.

Calculation of phonon dispersion in carbon nanotubes using a continuum-atomistic finite element approach

Directory of Open Access Journals (Sweden)

Michael J. Leamy

2011-12-01

Full Text Available Dispersion calculations are presented for cylindrical carbon nanotubes using a manifold-based continuum-atomistic finite element formulation combined with Bloch analysis. The formulated finite elements allow any (n,m chiral nanotube, or mixed tubes formed by periodically-repeating heterojunctions, to be examined quickly and accurately using only three input parameters (radius, chiral angle, and unit cell length and a trivial structured mesh, thus avoiding the tedious geometry generation and energy minimization tasks associated with ab initio and lattice dynamics-based techniques. A critical assessment of the technique is pursued to determine the validity range of the resulting dispersion calculations, and to identify any dispersion anomalies. Two small anomalies in the dispersion curves are documented, which can be easily identified and therefore rectified. They include difficulty in achieving a zero energy point for the acoustic twisting phonon, and a branch veering in nanotubes with nonzero chiral angle. The twisting mode quickly restores its correct group velocity as wavenumber increases, while the branch veering is associated with a rapid exchange of eigenvectors at the veering point, which also lessens its impact. By taking into account the two noted anomalies, accurate predictions of acoustic and low-frequency optical branches can be achieved out to the midpoint of the first Brillouin zone.

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

Finite Element Modeling of Thermo Creep Processes Using Runge-Kutta Method

Directory of Open Access Journals (Sweden)

Yu. I. Dimitrienko

2015-01-01

Full Text Available Thermo creep deformations for most heat-resistant alloys, as a rule, nonlinearly depend on stresses and are practically non- reversible. Therefore, to calculate the properties of these materials the theory of plastic flow is most widely used. Finite-element computations of a stress-strain state of structures with account of thermo creep deformations up to now are performed using main commercial software, including ANSYS package. However, in most cases to solve nonlinear creep equations, one should apply explicit or implicit methods based on the Euler method of approximation of time-derivatives. The Euler method is sufficiently efficient in terms of random access memory in computations, however this method is cumbersome in computation time and does not always provide a required accuracy for creep deformation computations.The paper offers a finite-element algorithm to solve a three-dimensional problem of thermo creep based on the Runge-Kutta finite-difference schemes of different orders with respect to time. It shows a numerical test example to solve the problem on the thermo creep of a beam under tensile loading. The computed results demonstrate that using the Runge-Kutta method with increasing accuracy order allows us to obtain a more accurate solution (with increasing accuracy order by 1 a relative error decreases, approximately, by an order too. The developed algorithm proves to be efficient enough and can be recommended for solving the more complicated problems of thermo creep of structures.

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2012-09-20

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

Analysis of uncertainty and variability in finite element computational models for biomedical engineering:characterization and propagation

Directory of Open Access Journals (Sweden)

Nerea Mangado

2016-11-01

Full Text Available Computational modeling has become a powerful tool in biomedical engineering thanks to its potential to simulate coupled systems. However, real parameters are usually not accurately known and variability is inherent in living organisms. To cope with this, probabilistic tools, statistical analysis and stochastic approaches have been used. This article aims to review the analysis of uncertainty and variability in the context of finite element modeling in biomedical engineering. Characterization techniques and propagation methods are presented, as well as examples of their applications in biomedical finite element simulations. Uncertainty propagation methods, both non-intrusive and intrusive, are described. Finally, pros and cons of the different approaches and their use in the scientific community are presented. This leads us to identify future directions for research and methodological development of uncertainty modeling in biomedical engineering.

Analysis of Uncertainty and Variability in Finite Element Computational Models for Biomedical Engineering: Characterization and Propagation.

Science.gov (United States)

Mangado, Nerea; Piella, Gemma; Noailly, Jérôme; Pons-Prats, Jordi; Ballester, Miguel Ángel González

2016-01-01

Computational modeling has become a powerful tool in biomedical engineering thanks to its potential to simulate coupled systems. However, real parameters are usually not accurately known, and variability is inherent in living organisms. To cope with this, probabilistic tools, statistical analysis and stochastic approaches have been used. This article aims to review the analysis of uncertainty and variability in the context of finite element modeling in biomedical engineering. Characterization techniques and propagation methods are presented, as well as examples of their applications in biomedical finite element simulations. Uncertainty propagation methods, both non-intrusive and intrusive, are described. Finally, pros and cons of the different approaches and their use in the scientific community are presented. This leads us to identify future directions for research and methodological development of uncertainty modeling in biomedical engineering.

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

Engineering computation of structures the finite element method

CERN Document Server

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

2015-01-01

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

Finite Element Method Application in Areal Rainfall Estimation Case Study; Mashhad Plain Basin

Directory of Open Access Journals (Sweden)

M. Irani

2016-10-01

Full Text Available Introduction: The hydrological models are very important tools for planning and management of water resources. These models can be used for identifying basin and nature problems and choosing various managements. Precipitation is based on these models. Calculations of rainfall would be affected by displacement and region factor such as topography, etc. Estimating areal rainfall is one of the basic needs in meteorological, water resources and others studies. There are various methods for the estimation of rainfall, which can be evaluated by using statistical data and mathematical terms. In hydrological analysis, areal rainfall is so important because of displacement of precipitation. Estimating areal rainfall is divided to three methods: 1- graphical. 2-topographical. 3-numerical. This paper represented calculating mean precipitation (daily, monthly and annual using Galerkin’s method (numerical method and it was compared with other methods such as kriging, IDW, Thiessen and arithmetic mean. In this study, there were 42 actual gauges and thirteen dummies in Mashhad plain basin which is calculated by Galerkin’s method. The method included the use of interpolation functions, allowing an accurate representation of shape and relief of catchment with numerical integration performed by Gaussian quadrature and represented the allocation of weights to stations. Materials and Methods:The estimation of areal rainfall (daily, monthly,… is the basic need for meteorological project. In this field ,there are various methods that one of them is finite element method. Present study aimed to estimate areal rainfall with a 16-year period (1997-2012 by using Galerkin method ( finite element in Mashhad plain basin for 42 station. Therefore, it was compared with other usual methods such as arithmetic mean, Thiessen, Kriging and IDW. The analysis of Thiessen, Kriging and IDW were in ArcGIS10.0 software environment and finite element analysis did by using of Matlab

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

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

The nonlinear finite element analysis program NUCAS (NUclear Containment Analysis System) for reinforced concrete containment building

Energy Technology Data Exchange (ETDEWEB)

Lee, Sang Jin; Lee, Hong Pyo; Seo, Jeong Moon [Korea Atomic Energy Research Institute, Taejeon (Korea)

2002-03-01

The maim goal of this research is to develop a nonlinear finite element analysis program NUCAS to accurately predict global and local failure modes of containment building subjected to internal pressure. In this report, we describe the techniques we developed throught this research. An adequate model to the analysis of containment building such as microscopic material model is adopted and it applied into the development Reissner-Mindlin degenerated shell element. To avoid finite element deficiencies, the substitute strains based on the assumed strain method is used in the shell formulation. Arc-length control method is also adopted to fully trace the peak load-displacement path due to crack formation. In addition, a benchmark test suite is developed to investigate the performance of NUCAS and proposed as the future benchmark tests for nonlinear analysis of reinforced concrete. Finally, the input format of NUCAS and the examples of input/output file are described. 39 refs., 65 figs., 8 tabs. (Author)

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.

A vortex model for Darrieus turbine using finite element techniques

Energy Technology Data Exchange (ETDEWEB)

Ponta, Fernando L. [Universidad de Buenos Aires, Dept. de Electrotecnia, Grupo ISEP, Buenos Aires (Argentina); Jacovkis, Pablo M. [Universidad de Buenos Aires, Dept. de Computacion and Inst. de Calculo, Buenos Aires (Argentina)

2001-09-01

Since 1970 several aerodynamic prediction models have been formulated for the Darrieus turbine. We can identify two families of models: stream-tube and vortex. The former needs much less computation time but the latter is more accurate. The purpose of this paper is to show a new option for modelling the aerodynamic behaviour of Darrieus turbines. The idea is to combine a classic free vortex model with a finite element analysis of the flow in the surroundings of the blades. This avoids some of the remaining deficiencies in classic vortex models. The agreement between analysis and experiment when predicting instantaneous blade forces and near wake flow behind the rotor is better than the one obtained in previous models. (Author)

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

African Journals Online (AJOL)

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

A finite element for plates and shells

International Nuclear Information System (INIS)

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

1981-08-01

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

COMPUTER EXPERIMENTS WITH FINITE ELEMENTS OF HIGHER ORDER

Directory of Open Access Journals (Sweden)

Khomchenko A.

2017-12-01

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

Fourier analysis of finite element preconditioned collocation schemes

Science.gov (United States)

Deville, Michel O.; Mund, Ernest H.

1990-01-01

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

Stochastic Finite Elements in Reliability-Based Structural Optimization

DEFF Research Database (Denmark)

Sørensen, John Dalsgaard; Engelund, S.

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

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

Energy Technology Data Exchange (ETDEWEB)

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

2014-10-15

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

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

The accuracy of the Gaussian-and-finite-element-Coulomb (GFC) method for the calculation of Coulomb integrals.

Science.gov (United States)

Przybytek, Michal; Helgaker, Trygve

2013-08-07

We analyze the accuracy of the Coulomb energy calculated using the Gaussian-and-finite-element-Coulomb (GFC) method. In this approach, the electrostatic potential associated with the molecular electronic density is obtained by solving the Poisson equation and then used to calculate matrix elements of the Coulomb operator. The molecular electrostatic potential is expanded in a mixed Gaussian-finite-element (GF) basis set consisting of Gaussian functions of s symmetry centered on the nuclei (with exponents obtained from a full optimization of the atomic potentials generated by the atomic densities from symmetry-averaged restricted open-shell Hartree-Fock theory) and shape functions defined on uniform finite elements. The quality of the GF basis is controlled by means of a small set of parameters; for a given width of the finite elements d, the highest accuracy is achieved at smallest computational cost when tricubic (n = 3) elements are used in combination with two (γ(H) = 2) and eight (γ(1st) = 8) Gaussians on hydrogen and first-row atoms, respectively, with exponents greater than a given threshold (αmin (G)=0.5). The error in the calculated Coulomb energy divided by the number of atoms in the system depends on the system type but is independent of the system size or the orbital basis set, vanishing approximately like d(4) with decreasing d. If the boundary conditions for the Poisson equation are calculated in an approximate way, the GFC method may lose its variational character when the finite elements are too small; with larger elements, it is less sensitive to inaccuracies in the boundary values. As it is possible to obtain accurate boundary conditions in linear time, the overall scaling of the GFC method for large systems is governed by another computational step-namely, the generation of the three-center overlap integrals with three Gaussian orbitals. The most unfavorable (nearly quadratic) scaling is observed for compact, truly three-dimensional systems

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.

Field Strain Measurement on the Fiber Scale in Carbon Fiber Reinforced Polymers Using Global Finite-Element Based Digital Image Correlation

KAUST Repository

Tao, Ran

2015-01-01

is aimed to accurately measure the displacement and strain fields at the fiber-matrix scale in a cross-ply composite. First, the theories of both local subset-based digital image correlation (DIC) and global finite-element based DIC are outlined. Second, in

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

Comparison of Finite Element Predictions to Measurements from the Sandia Microslip Experiment

Energy Technology Data Exchange (ETDEWEB)

LOBITZ,DONALD W.; GREGORY,DANNY LYNN; SMALLWOOD,DAVID O.

2000-11-09

When embarking on an experimental program for purposes of discovery and understanding, it is only prudent to use appropriate analysis tools to aid in the discovery process. Due to the limited scope of experimental measurement analytical results can significantly complement the data after a reasonable validation process has occurred. In this manner the analytical results can help to explain certain measurements, suggest other measurements to take and point to possible modifications to the experimental apparatus. For these reasons it was decided to create a detailed nonlinear finite element model of the Sandia Microslip Experiment. This experiment was designed to investigate energy dissipation due to microslip in bolted joints and to identify the critical parameters involved. In an attempt to limit the microslip to a single interface a complicated system of rollers and cables was devised to clamp the two slipping members together with a prescribed normal load without using a bolt. An oscillatory tangential load is supplied via a shaker. The finite element model includes the clamping device in addition to the sequence of steps taken in setting up the experiment. The interface is modeled using Coulomb friction requiring a modest validation procedure for estimating the coefficient of friction. Analysis results have indicated misalignment problems in the experimental procedure, identified transducer locations for more accurate measurements, predicted complex interface motions including the potential for galling, identified regions where microslip occurs and during which parts of the loading cycle it occurs, all this in addition to the energy dissipated per cycle. A number of these predictions have been experimentally corroborated in varying degrees and are presented in the paper along with the details of the finite element model.

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

Measurement and Finite Element Model Validation of Immature Porcine Brain-Skull Displacement during Rapid Sagittal Head Rotations.

Science.gov (United States)

Pasquesi, Stephanie A; Margulies, Susan S

2018-01-01

Computational models are valuable tools for studying tissue-level mechanisms of traumatic brain injury, but to produce more accurate estimates of tissue deformation, these models must be validated against experimental data. In this study, we present in situ measurements of brain-skull displacement in the neonatal piglet head ( n  = 3) at the sagittal midline during six rapid non-impact rotations (two rotations per specimen) with peak angular velocities averaging 51.7 ± 1.4 rad/s. Marks on the sagittally cut brain and skull/rigid potting surfaces were tracked, and peak values of relative brain-skull displacement were extracted and found to be significantly less than values extracted from a previous axial plane model. In a finite element model of the sagittally transected neonatal porcine head, the brain-skull boundary condition was matched to the measured physical experiment data. Despite smaller sagittal plane displacements at the brain-skull boundary, the corresponding finite element boundary condition optimized for sagittal plane rotations is far less stiff than its axial counterpart, likely due to the prominent role of the boundary geometry in restricting interface movement. Finally, bridging veins were included in the finite element model. Varying the bridging vein mechanical behavior over a previously reported range had no influence on the brain-skull boundary displacements. This direction-specific sagittal plane boundary condition can be employed in finite element models of rapid sagittal head rotations.

Measurement and Finite Element Model Validation of Immature Porcine Brain–Skull Displacement during Rapid Sagittal Head Rotations

Science.gov (United States)

Pasquesi, Stephanie A.; Margulies, Susan S.

2018-01-01

Computational models are valuable tools for studying tissue-level mechanisms of traumatic brain injury, but to produce more accurate estimates of tissue deformation, these models must be validated against experimental data. In this study, we present in situ measurements of brain–skull displacement in the neonatal piglet head (n = 3) at the sagittal midline during six rapid non-impact rotations (two rotations per specimen) with peak angular velocities averaging 51.7 ± 1.4 rad/s. Marks on the sagittally cut brain and skull/rigid potting surfaces were tracked, and peak values of relative brain–skull displacement were extracted and found to be significantly less than values extracted from a previous axial plane model. In a finite element model of the sagittally transected neonatal porcine head, the brain–skull boundary condition was matched to the measured physical experiment data. Despite smaller sagittal plane displacements at the brain–skull boundary, the corresponding finite element boundary condition optimized for sagittal plane rotations is far less stiff than its axial counterpart, likely due to the prominent role of the boundary geometry in restricting interface movement. Finally, bridging veins were included in the finite element model. Varying the bridging vein mechanical behavior over a previously reported range had no influence on the brain–skull boundary displacements. This direction-specific sagittal plane boundary condition can be employed in finite element models of rapid sagittal head rotations. PMID:29515995

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

Poisson-Nernst-Planck Equations for Simulating Biomolecular Diffusion-Reaction Processes I: Finite Element Solutions.

Science.gov (United States)

Lu, Benzhuo; Holst, Michael J; McCammon, J Andrew; Zhou, Y C

2010-09-20

In this paper we developed accurate finite element methods for solving 3-D Poisson-Nernst-Planck (PNP) equations with singular permanent charges for electrodiffusion in solvated biomolecular systems. The electrostatic Poisson equation was defined in the biomolecules and in the solvent, while the Nernst-Planck equation was defined only in the solvent. We applied a stable regularization scheme to remove the singular component of the electrostatic potential induced by the permanent charges inside biomolecules, and formulated regular, well-posed PNP equations. An inexact-Newton method was used to solve the coupled nonlinear elliptic equations for the steady problems; while an Adams-Bashforth-Crank-Nicolson method was devised for time integration for the unsteady electrodiffusion. We numerically investigated the conditioning of the stiffness matrices for the finite element approximations of the two formulations of the Nernst-Planck equation, and theoretically proved that the transformed formulation is always associated with an ill-conditioned stiffness matrix. We also studied the electroneutrality of the solution and its relation with the boundary conditions on the molecular surface, and concluded that a large net charge concentration is always present near the molecular surface due to the presence of multiple species of charged particles in the solution. The numerical methods are shown to be accurate and stable by various test problems, and are applicable to real large-scale biophysical electrodiffusion problems.

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

Comparison of 3-D finite elements for incompressible fluid flow

International Nuclear Information System (INIS)

Robichaud, M.; Tanguy, P.A.

1985-01-01

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

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

Science.gov (United States)

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

2014-01-01

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

Error-controlled adaptive finite elements in solid mechanics

National Research Council Canada - National Science Library

Stein, Erwin; Ramm, E

2003-01-01

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

Solution of the Lambda modes problem of a nuclear power reactor using an h–p finite element method

International Nuclear Information System (INIS)

Vidal-Ferrandiz, A.; Fayez, R.; Ginestar, D.; Verdú, G.

2014-01-01

Highlights: • An hp finite element method is proposed for the Lambda modes problem of a nuclear reactor. • Different strategies can be implemented for increasing the accuracy of the solutions. • 2D and 3D benchmarks have been studied obtaining accurate results. - Abstract: Lambda modes of a nuclear power reactor have interest in reactor physics since they have been used to develop modal methods and to study BWR reactor instabilities. An h–p-Adaptation finite element method has been implemented to compute the dominant modes the fundamental mode and the next subcritical modes of a nuclear reactor. The performance of this method has been studied in three benchmark problems, a homogeneous 2D reactor, the 2D BIBLIS reactor and the 3D IAEA reactor

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

Comparisons of Particle Tracking Techniques and Galerkin Finite Element Methods in Flow Simulations on Watershed Scales

Science.gov (United States)

Shih, D.; Yeh, G.

2009-12-01

This paper applies two numerical approximations, the particle tracking technique and Galerkin finite element method, to solve the diffusive wave equation in both one-dimensional and two-dimensional flow simulations. The finite element method is one of most commonly approaches in numerical problems. It can obtain accurate solutions, but calculation times may be rather extensive. The particle tracking technique, using either single-velocity or average-velocity tracks to efficiently perform advective transport, could use larger time-step sizes than the finite element method to significantly save computational time. Comparisons of the alternative approximations are examined in this poster. We adapt the model WASH123D to examine the work. WASH123D is an integrated multimedia, multi-processes, physics-based computational model suitable for various spatial-temporal scales, was first developed by Yeh et al., at 1998. The model has evolved in design capability and flexibility, and has been used for model calibrations and validations over the course of many years. In order to deliver a locally hydrological model in Taiwan, the Taiwan Typhoon and Flood Research Institute (TTFRI) is working with Prof. Yeh to develop next version of WASH123D. So, the work of our preliminary cooperationx is also sketched in this poster.

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

Accurate kinematic measurement at interfaces between dissimilar materials using conforming finite-element-based digital image correlation

KAUST Repository

Tao, Ran

2016-02-11

Digital image correlation (DIC) is now an extensively applied full-field measurement technique with subpixel accuracy. A systematic drawback of this technique, however, is the smoothening of the kinematic field (e.g., displacement and strains) across interfaces between dissimilar materials, where the deformation gradient is known to be large. This can become an issue when a high level of accuracy is needed, for example, in the interfacial region of composites or joints. In this work, we described the application of global conforming finite-element-based DIC technique to obtain precise kinematic fields at interfaces between dissimilar materials. Speckle images from both numerical and actual experiments processed by the described global DIC technique better captured sharp strain gradient at the interface than local subset-based DIC. © 2016 Elsevier Ltd. All rights reserved.

Parallelized Three-Dimensional Resistivity Inversion Using Finite Elements And Adjoint State Methods

Science.gov (United States)

Schaa, Ralf; Gross, Lutz; Du Plessis, Jaco

2015-04-01

The resistivity method is one of the oldest geophysical exploration methods, which employs one pair of electrodes to inject current into the ground and one or more pairs of electrodes to measure the electrical potential difference. The potential difference is a non-linear function of the subsurface resistivity distribution described by an elliptic partial differential equation (PDE) of the Poisson type. Inversion of measured potentials solves for the subsurface resistivity represented by PDE coefficients. With increasing advances in multichannel resistivity acquisition systems (systems with more than 60 channels and full waveform recording are now emerging), inversion software require efficient storage and solver algorithms. We developed the finite element solver Escript, which provides a user-friendly programming environment in Python to solve large-scale PDE-based problems (see https://launchpad.net/escript-finley). Using finite elements, highly irregular shaped geology and topography can readily be taken into account. For the 3D resistivity problem, we have implemented the secondary potential approach, where the PDE is decomposed into a primary potential caused by the source current and the secondary potential caused by changes in subsurface resistivity. The primary potential is calculated analytically, and the boundary value problem for the secondary potential is solved using nodal finite elements. This approach removes the singularity caused by the source currents and provides more accurate 3D resistivity models. To solve the inversion problem we apply a 'first optimize then discretize' approach using the quasi-Newton scheme in form of the limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) method (see Gross & Kemp 2013). The evaluation of the cost function requires the solution of the secondary potential PDE for each source current and the solution of the corresponding adjoint-state PDE for the cost function gradients with respect to the subsurface

The effect of loading time on flexible pavement dynamic response: a finite element analysis

Science.gov (United States)

Yin, Hao; Solaimanian, Mansour; Kumar, Tanmay; Stoffels, Shelley

2007-12-01

Dynamic response of asphalt concrete (AC) pavements under moving load is a key component for accurate prediction of flexible pavement performance. The time and temperature dependency of AC materials calls for utilizing advanced material characterization and mechanistic theories, such as viscoelasticity and stress/strain analysis. In layered elastic analysis, as implemented in the new Mechanistic-Empirical Pavement Design Guide (MEPDG), the time dependency is accounted for by calculating the loading times at different AC layer depths. In this study, the time effect on pavement response was evaluated by means of the concept of “pseudo temperature.” With the pavement temperature measured from instrumented thermocouples, the time and temperature dependency of AC materials was integrated into one single factor, termed “effective temperature.” Via this effective temperature, pavement responses under a transient load were predicted through finite element analysis. In the finite element model, viscoelastic behavior of AC materials was characterized through relaxation moduli, while the layers with unbound granular material were assumed to be in an elastic mode. The analysis was conducted for two different AC mixtures in a simplified flexible pavement structure at two different seasons. Finite element analysis results reveal that the loading time has a more pronounced impact on pavement response in the summer for both asphalt types. The results indicate that for reasonable prediction of dynamic response in flexible pavements, the effect of the depth-dependent loading time on pavement temperature should be considered.

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

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

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

International Nuclear Information System (INIS)

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

1983-01-01

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

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

International Nuclear Information System (INIS)

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

1980-07-01

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

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.

Well-Balanced Second-Order Approximation of the Shallow Water Equations With Friction via Continuous Galerkin Finite Elements

Science.gov (United States)

Quezada de Luna, M.; Farthing, M.; Guermond, J. L.; Kees, C. E.; Popov, B.

2017-12-01

The Shallow Water Equations (SWEs) are popular for modeling non-dispersive incompressible water waves where the horizontal wavelength is much larger than the vertical scales. They can be derived from the incompressible Navier-Stokes equations assuming a constant vertical velocity. The SWEs are important in Geophysical Fluid Dynamics for modeling surface gravity waves in shallow regimes; e.g., in the deep ocean. Some common geophysical applications are the evolution of tsunamis, river flooding and dam breaks, storm surge simulations, atmospheric flows and others. This work is concerned with the approximation of the time-dependent Shallow Water Equations with friction using explicit time stepping and continuous finite elements. The objective is to construct a method that is at least second-order accurate in space and third or higher-order accurate in time, positivity preserving, well-balanced with respect to rest states, well-balanced with respect to steady sliding solutions on inclined planes and robust with respect to dry states. Methods fulfilling the desired goals are common within the finite volume literature. However, to the best of our knowledge, schemes with the above properties are not well developed in the context of continuous finite elements. We start this work based on a finite element method that is second-order accurate in space, positivity preserving and well-balanced with respect to rest states. We extend it by: modifying the artificial viscosity (via the entropy viscosity method) to deal with issues of loss of accuracy around local extrema, considering a singular Manning friction term handled via an explicit discretization under the usual CFL condition, considering a water height regularization that depends on the mesh size and is consistent with the polynomial approximation, reducing dispersive errors introduced by lumping the mass matrix and others. After presenting the details of the method we show numerical tests that demonstrate the well

A layer-wise MITC9 finite element for the free-vibration analysis of plates with piezo-patches

Directory of Open Access Journals (Sweden)

Maria Cinefra

2015-04-01

Full Text Available The present article considers the free-vibration analysis of plate structures with piezoelectric patches by means of a plate finite element with variable through-the-thickness layer-wise kinematic. The refined models used are derived from Carrera’s Unified Formulation (CUF and they permit the vibration modes along the thickness to be accurately described. The finite-element method is employed and the plate element implemented has nine nodes, and the mixed interpolation of tensorial component (MITC method is used to contrast the membrane and shear locking phenomenon. The related governing equations are derived from the principle of virtual displacement, extended to the analysis of electromechanical problems. An isotropic plate with piezoelectric patches is analyzed, with clamped-free boundary conditions and subjected to open- and short-circuit configurations. The results, obtained with different theories, are compared with the higher-order type solutions given in the literature. The conclusion is reached that the plate element based on the CUF is more suitable and efficient compared to the classical models in the study of multilayered structures embedding piezo-patches.

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

Efficient Computation of Info-Gap Robustness for Finite Element Models

International Nuclear Information System (INIS)

Stull, Christopher J.; Hemez, Francois M.; Williams, Brian J.

2012-01-01

A recent research effort at LANL proposed info-gap decision theory as a framework by which to measure the predictive maturity of numerical models. Info-gap theory explores the trade-offs between accuracy, that is, the extent to which predictions reproduce the physical measurements, and robustness, that is, the extent to which predictions are insensitive to modeling assumptions. Both accuracy and robustness are necessary to demonstrate predictive maturity. However, conducting an info-gap analysis can present a formidable challenge, from the standpoint of the required computational resources. This is because a robustness function requires the resolution of multiple optimization problems. This report offers an alternative, adjoint methodology to assess the info-gap robustness of Ax = b-like numerical models solved for a solution x. Two situations that can arise in structural analysis and design are briefly described and contextualized within the info-gap decision theory framework. The treatments of the info-gap problems, using the adjoint methodology are outlined in detail, and the latter problem is solved for four separate finite element models. As compared to statistical sampling, the proposed methodology offers highly accurate approximations of info-gap robustness functions for the finite element models considered in the report, at a small fraction of the computational cost. It is noted that this report considers only linear systems; a natural follow-on study would extend the methodologies described herein to include nonlinear systems.

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

Towards a unified solution of localization failure with mixed finite elements

Science.gov (United States)

Benedetti, Lorenzo; Cervera, Miguel; Chiumenti, Michele; Zeidler, Antonia; Fischer, Jan-Thomas

2015-04-01

interpolations. As a fundamental enhancement of the displacement-pressure formulation above mentioned, this kind of formulation benefits of the following advantages: it provides enhanced rate of convergence for the strain (and stress) and it is able to deal with incompressible situations. The method is completed with constitutive laws from Von Mises and Drucker-Prager local plasticity models with nonlinear strain softening. Moreover, global and local error norms are discussed to support the advantages of the proposed method. Then, numerical examples of stability analysis of slopes are presented to demonstrate the capability of the method. It will be shown that not only soil slopes can be modeled but also snow avalanche release and their weak layer fracture can be similarly treated. Consequently, this formulation appears to be a general and accurate tool for the solution of mechanical problem involving failure with localization bands [3,4]. References [1] Y.R. Rashid, 'Ultimate strength analysis of prestressed concrete pressure vessels', Nuclear Engineering and Design, Volume 7, Issue 4, April, Pages 334-344, 1968. [2] M. Cervera, M. Chiumenti, D. Di Capua. 'Benchmarking on bifurcation and localization in J 2 plasticity for plane stress and plane strain conditions.' Computer Methods in Applied Mechanics and Engineering, Vol. 241-244, Pages 206-224, 2012. [3] L. Benedetti, M. Cervera, M. Chiumenti. 'Stress-accurate mixed FEM for soil failure under shallow foundations involving strain localization in plasticity' Computers and Geotechnics, Vol. 64, pp. 32-47, 2015. [4] Cervera, M., Chiumenti, M., Benedetti, L., Codina, R. 'Mixed stabilized finite element methods in nonlinear solid mechanics. Part III: Compressible and incompressible plasticity' Computer Methods in Applied Mechanics and Engineering, to appear, 2015.

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

Moving finite element method for ICF target implosion

Science.gov (United States)

Furuta, J.; Kawata, S.; Niu, K.

1985-03-01

One dimensional hydrodynamic codes for the analysis of internal confinement fusion (ICF) target implosion which include various effects were developed, but most of them utilize the artificial viscosity (e.g., Von Neumann's viscosity) which cannot reveal accurately the shock waves. A gain of ICF target implosion is much due to the dissipation at the shock fronts, so it is necessary to express correctly the shock waves which are affected by the viscosity. The width of the shock waves is usually a few times as large as the length of mean free path, therefore the meshes for the shock waves must be set to about 10 to the 4th to 10 to the 5th power. It is a serious problem because of the computational memories or CPU time. In the moving finite element (MPE) method, both nodal amplitudes and nodal positions move continuously with time in such a way as to satisfy simultaneous ordinary differential equations (OPDs) which minimize partial differential equation (PDE) residuals.

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

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.

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.

Non-linear finite element analyses of wide plate fracture mechanics experiments

International Nuclear Information System (INIS)

Harrop, L.P.; Gibson, S.

1988-06-01

A series of centre-cracked, wide plate fracture mechanics tests is being conducted with plates made from 0.36% carbon steel. This report gives an account of post-test finite element analyses performed to compare with the results of one of these tests (designated CSTP4) and a pre-test analysis of the next test which has a slightly different geometry (CSTP5). The plates are relatively thick (75mm) and have a width of 1.62m. The finite element analyses use a two-dimensional plane stress mesh. The work shows good agreement between the post-test analysis results and the overall experimental results for CSTP4. It is not expected that the analysis results will be accurate within the dimensions of the process zone ahead of the crack tip; the mesh is not sufficient for this. A vital ingredient in attaining the good overall agreement is the representation of the actual stress-strain curve of the material. The predicted response of test CSTP5 is markedly different from that of CSTP4 even though the only change is the increase in the height of the plate. In particular the shape and size of the plastic zone ahead of the crack tip is quite different in the two tests at the same nominal remote applied load. (author)

Finite Element Study on Continuous Rotating versus Reciprocating Nickel-Titanium Instruments.

Science.gov (United States)

El-Anwar, Mohamed I; Yousief, Salah A; Kataia, Engy M; El-Wahab, Tarek M Abd

2016-01-01

In the present study, GTX and ProTaper as continuous rotating endodontic files were numerically compared with WaveOne reciprocating file using finite element analysis, aiming at having a low cost, accurate/trustworthy comparison as well as finding out the effect of instrument design and manufacturing material on its lifespan. Two 3D finite element models were especially prepared for this comparison. Commercial engineering CAD/CAM package was used to model full detailed flute geometries of the instruments. Multi-linear materials were defined in analysis by using real strain-stress data of NiTi and M-Wire. Non-linear static analysis was performed to simulate the instrument inside root canal at a 45° angle in the apical portion and subjected to 0.3 N.cm torsion. The three simulations in this study showed that M-Wire is slightly more resistant to failure than conventional NiTi. On the other hand, both materials are fairly similar in case of severe locking conditions. For the same instrument geometry, M-Wire instruments may have longer lifespan than the conventional NiTi ones. In case of severe locking conditions both materials will fail similarly. Larger cross sectional area (function of instrument taper) resisted better to failure than the smaller ones, while the cross sectional shape and its cutting angles could affect instrument cutting efficiency.

A New Energy-Based Method for 3-D Finite-Element Nonlinear Flux Linkage computation of Electrical Machines

DEFF Research Database (Denmark)

Lu, Kaiyuan; Rasmussen, Peter Omand; Ritchie, Ewen

2011-01-01

This paper presents a new method for computation of the nonlinear flux linkage in 3-D finite-element models (FEMs) of electrical machines. Accurate computation of the nonlinear flux linkage in 3-D FEM is not an easy task. Compared to the existing energy-perturbation method, the new technique......-perturbation method. The new method proposed is validated using experimental results on two different permanent magnet machines....

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

Fast solution of neutron diffusion problem by reduced basis finite element method

International Nuclear Information System (INIS)

Chunyu, Zhang; Gong, Chen

2018-01-01

Highlights: •An extremely efficient method is proposed to solve the neutron diffusion equation with varying the cross sections. •Three orders of speedup is achieved for IAEA benchmark problems. •The method may open a new possibility of efficient high-fidelity modeling of large scale problems in nuclear engineering. -- Abstract: For the important applications which need carry out many times of neutron diffusion calculations such as the fuel depletion analysis and the neutronics-thermohydraulics coupling analysis, fast and accurate solutions of the neutron diffusion equation are demanding but necessary. In the present work, the certified reduced basis finite element method is proposed and implemented to solve the generalized eigenvalue problems of neutron diffusion with variable cross sections. The order reduced model is built upon high-fidelity finite element approximations during the offline stage. During the online stage, both the k eff and the spatical distribution of neutron flux can be obtained very efficiently for any given set of cross sections. Numerical tests show that a speedup of around 1100 is achieved for the IAEA two-dimensional PWR benchmark problem and a speedup of around 3400 is achieved for the three-dimensional counterpart with the fission cross-sections, the absorption cross-sections and the scattering cross-sections treated as parameters.

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

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.

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

Solution of three-dimensional energy equation using finite element method

International Nuclear Information System (INIS)

Bhasin, V.; Singh, R.K.; Dutta, B.K.; Kushwaha, H.S.

1993-01-01

In the present work an attempt has been made to formulate an efficient 3-D finite element program for solving coupled momentum-energy equation with unsymmetric frontal solver and a suitable upwinding scheme. Based on the above solution technique of energy equation it can be concluded that upwinding scheme can lead to fairly accurate and smooth results even with coarse mesh. Otherwise the mesh size requirement will be extremely stringent for most of the practical problems. With upwinding the additional computer time required is marginally more. This effort has resulted in getting practical solution for large size real life problems in nuclear industry. The program was used for computation of temperature field in heavy water moderator of Madras Atomic Power Station (MAPS) reactor, in new mode of operation. (author). 9 refs., 7 figs

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

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

Finite groups with three conjugacy class sizes of some elements

Indian Academy of Sciences (India)

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

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

International Nuclear Information System (INIS)

Ackroyd, R.T.

1987-01-01

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

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

Science.gov (United States)

Romero, Ignacio; Segurado, Javier; LLorca, Javier

2008-04-01

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

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

International Nuclear Information System (INIS)

Romero, Ignacio; Segurado, Javier; LLorca, Javier

2008-01-01

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

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

NARCIS (Netherlands)

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

2016-01-01

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

Three-dimensional Finite Element Modelling of Composite Slabs for High Speed Rails

Science.gov (United States)

Mlilo, Nhlanganiso; Kaewunruen, Sakdirat

2017-12-01

Currently precast steel-concrete composite slabs are being considered on railway bridges as a viable alternative replacement for timber sleepers. However, due to their nature and the loading conditions, their behaviour is often complex. Present knowledge of the behaviour of precast steel-concrete composite slabs subjected to rail loading is limited. FEA is an important tool used to simulate real life behaviour and is widely accepted in many disciples of engineering as an alternative to experimental test methods, which are often costly and time consuming. This paper seeks to detail FEM of precast steel-concrete slabs subjected to standard in-service loading in high-speed rail with focus on the importance of accurately defining material properties, element type, mesh size, contacts, interactions and boundary conditions that will give results representative of real life behaviour. Initial finite element model show very good results, confirming the accuracy of the modelling procedure

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.

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

International Nuclear Information System (INIS)

Lee Hae Sung.

1991-01-01

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

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.

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.

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.

2.5-D frequency-domain viscoelastic wave modelling using finite-element method

Science.gov (United States)

Zhao, Jian-guo; Huang, Xing-xing; Liu, Wei-fang; Zhao, Wei-jun; Song, Jian-yong; Xiong, Bin; Wang, Shang-xu

2017-10-01

2-D seismic modelling has notable dynamic information discrepancies with field data because of the implicit line-source assumption, whereas 3-D modelling suffers from a huge computational burden. The 2.5-D approach is able to overcome both of the aforementioned limitations. In general, the earth model is treated as an elastic material, but the real media is viscous. In this study, we develop an accurate and efficient frequency-domain finite-element method (FEM) for modelling 2.5-D viscoelastic wave propagation. To perform the 2.5-D approach, we assume that the 2-D viscoelastic media are based on the Kelvin-Voigt rheological model and a 3-D point source. The viscoelastic wave equation is temporally and spatially Fourier transformed into the frequency-wavenumber domain. Then, we systematically derive the weak form and its spatial discretization of 2.5-D viscoelastic wave equations in the frequency-wavenumber domain through the Galerkin weighted residual method for FEM. Fixing a frequency, the 2-D problem for each wavenumber is solved by FEM. Subsequently, a composite Simpson formula is adopted to estimate the inverse Fourier integration to obtain the 3-D wavefield. We implement the stiffness reduction method (SRM) to suppress artificial boundary reflections. The results show that this absorbing boundary condition is valid and efficient in the frequency-wavenumber domain. Finally, three numerical models, an unbounded homogeneous medium, a half-space layered medium and an undulating topography medium, are established. Numerical results validate the accuracy and stability of 2.5-D solutions and present the adaptability of finite-element method to complicated geographic conditions. The proposed 2.5-D modelling strategy has the potential to address modelling studies on wave propagation in real earth media in an accurate and efficient way.

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

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.

Ship vibration analysis by finite element technique. Pt. II: Vibration analysis / Analyse van scheepstrillingen door middel van de elementenmethode. Dl. II: Trillingsanalyse

NARCIS (Netherlands)

Hylarides, S.

1971-01-01

In the calculation of the natural frequencies of ships more accurate values are expected when the shell-like structure of ships is taken into account by the finite element technique, especially in the higher-node vibration modes. To avoid large matrix systems an elimination process has been

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 finite element method with overlapping meshes for free-boundary axisymmetric plasma equilibria in realistic geometries

Science.gov (United States)

Heumann, Holger; Rapetti, Francesca

2017-04-01

Existing finite element implementations for the computation of free-boundary axisymmetric plasma equilibria approximate the unknown poloidal flux function by standard lowest order continuous finite elements with discontinuous gradients. As a consequence, the location of critical points of the poloidal flux, that are of paramount importance in tokamak engineering, is constrained to nodes of the mesh leading to undesired jumps in transient problems. Moreover, recent numerical results for the self-consistent coupling of equilibrium with resistive diffusion and transport suggest the necessity of higher regularity when approximating the flux map. In this work we propose a mortar element method that employs two overlapping meshes. One mesh with Cartesian quadrilaterals covers the vacuum chamber domain accessible by the plasma and one mesh with triangles discretizes the region outside. The two meshes overlap in a narrow region. This approach gives the flexibility to achieve easily and at low cost higher order regularity for the approximation of the flux function in the domain covered by the plasma, while preserving accurate meshing of the geometric details outside this region. The continuity of the numerical solution in the region of overlap is weakly enforced by a mortar-like mapping.

A structurally detailed finite element human head model for simulation of transcranial magnetic stimulation.

Science.gov (United States)

Chen, Ming; Mogul, David Jeffery

2009-04-30

Computational studies of the head utilizing finite element models (FEMs) have been used to investigate a wide variety of brain-electromagnetic (EM) field interaction phenomena including magnetic stimulation of the head using transcranial magnetic stimulation (TMS), direct electric stimulation of the brain for electroconvulsive therapy, and electroencephalography source localization. However, no human head model of sufficient complexity for studying the biophysics under these circumstances has been developed which utilizes structures at both the regional and cellular levels and provides well-defined smooth boundaries between tissues of different conductivities and orientations. The main barrier for building such accurate head models is the complex modeling procedures that include 3D object reconstruction and optimized meshing. In this study, a structurally detailed finite element model of the human head was generated that includes details to the level of cerebral gyri and sulci by combining computed tomography and magnetic resonance images. Furthermore, cortical columns that contain conductive processes of pyramidal neurons traversing the neocortical layers were included in the head model thus providing structure at or near the cellular level. These refinements provide a much more realistic model to investigate the effects of TMS on brain electrophysiology in the neocortex.

Variable order spherical harmonic expansion scheme for the radiative transport equation using finite elements

International Nuclear Information System (INIS)

Surya Mohan, P.; Tarvainen, Tanja; Schweiger, Martin; Pulkkinen, Aki; Arridge, Simon R.

2011-01-01

Highlights: → We developed a variable order global basis scheme to solve light transport in 3D. → Based on finite elements, the method can be applied to a wide class of geometries. → It is computationally cheap when compared to the fixed order scheme. → Comparisons with local basis method and other models demonstrate its accuracy. → Addresses problems encountered n modeling of light transport in human brain. - Abstract: We propose the P N approximation based on a finite element framework for solving the radiative transport equation with optical tomography as the primary application area. The key idea is to employ a variable order spherical harmonic expansion for angular discretization based on the proximity to the source and the local scattering coefficient. The proposed scheme is shown to be computationally efficient compared to employing homogeneously high orders of expansion everywhere in the domain. In addition the numerical method is shown to accurately describe the void regions encountered in the forward modeling of real-life specimens such as infant brains. The accuracy of the method is demonstrated over three model problems where the P N approximation is compared against Monte Carlo simulations and other state-of-the-art methods.

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

Modelling of Conveyor Belt Passage by Driving Drum Using Finite Element Methods

Directory of Open Access Journals (Sweden)

Nikoleta Mikušová

2017-12-01

Full Text Available The finite element methods are used in many disciplines by the development of products, typically in mechanical engineering (for example in automotive industry, biomechanics, etc.. Some modern programs of the finite element's methods have specific tools (electromagnetic, fluid and structural simulations. The finite elements methods allow detailed presentation of structures by bending or torsion, complete design, testing and optimization before the prototype production. The aims of this paper were to the model of conveyor belt passage by driving drum. The model was created by the program Abaqus CAE. The created model presented data about forces, pressures, and deformation of the belt conveyor.

Finite-element pre-analysis for pressurized thermoshock tests

International Nuclear Information System (INIS)

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

1992-05-01

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

Finite element analysis of FRP-strengthened RC beams

Directory of Open Access Journals (Sweden)

Teeraphot Supaviriyakit

2004-05-01

Full Text Available This paper presents a non-linear finite element analysis of reinforced concrete beam strengthened with externally bonded FRP plates. The finite element modeling of FRP-strengthened beams is demonstrated. Concrete and reinforcing bars are modeled together as 8-node isoparametric 2D RC element. The FRP plate is modeled as 8-node isoparametric 2D elastic element. The glue is modeled as perfect compatibility by directly connecting the nodes of FRP with those of concrete since there is no failure at the glue layer. The key to the analysis is the correct material models of concrete, steel and FRP. Cracks and steel bars are modeled as smeared over the entire element. Stress-strain properties of cracked concrete consist of tensile stress model normal to crack, compressive stress model parallel to crack and shear stress model tangential to crack. Stressstrain property of reinforcement is assumed to be elastic-hardening to account for the bond between concrete and steel bars. FRP is modeled as elastic-brittle material. From the analysis, it is found that FEM can predict the load-displacement relation, ultimate load and failure mode of the beam correctly. It can also capture the cracking process for both shear-flexural peeling and end peeling modes similar to the experiment.

Application of the finite element groundwater model FEWA to the engineered test facility

International Nuclear Information System (INIS)

Craig, P.M.; Davis, E.C.

1985-09-01

A finite element model for water transport through porous media (FEWA) has been applied to the unconfined aquifer at the Oak Ridge National Laboratory Solid Waste Storage Area 6 Engineered Test Facility (ETF). The model was developed in 1983 as part of the Shallow Land Burial Technology - Humid Task (ONL-WL14) and was previously verified using several general hydrologic problems for which an analytic solution exists. Model application and calibration, as described in this report, consisted of modeling the ETF water table for three specialized cases: a one-dimensional steady-state simulation, a one-dimensional transient simulation, and a two-dimensional transient simulation. In the one-dimensional steady-state simulation, the FEWA output accurately predicted the water table during a long period in which there were no man-induced or natural perturbations to the system. The input parameters of most importance for this case were hydraulic conductivity and aquifer bottom elevation. In the two transient cases, the FEWA output has matched observed water table responses to a single rainfall event occurring in February 1983, yielding a calibrated finite element model that is useful for further study of additional precipitation events as well as contaminant transport at the experimental site

Finite Element Method Based Modeling of Resistance Spot-Welded Mild Steel

Directory of Open Access Journals (Sweden)

Miloud Zaoui

Full Text Available Abstract This paper deals with Finite Element refined and simplified models of a mild steel spot-welded specimen, developed and validated based on quasi-static cross-tensile experimental tests. The first model was constructed with a fine discretization of the metal sheet and the spot weld was defined as a special geometric zone of the specimen. This model provided, in combination with experimental tests, the input data for the development of the second model, which was constructed with respect to the mesh size used in the complete car finite element model. This simplified model was developed with coarse shell elements and a spring-type beam element was used to model the spot weld behavior. The global accuracy of the two models was checked by comparing simulated and experimental load-displacement curves and by studying the specimen deformed shapes and the plastic deformation growth in the metal sheets. The obtained results show that both fine and coarse finite element models permit a good prediction of the experimental tests.

Study on the influence of finite element formulation and equation of motion solution scheme on FEM analysis results based on the asymmetrically loaded plate problem

Directory of Open Access Journals (Sweden)

Marcin Krzeszowiec

2015-03-01

Full Text Available Computer simulations of physical phenomena are at the moment common both in science and industry. The possibility of finding approximate solutions for complicated systems of differential equations, mathematically describing issues in the fields of mechanics, physics or chemistry, allows for shorten design and research time, often significantly reducing the need for expensive experimental studies or costly production of prototypes. However, the mentioned prevalence of these methods, particularly the Finite Element Method, resulted in analysis outcomes to be often in advance regarded as accurate ones. The purpose of the article is to showcase, on a simple stress analysis problem, how parameters such as the density of the finite element mesh, finite element formulation or integration scheme significantly influence on the simulation results and how easy it is to end up with the results that do not hold any physical sense, despite the fact that all the basic assumptions of correct analysis (suitable boundary conditions, total system energy stored etc. have been met. The results of this study can serve as a warning against premature conclusion drawing from calculations carried out by means of FEM simulation.[b]Keywords[/b]: computational mechanics, finite element method, shell elements, numerical integration

Image segmentation with a finite element method

DEFF Research Database (Denmark)

Bourdin, Blaise

1999-01-01

regularization results, make possible to imagine a finite element resolution method.In a first time, the Mumford-Shah functional is introduced and some existing results are quoted. Then, a discrete formulation for the Mumford-Shah problem is proposed and its $\\Gamma$-convergence is proved. Finally, some...

Investigations on Actuator Dynamics through Theoretical and Finite Element Approach

Directory of Open Access Journals (Sweden)

Somashekhar S. Hiremath

2010-01-01

Full Text Available This paper gives a new approach for modeling the fluid-structure interaction of servovalve component-actuator. The analyzed valve is a precision flow control valve-jet pipe electrohydraulic servovalve. The positioning of an actuator depends upon the flow rate from control ports, in turn depends on the spool position. Theoretical investigation is made for No-load condition and Load condition for an actuator. These are used in finite element modeling of an actuator. The fluid-structure-interaction (FSI is established between the piston and the fluid cavities at the piston end. The fluid cavities were modeled with special purpose hydrostatic fluid elements while the piston is modeled with brick elements. The finite element method is used to simulate the variation of cavity pressure, cavity volume, mass flow rate, and the actuator velocity. The finite element analysis is extended to study the system's linearized response to harmonic excitation using direct solution steady-state dynamics. It was observed from the analysis that the natural frequency of the actuator depends upon the position of the piston in the cylinder. This is a close match with theoretical and simulation results. The effect of bulk modulus is also presented in the paper.

A multilevel correction adaptive finite element method for Kohn-Sham equation

Science.gov (United States)

Hu, Guanghui; Xie, Hehu; Xu, Fei

2018-02-01

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

Finite element modeling of trolling-mode AFM.

Science.gov (United States)

Sajjadi, Mohammadreza; Pishkenari, Hossein Nejat; Vossoughi, Gholamreza

2018-06-01

Trolling mode atomic force microscopy (TR-AFM) has overcome many imaging problems in liquid environments by considerably reducing the liquid-resonator interaction forces. The finite element model of the TR-AFM resonator considering the effects of fluid and nanoneedle flexibility is presented in this research, for the first time. The model is verified by ABAQUS software. The effect of installation angle of the microbeam relative to the horizon and the effect of fluid on the system behavior are investigated. Using the finite element model, frequency response curve of the system is obtained and validated around the frequency of the operating mode by the available experimental results, in air and liquid. The changes in the natural frequencies in the presence of liquid are studied. The effects of tip-sample interaction on the excitation of higher order modes of the system are also investigated in air and liquid environments. Copyright © 2018 Elsevier B.V. All rights reserved.

Finite element simulations with ANSYS workbench 16

CERN Document Server

Lee , Huei-Huang

2015-01-01

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

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

KAUST Repository

Bangerth, Wolfgang

2011-12-01

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

Rigid body formulation in a finite element context with contact interaction

Science.gov (United States)

Refachinho de Campos, Paulo R.; Gay Neto, Alfredo

2018-03-01

The present work proposes a formulation to employ rigid bodies together with flexible bodies in the context of a nonlinear finite element solver, with contact interactions. Inertial contributions due to distribution of mass of a rigid body are fully developed, considering a general pole position associated with a single node, representing a rigid body element. Additionally, a mechanical constraint is proposed to connect a rigid region composed by several nodes, which is useful for linking rigid/flexible bodies in a finite element environment. Rodrigues rotation parameters are used to describe finite rotations, by an updated Lagrangian description. In addition, the contact formulation entitled master-surface to master-surface is employed in conjunction with the rigid body element and flexible bodies, aiming to consider their interaction in a rigid-flexible multibody environment. New surface parameterizations are presented to establish contact pairs, permitting pointwise interaction in a frictional scenario. Numerical examples are provided to show robustness and applicability of the methods.

Discrete maximum principle for the P1 - P0 weak Galerkin finite element approximations