Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear operator theory”

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 181 records · Page 10

Fabrication and evaluation of polymeric early-warning fire-alarm devices

The electrical resistivities were investigated of some polymers known to be enhanced by the presence of certain gases. This was done to make a device capable of providing early warning to fire through its response with the gases produced in the early phases of combustion. Eight polymers were investigated: poly(phenyl acetylene), poly(p-aminophenyl acetylene), poly(p-nitrophenyl acetylene), poly(p-formamidophenyl acetylene), poly(ethynyl ferrocene), poly(ethynyl carborane), poly(ethynyl pyridine), and the polymer made from 1,2,3,6 tetramethyl pyridazine. A total of 40 usable thin-film sandwich devices and a total of 70 usable interdigitated-electrode lock-and-key devices were fabricated. The sandwich devices were used for measurements of contact linearity, polymer conductivity, and polymer dielectric constant. The lock-and-key devices were used to determine the response of the polymers to a spectrum of gases that included ammonia, carbon nonoxide, carbon dioxide, sulfur dioxide, ethylene, acrolein, water vapor, and normal laboratory air. Strongest responses were to water vapor, ammonia, and acrolein, and depending on the polymer, weaker responses to carbon dioxide, sulfur dioxide, and carbon monoxide were observed. A quantitative theory of device operation, capable of accounting for observed device leakage current and sensitivity, was developed. A prototype detection/alarm system was designed and built for use in demonstrating sensor performance.

Senturia, S. D.↗

Numerical analysis of right-half plane zeros for a single-link manipulator

The purpose of this research is to further develop an understanding of how nonminimum phase zero location is affected by structural link design. As the demand for light-weight robots that can operate in a large workspace increases, the structural flexibility of the links become more of an issue in controls problems. When the objective is to accurately position the tip while the robot is actuated at the base, the system is nonminimum phase. One important characteristic of nonminimum phase systems is system zeros in the right half of the Laplace plane. The ability to pick the location of these nonminimum phase zeros would give the designer a new freedom similar to pole placement. The research targets a single-link manipulator operating in the horizontal plane and modeled as a Euler-Bernoulli beam with pinned-free end conditions. Using transfer matrix theory, one can consider link designs that have variable cross-sections along the length of the beam. A FORTRAN program was developed to determine the location of poles and zeros given the system model. The program was used to confirm previous research on nonminimum phase systems, and develop a relationship for designing linearly tapered links. The method allows the designer to choose the location of the first pole and zero and then defines the appropriate taper to match the desired locations. With the pole and zero location fixes, the designer can independently change the link's moment of inertia about its axis of rotation by adjusting the height of the beam. These results can be applied to inverse dynamic algorithms currently under development at Georgia Tech.

Girvin, Douglas Lynn↗

Numerical analysis of nonminimum phase zero for nonuniform link design

As the demand for light-weight robots that can operate in a large workspace increases, the structural flexibility of the links becomes more of an issue in control. When the objective is to accurately position the tip while the robot is actuated at the base, the system is nonminimum phase. One important characteristic of nonminimum phase systems is system zeros in the right half of the Laplace plane. The ability to pick the location of these nonminimum phase zeros would give the designer a new freedom similar to pole placement. This research targets a single-link manipulator operating in the horizontal plane and modeled as a Euler-Bernoulli beam with pinned-free end conditions. Using transfer matrix theory, one can consider link designs that have variable cross-sections along the length of the beam. A FORTRAN program was developed to determine the location of poles and zeros given the system model. The program was used to confirm previous research on nonminimum phase systems, and develop a relationship for designing linearly tapered links. The method allows the designer to choose the location of the first pole and zero and then defines the appropriate taper to match the desired locations. With the pole and zero location fixed, the designer can independently change the link's moment of inertia about its axis of rotation by adjusting the height of the beam. These results can be applied to the inverse dynamic algorithms that are currently under development.

