Engineering PapersSearch

SEARCH · Engineering Papers

Results for “QUADRATIC EQUATION”

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 19 records

Optimization of structures to satisfy a flutter velocity constraint by use of quadratic equation fitting

Using the first and the second derivative of flutter velocity with respect to the parameters, the velocity hypersurface is made quadratic. This greatly simplifies the numerical procedure developed for determining the values of the design parameters such that a specified flutter velocity constraint is satisfied and the total structural mass is near a relative minimum. A search procedure is presented utilizing two gradient search methods and a gradient projection method. The procedure is applied to the design of a box beam, using finite-element representation. The results indicate that the procedure developed yields substantial design improvement satisfying the specified constraint and does converge to near a local optimum.

Motiwalla, S. K.

Numerical solution of quadratic matrix equations for free vibration analysis of structures

This paper is concerned with the efficient and accurate solution of the eigenvalue problem represented by quadratic matrix equations. Such matrix forms are obtained in connection with the free vibration analysis of structures, discretized by finite 'dynamic' elements, resulting in frequency-dependent stiffness and inertia matrices. The paper presents a new numerical solution procedure of the quadratic matrix equations, based on a combined Sturm sequence and inverse iteration technique enabling economical and accurate determination of a few required eigenvalues and associated vectors. An alternative procedure based on a simultaneous iteration procedure is also described when only the first few modes are the usual requirement. The employment of finite dynamic elements in conjunction with the presently developed eigenvalue routines results in a most significant economy in the dynamic analysis of structures.

Gupta, K. K.

Solution of quadratic matrix equations for free vibration analysis of structures.

An efficient digital computer procedure and the related numerical algorithm are presented herein for the solution of quadratic matrix equations associated with free vibration analysis of structures. Such a procedure enables accurate and economical analysis of natural frequencies and associated modes of discretized structures. The numerically stable algorithm is based on the Sturm sequence method, which fully exploits the banded form of associated stiffness and mass matrices. The related computer program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be substantially more accurate and economical than other existing procedures of such analysis. Numerical examples are presented for two structures - a cantilever beam and a semicircular arch.

Gupta, K. K.

Inverting x,y grid coordinates to obtain latitude and longitude in the vanderGrinten projection

The latitude and longitude of a point on the Earth's surface are found from its x,y grid coordinates in the vanderGrinten projection. The latitude is a solution of a cubic equation and the longitude a solution of a quadratic equation. Also, the x,y grid coordinates of a point on the Earth's surface can be found if its latitude and longitude are known by solving two simultaneous quadratic equations.

Rubincam, D. P.

Ground-Support Algorithms for Simulation, Processing, and Calibration of Barnes Static Earth Sensor Measurements: Applications to Tropical Rainfall Measuring Mission Observatory

New algorithms are described covering the simulation, processing, and calibration of penetration angles of the Barnes static Earth sensor assembly (SESA) as implemented in the Goddard Space Flight Center Flight Dynamics Division ground support system for the Tropical Rainfall Measuring Mission (TRMM) Observatory. The new treatment involves a detailed analysis of the measurements by individual quadrants. It is shown that, to a good approximation, individual quadrant misalignments can be treated simply as penetration angle biases. Simple formulas suitable for real-time applications are introduced for computing quadrant-dependent effects. The simulator generates penetration angles by solving a quadratic equation with coefficients uniquely determined by the spacecraft's position and the quadrant's orientation in GeoCentric Inertial (GCI) coordinates. Measurement processing for attitude determination is based on linearized equations obtained by expanding the coefficients of the aforementioned quadratic equation as a Taylor series in both the Earth oblateness coefficient (alpha approx. 1/150) and the angle between the pointing axis and the geodetic nadir vector. A simple formula relating a measured value of the penetration angle to the deviation of the Earth-pointed axis from the geodetic nadir vector is derived. It is shown that even near the very edge of the quadrant's Field Of View (FOV), attitude errors resulting from quadratic effects are a few hundredths of a degree, which is small compared to the attitude determination accuracy requirement (0.18 degree, 3 sigma) of TRMM. Calibration of SESA measurements is complicated by a first-order filtering used in the TRMM onboard algorithm to compute penetration angles from raw voltages. A simple calibration scheme is introduced where these complications are avoided by treating penetration angles as the primary raw measurements, which are adjusted using biases and scale factors. In addition to three misalignment parameters, the calibration state vector contains only two average penetration angle biases (one per each pair of opposite quadrants) since, because of the very narrow sensor FOV (+/- 2.6 degrees), differences between biases of the penetration angles measured by opposite quadrants cannot be distinguished from roll and pitch sensor misalignments. After calibration, the estimated misalignments and average penetration angle biases are converted to the four penetration angle biases and to the yaw misalignment angle. The resultant biases and the estimated scale factors are finally used to update the coefficients necessary for onboard computations of penetration angles from measured voltages.

Natanson, G. A.

Improvements in Modeling 90 degree Bleed Holes for Supersonic Inlets

