Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear Dynamic 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 235 records · Page 13

Automating the parallel processing of fluid and structural dynamics calculations

The NASA Lewis Research Center is actively involved in the development of expert system technology to assist users in applying parallel processing to computational fluid and structural dynamic analysis. The goal of this effort is to eliminate the necessity for the physical scientist to become a computer scientist in order to effectively use the computer as a research tool. Programming and operating software utilities have previously been developed to solve systems of ordinary nonlinear differential equations on parallel scalar processors. Current efforts are aimed at extending these capabilities to systems of partial differential equations, that describe the complex behavior of fluids and structures within aerospace propulsion systems. This paper presents some important considerations in the redesign, in particular, the need for algorithms and software utilities that can automatically identify data flow patterns in the application program and partition and allocate calculations to the parallel processors. A library-oriented multiprocessing concept for integrating the hardware and software functions is described.

Arpasi, Dale J.↗

Automating the parallel processing of fluid and structural dynamics calculations

The NASA Lewis Research Center is actively involved in the development of expert system technology to assist users in applying parallel processing to computational fluid and structural dynamic analysis. The goal of this effort is to eliminate the necessity for the physical scientist to become a computer scientist in order to effectively use the computer as a research tool. Programming and operating software utilities have previously been developed to solve systems of ordinary nonlinear differential equations on parallel scalar processors. Current efforts are aimed at extending these capabilties to systems of partial differential equations, that describe the complex behavior of fluids and structures within aerospace propulsion systems. This paper presents some important considerations in the redesign, in particular, the need for algorithms and software utilities that can automatically identify data flow patterns in the application program and partition and allocate calculations to the parallel processors. A library-oriented multiprocessing concept for integrating the hardware and software functions is described.

Arpasi, Dale J.↗

Minimum fuel coplanar aeroassisted orbital transfer using collocation and nonlinear programming

The fuel optimal control problem arising in coplanar orbital transfer employing aeroassisted technology is addressed. The mission involves the transfer from high energy orbit (HEO) to low energy orbit (LEO) without plane change. The basic approach here is to employ a combination of propulsive maneuvers in space and aerodynamic maneuvers in the atmosphere. The basic sequence of events for the coplanar aeroassisted HEO to LEO orbit transfer consists of three phases. In the first phase, the transfer begins with a deorbit impulse at HEO which injects the vehicle into a elliptic transfer orbit with perigee inside the atmosphere. In the second phase, the vehicle is optimally controlled by lift and drag modulation to satisfy heating constraints and to exit the atmosphere with the desired flight path angle and velocity so that the apogee of the exit orbit is the altitude of the desired LEO. Finally, the second impulse is required to circularize the orbit at LEO. The performance index is maximum final mass. Simulation results show that the coplanar aerocapture is quite different from the case where orbital plane changes are made inside the atmosphere. In the latter case, the vehicle has to penetrate deeper into the atmosphere to perform the desired orbital plane change. For the coplanar case, the vehicle needs only to penetrate the atmosphere deep enough to reduce the exit velocity so the vehicle can be captured at the desired LEO. The peak heating rates are lower and the entry corridor is wider. From the thermal protection point of view, the coplanar transfer may be desirable. Parametric studies also show the maximum peak heating rates and the entry corridor width are functions of maximum lift coefficient. The problem is solved using a direct optimization technique which uses piecewise polynomial representation for the states and controls and collocation to represent the differential equations. This converts the optimal control problem into a nonlinear programming problem which is solved numerically by using a modified version of NPSOL. Solutions were obtained for the described problem for cases with and without heating constraints. The method appears to be more robust than other optimization methods. In addition, the method can handle complex dynamical constraints.

Shi, Yun Yuan↗

Vascular mechanics of the coronary artery

