Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bilinear control”

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 37 records · Page 2

Electromagnetic Monitoring and Control of a Plurality of Nanosatellites

A method for monitoring position of and controlling a second nanosatellite (NS) relative to a position of a first NS. Each of the first and second NSs has a rectangular or cubical configuration of independently activatable, current-carrying solenoids, each solenoid having an independent magnetic dipole moment vector, .mu.1 and .mu.2. A vector force F and a vector torque are expressed as linear or bilinear combinations of the first set and second set of magnetic moments, and a distance vector extending between the first and second NSs is estimated. Control equations are applied to estimate vectors, .mu.1 and .mu.2, required to move the NSs toward a desired NS configuration. This extends to control of N nanosatellites.

Soloway, Donald I.↗

The algorithms for rational spline interpolation of surfaces

Two algorithms for interpolating surfaces with spline functions containing tension parameters are discussed. Both algorithms are based on the tensor products of univariate rational spline functions. The simpler algorithm uses a single tension parameter for the entire surface. This algorithm is generalized to use separate tension parameters for each rectangular subregion. The new algorithm allows for local control of tension on the interpolating surface. Both algorithms are illustrated and the results are compared with the results of bicubic spline and bilinear interpolation of terrain elevation data.

Schiess, J. R.↗

Numerical simulation of fluid flow around a scramaccelerator projectile

Numerical simulations of the fluid motion and temperature distribution around a 'scramaccelerator' projectile are obtained for Mach numbers in the 5-10 range. A finite element method is used to solve the equations of motion for inviscid and viscous two-dimensional or axisymmetric compressible flow. The time-dependent equations are solved explicitly, using bilinear isoparametric quadrilateral elements, mass lumping, and a shock-capturing Petrov-Galerkin formulation. Computed results indicate that maintaining on-design performance for controlling and stabilizing oblique detonation waves is critically dependent on projectile shape and Mach number.

Pepper, Darrell W.↗

FPGA Implementation of Stereo Disparity with High Throughput for Mobility Applications

High speed stereo vision can allow unmanned robotic systems to navigate safely in unstructured terrain, but the computational cost can exceed the capacity of typical embedded CPUs. In this paper, we describe an end-to-end stereo computation co-processing system optimized for fast throughput that has been implemented on a single Virtex 4 LX160 FPGA. This system is capable of operating on images from a 1024 x 768 3CCD (true RGB) camera pair at 15 Hz. Data enters the FPGA directly from the cameras via Camera Link and is rectified, pre-filtered and converted into a disparity image all within the FPGA, incurring no CPU load. Once complete, a rectified image and the final disparity image are read out over the PCI bus, for a bandwidth cost of 68 MB/sec. Within the FPGA there are 4 distinct algorithms: Camera Link capture, Bilinear rectification, Bilateral subtraction pre-filtering and the Sum of Absolute Difference (SAD) disparity. Each module will be described in brief along with the data flow and control logic for the system. The system has been successfully fielded upon the Carnegie Mellon University's National Robotics Engineering Center (NREC) Crusher system during extensive field trials in 2007 and 2008 and is being implemented for other surface mobility systems at JPL.

Random access memory↗

A method for the analysis of nonlinearities in aircraft dynamic response to atmospheric turbulence

An analytical method is developed which combines the equivalent linearization technique for the analysis of the response of nonlinear dynamic systems with the amplitude modulated random process (Press model) for atmospheric turbulence. The method is initially applied to a bilinear spring system. The analysis of the response shows good agreement with exact results obtained by the Fokker-Planck equation. The method is then applied to an example of control-surface displacement limiting in an aircraft with a pitch-hold autopilot.

Sidwell, K.↗

GRID3D-v2: An updated version of the GRID2D/3D computer program for generating grid systems in complex-shaped three-dimensional spatial domains

In order to generate good quality systems for complicated three-dimensional spatial domains, the grid-generation method used must be able to exert rather precise controls over grid-point distributions. Several techniques are presented that enhance control of grid-point distribution for a class of algebraic grid-generation methods known as the two-, four-, and six-boundary methods. These techniques include variable stretching functions from bilinear interpolation, interpolating functions based on tension splines, and normalized K-factors. The techniques developed in this study were incorporated into a new version of GRID3D called GRID3D-v2. The usefulness of GRID3D-v2 was demonstrated by using it to generate a three-dimensional grid system in the coolent passage of a radial turbine blade with serpentine channels and pin fins.

Steinthorsson, E.↗

Optimal guidance law development for an advanced launch system

