SEARCH · Engineering Papers

Results for “Finite Elements”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS2023

A Finite Element Approach for Simplified 2D Nonlinear Dynamic Contact/Impact Analysis

In this paper, a simplified numerical approach for finite element dynamic analysis of an inelastic solid structure subjected to solid object impact is presented. The approach approximates the impacting solid as the selected multiple nodes, for which mass of the impactor is distributed. The node-to-segment contact formulation with the penalty constraint technique incorporated is employed to impose contact conditions between the nodes and the surface of the receiver structure. The node-to-segment algorithm is integrated into Newton–Raphson time integration scheme and the Lagrange multiplier technique is applied to enforce the identical displacements for the selected nodes throughout the analysis process. The approach is verified using two-dimensional plane strain models considering elastic-perfectly-plastic material behavior. The results obtained using the proposed approach are in a good agreement with those simulated using a commercial finite element code, ABAQUS dynamic/implicit, in terms of displacements and stress distribution fields. The proposed approach is shown to be computationally superior to general finite element method-based contact/impact analysis without significantly sacrificing the accuracy.

Seungwook Seok↗

NASA NTRS1988

Global/local finite element analysis of composite materials

The motivation for performing global/local finite element analysis in composite materials is described. An example of such an analysis of a composite plate with a central circular hole is presented. Deformed finite element grids and interlaminar normal stress distributions are presented to aid understanding of the plate response. Such distribution at the plate edge is shown to be basically unaffected, although transverse displacements of the edge were slightly different from an analysis of a similar plate with no hole.

Griffin, O. Hayden, Jr.↗

NASA NTRS1992

An assessment of four-noded plate finite elements based on a generalized third-order theory

Plate finite elements based on the generalized third-order theory of Reddy and the first-order shear deformation theory are analyzed and compared on the basis of thick and thin plate modeling behavior, distortion sensitivity, overall accuracy, reliability, and efficiency. In particular, several four-noded Reddy-type elements and the nine-noded Lagrangian and heterosis (Mindlin-type) plate elements are analyzed to assess their behavior in bending, vibration, and stability of isotropic and laminated composite plates. A four-noded Reddy-type element is identified which is free of all spurious stiffness and zero energy modes, computationally efficient, and suitable for use in any general-purpose finite element program.

Averill, R. C.↗

NASA NTRS1991

Analysis of superconducting electromagnetic finite elements based on a magnetic vector potential variational principle

Electromagnetic finite elements are extended based on a variational principle that uses the electromagnetic four potential as primary variable. The variational principle is extended to include the ability to predict a nonlinear current distribution within a conductor. The extension of this theory is first done on a normal conductor and tested on two different problems. In both problems, the geometry remains the same, but the material properties are different. The geometry is that of a 1-D infinite wire. The first problem is merely a linear control case used to validate the new theory. The second problem is made up of linear conductors with varying conductivities. Both problems perform well and predict current densities that are accurate to within a few ten thousandths of a percent of the exact values. The fourth potential is then removed, leaving only the magnetic vector potential, and the variational principle is further extended to predict magnetic potentials, magnetic fields, the number of charge carriers, and the current densities within a superconductor. The new element produces good results for the mean magnetic field, the vector potential, and the number of superconducting charge carriers despite a relatively high system condition number. The element did not perform well in predicting the current density. Numerical problems inherent to this formulation are explored and possible remedies to produce better current predicting finite elements are presented.

Schuler, James J.↗

NASA NTRS1995

Alaska Crustal Deformation: Finite Element Modeling Constrained by Geologic and Very Long Baseline Interferometry Data