This paper describes our research into the vascular mechanics of the coronary artery and plaque. The three sections describe the determination of arterial mechanical properties using intravascular ultrasound (IVUS), a constitutive relation for the arterial wall, and finite element method (FEM) models of the arterial wall and atheroma. METHODS: Inflation testing of porcine left anterior descending coronary arteries was conducted. The changes in the vessel geometry were monitored using IVUS, and intracoronary pressure was recorded using a pressure transducer. The creep and quasistatic stress/strain responses were determined. A Standard Linear Solid (SLS) was modified to reproduce the non-linear elastic behavior of the arterial wall. This Standard Non-linear Solid (SNS) was implemented into an axisymetric thick-walled cylinder numerical model. Finite element analysis models were created for five age groups and four levels of stenosis using the Pathobiological Determinants of Atherosclerosis Youth (PDAY) database. RESULTS: The arteries exhibited non-linear elastic behavior. The total tissue creep strain was epsilon creep = 0.082 +/- 0.018 mm/mm. The numerical model could reproduce both the non-linearity of the porcine data and time dependent behavior of the arterial wall found in the literature with a correlation coefficient of 0.985. Increasing age had a strong positive correlation with the shoulder stress level, (r = 0.95). The 30% stenosis had the highest shoulder stress due to the combination of a fully formed lipid pool and a thin cap. CONCLUSIONS: Studying the solid mechanics of the arterial wall and the atheroma provide important insights into the mechanisms involved in plaque rupture.

NASA Program Biomedical Research and Countermeasur↗

Computational Aeroelastic Analysis of the Semi-Span Super-Sonic Transport (S4T) Wind-Tunnel Model

A summary of the computational aeroelastic analysis for the Semi-Span Super-Sonic Transport (S4T) wind-tunnel model is presented. A broad range of analysis techniques, including linear, nonlinear and Reduced Order Models (ROMs) were employed in support of a series of aeroelastic (AE) and aeroservoelastic (ASE) wind-tunnel tests conducted in the Transonic Dynamics Tunnel (TDT) at NASA Langley Research Center. This research was performed in support of the ASE element in the Supersonics Program, part of NASA's Fundamental Aeronautics Program. The analysis concentrated on open-loop flutter predictions, which were in good agreement with experimental results. This paper is one in a series that comprise a special S4T technical session, which summarizes the S4T project.

Sanetrik, Mark D.↗

Rotor transient analysis

Undamped modes approximate dynamic behavior of rotors and bearings. Application of modal analysis to uncouple equations of motion simplifies stability, steady-state unbalance response, and transient response analysis of system; nonlinear stability is predicted from calculated frequency spectra. Analysis provides designers with complete information without involving large-scale computational costs. Programs are written in FORTRAN IV for use on CDC 6600 computer.

Allaire, P. E.↗

Parallel Computation of the Jacobian Matrix for Nonlinear Equation Solvers Using MATLAB

Demonstrating speedup for parallel code on a multicore shared memory PC can be challenging in MATLAB due to underlying parallel operations that are often opaque to the user. This can limit potential for improvement of serial code even for the so-called embarrassingly parallel applications. One such application is the computation of the Jacobian matrix inherent to most nonlinear equation solvers. Computation of this matrix represents the primary bottleneck in nonlinear solver speed such that commercial finite element (FE) and multi-body-dynamic (MBD) codes attempt to minimize computations. A timing study using MATLAB's Parallel Computing Toolbox was performed for numerical computation of the Jacobian. Several approaches for implementing parallel code were investigated while only the single program multiple data (spmd) method using composite objects provided positive results. Parallel code speedup is demonstrated but the goal of linear speedup through the addition of processors was not achieved due to PC architecture.

Rose, Geoffrey K.↗

Development of the reentry flight dynamics simulator for evaluation of space shuttle orbiter entry systems

A nonlinear, six degree of freedom, digital computer simulation of a vehicle which has constant mass properties and whose attitudes are controlled by both aerodynamic surfaces and reaction control system thrusters was developed. A rotating, oblate Earth model was used to describe the gravitational forces which affect long duration Earth entry trajectories. The program is executed in a nonreal time mode or connected to a simulation cockpit to conduct piloted and autopilot studies. The program guidance and control software used by the space shuttle orbiter for its descent from approximately 121.9 km to touchdown on the runway.

Rowell, L. F.↗

Calculated effects of varying Reynolds Number and dynamic pressure on flexible wings at transonic speeds