The objective of this research effort was to develop a real-time guidance approach for launch vehicles ascent to orbit injection. Various analytical approaches combined with a variety of model order and model complexity reduction have been investigated. Singular perturbation methods were first attempted and found to be unsatisfactory. The second approach based on regular perturbation analysis was subsequently investigated. It also fails because the aerodynamic effects (ignored in the zero order solution) are too large to be treated as perturbations. Therefore, the study demonstrates that perturbation methods alone (both regular and singular perturbations) are inadequate for use in developing a guidance algorithm for the atmospheric flight phase of a launch vehicle. During a second phase of the research effort, a hybrid analytic/numerical approach was developed and evaluated. The approach combines the numerical methods of collocation and the analytical method of regular perturbations. The concept of choosing intelligent interpolating functions is also introduced. Regular perturbation analysis allows the use of a crude representation for the collocation solution, and intelligent interpolating functions further reduce the number of elements without sacrificing the approximation accuracy. As a result, the combined method forms a powerful tool for solving real-time optimal control problems. Details of the approach are illustrated in a fourth order nonlinear example. The hybrid approach is then applied to the launch vehicle problem. The collocation solution is derived from a bilinear tangent steering law, and results in a guidance solution for the entire flight regime that includes both atmospheric and exoatmospheric flight phases.

Calise, Anthony J.↗

On the algebraic structure of bilinear systems.

It is shown that a particular bilinear model is both quite general and easy to work with. A basic structure theory is developed with the aid of previous results. Some preliminary ideas are discussed together with the system interconnection, the canonical form, questions of controllability, aspects of observability, and equivalent realizations. It is pointed out that in actually determining equivalent realizations for systems and in the classification of systems, the results available in the study of Lie algebras are of fundamental importance.

Brockett, R. W.↗

Sliding control of pointing and tracking with operator spline estimation

It is shown how a variable structure control technique could be implemented to achieve precise pointing and good tracking of a deformable structure subject to fast slewing maneuvers. The correction torque that has to be applied to the structure is based on estimates of upper bounds on the model errors. For a rapid rotation of the deformable structure, the elastic response can be modeled by oscillators driven by angular acceleration, and where stiffness and damping coefficients are also angular velocity and acceleration dependent. By transforming this slew-driven elastic dynamics into bilinear form (be regarding the vector made up of the angular velocity, squared angular velocity and angular acceleration components, which appear in the coefficients as the input to the deformation dynamics), an operator spline can be constructed, that gives a low order estimate of the induced disturbance. Moreover, a worst case error bound between the estimated deformation and the unknown exact deformation is also generated, which can be used where required in the sliding control correction.

Dwyer, Thomas A. W., III↗

Improved Fabrication of Ceramic Matrix Composite/Foam Core Integrated Structures

The use of hybridized carbon/silicon carbide (C/SiC) fabric to reinforce ceramic matrix composite face sheets and the integration of such face sheets with a foam core creates a sandwich structure capable of withstanding high-heatflux environments (150 W/cm2) in which the core provides a temperature drop of 1,000 C between the surface and the back face without cracking or delamination of the structure. The composite face sheet exhibits a bilinear response, which results from the SiC matrix not being cracked on fabrication. In addition, the structure exhibits damage tolerance under impact with projectiles, showing no penetration to the back face sheet. These attributes make the composite ideal for leading edge structures and control surfaces in aerospace vehicles, as well as for acreage thermal protection systems and in high-temperature, lightweight stiffened structures. By tailoring the coefficient of thermal expansion (CTE) of a carbon fiber containing ceramic matrix composite (CMC) face sheet to match that of a ceramic foam core, the face sheet and the core can be integrally fabricated without any delamination. Carbon and SiC are woven together in the reinforcing fabric. Integral densification of the CMC and the foam core is accomplished with chemical vapor deposition, eliminating the need for bond-line adhesive. This means there is no need to separately fabricate the core and the face sheet, or to bond the two elements together, risking edge delamination during use. Fibers of two or more types are woven together on a loom. The carbon and ceramic fibers are pulled into the same pick location during the weaving process. Tow spacing may be varied to accommodate the increased volume of the combined fiber tows while maintaining a target fiber volume fraction in the composite. Foam pore size, strut thickness, and ratio of face sheet to core thickness can be used to tailor thermal and mechanical properties. The anticipated CTE for the hybridized composite is managed by the choice of constituents, varying fiber tow sizes and constituent part ratios. This structural concept provides high strength and stiffness at low density 1.06 g/cm3 in panels tested. Varieties of face sheet constructions are possible, including variations in fiber type and weave geometry. The integrated structures possible with this composite could eliminate the need for non-load-bearing thermal protection systems on top of a structural component. The back sheet can readily be integrated to substructures through the incorporation of ribs. This would eliminate weight and cost for aerospace missions.

