Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite-element”

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 433 records · Page 24

Integration of dynamic, aerodynamic and structural optimization of helicopter rotor blades

The purpose of the research is to study the integration of structural, dynamic, and aerodynamic considerations in the design-optimization process for helicopter rotorblades. This is to be done in three phases. Task 1 is to bring on-line computer codes that could perform the finite-element frequency analyses of rotor blades. The major features of this program are summarized. The second task was to bring on-line an optimization code for the work. Several were tried and it was decided to use CONMIN. Explicit volume constraints on the thicknesses and lumped masses used in the optimization were added. The specific aeroelastic constraint that the center of mass must be forward of the quarter chord in order to prevent flutter was applied. The bending-torsion coupling due to cg-ea offset within the blade cross section was included. Also included were some very simple stress constraints. The first three constraints are completed, and the fourth constraint is being completed.

Peters, David A.↗

Analytical modeling of helicopter static and dynamic induced velocity in GRASP

The methodology used by the General Rotorcraft Aeromechanical Stability Program (GRASP) to model the characteristics of the flow through a helicopter rotor in hovering or axial flight is described. Since the induced flow plays a significant role in determining the aeroelastic properties of rotorcraft, the computation of the induced flow is an important aspect of the program. Because of the combined finite-element/multibody methodology used as the basis for GRASP, the implementation of induced velocity calculations presented an unusual challenge to the developers. To preserve the modelling flexibility and generality of the code, it was necessary to depart from the traditional methods of computing the induced velocity. This is accomplished by calculating the actuator disc contributions to the rotor loads in a separate element called the air mass element, and then performing the calculations of the aerodynamic forces on individual blade elements within the aeroelastic beam element.

Kunz, Donald L.↗

Constitutive modeling for isotropic materials

The objective is to develop a unified constitutive model for finite element structural analysis of turbine engine hot-section components. This effort constitutes a different approach for non-linear finite-element computer codes which have heretofore been based on classical inelastic methods. The unified constitutive theory to be developed will avoid the simplifying assumptions of classical theory and should more accurately represent the behavior of superalloy materials under cyclic loading conditions and high temperature environments. During the first two years of the program, extensive experimental correlations were made with two representative unified models. The experiments were both uniaxial and biaxial at temperatures up to 1093 C (2000 F). In addition, the unified models were adopted to the MARC finite element code and used for stress analysis of notched bar and turbine blade geometries.

Lindholm, Ulric S.↗

Comparison of measured temperatures, thermal stresses and creep residues with predictions on a built-up titanium structure

Temperature, thermal stresses, and residual creep stresses were studied by comparing laboratory values measured on a built-up titanium structure with values calculated from finite-element models. Several such models were used to examine the relationship between computational thermal stresses and thermal stresses measured on a built-up structure. Element suitability, element density, and computational temperature discrepancies were studied to determine their impact on measured and calculated thermal stress. The optimum number of elements is established from a balance between element density and suitable safety margins, such that the answer is acceptably safe yet is economical from a computational viewpoint. It is noted that situations exist where relatively small excursions of calculated temperatures from measured values result in far more than proportional increases in thermal stress values. Measured residual stresses due to creep significantly exceeded the values computed by the piecewise linear elastic strain analogy approach. The most important element in the computation is the correct definition of the creep law. Computational methodology advances in predicting residual stresses due to creep require significantly more viscoelastic material characterization.

Jenkins, Jerald M.↗

Composite mechanics for engine structures

Recent research activities and accomplishments at Lewis Research Center on composite mechanics for engine structures are summarized. The activities focused mainly on developing procedures for the computational simulation of composite intrinsic and structural behavior. The computational simulation encompasses all aspects of composite mechanics, advanced three-dimensional finite-element methods, damage tolerance, composite structural and dynamic response, and structural tailoring and optimization.

Chamis, Christos C.↗

Steady and transient least square solvers for thermal problems

This paper develops a hierarchical least square solution algorithm for highly nonlinear heat transfer problems. The methodology's capability is such that both steady and transient implicit formulations can be handled. This includes problems arising from highly nonlinear heat transfer systems modeled by either finite-element or finite-difference schemes. The overall procedure developed enables localized updating, iteration, and convergence checking as well as constraint application. The localized updating can be performed at a variety of hierarchical levels, i.e., degree of freedom, substructural, material-nonlinear groups, and/or boundary groups. The choice of such partitions can be made via energy partitioning or nonlinearity levels as well as by user selection. Overall, this leads to extremely robust computational characteristics. To demonstrate the methodology, problems are drawn from nonlinear heat conduction. These are used to quantify the robust capabilities of the hierarchical least square scheme.

Padovan, Joe↗

Numerical methods for one-dimensional reaction-diffusion equations arising in combustion theory