Girvin, Douglas L.↗

Anomaly Detection in Test Equipment via Sliding Mode Observers

Nonlinear observers were originally developed based on the ideas of variable structure control, and for the purpose of detecting disturbances in complex systems. In this anomaly detection application, these observers were designed for estimating the distributed state of fluid flow in a pipe described by a class of advection equations. The observer algorithm uses collected data in a piping system to estimate the distributed system state (pressure and velocity along a pipe containing liquid gas propellant flow) using only boundary measurements. These estimates are then used to further estimate and localize possible anomalies such as leaks or foreign objects, and instrumentation metering problems such as incorrect flow meter orifice plate size. The observer algorithm has the following parts: a mathematical model of the fluid flow, observer control algorithm, and an anomaly identification algorithm. The main functional operation of the algorithm is in creating the sliding mode in the observer system implemented as software. Once the sliding mode starts in the system, the equivalent value of the discontinuous function in sliding mode can be obtained by filtering out the high-frequency chattering component. In control theory, "observers" are dynamic algorithms for the online estimation of the current state of a dynamic system by measurements of an output of the system. Classical linear observers can provide optimal estimates of a system state in case of uncertainty modeled by white noise. For nonlinear cases, the theory of nonlinear observers has been developed and its success is mainly due to the sliding mode approach. Using the mathematical theory of variable structure systems with sliding modes, the observer algorithm is designed in such a way that it steers the output of the model to the output of the system obtained via a variety of sensors, in spite of possible mismatches between the assumed model and actual system. The unique properties of sliding mode control allow not only control of the model internal states to the states of the real-life system, but also identification of the disturbance or anomaly that may occur.

Solano, Wanda M.↗

A Threshold Theory for Phase-Locked Loops

A model of a phase-locked loop has been developed which is valid for all signal-to-noise ratios. The model is in the form of a nonlinear feedback system with randomly time-varying parameters. The analysis considers two operating regions. In low signal-to-noise ratio regions, the important consideration is stability. We want to study the asymptotic stability in the mean of a nonlinear system. It follows directly that a necessary condition for asymptotic stability of any nonlinear system is that a linearized model about some equilibrium point be asymptotically stable. By considering all possible equilibrium points, we can find an upper bound on the value of noise density which makes the system unstable. This upper bound represents a threshold value for system operation. In high signal-to-noise ratio regions, our results provide an exact statistical. description of system behavior. Therefore, knowledge of the spectrum of the signal and noise may be used to optimize the system configuration.

Van Trees, H. L.↗

A Koopman-based residual modeling approach for the control of a soft robot arm

Soft robots are challenging to model and control due to their poorly defined kinematics and nonlinear dynamics. Recently, Koopman operator theory has been shown capable of constructing control-oriented soft robot models from data. However, building these models requires extensive data collection and they do not necessarily generalize well outside of the training observations. This paper presents a more data-efficient and generalizable approach to soft robot modeling that first identifies a physics-based Koopman model then supplements it with a data-driven residual Koopman model. The resulting combined model is linear and thus compatible with real-time model-based control techniques such as Model Predictive Control (MPC). The efficacy of the approach is demonstrated on several simulated systems and on a real soft robot arm, where it is shown to generate models that are more accurate than purely physics-based models and require less data to construct than purely data-driven models. Using a model-based controller, the soft arm is able to successfully track end effector trajectories, perform a pick-and-place task, and write on a dry-erase board, showcasing the applicability of this framework to increase the capabilities of soft robotic systems.

Daniel Bruder↗

A flight-dynamic helicopter mathematical model with a single flap-lag-torsion main rotor

