Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear programming”

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 1,099 records · Page 61

Design of suboptimal regulators for nonlinear systems

An optimal feedback control law is preferred for the regulation of a deterministic nonlinear system. In this paper, a practical, iterative design method leading to a sequence of suboptimal control laws with successively improved performance is presented. The design method requires the determination of an upper bound to the performance of each successive control law. This is obtained by solving a partial differential inequality by means of a linear programming technique. Robustness properties and the application of the design method to the control of a robot manipulator arm are also presented.

Balaram, J.↗

Application of triggered lightning numerical models to the F106B and extension to other aircraft

The goal of the F106B Thunderstorm Research Program is to characterize the lightning environment for aircraft in flight. This report describes the application of numerical electromagnetic models to this problem. Topics include: (1) Extensive application of linear triggered lightning to F106B data; (2) Electrostatic analysis of F106B field mill data; (3) Application of subgrid modeling to F106B nose region, including both static and nonlinear models; (4) Extension of F106B results to other aircraft of varying sizes and shapes; and (5) Application of nonlinear model to interaction of F106B with lightning leader-return stroke event.

Ng, Poh H.↗

Structural Reliability Using Probability Density Estimation Methods Within NESSUS

A reliability analysis studies a mathematical model of a physical system taking into account uncertainties of design variables and common results are estimations of a response density, which also implies estimations of its parameters. Some common density parameters include the mean value, the standard deviation, and specific percentile(s) of the response, which are measures of central tendency, variation, and probability regions, respectively. Reliability analyses are important since the results can lead to different designs by calculating the probability of observing safe responses in each of the proposed designs. All of this is done at the expense of added computational time as compared to a single deterministic analysis which will result in one value of the response out of many that make up the density of the response. Sampling methods, such as monte carlo (MC) and latin hypercube sampling (LHS), can be used to perform reliability analyses and can compute nonlinear response density parameters even if the response is dependent on many random variables. Hence, both methods are very robust; however, they are computationally expensive to use in the estimation of the response density parameters. Both methods are 2 of 13 stochastic methods that are contained within the Numerical Evaluation of Stochastic Structures Under Stress (NESSUS) program. NESSUS is a probabilistic finite element analysis (FEA) program that was developed through funding from NASA Glenn Research Center (GRC). It has the additional capability of being linked to other analysis programs; therefore, probabilistic fluid dynamics, fracture mechanics, and heat transfer are only a few of what is possible with this software. The LHS method is the newest addition to the stochastic methods within NESSUS. Part of this work was to enhance NESSUS with the LHS method. The new LHS module is complete, has been successfully integrated with NESSUS, and been used to study four different test cases that have been proposed by the Society of Automotive Engineers (SAE). The test cases compare different probabilistic methods within NESSUS because it is important that a user can have confidence that estimates of stochastic parameters of a response will be within an acceptable error limit. For each response, the mean, standard deviation, and 0.99 percentile, are repeatedly estimated which allows confidence statements to be made for each parameter estimated, and for each method. Thus, the ability of several stochastic methods to efficiently and accurately estimate density parameters is compared using four valid test cases. While all of the reliability methods used performed quite well, for the new LHS module within NESSUS it was found that it had a lower estimation error than MC when they were used to estimate the mean, standard deviation, and 0.99 percentile of the four different stochastic responses. Also, LHS required a smaller amount of calculations to obtain low error answers with a high amount of confidence than MC. It can therefore be stated that NESSUS is an important reliability tool that has a variety of sound probabilistic methods a user can employ and the newest LHS module is a valuable new enhancement of the program.

Chamis, Chrisos C.↗

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗

Design and evaluation of a sensor fail-operational control system for a digitally controlled turbofan engine

A self-learning, sensor fail-operational, control system for the TF30-P-3 afterburning turbofan engine was designed and evaluated. The sensor fail-operational control system includes a digital computer program designed to operate in conjunction with the standard TF30-P-3 bill-of-materials control. Four engine measurements and two compressor face measurements are tested. If any engine measurements are found to have failed, they are replaced by values synthesized from computer-stored information. The control system was evaluated by using a realtime, nonlinear, hybrid computer engine simulation at sea level static condition, at a typical cruise condition, and at several extreme flight conditions. Results indicate that the addition of such a system can improve the reliability of an engine digital control system.