A computational method is described that includes the effects of static aeroelastic wing deflections in steady transonic aerodynamic calculations. This method, known as the Transonic Aero-elastic Program System (TAPS), interacts a 3D transonic computer code with boundary layer and a linear finite element structural analysis codes to calculate wing pressures and deflections. The nonlinear nature of the transonic flow makes it necessary to couple the aerodynamic and structures codes in an iterative manner. TAPS has been arranged in a modular fashion so that different aerodynamic or structures programs may be used with a minimum of coding changes required. Results obtained using two different aerodynamic codes in TAPS are given, and those results are correlated with experimental data.

Campbell, R. L.↗

Symbolic modeling of flexible manipulators

An algorithm is presented to symbolically derive the full nonlinear dynamic equations of motion of multilink flexible manipulators. Lagrange's assumed modes method is the basis of new algorithm and adapted in a way suitable for symbolic manipulation by digital computers. It is applied to model a two-link flexible arm by means of a commercially available symbolic manipulation program. The advantages of the algorithm and simulation results are discussed.

Cetinkunt, Sabri↗

An Overview of Preliminary Computational and Experimental Results for the Semi-Span Super-Sonic Transport (S4T) Wind-Tunnel Model

A summary of computational and experimental aeroelastic and aeroservoelastic (ASE) results for the Semi-Span Super-Sonic Transport (S4T) wind-tunnel model is presented. A broad range of analyses and multiple ASE wind-tunnel tests of the S4T have been performed in support of the ASE element in the Supersonics Program, part of NASA's Fundamental Aeronautics Program. The computational results to be presented include linear aeroelastic and ASE analyses, nonlinear aeroelastic analyses using an aeroelastic CFD code, and rapid aeroelastic analyses using CFD-based reduced-order models (ROMs). Experimental results from two closed-loop wind-tunnel tests performed at NASA Langley's Transonic Dynamics Tunnel (TDT) will be presented as well.

Silva, Walter A.↗

NASA/FAA general aviation crash dynamics program - An update

Work in progress in the NASA/FAA General Aviation Crash Dynamics Program for the development of technology for increased crash-worthiness and occupant survivability of general aviation aircraft is presented. Full-scale crash testing facilities and procedures are outlined, and a chronological summary of full-scale tests conducted and planned is presented. The Plastic and Large Deflection Analysis of Nonlinear Structures and Modified Seat Occupant Model for Light Aircraft computer programs which form part of the effort to predict nonlinear geometric and material behavior of sheet-stringer aircraft structures subjected to large deformations are described, and excellent agreement between simulations and experiments is noted. The development of structural concepts to attenuate the load transmitted to the passenger through the seats and subfloor structure is discussed, and an apparatus built to test emergency locator transmitters in a realistic environment is presented.

Hayduk, R. J.↗

An application of nonlinear programming to the design of regulators of a linear-quadratic formulation

A design technique is proposed for linear regulators in which a feedback controller of fixed structure is chosen to minimize an integral quadratic objective function subject to the satisfaction of integral quadratic constraint functions. Application of a nonlinear programming algorithm to this mathematically tractable formulation results in an efficient and useful computer aided design tool. Particular attention is paid to computational efficiency and various recommendations are made. Two design examples illustrate the flexibility of the approach and highlight the special insight afforded to the designer. One concerns helicopter longitudinal dynamics and the other the flight dynamics of an aerodynamically unstable aircraft.

Fleming, P.↗

Analysis of a transport fuselage section drop test

Transport fuselage section drop tests provided useful information about the crash behavior of metal aircraft in preparation for a full-scale Boeing 720 controlled impact demonstration (CID). The fuselage sections have also provided an operational test environment for the data acquisition system designed for the CID test, and data for analysis and correlation with the DYCAST nonlinear finite-element program. The correlation of the DYCAST section model predictions was quite good for the total fuselage crushing deflection (22 to 24 inches predicted versus 24 to 26 inches measured), floor deformation, and accelerations for the floor and fuselage. The DYCAST seat and occupant model was adequate to approximate dynamic loading to the floor, but a more sophisticated model would be required for good correlation with dummy accelerations. Although a full-section model using only finite elements for the subfloor was desirable, constraints of time and computer resources limited the finite-element subfloor model to a two-frame model. Results from the two-frame model indicate that DYCAST can provide excellent correlation with experimental crash behavior of fuselage structure with a minimum of empirical force-deflection data representing structure in the analytical model.

E. L. Fasanella↗

