Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear equations”

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 487 records · Page 27

Iterative procedures for space shuttle main engine performance models

Performance models of the Space Shuttle Main Engine (SSME) contain iterative strategies for determining approximate solutions to nonlinear equations reflecting fundamental mass, energy, and pressure balances within engine flow systems. Both univariate and multivariate Newton-Raphson algorithms are employed in the current version of the engine Test Information Program (TIP). Computational efficiency and reliability of these procedures is examined. A modified trust region form of the multivariate Newton-Raphson method is implemented and shown to be superior for off nominal engine performance predictions. A heuristic form of Broyden's Rank One method is also tested and favorable results based on this algorithm are presented.

Santi, L. Michael↗

Simulation studies using multibody dynamics code DART

DART is a multibody dynamics code developed by Photon Research Associates for the Air Force Astronautics Laboratory (AFAL). The code is intended primarily to simulate the dynamics of large space structures, particularly during the deployment phase of their missions. DART integrates nonlinear equations of motion numerically. The number of bodies in the system being simulated is arbitrary. The bodies' interconnection joints can have an arbitrary number of degrees of freedom between 0 and 6. Motions across the joints can be large. Provision for simulating on-board control systems is provided. Conservation of energy and momentum, when applicable, are used to evaluate DART's performance. After a brief description of DART, studies made to test the program prior to its delivery to AFAL are described. The first is a large angle reorientating of a flexible spacecraft consisting of a rigid central hub and four flexible booms. Reorientation was accomplished by a single-cycle sine wave shape torque input. In the second study, an appendage, mounted on a spacecraft, was slewed through a large angle. Four closed-loop control systems provided control of this appendage and of the spacecraft's attitude. The third study simulated the deployment of the rim of a bicycle wheel configuration large space structure. This system contained 18 bodies. An interesting and unexpected feature of the dynamics was a pulsing phenomena experienced by the stays whole playout was used to control the deployment. A short description of the current status of DART is given.

Keat, James E.↗

The dynamics and control of large flexible space structures - 13

The optimal control of three-dimensional large angle maneuvers and vibrations of a Shuttle-mast-reflector system is considered. The nonlinear equations of motion are formulated by using Lagrange's formula, with the mast modeled as a continuous beam subject to three-dimensional deformations. Pontryagin's Maximum Principle is applied to the slewing problem, to derive the necessary conditions for the optimal controls, which are bounded by given saturation levels. The resulting two point boundary value problem is then solved by using the quasilinearization algorithm and the method of particular solutions. The study of the large angle maneuvering of the Shuttle-beam-reflector spacecraft in the plane of a circular earth orbit is extended to consider the effects of the structural offset connection, the axial shortening, and the gravitational torque on the slewing motion. Finally the effect of additional design parameters (such as related to additional payload requirement) on the linear quadratic regulator based design of an orbiting control/structural system is examined.

Bainum, Peter M.↗

Analytic theory for the selection of a two-dimensional needle crystal at arbitrary Peclet number

An accurate analytic theory is presented for the velocity selection of a two-dimensional needle crystal for arbitrary Peclet number for small values of the surface tension parameter. The velocity selection is caused by the effect of transcendentally small terms which are determined by analytic continuation to the complex plane and analysis of nonlinear equations. The work supports the general conclusion of previous small Peclet number analytical results of other investigators, though there are some discrepancies in details. It also addresses questions raised on the validity of selection theory owing to assumptions made on shape corrections at large distances from the tip.

Tanveer, S.↗

Assessment of structural analysis technology for elastic shell collapse problems

The prediction of the ultimate load carrying capability for compressively loaded shell structures is a challenging nonlinear analysis problem. Selected areas of finite element technology research and nonlinear solution technology are assessed. Herein, a finite element analysis procedure is applied to four shell collapse problems which have been used by computational structural mechanics researchers in the past. This assessment will focus on a number of different shell element formulations and on different approaches used to account for geometric nonlinearities. The results presented confirm that these aspects of nonlinear shell analysis can have a significant effect on the predicted nonlinear structural response. All analyses were performed using the CSM Testbed software system which allowed a convenient assessment of different element formulations with a consistent approach to solving the discretized nonlinear equations.

Knight, N. F., Jr.↗

Dynamic Analysis and Control of Lightweight Manipulators with Flexible Parallel Link Mechanisms

The objective is the theoretical analysis and the experimental verification of dynamics and control of a two link flexible manipulator with a flexible parallel link mechanism. Nonlinear equations of motion of the lightweight manipulator are derived by the Lagrangian method in symbolic form to better understand the structure of the dynamic model. The resulting equation of motion have a structure which is useful to reduce the number of terms calculated, to check correctness, or to extend the model to higher order. A manipulator with a flexible parallel link mechanism is a constrained dynamic system whose equations are sensitive to numerical integration error. This constrained system is solved using singular value decomposition of the constraint Jacobian matrix. Elastic motion is expressed by the assumed mode method. Mode shape functions of each link are chosen using the load interfaced component mode synthesis. The discrepancies between the analytical model and the experiment are explained using a simplified and a detailed finite element model.