Hurwitz, Frances I.↗

Inert Welding/Brazing Gas Filters and Dryers

The use of hybridized carbon/silicon carbide (C/SiC) fabric to reinforce ceramic matrix composite face sheets and the integration of such face sheets with a foam core creates a sandwich structure capable of withstanding high-heat-flux environments (150 W/sq cm) in which the core provides a temperature drop of 1,000 C between the surface and the back face without cracking or delamination of the structure. The composite face sheet exhibits a bilinear response, which results from the SiC matrix not being cracked on fabrication. In addition, the structure exhibits damage tolerance under impact with projectiles, showing no penetration to the back face sheet. These attributes make the composite ideal for leading-edge structures and control surfaces in aerospace vehicles, as well as for acreage thermal protection systems and in high-temperature, lightweight stiffened structures. By tailoring the coefficient of thermal expansion (CTE) of a carbon fiber containing ceramic matrix composite (CMC) face sheet to match that of a ceramic foam core, the face sheet and the core can be integrally fabricated without any delamination. Carbon and SiC are woven together in the reinforcing fabric. Integral densification of the CMC and the foam core is accomplished with chemical vapor deposition, eliminating the need for bond-line adhesive. This means there is no need to separately fabricate the core and the face sheet, or to bond the two elements together, risking edge delamination during use. Fibers of two or more types are woven together on a loom. The carbon and ceramic fibers are pulled into the same "pick" location during the weaving process. Tow spacing may be varied to accommodate the increased volume of the combined fiber tows while maintaining a target fiber volume fraction in the composite. Foam pore size, strut thickness, and ratio of face sheet to core thickness can be used to tailor thermal and mechanical properties. The anticipated CTE for the hybridized composite is managed by the choice of constituents, varying fiber tow sizes and constituent part ratios. This structural concept provides high strength and stiffness at low density 1.06 g/cu cm in panels tested. Varieties of face sheet constructions are possible, including variations in fiber type and weave geometry. The integrated structures possible with this composite could eliminate the need for non-load-bearing thermal protection systems on top of a structural component. The back sheet can readily be integrated to substructures through the incorporation of ribs. This would eliminate weight and cost for aerospace missions.

Goudy, Jerry↗

A methodology for using nonlinear aerodynamics in aeroservoelastic analysis and design

A methodology is presented for using the Volterra-Wiener theory of nonlinear systems in aeroservoelastic (ASE) analyses and design. The theory is applied to the development of nonlinear aerodynamic response models that can be defined in state-space form and are, therefore, appropriate for use in modern control theory. The theory relies on the identification of nonlinear kernels that can be used to predict the response of a nonlinear system due to an arbitrary input. A numerical kernel identification technique, based on unit impulse responses, is presented and applied to a simple bilinear, single-input single-output (SISO) system. The linear kernel (unit impulse response) and the nonlinear second-order kernel of the system are numerically-identified and compared with the exact, analytically-defined and linear and second-order kernels. This kernel identification technique is then applied to the CAP-TSD (Computational Aeroelasticity Program-Transonic Small Disturbance) code for identification of the linear and second-order kernels of a NACA64A010 rectangular wing undergoing pitch at M = 0.5, M = 8.5 (transonic), and M = 0.93 (transonic). Results presented demonstrate the feasibility of this approach for use with nonlinear, unsteady aerodynamic responses.

Silva, Walter A.↗

A methodology for using nonlinear aerodynamics in aeroservoelastic analysis and design

A methodology is presented for using the Volterra-Wiener theory of nonlinear systems in aeroservoelastic (ASE) analyses and design. The theory is applied to the development of nonlinear aerodynamic response models that can be defined in state-space form and are, therefore, appropriate for use in modern control theory. The theory relies on the identification of nonlinear kernels that can be used to predict the response of a nonlinear system due to an arbitrary input. A numerical kernel identification technique, based on unit impulse responses, is presented and applied to a simple bilinear, single-input-single-output system. The linear kernel (unit impulse response) and the nonlinear second-order kernel of the system are numerically-identified and compared with the exact, analytically-defined linear and second-order kernels. This kernel identification technique is then applied to the CAP-TSD code for identification of the linear and second-order kernels of a NACA64A010 rectangular wing undergoing pitch at M = 0.5, M = 0.85 (transonic), and M = 0.93 (transonic). Results presented demonstrate the feasibility of this approach for use with nonlinear, unsteady aerodynamic responses.

Silva, Walter A.↗