A review of numerical methods for one-dimensional reaction-diffusion equations arising in combustion theory is presented. The methods reviewed include explicit, implicit, quasi-linearization, time linearization, operator-splitting, random walk and finite-element techniques and methods of lines. Adaptive and nonadaptive procedures are also reviewed. These techniques are applied first to solve two model problems which have exact traveling wave solutions with which the numerical results can be compared. This comparison is performed in terms of both the wave profile and computed wave speed. It is shown that the computed wave speed is not a good indicator of the accuracy of a particular method. A fourth-order time-linearized, Hermitian compact operator technique is found to be the most accurate method for a variety of time and space sizes.

Ramos, J. I.↗

All-theoretical prediction of cabin noise due to impingement of propeller vortices on a wing structure

Reported calculations of structure-borne cabin noise for a small twin engine aircraft powered by tractor propellers rely on the following three-stage methodological breakup of the problem: (1) the unsteady-aerodynamic prediction of wing lift harmonics caused by the whipping action of the vortex system trailed from each propeller; (2) the associated wing/fuselage structural response; (3) the cabin noise field for the computed wall vibration. The first part--the estimate of airloads--skirts a full-fledged aeroelastic situation by assuming the wing to be fixed in space while cancelling the downwash field of the cutting vortices. The model is based on an approximate high-frequency lifting-surface theory justified by the blade rate and flight Mach number of application. Its results drive a finite-element representation of the wing accounting for upper and lower skin surfaces, spars, ribs, and the presence of fuel. The fuselage, modeled as a frame-stiffened cylindrical shell, is bolted to the wing.

Martinez, R.↗

Analysis of eccentric annular incompressible seals. II - Effects of eccentricity on rotordynamic coefficients

A new analysis procedure has been presented which solves for the flow variables of an annular pressure seal in which the rotor has a large static displacement (eccentricity) from the centered position. The present paper incorporates the solutions to investigate the effect of eccentricity on the rotordynamic coefficients. The analysis begins with a set of governing equations based on a turbulent bulk-flow model and Moody's friction factor equation. Perturbations of the flow variables yields a set of zeroth- and first-order equations. After integration of the zeroth-order equations, the resulting zeroth-order flow variables are used as input in the solution of the first-order equations. Further integration of the first order pressures yields the eccentric rotordynamic coefficients. The results from this procedure compare well with available experimental and theoretical data, with accuracy just as good or slightly better than the predictions based on a finite-element model.

Nelson, C. C.↗

Convergence properties of modal costs for certain distributed parameter systems

A complete modal cost analysis is presented for the vibration of various kinds of simple continua with different boundary conditions. Explicit formulas for the norm of the response (called the 'cost') are derived for these distributed parameter systems. The convergence theorems developed are useful in the model reduction of the equivalent continuum models of large space structures as well as in the selection of finite-element code for control design.

Hu, A.↗

Isoparametric finite element methods to estimate discontinuous coefficients in two-dimensional elliptic equations

A summary of findings is presented on the problem of estimating discontinuous coefficients, including the locations of interfaces at which discontinuities occur, for elliptic equations on a two-dimensional domain. Finite-element approximations are constructed that are based on an underlying grid that adaptively moves to follow locations of interfaces as they are iteratively updated in the numerical estimation process.

Lamm, Patricia K.↗

Frequency optimization of repetitive lattice beam-like structures using a continuum model

A new method for obtaining the maximum frequency design of a beam-like repetitive lattice structure is presented. Using existing techniques, the lattice is first modeled as an equivalent anisotropic Timoshenko beam. The computation of the stiffness and inertial properties of the beam, determined by matching the strain and kinetic energies of the beam with those of the lattice, is facilitated by the repetitive nature of the lattice. The optimum design is obtained by maximizing Rayleigh's quotient using methods of variational calculus. For the problem selected, results show excellent agreement with those obtained by traditional finite-element methods. Moreover, unlike FE methods, cpu time is relatively unaffected by the size of the truss.

Reiss, Robert↗

Analytical approach to peel stresses in bonded composite stiffened panels

A closed-form solution was obtained for the stresses and displacements of two bonded beams. A system of two fourth-order and two second-order differential equations with the associated boundary equations was determined using a variational work approach. A FORTRAN computer program was devised to solve for the eigenvalues and eigenvectors of this system and to calculate the coefficients from the boundary conditions. The results were then compared with NASTRAN finite-element solutions and shown to agree closely.

Barkey, Derek A.↗

Constitutive modeling for isotropic materials

The first year of progress on a NASA-Lewis contract with the General Electric Co is documented. The purpose of this contract (NAS3-23927) is to develop and evaluate unified constitutive equations for applications to hot-path components of aircraft gas turbine engines such as high pressure turbine blades and vanes. To accomplish this goal, uniaxial, notched, and multiaxial specimens made of conventionally cast Rene 80 are being tested under conditions that simulate engine operating conditions. To reduce the raw data, automated data reduction techniques are being developed that produce computer files containing the information needed to analyze proposed constitutive theories. Described are the analytical methods being developed to determine the parameters for these nonlinear unified theories by using the reduced data files. In another activity, a dedicated finite-element computer code is being developed to use unified theories in the structural analysis of hot-section components. This code was extensively verified for one such theory by successfully predicting the strain histories measured experimentally at the notch root of complex specimens taken from complex laboratory specimens.