A mathematical model of a helicopter system with a single main rotor that includes rigid, hinge-restrained rotor blades with flap, lag, and torsion degrees of freedom is described. The model allows several hinge sequences and two offsets in the hinges. Quasi-steady Greenberg theory is used to calculate the blade-section aerodynamic forces, and inflow effects are accounted for by using three-state nonlinear dynamic inflow model. The motion of the rigid fuselage is defined by six degrees of freedom, and an optional rotor rpm degree of freedom is available. Empennage surfaces and the tail rotor are modeled, and the effect of main-rotor downwash on these elements is included. Model trim linearization, and time-integration operations are described and can be applied to a subset of the model in the rotating or nonrotating coordinate frame. A preliminary validation of the model is made by comparing its results with those of other analytical and experimental studies. This publication presents the results of research compiled in November 1989.

Takahashi, Marc D.↗

Recent Development in the Analytical Investigation of Shipboard Rotorcraft Engage and Disengage Operations

This paper presents an overview of recent developments in an effort to predict transient aeroelastic rotor response during shipboard engage and disengage sequences. The blade is modeled as an elastic beam undergoing in flap, lag, extension and torsion. The blade equations of motion are formulated using Hamilton's principle and they are spatially discretized using the finite element method. The discretized blade equations of motion are integrated for a specified rotor speed run-up or run-down profile. Blade element theory is used to calculate quasi-steady or unsteady aerodynamic loads in linear and nonlinear regimes. The analysis is capable of simulating both articulated, hingeless, and gimballed rotor systems. Validation of the rotor code is discussed, including correlation with droop stop impact tests and wind tunnel experiments. Predictions of safe engagement and disengagement envelopes, limited by excessive blade tip deflections or hub moments, are presented. Future directions of study are also discussed.

Smith, Edward C.↗

The role of modern control theory in the design of controls for aircraft turbine engines

The development, applications, and current research in modern control theory (MCT) are reviewed, noting the importance for fuel-efficient operation of turbines with variable inlet guide vanes, compressor stators, and exhaust nozzle area. The evolution of multivariable propulsion control design is examined, noting a basis in a matrix formulation of the differential equations defining the process, leading to state space formulations. Reports and papers which appeared from 1970-1982 which dealt with problems in MCT applications to turbine engine control design are outlined, including works on linear quadratic regulator methods, frequency domain methods, identification, estimation, and model reduction, detection, isolation, and accommodation, and state space control, adaptive control, and optimization approaches. Finally, NASA programs in frequency domain design, sensor failure detection, computer-aided control design, and plant modeling are explored

Zeller, J.↗

Tilt-rotor flutter control in cruise flight

Tilt-rotor flutter control under cruising operation is analyzed. The rotor model consists of a straight fixed wing, a pylon attached to the wingtip, and a three-blade rotor. The wing is cantilevered to the fuselage and is allowed to bend forward and upward. It also has a torsional degree of freedom about the elastic axis. Each rotor blade has two bending degrees of freedom. Feedback of wingtip velocity and acceleration to cyclic pitch is investigated for flutter control, using strip theory and linearized equations of motion. To determine the feedback gain, an eigenvalue analysis is performed. A second, independent, timewise calculation is conducted to evaluate the control law while employing more sophisticated aerodynamics. The effectiveness of flutter control by cyclic pitch change was confirmed.

Nasu, Ken-Ichi↗

Eigenfunction methods in magnetospheric radial-diffusion theory

Complete sets of orthonormal basis functions constructed according to a generalization of the quantum-mechanical WKB approximation can be used to generate a nearly-diagonal matrix representation of the radial-transport operator for ring-current ions in the presence of radial diffusion and charge exchange. The resulting eigenfunctions (constructed by weighting the basis functions in proportion to the respective components of the eigenvectors of the matrix representation) and eigenvalues provide a spatial and temporal description of the evolving phase-space density during and following a magnetospheric disturbance (e.g., a magnetic storm). A linear superposition of the basis functions can also be used to eliminate any discrepancy between the steady-state solution of the transport equation and the appropriate WKB approximation of this steady-state solution.

Schulz, Michael↗

Forced oscillations in quadratically damped systems