Lee, Jeh Won↗

The MHOST finite element program: 3-D inelastic analysis methods for hot section components. Volume 1: Theoretical manual

Formulations and algorithms implemented in the MHOST finite element program are discussed. The code uses a novel concept of the mixed iterative solution technique for the efficient 3-D computations of turbine engine hot section components. The general framework of variational formulation and solution algorithms are discussed which were derived from the mixed three field Hu-Washizu principle. This formulation enables the use of nodal interpolation for coordinates, displacements, strains, and stresses. Algorithmic description of the mixed iterative method includes variations for the quasi static, transient dynamic and buckling analyses. The global-local analysis procedure referred to as the subelement refinement is developed in the framework of the mixed iterative solution, of which the detail is presented. The numerically integrated isoparametric elements implemented in the framework is discussed. Methods to filter certain parts of strain and project the element discontinuous quantities to the nodes are developed for a family of linear elements. Integration algorithms are described for linear and nonlinear equations included in MHOST program.

Nakazawa, Shohei↗

Large Angle Transient Dynamics (LATDYN) demonstration problem manual

LATDYN is a computer code for modeling the Large Angle Transient Dynamics of structures. The objective in developing the code was to investigate new techniques for analyzing flexible deformation and control/structure interaction problems associated with large angular motions of spacecraft. Such motions may consist of pointing the entire spacecraft or articulation of individual components, events which occur frequently during construction, operation, and maintenance of large spacecraft. This type of analysis is beyond the routine capability of conventional analytical tools without simplifying assumptions. In some instances, the motion may be sufficiently slow and the spacecraft (or component) sufficiently rigid to simplify the analysis of dynamics and controls by making pseudo-static and/or rigid body assumptions. LATDYN introduces a new approach to the problem by combining finite element structural analysis, multi-body dynamics, and control system analysis in a single tool. It includes a new type of finite element that can deform and rotate through large angles at the same time, and which can be connected to other finite elements either rigidly or through mechanical joints. LATDYN also provides symbolic capabilities for modeling control systems which are interfaced directly with the finite element structural model. Thus, the nonlinear equations representing the structural model are integrated along with the equations representing sensors, processing, and controls as a coupled system.

Wu, Shih-Chin↗

Large Angle Transient Dynamics (LATDYN) user's manual

A computer code for modeling the large angle transient dynamics (LATDYN) of structures was developed to investigate techniques for analyzing flexible deformation and control/structure interaction problems associated with large angular motions of spacecraft. This type of analysis is beyond the routine capability of conventional analytical tools without simplifying assumptions. In some instances, the motion may be sufficiently slow and the spacecraft (or component) sufficiently rigid to simplify analyses of dynamics and controls by making pseudo-static and/or rigid body assumptions. The LATDYN introduces a new approach to the problem by combining finite element structural analysis, multi-body dynamics, and control system analysis in a single tool. It includes a type of finite element that can deform and rotate through large angles at the same time, and which can be connected to other finite elements either rigidly or through mechanical joints. The LATDYN also provides symbolic capabilities for modeling control systems which are interfaced directly with the finite element structural model. Thus, the nonlinear equations representing the structural model are integrated along with the equations representing sensors, processing, and controls as a coupled system.

Abrahamson, A. Louis↗

Spatial evolution of nonlinear acoustic mode instabilities on hypersonic boundary layers

The effects are considered of strong critical layer nonlinearity on the spatial evolution of an initially linear acoustic mode instability wave on a hypersonic flat plate boundary layer. The analysis shows that nonlinearity, which is initially confined to a thin critical layer, first becomes important when the amplitude of the pressure fluctuations become O(1/M exp 4 in M exp 2), where M is the free stream Mach number. The flow outside the critical layer is still determined by linear dynamics and therefore takes the form of a linear instability wave, but with its amplitude completely determined by the flow within the critical layer. The latter flow is determined by a coupled set of nonlinear equations, which were solved numerically.

Goldstein, M. E.↗

Sensitivity of liquid bridges subject to axial residual acceleration

It has become evident that the float zone crystal growth method and isothermal liquid bridges may be very sensitive to the residual acceleration environment of a spacelab. Using a slender-body approximation, the problem of determining the axisymmetric response of the shape of the free surface of a cylindrical liquid column bounded by two solid regions is modeled by a one-dimensional system of nonlinear equations. It is found that the sensitivity of the zone shape depends on the static Bond number, BO, aspect ratio, and viscosity, as well as the amplitude and frequency of the disturbance. The general trend is an increase in tolerable residual gravity with increasing frequency. At the eigenfrequencies of the zone, however, there are dramatic deviations from this trend. At these frequencies the tolerable residual gravity level can be two orders of magnitude lower at this frequency.

Zhang, Yiqiang↗

Nonlinear coupled rotor-fuselage helicopter vibration studies with higher harmonic control

