Methods of mechanical contact constraint enforcement using the mortar finite element method
Explore the source record for details and available documents.
Engineering - Mechanical and civil engineering↗
SEARCH · Engineering Papers
Results for “Finite Element Method”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Developing Procedures to Implement Geometric Imperfections Beyond Right Circular Cylindrical Shells in Finite Element Method Models
Analysis of aerospace structures is frequently conducted using nominal dimensions and frequently assumes ideal conditions in loading, contact, constraints, et cetera. Off-nominal dimensions and nonideal conditions, however, are present in all structures. These are the result of widely ranging causes from coefficient of thermal expansion mismatches, manufacturing tooling anomalies, to assembly procedures that inadvertently alter the structure. Specifically, geometric imperfections can have potentially significant influence on the response of a structural test article observed in an experiment versus the response given by a numerical simulation. The Python Tool for Implementing Geometric Imperfections in Reduced Structures (Py_TIGIRS) was previously presented as a set of Python scripts to calculate and implement as-manufactured geometric midsurface and thickness imperfections into finite element method (FEM) shell models of nominally right circular cylinders. By taking advantage of the simple shape of a right circular cylinder, interpolations of the measured data points were able to be performed along directions that aligned to the cylindrical coordinate system axes of the entire structure. By taking advantage of the shell representation of the real structure as opposed to modeling using a continuum representation, the thickness variation was able to be implemented by shell section definitions instead of having to modify the position of multiple nodes in the thickness direction. Py_TIGIRS is a useful tool that established a procedural example on how to implement geometric imperfections in right circular cylindrical shell structures. Three new procedures, each expanded from concepts established in Py_TIGIRS, are proposed for various test-article designs and are intended to broaden the range of structures that can be modeled with measured geometric imperfections in the structural analysis community. Each test-article design introduces new challenges to successfully implement geometric imperfections into a FEM model. The first test-article design consists of a carbon fiber reinforced polymer square plate with a hat-shaped stiffener co-cured on one side. This test-article design was for a novel seven-point bend test that was also previously presented. Manufacturing and cure-cycle imperfections are observed using digital image correlation (DIC) techniques. As thermal expansion coefficient mismatches between the plate and stiffener materials were anticipated, a thermal analysis study with continuum shell and solid elements was conducted to capture the global shape observed prior to testing. The second test-article design is of a similar hat-stiffened plate configuration, but with a side length ratio near 3:1 with elongation in the stiffener direction. The test article was used to characterize the response to uniaxial compressive loading in the direction of the stiffener. Due to differing manufacturing steps, a thermal analysis like the one developed for the seven-point bend configuration was unable to mimic the observed geometric imperfections. Instead, a strategy based on applying deformations directly to the structure during analysis was developed for continuum shell and solid element representation of a stiffened panel.
Geometric imperfections↗
Analysis and Development of Finite Element Methods for the Study of Nonlinear Thermomechanical Behavior of Structural Components
Underintegrated methods are investigated with respect to their stability and convergence properties. The focus was on identifying regions where they work and regions where techniques such as hourglass viscosity and hourglass control can be used. Results obtained show that underintegrated methods typically lead to finite element stiffness with spurious modes in the solution. However, problems exist (scalar elliptic boundary value problems) where underintegrated with hourglass control yield convergent solutions. Also, stress averaging in underintegrated stiffness calculations does not necessarily lead to stable or convergent stress states.
Oden, J. Tinsley↗
Finite Element Method for Thermal Analysis
A two- and three-dimensional, finite-element thermal-analysis program which handles conduction with internal heat generation, convection, radiation, specified flux, and specified temperature boundary conditions is presented. Elements used in the program are the triangle and tetrahedron for two- and three-dimensional analysis, respectively. The theory used in the program is developed, and several sample problems demonstrating the capability and reliability of the program are presented. A guide to using the program, description of the input cards, and program listing are included.
Heuser, J.↗
A finite element method for nonlinear forced vibrations of beams
Techniques for defining a finite element model (FEM) for analysis of nonlinear vibrations in beam structures subjected to harmonic excitation are presented. The resulting model covers longitudinal deformation and inertial effects. The nonlinear oscillations of a beam element under forced excitation are modeled by a harmonic force matrix based on first order approximations of the Jacobian elliptic forcing function. Harmonic force and nonlinear stiffness matrices are derived and the nonlinear forced responses of beams are calculated under various boundary conditions. The results of FEM computations for simply-supported and clamped beams show that midplane stretching caused by large deflections increases the nonlinearity. Axially-restrained beams experience only hardening nonlinearity, while axially-free beams have reduced nonlinearity in deformation and inertia and an increase in linearity due to large deflection.
Mei, C.↗
Simulation of Two-Fluid Flows by the Least-Squares Finite Element Method Using a Continuum Surface Tension Model
In this paper a numerical procedure for simulating two-fluid flows is presented. This procedure is based on the Volume of Fluid (VOF) method proposed by Hirt and Nichols and the continuum surface force (CSF) model developed by Brackbill, et al. In the VOF method fluids of different properties are identified through the use of a continuous field variable (color function). The color function assigns a unique constant (color) to each fluid. The interfaces between different fluids are distinct due to sharp gradients of the color function. The evolution of the interfaces is captured by solving the convective equation of the color function. The CSF model is used as a means to treat surface tension effect at the interfaces. Here a modified version of the CSF model, proposed by Jacqmin, is used to calculate the tension force. In the modified version, the force term is obtained by calculating the divergence of a stress tensor defined by the gradient of the color function. In its analytical form, this stress formulation is equivalent to the original CSF model. Numerically, however, the use of the stress formulation has some advantages over the original CSF model, as it bypasses the difficulty in approximating the curvatures of the interfaces. The least-squares finite element method (LSFEM) is used to discretize the governing equation systems. The LSFEM has proven to be effective in solving incompressible Navier-Stokes equations and pure convection equations, making it an ideal candidate for the present applications. The LSFEM handles all the equations in a unified manner without any additional special treatment such as upwinding or artificial dissipation. Various bench mark tests have been carried out for both two dimensional planar and axisymmetric flows, including a dam breaking, oscillating and stationary bubbles and a conical liquid sheet in a pressure swirl atomizer.
Wu, Jie↗
High-order finite element method for atomic structure calculations
We introduce featom, an open source code that implements a high-order finite element solver for the radial Schrödinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or exponential, and the convergence can be systematically controlled by increasing the number and/or polynomial order of the finite element basis functions. The Dirac equation is solved using a squared Hamiltonian approach to eliminate spurious states. Here, to address the slow convergence of the $κ=±1$ states due to divergent derivatives at the origin, we incorporate known asymptotic forms into the solutions. We achieve a high level of accuracy (10 -8 Hartree) for total energies and eigenvalues of heavy atoms such as uranium in both Schrödinger and Dirac Kohn-Sham solutions. We provide detailed convergence studies and computational parameters required to attain commonly required accuracies. Finally, we compare our results with known analytic results as well as the results of other methods. In particular, we calculate benchmark results for atomic numbers (Z) from 1 to 92, verifying current benchmarks. We demonstrate significant speedup compared to the state-of-the-art shooting solver dftatom. An efficient, modular Fortran 2008 implementation, is provided under an open source, permissive license, including examples and tests, wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Real time time-dependent density functional theory using higher order finite-element methods
We present a computationally efficient approach to solve the time-dependent Kohn-Sham equations in real time using higher order finite-element spatial discretization, applicable to both pseudopotential and all-electron calculations. To this end, we develop an a priori mesh adaption technique, based on the semidiscrete (discrete in space but continuous in time) error estimate on the time-dependent Kohn-Sham orbitals, to construct an efficient finite-element discretization. Subsequently, we obtain the full-discrete error estimate to guide our choice of the time step. We employ spectral finite elements along with special reduced order quadrature to render the overlap matrix diagonal, thereby simplifying the inversion of the overlap matrix that features in the evaluation of the discrete time-evolution operator. We use the second-order Magnus operator as the time-evolution operator, wherein the action of the discrete Magnus operator, expressed as exponential of a matrix, on the Kohn-Sham orbitals is obtained efficiently through an adaptive Lanczos iteration. We observe close to optimal rates of convergence of the dipole moment with respect to spatial and temporal discretization, for both pseudopotential and all-electron calculations. We demonstrate a staggering 100-fold reduction in the computational time afforded by higher order finite elements over linear finite elements, for both pseudopotential and all-electron calculations. Further, for similar level of accuracy, we obtain significant computational savings by our approach as compared to state-of-the-art finite-difference methods. Here, we also demonstrate the competence of higher order finite elements for all-electron benchmark systems. Lastly, we observe good parallel scalability of the proposed method on many hundreds of processors.
36 MATERIALS SCIENCE↗
Multi-material ALE remap with interface sharpening using high-order matrix-free finite element methods
The arbitrary Lagrangian-Eulerian (ALE) technique involves remapping field quantities from a Lagrangian mesh to an optimized mesh in a conservative, accurate and bounds-preserving manner. For methods based on arbitrary order finite elements, as described in a reference, material volume fractions are advected in pseudo-time using flux-corrected transport (FCT) without any form of interface reconstruction. In practice, this can lead to excessive propagation of small volume fractions throughout the domain. In addition, this method requires assembly of a global advection matrix to compute the bounds-preserving low-order FCT solution. In this work, we introduce a new approach for ALE remap using a high-order matrix-free technique which incorporates a flux modification to sharpen material interfaces in a conservative manner. Our approach begins with computing a bounds-preserving low-order solution to the ALE remap equations at the element level. We then compute a sharp interface solution (not guaranteed to be bounds-preserving) which comes from solving an augmented version of the ALE remap equations with a conservative flux modification which acts to sharpen material volume fractions based on their gradients and transport directions. Using the sharp interface solution, we make global corrections to the bounds-preserving solution while maintaining preservation of bounds. By blending with the sharpened solution at the global level we are able to globally conserve mass without hindering the remap pseudo-time step. This new interface-aware ALE remap method is based entirely on partial assembly techniques where globally assembled matrix operators are no longer needed, resulting in a globally matrix-free FCT method for multi-material, multi-field ALE remap with high performance on GPU architectures. We present results of our new remap method on 1D, 2D and 3D benchmarks and describe the algorithmic tailoring for GPU architectures that was developed.
Vargas, Arturo [Lawrence Livermore National Labora↗
Verification of a 2-D to 3-D global/local finite element method for symmetric laminates
A two-dimensional to three-dimensional global/local finite element analysis technique, as applied to both cross-ply and general symmetric composite laminates, is presented. In particular, the local areas of interest for the laminate are the straight free edge in proximity to a hole and curved free edge of a hole. Verification of this technique was accomplished by comparison to results from complete three-dimensional finite element analyses. Examination of the interlaminar stress fields and the amount of time for their determination indicated that the global/local method developed yielded accurate results and was more efficient.
Thompson, Danniella M.↗
A simple finite element method for linear hyperbolic problems
Not Available
97 MATHEMATICS AND COMPUTING↗
Finite element methods of analysis for 3D inviscid compressible flows
The applicants have developed a finite element based approach for the solution of three-dimensional compressible flows. The procedure enables flow solutions to be obtained on tetrahedral discretizations of computational domains of complex form. A further development was the incorporation of a solution adaptive mesh strategy in which the adaptivity is achieved by complete remeshing of the solution domain. During the previous year, the applicants were working with the Advanced Aerodynamics Concepts Branch at NASA Ames Research Center with an implementation of the basic meshing and solution procedure. The objective of the work to be performed over this twelve month period was the transfer of the adaptive mesh technology and also the undertaking of basic research into alternative flow algorithms for the Euler equations on unstructured meshes.
Peraire, Jaime↗
Interior energy error estimates for the weak Galerkin finite element method
Consider the Poisson equation in a polytopal domain ΩRC d (d=2,3) as the model problem. In this paper, we study interior energy error estimates for the weak Galerkin finite element approximation to elliptic boundary value problems. In particular, we show that the interior error in the energy norm is bounded by three components: the best local approximation error, the error in negative norms, and the trace error on the element boundaries. This implies that the interior convergence rate can be polluted by solution singularities on the domain boundary, even when the solution is smooth in the interior region. Numerical results are reported to support the theoretical findings. Lastly, to the best of our knowledge, this is the first local energy error analysis that applies to general meshes consisting of polytopal elements and hanging nodes.
97 MATHEMATICS AND COMPUTING↗
Modelling NDE pulse-echo inspection of misorientated planar rough defects using an elastic finite element method
Pulse-echo ultrasonic NDE examination of large pressure vessel forgings is a design and construction code requirement in the power generation industry. Such inspections aim to size and characterise potential defects that may have formed during the forging process. Typically these defects have a range of orientations and surface roughnesses which can greatly affect ultrasonic wave scattering behaviour. Ultrasonic modelling techniques can provide insight into defect response and therefore aid in characterisation. However, analytical approaches to solving these scattering problems can become inaccurate, especially when applied to increasingly complex defect geometries. To overcome these limitations a elastic Finite Element (FE) method has been developed to simulate pulse-echo inspections of embedded planar defects. The FE model comprises a significantly reduced spatial domain allowing for a Monte-Carlo based approach to consider multiple realisations of defect orientation and surface roughness. The results confirm that defects aligned perpendicular to the path of beam propagation attenuate ultrasonic signals according to the level of surface roughness. However, for defects orientated away from this plane, surface roughness can increase the magnitude of the scattered component propagating back along the path of the incident beam. This study therefore highlights instances where defect roughness increases the magnitude of ultrasonic scattered signals, as opposed to attenuation which is more often assumed.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Application of the finite element method to rotary-wing aeroelasticity
Recent research in rotary-wing aeroelasticity has indicated that all fundamental problems in this area are inherently nonlinear. The non-linearities in this problem are due to the inclusion of finite slopes, due to moderate deflections, in the structural, inertia and aerodynamic operators associated with this aeroelastic problem. In this paper the equations of motion, which are both time and space dependent, for the aeroelastic problem are first formulated in P.D.E. form. Next the equations are linearized about a suitable equilibrium position. The spatial dependence in these equations is discretized using a local Galerkin method of weighted residuals resulting in a finite element formulation of the aeroelastic problem. As an illustration the method is applied to the coupled flap-lag problem of a helicopter rotor blade in hover. Comparison of the solutions with previously published solutions establishes the convergence properties of the method. It is concluded that this formulation is a practical tool for solving rotary-wing aeroelastic stability or response problems.
Friedmann, P.↗
Study on prestressed concrete reactor vessel structures. II-5: Crack analysis by three dimensional finite elements method of 1/20 multicavity type PCRV subjected to internal pressure
A three-dimensional finite elements analysis is reported of the nonlinear behavior of PCRV subjected to internal pressure by comparing calculated results with test results. As the first stage, an analysis considering the nonlinearity of cracking in concrete was attempted. As a result, it is found possible to make an analysis up to three times the design pressure (50 kg/sqcm), and calculated results agree well with test results.
Source record↗
Finite element methods of analysis for high speed viscous flows
Over the past three years a finite element based procedure for the solution of high speed viscous compressible flows was developed. The approach followed was to compute steady state solutions via a false transient, using an explicit time stepping scheme, and to attempt to improve the solution quality by incorporating adaptive mesh procedures. The main thrust of the work was to continue on the extension of the approach to the solution of some realistic compressible viscous flows. When flows at high Reynolds number are investigated, is soon becomes apparent that explicit techniwues have to be supplemented if they are to deal effectively with the large variations in element size and aspect ratio which characterize the computational grids necessary for adequate resolutin of the primary flow features. For this reason, the Taylor-Galerkin solution algorithm was rewritten in an explicit/implicit form. Solutions were computed for the problems of a flow past a flat plate, M-3, Re-100; shock/boundary layer interaction, M-2, Re-296000; flow over a compression corner, M-11.68, Re246000; and unifrom flow past a circular cylinder, M-6.34, Re-39770. A summary of the results is included and demonstrates the numerical performance of the scheme.
Source record↗
Application of the p-version of the finite-element method to global-local problems
A brief survey is given of some recent developments in finite-element analysis technology which bear upon the three main research areas under consideration in this workshop: (1) analysis methods; (2) software testing and quality assurance; and (3) parallel processing. The variational principle incorporated in a finite-element computer program, together with a particular set of input data, determines the exact solution corresponding to that input data. Most finite-element analysis computer programs are based on the principle of virtual work. In the following, researchers consider only programs based on the principle of virtual work and denote the exact displacement vector field corresponding to some specific set of input data by vector u(EX). The exact solution vector u(EX) is independent of the design of the mesh or the choice of elements. Except for very simple problems, or specially constructed test problems, vector u(EX) is not known. Researchers perform a finite-element analysis (or any other numerical analysis) because they wish to make conclusions concerning the response of a physical system to certain imposed conditions, as if vector u(EX) were known.
Szabo, Barna A.↗