Bayliss (1975) has studied the question whether in the case of linear differential equations the relationship between the stability of the homogeneous equations and the existence of almost periodic solutions to the inhomogeneous equation is preserved by finite difference approximations. In the current investigation analogous properties are considered for the case in which the damping is quadratic rather than linear. The properties of the considered equation for arbitrary forcing terms are examined and the validity is proved of a theorem concerning the characteristics of the unique solution. By using the Lipschitz continuity of the mapping and the contracting mapping principle, almost periodic solutions can be found for perturbations of the considered equation. Attention is also given to the Lipschitz continuity of the solution operator and the results of numerical tests which have been conducted to test the discussed theory.

Bayliss, A.↗

Rotor Dynamic State and Parameter Identification from Hovering Transients

State and parameter identifications based on a form of the maximum likelihood method are applied to the problem of extracting linear perturbation models, including rotor dynamic inflow effects, from transient blade flapping measurements. The estimation method is first studied in computer simulations and then applied to cyclic pitch stirring transients generated with a four-bladed rotor model operating in hovering trim conditions. The analytical perturbation models extracted from the transient test results are compared with transient and frequency response tests not used in the state and parameter identification. The identified analytical perturbation model is also compared with a simple theory. The method that is applicable both to small scale and full scale dynamic rotor testing is being extended to perturbations from forward flight trim conditions.

Hohenemser, K. H.↗

One GHz leaky SAW velocity of metal layers and bilayers evaporated onto fused quartz

An acoustic microscope operating at 1 GHz was used to determine the surface acoustic wave velocities v of thin film metal layers and metal bilayers deposited onto fused quartz. V(0)'s influence was reduced by gating. This produced a calibrated accuracy of 3 percent. A program was constructed from the explicitly solved 6 x 6 theoretical determinant. Single film thicknesses were decided upon by using this theory to produce a v equal to a standard. Since the single film v's were linear with respect to their thickness, half the single film thicknesses were used for the bilayered films. The velocities for these bilayered films agreed with theory. This experimentally confirms the theoretical technique used here to examine bilayered systems, or a prototype composite interphase. No discrepancy was seen for gold films as others have reported. V(z) seemed insensitive to the formation of intermetallics or CuO. Some annealed and unannealed platinum films did not change the v from that of fused quartz. Two platinum films whose v's were in agreement with theory peeled off with the tape test.

Walikainen, Dale↗

Application of a linearized unsteady aerodynamic analysis to standard cascade configurations

A linearized potential flow analysis, which accounts for the effects of nonuniform steady flow phenomena on the linearized unsteady aerodynamic response to prescribed blade motions, has been applied to five cascade configurations. These include the first, fifth, eighth and ninth standard configurations proposed as a result of the Second International Symposium on Aeroelasticity in Turbomachines and a NASA Lewis flutter cascade. Selected results from this study, including comparisons between analytical predictions and the experimental measurements submitted for three of the foregoing configurations, are described. The correlation between theory and experiment for the first standard configuration (a compressor cascade operating at low Mach number and frequency) is quite good. Moreover, the predictions and measurements for the NASA Lewis cascade of symmetric biconvex airfoils show good qualitative agreement. However, wide discrepancies exist between the theoretical predictions and the experimental measurements for the fifth standard configuration (a subsonic transonic fan tip cascade). These can be partially attributed to conditions being imposed in the experiment which differ from those commonly used in unsteady aerodynamic analyses.

Verdon, J. M.↗

Three-Dimensional Gear Crack Propagation Studied