We compute crustal motions in Alaska by calculating the finite element solution for an elastic spherical shell problem. The method we use allows the finite element mesh to include faults and very long baseline interferometry (VLBI) baseline rates of change. Boundary conditions include Pacific-North American (PA-NA) plate motions. The solution is constrained by the oblique orientation of the Fairweather-Queen Charlotte strike-slip faults relative to the PA-NA relative motion direction and the oblique orientation from normal convergence of the eastern Aleutian trench fault systems, as well as strike-shp motion along the Denali and Totschunda fault systems. We explore the effects that a range of fault slip constraints and weighting of VLBI rates of change has on the solution. This allows us to test the motion on faults, such as the Denali fault, where there are conflicting reports on its present-day slip rate. We find a pattern of displacements which produce fault motions generally consistent with geologic observations. The motion of the continuum has the general pattern of radial movement of crust to the NE away from the Fairweather-Queen Charlotte fault systems in SE Alaska and Canada. This pattern of crustal motion is absorbed across the Mackenzie Mountains in NW Canada, with strike-slip motion constrained along the Denali and Tintina fault systems. In south central Alaska and the Alaska forearc oblique convergence at the eastern Aleutian trench and the strike-shp motion of the Denali fault system produce a counterclockwise pattern of motion which is partially absorbed along the Contact and related fault systems in southern Alaska and is partially extruded into the Bering Sea and into the forearc parallel the Aleutian trench from the Alaska Peninsula westward. Rates of motion and fault slip are small in western and northern Alaska, but the motions we compute are consistent with the senses of strike-slip motion inferred geologically along the Kaltag, Kobuk Trench, and Thompson Creek faults and with the normal faulting observed in NW Alaska near Nome. The nonrigid behavior of our finite element solution produces patterns of motion that would not have been expected from rigid block models: strike-slip faults can exist in a continuum that has motion mostly perpendicular to their strikes, and faults can exhibit along-strike differences in magnitudes and directions.

Lundgren, Paul↗

NASA NTRS1989

Thermal-structural finite element analysis using linear flux formulation

A linear flux approach is developed for a finite element thermal-structural analysis of steady-state thermal and structural problems. The element fluxes are assumed to vary linearly in the same form as the element unknown variables, and the finite element matrices are evaluated in closed form. Since numerical integration is avoided, significant computational time saving is achieved. Solution accuracy and computational speed improvements are demonstrated by solving several two-and three-dimensional thermal-structural examples.

Pandey, Ajay K.↗

NASA NTRS1990

Thermal-structural finite element analysis using linear flux formulation

A linear flux approach is developed for a finite element thermal-structural analysis of steady state thermal and structural problems. The element fluxes are assumed to vary linearly in the same form as the element unknown variables, and the finite element matrices are evaluated in closed form. Since numerical integration is avoided, significant computational time saving is achieved. Solution accuracy and computational speed improvements are demonstrated by solving several two and three dimensional thermal-structural examples.

Pandey, Ajay K.↗

NASA NTRS1991

A finite element conjugate gradient FFT method for scattering

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.

Collins, Jeffery D.↗

NASA NTRS1989

Plan, formulate, discuss and correlate a NASTRAN finite element vibrations model of the Boeing Model 360 helicopter airframe

Boeing Helicopter, together with other United States helicopter manufacturers, participated in a finite element applications program to emplace in the United States a superior capability to utilize finite element analysis models in support of helicopter airframe design. The activities relating to planning and creating a finite element vibrations model of the Boeing Model 36-0 composite airframe are summarized, along with the subsequent analytical correlation with ground shake test data.

Gabel, R.↗

NASA NTRS1991

Finite element solution of transient fluid-structure interaction problems

A finite element approach using NASTRAN is developed for solving time-dependent fluid-structure interaction problems, with emphasis on the transient scattering of acoustic waves from submerged elastic structures. Finite elements are used for modeling both structure and fluid domains to facilitate the graphical display of the wave motion through both media. For the liquid, the use of velocity potential as the fundamental unknown results in a symmetric matrix equation. The approach is illustrated for the problem of transient scattering from a submerged elastic spherical shell subjected to an incident tone burst. The use of an analogy between the equations of elasticity and the wave equation of acoustics, a necessary ingredient to the procedure, is summarized.

Everstine, Gordon C.↗

NASA NTRS1993

Finite element study of plate buckling induced by spatial temperature gradients