The modeling of porous bleed regions as boundary conditions in computational fluid dynamics (CFD) simulations of supersonic inlet flows has been improved through a scaling of sonic flow coefficient data for 90deg bleed holes. The scaling removed the Mach number as a factor in computing the sonic flow coefficient and allowed the data to be fitted with a quadratic equation, with the only factor being the ratio of the plenum static pressure to the surface static pressure. The implementation of the bleed model into the Wind-US CFD flow solver was simplified by no longer requiring the evaluation of the flow properties at the boundary-layer edge. The quadratic equation can be extrapolated to allow the modeling of small amounts of blowing, which can exist when recirculation of the bleed flow occurs within the bleed region. The improved accuracy of the bleed model was demonstrated through CFD simulations of bleed regions on a flat plate in supersonic flow with and without an impinging oblique shock. The bleed model demonstrated good agreement with experimental data and three-dimensional CFD simulations of bleed holes.

Slater, John W.

An analytic approach to two-fixed-impulse transfers between Keplerian orbits

Solutions are obtained for the two-impulse transfer of a vehicle between arbitrary inclined orbits in an inverse-square force field with the restriction that the magnitude of each of the two impulses has a fixed preassigned value. The two magnitudes need not be equal. The equations for the conservation of angular momentum and energy are augmented by the Laplace integral; these equations establish linear relationships between several of the variables. This set of linear equations and one of two quadratic equations constitute the analytically tractable part of the solution. The remaining part, consisting of finding the zeros (if they exist) of a single trigonometric function of one variable, is solved using numerical methods. Explicit lower bounds on each of the magnitudes of the two impulses are obtained by requiring the solution be real. Graphical results are presented to illustrate the solution.

Jezewski, D. J.

Effects of magnetometer calibration and maneuvers on accuracies of magnetometer-only attitude-and-rate determination

Two different algorithms - a deterministic magnetic-field-only algorithm and a Kalman filter for gyroless spacecraft - are used to estimate the attitude and rates of the Rossi X-Ray Timing Explorer (RXTE) using only measurements from a three-axis magnetometer. The performance of these algorithms is examined using in-flight data from various scenarios. In particular, significant enhancements in accuracies are observed when' the telemetered magnetometer data are accurately calibrated using a recently developed calibration algorithm. Interesting features observed in these studies of the inertial-pointing RXTE include a remarkable sensitivity of the filter to the numerical values of the noise parameters and relatively long convergence time spans. By analogy, the accuracy of the deterministic scheme is noticeably lower as a result of reduced rates of change of the body-fixed geomagnetic field. Preliminary results show the filter-per-axis attitude accuracies ranging between 0.1 and 0.5 deg and rate accuracies between 0.001 deg/sec and 0.005 deg./sec, whereas the deterministic method needs a more sophisticated techniques for smoothing time derivatives of the measured geomagnetic field to clearly distinguish both attitude and rate solutions from the numerical noise. Also included is a new theoretical development in the deterministic algorithm: the transformation of a transcendental equation in the original theory into an 8th-order polynomial equation. It is shown that this 8th-order polynomial reduces to quadratic equations in the two limiting cases-infinitely high wheel momentum, and constant rates-discussed in previous publications.

Challa, M.

A curve fitting method for solving the flutter equation

A curve fitting approach was developed to solve the flutter equation for the critical flutter velocity. The psi versus nu curves are approximated by cubic and quadratic equations. The curve fitting technique utilized the first and second derivatives of psi with respect to nu. The method was tested for two structures, one structure being six times the total mass of the other structure. The algorithm never showed any tendency to diverge from the solution. The average time for the computation of a flutter velocity was 3.91 seconds on an IBM Model 50 computer for an accuracy of five per cent. For values of nu close to the critical root of the flutter equation the algorithm converged on the first attempt. The maximum number of iterations for convergence to the critical flutter velocity was five with an assumed value of nu relatively distant from the actual crossover.

Cooper, J. L.

An experimental and analytical investigation of axisymmetric diffusers

A finite difference computer program for turbulent compressible flow was used to establish the performance of several diffuser shapes for experimental testing. The diffusers were designed to have a linear change in Mach number, a linear change in pressure, or a curvature fitted by a quadratic equation. Testing was performed with M = 0.1 to 0.9 with and without boundary layer bleed. Above M = 0.6, data were obtained with a normal shock upstream of the diffuser entrance. Peak static pressure recovery occurred with a diffuser inlet M0.75. The quadratic diffuser yielded the highest total pressure recovery.

Povinelli, L. A.

An experimental and analytical investigation of axisymmetric diffusers

A finite difference computer program for turbulent compressible flow was used to establish the performance of several diffuser shapes for experimental testing. The diffusers were designed to have a linear change in Mach number, a linear change in pressure or a curvature fitted by a quadratic equation. Testing was performed with M = 0.1 to 0.9 with and without boundary layer bleed. Above M = 0.6, data were obtained with a normal shock upstream of the diffuser entrance. Peak static pressure recovery occurred with a diffuser inlet M-0.75. The quadratic diffuser yielded the highest total pressure recovery.

Povinelli, L. A.

Matrix quadratic solutions.

Matrix quadratic equation solution derivation applied in finding steady state solutions of Riccati differential equations with constant coefficients

OPTIMAL CONTROL

Numerical solution of the state dependent noise problem

A numerical technique is given for solving the matrix quadratic equation that arises in the optimal stationary control of linear systems with state (and/or control) dependent noise. The technique exploits fully existing, efficient algorithms for the matrix Lyapunov and Ricatti equations. The computational requirements are discussed, with an associated example.

Kleinman, D. L.