Bearing capacity and rigidity of short plastic-concrete-tubal vertical columns under transverse load

Science.gov (United States)

Dolzhenko, A. V.; Naumov, A. E.; Shevchenko, A. E.

2018-03-01

The results of mathematical modeling in determining strain-stress distribution parameters of a short plastic-concrete-tubal vertical column under horizontal load as those in vertical constructions are described. Quantitative parameters of strain-stress distribution during vertical and horizontal loads and horizontal stiffness were determined by finite element modeling. The internal stress in the concrete column core was analyzed according to equivalent stress in Mohr theory of failure. It was determined that the bearing capacity of a short plastic- concrete-tubal vertical column is 25% higher in resistibility and 15% higher in rigidness than those of the caseless concrete columns equal in size. Cracks formation in the core of a short plastic-concrete-tubal vertical column happens under significantly bigger horizontal loads with less amount of concrete spent than that in caseless concrete columns. The significant increase of bearing capacity and cracking resistance of a short plastic-concrete-tubal vertical column under vertical and horizontal loads allows recommending them as highly effective and highly reliable structural wall elements in civil engineering.

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

International Nuclear Information System (INIS)

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

2014-01-01

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

Modeling and Analysis of Size-Dependent Structural Problems by Using Low- Order Finite Elements with Strain Gradient Plasticity

International Nuclear Information System (INIS)

Park, Moon Shik; Suh, Yeong Sung; Song, Seung

2011-01-01

An elasto-plastic finite element method using the theory of strain gradient plasticity is proposed to evaluate the size dependency of structural plasticity that occurs when the configuration size decreases to micron scale. For this method, we suggest a low-order plane and three-dimensional displacement-based elements, eliminating the need for a high order, many degrees of freedom, a mixed element, or super elements, which have been considered necessary in previous researches. The proposed method can be performed in the framework of nonlinear incremental analysis in which plastic strains are calculated and averaged at nodes. These strains are then interpolated and differentiated for gradient calculation. We adopted a strain-gradient-hardening constitutive equation from the Taylor dislocation model, which requires the plastic strain gradient. The developed finite elements are tested numerically on the basis of typical size-effect problems such as micro-bending, micro-torsion, and micro-voids. With respect to the strain gradient plasticity, i.e., the size effects, the results obtained by using the proposed method, which are simple in their calculation, are in good agreement with the experimental results cited in previously published papers

Finite Element Analysis and Die Design of Non-specific Engineering Structure of Aluminum Alloy during Extrusion

International Nuclear Information System (INIS)

Chen, D.-C.; Lu, Y.-Y.

2010-01-01

Aluminum extension applies to industrial structure, light load, framework rolls and conveyer system platform. Many factors must be controlled in processing the non-specific engineering structure (hollow shape) of the aluminum alloy during extrusion, to obtain the required plastic strain and desired tolerance values. The major factors include the forming angle of the die and temperature of billet and various materials. This paper employs rigid-plastic finite element (FE) DEFORM 3D software to investigate the plastic deformation behavior of an aluminum alloy (A6061, A5052, A3003) workpiece during extrusion for the engineering structure of the aluminum alloy. This work analyzes effective strain, effective stress, damage and die radius load distribution of the billet under various conditions. The analytical results confirm the suitability of the current finite element software for the non-specific engineering structure of aluminum alloy extrusion.

Low Reynolds number steady state flow through a branching network of rigid vessels: II. A finite element mixture model

NARCIS (Netherlands)

Huyghe, J.M.R.J.; Oomens, C.W.J.; Campen, van D.H.; Heethaar, R.M.

1989-01-01

This research aims at formulating and verifying a finite element mixture formulation for blood perfusion. The equations derived in a companion paper [3] are discretized according to the Galerkin method. A flow experiment in a rigid model of a vascular tree of about 500 vessels is performed in order

Probabilistic finite elements

Science.gov (United States)

Belytschko, Ted; Wing, Kam Liu

1987-01-01

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

Finite element implementation of strain-hardening Drucker–Prager plasticity model with application to tunnel excavation

Directory of Open Access Journals (Sweden)

K. Liu

2017-09-01

Full Text Available This paper presents a finite element implementation of a strain-hardening Drucker–Prager model and its application to tunnel excavation. The computational model was constructed based on the return mapping scheme, in which an elastic trial step was first executed, followed by plastic correction involving the Newton–Raphson method to return the predicted state of stresses to the supposed yield surface. By combining the plastic shear hardening rule and stress correction equations, the loading index for the strain-hardening Drucker–Prager model was solved. It is therefore possible to update the stresses, elastic and plastic strains, and slope of the yield locus at the end of each incremental step. As an illustrative example, an integration algorithm was incorporated into ABAQUS through the user subroutine UMAT to solve the tunnel excavation problem in strain-hardening Drucker–Prager rock formations. The obtained numerical results were found to be in excellent agreement with the available analytical solutions, thus indicating the validity and accuracy of the proposed UMAT code, as well as the finite element model.

Discussion of plastic instabilities of a pressure vessel by means of the finite element method

International Nuclear Information System (INIS)

Laemmer, H.; Ritter, B.; Tsakmakis, C.

1994-06-01

