Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element problems”

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.

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Vibration of thick composite laminates - Analytic theory and finite element approximations

Two versions of a (1,2) higher-order laminate plate theory are extended to elastodynamics, and their performance on free vibration problems for homogeneous and laminated plates is examined. Both versions are based on continuously varying through-the-thickness linear and parabolic distributions for the inplane and transverse displacements, respectively. They also use continuous parabolic distributions for the transverse shear strains. The theory is capable of providing an improved accuracy as compared to other approximate theories and accurately predicting the lowest thickness-stretch frequencies which cannot be obtained with shear-deformable theories. The theory is considered to be particularly well-suited for generating simple and efficient variationally-based displacement finite elements leading to C sup 0 elements with only five degrees of freedom per node.

Tessler, A.↗
DOE OSTI2024

Prediction of residual stresses in additively manufactured parts using lumped capacitance and classical lamination theory

Several industries are interested in Laser Powder Bed Fusion (L-PBF) Additively Manufactured (AM) metal parts because their designs can be made arbitrarily complex while retaining bulk-type material properties. However, the residual stresses (RS) and distortions caused by the heat gradients inherent to L-PBF processes are detrimental to the structural integrity of the parts and must be taken into consideration during the part design cycle. Predicting the state of stresses in as-built 3D printed parts is a difficult problem that is typically approached with the use of transient thermomechanical Finite Element Models (FEMs). However, the nonlinearities associated with AM processes are difficult to capture in these FEMs without increasing the computational cost of the simulation, limiting their ability to be incorporated into practical design cycles. This work presents a novel analytical framework that combines lumped capacitance nonlinear heat transfer with time dependent classical lamination theory to efficiently and accurately predict RS in as-built L-PBF parts without the need of FEMs. The simulation was compared to Neutron Diffraction (ND) residual strain measurements taken at Oak Ridge National Laboratories (ORNL) as well as Synchrotron X-ray Diffraction (XRD) strain data published by the National Institute of Standards and Technology (NIST). The simulation predictions and the experimental data showed excellent agreement for the in-plane strain directions, and general agreement for the out of plane strain component, highlighting an area where further development can be implemented.

42 ENGINEERING↗
DOE OSTI2020

A particle finite element-based model for droplet spreading analysis

Here, a particle finite element method-based model is proposed to analyze droplet dynamics problems, particularly droplet spreading on solid substrates (wetting). The model uses an updated Lagrangian framework to formulate the governing equations of the liquid. The curvature of the liquid surface is tracked accurately using a deforming boundary mesh. In order to predict the spreading rate of the droplet on the solid substrate and track the corresponding contact angle evolution, dissipative forces at the contact line are included in the formulation in addition to the Navier-slip boundary conditions at the solid–liquid interface. The inclusion of these boundary conditions makes it possible to account for the induced Young’s stress at the contact line and for the viscous dissipation along the solid–liquid interfacial region. These are found to be essential to obtain a mesh-independent physical solution. The temporal evolution of the contact angle and the contact line velocity of the proposed model are compared with spreading droplets and micro-sessile droplet injection experiments and are shown to be in good agreement.

42 ENGINEERING↗
DOE OSTI2022

Three-dimensional modelling of a self-sustained atmospheric pressure glow discharge

The atmospheric pressure glow discharge (APGD) is a relatively simple and versatile plasma source used in diverse applications. Stable APGD operation at high currents, generally a challenge due to instabilities leading to glow-to-arc transition, has been demonstrated using actively-controlled cathodic cooling. This article presents the computational modelling and simulation of a self-sustained direct-current APGD in helium within a 10 mm pin-to-plate inter-electrode gap for currents ranging from 4 to 40 mA. The APGD model is comprised of the conservation equations for total mass, chemical species, momentum, thermal energy of heavy-species and of free electrons, and electric charge. The model equations are discretized using a nonlinear variational multi-scale finite element method that has demonstrated superior accuracy in other plasma flow problems, on a temporal and three-dimensional computational domain suitable to unveil the potential occurrence of instabilities. Modelling results show good agreement with experimental measurements of voltage drop and the same trend but higher values of temperature. The higher temperatures obtained by the simulations appear to be due to the absence of a near-cathode heat dissipation model. Here, the results also reveal that the distribution of electron density and of the ratio of atomic helium ions to total ions transitions from monotonically increasing away from the cathode to presenting a minimum near the centre of the gap with increasing current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Langevin and Fokker-Planck analyses for diffusion-mediated passing of circular and discorectangular species in two-dimensional channels