Gears used in current helicopters and turboprops are designed for light weight, high margins of safety, and high reliability. However, unexpected gear failures may occur even with adequate tooth design. To design an extremely safe system, the designer must ask and address the question, "What happens when a failure occurs?" With gear-tooth bending fatigue, tooth or rim fractures may occur. A crack that propagates through a rim will be catastrophic, leading to disengagement of the rotor or propeller, loss of an aircraft, and possible fatalities. This failure mode should be avoided. A crack that propagates through a tooth may or may not be catastrophic, depending on the design and operating conditions. Also, early warning of this failure mode may be possible because of advances in modern diagnostic systems. One concept proposed to address bending fatigue fracture from a safety aspect is a splittooth gear design. The prime objective of this design would be to control crack propagation in a desired direction such that at least half of the tooth would remain operational should a bending failure occur. A study at the NASA Lewis Research Center analytically validated the crack-propagation failsafe characteristics of a split-tooth gear. It used a specially developed three-dimensional crack analysis program that was based on boundary element modeling and principles of linear elastic fracture mechanics. Crack shapes as well as the crack-propagation life were predicted on the basis of the calculated stress intensity factors, mixed-mode crack-propagation trajectory theories, and fatigue crack-growth theories. The preceding figures show the effect of the location of initial cracks on crack propagation. Initial cracks in the fillet of the teeth produced stress intensity factors of greater magnitude (and thus, greater crack growth rates) than those in the root or groove areas of the teeth. Crack growth was simulated in a case study to evaluate crack-propagation paths. Tooth fracture was predicted from the crackgrowth simulation for an initial crack in the tooth fillet region. This was the desired failure mode for an ultrasafe design. Lastly, tooth loads on the uncracked mesh of the split-tooth design were up to five times greater than those on the cracked mesh if equal deflections of the cracked and uncracked teeth were considered. This effect needs to be considered in the design of a split-tooth configuration. This work was done in-house at Lewis in support of the National Rotorcraft Technology Center project, Ultra-Safe Gear Design, with the Boeing Defense and Space Group. The crack-propagation package was developed by the Cornell Fracture Group at Cornell University. The reported results, which are the initial findings of Lewis gear-crackpropagation research for the Rotorcraft Base Program, will be further investigated to develop generalized gear design guidelines.

Lewicki, David G.↗

NASTRAN cyclic symmetry capability

A development for NASTRAN which facilitates the analysis of structures made up of identical segments symmetrically arranged with respect to an axis is described. The key operation in the method is the transformation of the degrees of freedom for the structure into uncoupled symmetrical components, thereby greatly reducing the number of equations which are solved simultaneously. A further reduction occurs if each segment has a plane of reflective symmetry. The only required assumption is that the problem be linear. The capability, as developed, will be available in level 16 of NASTRAN for static stress analysis, steady state heat transfer analysis, and vibration analysis. The paper includes a discussion of the theory, a brief description of the data supplied by the user, and the results obtained for two example problems. The first problem concerns the acoustic modes of a long prismatic cavity imbedded in the propellant grain of a solid rocket motor. The second problem involves the deformations of a large space antenna. The latter example is the first application of the NASTRAN Cyclic Symmetry capability to a really large problem.

Macneal, R. H.↗

The collective excitation of g-modes in the sun

Oscillations of the solar interior (linear g-modes) may be strongly driven by the collective influence of all the modes upon the nuclear reactions in the core. This heretofore neglected effect could couple the modes, reduce the effective amplitudes near the center, and spatially concentrate most of the oscillation energy into just a portion of the radiative interior. If operating at sufficient strength, this can reverse the conventional conclusion, drawn from single mode calculations, that almost all solar g-modes are damped. Furthermore, it would put the theory into rough harmony with three otherwise troubling observations: (1) the 'low' neutrino flux measured by Davis (1978), (2) the high correspondence found by Wolff (1976) between recurrence periods in solar activity and the rotational beat periods of g-modes, and (3) the fluctuations in the sun's diameter which imply g-mode activity at high angular harmonics (Hill and Caudell 1979). A nonlinear expression is derived for the local rate of work done on an array of oscillation modes by the nuclear reactions. Three additional tests of the model are suggested.

Wolff, C. L.↗