This paper addresses the problem of vibration prediction and vibration reduction in helicopters by means of active control methodologies. The nonlinear equations of a coupled rotor/flexible-fuselage system have been derived using computer algebra, thus relegating this tedious task to the computer. In the solution procedure the trim state and vibratory response of the helicopter are obtained in a single pass by using a harmonic balance technique and simultaneously satisfying the trim and the vibratory response of the helicopter in all the rotor and fuselage degrees of freedom. Using this solution procedure, the influence of the fuselage flexibility on the vibratory response is studied. In addition, it is shown that the conventional single frequency HHC is capable of reducing either the hub loads or only the fuselage vibrations but not both simultaneously. A new scheme called MHHC, having multiple higher harmonic pitch inputs, was used to accomplish this task of simultaneously reducing both the vibratory hub loads and fuselage vibratory response. In addition, the uniqueness of this MHHC scheme is explained in detail.

Friedmann, P. P.↗

The structure of a passive scalar field with a uniform mean gradient in rapidly sheared homogeneous turbulent flow

The effect of an arbitrarily oriented mean passive scalar gradient on one-point passive scalar statistics is studied in homogeneous turbulent shear flow in the limit of rapid shearing. By neglecting the nonlinear inertial transfer to small scales an analytical solution for individual Fourier modes is obtained for the case of unity Prandtl number. This solution is used to compute the development of one-point statistics for the velocity and scalar fields in the inviscid limit. Comparisons to direct numerical simulations of the full nonlinear equations for the same flow show that, in addition to describing the early time response to the imposed shear, the linear solution gives reasonable estimates of several correlation coefficients for the developed shear flow.

Rogers, Michael M.↗

Coupled rotor-flexible fuselage vibration reduction using open loop higher harmonic control

A fundamental study of vibration prediction and vibration reduction in helicopters using active controls was performed. The nonlinear equations of motion for a coupled rotor/flexible fuselage system have been derived using computer algebra on a special purpose symbolic computer facility. The trim state and vibratory response of the helicopter are obtained in a single pass by applying the harmonic balance technique and simultaneously satisfying the trim and the vibratory response of the helicopter for all rotor and fuselage degrees of freedom. The influence of the fuselage flexibility on the vibratory response is studied. It is shown that the conventional single frequency higher harmonic control is capable of reducing either the hub loads or only the fuselage vibrations but not both simultaneously. It is demonstrated that for simultaneous reduction of hub shears and fuselae vibrations a new scheme called multiple higher harmonic control is required.

Papavassiliou, I.↗

Generalized multi-point inverse airfoil design

In a rather general sense, inverse airfoil design can be taken to mean the problem of specifying a desired set of airfoil characteristics, such as the airfoil maximum thickness ratio, pitching moment, part of the velocity distribution or boundary-layer development, etc., then from this information determine the corresponding airfoil shape. This paper presents a method which approaches the design problem from this perspective. In particular, the airfoil is divided into segments along which, together with the design conditions, either the velocity distribution or boundary-layer development may be prescribed. In addition to these local desired distributions, single parameters like the airfoil thickness can be specified. The problem of finding the airfoil shape is determined by coupling an incompressible, inviscid, inverse airfoil design method with a direct integral boundary-layer analysis method and solving the resulting nonlinear equations via a multidimensional Newton iteration technique. The approach is fast and easily allows for interactive design. It is also flexible and could be adapted to solving compressible, inverse airfoil design problems.

Selig, Michael S.↗

Computer Program For Variable-Conductance Heat Pipes

VCHPDA provides accurate mathematical models of transient as well as steady-state performance of variable-conductance heat pipes over wide range of operating conditions. Applies to heat pipes with either cold, wicked or hot, nonwicked gas reservoirs and uses ideal-gas law and "flat-front" (negligible vapor diffusion) gas theory. Calculates length of gas-blocked region and temperature of vapor in active portion of heat pipe by solving set of nonlinear equations for conservation of energy and mass. Written in FORTRAN 77.

Antoniuk, D.↗

Dynamics And Control Of Flexible Manipulator

Nonlinear equations solved with linear approximations applicable during small increments of time. Report presents theoretical study of selected aspects of dynamics and control of robotic manipulator that includes multiple flexible links. Addresses problem of control primarily from perspective of inverse-dynamics problem.

Gawronski, Wodek K.↗

Projection methods for the numerical solution of Markov chain models

Projection methods for computing stationary probability distributions for Markov chain models are presented. A general projection method is a method which seeks an approximation from a subspace of small dimension to the original problem. Thus, the original matrix problem of size N is approximated by one of dimension m, typically much smaller than N. A particularly successful class of methods based on this principle is that of Krylov subspace methods which utilize subspaces of the form span(v,av,...,A(exp m-1)v). These methods are effective in solving linear systems and eigenvalue problems (Lanczos, Arnoldi,...) as well as nonlinear equations. They can be combined with more traditional iterative methods such as successive overrelaxation, symmetric successive overrelaxation, or with incomplete factorization methods to enhance convergence.

Saad, Youcef↗