The propensity for pairs of diffusing species to pass each other within narrow channels or pores is of basic interest as a first-passage-type problem. It is also of relevance for solution-phase transport in nanoporous materials, and in particular for catalytic conversion reactions where high yield requires that product species can efficiently pass reactant species to exit the pores. Here, we analyze a two-dimensional model with nonoverlapping circular and discorectangular species confined to a rectangular channel, and where passing is mediated by Brownian dynamics in an implicit solvent. For narrower channels where passing is still possible, the discorectangle must align with the channel to pass the circular species. Behavior of the passing propensity, 𝑃, can be assessed by strongly damped Langevin simulations, or within an equivalent Fokker-Planck equation (FPE) formalism. The latter corresponds to a diffusion problem in a “higher-dimensional channel” with a constriction. We assess the variation of the passing propensity, 𝑃, for a broad range of channel width including its scaling just above the threshold where passing is sterically blocked. Analysis of 𝑃 versus the rotational diffusion coefficient 𝐷 𝑟 of the discorectangle reveals a significant decrease in 𝑃 for lower 𝐷 𝑟 for moderate channel width. This prompts a direct analysis of the regime where 𝐷 𝑟 → 0 , for which the FPE can be reduced to a three-dimensional diffusion problem, precise analysis of which is facilitated by adaptive-mesh finite element methods. The dependence of 𝑃 on the aspect ratio of the discorectangle is also assessed.

