Engineering PapersSearch

Engineering topics

Desmarais, R. N.

Publications and source records attributed to Desmarais, R. N..

At least 19 records

A method for computing the kernel of the downwash integral equation for arbitrary complex frequencies

For the design of active controls to stabilize flight vehicles, which requires the use of unsteady aerodynamics that are valid for arbitrary complex frequencies, algorithms are derived for evaluating the nonelementary part of the kernel of the integral equation that relates unsteady pressure to downwash. This part of the kernel is separated into an infinite limit integral that is evaluated using Bessel and Struve functions and into a finite limit integral that is expanded in series and integrated termwise in closed form. The developed series expansions gave reliable answers for all complex reduced frequencies and executed faster than exponential approximations for many pressure stations.

Desmarais, R. N.

Generalization of the subsonic kernel function in the s-plane, with applications to flutter analysis

A generalized subsonic unsteady aerodynamic kernel function, valid for both growing and decaying oscillatory motions, is developed and applied in a modified flutter analysis computer program to solve the boundaries of constant damping ratio as well as the flutter boundary. Rates of change of damping ratios with respect to dynamic pressure near flutter are substantially lower from the generalized-kernel-function calculations than from the conventional velocity-damping (V-g) calculation. A rational function approximation for aerodynamic forces used in control theory for s-plane analysis gave rather good agreement with kernel-function results, except for strongly damped motion at combinations of high (subsonic) Mach number and reduced frequency.

Cunningham, H. J.

Planar-Wing Flutter Analysis

Flutter Analysis System, FAST, is group of five programs that perform flutter analysis of single planar wing. Each program performs certain portions of flutter analysis and may be run sequentially or individually. Modular approach is very versatile and flexible. FAST uses natural vibration modes as input data and performs conventional V-g solution.

Desmarais, R. N.

Steady, Oscillatory, and Unsteady Subsonic and Supersonic Aerodynamics

Computer program SOUSSA-P (Steady, Oscillatory, and Unsteady Subsonic and Supersonic Aerodynamics--Production Version) accurately and efficiently evaluates steady and unsteady aerodynamic loads on aircraft having arbitrary shapes and motions, including structural deformations.

Desmarais, R. N.

Subsonic aerodynamic and flutter characteristics of several wings calculated by the SOUSSA P1.1 panel method