Non-oscillatory and non-diffusive solution of convection problems by the iteratively reweighted least-squares finite element method

A comparative description is presented for the least-squares FEM (LSFEM) for 2D steady-state pure convection problems. In addition to exhibiting better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method, numerical convergence rates are obtained which show the LSFEM to be virtually optimal. The LSFEM is used as a framework for an iteratively reweighted LSFEM yielding nonoscillatory and nondiffusive solutions for problems with contact discontinuities; this method is shown to convect contact discontinuities without error when using triangular and bilinear elements.

Jiang, Bo-Nan↗

Programmer's manual for MMLE3, a general FORTRAN program for maximum likelihood parameter estimation

The MMLE3 is a maximum likelihood parameter estimation program capable of handling general bilinear dynamic equations of arbitrary order with measurement noise and/or state noise (process noise). The basic MMLE3 program is quite general and, therefore, applicable to a wide variety of problems. The basic program can interact with a set of user written problem specific routines to simplify the use of the program on specific systems. A set of user routines for the aircraft stability and control derivative estimation problem is provided with the program. The implementation of the program on specific computer systems is discussed. The structure of the program is diagrammed, and the function and operation of individual routines is described. Complete listings and reference maps of the routines are included on microfiche as a supplement. Four test cases are discussed; listings of the input cards and program output for the test cases are included on microfiche as a supplement.

Maine, R. E.↗

Interior Noise Reduction by Adaptive Feedback Vibration Control

The objective of this project is to investigate the possible use of adaptive digital filtering techniques in simultaneous, multiple-mode identification of the modal parameters of a vibrating structure in real-time. It is intended that the results obtained from this project will be used for state estimation needed in adaptive structural acoustics control. The work done in this project is basically an extension of the work on real-time single mode identification, which was performed successfully using a digital signal processor (DSP) at NASA, Langley. Initially, in this investigation the single mode identification work was duplicated on a different processor, namely the Texas Instruments TMS32OC40 DSP. The system identification results for the single mode case were very good. Then an algorithm for simultaneous two mode identification was developed and tested using analytical simulation. When it successfully performed the expected tasks, it was implemented in real-time on the DSP system to identify the first two modes of vibration of a cantilever aluminum beam. The results of the simultaneous two mode case were good but some problems were identified related to frequency warping and spurious mode identification. The frequency warping problem was found to be due to the bilinear transformation used in the algorithm to convert the system transfer function from the continuous-time domain to the discrete-time domain. An alternative approach was developed to rectify the problem. The spurious mode identification problem was found to be associated with high sampling rates. Noise in the signal is suspected to be the cause of this problem but further investigation will be needed to clarify the cause. For simultaneous identification of more than two modes, it was found that theoretically an adaptive digital filter can be designed to identify the required number of modes, but the algebra became very complex which made it impossible to implement in the DSP system used in this study. The on-line identification algorithm developed in this research will be useful in constructing a state estimator for feedback vibration control.

Lim, Tae W.↗

User's manual for MMLE3, a general FORTRAN program for maximum likelihood parameter estimation

A user's manual for the FORTRAN IV computer program MMLE3 is described. It is a maximum likelihood parameter estimation program capable of handling general bilinear dynamic equations of arbitrary order with measurement noise and/or state noise (process noise). The theory and use of the program is described. The basic MMLE3 program is quite general and, therefore, applicable to a wide variety of problems. The basic program can interact with a set of user written problem specific routines to simplify the use of the program on specific systems. A set of user routines for the aircraft stability and control derivative estimation problem is provided with the program.

Maine, R. E.↗

The Flux-integral Method for Multidimensional Convection and Diffusion

The flux-integral method is a procedure for constructing an explicit, single-step, forward-in-time, conservative, control volume update of the unsteady, multidimensional convection-diffusion equation. The convective plus diffusive flux at each face of a control-volume cell is estimated by integrating the transported variable and its face-normal derivative over the volume swept out by the convecting velocity field. This yields a unique description of the fluxes, whereas other conservative methods rely on nonunique, arbitrary pseudoflux-difference splitting procedures. The accuracy of the resulting scheme depends on the form of the subcell interpolation assumed, given cell-average data. Cellwise constant behavior results in a (very artificially diffusive) first-order convection scheme. Second-order convection-diffusion schemes correspond to cellwise linear (or bilinear) subcell interpolation. Cellwise quadratic subcell interpolants generate a highly accurate convection-diffusion scheme with excellent phase accuracy. Under constant-coefficient conditions, this is a uniformly third-order polynomial interpolation algorithm (UTOPIA).

Leonard, B. P.↗