In this report the influence of geometrical and thermal imperfections on the load-carrying capacity of pressure vessels has been studied by means of the ABAQUS finite element code. Using a thermo-plasticity model for finite strains failure due to plastic instability has been investigated. The parameters needed for this model were identified with uniaxial tensile tests. The parameters were fitted in such a way that necking of the tensile rod was described in a satisfactory manner. Starting from the tensile rod the influence of imperfections on necking of a structure is discussed. Axisymmetric and three dimensional calculations showed that for a sufficient size of the imperfection necking can be obtained. As expected, the maximum internal pressure at failure and the beginning of localization of plastic deformations depend on the size of the imperfection. (orig.) [de

Preliminary analysis of construction of the test drift in boom clay at Mol using plasticity solutions and finite elements

International Nuclear Information System (INIS)

Mair, R.J.; Taylor, R.N.; Higgins, K.G.; Potts, D.M.

1989-01-01

Analyses have been undertaken on an advancing tunnel heading at great depth in a clay formation corresponding to the test drift construction at Mol. Belgium. Simplifying assumptions enable plasticity solutions to be used to model the behaviour of a tunnel heading in a linear elastic-perfectly plastic soil. Finite element analysis with the same soil model has been undertaken of the test drift construction, assuming axisymmetric conditions. The results are compared with the plasticity solutions and with the measurements of lining stresses, soil movements and pore pressures by SCK/CEN. Good agreement is obtained between the plasticity solutions and finite element analysis. The measured immediate build-up of stress on the linings is well-predicted and reasonable agreement is obtained between predicted and measured soil movements. The measured pore-pressure changes are poorly predicted by the analyses

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

Finite element analysis of elasto-plastic tee joints

International Nuclear Information System (INIS)

Powell, G.H.

1974-09-01

The theory and computational procedures used in the computer program B169TJ/EP for the analysis of elasto-plastic tee joints are described, and detailed user's guide is presented. The program is particularly applicable to joints conforming to the ANSI B16.9 Manufacturing Standard, but can also be applied to other joint geometries. The joint may be loaded by internal pressure and by arbitrary combinations of applied forces and moments at the ends of the branch and run pipes, and the loading sequence may be arbitrary. The joint material is assumed to yield according to the von Mises criterion, and to exhibit either linear kinematic hardening or nonlinear isotropic hardening after yield. The program makes use of the finite element and mesh generation procedures previously applied in the elastic stress analysis program B16.9TJ/ SA, with minor modifications. (U.S.)

Robust and efficient handling of yield surface discontinuities in elasto-plastic finite element calculations

DEFF Research Database (Denmark)

Clausen, Johan Christian; Damkilde, Lars; Andersen, Lars Vabbersgaard

2015-01-01

Purpose – The purpose of this paper is to present several methods on how to deal with yield surface discontinuities. The explicit formulations, first presented by Koiter (1953), result in multisingular constitutive matrices which can cause numerical problems in elasto-plastic finite element...... documented in the literature all present “easy” calculation examples, e.g. low friction angles and few elements. The amendments presented in this paper result in robust elasto-plastic computations, making the solution of “hard” problems possible without introducing approximations in the yield surfaces...... calculations. These problems, however, are not documented in previous literature. In this paper an amendment to the Koiter formulation of the constitutive matrices for stress points located on discontinuities is proposed. Design/methodology/approach – First, a review of existing methods of handling yield...

A work-hardening rule for finite elastic-plastic deformation of metals at elevated temperatures

International Nuclear Information System (INIS)

Lee, L.H.N.; Horng, J.T.

1975-01-01

The paper is concerned with an extension of Prager-Ziegler's kinematic work-hardening rule for infinitesimal elastic-plastic deformation to a work-hardening rule for finite elastic-plastic deformation of a polycrystalline metal. It is shown that the finite work-hardening rule, which accounts for the Bauschinger and temperature effects within certain pressure and temperature ranges, satisfies certain invariant, continuity and thermodynamic requirements. A description of the kinematics of an elastic-plastic body is employed with reference to three separate configurations: initial, current and an intermediate configuration. The intermediate configuration is a conceptual, local configuration obtained by removing the stress and temperature changes in the neighborhood of an element. A rigid body rotation of the intermediate configuration is allowed. Piola-Kirchhoff stresses and Green deformation tensors referred to the initial and intermediate configurations are employed as stress and strain measures. The plastic deformation has been associated with the motion and production of dislocations. It has been observed that the motion of mobile dislocations usually occur in the narrow slip bands in each grain, leaving the basic lattice structure practically intact, so that the macroscopic elastic properties of the material are essentially independent of plastic deformation. Employing this fact and the thermodynamic laws, a simplified elastic stress-strain relationship of the plastically deformed material, which agrees with the results of Naghdi and Trapp, is obtained

Finite element analysis of large elasto-plastic deformation for sealing ring in nuclear pressure vessel

International Nuclear Information System (INIS)

Xiao Xuejian; Chen Ruxin

1995-02-01

Based on the R. Hills incremental virtual power principle and the elasto-plastic constitution equation for large deformation and by considering physical nonlinear, geometric nonlinear and thermal effects, a plane and axisymmetric finite element equation for thermal large elasto-plastic deformation has been established in the Euler description. The corresponding analysis program ATLEPD has been also complied for thermal large elasto-plastic deformation process of O-ring in RPV. The variations of stress, strain, contact specific pressure, mesh deformation and the aspects of spring back in upsetting and spring back process have been also investigated. Numerical results are fairly consistent with experimental ones. (5 figs., 4 tabs.)

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

Non-uniform plastic deformation of micron scale objects

DEFF Research Database (Denmark)

Niordson, Christian Frithiof; Hutchinson, J. W.

2003-01-01

Significant increases in apparent flow strength are observed when non-uniform plastic deformation of metals occurs at the scale ranging from roughly one to ten microns. Several basic plane strain problems are analyzed numerically in this paper based on a new formulation of strain gradient...... plasticity. The problems are the tangential and normal loading of a finite rectangular block of material bonded to rigid platens and having traction-free ends, and the normal loading of a half-space by a flat, rigid punch. The solutions illustrate fundamental features of plasticity at the micron scale...... that are not captured by conventional plasticity theory. These include the role of material length parameters in establishing the size dependence of strength and the elevation of resistance to plastic flow resulting from constraint on plastic flow at boundaries. Details of the finite element method employed...

Finite element modeling of penetration of rigid cylindrical bar impacting on a clamped circular plate, employing a strain rate dependent flow rate and a Gruneisen equation of state

International Nuclear Information System (INIS)

Tariq, M.; Khan, I.A.

2003-01-01

A time dependent Finite Element simulation of penetration of a rigid cylindrical bar impacting on a copper plate is conducted, to demonstrate how material behavior appears to change when Johnson-Cook plasticity rule is employed along with a Gruneisen, equation of state with cubic shock velocity-particle relationship, and defining pressure both for compressed and expanded materials, as compared to the behavior when only isotropic strain-hardening model is employed. The bar impacts the plate with a velocity of 1000 m/s, and penetrates the plate, a portion of it coming out of the other side. Results are obtained and compared taking both an isotropic strain-hardening model, and a model incorporating Johnson-Cook flow rule along with Gruneisen equation of state. (author)

Rigid-plastic seismic design of reinforced concrete structures

DEFF Research Database (Denmark)

Costa, Joao Domingues; Bento, R.; Levtchitch, V.

2007-01-01

structural strength with respect to a pre-defined performance parameter using a rigid-plastic response spectrum, which is characteristic of the ground motion alone. The maximum strength demand at any point is solely dependent on the intensity of the ground motion, which facilitates the task of distributing......In this paper a new seismic design procedure for Reinforced Concrete (R/C) structures is proposed-the Rigid-Plastic Seismic Design (RPSD) method. This is a design procedure based on Non-Linear Time-History Analysis (NLTHA) for systems expected to perform in the non-linear range during a lifetime...... earthquake event. The theoretical background is the Theory of Plasticity (Rigid-Plastic Structures). Firstly, a collapse mechanism is chosen and the corresponding stress field is made safe outside the regions where plastic behaviour takes place. It is shown that this allows the determination of the required...

Finite element analysis of cylindrical indentation for determining plastic properties of materials in small volumes

International Nuclear Information System (INIS)

Lu, Y Charles; Kurapati, Siva N V R K; Yang Fuqian

2008-01-01

The cylindrical indentation is analysed, using the finite element method, for determining the plastic properties of elastic-plastic materials and the effect of strain hardening. The results are compared with those obtained from spherical indentation, the commonly used technique for measuring plastic properties of materials in small volumes. The analysis shows that the deformation under a cylindrical indenter quickly reaches a fully plastic state and that the size (diameter) of the plastic zone remains constant during further indentation. The indentation load is proportional to the indentation depth at large indentation depth, from which the indentation pressure P m at the onset of yielding can be readily extrapolated. The analysis of cylindrical indentation suggests that it does not need parameters such as impression radius (a) and contact stiffness (S) for determining the plastic behaviour of materials. Thus, the cylindrical indentation can suppress the uncertainties in measuring material properties

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

Efficient Non-Linear Finite Element Implementation of Elasto-Plasticity for Geotechnical Problems

DEFF Research Database (Denmark)

Clausen, Johan

-Coulomb yield criterion and the corresponding plastic potential possess corners and an apex, which causes numerical difficulties. A simple, elegant and efficient solution to these problems is presented in this thesis. The solution is based on a transformation into principal stress space and is valid for all...... linear isotropic plasticity models in which corners and apexes are encountered. The validity and merits of the proposed solution are examined in relation to the Mohr-Coulomb and the Modified Mohr-Coulomb material models. It is found that the proposed method compares well with existing methods......-Brown material. The efficiency and validity are demonstrated by comparing the finite-element results with well-known solutions for simple geometries. A common geotechnical problem is the assessment of slope stability. For slopes with simple geometries and consisting of a linear Mohr-Coulomb material, this can...

Mechanical evolution of transpression zones affected by fault interactions: Insights from 3D elasto-plastic finite element models

Science.gov (United States)

Nabavi, Seyed Tohid; Alavi, Seyed Ahmad; Mohammadi, Soheil; Ghassemi, Mohammad Reza

2018-01-01

The mechanical evolution of transpression zones affected by fault interactions is investigated by a 3D elasto-plastic mechanical model solved with the finite-element method. Ductile transpression between non-rigid walls implies an upward and lateral extrusion. The model results demonstrate that a, transpression zone evolves in a 3D strain field along non-coaxial strain paths. Distributed plastic strain, slip transfer, and maximum plastic strain occur within the transpression zone. Outside the transpression zone, fault slip is reduced because deformation is accommodated by distributed plastic shear. With progressive deformation, the σ3 axis (the minimum compressive stress) rotates within the transpression zone to form an oblique angle to the regional transport direction (∼9°-10°). The magnitude of displacement increases faster within the transpression zone than outside it. Rotation of the displacement vectors of oblique convergence with time suggests that transpression zone evolves toward an overall non-plane strain deformation. The slip decreases along fault segments and with increasing depth. This can be attributed to the accommodation of bulk shortening over adjacent fault segments. The model result shows an almost symmetrical domal uplift due to off-fault deformation, generating a doubly plunging fold and a 'positive flower' structure. Outside the overlap zone, expanding asymmetric basins subside to 'negative flower' structures on both sides of the transpression zone and are called 'transpressional basins'. Deflection at fault segments causes the fault dip fall to less than 90° (∼86-89°) near the surface (∼1.5 km). This results in a pure-shear-dominated, triclinic, and discontinuous heterogeneous flow of the transpression zone.

Finite element evaluation of elasto-plastic accommodation energies during solid state transformations: Coherent, spherical precipitate in finite matrix

International Nuclear Information System (INIS)

Sen, S.; Balasubramaniam, R.; Sethuraman, R.

1996-01-01

The molar volume difference between the matrix and the precipitate phases in the case of solid state phase transformations results in the creation of stain energy in the system due to the misfit strains. A finite element model based on the initial strain approach is proposed to evaluate elasto-plastic accommodation energies during solid state transformation. The three-dimensional axisymmetric model has been used to evaluate energies as a function of transformation for α-β hydrogen transformations in the Nb-H system. The transformation has been analyzed for the cases of transformation progressing both from the center to surface and from the surface to center of the system. The effect of plastic deformation has been introduced to make the model realistic, specifically to the Nb-NbH phase transformation which involves a 4% linear misfit strain. It has been observed that plastic deformation reduces the strain energies compared to the linear elastic analysis

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.

Elastic-plastic analysis of an axi-symmetric problem by a finite element method

International Nuclear Information System (INIS)

Isozaki, Toshikuni

1984-06-01

Generally speaking, many structures are designed and fabricated on the basis of an axi-symmetric structure. Finite Element Method is the capable method to solve these axi-symmetric problems beyond the elastic limit. As the first step to solve these problems, the computer program for the elastic-plastic analysis of the axi-symmetric problem is composed. The basic program is based upon that described in Zienkiewicz's text book to solve the elastic plane stress problem, taking the plastic stress matrix by Yamada's method into consideration and it is converted to solve the axi-symmetric problem. For the verification of the program, the plane strain problem of a cylindrical tube under internal pressure was solved. The computed results were compared with those shown in ADINA's user's manual. They showed close agreement. (author)

Analysis of Elastic-Plastic J Integrals for 3-Dimensional Cracks Using Finite Element Alternating Method

International Nuclear Information System (INIS)

Park, Jai Hak

2009-01-01

SGBEM(Symmetric Galerkin Boundary Element Method)-FEM alternating method has been proposed by Nikishkov, Park and Atluri. In the proposed method, arbitrarily shaped three-dimensional crack problems can be solved by alternating between the crack solution in an infinite body and the finite element solution without a crack. In the previous study, the SGBEM-FEM alternating method was extended further in order to solve elastic-plastic crack problems and to obtain elastic-plastic stress fields. For the elastic-plastic analysis the algorithm developed by Nikishkov et al. is used after modification. In the algorithm, the initial stress method is used to obtain elastic-plastic stress and strain fields. In this paper, elastic-plastic J integrals for three-dimensional cracks are obtained using the method. For that purpose, accurate values of displacement gradients and stresses are necessary on an integration path. In order to improve the accuracy of stress near crack surfaces, coordinate transformation and partitioning of integration domain are used. The coordinate transformation produces a transformation Jacobian, which cancels the singularity of the integrand. Using the developed program, simple three-dimensional crack problems are solved and elastic and elastic-plastic J integrals are obtained. The obtained J integrals are compared with the values obtained using a handbook solution. It is noted that J integrals obtained from the alternating method are close to the values from the handbook

Finite element modelling of shot peening process: Prediction of the compressive residual stresses, the plastic deformations and the surface integrity

International Nuclear Information System (INIS)

Frija, M.; Hassine, T.; Fathallah, R.; Bouraoui, C.; Dogui, A.

2006-01-01

This paper presents a numerical simulation of the shot peening process using finite element method. The majority of the controlling parameters of the process have been taken into account. The shot peening loading has been characterised by using energy equivalence between the dynamic impact and a static indentation of a peening shot in the treated surface. The behaviour of the subjected material is supposed to be elastic plastic with damage. An integrated law of the damage proposed by Lemaitre and Chaboche has been used. The proposed model leads to obtain the residual stress, the plastic deformation profiles and the surface damage. An application on a shot peened Ni-based super alloy Waspaloy has been carried out. The comparison of the residual stresses, obtained by X-ray diffraction method and by finite element calculation, shows a good correlation. The in-depth profile of the plastic deformations and the superficial damage values are in good agreement with the experimental observations

Simulation of finite-strain inelastic phenomena governed by creep and plasticity

Science.gov (United States)

Li, Zhen; Bloomfield, Max O.; Oberai, Assad A.

2017-11-01

Inelastic mechanical behavior plays an important role in many applications in science and engineering. Phenomenologically, this behavior is often modeled as plasticity or creep. Plasticity is used to represent the rate-independent component of inelastic deformation and creep is used to represent the rate-dependent component. In several applications, especially those at elevated temperatures and stresses, these processes occur simultaneously. In order to model these process, we develop a rate-objective, finite-deformation constitutive model for plasticity and creep. The plastic component of this model is based on rate-independent J_2 plasticity, and the creep component is based on a thermally activated Norton model. We describe the implementation of this model within a finite element formulation, and present a radial return mapping algorithm for it. This approach reduces the additional complexity of modeling plasticity and creep, over thermoelasticity, to just solving one nonlinear scalar equation at each quadrature point. We implement this algorithm within a multiphysics finite element code and evaluate the consistent tangent through automatic differentiation. We verify and validate the implementation, apply it to modeling the evolution of stresses in the flip chip manufacturing process, and test its parallel strong-scaling performance.

An explanation for the shape of nanoindentation unloading curves based on finite element simulation

International Nuclear Information System (INIS)

Bolshakov, A.; Pharr, G.M.

1995-01-01

Current methods for measuring hardness and modulus from nanoindentation load-displacement data are based on Sneddon's equations for the indentation of an elastic half-space by an axially symmetric rigid punch. Recent experiments have shown that nanoindentation unloading data are distinctly curved in a manner which is not consistent with either the flat punch or the conical indenter geometries frequently used in modeling, but are more closely approximated by a parabola of revolution. Finite element simulations for conical indentation of an elastic-plastic material are presented which corroborate the experimental observations, and from which a simple explanation for the shape of the unloading curve is derived. The explanation is based on the concept of an effective indenter shape whose geometry is determined by the shape of the plastic hardness impression formed during indentation

Plastic Limit Loads for Slanted Circumferential Through-Wall Cracked Pipes Using 3D Finite-Element Limit Analyses

International Nuclear Information System (INIS)

Jang, Hyun Min; Cho, Doo Ho; Kim, Young Jin; Huh, Nam Su; Shim, Do Jun; Choi, Young Hwan; Park, Jung Soon

2011-01-01

On the basis of detailed 3D finite-element (FE) limit analyses, the plastic limit load solutions for pipes with slanted circumferential through-wall cracks (TWCs) subjected to axial tension, global bending, and internal pressure are reported. The FE model and analysis procedure employed in the present numerical study were validated by comparing the present FE results with existing solutions for plastic limit loads of pipes with idealized TWCs. For the quantification of the effect of slanted crack on plastic limit load, slant correction factors for calculating the plastic limit loads of pipes with slanted TWCs from pipes with idealized TWCs are newly proposed from extensive 3D FE calculations. These slant-correction factors are presented in tabulated form for practical ranges of geometry and for each set of loading conditions

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

International Nuclear Information System (INIS)

Miller, J.D.

1985-01-01

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

Hierarchical finite element modeling of SiCp/Al2124 T4 composites with dislocation plasticity and size dependent failure

International Nuclear Information System (INIS)

Suh, Yeong Sung; Kim, Yong Bae

2012-01-01

The strength of particle reinforced metal matrix composites is, in general, known to be increased by the geometrically necessary dislocations punched around a particle that form during cooling after consolidation because of coefficient of thermal expansion (CTE) mismatch between the particle and the matrix. An additional strength increase may also be observed, since another type of geometrically necessary dislocation can be formed during extensive deformation as a result of the strain gradient plasticity due to the elastic plastic mismatch between the particle and the matrix. In this paper, the magnitudes of these two types of dislocations are calculated based on the dislocation plasticity. The dislocations are then converted to the respective strengths and allocated hierarchically to the matrix around the particle in the axisymmetric finite element unit cell model. the proposed method is shown to be very effective by performing finite element strength analysis of SiC p /Al2124 T4 composites that included ductile in the matrix and particle matrix decohesion. The predicted results for different particle sizes and volume fractions show that the length scale effect of the particle size obviously affects the strength and failure behavior of the particle reinforced metal matrix composites

Assessment of stress-strain data suitable for finite-element elastic--plastic analysis of shipping containers

International Nuclear Information System (INIS)

Rack, H.J.; Knorovsky, G.A.

1978-09-01

Stress-strain data which describes the influence of strain rate and temperature on the mechanical response of materials presently being used for light water reactor fuel shipping containers have been assembled. Selection of data has been limited to that which is suitable for use in finite-element elastic--plastic analysis of shipping containers (e.g., they must include complete material history profiles). Based on this information, recommendations have been made for further work which is required to complete the necessary data base

Consistency between analytical and finite element predictions for safety of cylindrical pressure vessels at higher temperatures

International Nuclear Information System (INIS)

Iancu, Otto Theodor

2014-01-01

The prediction of the plastic collapse load of cylindrical pressure vessels is very often made by using expensive Finite Element computations. The calculation of the collapse load requires an elastic-plastic material model and the consideration of non-linear geometry effects. The plastic collapse load causes overall structural instability and cannot be determined directly from a Finite Element analysis. In the present paper the plastic collapse load for a cylindrical pressure vessel is determined by an analytical method based on a linear elastic perfectly plastic material model. When plasticity occurs the material is considered to be incompressible and the tensor of plastic strains to be parallel to the stress deviator tensor. In this case the finite stress-strain relationships of Henkel can be used for calculating the pressure for which plastic flow occurs. The analytical results are completely confirmed by Finite Element predictions. (orig.)

Regularized finite element modeling of progressive failure in soils within nonlocal softening plasticity

Science.gov (United States)

Huang, Maosong; Qu, Xie; Lü, Xilin

2017-11-01

By solving a nonlinear complementarity problem for the consistency condition, an improved implicit stress return iterative algorithm for a generalized over-nonlocal strain softening plasticity was proposed, and the consistent tangent matrix was obtained. The proposed algorithm was embodied into existing finite element codes, and it enables the nonlocal regularization of ill-posed boundary value problem caused by the pressure independent and dependent strain softening plasticity. The algorithm was verified by the numerical modeling of strain localization in a plane strain compression test. The results showed that a fast convergence can be achieved and the mesh-dependency caused by strain softening can be effectively eliminated. The influences of hardening modulus and material characteristic length on the simulation were obtained. The proposed algorithm was further used in the simulations of the bearing capacity of a strip footing; the results are mesh-independent, and the progressive failure process of the soil was well captured.

Experimental and finite element analyses of plastic deformation behavior in vortex extrusion

International Nuclear Information System (INIS)

Shahbaz, M.; Pardis, N.; Kim, J.G.; Ebrahimi, R.; Kim, H.S.

2016-01-01

Vortex extrusion (VE) is a single pass severe plastic deformation (SPD) technique which can impose high strain values with almost uniform distribution within cross section of the processed material. This technique needs no additional facilities for installation on any conventional extrusion equipment. In this study the deformation behavior of material during VE is investigated and the results are compared with those of conventional extrusion (CE). These investigations include finite element analysis, visioplasticity, and microstructural characterization of the processed samples. The results indicate that the VE process can accumulate a higher strain value by applying an additional torsional deformation. The role of this additional deformation mode on the microstructural evolution of the VE sample is discussed and compared with the results obtained on the CE samples.

Elastic-Plastic J-Integral Solutions or Surface Cracks in Tension Using an Interpolation Methodology. Appendix C -- Finite Element Models Solution Database File, Appendix D -- Benchmark Finite Element Models Solution Database File

Science.gov (United States)

Allen, Phillip A.; Wells, Douglas N.

2013-01-01

No closed form solutions exist for the elastic-plastic J-integral for surface cracks due to the nonlinear, three-dimensional nature of the problem. Traditionally, each surface crack must be analyzed with a unique and time-consuming nonlinear finite element analysis. To overcome this shortcoming, the authors have developed and analyzed an array of 600 3D nonlinear finite element models for surface cracks in flat plates under tension loading. The solution space covers a wide range of crack shapes and depths (shape: 0.2 less than or equal to a/c less than or equal to 1, depth: 0.2 less than or equal to a/B less than or equal to 0.8) and material flow properties (elastic modulus-to-yield ratio: 100 less than or equal to E/ys less than or equal to 1,000, and hardening: 3 less than or equal to n less than or equal to 20). The authors have developed a methodology for interpolating between the goemetric and material property variables that allows the user to reliably evaluate the full elastic-plastic J-integral and force versus crack mouth opening displacement solution; thus, a solution can be obtained very rapidly by users without elastic-plastic fracture mechanics modeling experience. Complete solutions for the 600 models and 25 additional benchmark models are provided in tabular format.

FINITE ELEMENT ANALYSIS OF CONCRETE FILLER INFLUENCE ON DYNAMIC RIGIDITY OF HEAVY MACHINE TOOL PORTAL

Directory of Open Access Journals (Sweden)

Yu. V. Vasilevich

2016-01-01

Full Text Available Virtual testing of portal machine tool has been carried out with the help of finite elements method (FEM. Static, modal and harmonic analyses have been made for a heavy planer. The paper reveals influence of concrete filler on machine tool dynamic flexibility. A peculiar feature of the simulation is concrete filling of a high-level transverse beam. Such approach oes look a typical one for machine-tool industry. Concrete has been considered as generalized material in two variants. It has been established that concrete application provides approximately 3-fold increase in machine tool rigidity per each coordinate. In this regard it is necessary to arrange closure of rigidity contour by filling all the cavities inside of the portal. Modal FEA makes it possible to determine that concrete increases comparatively weakly (1.3–1.4-fold frequencies of resonance modes. Frequency of the lowest mode rises only from 30.25 to 42.86 Hz. The following most active whole-machine eigenmodes have been revealed in the paper: “Portal pecking”, “Parallelogram” and “Traverse pecking”. In order to restrain the last mode it is necessary to carry out concrete filling of the traverse, in particular. Frequency-response characteristics and curves of dynamic rigidity for a spindle have been plotted for 0–150 Hz interval while using harmonic FEM. It has been determined that concrete increases dynamic machine tool rigidity by 2.5–3.5-fold. The effect is obtained even in the case when weakly damping concrete (2 % is used. This is due to distribution of vibrational energy flow along concrete and along cast iron as well. Thus energy density and vibration amplitudes must decrease. The paper shows acceptability for internal reinforcement of high-level machine tool parts (for example, portal traverses and fillers are applied for this purpose. Traverse weighting is compensated by additional torsional, shear and bending rigidity. The machine tool obtains the

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

International Nuclear Information System (INIS)

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

2008-01-01

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

Comparison of finite element and fast Fourier transform crystal plasticity solvers for texture prediction

International Nuclear Information System (INIS)

Liu, B; Raabe, D; Roters, F; Eisenlohr, P; Lebensohn, R A

2010-01-01

We compare two full-field formulations, i.e. a crystal plasticity fast Fourier transform-based (CPFFT) model and the crystal plasticity finite element model (CPFEM) in terms of the deformation textures predicted by both approaches. Plane-strain compression of a 1024-grain ensemble is simulated with CPFFT and CPFEM to assess the models in terms of their predictions of texture evolution for engineering applications. Different combinations of final textures and strain distributions are obtained with the CPFFT and CPFEM models for this 1024-grain polycrystal. To further understand these different predictions, the correlation between grain rotations and strain gradients is investigated through the simulation of plane-strain compression of bicrystals. Finally, a study of the influence of the initial crystal orientation and the crystallographic neighborhood on grain rotations and grain subdivisions is carried out by means of plane-strain compression simulations of a 64-grain cluster

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.

Analysis of polyethylene terephthalate PET plastic bottle jointing system using finite element method (FEM)

Science.gov (United States)

Zaidi, N. A.; Rosli, Muhamad Farizuan; Effendi, M. S. M.; Abdullah, Mohamad Hariri

2017-09-01

For almost all injection molding applications of Polyethylene Terephthalate (PET) plastic was analyzed the strength, durability and stiffness of properties by using Finite Element Method (FEM) for jointing system of wood furniture. The FEM was utilized for analyzing the PET jointing system for Oak and Pine as wood based material of furniture. The difference pattern design of PET as wood jointing furniture gives the difference value of strength furniture itself. The results show the wood specimen with grooves and eclipse pattern design PET jointing give lower global estimated error is 28.90%, compare to the rectangular and non-grooves wood specimen of global estimated error is 63.21%.

Finite element discretization of Darcy's equations with pressure dependent porosity

KAUST Repository

Girault, Vivette; Murat, Franç ois; Salgado, Abner

2010-01-01

We consider the flow of a viscous incompressible fluid through a rigid homogeneous porous medium. The permeability of the medium depends on the pressure, so that the model is nonlinear. We propose a finite element discretization of this problem and

A new dedicated finite element for push-over analysis of reinforced concrete shear wall systems

Directory of Open Access Journals (Sweden)

Delal Doğru ORMANCI

2016-06-01

Full Text Available In this study, a finite element which has been analyzed based on anisotropic behavior of reinforced shear walls is developed. Element stiffness matrices were varied based on whether the element is in the tension or the compression zone of the cross-section. Nonlinear behavior of reinforced shear wall model is investigated under horizontal loads. This behavior is defined with a similar approach to plastic hinge assumption in frame structures that the finite element behaves lineer elastic between joints and plastic deformations are concentrated on joints as vertical plastic displacements. According to this acceptance, plastic behavior of reinforced shear wall occurs when the vertical strain reaches elastic strain limit. In the definition of finite element, displacement functions are chosen considering that the partition of shear walls just at floor levels, are enough for solution. Results of this study are compared with the solution obtained from a different computer programme and experimental results.

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.

On crack interaction effects of in-plane surface cracks using elastic and elastic-plastic finite element analyses

International Nuclear Information System (INIS)

Kim, Jong Min; Huh, Nam Su

2010-01-01

The crack-tip stress fields and fracture mechanics assessment parameters for a surface crack, such as the elastic stress intensity factor or the elastic-plastic J-integral, can be affected significantly by the adjacent cracks. Such a crack interaction effect due to multiple cracks can alter the fracture mechanics assessment parameters significantly. There are many factors to be considered, for instance the relative distance between adjacent cracks, the crack shape, and the loading condition, to quantify the crack interaction effect on the fracture mechanics assessment parameters. Thus, the current assessment codes on crack interaction effects (crack combination rules), including ASME Sec. XI, BS7910, British Energy R6 and API 579-1/ASME FFS-1, provide different rules for combining multiple surface cracks into a single surface crack. The present paper investigates crack interaction effects by evaluating the elastic stress intensity factor and the elastic-plastic J-integral of adjacent in-plane surface cracks in a plate through detailed 3-dimensional elastic and elastic-plastic finite element analyses. The effects on the fracture mechanics assessment parameters of the geometric parameters, the relative distance between two cracks, and the crack shape are investigated systematically. As for the loading condition, an axial tension is considered. Based on the finite element results, the acceptability of the crack combination rules provided in the existing guidance was investigated, and the relevant recommendations on a crack interaction for in-plane surface cracks are discussed. The present results can be used to develop more concrete guidance on crack interaction effects for crack shape characterization to evaluate the integrity of defective components

A parallel finite element method for the analysis of crystalline solids

DEFF Research Database (Denmark)

Sørensen, N.J.; Andersen, B.S.

1996-01-01

A parallel finite element method suitable for the analysis of 3D quasi-static crystal plasticity problems has been developed. The method is based on substructuring of the original mesh into a number of substructures which are treated as isolated finite element models related via the interface...... conditions. The resulting interface equations are solved using a direct solution method. The method shows a good speedup when increasing the number of processors from 1 to 8 and the effective solution of 3D crystal plasticity problems whose size is much too large for a single work station becomes possible....

3D finite element model of elastoplastic contact on the double sinus rough surface

International Nuclear Information System (INIS)

Hagege, H; Bouvier, S; Mazeran, P-E; Bigerelle, M

2011-01-01

One of the objectives in the field of tribology is to solve the mechanical stress-displacement problem involved by rough contacts. In our approach, the surface chosen is a 256-256 μm 2 3D sinusoidal shape (amplitude 4.5μm, wavelength 50μm) with an elastoplastic constitutive behaviour. The constitutive law combines isotropic and kinematic hardening and is experimentally identified from 316L steel sheets. The FEM deformable surface is crushed then uncrushed by a rigid flat surface: stresses, contact pressure and plastic cumulated strain are computed. We investigate the results sensitivity with respect to the level of in-plane refinement. At last, we conclude on some guidelines for 3D finite elements modelling of rough surfaces.

Rigid body displacement fields of an in-plane-deformable curved beam based on conventional strain definition

International Nuclear Information System (INIS)

Moon, Won Joo; Min, Oak Key; Kim, Yong Woo

1998-01-01

To improve the convergence and the accuracy of a finite element, the finite element has to describe not only displacement and stress distributions in a static analysis but also rigid body displacements. In this paper, we consider the in-plane-deformable curved beam element to understand the descriptive capability of rigid body displacements of a finite element. We derive the rigid body displacement fields of a single finite element under various essential boundary conditions when the nodal displacements are caused by the rigid body displacement. We also examine the rigid body displacement fields of a quadratic curved beam element by employing the reduced minimization theory

Elastoplastic finite element analysis for wet multidisc brake during lasting braking

Directory of Open Access Journals (Sweden)

Ji Zhanling

2015-01-01

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

Computational strain gradient crystal plasticity

DEFF Research Database (Denmark)

Niordson, Christian Frithiof; Kysar, Jeffrey W.

2014-01-01

A numerical method for viscous strain gradient crystal plasticity theory is presented, which incorporates both energetic and dissipative gradient effects. The underlying minimum principles are discussed as well as convergence properties of the proposed finite element procedure. Three problems...... of plane crystal plasticity are studied: pure shear of a single crystal between rigid platens as well as plastic deformation around cylindrical voids in hexagonal close packed and face centered cubic crystals. Effective in-plane constitutive slip parameters for plane strain deformation of specifically...... oriented face centered cubic crystals are developed in terms of the crystallographic slip parameters. The effect on geometrically necessary dislocation structures introduced by plastic deformation is investigated as a function of the ratio of void radius to plasticity length scale....

A finite element-based algorithm for rubbing induced vibration prediction in rotors

Science.gov (United States)

Behzad, Mehdi; Alvandi, Mehdi; Mba, David; Jamali, Jalil

2013-10-01

In this paper, an algorithm is developed for more realistic investigation of rotor-to-stator rubbing vibration, based on finite element theory with unilateral contact and friction conditions. To model the rotor, cross sections are assumed to be radially rigid. A finite element discretization based on traditional beam theories which sufficiently accounts for axial and transversal flexibility of the rotor is used. A general finite element discretization model considering inertial and viscoelastic characteristics of the stator is used for modeling the stator. Therefore, for contact analysis, only the boundary of the stator is discretized. The contact problem is defined as the contact between the circular rigid cross section of the rotor and “nodes” of the stator only. Next, Gap function and contact conditions are described for the contact problem. Two finite element models of the rotor and the stator are coupled via the Lagrange multipliers method in order to obtain the constrained equation of motion. A case study of the partial rubbing is simulated using the algorithm. The synchronous and subsynchronous responses of the partial rubbing are obtained for different rotational speeds. In addition, a sensitivity analysis is carried out with respect to the initial clearance, the stator stiffness, the damping parameter, and the coefficient of friction. There is a good agreement between the result of this research and the experimental result in the literature.

Hybrid-finite-element analysis of some nonlinear and 3-dimensional problems of engineering fracture mechanics

Science.gov (United States)

Atluri, S. N.; Nakagaki, M.; Kathiresan, K.

1980-01-01

In this paper, efficient numerical methods for the analysis of crack-closure effects on fatigue-crack-growth-rates, in plane stress situations, and for the solution of stress-intensity factors for arbitrary shaped surface flaws in pressure vessels, are presented. For the former problem, an elastic-plastic finite element procedure valid for the case of finite deformation gradients is developed and crack growth is simulated by the translation of near-crack-tip elements with embedded plastic singularities. For the latter problem, an embedded-elastic-singularity hybrid finite element method, which leads to a direct evaluation of K-factors, is employed.

Identification of material parameters for plasticity models: A comparative study on the finite element model updating and the virtual fields method

Science.gov (United States)

Martins, J. M. P.; Thuillier, S.; Andrade-Campos, A.

2018-05-01

The identification of material parameters, for a given constitutive model, can be seen as the first step before any practical application. In the last years, the field of material parameters identification received an important boost with the development of full-field measurement techniques, such as Digital Image Correlation. These techniques enable the use of heterogeneous displacement/strain fields, which contain more information than the classical homogeneous tests. Consequently, different techniques have been developed to extract material parameters from full-field measurements. In this study, two of these techniques are addressed, the Finite Element Model Updating (FEMU) and the Virtual Fields Method (VFM). The main idea behind FEMU is to update the parameters of a constitutive model implemented in a finite element model until both numerical and experimental results match, whereas VFM makes use of the Principle of Virtual Work and does not require any finite element simulation. Though both techniques proved their feasibility in linear and non-linear constitutive models, it is rather difficult to rank their robustness in plasticity. The purpose of this work is to perform a comparative study in the case of elasto-plastic models. Details concerning the implementation of each strategy are presented. Moreover, a dedicated code for VFM within a large strain framework is developed. The reconstruction of the stress field is performed through a user subroutine. A heterogeneous tensile test is considered to compare FEMU and VFM strategies.

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

International Nuclear Information System (INIS)

Soba, Alejandro; Denis, Alicia C.

1999-01-01

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

Analysis of elastic-plastic problems using edge-based smoothed finite element method

International Nuclear Information System (INIS)

Cui, X.Y.; Liu, G.R.; Li, G.Y.; Zhang, G.Y.; Sun, G.Y.

2009-01-01

In this paper, an edge-based smoothed finite element method (ES-FEM) is formulated for stress field determination of elastic-plastic problems using triangular meshes, in which smoothing domains associated with the edges of the triangles are used for smoothing operations to improve the accuracy and the convergence rate of the method. The smoothed Galerkin weak form is adopted to obtain the discretized system equations, and the numerical integration becomes a simple summation over the edge-based smoothing domains. The pseudo-elastic method is employed for the determination of stress field and Hencky's total deformation theory is used to define effective elastic material parameters, which are treated as field variables and considered as functions of the final state of stress fields. The effective elastic material parameters are then obtained in an iterative manner based on the strain controlled projection method from the uniaxial material curve. Some numerical examples are investigated and excellent results have been obtained demonstrating the effectivity of the present method.

Static Buckling Model Tests and Elasto-plastic Finite Element Analysis of a Pile in Layers with Various Thicknesses

Science.gov (United States)

Okajima, Kenji; Imai, Junichi; Tanaka, Tadatsugu; Iida, Toshiaki

Damage to piles in the liquefied ground is frequently reported. Buckling by the excess vertical load could be one of the causes of the pile damage, as well as the lateral flow of the ground and the lateral load at the pile head. The buckling mechanism is described as a complicated interaction between the pile deformation by the vertical load and the earth pressure change cased by the pile deformation. In this study, series of static buckling model tests of a pile were carried out in dried sand ground with various thickness of the layer. Finite element analysis was applied to the test results to verify the effectiveness of the elasto-plastic finite element analysis combining the implicit-explicit mixed type dynamic relaxation method with the return mapping method to the pile buckling problems. The test results and the analysis indicated the possibility that the buckling load of a pile decreases greatly where the thickness of the layer increases.

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.

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

International Nuclear Information System (INIS)

Rudolph, Juergen; Goetz, Andreas; Hilpert, Roland

2012-01-01

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

Stochastic higher order finite element elasto-plastic analysis of the necking phenomenon

Directory of Open Access Journals (Sweden)

Strąkowski Michał

2017-01-01

Full Text Available The principal goal of this work is to investigate an application of the stochastic perturbation technique of the 10th order in coupled thermo-elasto-plastic analysis of tension of the steel elastic bar exposed to fire with thermally dependent material characteristics. An ambient temperature, calculated from the fire curve after ISO 834-1, equivalent to the fire exposure of the steel structure is treated here as the input Gaussian random variable. It is uniquely defined by the constant mean value at outer surfaces of this element, where material parameters of the steel as Young modulus, yield strength, heat conductivity, capacity and thermal elongation are considered all as highly temperature-dependent. Computational implementation known as the Stochastic Finite Element Method is carried out with the use of the FEM system ABAQUS and computer algebra system MAPLE. It uses both polynomial and non-polynomial local response functions of stresses and displacements. The basic probabilistic characteristics of time-dependent structural response are determined (expectations, coefficients of variation, skewness and kurtosis and verified with classical Monte-Carlo simulation scheme and semi-analytical technique for input coefficient of variation not larger than 0.20. Finally, probabilistic convergence of all three methods versus increasing input uncertainty level is investigated.

Finite element modelling of fibre-reinforced brittle materials

NARCIS (Netherlands)

Kullaa, J.

1997-01-01

The tensile constitutive behaviour of fibre-reinforced brittle materials can be extended to two or three dimensions by using the finite element method with crack models. The three approaches in this study include the smeared and discrete crack concepts and a multi-surface plasticity model. The

Crystal Plasticity Finite Element Analysis of Loading-Unloading Behaviour in Magnesium Alloy Sheet

International Nuclear Information System (INIS)

Hama, Takayuki; Fujimoto, Hitoshi; Takuda, Hirohiko

2010-01-01

Magnesium alloy sheets exhibit strong inelastic response during unloading. In this study crystal plasticity finite element analysis of loading-unloading behaviour during uniaxial tension in a rolled magnesium alloy sheet was carried out, and the mechanism of this inelastic response was examined in detail in terms of macroscopic and mesoscopic deformations. The unloading behaviour obtained by the simulation was in good agreement with the experiment in terms of variation with stress of instantaneous tangent modulus during unloading. Variations of activities of each family of slip systems during the deformation showed that the activation of basal slip systems is the largest during unloading, and the slip direction during unloading is opposite from during loading. These results indicated that one of the factors of the inelastic behaviour during unloading is the fact that the basal slip systems are easily activated during unloading because of their low strengths.

Rigid-Plastic Post-Buckling Analysis of Columns and Quadratic Plates

DEFF Research Database (Denmark)

Jönsson, Jeppe

2008-01-01

the compressive load as a function of the transverse displacement. An estimate of the magnitude of the transverse displacement prior to the forming of the collapse mechanism is introduced into the compressive load function, determined by the virtual work equation, thereby revealing a qualified estimate...... yield lines accommodate differential rotations of rigid parts and the area “collapse” yield lines accommodate local area changes of the rigid parts thereby preserving compatibility of the rigid parts of a plate. The approach will be illustrated for rigid plastic column analysis and for a quadratic plate...

Strain gradient effects on cyclic plasticity

DEFF Research Database (Denmark)

Niordson, Christian Frithiof; Legarth, Brian Nyvang

2010-01-01

Size effects on the cyclic shear response are studied numerically using a recent higher order strain gradient visco-plasticity theory accounting for both dissipative and energetic gradient hardening. Numerical investigations of the response under cyclic pure shear and shear of a finite slab between...... rigid platens have been carried out, using the finite element method. It is shown for elastic–perfectly plastic solids how dissipative gradient effects lead to increased yield strength, whereas energetic gradient contributions lead to increased hardening as well as a Bauschinger effect. For linearly...... hardening materials it is quantified how dissipative and energetic gradient effects promote hardening above that of conventional predictions. Usually, increased hardening is attributed to energetic gradient effects, but here it is found that also dissipative gradient effects lead to additional hardening...

Plastic collapse of tubes submitted to a ring load by optimization

International Nuclear Information System (INIS)

Zouain, N.

1982-05-01

The limit analysis of a tube with finite lenght, made of a rigid - plastic material, is considered for the case of an internal load uniformely distributed in a cross section of the tube. The exact creep law is calculated for several qualitatively differents cases, namely different tube lenghts. The corresponding stress and collapse mechanisms are given so that they can be compared to the approximations developed here. The static and kinematic theorems on plastic collapse are used for establishing two numerical methods of resolution, specifically mathematical programming and finite element method. These mathematical methods are applied to collapse load for the considered tube. (E.G.) [pt

PLANS; a finite element program for nonlinear analysis of structures. Volume 2: User's manual

Science.gov (United States)

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

1977-01-01

The PLANS system, rather than being one comprehensive computer program, is a collection of finite element programs used for the nonlinear analysis of structures. This collection of programs evolved and is based on the organizational philosophy in which classes of analyses are treated individually based on the physical problem class to be analyzed. Each of the independent finite element computer programs of PLANS, with an associated element library, can be individually loaded and used to solve the problem class of interest. A number of programs have been developed for material nonlinear behavior alone and for combined geometric and material nonlinear behavior. The usage, capabilities, and element libraries of the current programs include: (1) plastic analysis of built-up structures where bending and membrane effects are significant, (2) three dimensional elastic-plastic analysis, (3) plastic analysis of bodies of revolution, and (4) material and geometric nonlinear analysis of built-up structures.

An efficient structural finite element for inextensible flexible risers

Science.gov (United States)

Papathanasiou, T. K.; Markolefas, S.; Khazaeinejad, P.; Bahai, H.

2017-12-01

A core part of all numerical models used for flexible riser analysis is the structural component representing the main body of the riser as a slender beam. Loads acting on this structural element are self-weight, buoyant and hydrodynamic forces, internal pressure and others. A structural finite element for an inextensible riser with a point-wise enforcement of the inextensibility constrain is presented. In particular, the inextensibility constraint is applied only at the nodes of the meshed arc length parameter. Among the virtues of the proposed approach is the flexibility in the application of boundary conditions and the easy incorporation of dissipative forces. Several attributes of the proposed finite element scheme are analysed and computation times for the solution of some simplified examples are discussed. Future developments aim at the appropriate implementation of material and geometric parameters for the beam model, i.e. flexural and torsional rigidity.

Finite Element Simulation of Temperature and Strain Distribution during Friction Stir Welding of AA2024 Aluminum Alloy

Science.gov (United States)

Jain, Rahul; Pal, Surjya Kanta; Singh, Shiv Brat

2017-02-01

Friction Stir Welding (FSW) is a solid state joining process and is handy for welding aluminum alloys. Finite Element Method (FEM) is an important tool to predict state variables of the process but numerical simulation of FSW is highly complex due to non-linear contact interactions between tool and work piece and interdependency of displacement and temperature. In the present work, a three dimensional coupled thermo-mechanical method based on Lagrangian implicit method is proposed to study the thermal history, strain distribution and thermo-mechanical process in butt welding of Aluminum alloy 2024 using DEFORM-3D software. Workpiece is defined as rigid-visco plastic material and sticking condition between tool and work piece is defined. Adaptive re-meshing is used to tackle high mesh distortion. Effect of tool rotational and welding speed on plastic strain is studied and insight is given on asymmetric nature of FSW process. Temperature distribution on the workpiece and tool is predicted and maximum temperature is found in workpiece top surface.

Finite element modelling of creep process - steady state stresses and strains

Directory of Open Access Journals (Sweden)

Sedmak Aleksandar S.

2014-01-01

Full Text Available Finite element modelling of steady state creep process has been described. Using an analogy of visco-plastic problem with a described procedure, the finite element method has been used to calculate steady state stresses and strains in 2D problems. An example of application of such a procedure have been presented, using real life problem - cylindrical pipe with longitudinal crack at high temperature, under internal pressure, and estimating its residual life, based on the C*integral evaluation.

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.

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

International Nuclear Information System (INIS)

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

2010-01-01

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

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%

Basic Strain Gradient Plasticity Theories with Application to Constrained Film Deformation

DEFF Research Database (Denmark)

Niordson, Christian Frithiof; Hutchinson, John W.

2011-01-01

films: the compression or extension of a finite layer joining rigid platens. Full elastic-plastic solutions are obtained for the same problem based on a finite element method devised for the new class of flow theories. Potential difficulties and open issues associated with the new class of flow theories......A family of basic rate-independent strain gradient plasticity theories is considered that generalize conventional J(2) deformation and flow theories of plasticity to include a dependence on strain gradients in a simple way. The theory builds on three recent developments: the work of Gudmundson (J....... Mech. Phys. Solids 52 (2004), 1379-1406) and Gurtin and Anand (J. Mech. Phys. Solids 57 (2009), 405-421), proposing constitutive relations for flow theories consistent with requirements of positive plastic dissipation; the work of Fleck and Willis (J. Mech. Phys. Solids 57 (2009), 161-177 and 1045...

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

ACCEPT: a three-dimensional finite element program for large deformation elastic-plastic-creep analysis of pressurized tubes (LWBR/AWBA Development Program)

International Nuclear Information System (INIS)

Hutula, D.N.; Wiancko, B.E.

1980-03-01

ACCEPT is a three-dimensional finite element computer program for analysis of large-deformation elastic-plastic-creep response of Zircaloy tubes subjected to temperature, surface pressures, and axial force. A twenty-mode, tri-quadratic, isoparametric element is used along with a Zircaloy materials model. A linear time-incremental procedure with residual force correction is used to solve for the time-dependent response. The program features an algorithm which automatically chooses the time step sizes to control the accuracy and numerical stability of the solution. A contact-separation capability allows modeling of interaction of reactor fuel rod cladding with fuel pellets or external supports

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.

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

Science.gov (United States)

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

1998-02-01

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

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.

Finite element analysis of the cyclic indentation of bilayer enamel

International Nuclear Information System (INIS)

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

2014-01-01

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

Finite element analysis of the cyclic indentation of bilayer enamel

Science.gov (United States)

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

2014-04-01

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

Experimental and Finite Element Analysis of Asymmetric Rolling of 6061 Aluminum Alloy Using Two-Scale Elasto-Plastic Constitutive Relation

Directory of Open Access Journals (Sweden)

Wronski M.

2017-12-01

Full Text Available The goal of this work was theoretical and experimental study of micro- and macroscopic mechanical fields of 6061 aluminum alloy induced by the asymmetric rolling process. Two-scale constitutive law was used by implementing an elasto-plastic self-consistent scheme into the Finite Element code (ABAQUS/Explicit. The model was applied to study the asymmetric rolling. Such a deformation process induces heterogeneous mechanical fields that were reproduced by the model thanks to the crystallographic nature of constitutive law used. The studied material was processed, at room temperature, in one rolling pass to 36% reduction. The resulting material modifications were compared with predictions of the two-scale model. Namely, the calculated textures were compared with experimental ones determined by X-ray diffraction. Especially, detailed quantitative analysis of texture variation across the sample thickness was done. The influence of this texture variation on plastic anisotropy was studied. The advantages of asymmetric rolling process over symmetric one were identified. The main benefits are a nearly homogeneous crystallographic texture, reduced rolling normal forces and homogenization of plastic anisotropy through the sample thickness.

Finite element historical deformation analysis in piecewise linear plasticity by mathematical programming

International Nuclear Information System (INIS)

De Donato, O.; Parisi, M.A.

1977-01-01

When loads increase proportionally beyond the elastic limit in the presence of elastic-plastic piecewise-linear constitutive laws, the problem of finding the whole evolution of the plastic strain and displacements of structures was recently shown to be amenable to a parametric linear complementary problem (PLCP) in which the parameter is represented by the load factor, the matrix is symmetric positive definite or at least semi-definite (for perfect plasticity) and the variables with a direct mechanical meaning are the plastic multipliers. With reference to plane trusses and frames with elastic-plastic linear work-hardening material behaviour numerical solutions were also fairly efficiently obtained using a recent mathematical programming algorithm (due to R.W. Cottle) which is able to provide the whole deformation history of the structure and, at the same time to rule out local unloadings along the given proportional loading process by means of 'a priori' checks carried out before each pivotal step of the procedure. Hence it becomes possible to use the holonomic (reversible, path-independent) constitutive laws in finite terms and to benefit by all the relevant numerical and computational advantages despite the non-holonomic nature of plastic behaviour. In the present paper the method of solution is re-examined in view to overcome an important drawback of the algorithm deriving from the size of PLCP fully populated matrix when structural problems with large number of variables are considered and, consequently, the updating, the storing or, generally, the handling of the current tableau may become prohibitive. (Auth.)

APPLICATION OF RIGID LINKS IN STRUCTURAL DESIGN MODELS

Directory of Open Access Journals (Sweden)

Sergey Yu. Fialko

2017-09-01

Full Text Available A special finite element modelling rigid links is proposed for the linear static and buckling analysis. Unlike the classical approach based on the theorems of rigid body kinematics, the proposed approach preserves the similarity between the adjacency graph for a sparse matrix and the adjacency graph for nodes of the finite element model, which allows applying sparse direct solvers more effectively. Besides, the proposed approach allows significantly reducing the number of nonzero entries in the factored stiffness matrix in comparison with the classical one, which greatly reduces the duration of the solution. For buckling problems of structures containing rigid bodies, this approach gives correct results. Several examples demonstrate its efficiency.

Generalized finite elements

International Nuclear Information System (INIS)

Wachspress, E.

2009-01-01

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

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

Tube Bulge Process : Theoretical Analysis and Finite Element Simulations

International Nuclear Information System (INIS)

Velasco, Raphael; Boudeau, Nathalie

2007-01-01

This paper is focused on the determination of mechanics characteristics for tubular materials, using tube bulge process. A comparative study is made between two different models: theoretical model and finite element analysis. The theoretical model is completely developed, based first on a geometrical analysis of the tube profile during bulging, which is assumed to strain in arc of circles. Strain and stress analysis complete the theoretical model, which allows to evaluate tube thickness and state of stress, at any point of the free bulge region. Free bulging of a 304L stainless steel is simulated using Ls-Dyna 970. To validate FE simulations approach, a comparison between theoretical and finite elements models is led on several parameters such as: thickness variation at the free bulge region pole with bulge height, tube thickness variation with z axial coordinate, and von Mises stress variation with plastic strain. Finally, the influence of geometrical parameters deviations on flow stress curve is observed using analytical model: deviations of the tube outer diameter, its initial thickness and the bulge height measurement are taken into account to obtain a resulting error on plastic strain and von Mises stress

Crystal plasticity finite element analysis of deformation behaviour in SAC305 solder joint

Science.gov (United States)

Darbandi, Payam

Due to the awareness of the potential health hazards associated with the toxicity of lead (Pb), actions have been taken to eliminate or reduce the use of Pb in consumer products. Among those, tin (Sn) solders have been used for the assembly of electronic systems. Anisotropy is of significant importance in all structural metals, but this characteristic is unusually strong in Sn, making Sn based solder joints one of the best examples of the influence of anisotropy. The effect of anisotropy arising from the crystal structure of tin and large grain microstructure on the microstructure and the evolution of constitutive responses of microscale SAC305 solder joints is investigated. Insights into the effects of key microstructural features and dominant plastic deformation mechanisms influencing the measured relative activity of slip systems in SAC305 are obtained from a combination of optical microscopy, orientation imaging microscopy (OIM), slip plane trace analysis and crystal plasticity finite element (CPFE) modeling. Package level SAC305 specimens were subjected to shear deformation in sequential steps and characterized using optical microscopy and OIM to identify the activity of slip systems. X-ray micro Laue diffraction and high energy monochromatic X-ray beam were employed to characterize the joint scale tensile samples to provide necessary information to be able to compare and validate the CPFE model. A CPFE model was developed that can account for relative ease of activating slip systems in SAC305 solder based upon the statistical estimation based on correlation between the critical resolved shear stress and the probability of activating various slip systems. The results from simulations show that the CPFE model developed using the statistical analysis of activity of slip system not only can satisfy the requirements associated with kinematic of plastic deformation in crystal coordinate systems (activity of slip systems) and global coordinate system (shape changes

21 CFR 178.3790 - Polymer modifiers in semirigid and rigid vinyl chloride plastics.

Science.gov (United States)

2010-04-01

... 21 Food and Drugs 3 2010-04-01 2009-04-01 true Polymer modifiers in semirigid and rigid vinyl...: ADJUVANTS, PRODUCTION AIDS, AND SANITIZERS Certain Adjuvants and Production Aids § 178.3790 Polymer modifiers in semirigid and rigid vinyl chloride plastics. The polymers identified in paragraph (a) of this...

Elasto-plastic impact of hemispherical shell impacting on hard rigid sphere

Science.gov (United States)

Raftopoulos, D. D.; Spicer, A. L.

1976-01-01

An analysis of plastic stress waves for cylindrical metallic projectile in impact is extended to an analysis of a hemispherical shell suffereing plastic deformation during the process of impact. It is assumed that the hemispherical shell with a prescribed launch velocity impinges a fixed rigid sphere of diameter equal to the internal diameter of the shell. The dynamic biaxial state of stress present in the shell during deformation is investigated. The analysis is valuable for studying the state of stress during large plastic deformation of a hemispherical shell.

Simulation of plastic displacement distributions for multistory shear frames excited by Gaussian white noise

DEFF Research Database (Denmark)

Lazarov, Boyan Stefanov; Ditlevsen, Ove Dalager

2003-01-01

The object of study is a stationary Gaussian white noise excited MDOF linear elastic, ideal plastic, linearly damped, statically determinate oscillator with several potential elements of ideal plastic yielding. Specifically the study is exemplified for a plane multistory shear frame with rigid...... traverses where all the connecting columns except the columns in one or more of the bottom floors have finite symmetrical yield limits. The white noise excitation acts on the mass of the first floor making the movement of the elastic bottom floors simulate a ground motion that interacts with the structure...

Simulation by Slepian method of plastic displacements of Gaussian process excited multistory shear frame

DEFF Research Database (Denmark)

Lazarov, Boyan Stefanov; Ditlevsen, Ove Dalager

2004-01-01

The object of study is a stationary Gaussian white noise excited multi-degree-of-freedom (MDOF) linear elastic, ideal plastic, linearly damped, statically determinate oscillator with several potential elements of ideal plastic yielding. Specifically the study is exemplified for a plane multistory...... shear frame with rigid traverses where all the connecting columns except the columns in one or more of the bottom floors have finite symmetrical yield limits. The white noise excitation acts on the mass of the first floor making the movement of the elastic bottom floors simulate a ground motion...

Stability Analysis of Anchored Soil Slope Based on Finite Element Limit Equilibrium Method

Directory of Open Access Journals (Sweden)

Rui Zhang

2016-01-01

Full Text Available Under the condition of the plane strain, finite element limit equilibrium method is used to study some key problems of stability analysis for anchored slope. The definition of safe factor in slices method is generalized into FEM. The “true” stress field in the whole structure can be obtained by elastic-plastic finite element analysis. Then, the optimal search for the most dangerous sliding surface with Hooke-Jeeves optimized searching method is introduced. Three cases of stability analysis of natural slope, anchored slope with seepage, and excavation anchored slope are conducted. The differences in safety factor quantity, shape and location of slip surface, anchoring effect among slices method, finite element strength reduction method (SRM, and finite element limit equilibrium method are comparatively analyzed. The results show that the safety factor given by the FEM is greater and the unfavorable slip surface is deeper than that by the slice method. The finite element limit equilibrium method has high calculation accuracy, and to some extent the slice method underestimates the effect of anchor, and the effect of anchor is overrated in the SRM.

Finite element method for starved hydrodynamic lubrication with film separation and free surface effects

DEFF Research Database (Denmark)

Poulios, Konstantinos; Vølund, Anders; Klit, Peder

2018-01-01

This paper proposes a numerical method for determining the evolution of lubricant film thickness and pressure in partially and fully flooded regions of a hydrodynamic contact between two non-conformal rigid surfaces. The proposed method accounts for the classical Reynolds equation within the fully...... thickness is zero. Both pressure and film thickness fields are considered as unknowns to solve for in each time step and they are approximated through quadratic B-spline finite elements. The geometry of the gap between the rigid surfaces delimiting the lubricant is accounted for in the form of a unilateral...

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

International Nuclear Information System (INIS)

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

1981-01-01

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

Rigid-Plastic Approximations for Predicting Plastic Deformation of Cylindrical Shells Subject to Dynamic Loading

Directory of Open Access Journals (Sweden)

Michelle S. Hoo Fatt

1996-01-01

Full Text Available A theoretical approach was developed for predicting the plastic deformation of a cylindrical shell subject to asymmetric dynamic loads. The plastic deformation of the leading generator of the shell is found by solving for the transverse deflections of a rigid-plastic beam/string-on-foundation. The axial bending moment and tensile force in the beam/string are equivalent to the longitudinal bending moments and membrane forces of the shell, while the plastic foundation force is equivalent to the shell circumferential bending moment and membrane resistances. Closed-form solutions for the transient and final deformation profile of an impulsive loaded shell when it is in a “string” state were derived using the eigenfunction expansion method. These results were compared to DYNA 3D predictions. The analytical predictions of the transient shell and final centerline deflections were within 25% of the DYNA 3D results.

Computational Strain Gradient Crystal Plasticity

DEFF Research Database (Denmark)

Niordson, Christian Frithiof; Kysar, Jeffrey W.

2011-01-01

A model for strain gradient crystal visco-plasticity is formulated along the lines proposed by Fleck andWillis (2009) for isotropic plasticity. Size-effects are included in the model due to the addition of gradient terms in both the free energy as well as through a dissipation potential. A finite...... element solution method is presented, which delivers the slip-rate field and the velocity-field based on two minimum principles. Some plane deformation problems relevant for certain specific orientations of a face centered cubic crystal under plane loading conditions are studied, and effective in......-plane parameters are developed based on the crystallographic properties of the material. The problem of cyclic shear of a single crystal between rigid platens is studied as well as void growth of a cylindrical void....

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.

Stability analysis and finite element simulations of superplastic forming in the presence of hydrostatic pressure

Science.gov (United States)

Nazzal, M. A.

2018-04-01

It is established that some superplastic materials undergo significant cavitation during deformation. In this work, stability analysis for the superplastic copper based alloy Coronze-638 at 550 °C based on Hart's definition of stable plastic deformation and finite element simulations for the balanced biaxial loading case are carried out to study the effects of hydrostatic pressure on cavitation evolution during superplastic forming. The finite element results show that imposing hydrostatic pressure yields to a reduction in cavitation growth.

49 CFR 178.925 - Standards for rigid plastic Large Packagings.

Science.gov (United States)

2010-10-01

... manufacture of the tested design type, retesting may be omitted if changes in the carbon black content, the... or chemical properties of the material of construction. (3) No used material other than production residues or regrind from the same manufacturing process may be used in the manufacture of rigid plastic...

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

A finite difference method for off-fault plasticity throughout the earthquake cycle

Science.gov (United States)

Erickson, Brittany A.; Dunham, Eric M.; Khosravifar, Arash

2017-12-01

We have developed an efficient computational framework for simulating multiple earthquake cycles with off-fault plasticity. The method is developed for the classical antiplane problem of a vertical strike-slip fault governed by rate-and-state friction, with inertial effects captured through the radiation-damping approximation. Both rate-independent plasticity and viscoplasticity are considered, where stresses are constrained by a Drucker-Prager yield condition. The off-fault volume is discretized using finite differences and tectonic loading is imposed by displacing the remote side boundaries at a constant rate. Time-stepping combines an adaptive Runge-Kutta method with an incremental solution process which makes use of an elastoplastic tangent stiffness tensor and the return-mapping algorithm. Solutions are verified by convergence tests and comparison to a finite element solution. We quantify how viscosity, isotropic hardening, and cohesion affect the magnitude and off-fault extent of plastic strain that develops over many ruptures. If hardening is included, plastic strain saturates after the first event and the response during subsequent ruptures is effectively elastic. For viscoplasticity without hardening, however, successive ruptures continue to generate additional plastic strain. In all cases, coseismic slip in the shallow sub-surface is diminished compared to slip accumulated at depth during interseismic loading. The evolution of this slip deficit with each subsequent event, however, is dictated by the plasticity model. Integration of the off-fault plastic strain from the viscoplastic model reveals that a significant amount of tectonic offset is accommodated by inelastic deformation ( ∼ 0.1 m per rupture, or ∼ 10% of the tectonic deformation budget).

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

Application of the Single Hardening Model in the Finite Element Program ABAQUS

DEFF Research Database (Denmark)

Jakobsen, Kim Parsberg

model, developed by Lade and Kim (Kim & Lade 1988, Lade & Kim 1988a, Lade & Kim 1988b) is implemented as a user defined material module, UMAT, in the commercial finite element program, ABAQUS. The advantages of the Single Hardening Model Iie in its ability to predict elastic and plastic displacements...

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

International Nuclear Information System (INIS)

Nakamachi, Eiji

2005-01-01

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

Three-dimensional finite element analysis of residual magnetic field for ferromagnets under early damage

International Nuclear Information System (INIS)

Yao, Kai; Shen, Kai; Wang, Zheng-Dao; Wang, Yue-Sheng

2014-01-01

In this study, 3D finite element analysis is presented by calculating the residual magnetic field signals of ferromagnets under the plastic deformation. The contour maps of tangential and normal RMF gradients are given, and the 3D effect is discussed. The results show that the tangential peak–peak amplitude and normal peak–vale amplitude are remarkably different in 2D and 3D simulations, but the tangential peak–peak width and normal peak–vale width are similar. Moreover, some key points are capable of capturing the plastic-zone shape, especially when the lift-off is small enough. The present study suggests an effective defect identification method with Metal magnetic memory (MMM) technique. - Highlights: • Three-dimensional (3D) finite element analysis is presented by calculating the residual magnetic field signals of ferromagnets under the plastic deformation. • The contour maps of gradients of the tangential and normal residual magnetic fields are given, and the 3D effect is discussed. • The present study suggests an effective defect identification method with metal magnetic memory technique

SAFE-PLANE, Stress Analysis of Planar Structure by Finite Elements Method

International Nuclear Information System (INIS)

Cornell, D.C.; Reich, Morris

1967-01-01

1 - Description of problem or function: SAFE-PLANE is applied to two- dimensional structures of arbitrary geometry under in-plane loads. Either plane stress or plane strain conditions may be imposed. Mechanical and thermal loads are permitted. 2 - Method of solution: The finite-element method is used to construct a mathematical model by assembling discrete elements. The total potential energy of the structure is determined and subsequently minimized by iteration on components of the displacement field until static equilibrium of the structure is attained. Strains and stresses are computed from the resulting displacements. 3 - Restrictions on the complexity of the problem: Multi-material structures with varying rigidities converge very slowly. Not valid for incompressible materials. Maximum number of nodal points = 675. Maximum number of elements = 1350

Rigid inclusions-Comparison between analytical and numerical methods

International Nuclear Information System (INIS)

Gomez Perez, R.; Melentijevic, S.

2014-01-01

This paper compares different analytical methods for analysis of rigid inclusions with finite element modeling. First of all, the load transfer in the distribution layer is analyzed for its different thicknesses and different inclusion grids to define the range between results obtained by analytical and numerical methods. The interaction between the soft soil and the inclusion in the estimation of settlements is studied as well. Considering different stiffness of the soft soil, settlements obtained analytical and numerically are compared. The influence of the soft soil modulus of elasticity on the neutral point depth was also performed by finite elements. This depth has a great importance for the definition of the total length of rigid inclusion. (Author)

Steam generator tube rupture simulation using extended finite element method

Energy Technology Data Exchange (ETDEWEB)

Mohanty, Subhasish, E-mail: smohanty@anl.gov; Majumdar, Saurin; Natesan, Ken

2016-08-15

Highlights: • Extended finite element method used for modeling the steam generator tube rupture. • Crack propagation is modeled in an arbitrary solution dependent path. • The FE model is used for estimating the rupture pressure of steam generator tubes. • Crack coalescence modeling is also demonstrated. • The method can be used for crack modeling of tubes under severe accident condition. - Abstract: A steam generator (SG) is an important component of any pressurized water reactor. Steam generator tubes represent a primary pressure boundary whose integrity is vital to the safe operation of the reactor. SG tubes may rupture due to propagation of a crack created by mechanisms such as stress corrosion cracking, fatigue, etc. It is thus important to estimate the rupture pressures of cracked tubes for structural integrity evaluation of SGs. The objective of the present paper is to demonstrate the use of extended finite element method capability of commercially available ABAQUS software, to model SG tubes with preexisting flaws and to estimate their rupture pressures. For the purpose, elastic–plastic finite element models were developed for different SG tubes made from Alloy 600 material. The simulation results were compared with experimental results available from the steam generator tube integrity program (SGTIP) sponsored by the United States Nuclear Regulatory Commission (NRC) and conducted at Argonne National Laboratory (ANL). A reasonable correlation was found between extended finite element model results and experimental results.

Steam generator tube rupture simulation using extended finite element method

International Nuclear Information System (INIS)

Mohanty, Subhasish; Majumdar, Saurin; Natesan, Ken

2016-01-01

Highlights: • Extended finite element method used for modeling the steam generator tube rupture. • Crack propagation is modeled in an arbitrary solution dependent path. • The FE model is used for estimating the rupture pressure of steam generator tubes. • Crack coalescence modeling is also demonstrated. • The method can be used for crack modeling of tubes under severe accident condition. - Abstract: A steam generator (SG) is an important component of any pressurized water reactor. Steam generator tubes represent a primary pressure boundary whose integrity is vital to the safe operation of the reactor. SG tubes may rupture due to propagation of a crack created by mechanisms such as stress corrosion cracking, fatigue, etc. It is thus important to estimate the rupture pressures of cracked tubes for structural integrity evaluation of SGs. The objective of the present paper is to demonstrate the use of extended finite element method capability of commercially available ABAQUS software, to model SG tubes with preexisting flaws and to estimate their rupture pressures. For the purpose, elastic–plastic finite element models were developed for different SG tubes made from Alloy 600 material. The simulation results were compared with experimental results available from the steam generator tube integrity program (SGTIP) sponsored by the United States Nuclear Regulatory Commission (NRC) and conducted at Argonne National Laboratory (ANL). A reasonable correlation was found between extended finite element model results and experimental results.

Strain Localization during Equal-Channel Angular Pressing Analyzed by Finite Element Simulations

Directory of Open Access Journals (Sweden)

Tobias Daniel Horn

2018-01-01

Full Text Available Equal-Channel Angular Pressing (ECAP is a method used to introduce severe plastic deformation into a metallic billet without changing its geometry. In special cases, strain localization occurs and a pattern consisting of regions with high and low deformation (so-called shear and matrix bands can emerge. This paper studies this phenomenon numerically adopting two-dimensional finite element simulations of one ECAP pass. The mechanical behavior of aluminum is modeled using phenomenological plasticity theory with isotropic or kinematic hardening. The effects of the two different strain hardening types are investigated numerically by systematic parameter studies: while isotropic hardening only causes minor fluctuations in the plastic strain fields, a material with high initial hardening rate and sufficient strain hardening capacity can exhibit pronounced localized deformation after ECAP. The corresponding finite element simulation results show a regular pattern of shear and matrix bands. This result is confirmed experimentally by ECAP-processing of AA6060 material in a severely cold worked condition, where microstructural analysis also reveals the formation of shear and matrix bands. Excellent agreement is found between the experimental and numerical results in terms of shear and matrix band width and length scale. The simulations provide additional insights regarding the evolution of the strain and stress states in shear and matrix bands.

Analytical Study of Common Rigid Steel Connections under the Effect of Heat

Directory of Open Access Journals (Sweden)

Rohola Rahnavard

2014-01-01

Full Text Available One of the most important members of steel structure’s connection region is beam-to-column connection. Rigid connection in steel moment frame has special role in the behavior of these structures and the fire resistance of these connections can be important. In this paper the behaviors of three common types of rigid connections in Iran under the effect of heat were studied by the use of numerical finite element methods through ABAQUS software. The models were verified by the use of an experimental model through elastic and plastic amplitudes up to collapse and during numerical results, and the effect of large deformation in the nonlinear region has also been considered. The results show that the connection with the end plate had a better performance against heat than other connections. Also reduced stiffness and lateral buckling in this connection were less than other connections.

Model tests and elasto-plastic finite element analysis on multicavity type PCRV

International Nuclear Information System (INIS)

Nojiri, Y.; Yamazaki, M.; Kotani, K.; Matsuzaki, Y.

1978-01-01

Multicavity type PCRV models were tested to investigate elastic stress distributions, cracking and failure mode of the models, and to determine the adequacy and relative accuracy of finite element structural analyses. The behavior of the models under pressure was investigated, and it was found that the predictions of the analyses showed a good agreement with the test results

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

Directory of Open Access Journals (Sweden)

M.A. Najafgholipour

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

Finite 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

Investigation of Apple Vibration Characteristics Using Finite Element Modal Analysis

Directory of Open Access Journals (Sweden)

R Mirzaei

2013-02-01

Full Text Available The most important quality indicator of fruits is the flesh firmness which is well correlated to their young’s modulus. In this research variation of vibration characteristics (shape modes, natural frequency of apple due to change of material characteristics (density, young's models, Poisson ratio and apple volume was investigated using Finite Element simulation. An image processing technique was used to obtain an unsymmetrical and non-spherical geometric model of apple. The exact three-dimensional shape of the fruit was created by determining the coordinates of apple surface and forming uneven rotational curvatures. Modal analysis with no boundary constraints has been applied. The first 20 Eigen frequencies and the corresponding mode shape were determined. Six rigid body modes possess zero resonant frequency which is related to the degree of freedom of a rigid body in space indicated the validity of finite element model. The modal analysis results showed that resonant frequency increased by increasing young's modulus of the fruit, while it decreased by increasing apple density. First mode torsion has a mean resonant frequency of 584 Hz. Variations of natural frequency due to change in young's modulus, density, and Poisson ratio were 80%, 11% and 4%, respectively. Coefficient of variation of resonant frequency in response to changing young's modulus was 2-3 times of that of density which shows the greatest effect of young modulus changes on natural frequency of fruits. Consequently with determination of fruits' natural frequency, their young modulus and firmness can be estimated.

On fracture in finite strain gradient plasticity

DEFF Research Database (Denmark)

Martínez Pañeda, Emilio; Niordson, Christian Frithiof

2016-01-01

In this work a general framework for damage and fracture assessment including the effect of strain gradients is provided. Both mechanism-based and phenomenological strain gradient plasticity (SGP) theories are implemented numerically using finite deformation theory and crack tip fields are invest......In this work a general framework for damage and fracture assessment including the effect of strain gradients is provided. Both mechanism-based and phenomenological strain gradient plasticity (SGP) theories are implemented numerically using finite deformation theory and crack tip fields...... are investigated. Differences and similarities between the two approaches within continuum SGP modeling are highlighted and discussed. Local strain hardening promoted by geometrically necessary dislocations (GNDs) in the vicinity of the crack leads to much higher stresses, relative to classical plasticity...... in the multiple parameter version of the phenomenological SGP theory. Since this also dominates the mechanics of indentation testing, results suggest that length parameters characteristic of mode I fracture should be inferred from nanoindentation....

A finite deformation theory of higher-order gradient crystal plasticity

DEFF Research Database (Denmark)

Kuroda, Mitsutoshi; Tvergaard, Viggo

2008-01-01

crystal plasticity that is based on an assumption of the existence of higher-order stresses. Furthermore, a boundary-value problem for simple shear of a constrained thin strip is studied numerically, and some characteristic features of finite deformation are demonstrated through a comparison to a solution......For higher-order gradient crystal plasticity, a finite deformation formulation is presented. The theory does not deviate much from the conventional crystal plasticity theory. Only a back stress effect and additional differential equations for evolution of the geometrically necessary dislocation...

Frost heave modelling of buried pipelines using non-linear Fourier finite elements

International Nuclear Information System (INIS)

Wan, R. G.; You, R.

1998-01-01

Numerical analysis of the response of a three-dimensional soil-pipeline system in a freezing environment using non-linear Fourier finite elements was described as an illustration of the effectiveness of this technique in analyzing plasticity problems. Plastic deformations occur when buried pipeline is under the action of non-uniform frost heave. The three-dimensional frost heave which develops over time including elastoplastic deformations of the soil and pipe are computed. The soil heave profile obtained in the numerical analysis was consistent with experimental findings for similar configurations. 8 refs., 8 figs

Finite element analysis of the combined fine blanking and extrusion process

Science.gov (United States)

Zheng, Peng-Fei

The combined fine blanking and extrusion process is such a metal forming process that fine blanking and forward extrusion are carried out on sheet metal material at the same time. There are two typical characteristics in this process, one is the fine blanking whose deformation mechanism is different from conventional blanking; the other is the sheet metal extrusion, which is different from the conventional extrusion. Even though fine blanking has been used in industry for many years, only limited literature can be found which deals with the theoretical analysis of it. On the other hand, no publications on the theoretical analysis of the sheet metal extrusion have been found. Intensive work should be carried out to reveal the mechanism of both fine blanking process and sheet metal extrusion process, and further the combined fine blanking and extrusion process. The scope of this thesis is to study the mechanics of fine blanking, sheet metal extrusion, and combined fine blanking and extrusion process one by one with the rigid-plastic finite element method. All of above processes are typical unsteady ones, especially the fine blanking process in which extremely severe and localized deformation occurs. Therefore, commercial programs can not be used to solve these problems up till now. Owing to this reason, a rigid-plastic finite element program was developed for simulating these processes where remeshing and mesh tracing techniques as well as the golden section method were adopted according to the characteristics of these processes in this thesis. Moreover, a permissible kinematic velocity field was adopted as the initial velocity field for simulating extrusion process successfully. Results from the simulation included the distorted mesh, the field of material flow, the stress and the strain distributions at various moments of deformation. Results under different deformation conditions such as different blanking clearances, different diameters of the extrusion punch and

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

Global finite-time attitude stabilization for rigid spacecraft in the exponential coordinates

Science.gov (United States)

Shi, Xiao-Ning; Zhou, Zhi-Gang; Zhou, Di

2018-06-01

This paper addresses the global finite-time attitude stabilisation problem on the special orthogonal group (SO(3)) for a rigid spacecraft via homogeneous feedback approach. Considering the topological and geometric properties of SO(3), the logarithm map is utilised to transform the stabilisation problem on SO(3) into the one on its associated Lie algebra (?). A model-independent discontinuous state feedback plus dynamics compensation scheme is constructed to achieve the global finite-time attitude stabilisation in a coordinate-invariant way. In addition, to address the absence of angular velocity measurements, a sliding mode observer is proposed to reconstruct the unknown angular velocity information within finite time. Then, an observer-based finite-time output feedback control strategy is obtained. Numerical simulations are finally performed to demonstrate the effectiveness of the proposed finite-time controllers.

SAFE-AXISYM, Stress Analysis of Axisymmetric Composite Structure by Finite Elements Method

International Nuclear Information System (INIS)

Cornell, D.C.

1967-01-01

1 - Nature of physical problem solved: SAFE-AXISYM is a program for the analysis of multi-material axisymmetric composite structures. It is designed for the analysis of heterogeneous structures such as reinforced and/or prestressed concrete vessels. The structure is assumed to be linearly elastic, and only bodies of revolution subjected to axisymmetric loading can be treated. 2 - Method of solution: SAFE-AXISYM uses a finite element method with a modified Gauss-Seidel iteration scheme. A reference grid subdivides the structure into ring-like small, finite elements, the vertices of which are called nodes. The grid may be generated by hand, by the computer or by a combination of the two methods. Each node has two degrees of freedom, translation in the and in the axial direction. Both zero and non-zero fixed displacement constraints may be assumed, and the loading condition may be mechanical and/or thermal. 3 - Restrictions on the complexity of the problem: Multi-material structures with varying rigidities converge very slowly. Not valid for incompressible materials. Maximum number of nodes = 475. Maximum number of elements = 1100

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.

Rotational Stiffness of Precast Beam-Column Connection using Finite Element Method

Science.gov (United States)

Hashim, N.; Agarwal, J.

2018-04-01

Current design practice in structural analysis is to assume the connection as pinned or rigid, however this cannot be relied upon for safety against collapse because during services the actual connection reacts differently where the connection has rotated in relevance. This situation may lead to different reactions and consequently affect design results and other frame responses. In precast concrete structures, connections play an important part in ensuring the safety of the whole structure. Thus, investigates on the actual connection behavior by construct the moment-rotation relationship is significant. Finite element (FE) method is chosen for modeling a 3-dimensional beam-column connection. The model is built in symmetry to reduce analysis time. Results demonstrate that precast billet connection is categorized as semi-rigid connection with Sini of 23,138kNm/rad. This is definitely different from the assumption of pinned or rigid connection used in design practice. Validation were made by comparing with mathematical equation and small differences were achieved that led to the conclusion where precast billet connection using FE method is acceptable.

Hybrid finite difference/finite element immersed boundary method.

Science.gov (United States)

E Griffith, Boyce; Luo, Xiaoyu

2017-12-01

The immersed boundary method is an approach to fluid-structure interaction that uses a Lagrangian description of the structural deformations, stresses, and forces along with an Eulerian description of the momentum, viscosity, and incompressibility of the fluid-structure system. The original immersed boundary methods described immersed elastic structures using systems of flexible fibers, and even now, most immersed boundary methods still require Lagrangian meshes that are finer than the Eulerian grid. This work introduces a coupling scheme for the immersed boundary method to link the Lagrangian and Eulerian variables that facilitates independent spatial discretizations for the structure and background grid. This approach uses a finite element discretization of the structure while retaining a finite difference scheme for the Eulerian variables. We apply this method to benchmark problems involving elastic, rigid, and actively contracting structures, including an idealized model of the left ventricle of the heart. Our tests include cases in which, for a fixed Eulerian grid spacing, coarser Lagrangian structural meshes yield discretization errors that are as much as several orders of magnitude smaller than errors obtained using finer structural meshes. The Lagrangian-Eulerian coupling approach developed in this work enables the effective use of these coarse structural meshes with the immersed boundary method. This work also contrasts two different weak forms of the equations, one of which is demonstrated to be more effective for the coarse structural discretizations facilitated by our coupling approach. © 2017 The Authors International  Journal  for  Numerical  Methods  in  Biomedical  Engineering Published by John Wiley & Sons Ltd.

Finite element analysis of the biaxial cyclic tensile loading of the elastoplastic plate with the central hole: asymptotic regimes

Science.gov (United States)

Turkova, Vera; Stepanova, Larisa

2018-03-01

For elastistoplastic structure elements under cyclic loading three types of asymptotic behavior are well known: shakedown, cyclic plasticity or ratcheting. In structure elements operating in real conditions ratcheting must always be excluded since it caused the incremental fracture of structure by means of the accumulation of plastic strains. In the present study results of finite-element (FEM) calculations of the asymptotical behavior of an elastoplastic plate with the central circular and elliptic holes under the biaxial cyclic loading for three different materials are presented. Incremental cyclic loading of the sample with stress concentrator (the central hole) is performed in the multifunctional finite-element package SIMULIA Abaqus. The ranges of loads found for shakedown, cyclic plasticity and ratcheting are presented. The results obtained are generalized and analyzed. Convenient normalization is suggested. The chosen normalization allows us to present all computed results, corresponding to separate materials, within one common curve with minimum scattering of the points. Convenience of the generalized diagram consists in a possibility to find an asymptotical behavior of an inelastic structure for materials for which computer calculations were not made.

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.

An inverse method based on finite element model to derive the plastic flow properties from non-standard tensile specimens of Eurofer97 steel

Directory of Open Access Journals (Sweden)

S. Knitel

2016-12-01

Full Text Available A new inverse method was developed to derive the plastic flow properties of non-standard disk tensile specimens, which were so designed to fit irradiation rods used for spallation irradiations in SINQ (Schweizer Spallations Neutronen Quelle target at Paul Scherrer Institute. The inverse method, which makes use of MATLAB and the finite element code ABAQUS, is based upon the reconstruction of the load-displacement curve by a succession of connected small linear segments. To do so, the experimental engineering stress/strain curve is divided into an elastic and a plastic section, and the plastic section is further divided into small segments. Each segment is then used to determine an associated pair of true stress/plastic strain values, representing the constitutive behavior. The main advantage of the method is that it does not rely on a hypothetic analytical expression of the constitutive behavior. To account for the stress/strain gradients that develop in the non-standard specimen, the stress and strain were weighted over the volume of the deforming elements. The method was validated with tensile tests carried out at room temperature on non-standard flat disk tensile specimens as well as on standard cylindrical specimens made of the reduced-activation tempered martensitic steel Eurofer97. While both specimen geometries presented a significant difference in terms of deformation localization during necking, the same true stress/strain curve was deduced from the inverse method. The potential and usefulness of the inverse method is outlined for irradiated materials that suffer from a large uniform elongation reduction.

A 2D finite element implementation of the Fleck–Willis strain-gradient flow theory

DEFF Research Database (Denmark)

Nielsen, Kim Lau; Niordson, Christian Frithiof

2013-01-01

The lay-out of a numerical solution procedure for the strain gradient flow (rate-independent) theory by Fleck and Willis [A mathematical basis for strain-gradient theory – Part II: Tensorial plastic multiplier, 57:1045–1057; 2009, JMPS] has been an open issue, and its finite element implementation...

Determination of Weight Suspension Rigidity in the Transport-Erector Aggregates

Directory of Open Access Journals (Sweden)

V. A. Zverev

2016-01-01

Full Text Available The aim is to determine weight suspension rigidity in aggregates designed to perform technological transport-erector operations at the miscellaneous launch complexes.We consider the weight suspension comprising the following distinctive structural components: the executive weight-lowering mechanism, polyspast mechanism, rope, traverse, and rods. A created structural dynamic model of suspension allowed us to define weight suspension rigidity. Within the framework of design analysis of a dynamic model we determined the rigidity of its structural units, i.e. traverse, rope, and polyspast.Known analytical relationships were used to calculate the rope rigidity. To determine rigidity of polyspast and traverse have been created special models based on the finite element method. For each model deformation in the specific points under the test load have been defined. Data obtained were used to determine trigidity of traverses and polyspast, and also rigidity of suspension in total. The rigidity models of polispast mechanism and traverse have been developed and calculated using the software complex "Zenit-95".As the research results, the paper presents a dynamic model of the weight suspension of the transport-erector aggregate, the finite element models of the polispast mechanism and traverse, an algorithm for determining the weight suspension rigidity and relevant analytical relationships.Independent calculation of weight suspension rigidity enables us to simplify further dynamic calculation of the aggregate-weight system because it allows attaining a simpler model of the aggregate-weight system that uses the weight suspension model as an element of equivalent rigidity. Despite this simplification the model allows us to determine correctly weight movement parameters and overloads in the aggregate-weight system in the process of technical operations.

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.

Finite element simulation of nanoindentation tests using a macroscopic computational model

International Nuclear Information System (INIS)

Khelifa, Mourad; Fierro, Vanessa; Celzard, Alain

2014-01-01

The aim of this work was to develop a numerical procedure to simulate nanoindentation tests using a macroscopic computational model. Both theoretical and numerical aspects of the proposed methodology, based on the coupling of isotropic elasticity and anisotropic plasticity described with the quadratic criterion of Hill are presented to model this behaviour. The anisotropic plastic behaviour accounts for the mixed nonlinear hardening (isotropic and kinematic) under large plastic deformation. Nanoindentation tests were simulated to analyse the nonlinear mechanical behaviour of aluminium alloy. The predicted results of the finite element (FE) modelling are in good agreement with the experimental data, thereby confirming the accuracy level of the suggested FE method of analysis. The effects of some technological and mechanical parameters known to have an influence during the nanoindentation tests were also investigated.

Finite element methods in a simulation code for offshore wind turbines

Science.gov (United States)

Kurz, Wolfgang

1994-06-01

Offshore installation of wind turbines will become important for electricity supply in future. Wind conditions above sea are more favorable than on land and appropriate locations on land are limited and restricted. The dynamic behavior of advanced wind turbines is investigated with digital simulations to reduce time and cost in development and design phase. A wind turbine can be described and simulated as a multi-body system containing rigid and flexible bodies. Simulation of the non-linear motion of such a mechanical system using a multi-body system code is much faster than using a finite element code. However, a modal representation of the deformation field has to be incorporated in the multi-body system approach. The equations of motion of flexible bodies due to deformation are generated by finite element calculations. At Delft University of Technology the simulation code DUWECS has been developed which simulates the non-linear behavior of wind turbines in time domain. The wind turbine is divided in subcomponents which are represented by modules (e.g. rotor, tower etc.).

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

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

International Nuclear Information System (INIS)

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

2007-01-01

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

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

Rigid missiles impact on reinforced concrete structures: analysis by discrete element method

International Nuclear Information System (INIS)

Shiu, W.J.

2008-10-01

The constructions likely to be subjected to some extreme loadings like reactor containment buildings have to be dimensioned accordingly. As a part of study of concrete structures, this thesis focuses on numerical modelling of rigid missile impacts against a rigid reinforced concrete slab. Based on some experiment tests data, an elasto-plastic-damaged constitutive law has been implanted into a discrete element numerical code. To calibrate certain parameters of the numerical model, some quasi static tests have been first simulated. Once the model calibration was done, some missile impact simulation tests have then been carried out. The numerical results are well agree with these provided by French Atomic Energy Agency (Cea) and the French Electrical power Company (EDF) in terms of the trajectory of the missile. We were able to show the need of a constitutive law taking into account the compaction behaviour of the concrete when the predictions of penetration and perforation of a thick slab was demanded. Finally, a parametric study confirmed that the numerical model can be used the way predictive as well as the empirical prediction law, while the first can provide additional significant mechanical description. (author)

The Numerical Simulation of the Crack Elastoplastic Extension Based on the Extended Finite Element Method

Directory of Open Access Journals (Sweden)

Xia Xiaozhou

2013-01-01

Full Text Available In the frame of the extended finite element method, the exponent disconnected function is introduced to reflect the discontinuous characteristic of crack and the crack tip enrichment function which is made of triangular basis function, and the linear polar radius function is adopted to describe the displacement field distribution of elastoplastic crack tip. Where, the linear polar radius function form is chosen to decrease the singularity characteristic induced by the plastic yield zone of crack tip, and the triangle basis function form is adopted to describe the displacement distribution character with the polar angle of crack tip. Based on the displacement model containing the above enrichment displacement function, the increment iterative form of elastoplastic extended finite element method is deduced by virtual work principle. For nonuniform hardening material such as concrete, in order to avoid the nonsymmetry characteristic of stiffness matrix induced by the non-associate flowing of plastic strain, the plastic flowing rule containing cross item based on the least energy dissipation principle is adopted. Finally, some numerical examples show that the elastoplastic X-FEM constructed in this paper is of validity.

Stresses in Circular Plates with Rigid Elements

Science.gov (United States)

Velikanov, N. L.; Koryagin, S. I.; Sharkov, O. V.

2018-05-01

Calculations of residual stress fields are carried out by numerical and static methods, using the flat cross-section hypothesis. The failure of metal when exposed to residual stresses is, in most cases, brittle. The presence in the engineering structures of rigid elements often leads to the crack initiation and structure failure. This is due to the fact that rigid elements under the influence of external stresses are stress concentrators. In addition, if these elements are fixed by welding, the residual welding stresses can lead to an increase in stress concentration and, ultimately, to failure. The development of design schemes for such structures is a very urgent task for complex technical systems. To determine the stresses in a circular plate with a welded circular rigid insert under the influence of an external load, one can use the solution of the plane stress problem for annular plates in polar coordinates. The polar coordinates of the points are the polar radius and the polar angle, and the stress state is determined by normal radial stresses, tangential and shearing stresses. The use of the above mentioned design schemes, formulas, will allow more accurate determination of residual stresses in annular welded structures. This will help to establish the most likely directions of failure and take measures at the stages of designing, manufacturing and repairing engineering structures to prevent these failures. However, it must be taken into account that the external load, the presence of insulation can lead to a change in the residual stress field.

Thermo-plastic finite element analysis for metal honeycomb structure

Directory of Open Access Journals (Sweden)

Ji Zhanling

2013-01-01

Full Text Available This paper deals with thermal-plastic analysis for the metal honeycomb structure. The heat transfer equation and thermal elastoplastic constitutive equation of a multilayer panel are established and studied numerically using ANSYS software. The paper elucidates that only the outer skin produces easily plastic deformation, and the outer skin still exists some residual stress and residual deformation after cooling. The dynamic evolution of plastic deformation and material performance degradation under high energy thermal load are revealed.

One-Dimensional Finite Elements An Introduction to the FE Method

CERN Document Server

Öchsner, Andreas

2013-01-01

 This textbook presents finite element methods using exclusively  one-dimensional elements. The aim is to present the complex methodology in  an easily understandable but mathematically correct fashion. The approach of  one-dimensional elements enables the reader to focus on the understanding of  the principles of basic and advanced mechanical problems. The reader easily  understands the assumptions and limitations of mechanical modeling as well  as the underlying physics without struggling with complex mathematics. But  although the description is easy it remains scientifically correct.   The approach using only one-dimensional elements covers not only standard  problems but allows also for advanced topics like plasticity or the  mechanics of composite materials. Many examples illustrate the concepts and  problems at the end of every chapter help to familiarize with the topics.

Evaluation and optimization of footwear comfort parameters using finite element analysis and a discrete optimization algorithm

Science.gov (United States)

Papagiannis, P.; Azariadis, P.; Papanikos, P.

2017-10-01

Footwear is subject to bending and torsion deformations that affect comfort perception. Following review of Finite Element Analysis studies of sole rigidity and comfort, a three-dimensional, linear multi-material finite element sole model for quasi-static bending and torsion simulation, overcoming boundary and optimisation limitations, is described. Common footwear materials properties and boundary conditions from gait biomechanics are used. The use of normalised strain energy for product benchmarking is demonstrated along with comfort level determination through strain energy density stratification. Sensitivity of strain energy against material thickness is greater for bending than for torsion, with results of both deformations showing positive correlation. Optimization for a targeted performance level and given layer thickness is demonstrated with bending simulations sufficing for overall comfort assessment. An algorithm for comfort optimization w.r.t. bending is presented, based on a discrete approach with thickness values set in line with practical manufacturing accuracy. This work illustrates the potential of the developed finite element analysis applications to offer viable and proven aids to modern footwear sole design assessment and optimization.

Prediction of plastic deformation under contact condition by quasi-static and dynamic simulations using explicit finite element analysis

International Nuclear Information System (INIS)

Siswanto, W. A.; Nagentrau, M.; Tobi, A. L. Mohd; Tamin, M. N.

2016-01-01

We compared the quasi-static and dynamic simulation responses on elastic-plastic deformation of advanced alloys using Finite element (FE) method with an explicit numerical algorithm. A geometrical model consisting of a cylinder-on-flat surface contact under a normal load and sliding motion was examined. Two aeroengine materials, Ti-6Al-4V and Super CMV (Cr-Mo-V) alloy, were employed in the FE analysis. The FE model was validated by comparative magnitudes of the FE-predicted maximum contact pressure variation along the contact half-width length with the theoretical Hertzian contact solution. Results show that the (compressive) displacement of the initial contact surface steadily increases for the quasi-static load case, but accumulates at an increasing rate to the maximum level for the dynamic loading. However, the relatively higher stiffness and yield strength of the Super CMV alloy resulted in limited deformation and low plastic strain when compared to the Ti-6Al-4V alloy. The accumulated equivalent plastic strain of the material point at the initial contact position was nearly a thousand times higher for the dynamic load case (for example, 6.592 for Ti-6Al-4V, 1.0 kN) when compared to the quasi-static loading (only 0.0072). During the loading step, the von Mises stress increased with a decreasing and increasing rate for the quasi-static and dynamic load case, respectively. A sudden increase in the stress magnitude to the respective peak value was registered due to the additional constraint to overcome the static friction of the mating surfaces during the sliding step

Prediction of plastic deformation under contact condition by quasi-static and dynamic simulations using explicit finite element analysis

Energy Technology Data Exchange (ETDEWEB)

Siswanto, W. A.; Nagentrau, M.; Tobi, A. L. Mohd [Faculty of Mechanical and Manufacturing Engineering, Universiti Tun Hussein Onn Malaysia, Batu Pahat (Malaysia); Tamin, M. N. [Faculty of Mechanical Engineering, Universiti Teknologi Malaysia, Johor Bahru (Malaysia)

2016-11-15

We compared the quasi-static and dynamic simulation responses on elastic-plastic deformation of advanced alloys using Finite element (FE) method with an explicit numerical algorithm. A geometrical model consisting of a cylinder-on-flat surface contact under a normal load and sliding motion was examined. Two aeroengine materials, Ti-6Al-4V and Super CMV (Cr-Mo-V) alloy, were employed in the FE analysis. The FE model was validated by comparative magnitudes of the FE-predicted maximum contact pressure variation along the contact half-width length with the theoretical Hertzian contact solution. Results show that the (compressive) displacement of the initial contact surface steadily increases for the quasi-static load case, but accumulates at an increasing rate to the maximum level for the dynamic loading. However, the relatively higher stiffness and yield strength of the Super CMV alloy resulted in limited deformation and low plastic strain when compared to the Ti-6Al-4V alloy. The accumulated equivalent plastic strain of the material point at the initial contact position was nearly a thousand times higher for the dynamic load case (for example, 6.592 for Ti-6Al-4V, 1.0 kN) when compared to the quasi-static loading (only 0.0072). During the loading step, the von Mises stress increased with a decreasing and increasing rate for the quasi-static and dynamic load case, respectively. A sudden increase in the stress magnitude to the respective peak value was registered due to the additional constraint to overcome the static friction of the mating surfaces during the sliding step.

Stress categorization in nozzle to pressure vessel connections finite elements models

International Nuclear Information System (INIS)

Albuquerque, Levi Barcelos de

1999-01-01

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

Finite Element Based Design Optimization of WENDELSTEIN 7-X Divertor Targets

International Nuclear Information System (INIS)

Plankensteiner, A.; Leuprecht, A.; Schedler, B.; Scheiber, K.; Greuner, H.

2006-01-01

In the fusion experiment WENDELSTEIN 7-X divertor plasma facing components have to withstand severe loading conditions. In general thermally induced mechanical stressing turns out to be most critical with respect to life time predictions of the component. In the specific case flat tiles of CFC grade NB31 are joined to the precipitation hardened CuCrZr heat sink by employing an active metal cast (AMC)-Cu as an interlayer between CFC and CuCrZr. Residual stresses resulting from the manufacturing process act as initial stresses in the subsequent operational heat flux loading. For the latter loading regime these stresses intrinsically are generated due to the large contrast in the CTE for CFC and Cu. Different design variants of those CFC flat tile armoured target elements have been analysed via the finite element package ABAQUS aiming at derivation of an optimized component design. The numerical study comprises variants with different degrees of tessellation of the CFC flat tile section, orientation of the CFC, lamellar design of the AMC-interlayer, and different designs of the cooling channels. The thermo-mechanical material characteristics are accounted for the finite element models with elastic-plastic properties being assigned to the metallic sections CuCrZr and AMC-Cu, respectively, and orthotropic nonlinear-elastic properties being used to the CFC section. The latter has been realized in form of a user-defined material subroutine that is used at the integration point level of the finite element model. In particular, twelve scalar-type damage parameters obeying their own evolution equations with respect to the loading history account for specific stress-strain relationships in the three principal material directions and planes with six damage parameters being used for normal loading under tensile and compressive stress states, respectively, and six parameters being used for shear loading. For the aim of model verification calculated surface temperatures, global

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

International Nuclear Information System (INIS)

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

2004-01-01

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

Analysis of Switched-Rigid Floating Oscillator

Directory of Open Access Journals (Sweden)

Prabhakar R. Marur

2009-01-01

Full Text Available In explicit finite element simulations, a technique called deformable-to-rigid (D2R switching is used routinely to reduce the computation time. Using the D2R option, the deformable parts in the model can be switched to rigid and reverted back to deformable when needed during the analysis. The time of activation of D2R however influences the overall dynamics of the system being analyzed. In this paper, a theoretical basis for the selection of time of rigid switching based on system energy is established. A floating oscillator problem is investigated for this purpose and closed-form analytical expressions are derived for different phases in rigid switching. The analytical expressions are validated by comparing the theoretical results with numerical computations.

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

Energy Technology Data Exchange (ETDEWEB)

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

2016-07-01

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

Strain-based finite elements for the analysis of cylinders with holes and normally intersecting cylinders

International Nuclear Information System (INIS)

Sabir, A.B.

1983-01-01

A finite element solution to the problems of stress distribution for cylindrical shells with circular and elliptical holes and also for normally intersecting thin elastic cylindrical shells is given. Quadrilateral and triangular curved finite elements are used in the analysis. The elements are of a new class, based on simple independent generalised strain functions insofar as this is allowed by the compatibility equations. The elements also satisfy exactly the requirements of strain-free-rigid body displacements and uses only the external 'geometrical' nodal degrees of freedom to avoid the difficulties associated with unnecessary internal degrees of freedom. We first develop strain based quadrilateral and triangular elements and apply them to the solution of the problem of stress concentrations in the neighbourhood of small and large circular and elliptical holes when the cylinders are subjected to a uniform axial tension. These results are compared with analytical solutions based on shallow shell approximations and show that the use of these strain based elements obviates the need for using an inordinately large number of elements. Normally intersecting cylinders are common configurations in structural components for nuclear reactor systems and design information for such configurations are generally lacking. The opportunity is taken in the present paper to provide a finite element solution to this problem. A method of substructing will be introduced to enable a solution to the large number of non banded set of simultaneous equations encountered. (orig./HP)

Determination of the ultimate load in concrete slabs by the yield line finite element method

International Nuclear Information System (INIS)

Vaz, L.E.; Feijo, B.; Martha, L.F.R.; Lopes, M.M.

1984-01-01

A method for calculating the ultimate load in reinforced concrete slabs is proposed. The method follows the finite element aproach representating the continuum slab as an assembly of rigid triangular plates connected along their sides through yield line elements. This approach leads to the definition of the displacement configuration of the plate only as a function of the transversal displacement at the nodes of the mesh (1 DOF per node) reducing significantly the number of DOF's in relation to the conventional formulation by means of the finite element method (minimum of 3 DOF per node). Nonlinear behaviour of the reinforced concrete section is considered in the definition of the moment rotation curve of the yield lines. The effect of the in plane forces acting in the middle surface of the plate is also taken into account. The validity of the model is verified comparing the numerical solutions with the results of the classical yield line theory. (Author) [pt

Active earth pressure model tests versus finite element analysis

Science.gov (United States)

Pietrzak, Magdalena

2017-06-01

The purpose of the paper is to compare failure mechanisms observed in small scale model tests on granular sample in active state, and simulated by finite element method (FEM) using Plaxis 2D software. Small scale model tests were performed on rectangular granular sample retained by a rigid wall. Deformation of the sample resulted from simple wall translation in the direction `from the soil" (active earth pressure state. Simple Coulomb-Mohr model for soil can be helpful in interpreting experimental findings in case of granular materials. It was found that the general alignment of strain localization pattern (failure mechanism) may belong to macro scale features and be dominated by a test boundary conditions rather than the nature of the granular sample.

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.

Theoretical basis for a transient thermal elastic-plastic stress analysis of nuclear reactor fuel elements

International Nuclear Information System (INIS)

Hsu, T.R.; Bertels, A.W.M.; Banerjee, S.; Harrison, W.C.

1976-07-01

This report presents the theoretical basis for a transient thermal elastic-plastic stress analysis of a nuclear reactor fuel element subject to severe transient thermo-mechanical loading. A finite element formulation is used for both the non-linear stress analysis and thermal analysis. These two major components are linked together to form an integrated program capable of predicting fuel element transient behaviour in two dimensions. Specific case studies are presented to illustrate capabilities of the analysis. (author)

Friction welding; Magnesium; Finite element; Shear test.

Directory of Open Access Journals (Sweden)

Leonardo Contri Campanelli

2013-06-01

Full Text Available Friction spot welding (FSpW is one of the most recently developed solid state joining technologies. In this work, based on former publications, a computer aided draft and engineering resource is used to model a FSpW joint on AZ31 magnesium alloy sheets and subsequently submit the assembly to a typical shear test loading, using a linear elastic model, in order to conceive mechanical tests results. Finite element analysis shows that the plastic flow is concentrated on the welded zone periphery where yield strength is reached. It is supposed that “through the weld” and “circumferential pull-out” variants should be the main failure behaviors, although mechanical testing may provide other types of fracture due to metallurgical features.

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.

Nonlinear finite element analyses: advances and challenges in dental applications.

Science.gov (United States)

Wakabayashi, N; Ona, M; Suzuki, T; Igarashi, Y

2008-07-01

To discuss the development and current status of application of nonlinear finite element method (FEM) in dentistry. The literature was searched for original research articles with keywords such as nonlinear, finite element analysis, and tooth/dental/implant. References were selected manually or searched from the PUBMED and MEDLINE databases through November 2007. The nonlinear problems analyzed in FEM studies were reviewed and categorized into: (A) nonlinear simulations of the periodontal ligament (PDL), (B) plastic and viscoelastic behaviors of dental materials, (C) contact phenomena in tooth-to-tooth contact, (D) contact phenomena within prosthodontic structures, and (E) interfacial mechanics between the tooth and the restoration. The FEM in dentistry recently focused on simulation of realistic intra-oral conditions such as the nonlinear stress-strain relationship in the periodontal tissues and the contact phenomena in teeth, which could hardly be solved by the linear static model. The definition of contact area critically affects the reliability of the contact analyses, especially for implant-abutment complexes. To predict the failure risk of a bonded tooth-restoration interface, it is essential to assess the normal and shear stresses relative to the interface. The inclusion of viscoelasticity and plastic deformation to the program to account for the time-dependent, thermal sensitive, and largely deformable nature of dental materials would enhance its application. Further improvement of the nonlinear FEM solutions should be encouraged to widen the range of applications in dental and oral health science.

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)

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.

Application of a general purpose finite element program system in pressure vessel technology

International Nuclear Information System (INIS)

Aamodt, B.; Sandsmark, N.; Medonos, S.

1977-01-01

Main advantages of using general purpose finite element program systems in structural analysis are summarized. Several illustrative applications of the program system SESAM-69 to pressure vessel problems are described. The first example is a dynamic analysis of the motor housing of the internal main circulation pump of a BWR nuclear reactor. The next example is a transient heat conduction and stress analysis of deflector of feeding nozzle of PWR nuclear reactor. Then, numerical calculations of stress intensity factors and fatigue crack growth of semi-elliptical surface cracks are discussed. And finally, an elasto-plastic analysis of a thick plate with edge-cracks is considered. It is concluded that due to the fact that general purpose finite element program systems are general and user-orientated, they will gain increasingly higher popularity in the years ahead

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.

Finite element discretization of Darcy's equations with pressure dependent porosity

KAUST Repository

Girault, Vivette

2010-02-23

We consider the flow of a viscous incompressible fluid through a rigid homogeneous porous medium. The permeability of the medium depends on the pressure, so that the model is nonlinear. We propose a finite element discretization of this problem and, in the case where the dependence on the pressure is bounded from above and below, we prove its convergence to the solution and propose an algorithm to solve the discrete system. In the case where the dependence on the pressure is exponential, we propose a splitting scheme which involves solving two linear systems, but parts of the analysis of this method are still heuristic. Numerical tests are presented, which illustrate the introduced methods. © 2010 EDP Sciences, SMAI.

Bending behaviors of fully covered biodegradable polydioxanone biliary stent for human body by finite element method.

Science.gov (United States)

Liu, Yanhui; Zhu, Guoqing; Yang, Huazhe; Wang, Conger; Zhang, Peihua; Han, Guangting

2018-01-01

This paper presents a study of the bending flexibility of fully covered biodegradable polydioxanone biliary stents (FCBPBs) developed for human body. To investigate the relationship between the bending load and structure parameter (monofilament diameter and braid-pin number), biodegradable polydioxanone biliary stents derived from braiding method were covered with membrane prepared via electrospinning method, and nine FCBPBSs were then obtained for bending test to evaluate the bending flexibility. In addition, by the finite element method, nine numerical models based on actual biliary stent were established and the bending load was calculated through the finite element method. Results demonstrate that the simulation and experimental results are in good agreement with each other, indicating that the simulation results can be provided a useful reference to the investigation of biliary stents. Furthermore, the stress distribution on FCBPBSs was studied, and the plastic dissipation analysis and plastic strain of FCBPBSs were obtained via the bending simulation. Copyright © 2017 Elsevier Ltd. All rights reserved.

Patient-specific non-linear finite element modelling for predicting soft organ deformation in real-time: application to non-rigid neuroimage registration.

Science.gov (United States)

Wittek, Adam; Joldes, Grand; Couton, Mathieu; Warfield, Simon K; Miller, Karol

2010-12-01

Long computation times of non-linear (i.e. accounting for geometric and material non-linearity) biomechanical models have been regarded as one of the key factors preventing application of such models in predicting organ deformation for image-guided surgery. This contribution presents real-time patient-specific computation of the deformation field within the brain for six cases of brain shift induced by craniotomy (i.e. surgical opening of the skull) using specialised non-linear finite element procedures implemented on a graphics processing unit (GPU). In contrast to commercial finite element codes that rely on an updated Lagrangian formulation and implicit integration in time domain for steady state solutions, our procedures utilise the total Lagrangian formulation with explicit time stepping and dynamic relaxation. We used patient-specific finite element meshes consisting of hexahedral and non-locking tetrahedral elements, together with realistic material properties for the brain tissue and appropriate contact conditions at the boundaries. The loading was defined by prescribing deformations on the brain surface under the craniotomy. Application of the computed deformation fields to register (i.e. align) the preoperative and intraoperative images indicated that the models very accurately predict the intraoperative deformations within the brain. For each case, computing the brain deformation field took less than 4 s using an NVIDIA Tesla C870 GPU, which is two orders of magnitude reduction in computation time in comparison to our previous study in which the brain deformation was predicted using a commercial finite element solver executed on a personal computer. Copyright © 2010 Elsevier Ltd. All rights reserved.

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.

Concrete model for finite element analysis of structures subjected to severe damages

International Nuclear Information System (INIS)

Jamet, Ph.; Millard, A.; Hoffmann, A.; Nahas, G.; Barbe, B.

1984-01-01

A specific concrete model has been developed, in order to perform mechanical analysis of civil engineering structures, when subjected to accidental loadings, leading to severe damages. Its formulation is based on the physical mechanisms, which have been observed on laboratory specimens. The model has been implemented into the CASTEM finite element system, and the case of a concrete slab perforation by a rigid missile has been considered. The qualitative behaviour of the structure is well predicted by the model. Comparison between numerical and experimental results is also performed, using two main curves: missile velocity versus penetration depth; reaction forces versus time. (Author) [pt

A new approach for elasto-plastic finite strain analysis of cantilever ...

Indian Academy of Sciences (India)

A new approach for elasto-plastic finite strain analysis of cantilever beams subjected to uniform bending moment ... Curvature; deflection curve; cantilever beam; elasto-plastic analysis; tapered beam subjected to tipmoment; ... Sadhana | News.

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.

Investigation of High-Speed Cryogenic Machining Based on Finite Element Approach

Directory of Open Access Journals (Sweden)

Pooyan Vahidi Pashaki

Full Text Available Abstract The simulation of cryogenic machining process because of using a three-dimensional model and high process duration time in the finite element method, have been studied rarely. In this study, to overcome this limitation, a 2.5D finite element model using the commercial finite element software ABAQUS has been developed for the cryogenic machining process and by considering more realistic assumptions, the chip formation procedure investigated. In the proposed method, the liquid nitrogen has been used as a coolant. At the modeling of friction during the interaction of tools - chip, the Coulomb law has been used. In order to simulate the behavior of plasticity and failure criterion, Johnson-Cook model was used, and unlike previous investigations, thermal and mechanical properties of materials as a function of temperature were applied to the software. After examining accuracy of the model with present experimental data, the effect of parameters such as rake angle and the cutting speed as well as dry machining of aluminum alloy by the use of coupled dynamic temperature solution has been studied. Results indicated that at the cutting velocity of 10 m/s, cryogenic cooling has caused into decreasing 60 percent of tools temperature in comparison with the dry cooling. Furthermore, a chip which has been made by cryogenic machining were connected and without fracture in contrast to dry machining.

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.

Simulation of micromechanical behavior of polycrystals: finite elements versus fast Fourier transforms

International Nuclear Information System (INIS)

Prakash, A; Lebensohn, R A

2009-01-01

In this work, we compare finite element and fast Fourier transform approaches for the prediction of the micromechanical behavior of polycrystals. Both approaches are full-field approaches and use the same visco-plastic single crystal constitutive law. We investigate the texture and the heterogeneity of the inter- and intragranular stress and strain fields obtained from the two models. Additionally, we also look into their computational performance. Two cases—rolling of aluminum and wire drawing of tungsten—are used to evaluate the predictions of the two models. Results from both the models are similar, when large grain distortions do not occur in the polycrystal. The finite element simulations were found to be highly computationally intensive, in comparison with the fast Fourier transform simulations. Figure 9 was corrected in this article on the 25 August 2009. The corrected electronic version is identical to the print version

The Guidelines for Modelling the Preloading Bolts in the Structural Connection Using Finite Element Methods

OpenAIRE

Paulina Krolo; Davor Grandić; Mladen Bulić

2016-01-01

The aim of this paper is the development of the two different numerical techniques for the preloading of bolts by the finite element method using the software Abaqus Standard. Furthermore, this paper gave detailed guidelines for modelling contact, method for solving the numerical error problems such as numerical singularity error and negative eigenvalues due to rigid body motion or the problem of the extensive elongation of bolts after pretension which is occurring during the analysis. The be...

Finite-difference analysis of shells impacting rigid barriers

International Nuclear Information System (INIS)

Pirotin, S.D.; Witmer, E.A.

1977-01-01

Nuclear power plants must be protected from the adverse effects of missile impacts. A significant category of missile impact involves deformable structures (pressure vessel components, whipping pipes) striking relatively rigid targets (concrete walls, bumpers) which act as protective devices. The response and interaction of these structures is needed to assess the adequacy of these barriers for protecting vital safety related equipment. The present investigation represents an initial attempt to develop an efficient numerical procedure for predicting the deformations and impact force time-histories of shells which impact upon a rigid target. The general large-deflection equations of motion of the shell are expressed in finite-difference form in space and integrated in time through application of the central-difference temporal operator. The effect of material nonlinearities is treated by a mechanical sublayer material model which handles the strain-hardening, Bauschinger, and strain-rate effects. The general adequacy of this shell treatment has been validated by comparing predictions with the results of various experiments in which structures have been subjected to well-defined transient forcing functions (typically high-explosive impulse loading). The 'new' ingredient addressed in the present study involves an accounting for impact interaction and response of both the target structure and the attacking body. (Auth.)

An efficient finite element solution for gear dynamics

International Nuclear Information System (INIS)

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

2010-01-01

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

Finite element analysis theory and application with ANSYS

CERN Document Server

Moaveni, Saeed

2015-01-01

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

Finite element application to global reactor analysis

International Nuclear Information System (INIS)

Schmidt, F.A.R.

1981-01-01

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

Finite-Element Software for Conceptual Design

DEFF Research Database (Denmark)

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

2010-01-01

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

Books and monographs on finite element technology

Science.gov (United States)

Noor, A. K.

1985-01-01

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

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.

Plasticity and beyond microstructures, crystal-plasticity and phase transitions

CERN Document Server

Hackl, Klaus

2014-01-01

The book presents the latest findings in experimental plasticity, crystal plasticity, phase transitions, advanced mathematical modeling of finite plasticity and multi-scale modeling. The associated algorithmic treatment is mainly based on finite element formulations for standard (local approach) as well as for non-standard (non-local approach) continua and for pure macroscopic as well as for directly coupled two-scale boundary value problems. Applications in the area of material design/processing are covered, ranging from grain boundary effects in polycrystals and phase transitions to deep-drawing of multiphase steels by directly taking into account random microstructures.

Finite element analysis for prediction of the residual stresses induced by shot peening II

International Nuclear Information System (INIS)

Kim, Cheol; Seok, Chang Sung; Yang, Won Ho; Ryu, Myung Hai

2002-01-01

Shot peening is a surface impact treatment widely used to improve the performance of metal parts and welded details subjected to fatigue loading, contact fatigue, stress corrosion and other damage mechanisms. The better performance of the peened parts is mainly due to the residual stresses resulting from the plastic deformation of the surface layers of the material caused by the impact of the shot. In this paper the simulation technique is applied to predict the magnitude and distribution of the residual stress and plastic deformation caused by shot peening with the help of finite element analysis

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.

Extensions to a nonlinear finite element axisymmetric shell model based on Reissner's shell theory

International Nuclear Information System (INIS)

Cook, W.A.

1981-01-01

A finite element shell-of-revolution model has been developed to analyze shipping containers under severe impact conditions. To establish the limits for this shell model, I studied the basic assumptions used in its development; these are listed in this paper. Several extensions were evident from the study of these limits: a thick shell, a plastic hinge, and a linear normal stress. (orig./HP)

Finite element based design optimization of WENDELSTEIN 7-X divertor components under high heat flux loading

International Nuclear Information System (INIS)

Plankensteiner, A.; Leuprecht, A.; Schedler, B.; Scheiber, K.-H.; Greuner, H.

2007-01-01

In the divertor of the nuclear fusion experiment WENDELSTEIN 7-X (W7-X) plasma facing high heat flux target elements have to withstand severe loading conditions. The thermally induced mechanical stressing turns out to be most critical with respect to lifetime predictions of the target elements. Therefore, different design variants of those CFC flat tile armoured high heat flux components have been analysed via the finite element package ABAQUS aiming at derivation of an optimized component design under high heat flux conditions. The investigated design variants comprise also promising alterations in the cooling channel design and castellation of the CFC flat tiles which, however, from a system integration and manufacturing standpoint of view, respectively, are evaluated to be critical. Therefore, the numerical study as presented here mainly comprises a reference variant that is comparatively studied with a variant incorporating a bi-layer-type AMC-Cu/OF-Cu interlayer at the CFC/Cu-interface. The thermo-mechanical material characteristics are accounted for in the finite element models with elastic-plastic properties being assigned to the metallic sections CuCrZr, AMC-Cu and OF-Cu, respectively, and orthotropic nonlinear-elastic properties being used for the CFC sections. The calculated temporal and spatial evolution of temperatures, stresses, and strains for the individual design variants are evaluated with special attention being paid to stress measures, plastic strains, and damage parameters indicating the risk of failure of CFC and the CFC/Cu-interface, respectively. This way the finite element analysis allows to numerically derive an optimized design variant within the framework of expected operating conditions in W7-X

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.

Study of 3-D stress development in parent and twin pairs of a hexagonal close-packed polycrystal: Part II - Crystal plasticity finite element modeling

DEFF Research Database (Denmark)

Abdolvand, Hamidreza; Majkut, Marta; Oddershede, Jette

2015-01-01

for each grain from the 3DXRD experiment is within the stress variation zone of the grain modeled in the CPFE simulation. Also, the CPFE average stress calculation for each grain is in good agreement with the measured average stress values. It is shown that upon considering the stress variations within......Stress heterogeneity within each individual grain of polycrystalline Zircaloy-2 is studied using a crystal plasticity finite element (CPFE) model. For this purpose, the weighted Voronoi tessellation method is used to construct 3D geometries of more than 2600 grains based on their center......-of-mass positions and volumes as measured by three-dimensional X-ray diffraction (3DXRD) microscopy. The constructed microstructure is meshed with different element densities and for different numbers of grains. Then a selected group of twin and parent pairs are studied. It is shown that the measured average stress...

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

International Nuclear Information System (INIS)

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

1980-01-01

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

Calculation of foundation response to spatially varying ground motion by finite element method

International Nuclear Information System (INIS)

Wang, F.; Gantenbein, F.

1995-01-01

This paper presents a general method to compute the response of a rigid foundation of arbitrary shape resting on a homogeneous or multilayered elastic soil when subjected to a spatially varying ground motion. The foundation response is calculated from the free-field ground motion and the contact tractions between the foundation and the soil. The spatial variation of ground motion in this study is introduced by a coherence function and the contact tractions are obtained numerically using the Finite Element Method in the process of calculating the dynamic compliance of the foundation. Applications of this method to a massless rigid disc supported on an elastic half space and to that founded on an elastic medium consisting of a layer of constant thickness supported on an elastic half space are described. The numerical results obtained are in very good agreement with analytical solutions published in the literature. (authors). 5 refs., 8 figs

Finite element simulations and experiments to determine the residual damage of a CFRP composite material after ballistic impacts

NARCIS (Netherlands)

Herlaar, K.; Jagt-Deutekom, M. van der

2005-01-01

The use of lightweight systems is essential for future combat systems. More and more steel structures are replaced by composite structures. This also influences the vulnerability of the platform. A finite element material model is created in Autodyn of the carbon fiber reinforced plastic (CFRP)

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

Two-Dimensional Nonlinear Finite Element Analysis of CMC Microstructures

Science.gov (United States)

Mital, Subodh K.; Goldberg, Robert K.; Bonacuse, Peter J.

2012-01-01

A research program has been developed to quantify the effects of the microstructure of a woven ceramic matrix composite and its variability on the effective properties and response of the material. In order to characterize and quantify the variations in the microstructure of a five harness satin weave, chemical vapor infiltrated (CVI) SiC/SiC composite material, specimens were serially sectioned and polished to capture images that detailed the fiber tows, matrix, and porosity. Open source quantitative image analysis tools were then used to isolate the constituents, from which two dimensional finite element models were generated which approximated the actual specimen section geometry. A simplified elastic-plastic model, wherein all stress above yield is redistributed to lower stress regions, is used to approximate the progressive damage behavior for each of the composite constituents. Finite element analyses under in-plane tensile loading were performed to examine how the variability in the local microstructure affected the macroscopic stress-strain response of the material as well as the local initiation and progression of damage. The macroscopic stress-strain response appeared to be minimally affected by the variation in local microstructure, but the locations where damage initiated and propagated appeared to be linked to specific aspects of the local microstructure.

The problem solution on wedge penetration in an initially anisotropic medium within the rigid-plastic scheme

Science.gov (United States)

Chanyshev, AI; Abdulin, IM

2018-03-01

Two problems are solved in the paper: on ultimate loads in the initial stage of indentation of an absolutely rigid smooth wedge into a layer of an initially anisotropic plastic medium and in the final stage when the tool penetrates through the layer. The problems are solved with Chanyshev’s constitutive relations of plasticity of the initially anisotropic medium based on use of the eigen elasticity tensors.

Analysis of the Behaviour of Semi Rigid Steel End Plate Connections

Directory of Open Access Journals (Sweden)

Bahaz A.

2018-01-01

Full Text Available The analysis of steel-framed building structures with full strength beam to column joints is quite standard nowadays. Buildings utilizing such framing systems are widely used in design practice. However, there is a growing recognition of significant benefits in designing joints as partial strength/semi-rigid. The design of joints within this partial strength/semi-rigid approach is becoming more and more popular. This requires the knowledge of the full nonlinear moment-rotation behaviour of the joint, which is also a design parameter. The rotational behaviour of steel semi rigid connections can be studied using the finite element method for the following three reasons: i such models are inexpensive; ii they allow the understanding of local effects, which are difficult to measure accurately physically, and iii they can be used to generate extensive parametric studies. This paper presents a three-dimensional finite element model using ABAQUS software in order to identify the effect of different parameters on the behaviour of semi rigid steel beam to column end plate connections. Contact and sliding between different elements, bolt pretension and geometric and material non-linearity are included in this model. A parametric study is conducted using a model of two end-plate configurations: flush and extended end plates. The studied parameters were as follows: bolts type, end plate thickness and column web stiffener. Then, the model was calibrated and validated with experimental results taken from the literature and with the model proposed by Eurocode3. The procedure for determining the moment–rotation curve using finite element analysis is also given together with a brief explanation of how the design moment resistance and the initial rotational stiffness of the joint are obtained.

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.

Cyclic behavior of 316L steel predicted by means of finite element computations

International Nuclear Information System (INIS)

Liu, J.; Sauzay, M.; Robertson, C.; Liu, J.

2011-01-01

The cyclic behavior of 316L steels is predicted based on crystalline elastoplastic constitutive laws. Calculations are performed with the finite element software CAST3M, using a polycrystalline mesh where the individual grains are modeled as cubes, having random crystallographic orientations. At the grain scale, the constitutive law parameters are adjusted using single crystal cyclic stress strain curves (CSSCs) from literature. Calculations are performed for different loading conditions (uniaxial tension-compression, biaxial tension-compression and alternated torsion) and a large range of three remote plastic strain amplitudes. We obtained 3 close macroscopic CSSCs. Somewhat lower stresses are obtained in torsion, particularly at high plastic strain amplitude. Our results are in agreement with all the published experimental data. The mean plastic strain is computed in each grain, yielding a particular polycrystalline mean grain plastic strain distribution for each loading condition and remote plastic strain. The plastic strain scatter increases for decreasing macroscopic strains. The number of cycles to the first micro-crack initiation corresponding to the aforesaid plastic strain distributions is then calculated using a surface roughness based initiation criterion. The effect of the different loading conditions is finally discussed. (authors)

Surgery simulation using fast finite elements

DEFF Research Database (Denmark)

Bro-Nielsen, Morten

1996-01-01

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

Using Finite Element Method

Directory of Open Access Journals (Sweden)

M.H.R. Ghoreishy

2008-02-01

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

Precise magnetostatic field using the finite element method

International Nuclear Information System (INIS)

Nascimento, Francisco Rogerio Teixeira do

2013-01-01

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

Finite element elasto-plastic analysis of thin walled structures of reinforced concrete as applied to reactor facilities

International Nuclear Information System (INIS)

Fujita, F.; Tsuboi, Y.

1981-01-01

The authors developed a new program of elasto-plastic analysis of reinforced concrete shells, in which the simplest model of shell element and an orthotropic constitutive relation are adopted, and verified its validity with reference to the results of model experiments of containers and box-wall structures with various loading conditions. For the two-dimensional stress-strain relationship of concrete, an orthotropic nonlinear formula proposed by one of the authors was adopted. For concrete, the octahedral shear failure and tension cut-off criteria were also imposed. The Kirchhoff-Love's assumptions were assumed to be valid for the whole range of the analysis and the layered approach of elasto-plastic stiffness evaluation. Derivation of the shell element is outlined with examination of its accuracy in elastic range and the assumption of elasto-plastic material property and the procedure of nonlinear analysis are described. As examples, the method is applied to the analysis of a cylindrical container and a box-wall structure. Comparison of the computed results with the corresponding experimental data indicates the applicability of the proposed method. (orig./HP)

Vector form Intrinsic Finite Element Method for the Two-Dimensional Analysis of Marine Risers with Large Deformations

Science.gov (United States)

Li, Xiaomin; Guo, Xueli; Guo, Haiyan

2018-06-01

Robust numerical models that describe the complex behaviors of risers are needed because these constitute dynamically sensitive systems. This paper presents a simple and efficient algorithm for the nonlinear static and dynamic analyses of marine risers. The proposed approach uses the vector form intrinsic finite element (VFIFE) method, which is based on vector mechanics theory and numerical calculation. In this method, the risers are described by a set of particles directly governed by Newton's second law and are connected by weightless elements that can only resist internal forces. The method does not require the integration of the stiffness matrix, nor does it need iterations to solve the governing equations. Due to these advantages, the method can easily increase or decrease the element and change the boundary conditions, thus representing an innovative concept of solving nonlinear behaviors, such as large deformation and large displacement. To prove the feasibility of the VFIFE method in the analysis of the risers, rigid and flexible risers belonging to two different categories of marine risers, which usually have differences in modeling and solving methods, are employed in the present study. In the analysis, the plane beam element is adopted in the simulation of interaction forces between the particles and the axial force, shear force, and bending moment are also considered. The results are compared with the conventional finite element method (FEM) and those reported in the related literature. The findings revealed that both the rigid and flexible risers could be modeled in a similar unified analysis model and that the VFIFE method is feasible for solving problems related to the complex behaviors of marine risers.

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

Global and Local Mechanical Responses for Necking of Rectangular Bars Using Updated and Total Lagrangian Finite Element Formulations

Directory of Open Access Journals (Sweden)

Claudio A. Careglio

2016-01-01

Full Text Available In simulations of forged and stamping processes using the finite element method, load displacement paths and three-dimensional stress and strains states should be well and reliably represented. The simple tension test is a suitable and economical tool to calibrate constitutive equations with finite strains and plasticity for those simulations. A complex three-dimensional stress and strain states are developed when this test is done on rectangular bars and the necking phenomenon appears. In this work, global and local numerical results of the mechanical response of rectangular bars subjected to simple tension test obtained from two different finite element formulations are compared and discussed. To this end, Updated and Total Lagrangian formulations are used in order to get the three-dimensional stress and strain states. Geometric changes together with strain and stress distributions at the cross section where necking occurs are assessed. In particular, a detailed analysis of the effective plastic strain, stress components in axial and transverse directions and pressure, and deviatoric stress components is presented. Specific numerical results are also validated with experimental measurements comparing, in turn, the performance of the two numerical approaches used in this study.

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

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

Energy Technology Data Exchange (ETDEWEB)

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

1979-06-01

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

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

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

Formulation of stiffness equation for a three-dimensional isoparametric element with elastic-plastic material and large deformation

International Nuclear Information System (INIS)

Chang, T.Y.; Prachuktam, S.; Reich, M.

1975-01-01

The formulation of the stiffness equation for an 8 to 21 node isoparametric element with elastic-plastic material and large deformation is presented. The formulation has been implemented in a nonlinear finite element program for the analysis of three-dimensional continuums. To demonstrate the utility of the formulation, a thick-walled cylinder was analyzed and the results are compared favorably with a known solution. The element type presented can be applied not only to 3-D continuums, but also to plate or shell structures, for which degenerated isoparametric elements may be used

Finite element modeling of multilayered structures of fish scales.

Science.gov (United States)

Chandler, Mei Qiang; Allison, Paul G; Rodriguez, Rogie I; Moser, Robert D; Kennedy, Alan J

2014-12-01

The interlinked fish scales of Atractosteus spatula (alligator gar) and Polypterus senegalus (gray and albino bichir) are effective multilayered armor systems for protecting fish from threats such as aggressive conspecific interactions or predation. Both types of fish scales have multi-layered structures with a harder and stiffer outer layer, and softer and more compliant inner layers. However, there are differences in relative layer thickness, property mismatch between layers, the property gradations and nanostructures in each layer. The fracture paths and patterns of both scales under microindentation loads were different. In this work, finite element models of fish scales of A. spatula and P. senegalus were built to investigate the mechanics of their multi-layered structures under penetration loads. The models simulate a rigid microindenter penetrating the fish scales quasi-statically to understand the observed experimental results. Study results indicate that the different fracture patterns and crack paths observed in the experiments were related to the different stress fields caused by the differences in layer thickness, and spatial distribution of the elastic and plastic properties in the layers, and the differences in interface properties. The parametric studies and experimental results suggest that smaller fish such as P. senegalus may have adopted a thinner outer layer for light-weighting and improved mobility, and meanwhile adopted higher strength and higher modulus at the outer layer, and stronger interface properties to prevent ring cracking and interface cracking, and larger fish such as A. spatula and Arapaima gigas have lower strength and lower modulus at the outer layers and weaker interface properties, but have adopted thicker outer layers to provide adequate protection against ring cracking and interface cracking, possibly because weight is less of a concern relative to the smaller fish such as P. senegalus. Published by Elsevier Ltd.

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

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

International Nuclear Information System (INIS)

Peterson, A.

1984-01-01

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

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

International Nuclear Information System (INIS)

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

1976-01-01

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

Elastic properties of rigid fiber-reinforced composites

Science.gov (United States)

Chen, J.; Thorpe, M. F.; Davis, L. C.

1995-05-01

We study the elastic properties of rigid fiber-reinforced composites with perfect bonding between fibers and matrix, and also with sliding boundary conditions. In the dilute region, there exists an exact analytical solution. Around the rigidity threshold we find the elastic moduli and Poisson's ratio by decomposing the deformation into a compression mode and a rotation mode. For perfect bonding, both modes are important, whereas only the compression mode is operative for sliding boundary conditions. We employ the digital-image-based method and a finite element analysis to perform computer simulations which confirm our analytical predictions.

Nonlinear finite element modeling of corrugated board

Science.gov (United States)

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

1999-01-01

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

Peridynamic Multiscale Finite Element Methods

Energy Technology Data Exchange (ETDEWEB)

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

2015-12-01

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

Contribution to finite element modelling of airfoil aeroelastic instabilities

Directory of Open Access Journals (Sweden)

Horáček J.

2007-10-01

Full Text Available Nonlinear equations of motion for a flexibly supported rigid airfoil with additional degree of freedom for controlling of the profile motion by a trailing edge flap are derived for large vibration amplitudes. Preliminary results for numerical simulation of flow-induced airfoil vibrations in a laminar incompressible flow are presented for the NACA profile 0012 with three-degrees of freedom (vertical translation, rotation around the elastic axis and rotation of the flap. The developed numerical solution of the Navier – Stokes equations and the Arbitrary Eulerian-Lagrangian approach enable to consider the moving grid for the finite element modelling of the fluid flow around the oscillating airfoil. A sequence of numerical simulation examples is presented for Reynolds numbers up to about Re~10^5, when the system loses the aeroelastic stability, and when the large displacements of the profile and a post-critical behaviour of the system take place.

Bibliography for finite elements. [2200 references

Energy Technology Data Exchange (ETDEWEB)

Whiteman, J R [comp.

1975-01-01

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

Complex finite element sensitivity method for creep analysis

International Nuclear Information System (INIS)

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

2015-01-01

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

Finite strain logarithmic hyperelasto-plasticity with softening: a strongly non-local implicit gradient framework

NARCIS (Netherlands)

Geers, M.G.D.

2004-01-01

This paper addresses the extension of a Eulerian logarithmic finite strain hyperelasto-plasticity model in order to incorporate an isotropic plastic damage variable that leads to softening and failure of the plastic material. It is shown that a logarithmic elasto-plastic model with a strongly

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

Energy Technology Data Exchange (ETDEWEB)

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

1998-03-01

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

Application of Mass Lumped Higher Order Finite Elements

International Nuclear Information System (INIS)

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

2005-01-01

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

The finite element method its basis and fundamentals

CERN Document Server

Zienkiewicz, Olek C; Zhu, JZ

2013-01-01

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

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

International Nuclear Information System (INIS)

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

1998-03-01

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

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

OpenAIRE

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

2004-01-01

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

A 3D finite element simulation model for TBM tunnelling in soft ground

Science.gov (United States)

Kasper, Thomas; Meschke, Günther

2004-12-01

A three-dimensional finite element simulation model for shield-driven tunnel excavation is presented. The model takes into account all relevant components of the construction process (the soil and the ground water, the tunnel boring machine with frictional contact to the soil, the hydraulic jacks, the tunnel lining and the tail void grouting). The paper gives a detailed description of the model components and the stepwise procedure to simulate the construction process. The soil and the grout material are modelled as saturated porous media using a two-field finite element formulation. This allows to take into account the groundwater, the grouting pressure and the fluid interaction between the soil and slurry at the cutting face and between the soil and grout around the tail void. A Cam-Clay plasticity model is used to describe the material behaviour of cohesive soils. The cementitious grouting material in the tail void is modelled as an ageing elastic material with time-dependent stiffness and permeability. To allow for an automated computation of arbitrarily long and also curvilinear driving paths with suitable finite element meshes, the simulation procedure has been fully automated. The simulation of a tunnel advance in soft cohesive soil below the ground water table is presented and the results are compared with measurements taken from the literature. Copyright

Finite element formulation for a digital image correlation method

International Nuclear Information System (INIS)

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

2005-01-01

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

Numerical Investigation of Slab-Column Connection by Finite Element Method

International Nuclear Information System (INIS)

Akram, T.; Shaikh, M.A.; Memon, A.A.

2007-01-01

The flat slab-on-column construction subjected to high transverse stresses concentrated at the slab-column connection can lead to a non-ductile, sudden and brittle punching failure and results in the accidental collapse of flat slab buildings. The major parameters affecting the slab-column connection are the concrete strength, slab thickness, slab reinforcement and aspect ratio of column. The application of numerical methods based on the finite element theory for solving practical tasks allow to perform virtual testing of structures and explore their behavior under load and other effects in different conditions taking into account the elastic and plastic behavior of materials, appearance and development of cracks and other damages (disintegrations), and finally to simulate the failure mechanism and its consequences. In this study, the models are developed to carry out the finite element analysis of slab- column connection using ADINA (Automatic Dynamic Incremental Nonlinear Analysis) by varying the slab thickness and slab confining reinforcement and to investigate their effect on the deflection and load carrying capacity. Test results indicate that by increasing the slab thickness, the deflection and the load carrying capacity of slab-column connection increases, more over, by increasing the slab confining reinforcement, the deflection decreases where as the load carrying capacity increases. (author)

Finite element formulation for thermal stress analysis of thin reactor structures

International Nuclear Information System (INIS)

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

1978-01-01

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

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

Study of 3-D stress development in parent and twin pairs of a hexagonal close-packed polycrystal: Part II – crystal plasticity finite element modeling

International Nuclear Information System (INIS)

Abdolvand, Hamidreza; Majkut, Marta; Oddershede, Jette; Wright, Jonathan P.; Daymond, Mark R.

2015-01-01

Stress heterogeneity within each individual grain of polycrystalline Zircaloy-2 is studied using a crystal plasticity finite element (CPFE) model. For this purpose, the weighted Voronoi tessellation method is used to construct 3D geometries of more than 2600 grains based on their center-of-mass positions and volumes as measured by three-dimensional X-ray diffraction (3DXRD) microscopy. The constructed microstructure is meshed with different element densities and for different numbers of grains. Then a selected group of twin and parent pairs are studied. It is shown that the measured average stress for each grain from the 3DXRD experiment is within the stress variation zone of the grain modeled in the CPFE simulation. Also, the CPFE average stress calculation for each grain is in good agreement with the measured average stress values. It is shown that upon considering the stress variations within each grain, stresses in the parent and twin are quite different if they are plotted in the global coordinate system. However, if the stress tensor is rotated into the local coordinate system of the twin habit plane, all the stress components averaged over the presented population are close, except for the shear acting on the twin plane and the transverse stress. This result is significant as it provides information needed to model such parent-twin interactions in crystal plasticity codes

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

The development and use of a piece-wise continuous finite element for plate and shell analysis

International Nuclear Information System (INIS)

Jobson, D.A.; Knowles, J.A.

1975-01-01

The implementation of general purpose programs for the numerical analysis of plate and shell structures calls for the adoption of finite element stiffness expressions which take into account of both lateral distortion and membrane action. It is important that design-oriented programs of the above kind be perfectly general. In particular the element behaviour must be independent of the choice of base axes, and not prone either to singularities or to doubts over convergence with successive mesh refinement. The basic elements should also be mathematically isotropic and the imposition of rigid body displacements should not cause self-straining. Ideally the program should allow the assembly of a wide variety of elements, oriented in any conceivable way and of a freely chosen shape. The present paper documents a procedure for synthesising the latter from a three node primary element which satisfies the above requirements. (Auth.)

Measurement of shear strength resistance in flexion test on PUR sandwich panels: analysis of difficulties and finite element method simulation

Directory of Open Access Journals (Sweden)

J. M. Chillón Moreno

2017-09-01

Full Text Available The use of the composite construction elements combining different materials with complementary characteristics, it has spread increasingly. The sandwich panels composed by external metallic sheets that they contributes resistance and core of rigid insulating thermal foam, that provides qualities that improve the thermal comfort inside all kinds of constructions. They are in use in closings and covers that shape the surrounding one of the buildings. Of the different quality controls to which they have to surrender. In this article one proposes an improvement to the indicated one in the procedure for the determination of the resistance to the shear strength, Managing to avoid many anomalous results obtained by the utilization of rigid plates in the supports. Finally, so much the problem observed as the proposed solution, they are modeled and simulate by means of the method of finite elements.

Finite-element analysis of dynamic fracture

Science.gov (United States)

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

1976-01-01

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

The finite element method in electromagnetics

CERN Document Server

Jin, Jianming

2014-01-01

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

Finite elements for analysis and design

CERN Document Server

Akin, J E; Davenport, J H

1994-01-01

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

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

International Nuclear Information System (INIS)

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

1977-01-01

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

Shakedown analysis by finite element incremental procedures

International Nuclear Information System (INIS)

Borkowski, A.; Kleiber, M.

1979-01-01

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

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

Inside finite elements

CERN Document Server

Weiser, Martin

2016-01-01

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

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

Advanced Pavement Design: Finite Element Modeling for Rigid Pavement Joints, Report II: Model Development

National Research Council Canada - National Science Library

Hammons, Michael

1998-01-01

.... The objective of this research was to obtain data on the response of the ng'id pavement slab-joint-foundation system by conducting laboratory-scale experiments on jointed rigid pavement models...

Implicit three-dimensional finite-element formulation for the nonlinear structural response of reactor components

International Nuclear Information System (INIS)

Kulak, R.F.; Belytschko, T.B.

1975-09-01

The formulation of a finite-element procedure for the implicit transient and static analysis of plate/shell type structures in three-dimensional space is described. The triangular plate/shell element can sustain both membrane and bending stresses. Both geometric and material nonlinearities can be treated, and an elastic-plastic material law has been incorporated. The formulation permits the element to undergo arbitrarily large rotations and translations; but, in its present form it is restricted to small strains. The discretized equations of motion are obtained by a stiffness method. An implicit integration algorithm based on trapezoidal integration formulas is used to integrate the discretized equations of motion in time. To ensure numerical stability, an iterative solution procedure with equilibrium checks is used

A finite element model of rigid body structures actuated by dielectric elastomer actuators

Science.gov (United States)

Simone, F.; Linnebach, P.; Rizzello, G.; Seelecke, S.

2018-06-01

This paper presents on finite element (FE) modeling and simulation of dielectric elastomer actuators (DEAs) coupled with articulated structures. DEAs have proven to represent an effective transduction technology for the realization of large deformation, low-power consuming, and fast mechatronic actuators. However, the complex dynamic behavior of the material, characterized by nonlinearities and rate-dependent phenomena, makes it difficult to accurately model and design DEA systems. The problem is further complicated in case the DEA is used to activate articulated structures, which increase both system complexity and implementation effort of numerical simulation models. In this paper, we present a model based tool which allows to effectively implement and simulate complex articulated systems actuated by DEAs. A first prototype of a compact switch actuated by DEA membranes is chosen as reference study to introduce the methodology. The commercially available FE software COMSOL is used for implementing and coupling a physics-based dynamic model of the DEA with the external structure, i.e., the switch. The model is then experimentally calibrated and validated in both quasi-static and dynamic loading conditions. Finally, preliminary results on how to use the simulation tool to optimize the design are presented.

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.

Stability Analysis Of 3-d Conventional Pallet Rack Structures With Semi-rigid Connections

Directory of Open Access Journals (Sweden)

Kamal M. Bajoria

2009-12-01

Full Text Available This paper describe the three dimensional finite element modeling and buckling analysis of conventional pallet racking system with semi rigid connection. In this study three dimensional models of conventional pallet racking system were prepared using the finiteelement program ANSYS and finite element analysis carried out on conventional pallet racks with the 18 types of column sections developed along with semi-rigid connections. A parametric study was carried out to compare the effective length approach and the finiteelement method for accuracy and appropriateness for cold-formed steel frame design. Numerous frame elastic buckling analyses were carried out to evaluate the alignment chart and the AISI torsional-flexural buckling provisions. The parameters that influence the valueof Kx for column flexural buckling were examined in this study. The alignment chart and the AISI torsional-flexural buckling provisions, used to obtain the effective lengths and elastic buckling load of members were also evaluated. Results showed that the elastic buckling load obtained from the AISI torsional-flexural buckling provisions is generally conservative compared to the results obtained from performing frame elastic buckling analysis. Results also showed that, the effective length approach is more conservative than the finite element approach.

Investigating ASME allowable loads with finite element analyses

International Nuclear Information System (INIS)

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

1995-01-01

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

Finite element analysis of interface stress between neutron absorption coating and chop disk

International Nuclear Information System (INIS)

Tang Changliang; Zhang Xiaozhang; Jiang Lei; Dai Xingjian

2012-01-01

The performance of disk chopper is directly affected by bond strength between neutron absorption coating and chop disk. Based on the finite element analysis software ANSYS, the interface stress distribution under high speed centrifugal load was calculated, which was to investigate the effects of coating's elastic modulus, poisson ratio and coating thickness on the interfacial stress distribution. The results show that soft and tough coating can reduce the peak stress effectively, and coating thickness reducing is helpful to avoid the plastic failure of opening in the disk under high speed centrifugal load. (authors)

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

International Nuclear Information System (INIS)

Baker, A.R.

1982-07-01

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

A geometric toolbox for tetrahedral finite element partitions

NARCIS (Netherlands)

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

2011-01-01

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

Mixed finite element simulations in two-dimensional groundwater flow problems

International Nuclear Information System (INIS)

Kimura, Hideo

1989-01-01

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

3D Finite Element Simulation of Pellet-Cladding Mechanical Interaction

International Nuclear Information System (INIS)

Seo, Sang Kyu; Lee, Sung Uk; Lee, Eun Ho; Yang, Dong Yol; Kim, Hyo Chan; Yang, Dong Yol

2016-01-01

In a nuclear power plant, the fuel assembly, which is composed of fuel rods, burns, and the high temperature can generate power. The fuel rod consists of pellets and a cladding that covers the pellets. It is important to understand the pellet-cladding mechanical interaction with regard to nuclear safety. This paper proposes simulation of the PCMI. The gap between the pellets and the cladding, and the contact pressure are very important for conducting thermal analysis. Since the gap conductance is not known, it has to be determined by a suitable method. This paper suggests a solution. In this study, finite element (FE) contact analysis is conducted considering thermal expansion of the pellets. As the contact causes plastic deformation, this aspect is considered in the analysis. A 3D FE module is developed to analyze the PCMI using FORTRAN 90. The plastic deformation due to the contact between the pellets and the cladding is the major physical phenomenon. The simple analytical solution of a cylinder is proposed and compared with the fuel rod performance code results

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

International Nuclear Information System (INIS)

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

1977-01-01

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

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.

Coupled Vortex-Lattice Flight Dynamic Model with Aeroelastic Finite-Element Model of Flexible Wing Transport Aircraft with Variable Camber Continuous Trailing Edge Flap for Drag Reduction

Science.gov (United States)

Nguyen, Nhan; Ting, Eric; Nguyen, Daniel; Dao, Tung; Trinh, Khanh

2013-01-01

This paper presents a coupled vortex-lattice flight dynamic model with an aeroelastic finite-element model to predict dynamic characteristics of a flexible wing transport aircraft. The aircraft model is based on NASA Generic Transport Model (GTM) with representative mass and stiffness properties to achieve a wing tip deflection about twice that of a conventional transport aircraft (10% versus 5%). This flexible wing transport aircraft is referred to as an Elastically Shaped Aircraft Concept (ESAC) which is equipped with a Variable Camber Continuous Trailing Edge Flap (VCCTEF) system for active wing shaping control for drag reduction. A vortex-lattice aerodynamic model of the ESAC is developed and is coupled with an aeroelastic finite-element model via an automated geometry modeler. This coupled model is used to compute static and dynamic aeroelastic solutions. The deflection information from the finite-element model and the vortex-lattice model is used to compute unsteady contributions to the aerodynamic force and moment coefficients. A coupled aeroelastic-longitudinal flight dynamic model is developed by coupling the finite-element model with the rigid-body flight dynamic model of the GTM.

Impact of new computing systems on finite element computations

International Nuclear Information System (INIS)

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

1983-01-01

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

Finite size effects of a pion matrix element

International Nuclear Information System (INIS)

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

2004-01-01

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

Finite elements in fracture mechanics theory, numerics, applications

CERN Document Server

Kuna, Meinhard

2013-01-01

Fracture mechanics has established itself as an important discipline of growing interest to those working to assess the safety, reliability and service life of engineering structures and materials. In order to calculate the loading situation at cracks and defects, nowadays numerical techniques like finite element method (FEM) have become indispensable tools for a broad range of applications. The present monograph provides an introduction to the essential concepts of fracture mechanics, its main goal being to procure the special techniques for FEM analysis of crack problems, which have to date only been mastered by experts. All kinds of static, dynamic and fatigue fracture problems are treated in two- and three-dimensional elastic and plastic structural components. The usage of the various solution techniques is demonstrated by means of sample problems selected from practical engineering case studies. The primary target group includes graduate students, researchers in academia and engineers in practice.

Evaluation of radiation damping using 3-D finite element models

International Nuclear Information System (INIS)

Vaughan, D.K.; Isenberg, J.

1983-01-01

The paper presents an analytic approach which is being used to quantify the contribution of radiation damping to overall system damping. The approach uses three-dimensional finite element techniques and can easily include details of site geology, foundation shape, and embedment depth. The approach involves performing free vibration response analyses for each soil-structure interaction (SSI) mode of interest. The structural model is specified without damping and, consequently, amplitude decay of the structure's free vibration response is a measure of the radiation damping characteristics of the soil-structure system for the particular deformational mode being investigated. The computational approach developed is highly efficient in order to minimize the impact of including three-dimensional geometry within the model. A new finite element code, FLEX, has been developed to represent the soil continuum. FLEX uses a highly optimized explicit time integration algorithm which takes advantage of parallel processing on vector machines, such as the CRAY 1 computer. A modal representation of the superstructure is used in combination with a substructuring approach to solve for the coupled response of the soil-structure system. This requires solving for numerical Green's functions for each degree-of-freedom of the foundation (assumed rigid). Once computed for a particular site and foundation, these Green's functions may be used within a convolution integral to represent the continuum forces on the foundation for any free vibration SSI response computation of any superstructure model. This analytic approach is applied to an investigation of the radiation damping coefficients for the first two fundamental SSI modes of the HDR containment structure. (orig./HP)

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.

Indentation of elastically soft and plastically compressible solids

DEFF Research Database (Denmark)

Needleman, A.; Tvergaard, Viggo; Van der Giessen, E.

2015-01-01

rapidly for small deviations from plastic incompressibility and then decreases rather slowly for values of the plastic Poisson's ratio less than 0.25. For both soft elasticity and plastic compressibility, the main reason for the lower values of indentation hardness is related to the reduction......The effect of soft elasticity, i.e., a relatively small value of the ratio of Young's modulus to yield strength and plastic compressibility on the indentation of isotropically hardening elastic-viscoplastic solids is investigated. Calculations are carried out for indentation of a perfectly sticking...... rigid sharp indenter into a cylinder modeling indentation of a half space. The material is characterized by a finite strain elastic-viscoplastic constitutive relation that allows for plastic as well as elastic compressibility. Both soft elasticity and plastic compressibility significantly reduce...

Verification of Orthogrid Finite Element Modeling Techniques

Science.gov (United States)

Steeve, B. E.

1996-01-01

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

Review of Tomographic Imaging using Finite Element Method

Directory of Open Access Journals (Sweden)

Mohd Fua’ad RAHMAT

2011-12-01

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

Magnetic materials and 3D finite element modeling

CERN Document Server

Bastos, Joao Pedro A

2014-01-01

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

New mixed finite-element methods

International Nuclear Information System (INIS)

Franca, L.P.

1987-01-01

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

Evaluation of Strip Footing Bearing Capacity Built on the Anthropogenic Embankment by Random Finite Element Method

Science.gov (United States)

Pieczynska-Kozlowska, Joanna

2014-05-01

One of a geotechnical problem in the area of Wroclaw is an anthropogenic embankment layer delaying to the depth of 4-5m, arising as a result of historical incidents. In such a case an assumption of bearing capacity of strip footing might be difficult. The standard solution is to use a deep foundation or foundation soil replacement. However both methods generate significant costs. In the present paper the authors focused their attention on the influence of anthropogenic embankment variability on bearing capacity. Soil parameters were defined on the basis of CPT test and modeled as 2D anisotropic random fields and the assumption of bearing capacity were made according deterministic finite element methods. Many repeated of the different realizations of random fields lead to stable expected value of bearing capacity. The algorithm used to estimate the bearing capacity of strip footing was the random finite element method (e.g. [1]). In traditional approach of bearing capacity the formula proposed by [2] is taken into account. qf = c'Nc + qNq + 0.5γBN- γ (1) where: qf is the ultimate bearing stress, cis the cohesion, qis the overburden load due to foundation embedment, γ is the soil unit weight, Bis the footing width, and Nc, Nq and Nγ are the bearing capacity factors. The method of evaluation the bearing capacity of strip footing based on finite element method incorporate five parameters: Young's modulus (E), Poisson's ratio (ν), dilation angle (ψ), cohesion (c), and friction angle (φ). In the present study E, ν and ψ are held constant while c and φ are randomized. Although the Young's modulus does not affect the bearing capacity it governs the initial elastic response of the soil. Plastic stress redistribution is accomplished using a viscoplastic algorithm merge with an elastic perfectly plastic (Mohr - Coulomb) failure criterion. In this paper a typical finite element mesh was assumed with 8-node elements consist in 50 columns and 20 rows. Footings width B

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

A finite element elastic-plastic analysis of residual stresses due to clad welding in reactor vessels

International Nuclear Information System (INIS)

Buchalet, C.; Riccardella, P.C.

1972-01-01

Residual stresses due to weld deposited cladding on the inside of a typical Westinghouse pressurized water reactor vessel are investigated using an axisymmetric finite element elastic-plastic analysis. At the beginning of the analysis, one head of the weld cladding is assumed to lie on the reactor vessel wall at melting temperature (2600degF), but in the solid phase, while the vessel remains at 300degF (preheat temperature). All material properties used in the calculations are taken as temperature-dependent. Temperature profiles are obtained in the cladding and base metal at several discrete time intervals. These temperatures profiles are used to obtain the stress distribution for the same time intervals. Residual hoop tensile stresses of approximately 25 ksi were found to exist in the cladding. Peak tensile stresses in the hoop direction occur in the base metal near the cladding interface and reach a value of 60 ksi at the end of the transient. The tensile stress decreases very rapidly through the thickness of the base metal and becomes insignificant at about two inches from the inside surface. In order to lower residual stresses, a post-weld heat treatment is performed by uniformly heating the vessel to 1100degF, holding at that temperature for a specified period of time and then cooling slowly. The analysis shows that after this treatment, the peak stresses in the base metal decrease from 60 ksi to 32 ksi, while the stress in the cladding does not change significantly. (author)

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

International Nuclear Information System (INIS)

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

1993-01-01

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

Finite element simulation of the T-shaped ECAP processing of round samples

Science.gov (United States)

Shaban Ghazani, Mehdi; Fardi-Ilkhchy, Ali; Binesh, Behzad

2018-05-01

Grain refinement is the only mechanism that increases the yield strength and toughness of the materials simultaneously. Severe plastic deformation is one of the promising methods to refine the microstructure of materials. Among different severe plastic deformation processes, the T-shaped equal channel angular pressing (T-ECAP) is a relatively new technique. In the present study, finite element analysis was conducted to evaluate the deformation behavior of metals during T-ECAP process. The study was focused mainly on flow characteristics, plastic strain distribution and its homogeneity, damage development, and pressing force which are among the most important factors governing the sound and successful processing of nanostructured materials by severe plastic deformation techniques. The results showed that plastic strain is localized in the bottom side of sample and uniform deformation cannot be possible using T-ECAP processing. Friction coefficient between sample and die channel wall has a little effect on strain distributions in mirror plane and transverse plane of deformed sample. Also, damage analysis showed that superficial cracks may be initiated from bottom side of sample and their propagation will be limited due to the compressive state of stress. It was demonstrated that the V shaped deformation zone are existed in T-ECAP process and the pressing load needed for execution of deformation process is increased with friction.

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

MEL finite element analysis of water-shell interactions in the context of a PWR-LOCA

International Nuclear Information System (INIS)

Verbiese, S.; Vrije Universiteit Brussels; Goethem, G. van

1979-01-01

In the framework of the computational effort engaged towards and understanding of the transient dynamic fluid-structure phenomena taking place in the very first instants of the PWR loss-of-coolant-accident, before the ebullition crisis and the subsequent two-phase flow, two finite element programs have been selected and coupled to describe this class of events in pressure vessels undergoing moderate plastic deformations. Water is modeled by a compressible inviscid Eulerian (bulk of the fluid) - mixed Eulerian-Lagrangian MEL (boundary elements in contact with the moving structure) program. For the shells a convected coordinates elastic-plastic structural code (EURDYN I) is used. A 1-D discussion on the MEL integration scheme is presented, as well as a flow chart of the combined program. Emphases is placed, during the present calculations limited to very simple axisymmetric configurations, upon the computational aspects in dealing with the interaction of both media at the fluid-structure interface, such as weak code coupling, subcycling and pressure relaxation. (orig.)

An introduction to the UNCLE finite element scheme

International Nuclear Information System (INIS)

Enderby, J.A.

1983-01-01

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

An introduction to the UNCLE finite element scheme

Energy Technology Data Exchange (ETDEWEB)

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

1983-05-01

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

A particle finite element method for machining simulations

Science.gov (United States)

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

2014-07-01

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

A finite element primer for beginners the basics

CERN Document Server

Zohdi, Tarek I

2014-01-01

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

On the reliability of finite element solutions

International Nuclear Information System (INIS)

Prasad, K.S.R.K.

1975-01-01

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

An evaluation on the effect of reversed plastic zone on the fatigue crack opening behavior under 2-D plane stress

International Nuclear Information System (INIS)

Choi, Hyeon Chang

2005-01-01

The relationship between fatigue crack opening behavior and the reversed plastic zone sizes is studied. An elastic-plastic Finite Element Analysis (FEA) is performed to examine the opening behavior of fatigue crack, where the contact elements are used in the mesh of the crack tip area. The smaller element size than reversed plastic zone size is used for evaluating the distribution of reversed plastic zone. In the author's previous results the FEA could predict the crack opening level, which crack tip elements were in proportion to the theoretical reversed plastic zone size. It is found that the calculated reversed plastic zone size is related to the theoretical reversed plastic zone size and crack opening level. The calculated reversed plastic zone sizes are almost equal to the reversed plastic zone considering crack opening level obtained by experimental results. It can be possible to predict the crack opening level from the reversed plastic zone size calculated by finite element method. We find that the experimental crack opening levels correspond with the opening values of contact nodes on the calculated reversed plastic zone of finite element simulation

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

International Nuclear Information System (INIS)

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

2011-01-01

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

Automation of finite element methods

CERN Document Server

Korelc, Jože

2016-01-01

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

Finite Element Based Design and Optimization for Piezoelectric Accelerometers

DEFF Research Database (Denmark)

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

1998-01-01

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

Rigid inclusions-Comparison between analytical and numerical methods; Inclusiones rigidas-Analisis comparativo de los metodos analiticos y numericos

Energy Technology Data Exchange (ETDEWEB)

Gomez Perez, R.; Melentijevic, S.

2014-07-01

This paper compares different analytical methods for analysis of rigid inclusions with finite element modeling. First of all, the load transfer in the distribution layer is analyzed for its different thicknesses and different inclusion grids to define the range between results obtained by analytical and numerical methods. The interaction between the soft soil and the inclusion in the estimation of settlements is studied as well. Considering different stiffness of the soft soil, settlements obtained analytical and numerically are compared. The influence of the soft soil modulus of elasticity on the neutral point depth was also performed by finite elements. This depth has a great importance for the definition of the total length of rigid inclusion. (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...

Finite Element Simulation and Experimental Verification of Internal Stress of Quenched AISI 4140 Cylinders

Science.gov (United States)

Liu, Yu; Qin, Shengwei; Hao, Qingguo; Chen, Nailu; Zuo, Xunwei; Rong, Yonghua

2017-03-01

The study of internal stress in quenched AISI 4140 medium carbon steel is of importance in engineering. In this work, the finite element simulation (FES) was employed to predict the distribution of internal stress in quenched AISI 4140 cylinders with two sizes of diameter based on exponent-modified (Ex-Modified) normalized function. The results indicate that the FES based on Ex-Modified normalized function proposed is better consistent with X-ray diffraction measurements of the stress distribution than FES based on normalized function proposed by Abrassart, Desalos and Leblond, respectively, which is attributed that Ex-Modified normalized function better describes transformation plasticity. Effect of temperature distribution on the phase formation, the origin of residual stress distribution and effect of transformation plasticity function on the residual stress distribution were further discussed.

Biomechanical Analysis of Human Abdominal Impact Responses and Injuries through Finite Element Simulations of a Full Human Body Model.

Science.gov (United States)

Ruan, Jesse S; El-Jawahri, Raed; Barbat, Saeed; Prasad, Priya

2005-11-01

Human abdominal response and injury in blunt impacts was investigated through finite element simulations of cadaver tests using a full human body model of an average-sized adult male. The model was validated at various impact speeds by comparing model responses with available experimental cadaver test data in pendulum side impacts and frontal rigid bar impacts from various sources. Results of various abdominal impact simulations are presented in this paper. Model-predicted abdominal dynamic responses such as force-time and force-deflection characteristics, and injury severities, measured by organ pressures, for the simulated impact conditions are presented. Quantitative results such as impact forces, abdominal deflections, internal organ stresses have shown that the abdomen responded differently to left and right side impacts, especially in low speed impact. Results also indicated that the model exhibited speed sensitive response characteristics and the compressibility of the abdomen significantly influenced the overall impact response in the simulated impact conditions. This study demonstrates that the development of a validated finite element human body model can be useful for abdominal injury assessment. Internal organ injuries, which are difficult to detect in experimental studies with human cadavers due to the difficulty of instrumentation, may be more easily identified with a validated finite element model through stress-strain analysis.

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.

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

International Nuclear Information System (INIS)

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

2011-01-01

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

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

Science.gov (United States)

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

2004-01-01

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

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

African Journals Online (AJOL)

user

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

Finite element analysis of tibial fractures

DEFF Research Database (Denmark)

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

2010-01-01

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

Generalized multiscale finite element method. Symmetric interior penalty coupling

KAUST Repository

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

2013-01-01

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

Generalized multiscale finite element method. Symmetric interior penalty coupling

KAUST Repository

Efendiev, Yalchin R.

2013-12-01

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

Finite element analysis of fatigue crack closure under plane strain state

International Nuclear Information System (INIS)

Lee, Hak Joo; Kang, Jae Youn; Song, Ji Ho

2004-01-01

An elastic-plastic finite element analysis of fatigue crack closure is performed for plane strain conditions. The stabilization behavior of crack opening level and the effect of mesh size on the crack opening stress are investigated. In order to obtain a stabilized crack opening level for plane strain conditions, the crack must be advanced through approximately four times the initial monotonic plastic zone. The crack opening load tends to increase with the decrease of mesh size. The mesh size nearly equal to the theoretical plane strain cyclic plastic zone size may provide reasonable numerical results comparable with experimental crack opening data. The crack opening behavior is influenced by the crack growth increment and discontinuous opening behavior is observed. A procedure to predict the most appropriate mesh size for different stress ratio is suggested. Crack opening loads predicted by the FE analysis based on the procedure suggested resulted in good agreement with experimental ones within the error of 5 %. Effect of the distance behind the crack tip on the crack opening load determined by the ASTM compliance offset method based on the load-displacement relation and by the rotational offset method based on the load-differential displacement relation is investigated. Optimal gage location and method to determine the crack opening load is suggested

On the finite element modeling of the asymmetric cracked rotor

Science.gov (United States)

AL-Shudeifat, Mohammad A.

2013-05-01

The advanced phase of the breathing crack in the heavy duty horizontal rotor system is expected to be dominated by the open crack state rather than the breathing state after a short period of operation. The reason for this scenario is the expected plastic deformation in crack location due to a large compression stress field appears during the continuous shaft rotation. Based on that, the finite element modeling of a cracked rotor system with a transverse open crack is addressed here. The cracked rotor with the open crack model behaves as an asymmetric shaft due to the presence of the transverse edge crack. Hence, the time-varying area moments of inertia of the cracked section are employed in formulating the periodic finite element stiffness matrix which yields a linear time-periodic system. The harmonic balance method (HB) is used for solving the finite element (FE) equations of motion for studying the dynamic behavior of the system. The behavior of the whirl orbits during the passage through the subcritical rotational speeds of the open crack model is compared to that for the breathing crack model. The presence of the open crack with the unbalance force was found only to excite the 1/2 and 1/3 of the backward critical whirling speed. The whirl orbits in the neighborhood of these subcritical speeds were found to have nearly similar behavior for both open and breathing crack models. While unlike the breathing crack model, the subcritical forward whirling speeds have not been observed for the open crack model in the response to the unbalance force. As a result, the behavior of the whirl orbits during the passage through the forward subcritical rotational speeds is found to be enough to distinguish the breathing crack from the open crack model. These whirl orbits with inner loops that appear in the neighborhood of the forward subcritical speeds are then a unique property for the breathing crack model.

Finite Element Simulation of Shot Peening: Prediction of Residual Stresses and Surface Roughness

Science.gov (United States)

Gariépy, Alexandre; Perron, Claude; Bocher, Philippe; Lévesque, Martin

Shot peening is a surface treatment that consists of bombarding a ductile surface with numerous small and hard particles. Each impact creates localized plastic strains that permanently stretch the surface. Since the underlying material constrains this stretching, compressive residual stresses are generated near the surface. This process is commonly used in the automotive and aerospace industries to improve fatigue life. Finite element analyses can be used to predict residual stress profiles and surface roughness created by shot peening. This study investigates further the parameters and capabilities of a random impact model by evaluating the representative volume element and the calculated stress distribution. Using an isotropic-kinematic hardening constitutive law to describe the behaviour of AA2024-T351 aluminium alloy, promising results were achieved in terms of residual 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.

Parametric study of unconstrained high-pressure torsion- Finite element analysis

International Nuclear Information System (INIS)

Halloumi, A; Busquet, M; Descartes, S

2014-01-01

The high-pressure torsion (HPT) experiments have been investigated numerically. An axisymmetric model with twist was developed with commercial finite element software (Abaqus) to study locally the specificity of the stress and strain history within the transformed layers produced during HPT processing. The material local behaviour law in the plastic domain was modelled. A parametric study highlights the role of the imposed parameters (friction coefficient at the interfaces anvil surfaces/sample, imposed pressure) on the stress/strain distribution in the sample bulk for two materials: ultra-high purity iron and steel grade R260. The present modelling provides a tool to investigate and to analyse the effect of pressure and friction on the local stress and strain history during the HPT process and to couple with experimental results

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

Introduction to Computational Plasticity

International Nuclear Information System (INIS)

Hartley, P

2006-01-01

matrices with a clear explanation of their development, is a recurring, and commendable, feature of the book, which provides an invaluable introduction for those new to the subject. The chapter moves on from time-independent behaviour to introduce viscoplasticity and creep. Chapter three takes the theories of deformation another stage further to consider the problems associated with large deformation in which an important concept is the separation of the phenomenon into material stretch and rotation. The latter is crucial to allow correct measures of strain and stress to be developed in which the effects of rigid body rotation do not contribute to these variables. Hence, the introduction of 'objective' measures for stress and strain. These are described with reference to deformation gradients, which are clearly explained; however, the introduction of displacement gradients passes with little comment, although velocity gradients appear later in the chapter. The interpretation of different strain measures, e.g. Green--Lagrange and Almansi, is covered briefly, followed by a description of the spin tensor and its use in developing the objective Jaumann rate of stress. It is tempting here to suggest that a more complete description should be given together with other measures of strain and stress, of which there are several, but there would be a danger of changing the book from an 'introduction' to a more comprehensive text, and examples of such exist already. Chapter four begins the process of developing the plasticity theories into a form suitable for inclusion in the finite-element method. The starting point is Hamilton's principle for equilibrium of a dynamic system. A very brief introduction to the finite-element method is then given, followed by the finite-element equilibrium equations and a description of how they are incorporated into Hamilton's principle. A useful clarification is provided by comparing tensor notation and the form normally used in finite-element

Elastic–plastic adhesive contact of non-Gaussian rough surfaces

Indian Academy of Sciences (India)

Adhesion; asymmetric roughness; elastic–plastic contact; non-Gaussian rough surfaces. ... model of contact deformation that is based on accurate Finite Element Analysis (FEA) of an elastic–plastic single asperity contact. ... Sadhana | News.

Modeling plasticity by non-continuous deformation

Science.gov (United States)

Ben-Shmuel, Yaron; Altus, Eli

2017-10-01

Plasticity and failure theories are still subjects of intense research. Engineering constitutive models on the macroscale which are based on micro characteristics are very much in need. This study is motivated by the observation that continuum assumptions in plasticity in which neighbour material elements are inseparable at all-time are physically impossible, since local detachments, slips and neighbour switching must operate, i.e. non-continuous deformation. Material microstructure is modelled herein by a set of point elements (particles) interacting with their neighbours. Each particle can detach from and/or attach with its neighbours during deformation. Simulations on two- dimensional configurations subjected to uniaxial compression cycle are conducted. Stochastic heterogeneity is controlled by a single "disorder" parameter. It was found that (a) macro response resembles typical elasto-plastic behaviour; (b) plastic energy is proportional to the number of detachments; (c) residual plastic strain is proportional to the number of attachments, and (d) volume is preserved, which is consistent with macro plastic deformation. Rigid body displacements of local groups of elements are also observed. Higher disorder decreases the macro elastic moduli and increases plastic energy. Evolution of anisotropic effects is obtained with no additional parameters.

Stress analysis of disconnected structures in contact through finite element gaps

International Nuclear Information System (INIS)

Stadter, J.T.; Weiss, R.O.

1976-07-01

A numerical procedure is presented for analyzing thermal stress problems of disconnected structures in contact across separations or gaps. The new procedure is called SAASGAPS, an adaptation of the basic SAAS III computer program. The SAAS program uses the finite element method and allows analyses of plane and axisymmetric bodies with temperature dependent material properties, subject to thermal and mechanical loads. A secant modulus approach with a bilinear stress-strain curve is used for elastic-plastic problems. The SAASGAPS version contains all of the features of the original SAAS program. A special gap element is used together with a stress invariance principle to model the contact process. The iterative procedure implemented in SAASGAPS is described. Results are discussed for five problems involving frictionless contact. Two of these problems are associated with the thermal stress analysis of the heat shield for the Multi-Hundred Watt Radioisotope Thermoelectric Generator. Input instructions for the program are described in an appendix

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.

Inelastic analysis of finite length and depth cracked tubes

International Nuclear Information System (INIS)

Reich, M.; Gardner, D.; Prachuktam, S.; Chang, T.Y.

1977-01-01

Steam generator tube failure can at times result in reactor safety problems and subsequent premature reactor shutdown. This paper concerns itself with the prediction of the failure pressures for typical PWR steam generator tubes with longitudinal finite length and finite depth cracks. Only local plastic overload failure is considered since the material is non-notch sensitive. Non-linear finite element analyses are carried out to determine the burst pressures of steam generator tubes containing longitudinal cracks located on the outer surface of the tubes. The non-linearities considered herein include elastic-plastic material behaviour and large deformations. A non-proprietary general purpose non-linear finite element program, NFAP was adopted for the analysis. Due to the asymmetric nature of the cracks, two-dimensional as well as three-dimensional finite element analyses, were performed. The analysis clearly shows that for short cracks axial effects play a significant role. For long cracks, they are not important since two-dimensional conditions predominate and failure is governed by circumferential or hoop stress conditions. (Auth.)

A 3D elasto-plastic soil model for lateral buckling analysis

DEFF Research Database (Denmark)

Hededal, Ole; Strandgaard, Torsten

2008-01-01

Modeling the lay-down of pipelines and subsequently the in- service conditions for a pipeline involves definition of a pipe-soil interaction model. A generalized true 3D elasto-plastic spring element based on an anisotropic hardening/degradation model for sliding is presented. The basis...... for the model is the elasto-plastic framework. A generic format is selected, allowing different yield criteria and flow rules to be implemented in a simple way. The model complies to a finite element format allowing it to be directly implemented into a standard finite element code. Examples demonstrating...

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

Fracture toughness evaluation of elastic-plastic J-integral for high temperature components of gas turbine in power plants

International Nuclear Information System (INIS)

Chung, Nam Yong; Kim, Moon Young; Kim, Jong Woo

1999-01-01

In the study, the analysis of elastic-plastic J-integral was performed in high temperature components for gas turbine based on elastic-plastic fracture mechanics. It had been operated on the range of about 700 deg C and degraded by high temperature. It was tested for material properties of used component because of material properties changing at high temperature condition. The elastic-plastic fracture mechanics parameter, J is obtained with finite element method. A method is suggested which determines J Ic applying analysis of elastic-plastic finite element method and results of experimental load-displacements with CT specimen. It is also investigated that J-integral is applied for the elastic-plastic analysis in high temperature components. The elastic-plastic fracture toughness. J Ic determined by finite element was obtained with high accuracy using the experimental method.=20

Elastic-plastic analysis of AS4/PEEK composite laminate using a one-parameter plasticity model

Science.gov (United States)

Sun, C. T.; Yoon, K. J.

1992-01-01

A one-parameter plasticity model was shown to adequately describe the plastic deformation of AS4/PEEK (APC-2) unidirectional thermoplastic composite. This model was verified further for unidirectional and laminated composite panels with and without a hole. The elastic-plastic stress-strain relations of coupon specimens were measured and compared with those predicted by the finite element analysis using the one-parameter plasticity model. The results show that the one-parameter plasticity model is suitable for the analysis of elastic-plastic deformation of AS4/PEEK composite laminates.

Numerical experiment on finite element method for matching data

International Nuclear Information System (INIS)

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

1993-03-01

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

Analytical and finite element modeling of grounding systems

Energy Technology Data Exchange (ETDEWEB)

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

2007-07-01

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

Nonlinear finite element analysis of liquid sloshing in complex vehicle motion scenarios

Science.gov (United States)

Nicolsen, Brynne; Wang, Liang; Shabana, Ahmed

2017-09-01

The objective of this investigation is to develop a new total Lagrangian continuum-based liquid sloshing model that can be systematically integrated with multibody system (MBS) algorithms in order to allow for studying complex motion scenarios. The new approach allows for accurately capturing the effect of the sloshing forces during curve negotiation, rapid lane change, and accelerating and braking scenarios. In these motion scenarios, the liquid experiences large displacements and significant changes in shape that can be captured effectively using the finite element (FE) absolute nodal coordinate formulation (ANCF). ANCF elements are used in this investigation to describe complex mesh geometries, to capture the change in inertia due to the change in the fluid shape, and to accurately calculate the centrifugal forces, which for flexible bodies do not take the simple form used in rigid body dynamics. A penalty formulation is used to define the contact between the rigid tank walls and the fluid. A fully nonlinear MBS truck model that includes a suspension system and Pacejka's brush tire model is developed. Specified motion trajectories are used to examine the vehicle dynamics in three different scenarios - deceleration during straight-line motion, rapid lane change, and curve negotiation. It is demonstrated that the liquid sloshing changes the contact forces between the tires and the ground - increasing the forces on certain wheels and decreasing the forces on other wheels. In cases of extreme sloshing, this dynamic behavior can negatively impact the vehicle stability by increasing the possibility of wheel lift and vehicle rollover.

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

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

Science.gov (United States)

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

2017-12-01

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

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.

3D finite element analysis on crack-tip plastic zone

African Journals Online (AJOL)

user

International Journal of Engineering, Science and Technology ... directed in this investigation to achieve a few numerical estimations of 3D plastic zone shape ..... crack front, using domain integral method the software (ABAQUS) automatically.

Inelastic analysis of finite length and depth cracked tubes

International Nuclear Information System (INIS)

Reich, M.; Gardner, D.; Prachuktam, S.; Chang, T.Y.

1977-01-01

Steam generator tube failure can at times result in reactor safety problems and subsequent premature reactor shutdowns. This paper concerns itself with the prediction of the failure pressures for typical PWR steam generator tubes with longitudinal finite length and finite depth cracks. Only local plastic overload failure is considered since the material is non-notch sensitive. Non-linear finite element analyses are carried out to determine the burst pressures of steam generator tubes containing longitudinal cracks located on the outer surface of the tubes. The non-linearities considered herein include elastic-plastic material behavior and large deformations. A non-proprietary general purpose non-linear finite element program, NFAP was adopted for the analysis. Due to the asymmetric nature of the cracks, two-dimensional, as well as three-dimensional finite element analyses, were performed. The two-dimensional element and its formulations are similar to those of NONSAP. The three-dimensional isoparametric element with elastic-plastic material characteristics together with the large deformation formulations used in NFAP are described in the Report BNL-20684. The numerical accuracy of the program was investigated and checked with known solutions of benchmark problems. In addition to the three-dimensional element which was specifically inserted into NFAP for this problem, other features such as direct pressure inputs for isoparametric elements, automatic load increment adjustments for convergent non-linear solutions, and automatic bandwidth reduction schemes are incorporated into the program thus allowing for a more economical evaluation of three-dimensional inelastic analysis. In summary the analysis clearly shows that for short cracks axial effects play a significant role. For long cracks, they are not important since two-dimensional conditions predominate and failure is governed by circumferential or hoop stress conditions

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

Energy Technology Data Exchange (ETDEWEB)

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

1999-03-01

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

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

International Nuclear Information System (INIS)

Ishiguro, Misako; Higuchi, Kenji

1983-01-01

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

A finite element analysis of novel vented dental abutment geometries for cement-retained crown restorations.

Science.gov (United States)

Rodriguez, Lucas C; Saba, Juliana N; Meyer, Clark A; Chung, Kwok-Hung; Wadhwani, Chandur; Rodrigues, Danieli C

2016-11-01

Recent literature indicates that the long-term success of dental implants is, in part, attributed to how dental crowns are attached to their associated implants. The commonly utilized method for crown attachment - cementation, has been criticized because of recent links between residual cement and peri-implant disease. Residual cement extrusion from crown-abutment margins post-crown seating is a growing concern. This study aimed at (1) identifying key abutment features, which would improve dental cement flow characteristics, and (2) understanding how these features would impact the mechanical stability of the abutment under functional loads. Computational fluid dynamic modeling was used to evaluate cement flow in novel abutment geometries. These models were then evaluated using 3D-printed surrogate models. Finite element analysis also provided an understanding of how the mechanical stability of these abutments was altered after key features were incorporated into the geometry. The findings demonstrated that the key features involved in improved venting of the abutment during crown seating were (1) addition of vents, (2) diameter of the vents, (3) location of the vents, (4) addition of a plastic screw insert, and (5) thickness of the abutment wall. This study culminated in a novel design for a vented abutment consisting of 8 vents located radially around the abutment neck-margin plus a plastic insert to guide the cement during seating and provide retrievability to the abutment system.Venting of the dental abutment has been shown to decrease the risk of undetected residual dental cement post-cement-retained crown seating. This article will utilize a finite element analysis approach toward optimizing dental abutment designs for improved dental cement venting. Features investigated include (1) addition of vents, (2) diameter of vents, (3) location of vents, (4) addition of plastic screw insert, and (5) thickness of abutment wall.

A finite element analysis of novel vented dental abutment geometries for cement‐retained crown restorations

Science.gov (United States)

Rodriguez, Lucas C.; Saba, Juliana N.; Meyer, Clark A.; Chung, Kwok‐Hung; Wadhwani, Chandur

2016-01-01

Abstract Recent literature indicates that the long‐term success of dental implants is, in part, attributed to how dental crowns are attached to their associated implants. The commonly utilized method for crown attachment – cementation, has been criticized because of recent links between residual cement and peri‐implant disease. Residual cement extrusion from crown‐abutment margins post‐crown seating is a growing concern. This study aimed at (1) identifying key abutment features, which would improve dental cement flow characteristics, and (2) understanding how these features would impact the mechanical stability of the abutment under functional loads. Computational fluid dynamic modeling was used to evaluate cement flow in novel abutment geometries. These models were then evaluated using 3D‐printed surrogate models. Finite element analysis also provided an understanding of how the mechanical stability of these abutments was altered after key features were incorporated into the geometry. The findings demonstrated that the key features involved in improved venting of the abutment during crown seating were (1) addition of vents, (2) diameter of the vents, (3) location of the vents, (4) addition of a plastic screw insert, and (5) thickness of the abutment wall. This study culminated in a novel design for a vented abutment consisting of 8 vents located radially around the abutment neck‐margin plus a plastic insert to guide the cement during seating and provide retrievability to the abutment system.Venting of the dental abutment has been shown to decrease the risk of undetected residual dental cement post‐cement‐retained crown seating. This article will utilize a finite element analysis approach toward optimizing dental abutment designs for improved dental cement venting. Features investigated include (1) addition of vents, (2) diameter of vents, (3) location of vents, (4) addition of plastic screw insert, and (5) thickness of abutment wall. PMID

Finite element approximation to the even-parity transport equation

International Nuclear Information System (INIS)

Lewis, E.E.

1981-01-01

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

Linear finite element method for one-dimensional diffusion problems

Energy Technology Data Exchange (ETDEWEB)

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

2011-07-01

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

Quadrature representation of finite element variational forms

DEFF Research Database (Denmark)

Ølgaard, Kristian Breum; Wells, Garth N.

2012-01-01

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

Effects of the Strain Rate Sensitivity and Strain Hardening on the Saturated Impulse of Plates

Directory of Open Access Journals (Sweden)

Ling Zhu

Full Text Available Abstract This paper studies the stiffening effects of the material strain rate sensitivity and strain hardening on the saturated impulse of elastic, perfectly plastic plates. Finite element (FE code ABAQUS is employed to simulate the elastoplastic response of square plates under rectangular pressure pulse. Rigid-plastic analyses for saturated impulse, which consider strain rate sensitivity and strain hardening, are conducted. Satisfactory agreement between the finite element models (FEM and predictions of the rigid-plastic analysis is obtained, which verifies that the proposed rigid-plastic methods are effective to solve the problem including strain rate sensitivity and strain hardening. The quantitative results for the scale effect of the strain rate sensitivity are given. The results for the stiffening effects suggest that two general stiffening factors n 1 and n 2, which characterizes the strain rate sensitivity and strain hardening effect, respectively can be defined. The saturated displacement is inversely proportional to the stiffening factors (i.e. n 1 and n 2 and saturated impulse is inversely proportional to the square roots of the stiffening factors (i.e. n 1 and n 2. Formulae for displacement and saturated impulse are proposed based on the empirical analysis.

Residual Stress Analysis of Severe Plastic Deformed Materials using the Finite Element Method and the Neutron Diffraction Method

International Nuclear Information System (INIS)

Kang, Mi Hyun; Seong, Back Suck; Kim, Hyoung Seop

2009-01-01

Severe plastic deformation (SPD) is one of the most promising top-down techniques, moving towards industrialization to fabricate bulk ultrafine grain materials. The strain distribution and deformation behavior during the ECAP (equal channel angular pressing), influenced by tool angles, friction and material behavior, was studied through experimental and numerical analyses. The residual stress of work piece which was straight before ECAP produces many serious problems in the next processing e.g. input of the work piece for the next ECAP. The bent work piece needs additional straightening or surface polishing even if the amount of bending is small, and residual stress need to be released before service applications. Residual stress, particularly tensile residual stress can be a very important factor in affecting the reliability and integrity of working parts. The formation of tensile residual stress may result in initiation of fatigue cracks, stress corrosion cracking, or other types of fracture. Hence, residual stress and resulting bending need to be controlled during ECAP. Thus, in current study the bending behavior and the residual stress of the work piece in ECAP are analyzed through experimental and finite element analyses by considering the effects of material, geometric, and processing parameters individually. The stress states in the ECAP processed work piece were measured by the non-destructive way using neutron diffraction. Efforts were made to suggest the alternate routes to reduce the residual stress and bending of work piece in ECAP

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.

Modelling of transport and collisions between rigid bodies to simulate the jam formation in urban flows

Directory of Open Access Journals (Sweden)

S Hadji

2008-09-01

Full Text Available This study deals with the simulation of transport and interaction betweenbodies considered as a rectangular shape particles, in urban flow. We usedan hydrodynamic two-dimensional finite elements model coupled to theparticles model based on Maxey-Riley equations, and taking into accountof contact between bodies. The finite element discretization is based onthe velocity field richer than pressure field, and the particles displacementsare computed by using a rigid body motion method. A collision strategy isalso developed to handle cases in which bodies touch.

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

Extensions to a nonlinear finite-element axisymmetric shell model based on Reissner's shell theory

International Nuclear Information System (INIS)

Cook, W.A.

1981-01-01

Extensions to shell analysis not usually associated with shell theory are described in this paper. These extensions involve thick shells, nonlinear materials, a linear normal stress approximation, and a changing shell thickness. A finite element shell-of-revolution model has been developed to analyze nuclear material shipping containers under severe impact conditions. To establish the limits for this shell model, the basic assumptions used in its development were studied; these are listed in this paper. Several extensions were evident from the study of these limits: a thick shell, a plastic hinge, and a linear normal stress

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.

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.

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.

An analysis of heat field of metal sheet during elastic-plastic deformation

International Nuclear Information System (INIS)

Li, S.X.; Huang, Y.; Shih, C.H.

1985-08-01

This paper describes the application of the finite element analysis to calculate the temperature distribution generated during the process of elastic-plastic deformation. A better agreement is found between the results of heat field computed by use of the finite element analysis and that measured by use of an infrared camera. The results indicate that the method of finite element analysis used for heat field evaluation is reliable. (author)

NON-LINEAR FINITE ELEMENT MODELING OF DEEP DRAWING PROCESS

Directory of Open Access Journals (Sweden)

Hasan YILDIZ

2004-03-01

Full Text Available Deep drawing process is one of the main procedures used in different branches of industry. Finding numerical solutions for determination of the mechanical behaviour of this process will save time and money. In die surfaces, which have complex geometries, it is hard to determine the effects of parameters of sheet metal forming. Some of these parameters are wrinkling, tearing, and determination of the flow of the thin sheet metal in the die and thickness change. However, the most difficult one is determination of material properties during plastic deformation. In this study, the effects of all these parameters are analyzed before producing the dies. The explicit non-linear finite element method is chosen to be used in the analysis. The numerical results obtained for non-linear material and contact models are also compared with the experiments. A good agreement between the numerical and the experimental results is obtained. The results obtained for the models are given in detail.

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

BEAM: A Finite Element Program for the Collapse Analysis of Vehicle Structures

Science.gov (United States)

1994-06-01

deflects a latera: d&stance 8, its bending stresses are increased. Nor can BEAM account for the reduction of plastic moment capacity due to axial loads...Figure 9: The load -displacement curve for Frame 4, comparing elastic-, rigid plastuc and Sttq’ BI-Step analyses with experimental results. The

On the stress assessment and verification of 3D finite element models

International Nuclear Information System (INIS)

Neto, M.M.; Miranda, C.A.J.; Cruz, J.R.B.; Bezerra, L.M.

1995-01-01

The evaluation of components using three dimensional (3D) finite element analysis (FEA) does not generally fall into the shell type verification. Consequently, the demonstration that the modes of failure are avoided sometimes is not straightforward. Elastic rules, developed by limit load theory, require the computation of the shell type through wall membrane and bending stresses. How to calculate these stresses from 3D FEA is not necessarily self-evident. One approach to be considered is to develop recommendations in a case-by-case basis for the most common pressure vessel geometries and loads based on comparison between the results of elastic and also plastic FEA. In this paper the case of a complex geometry -- lugs attached to a cylindrical pressure vessel wall -- is examined and discussed. This case is typically a three-dimensional (3D) configuration where it is not a simple task to check the requirements of the ASME code. From the comparison of the results of 3D elastic and elastic-plastic FEA some conclusions are addressed

A Finite Element Model to Simulate Defect Formation during Friction Stir Welding

Directory of Open Access Journals (Sweden)

Zhi Zhu

2017-07-01

Full Text Available In this study, a 3D coupled thermo-mechanical finite element model is developed to predict and analyze the defect formation during friction stir welding based on coupled Eulerian Lagrangian method. The model is validated by comparing the estimated welding temperature, processed zone shape and void size with those obtained experimentally. The results compared indicate that the simulated temperature and the data measured are in good agreement with each other. In addition, the model can predict the plasticized zone shape and the presence of a void in the weld quite accurately. However, the void size is overestimated. The effects of welding parameters and tool pin profile are also analyzed. The results reveal that welding at low welding speed or high tool rotational speed could produce a smaller void. Moreover, compared to a smooth tool pin, a featured tool pin can enhance plastic flow in the weld and achieve defect-free weldment. The results are helpful for the optimization of the welding process and the design of welding tools.

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.

Stress Concentration Factor of Expanded Aluminum Tubes Using Finite Element Modeling

Directory of Open Access Journals (Sweden)

L Mhamdi

2013-06-01

Full Text Available This paper discusses the development of semi-empirical relations for the maximum stress concentration factor (SCF around circular holes embedded in aluminum tubes under various expansion ratios and mandrel angles. Finite element models were developed to study the expansion of a typical aluminum tube with embedded holes of various sizes. An elastic perfectly-plastic material behaviour was used to describe the structural response of the tubes under expansion. Various hole-diameter-to-tubewall- thickness ratios, tube expansion ratios, and mandrel angles were considered to determine the stress state around the hole at zero and 90 degree locations from which the maximum SCF was determined. Semi-empirical relations for the maximum SCF using the Lagrange interpolation formulation were developed. The developed relations were found to predict the SCFs accurately.

Finite element method programs to analyze irradiation behavior of fuel pellets

International Nuclear Information System (INIS)

Yamada, Rayji; Harayama, Yasuo; Ishibashi, Akihiro; Ono, Masao.

1979-09-01

For the safety assessment of reactor fuel, it is important to grasp local changes of fuel pins due to irradiation in a reactor. Such changes of fuel result mostly from irradiation of fuel pellets. Elasto-plastic analysis programs based on the finite element method were developed to analyze these local changes. In the programs, emphasis is placed on the analysis of cracks in pellets; the interaction between cracked-pellets and cladding is not taken into consideration. The two programs developed are FEMF3 based on a two-dimensional axially symmetric model (r-z system) and FREB4 on a two-dimensional plane model (r-theta system). It is discussed in this report how the occurrence and distribution of cracks depend on heat rate of the fuel pin. (author)

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

DEFF Research Database (Denmark)

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

1999-01-01

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

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

Finite Element Modeling of Material Fatigue and Cracking Problems for Steam Power System HP Devices Exposed to Thermal Shocks

Directory of Open Access Journals (Sweden)

Pawlicki Jakub

2016-09-01

Full Text Available The paper presents a detailed analysis of the material damaging process due to low-cycle fatigue and subsequent crack growth under thermal shocks and high pressure. Finite Element Method (FEM model of a high pressure (HP by-pass valve body and a steam turbine rotor shaft (used in a coal power plant is presented. The main damaging factor in both cases is fatigue due to cycles of rapid temperature changes. The crack initiation, occurring at a relatively low number of load cycles, depends on alternating or alternating-incremental changes in plastic strains. The crack propagation is determined by the classic fracture mechanics, based on finite element models and the most dangerous case of brittle fracture. This example shows the adaptation of the structure to work in the ultimate conditions of high pressure, thermal shocks and cracking.

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

Effect of material variation on the biomechanical behaviour of orthodontic fixed appliances: a finite element analysis.

Science.gov (United States)

Papageorgiou, Spyridon N; Keilig, Ludger; Hasan, Istabrak; Jäger, Andreas; Bourauel, Christoph

2016-06-01

Biomechanical analysis of orthodontic tooth movement is complex, as many different tissues and appliance components are involved. The aim of this finite element study was to assess the relative effect of material alteration of the various components of the orthodontic appliance on the biomechanical behaviour of tooth movement. A three-dimensional finite element solid model was constructed. The model consisted of a canine, a first, and a second premolar, including the surrounding tooth-supporting structures and fixed appliances. The materials of the orthodontic appliances were alternated between: (1) composite resin or resin-modified glass ionomer cement for the adhesive, (2) steel, titanium, ceramic, or plastic for the bracket, and (3) β-titanium or steel for the wire. After vertical activation of the first premolar by 0.5mm in occlusal direction, stress and strain calculations were performed at the periodontal ligament and the orthodontic appliance. The finite element analysis indicated that strains developed at the periodontal ligament were mainly influenced by the orthodontic wire (up to +63 per cent), followed by the bracket (up to +44 per cent) and the adhesive (up to +4 per cent). As far as developed stresses at the orthodontic appliance are concerned, wire material had the greatest influence (up to +155 per cent), followed by bracket material (up to +148 per cent) and adhesive material (up to +8 per cent). The results of this in silico study need to be validated by in vivo studies before they can be extrapolated to clinical practice. According to the results of this finite element study, all components of the orthodontic fixed appliance, including wire, bracket, and adhesive, seem to influence, to some extent, the biomechanics of tooth movement. © The Author 2015. Published by Oxford University Press on behalf of the European Orthodontic Society. All rights reserved. For permissions, please email: journals.permissions@oup.com.

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

Continuum mechanics elasticity, plasticity, viscoelasticity

CERN Document Server

Dill, Ellis H

2006-01-01

FUNDAMENTALS OF CONTINUUM MECHANICSMaterial ModelsClassical Space-TimeMaterial BodiesStrainRate of StrainCurvilinear Coordinate SystemsConservation of MassBalance of MomentumBalance of EnergyConstitutive EquationsThermodynamic DissipationObjectivity: Invariance for Rigid MotionsColeman-Mizel ModelFluid MechanicsProblems for Chapter 1BibliographyNONLINEAR ELASTICITYThermoelasticityMaterial SymmetriesIsotropic MaterialsIncompressible MaterialsConjugate Measures of Stress and StrainSome Symmetry GroupsRate Formulations for Elastic MaterialsEnergy PrinciplesGeometry of Small DeformationsLinear ElasticitySpecial Constitutive Models for Isotropic MaterialsMechanical Restrictions on the Constitutive RelationsProblems for Chapter 2BibliographyLINEAR ELASTICITYBasic EquationsPlane StrainPlane StressProperties of SolutionsPotential EnergySpecial Matrix NotationThe Finite Element Method of SolutionGeneral Equations for an Assembly of ElementsFinite Element Analysis for Large DeformationsProblems for Chapter 3Bibliograph...

Engineering computation of structures the finite element method

CERN Document Server

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

2015-01-01

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

Modelling drawbeads with finite elements and verification

NARCIS (Netherlands)

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

1994-01-01

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

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

Finite element simulation of lower limb injuries to the driver in minibus frontal collisions.

Science.gov (United States)

Shi, Liang-Liang; Lei, Chen; Li, Kui; Fu, Shuo-Zhen; Wu, Zheng-Wei; Yin, Zhi-Yong

2016-06-01

This study aims to explore the biomechanical mechanism of lower limb injuries to the driver by establishing a finite element (FE) simulation model of collisions. First a minibus FE model was integrated with a seat belt system. Then it was used to rebuild two collisions together with the total human model for safety (THUMS) provided by Toyota Motor Corporation: a rear-end collision between a minibus and a truck and a head-on collision of a minibus to a rigid wall. The impact velocities of both collisions were set at 56 km/h. The vehicle dynamic response, vehicle deceleration, and dashboard intrusion in the two collisions were compared. In the minibus rear-end truck collision, the peak values of the von Mises equivalent stress at the tibia and the femur were 133 MPa and 126 MPa respectively; while in the minibus head-on rigid wall collision, the data were 139 MPa and 99 MPa. Compared with the minibus head-on rigid wall collision, the vehicle deceleration was smaller and the dashboard intrusion was larger in the minibus rear-end truck collision. The results illustrate that a longer dashboard incursion distance corresponds to a higher von Mises equivalent stress at the femur. The simulation results are consistent with the driver's autopsy report on lower limbs injuries. These findings verify that FE simulation method is reliable and useful to analyze the mechanisms of lower limb injuries to the driver in minibus frontal collisions.

Effect of Punch Stroke on Deformation During Sheet Forming Through Finite Element

Science.gov (United States)

Akinlabi, Stephen; Akinlabi, Esther

2017-08-01

Forming is one of the traditional methods of making shapes, bends and curvature in metallic components during a fabrication process. Mechanical forming, in particular, employs the use of a punch, which is pressed against the sheet material to be deformed into a die by the application of an external force. This study reports on the finite element analysis of the effects of punch stroke on the resulting sheet deformation, which is directly a function of the structural integrity of the formed components for possible application in the automotive industry. The results show that punch stroke is directly proportional to the resulting bend angle of the formed components. It was further revealed that the developed plastic strain increases as the punch stroke increases.

Numerical investigation of a shot peening process by a finite element approach

DEFF Research Database (Denmark)

Liu, Hongsheng; Zhang, Xiaodan; Hansen, Niels

2014-01-01

Shot peening is a surface impact treatment widely used to improve the performance of a metal or a component. The better performance of the shot peened part is controlled by compressive residual stresses resulting from the plastic deformation of the surface layers by impacts of shot. The compressive...... residual stress is generally measured by X-ray diffraction. However, considerable cost and time are needed for such measurements. For this reason, in this work a 3D finite element (FE) model is introduced for a shot peening process. Through the FE simulations, the effect of process parameters...... such as damping ratio of material, friction coefficient, shot velocity and shot angle on the magnitude and distribution of the compressive residual stress is examined....

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

Elastic-plastic failure analysis of pressure burst tests of thin toroidal shells

International Nuclear Information System (INIS)

Jones, D.P.; Holliday, J.E.; Larson, L.D.

1998-07-01

This paper provides a comparison between test and analysis results for bursting of thin toroidal shells. Testing was done by pressurizing two toroidal shells until failure by bursting. An analytical criterion for bursting is developed based on good agreement between structural instability predicted by large strain-large displacement elastic-plastic finite element analysis and observed burst pressure obtained from test. The failures were characterized by loss of local stability of the membrane section of the shells consistent with the predictions from the finite element analysis. Good agreement between measured and predicted burst pressure suggests that incipient structural instability as calculated by an elastic-plastic finite element analysis is a reasonable way to calculate the bursting pressure of thin membrane structures

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