Ramaswamy, V. G.↗

Recent developments in the experimental identification of the dynamics of a highly flexible grid

Control effectiveness tests of reaction wheel actuators attached to a highly flexible grid are reported. Analytic determination of actuator control effectiveness is accomplished with finite-element modeling. Experimental determination is done with a least-square parameter identification algorithm that identifies the control coefficients of the second-order difference equation model of each vibration mode. The algorithm assumes a model with frequency and damping predetermined from free-decay tests for each mode. Accounting for the difference in forced and resonant frequency was necessary to produce control effectiveness estimates that are in reasonable agreement with the analytic predictions. The average error for control effectiveness coefficients greater than 5/sq sec was 6.384 percent.

Montgomery, Raymond C.↗

Constitutive modeling for isotropic materials (HOST)

The results of the third year of work on a program which is part of the NASA Hot Section Technology program (HOST) are presented. The goals of this program are: (1) the development of unified constitutive models for rate dependent isotropic materials; and (2) the demonstration of the use of unified models in structural analyses of hot section components of gas turbine engines. The unified models selected for development and evaluation are those of Bodner-Partom and of Walker. A test procedure was developed for assisting the generation of a data base for the Bodner-Partom model using a relatively small number of specimens. This test procedure involved performing a tensile test at a temperature of interest that involves a succession of strain-rate changes. The results for B1900+Hf indicate that material constants related to hardening and thermal recovery can be obtained on the basis of such a procedure. Strain aging, thermal recovery, and unexpected material variations, however, preluded an accurate determination of the strain-rate sensitivity parameter is this exercise. The effects of casting grain size on the constitutive behavior of B1900+Hf were studied and no particular grain size effect was observed. A systematic procedure was also developed for determining the material constants in the Bodner-Partom model. Both the new test procedure and the method for determining material constants were applied to the alternate material, Mar-M247 . Test data including tensile, creep, cyclic and nonproportional biaxial (tension/torsion) loading were collected. Good correlations were obtained between the Bodner-Partom model and experiments. A literature survey was conducted to assess the effects of thermal history on the constitutive behavior of metals. Thermal history effects are expected to be present at temperature regimes where strain aging and change of microstructure are important. Possible modifications to the Bodner-Partom model to account for these effects are outlined. The use of a unified constitutive model for hot section component analyses was demonstrated by applying the Walker model and the MARC finite-element code to a B1900+Hf airfoil problem.

Chan, Kwai S.↗

The Effects of Variable Mass and Geometry, Pretwist, Shear Deformation and Rotatory Inertia on the Resonant Frequencies of Intact Long Bones: A Finite Element Model Analysis

The influence of pretwist, nonuniformities in mass and flexural stiffness, rotatory inertia and shear deformation on the natural frequencies of intact bones is evaluated by means of a linear elastic, finite-element model which has been programmed for solution on the digital computer. Theoretical results are compared to the results on the forced vibration of intact canine radii obtained experimentally by Thompson. Surprisingly, inclusion of fairly large pretwist angles (from -14 to 12 deg for one specimen) had little affect on the first three frequencies of transverse vibration in either the cranial or lateral directions. Inclusion of shear deformation reduced the third-mode frequency in the stiffest (lateral) direction by about six percent, otherwise shear deformation played a minor role in determining natural frequencies. Similarly. rotatory inertia had negligible influence up to the third natural frequency. The predominant influence on the first three natural frequencies of transverse vibration could be attributed to the variations in mass and flexural stiffness along the length of the test specimens. Different effective moduli of elasticity are required to yield correct absolute values for the frequencies which correspond to experimental findings. thus implying the presence of some inhomogeneities in material properties around the bone cross-section and/or along its length.

Young, Donald R.↗

Subtle alternating electrocardiographic morphology as an indicator of decreased cardiac electrical stability

Observations from finite-element computer models, together with analytic developments based on percolation theory have suggested that subtle fluctuations of ECG morphology might serve as an indicator diminished cardiac electrical stability. With fixed-rate atrial pacing in canines, we have previously observed a pattern of alternation in T wave energy which correlated with cardiac electrical stability. We report here on a series of 20 canine experiments in which cardiac electrical stability (measured via Ventricular Fibrillation Threshold determination) was compared to a non-degenerate, multidimensional measurement of the degree of alternating activity present in the ECG complex morphology. The decrease in cardiac electrical stability brought on by both coronary artery occlusion and systemic hypothermia was consistently accompanied by subtle alternation in ECG morphology, with the absolute degree of alternating activity being significantly (negatively) correlated with cardiac electrical stability.

Ventricular Fibrillation↗