The MHOST finite element program: 3-D inelastic analysis methods for hot section components. Volume 1: Theoretical manual

Formulations and algorithms implemented in the MHOST finite element program are discussed. The code uses a novel concept of the mixed iterative solution technique for the efficient 3-D computations of turbine engine hot section components. The general framework of variational formulation and solution algorithms are discussed which were derived from the mixed three field Hu-Washizu principle. This formulation enables the use of nodal interpolation for coordinates, displacements, strains, and stresses. Algorithmic description of the mixed iterative method includes variations for the quasi static, transient dynamic and buckling analyses. The global-local analysis procedure referred to as the subelement refinement is developed in the framework of the mixed iterative solution, of which the detail is presented. The numerically integrated isoparametric elements implemented in the framework is discussed. Methods to filter certain parts of strain and project the element discontinuous quantities to the nodes are developed for a family of linear elements. Integration algorithms are described for linear and nonlinear equations included in MHOST program.

Nakazawa, Shohei↗

A study of ignition and simulation circuits for arcjet thrusters, part 1

A 1 kW electronic load was programmed to simulate the nonlinear i-v (volt-ampere) characteristics of an arcjet, both ignited and unignited. The simulator was tested and found to closely resemble an arcjet both for large transients and small perturbances up to about 40 kHz. No attempt was made to simulate the ignition process itself. The dynamic behavior of the arcjet (and the simulator) was shown to differ significantly from that of a resistor bank. Previous research led to the design and construction of a 1 kW arcjet power supply. A high voltage ignition circuit was added to this hardware, and tests on a 1 kW arcjet were performed at NASA-Lewis. All tests were successful and no ignition failures were observed. Circuit documentation and test results are included.

Stuart, Thomas A.↗

The effects of pressure sensor acoustics on airdata derived from a high-angle-of-attack flush airdata sensing (HI-FADS) system

The accuracy of a prototype nonintrusive airdata system derived for high-angle-of-attack measurements was demonstrated for quasi-steady maneuvers as great as 55 degrees during phase one of the F-18 high alpha research vehicle flight test program. This system consists of a matrix of nine pressure ports arranged in annular rings on the aircraft nose, and estimates the complete airdata set utilizing flow modeling and nonlinear regression. Particular attention is paid to the effects of acoustical distortions within the individual pressure sensors of the HI-FADS pressure matrix. A dynamic model to quantify these effects which describes acoustical distortion is developed and solved in closed form for frequency response.

Whitmore, Stephen R.↗

Orion Parachute Riser Cutter Development

This paper presents the tests and analytical approach used on the development of a steel riser cutter for the CEV Parachute Assembly System (CPAS) used on the Orion crew module. Figure 1 shows the riser cutter and the steel riser bundle which consists of six individual cables. Due to the highly compressed schedule, initial unavailability of the riser material and the Orion Forward Bay mechanical constraints, JSC primarily relied on a combination of internal ballistics analysis and LS-DYNA simulation for this project. Various one dimensional internal ballistics codes that use standard equation of state and conservation of energy have commonly used in the development of CAD devices for initial first order estimates and as an enhancement to the test program. While these codes are very accurate for propellant performance prediction, they usually lack a fully defined kinematic model for dynamic predictions. A simple piston device can easily and accurately be modeled using an equation of motion. However, the accuracy of analytical models is greatly reduced on more complicated devices with complex external loads, nonlinear trajectories or unique unlocking features. A 3D finite element model of CAD device with all critical features included can vastly improve the analytical ballistic predictions when it is used as a supplement to the ballistic code. During this project, LS-DYNA structural 3D model was used to predict the riser resisting load that was needed for the ballistic code. A Lagrangian model with eroding elements shown in Figure 2 was used for the blade, steel riser and the anvil. The riser material failure strain was fine tuned by matching the dent depth on the anvil with the actual test data. LS-DYNA model was also utilized to optimize the blade tip design for the most efficient cut. In parallel, the propellant type and the amount were determined by using CADPROG internal ballistics code. Initial test results showed a good match with LS-DYNA and CADPROG simulations. Final paper will present a detailed roadmap from initial ballistic modeling and LS-DYNA simulation to the performance testing. Blade shape optimization study will also be presented.

Oguz, Sirri↗