Finite element analyses of thermal buckling of thin metallic plates with prescribed spatial temperature distributions are described. Thermally induced compressive membrane stresses and transverse plate displacement imperfections initiate plates buckling. A finite element formulation based on von Karman plate theory is presented. The resulting nonlinear equations are solved for incremental temperature increases by Newton-Raphson iteration. The computational method is used to investigate the buckling response of rectangular plates with steady and unsteady spatially varying temperature distributions. The role of initial plate imperfections and temperature distributions on the nonlinear response of plate displacements and stresses is described. The relatively high levels of stress induced by spatial temperature gradients should be considered carefully in the postbuckling design of panels for aerospace vehicles subjected to combined mechanical and thermal loads.

Thornton, Earl A.↗

NASA NTRS1986

Anisotropic constitutive model for nickel base single crystal alloys: Development and finite element implementation

A tool for the mechanical analysis of nickel base single crystal superalloys, specifically Rene N4, used in gas turbine engine components is developed. This is achieved by a rate dependent anisotropic constitutive model implemented in a nonlinear three dimensional finite element code. The constitutive model is developed from metallurigical concepts utilizing a crystallographic approach. A non Schmid's law formulation is used to model the tension/compression asymmetry and orientation dependence in octahedral slip. Schmid's law is a good approximation to the inelastic response of the material in cube slip. The constitutive equations model the tensile behavior, creep response, and strain rate sensitivity of these alloys. Methods for deriving the material constants from standard tests are presented. The finite element implementation utilizes an initial strain method and twenty noded isoparametric solid elements. The ability to model piecewise linear load histories is included in the finite element code. The constitutive equations are accurately and economically integrated using a second order Adams-Moulton predictor-corrector method with a dynamic time incrementing procedure. Computed results from the finite element code are compared with experimental data for tensile, creep and cyclic tests at 760 deg C. The strain rate sensitivity and stress relaxation capabilities of the model are evaluated.

Dame, L. T.↗

NASA NTRS1987

Dynamic and thermal response finite element models of multi-body space structural configurations

Presented is structural dynamics modeling of two multibody space structural configurations. The first configuration is a generic space station model of a cylindrical habitation module, two solar array panels, radiator panel, and central connecting tube. The second is a 15-m hoop-column antenna. Discussed is the special joint elimination sequence used for these large finite element models, so that eigenvalues could be extracted. The generic space station model aided test configuration design and analysis/test data correlation. The model consisted of six finite element models, one of each substructure and one of all substructures as a system. Static analysis and tests at the substructure level fine-tuned the finite element models. The 15-m hoop-column antenna is a truss column and structural ring interconnected with tension stabilizing cables. To the cables, pretensioned mesh membrane elements were attached to form four parabolic shaped antennae, one per quadrant. Imposing thermal preloads in the cables and mesh elements produced pretension in the finite element model. Thermal preload variation in the 96 control cables was adjusted to maintain antenna shape within the required tolerance and to give pointing accuracy.

Edighoffer, Harold H.↗

NASA NTRS1985

A New Finite Element Approach for Prediction of Aerothermal Loads: Progress in Inviscid Flow Computations

Recent progress in the development of finite element methodology for the prediction of aerothermal loads is described. Two dimensional, inviscid computations are presented, but emphasis is placed on development of an approach extendable to three dimensional viscous flows. Research progress is described for: (1) utilization of a commerically available program to construct flow solution domains and display computational results, (2) development of an explicit Taylor-Galerkin solution algorithm, (3) closed form evaluation of finite element matrices, (4) vector computer programming strategies, and (5) validation of solutions. Two test problems of interest to NASA Langley aerothermal research are studied. Comparisons of finite element solutions for Mach 6 flow with other solution methods and experimental data validate fundamental capabilities of the approach for analyzing high speed inviscid compressible flows.

Bey, K. S.↗

NASA NTRS1985

A new finite element approach for prediction of aerothermal loads - Progress in inviscid flow computations

Recent progress in the development of finite element methodology for the prediction of aerothermal loads is described. Two dimensional, inviscid computations are presented, but emphasis is placed on development of an approach extendable to three dimensional viscous flows. Research progress is described for: (1) utilization of a commercially available program to construct flow solution domains and display computational results, (2) development of an explicit Taylor-Galerkin solution algorithm, (3) closed form evaluation of finite element matrices, (4) vector computer programming strategies, and (5) validation of solutions. Two test problems of interest to NASA Langley aerothermal research are studied. Comparisons of finite element solutions for Mach 6 flow with other solution methods and experimental data validate fundamental capabilities of the approach for analyzing high speed inviscid compressible flows.