Hrach, F. J.↗

Comprehensive missile aerodynamics programs for preliminary design

Two different classes of missile aeroprediction programs have been recently developed. The first class of programs provides rapid engineering predictions and includes MISSILE1 and MISSILE2 applicable to missile configurations with axisymmetric bodies. The second class of programs consists of the DEMON series, including a simplified version NSWCDM, designed to calculate detailed loadings acting on supersonic missiles which may have non-circular body cross sections. Both classes account for high angles of attack and track vortices from canard or wing section to the tail section. Extensive comparisons with experimental data are presented including nonlinear effects of canard control.

Dillenius, M. F. E.↗

Unified aeroacoustics analysis for high speed turboprop aerodynamics and noise. Volume 3: Application of theory for blade loading, wakes, noise, and wing shielding

Results of the program for the generation of a computer prediction code for noise of advanced single rotation, turboprops (prop-fans) such as the SR3 model are presented. The code is based on a linearized theory developed at Hamilton Standard in which aerodynamics and acoustics are treated as a unified process. Both steady and unsteady blade loading are treated. Capabilities include prediction of steady airload distributions and associated aerodynamic performance, unsteady blade pressure response to gust interaction or blade vibration, noise fields associated with thickness and steady and unsteady loading, and wake velocity fields associated with steady loading. The code was developed on the Hamilton Standard IBM computer and has now been installed on the Cray XMP at NASA-Lewis. The work had its genesis in the frequency domain acoustic theory developed at Hamilton Standard in the late 1970s. It was found that the method used for near field noise predictions could be adapted as a lifting surface theory for aerodynamic work via the pressure potential technique that was used for both wings and ducted turbomachinery. In the first realization of the theory for propellers, the blade loading was represented in a quasi-vortex lattice form. This was upgraded to true lifting surface loading. Originally, it was believed that a purely linear approach for both aerodynamics and noise would be adequate. However, two sources of nonlinearity in the steady aerodynamics became apparent and were found to be a significant factor at takeoff conditions. The first is related to the fact that the steady axial induced velocity may be of the same order of magnitude as the flight speed and the second is the formation of leading edge vortices which increases lift and redistribute loading. Discovery and properties of prop-fan leading edge vortices were reported in two papers. The Unified AeroAcoustic Program (UAAP) capabilites are demonstrated and the theory verified by comparison with the predictions with data from tests at NASA-Lewis. Steady aerodyanmic performance, unsteady blade loading, wakes, noise, and wing and boundary layer shielding are examined.

Hanson, D. B.↗

Substructure analysis using NICE/SPAR and applications of force to linear and nonlinear structures

Parallel computing studies are presented for a variety of structural analysis problems. Included are the substructure planar analysis of rectangular panels with and without a hole, the static analysis of space mast, using NICE/SPAR and FORCE, and substructure analysis of plane rigid-jointed frames using FORCE. The computations are carried out on the Flex/32 MultiComputer using one to eighteen processors. The NICE/SPAR runstream samples are documented for the panel problem. For the substructure analysis of plane frames, a computer program is developed to demonstrate the effectiveness of a substructuring technique when FORCE is enforced. Ongoing research activities for an elasto-plastic stability analysis problem using FORCE, and stability analysis of the focus problem using NICE/SPAR are briefly summarized. Speedup curves for the panel, the mast, and the frame problems provide a basic understanding of the effectiveness of parallel computing procedures utilized or developed, within the domain of the parameters considered. Although the speedup curves obtained exhibit various levels of computational efficiency, they clearly demonstrate the excellent promise which parallel computing holds for the structural analysis problem. Source code is given for the elasto-plastic stability problem and the FORCE program.

Razzaq, Zia↗

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.↗

Development of welded metal bellows having minimum effective diameter change