The SOUSSA (steady, oscillatory, and unsteady subsonic and supersonic aerodynamics) program is the computational implementation of a general potential flow analysis (by the Green's function method) that can generate pressure distributions on complete aircraft having arbitrary shapes, motions and deformations. Some applications of the initial release version of this program to several wings in steady and oscillatory motion, including flutter are presented. The results are validated by comparisons with other calculations and experiments. Experiences in using the program as well as some recent improvements are described.

Yates, E. C., Jr.

Subsonic aerodynamic and flutter characteristics of several wings calculated by the SOUSSA P1.1 panel method

Several applications of the steady, oscillatory, and unsteady subsonic and supersonic aerodynamics (SOUSSA) computer program to wings with steady and oscillatory motion, including flutter, are discussed. The program employs a generalized Green's function to the full, time-dependent potential-flow equation to obtain an integral equation for the velocity potential at any point in a flow, even points on a body or whole bodies in a flow. Aerodynamic calculations are provided for two rectangular wings, a clipped-tip delta wing, and two swept wings with and without a fuselage. The number and distribution of the finite element panels are varied in order to demonstrate the convergence of the results. The results are shown to be close to those of lifting-surface theory, and further applications with bodies having deformities, arbitrary shapes, motions, and deformations are indicated.

Yates, E. C., Jr.

An accurate method for evaluating the kernel of the integral equation relating lift to downwash in unsteady potential flow

The method is capable of generating approximations of arbitrary accuracy. It is based on approximating the algebraic part of the nonelementary integrals in the kernel by exponential functions and then integrating termwise. The exponent spacing in the approximation is a geometric sequence. The coefficients and exponent multiplier of the exponential approximation are computed by least squares so the method is completely automated. Exponential approximates generated in this manner are two orders of magnitude more accurate than the exponential approximation that is currently most often used for this purpose. The method can be used to generate approximations to attain any desired trade-off between accuracy and computing cost.

Desmarais, R. N.

An accurate and efficient method for evaluating the kernel of the integral equation relating pressure to normalwash in unsteady potential flow

This paper describes an accurate economical method for generating approximations to the kernel of the integral equation relating unsteady pressure to normalwash in nonplanar flow. The method is capable of generating approximations of arbitrary accuracy. It is based on approximating the algebraic part of the non elementary integrals in the kernel by exponential approximations and then integrating termwise. The exponent spacing in the approximation is a geometric sequence. The coefficients and exponent multiplier of the exponential approximation are computed by least squares so the method is completely automated. Exponential approximates generated in this manner are two orders of magnitude more accurate than the exponential approximation that is currently most often used for this purpose. Coefficients for 8, 12, 24, and 72 term approximations are tabulated in the report. Also, since the method is automated, it can be used to generate approximations to attain any desired trade-off between accuracy and computing cost.

Desmarais, R. N.

Prediction of transonic flutter for a supercritical wing by modified strip analysis and comparison with experiment

Use of a supercritical airfoil can adversely affect wing flutter speeds in the transonic range. As adequate theories for three dimensional unsteady transonic flow are not yet available, the modified strip analysis was used to predict the transonic flutter boundary for the supercritical wing. The steady state spanwise distributions of section lift curve slope and aerodynamic center, required as input for the flutter calculations, were obtained from pressure distributions. The calculated flutter boundary is in agreement with experiment in the subsonic range. In the transonic range, a transonic bucket is calculated which closely resembles the experimental one with regard to both shape and depth, but it occurs at about 0.04 Mach number lower than the experimental one.

Yates, E. C., Jr.

Prediction of transonic flutter for a supercritical wing by modified strip analysis and comparison with experiment

The experiments of Farmer, Hanson, and Wynne (1976) showed that use of a supercritical airfoil can adversely affect wing flutter speeds in the transonic range. Inasmuch as adequate theories for three-dimensional unsteady transonic flow are not yet available, the modified strip analysis published by Yates in 1958 has been used to predict the transonic flutter boundary for the supercritical wing tested by Farmer, Hanson, and Wynne. The steady-state spanwise distributions of section lift-curve slope and aerodynamic center, required as input for the flutter calculations, were obtained from pressure distributions measured by Harris (1972). The calculated flutter boundary is in excellent agreement with experiment in the subsonic range. In the transonic range a 'transonic bucket' is calculated which closely resembles the experimental one with regard to both shape and depth, but it occurs at about 0.04 Mach number lower than the experimental one.

Yates, E. C., Jr.

Wing/store flutter with nonlinear pylon stiffness

Recent wind tunnel tests and analytical studies show that a store mounted on a pylon with soft pitch stiffness provides substantial increase in flutter speed of fighter aircraft and reduces dependency of flutter on mass and inertia of the store. This concept, termed the decoupler pylon, utilizes a low frequency control system to maintain pitch alignment of the store during maneuvers and changing flight conditions. Under rapidly changing transient loads, however, the alignment control system may allow the store to momentarily bottom against a relatively stiff backup structure in which case the pylon stiffness acts as a hardening nonlinear spring. Such structural nonlinearities are known to affect not only the flutter speed but also the basic behavior of the instability. The influence of pylon stiffness nonlinearities or the flutter characteristics of wing mounted external stores is examined.

Desmarais, R. N.

Wing/store flutter with nonlinear pylon stiffness

Recent wind tunnel tests and analytical studies show that a store mounted on a pylon with 'soft' pitch stiffness provides substantial increase in flutter speed of fighter aircraft and reduces dependency of flutter on mass and inertia of the store. This concept, termed the decoupler pylon, utilizes a low-frequency control system to maintain pitch alignment of the store during maneuvers and changing flight conditions. Under rapidly changing transient loads, however, the alignment control system may allow the store to momentarily bottom against a relatively stiff backup structure in which case the pylon stiffness acts as a hardening nonlinear spring. Such structural nonlinearities are known to affect not only the flutter speed but also the basic behavior of the instability. This paper examines the influence of pylon stiffness nonlinearities on the flutter characteristics of wing-mounted external stores.

Desmarais, R. N.

Accurate numerical evaluation of modified Struve functions occurring in unsteady aerodynamics

A method is developed for computing the modified Struve functions that occur in unsteady aerodynamics. The method uses a rational approximation supplemented by an asymptotic series for large argument. Simple recursive formulas for generating the coefficients are derived. The method is capable of generating results of arbitrary accuracy. It can also be used for complex argument and order. For greater computing speed, a method is presented that uses the rational and asymptotic approximations to generate Chebyshev coefficients.

Desmarais, R. N.

User's Guide for a Modular Flutter Analysis Software System (Fast Version 1.0)

The use and operation of a group of computer programs to perform a flutter analysis of a single planar wing are described. This system of programs is called FAST for Flutter Analysis System, and consists of five programs. Each program performs certain portions of a flutter analysis and can be run sequentially as a job step or individually. FAST uses natural vibration modes as input data and performs a conventional V-g type of solution. The unsteady aerodynamics programs in FAST are based on the subsonic kernel function lifting-surface theory although other aerodynamic programs can be used. Application of the programs is illustrated by a sample case of a complete flutter calculation that exercises each program.

Desmarais, R. N.

Curve fitting of aeroelastic transient response data with exponential functions

The extraction of frequency, damping, amplitude, and phase information from unforced transient response data is considered. These quantities are obtained from the parameters determined by fitting the digitized time-history data in a least-squares sense with complex exponential functions. The highlights of the method are described, and the results of several test cases are presented. The effects of noise are considered both by using analytical examples with random noise and by estimating the standard deviation of the parameters from maximum-likelihood theory.

Bennett, R. M.

Programs for computing abscissas and weights for classical and nonclassical Gaussian quadrature formulas

Computer programs for computing Gaussian quadrature abscissas and weights are described. For the classical case the programs use Laguerre iteration to compute abscissas as zeros of orthogonal polynomials. The polynomials are evaluated from known recursion coefficients. The nonclassical case is handled similarly except that the recursion coefficients are computed by numerical integration. A sample problem, with input and output, is presented to illustrate the use of the programs. It computes the quadrature abscissas and weights associated with the weight function over the interval (0,1) for quadrature orders from 16 to 96 in increments of 8.

Desmarais, R. N.