Bey, K. S.↗

NASA NTRS1989

Mesh refinement in finite element analysis by minimization of the stiffness matrix trace

Most finite element packages provide means to generate meshes automatically. However, the user is usually confronted with the problem of not knowing whether the mesh generated is appropriate for the problem at hand. Since the accuracy of the finite element results is mesh dependent, mesh selection forms a very important step in the analysis. Indeed, in accurate analyses, meshes need to be refined or rezoned until the solution converges to a value so that the error is below a predetermined tolerance. A-posteriori methods use error indicators, developed by using the theory of interpolation and approximation theory, for mesh refinements. Some use other criterions, such as strain energy density variation and stress contours for example, to obtain near optimal meshes. Although these methods are adaptive, they are expensive. Alternatively, a priori methods, until now available, use geometrical parameters, for example, element aspect ratio. Therefore, they are not adaptive by nature. An adaptive a-priori method is developed. The criterion is that the minimization of the trace of the stiffness matrix with respect to the nodal coordinates, leads to a minimization of the potential energy, and as a consequence provide a good starting mesh. In a few examples the method is shown to provide the optimal mesh. The method is also shown to be relatively simple and amenable to development of computer algorithms. When the procedure is used in conjunction with a-posteriori methods of grid refinement, it is shown that fewer refinement iterations and fewer degrees of freedom are required for convergence as opposed to when the procedure is not used. The mesh obtained is shown to have uniform distribution of stiffness among the nodes and elements which, as a consequence, leads to uniform error distribution. Thus the mesh obtained meets the optimality criterion of uniform error distribution.

Kittur, Madan G.↗

NASA NTRS1994

Manual for automatic generation of finite element models of spiral bevel gears in mesh

The goal of this research is to develop computer programs that generate finite element models suitable for doing 3D contact analysis of faced milled spiral bevel gears in mesh. A pinion tooth and a gear tooth are created and put in mesh. There are two programs: Points.f and Pat.f to perform the analysis. Points.f is based on the equation of meshing for spiral bevel gears. It uses machine tool settings to solve for an N x M mesh of points on the four surfaces, pinion concave and convex, and gear concave and convex. Points.f creates the file POINTS.OUT, an ASCI file containing N x M points for each surface. (N is the number of node points along the length of the tooth, and M is nodes along the height.) Pat.f reads POINTS.OUT and creates the file tl.out. Tl.out is a series of PATRAN input commands. In addition to the mesh density on the tooth face, additional user specified variables are the number of finite elements through the thickness, and the number of finite elements along the tooth full fillet. A full fillet is assumed to exist for both the pinion and gear.

Bibel, G. D.↗

NASA NTRS1990

Plan, formulate, and discuss a NASTRAN finite element model of the UH-60A helicopter airframe

Under a rotorcraft structural dynamics program sponsored by the NASA Langley Research Center, Sikorsky Aircraft, together with the other major helicopter airframe manufacturers, is engaged in a study to improve the use of finite element analysis to predict the dynamic behavior of helicopter airframes. This program, which was designated DAMVIBS (Design Analysis Methods for VIBrationS), includes activities in the areas of: planning, creating, and documenting finite element models of helicopter airframes; the performance of ground vibration tests; and the correlation of test and analysis. The work performed at Sikorsky Aircraft for planning, creating, and documenting a finite element model of the UH-60A BLACK HAWK helicopter airframe is summarized. A complete description of the components of the helicopter which are to be represented in the model is presented and includes: the structural arrangement, the identification of primary and secondary structure, the components of the drive and power trains, and the attachment of large weight items to the structure. Also presented are the techniques which were used to formulate the structural finite element model for static analysis, for forming the mass and vibration models for dynamic analysis, and the procedures which were used to check out and verify the integrity of the model. Initial predictions for the vibration modes for the helicopter are included.

Dinyovszky, P.↗