A program of analysis, design, and fabrication was conducted to develop welded metal bellows having a minimum change in effective diameter for cryogenic turbomachinery face seal applications. Linear analysis of the principle types of bellows provided identification of concepts capable of meeting basic operation requirements. For the 6-inch (.152 m) mean diameter, 1.5-inch free length bellows studied, nonlinear analysis showed that opposed and nested toroidal bellows plates stiffened by means of alternating stiffener rings were capable of maintaining constant effective diameter within 0.3% and 0.1% respectively under the operating conditions of interest. Changes in effective diameter were due principally to bellows axial deflection with pressure differential having a lesser influence. Fabrication problems associated with joining the thin bellows plates to the relatively heavy stiffener rings were encountered and precluded assembly and testing of a bellows core. Fabrication problems are summarized and recommended fabrication methods for future effort are presented.

Henschel, J. K.↗

Development of an advanced pitch active control system and a reduced area horizontal tail for a wide-body jet aircraft

The development of an advanced pitch active control system (PACS) and a reduced area horizontal tail for a wide-body jet transport (L-1011) with a flying horizontal stabilizer is discussed. The advanced PACS control law design objectives were to provide satisfactory handling qualities for aft c.g. flight conditions to negative static stability margins of 10 percent and to provide good maneuver control column force gradients for nonlinear stability flight conditions. Validity of the control laws were demonstrated by piloted flight simulation tests on the NASA Langley Visual Motion Simulator. Satisfactory handling qualities were actually demonstrated to a negative 20 percent static stability margin. The PACS control laws were mechanized to provide the system architecture that would be suitable for an L-1011 flight test program to a negative stability margin of 3 percent which represents the aft c.g. limits of the aircraft. Reduced area horizontal tail designs of 30 and 38 percent with respect to the L-1011 standard tail were designed, fabricated and wind tunnel tested. Drag reductions and weight savings of the 30 percent smaller tail would provide an L/D benefit of about 2% and the 38% small tail L/D benefit would be about 3 percent. However, forward c.g. limitations would have to be imposed on the aircraft because the maximum horizontal tail lift goal was not achieved and sufficient aircraft nose-up control authority was not available. This limitation would not be required for a properly designed new aircraft.

Guinn, Wiley A.↗

Application of decomposition techniques to the preliminary design of a transport aircraft

A multifaceted decomposition of a nonlinear constrained optimization problem describing the preliminary design process for a transport aircraft has been made. Flight dynamics, flexible aircraft loads and deformations, and preliminary structural design subproblems appear prominently in the decomposition. The use of design process decomposition for scheduling design projects, a new system integration approach to configuration control, and the application of object-centered programming to a new generation of design tools are discussed.

Rogan, J. E.↗

Analysis of structures with rotating, flexible substructures applied to rotorcraft aeroelasticity in GRASP

Application to the General Rotorcraft Aeromechanical Stability Program (GRASP) of new methodology for structural dynamic analysis, including substructuring, frames of reference, nodes, finite elements, and constraints, is discussed. The structure is decomposed into a hierarchy of substructures, and discrete relative motion between substructures is analyzed exactly. The finite element method is used to treat deformation of continua, and the library of finite elements includes a nonlinear beam element incorporating aeroelastic effects. Analytical bases for the aeroelastic beam element and the screw constraint are considered, and the important role of geometric stiffness in the formulation is shown.

Hodges, Dewey H.↗

Using NASTRAN To Analyze Vibrations Of Rotor Blades

Report gives information on use of NASTRAN computer program in finite-element analysis of rotating flexible blades like in compressors and on turboprop engines. Predicts steady-state components of deflections and stresses under centrifugal forces, generates data for plots of natural frequency versus rotational speed, and provides vibration-mode data for calculations of flutter. Describes use of NASTRAN solution sequence 64 for geometrical nonlinear analysis and solution sequence 63 for determination of frequencies and vibrational-mode shapes. Includes sample problem with NASTRAN input data. Emphasizes key factors in analysis of rotating blades, such as setting angle and centrifugal softening effects. Combined analyses of solution sequences 64 and 63 reduces computer time and number of output listings, in comparison with separate analyses. In central-processing-unit time cut in half.