Rahman, Md Khaledur [Ames Laboratory (AMES), Ames,↗
DOE OSTI2023

An $\bar{F}$ Meshfree Treatment for Nearly Incompressible Materials

The computational modeling of nearly incompressible materials is a difficult task for many numerical methods, and even after several decades of investigation, it is still an active research area. This report seeks to address the treatment of incompressible materials in meshfree methods using a synergistic combination of two treatments. The first treatment is an $\bar{F}$ method, where the decomposed dilatational and deviatoric parts are calculated over different smoothing domains. The second treatment “activates” additional nodes throughout the domain to increase the flexibility of the model. We implement this synergistic combination in the context of the reproducing kernel particle method (RKPM) and present results for the Cook’s membrane benchmark problem. The results are compared with those using the composite tet10 finite element with a volume-averaged J formulation. We show that the combined treatment is an effective way to deal with nearly incompressible materials in a meshfree framework and compares well with other highly-effective treatments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Stability of Panels of a Thermoelastic Shell Heated at the Edge by a Mobile Point Source

The thermostability of thin-walled structures is considered by means of a thermoelasticity problem in a three-dimensional quasi-static formulation. Finite-element models of panels in the form of a circular cylinder and a plate reinforced by ribs are constructed. Numerical solutions are obtained for the stability loss of thermoelastic thin-walled reinforced structures with intense edge heating by a mobile heat source.

42 ENGINEERING↗
NASA NTRS1972

Structural-Thermal-Optical Program (STOP)

A structural thermal optical computer program is developed which uses a finite element approach and applies the Ritz method for solving heat transfer problems. Temperatures are represented at the vertices of each element and the displacements which yield deformations at any point of the heated surface are interpolated through grid points.

Lee, H. P.↗
NASA NTRS1971

The application of a curved, mixed-type shell element.

The variational condition given by Herrmann (1967) for finite element bending analyses of plates is extended to general thin shell problems. Emphasis is laid upon shell behavior and, in particular, on boundary layers and inextensional bending. A study of the boundary layer problem on a cylindrical shell compares the refined element and the Herrmann-Campbell (1968) element. Some numerical results on cross-tube and mitred-bend problems are given.

Visser, W.↗
NASA NTRS1976

Finite element solutions for crack-tip behavior in small-scale yielding

The subject considered is the stress and deformation fields in a cracked elastic-plastic power law hardening material under plane strain tensile loading. An incremental plasticity finite element formulation is developed for accurate analysis of the complete field problem including the extensively deformed near tip region, the elastic-plastic region, and the remote elastic region. The formulation has general applicability and was used to solve the small scale yielding problem for a set of material hardening exponents. Distributions of stress, strain, and crack opening displacement at the crack tip and through the elastic-plastic zone are presented as a function of the elastic stress intensity factor and material properties.

Tracey, D. M.↗
NASA NTRS1976

The Boeing plastic analysis capability for engines

The current BOPACE program is described as a nonlinear stress analysis program, which is based on a family of isoparametric finite elements. The theoretical, user, programmer, preprocessing aspects are discussed, and example problems are included. New features in the current program version include substructuring, an out-of-core Gauss wavefront equation solver, multipoint constraints, combined material and geometric nonlinearities, automatic calculation of inertia effects, provision for distributed as well as concentrated mechanical loads, follower forces, singular crack-tip elements, the SAIL automatic generation capability, and expanded user control over input quantity definition, output selection, and program execution. BOPACE is written in FORTRAN 4 and is currently available for both the IBM 360/370 and the UNIVAC 1108 machines.

Vos, R. G.↗
NASA NTRS1977

Numerical analysis of free vibrations of damped rotating structures

This paper is concerned with the efficient numerical solution of damped and undamped free vibration problems of rotating structures. While structural discretization is achieved by the finite element method, the associated eigenproblem solution is effected by a combined Sturm sequence and inverse iteration technique that enables the computation of a few required roots only without having to compute any other. For structures of complex configurations, a modal synthesis technique is also presented, which is based on appropriate combinations of eigenproblem solution of various structural components. Such numerical procedures are general in nature, which fully exploit matrix sparsity inherent in finite element discretizations, and prove to be most efficient for the vibration analysis of any damped rotating structure, such as rotating machineries, helicopter and turbine blades, spinning space stations, among others.

Gupta, K. K.↗
NASA NTRS1981

Application of 'steady' state finite element and transient finite difference theory to sound propagation in a variable duct - A comparison with experiment

Experimental data are presented for sound propagation in a simulated infinite hard wall duct with a large change in duct cross sectional area. The data are conveniently tabulated for further use. The 'steady' state finite element theory of Astley and Eversman (1981) and the transient finite difference theory of White (1981) are in good agreement with the data for both the axial and transverse pressure profiles and the axial phase angle. Therefore, numerical finite difference and finite element theories appear to be ideally suited for handling duct propagation problems which encounter large axial gradients in acoustic parameters. The measured energy reflection coefficient agrees with the values from the Astley-Eversman modal coupling model.

Baumeister, K. J.↗
NASA NTRS1983

Modeling global structural damping in trusses using simple continuum models

Truss beams with members having viscous damping are modeled as continuum Timoshenko beams. Procedures for deriving the equivalent beam stiffnesses and damping are presented. The global damping for the continuum beam is explicitly expressed in terms of the damping coefficients of the individual truss members. The continuum beam model is used to study transient vibration problems and the solutions are compared well with the full scale finite element solutions. The gradient method is used for parameter estimations in conjunction with the Timoshenko beam model. It is shown that the Timoshenko beam model can be easily updated with measured data and the updated model can yield very accurate transient solutions.

Sun, C. T.↗
NASA NTRS1984

A shock-capturing finite element method

Finite element methods are developed for a one-dimensional singular perturbation problem involving a nonlinear flux function. In the limit of vanishing dissipation, the numerical flux is found to be almost identical to that due to Engquist and Osher, the only difference being in the treatment of transonic compression. This facilitates sharper resolution of stationary shocks, but is accompanied by loss of smoothness and monotonicity.

Hughes, T. J. R.↗
NASA NTRS1989

Static aeroelastic analysis and tailoring of missile control fins

A concept for enhancing the design of control fins for supersonic tactical missiles is described. The concept makes use of aeroelastic tailoring to create fin designs (for given planforms) that limit the variations in hinge moments that can occur during maneuvers involving high load factors and high angles of attack. It combines supersonic nonlinear aerodynamic load calculations with finite-element structural modeling, static and dynamic structural analysis, and optimization. The problem definition is illustrated. The fin is at least partly made up of a composite material. The layup is fixed, and the orientations of the material principal axes are allowed to vary; these are the design variables. The objective is the magnitude of the difference between the chordwise location of the center of pressure and its desired location, calculated for a given flight condition. Three types of constraints can be imposed: upper bounds on static displacements for a given set of load conditions, lower bounds on specified natural frequencies, and upper bounds on the critical flutter damping parameter at a given set of flight speeds and altitudes. The idea is to seek designs that reduce variations in hinge moments that would otherwise occur. The block diagram describes the operation of the computer program that accomplishes these tasks. There is an option for a single analysis in addition to the optimization.

Mcintosh, S. C., Jr.↗
NASA NTRS1988

Probabilistic structural analysis methods and applications

An advanced algorithm for simulating the probabilistic distribution of structural responses due to statistical uncertainties in loads, geometry, material properties, and boundary conditions is reported. The method effectively combines an advanced algorithm for calculating probability levels for multivariate problems (fast probability integration) together with a general-purpose finite-element code for stress, vibration, and buckling analysis. Application is made to a space propulsion system turbine blade for which the geometry and material properties are treated as random variables.

Cruse, T. A.↗