SEARCH · Engineering Papers
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.
The evaluation of the dynamic behavior of a rotating system is possible by means of modal parameters (Eigenvalues and Eigenvectors). A mixed analytical and experimetal approach is used to identify the modal parameters of a specially designed test rig. The modal identification is done both for nonrotating as well as rotating systems. These modal parameters are used to validate a developed Finite Element Model.
Influence of nonlinear factors in fluid lubricated bearings on rotor vibrations is discussed as are computation of threshold of stability for systems with multiple degrees of freedom. The reason for this investigation was to learn more about the influence of tilting pad bearings on vertical rotors. The eigenvalues, the stability threshold, the response on the Heaviside-function and the unbalance response were calculated for a rigid rotor, a Jeffcott rotor and for the finite element model of a real rotor. The comparison of the linear and nonlinear calculations shows the limit of the linear calculations. The computation is supposed to show, that the model of the Jeffcott rotor is not sufficient in most practical cases to determine the stability threshold. This was known for some years, when computed results were compared with the behavior of real machines. That's why a program was developed, which is roughly described.
Efficient solution of the equations from spectral discretizations is essential if the high-order accuracy of these methods is to be realized. Direct solution of these equations is rarely feasible, thus iterative techniques are required. A preconditioning scheme for first-order Chebyshev collocation operators is proposed herein, in which the central finite difference mesh is finer than the collocation mesh. Details of the proper techniques for transferring information between the meshes are given here, and the scheme is analyzed by examination of the eigenvalue spectra of the preconditioned operators. The effect of artificial viscosity required in the inversion of the finite difference operator is examined. A second preconditioning scheme, involving a high-order upwind finite difference operator of the van Leer type is also analyzed to provide a comparison with the present scheme. Finally, the performance of the present scheme is verified by application to several test problems.
This paper describes a method for the identification of the parameters entering into the equations of motion of distributed systems. Because the motion of distributed systems is described in terms of partial differential equations, these parameters are in general continuous functions of the spatial variables. For vibrating systems, these parameters ordinarily represent the mass, stiffness and damping distributions. In this paper, these distributions are expanded in terms of finite series of known functions of the spatial variables multiplied by undetermined coefficients. It is assumed that the nature of the equations of motion is known and that a limited number of eigenvalues and eigenfunctions is identified in advance. Use is then made of the least squares method, in conjunction with the eigenfunctions' orthogonality, to compute the undetermined coefficients, thus identifying the system distributed parameters. A method for the identification of the eigensolution is also presented. The procedure for the identification of the eigensolution and of the system parameters is demonstrated via a numerical example.
To permit investigation of the relative effects of higher-order terms in the potential function and allow computation of the gravity gradient torques consistent with the gravitational acceleration for a given potential function, a simple formula is derived for computing the gravity gradient torque on a vehicle in terms of the eigenvalues and eigenvectors of the gravity gradient matrix and the inertia properties of the vehicle. The formulation is applicable to any potential function for which the gravity matrix can be computed for a given point in the potential field.
The present paper is concerned with an analysis of the bonding in Fe2(CO)9, taking into account the results of self-consistent-field (SCF) calculations. The basis sets used are the same as those employed in the study of Fe(CO)5 conducted by Bauschlicher and Bagus (1985). The Fe basis set starts with the 14s9p5d primitive set optimized by Wachters (1970), which is contracted to 8s4p5d. An SCF wave function, consisting of only the core electrons, was determined, and then electrons were slowly added to the virtual orbitals with the lowest eigenvalues. The conducted calculations show no direct Fe-Fe bonding, but rather bonding resulting from the bridging CO's. The bonding is found to involve both the 5 sigma and 2 pi(asterisk) orbitals of the bridging CO's.
A fully coupled, shock capturing technique is presented for chemically reacting flows at high Mach numbers. The technique makes use of a total variation diminishing (TVD) dissipation operator which results in sharp, crisp shocks. The eigenvalues and eigenvectors of the fully coupled system, which includes species conversion equations in addition to the gas dynamics equations, are analytically derived for a general reacting gas. Species production terms for a model dissociating gas are introduced and are included in the algorithm. The convective terms are solved using a first-order TVD scheme while the source terms are solved using a fourth-order Runge-Kutta scheme to enhance stability. Results from one-dimensional numerical experiments are shown for a two species and a three species gas.
This investigation represents an extension of a study of Roussos (1985) who considered the noise transmission loss of a rectangular plate in an infinite baffle. Roussos, who employed an analytical formulation, considered an unstiffened plate. While it is difficult to consider stiffeners by means of analytical methods, the difficulties can be avoided by employing a finite element procedure. For this reason, the present study is concerned with the implementation of a finite element method. The representation of the panel transmission loss is discussed, and the determination of the panel motion by means of the finite element technique is described, taking into account an isotropic flat panel, the exciting force, an eigenvalue problem, the radiation pressure, a plate element, and a cylindrical shell element. Numerical results are considered for a flat panel, a curved panel, and a stiffened flat panel.
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.
The flow of a fluid between concentric rotating cylinders (the Taylor problem) is studied by exploiting the symmetries of the system. The Navier-Stokes equations, linearized about Couette flow, possess two zero and four purely imaginary eigenvalues at a suitable value of the speed of rotation of the outer cylinder. There is thus a reduced bifurcation equation on a six-dimensonal space which can be shown to commute with an action of the symmetry group 0(2) x S0(2). The group structure is used to analyze this bifurcation equation in the simplest (nondegenerate) case, and to compute the stabilities of solutions. In particular, when the outer cylinder is counterrotated, transitions which seem to agree with recent experiments of Andereck, Liu, and Swinney (1984) are obtained. It is also possible to obtain the 'main sequence' in this model. This sequence is normally observed in experiments when the outer cylinder is held fixed.
A new interactive computer software package for eigenvalue/eigenvector assignment using constant state feedback is described. The package consists of ten subprograms, each associated with a specific design objective, accessible from a main control program. Using this package, primary design objectives of assigning eigenvalues and approximating eigenvectors are first achieved. Then secondary design objectives, including modification of specific eigenvector components, reduction in specified elements of the feedback gain matrix, and reduction in eigensystem sensitivity to changes in plant parameters, are addressed. These secondary objectives are achieved by a systematic modification of the assigned eigenvectors in a small region about the initial assignment. In addition, the program implements the use of spectral assignment procedures with reduced-order system models. Program modes are described and illustrated by numerical examples.
Eigenspace techniques allow the control system designer to use feedback to place eigenvalues and shape eigenvectors so as to modify closed loop dynamic response characteristics to achieve performance objectives. In this paper the theory of eigenspace design is reviewed, extended and applied to several aircraft control problems.
The design of a maneuver autopilot for flight test trajectory control using constrained eigenvalue/eigenvector assignment is examined. The aircraft considered was a high-performance fighter with a command augmentation system engaged in all three axes. Attention is given to difficulties encountered in the generation of the desired eigenvalues and eigenvectors. It is found that this approach demands several iterations to converge to a satisfactory result, and does not appear to easily yield suitable insight for the output feedback design of high-order multivariable systems which will be used at other operating points. It is concluded that this technique could be made more attractive by generating gradients of the eigensystem between flight conditions, and including this information in the single-point design technique.
The requirements of complex aerospace vehicles combined with the age of structural analysis systems enhance the need to advance technology toward a new generation of structural analysis capability. Recent and impeding advances in parallel and supercomputers provide the opportunity to significantly improve these structural analysis capabilities for large order finite element problems. Long-term research in parallel computing, associated with the NASA Finite Element Machine project, is discussed. The results show the potential of parallel computers to provide substantial increases in computation speed over sequential computers. Results are given for sample problems in the areas of eigenvalue analysis and transient response.
The growth of the momentum thickness and the modal disturbance energies are examined to study the nature and onset of nonlinearity in a temporally growing free shear layer. A shooting technique is used to find solutions to the linearized eigenvalue problem, and pseudospectral weakly nonlinear simulations of this flow are obtained for comparison. The roll-up of a fundamental disturbance follows linear theory predictions even with a 20 percent disturbance amplitude. A weak nonlinear interaction of the disturbance creates a finite-amplitude mean shear stress which dominates the growth of the layer momentum thickness, and the disturbance growth rate changes until the fundamental disturbance dominates. The fundamental then becomes an energy source for the harmonic, resulting in an increase in the growth rate of the subharmonic over the linear prediction even when the fundamental has no energy to give. Also considered are phase relations and the wall influence.
The exact FORTRAN computer program BUNVIS-RG which uses an exact stiffness matrix method to find eigenvalues and modes of three-dimensional planes, is presented in detail. It is shown that BUNVIS-RG solution times can be predicted fairly accurately from two computer dependent constants, the number of nodes and members in the frame, and the maximum node number difference of any pair of connected nodes. The program's use of repetitive geometry, stayed column substructures and nonuniform member options was found to enhance efficiency and range of application.
Methods for calculating sensitivity derivatives for discrete structural systems are surveyed, primarily covering literature published during the past two decades. Methods are described for calculating derivatives of static displacements and stresses, eigenvalues and eigenvectors, transient structural response, and derivatives of optimum structural designs with respect to problem parameters. The survey is focused on publications developed in nonstructural fields such as electronics, controls, and physical chemistry which are directly applicable to structural problems. Most notable among the nonstructural-based methods are the adjoint variable technique from control theory, and the Green's function and FAST methods from physical chemistry.
A general relationship is established between wave-packet evolution and quantization in classically integrable systems. Because of wave-packet spreading, one cannot simply take the Fourier transform of the time evolution of a wave-packet. Instead, one must propagate the wave-packet using the actions as Hamiltonians. The energy eigenvalues which result are the Einstein-Brillouin-Keller values, and new forms for the eigenfunctions appear. These are free of caustic singularities, and represent averages of wave packets over the invariant torus.