Lawrence, Charles↗

Investigation of difficult component effects on FEM vibration prediction for the AH-1G helicopter

Under the NASA-sponsored Design Analysis Methods for Vibrations program, a series of ground vibration tests and NASTRAN finite element model correlations were conducted on the Bell AH-1G helicopter gunship to investigate the effects of difficult components on the vibration response of the airframe. Secondary structure and damping were found to have significant effects on the frequency response of the airframe above 15 Hz. The nonlinear effects of thrust stiffening and elastomeric mounts on the low-frequency pylon modes below the main rotor were also significant.

Dompka, Robert V.↗

Thermoviscoelastic characterization and prediction of Kevlar/epoxy composite laminates

The thermoviscoelastic characterization of Kevlar 49/Fiberite 7714A epoxy composite lamina and the development of a numerical procedure to predict the viscoelastic response of any general laminate constructed from the same material were studied. The four orthotropic material properties, S sub 11, S sub 12, S sub 22, and S sub 66, were characterized by 20 minute static creep tests on unidirectional (0) sub 8, (10) sub 8, and (90) sub 16 lamina specimens. The Time-Temperature Superposition-Principle (TTSP) was used successfully to accelerate the characterization process. A nonlinear constitutive model was developed to describe the stress dependent viscoelastic response for each of the material properties. A numerical procedure to predict long term laminate properties from lamina properties (obtained experimentally) was developed. Numerical instabilities and time constraints associated with viscoelastic numerical techniques were discussed and solved. The numerical procedure was incorporated into a user friendly microcomputer program called Viscoelastic Composite Analysis Program (VCAP), which is available for IBM PC type computers. The program was designed for ease of use. The final phase involved testing actual laminates constructed from the characterized material, Kevlar/epoxy, at various temperatures and load level for 4 to 5 weeks. These results were compared with the VCAP program predictions to verify the testing procedure and to check the numerical procedure used in the program. The actual tests and predictions agreed for all test cases which included 1, 2, 3, and 4 fiber direction laminates.

Gramoll, K. C.↗

Crack Growth Simulation and Residual Strength Prediction in Airplane Fuselages

The objectives were to create a capability to simulate curvilinear crack growth and ductile tearing in aircraft fuselages subjected to widespread fatigue damage and to validate with tests. Analysis methodology and software program (FRANC3D/STAGS) developed herein allows engineers to maintain aging aircraft economically, while insuring continuous airworthiness, and to design more damage-tolerant aircraft for the next generation. Simulations of crack growth in fuselages were described. The crack tip opening angle (CTOA) fracture criterion, obtained from laboratory tests, was used to predict fracture behavior of fuselage panel tests. Geometrically nonlinear, elastic-plastic, thin shell finite element crack growth analyses were conducted. Comparisons of stress distributions, multiple stable crack growth history, and residual strength between measured and predicted results were made to assess the validity of the methodology. Incorporation of residual plastic deformations and tear strap failure was essential for accurate residual strength predictions. Issue related to predicting crack trajectory in fuselages were also discussed. A directional criterion, including T-stress and fracture toughness orthotropy, was developed. Curvilinear crack growth was simulated in coupon and fuselage panel tests. Both T-stress and fracture toughness orthotropy were essential to predict the observed crack paths. Flapping of fuselages were predicted. Measured and predicted results agreed reasonable well.

Chen, Chuin-Shan↗

Recent research on aerodynamic characteristics of fighter configurations during spins

The NASA Langley Research Center is currently conducting a stall/spin research program to define fighter aerodynamics applicable during developed spins and to develop analytical methods to use such measured aerodynamics for theoretically calculating spin motions. Some static, forced-oscillation and continuous rotation aerodynamic data have been measured for several current fighter models at developed spin angles of attack. The paper discusses these aerodynamic data and illustrates both the extremely nonlinear dependence of such data on several variables and the correlation that exists between the three types of measured aerodynamics. The current analytical methods for using these aerodynamics to calculate spin motions are discussed and correlated with experimentally obtained spins.

Anglin, E. L.↗