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,423 records · Page 79

Vector algorithms for geometrically nonlinear 3D finite element analysis

Algorithms for geometrically nonlinear finite element analysis are presented which exploit the vector processing capability of the VPS-32, which is closely related to the CYBER 205. By manipulating vectors (which are long lists of numbers) rather than individual numbers, very high processing speeds are obtained. Long vector lengths are obtained without extensive replication or reordering by storage of intermediate results in strategic patterns at all stages of the computations. Comparisons of execution times with those from programs using either scalar or other vector programming techniques indicate that the algorithms presented are quite efficient.

Whitcomb, John D.↗

Time-domain finite elements in optimal control with application to launch-vehicle guidance

A time-domain finite element method is developed for optimal control problems. The theory derived is general enough to handle a large class of problems including optimal control problems that are continuous in the states and controls, problems with discontinuities in the states and/or system equations, problems with control inequality constraints, problems with state inequality constraints, or problems involving any combination of the above. The theory is developed in such a way that no numerical quadrature is necessary regardless of the degree of nonlinearity in the equations. Also, the same shape functions may be employed for every problem because all strong boundary conditions are transformed into natural or weak boundary conditions. In addition, the resulting nonlinear algebraic equations are very sparse. Use of sparse matrix solvers allows for the rapid and accurate solution of very difficult optimization problems. The formulation is applied to launch-vehicle trajectory optimization problems, and results show that real-time optimal guidance is realizable with this method. Finally, a general problem solving environment is created for solving a large class of optimal control problems. The algorithm uses both FORTRAN and a symbolic computation program to solve problems with a minimum of user interaction. The use of symbolic computation eliminates the need for user-written subroutines which greatly reduces the setup time for solving problems.

Bless, Robert R.↗

Delamination and Stitched Failure in Stitched Composite Joints

The effect of stitches on the failure of a single lap joint configuration was determined in a combined experimental and finite element study. The experimental program was conducted to determine debond growth under static monotonic loading. The stitches were shown to delay the initiation of the debond and provide load transfer beyond the load necessary to completely debond the stitched lap joint. The experimentally determined debond length vs. applied load was used as an input parameter in the finite element analysis of both configurations. The strain energy release rates at the debond from were calculated using plate finite elements. Nonlinear fastener elements were used to model the stitches and multipoint constraints were used to model the contact problem. Models of the unstitched configuration showed significant values of modes I and II across the width of the joint and showed that mode III is zero at the centerline but increases near the free edge. Models of the stitched configuration showed that the stitches were effective in reducing mode I to zero, but had less of an effect on modes II and III.

Glaessgen, E. H.↗

Buckling Testing and Analysis of Honeycomb Sandwich Panel Arc Segments of a Full-Scale Fairing Barrel: Comparison of In- and Out-of-Autoclave Facesheet Configurations

Four honeycomb sandwich panels, representing 1/16th arc segments of a 10-m diameter barrel section of the Heavy Lift Launch Vehicle, were manufactured and tested under the NASA Composites for Exploration and the NASA Constellation Ares V programs. Two configurations were chosen for the panels: 6-ply facesheets with 1.125 in. honeycomb core and 8-ply facesheets with 1.0 in. honeycomb core. Additionally, two separate carbon fiber/epoxy material systems were chosen for the facesheets: in-autoclave IM7/977-3 and out-of-autoclave T40-800b/5320-1. Smaller 3 ft. by 5 ft. panels were cut from the 1/16th barrel sections and tested under compressive loading. Furthermore, linear eigenvalue and geometrically nonlinear finite element analyses were performed to predict the compressive response of each 3 ft. by 5 ft. panel. To improve the robustness of the geometrically nonlinear finite element model, measured surface imperfections were included in the geometry of the model. Both the linear and nonlinear models yielded good qualitative and quantitative predictions. Additionally, it was correctly predicted that the panel would fail in buckling prior to failing in strength. Furthermore, several imperfection studies were performed to investigate the influence of geometric imperfections, fiber angle misalignments, and three-dimensional effects on the compressive response of the panel.

epoxy matrix composites↗

Simulation studies using multibody dynamics code DART

DART is a multibody dynamics code developed by Photon Research Associates for the Air Force Astronautics Laboratory (AFAL). The code is intended primarily to simulate the dynamics of large space structures, particularly during the deployment phase of their missions. DART integrates nonlinear equations of motion numerically. The number of bodies in the system being simulated is arbitrary. The bodies' interconnection joints can have an arbitrary number of degrees of freedom between 0 and 6. Motions across the joints can be large. Provision for simulating on-board control systems is provided. Conservation of energy and momentum, when applicable, are used to evaluate DART's performance. After a brief description of DART, studies made to test the program prior to its delivery to AFAL are described. The first is a large angle reorientating of a flexible spacecraft consisting of a rigid central hub and four flexible booms. Reorientation was accomplished by a single-cycle sine wave shape torque input. In the second study, an appendage, mounted on a spacecraft, was slewed through a large angle. Four closed-loop control systems provided control of this appendage and of the spacecraft's attitude. The third study simulated the deployment of the rim of a bicycle wheel configuration large space structure. This system contained 18 bodies. An interesting and unexpected feature of the dynamics was a pulsing phenomena experienced by the stays whole playout was used to control the deployment. A short description of the current status of DART is given.

Keat, James E.↗

Loads Model Development and Analysis for the F/A-18 Active Aeroelastic Wing Airplane

The Active Aeroelastic Wing airplane was successfully flight-tested in March 2005. During phase 1 of the two-phase program, an onboard excitation system provided independent control surface movements that were used to develop a loads model for the wing structure and wing control surfaces. The resulting loads model, which was used to develop the control laws for phase 2, is described. The loads model was developed from flight data through the use of a multiple linear regression technique. The loads model input consisted of aircraft states and control surface positions, in addition to nonlinear inputs that were calculated from flight-measured parameters. The loads model output for each wing consisted of wing-root bending moment and torque, wing-fold bending moment and torque, inboard and outboard leading-edge flap hinge moment, trailing-edge flap hinge moment, and aileron hinge moment. The development of the Active Aeroelastic Wing loads model is described, and the ability of the model to predict loads during phase 2 research maneuvers is demonstrated. Results show a good match to phase 2 flight data for all loads except inboard and outboard leading-edge flap hinge moments at certain flight conditions. The average load prediction errors for all loads at all flight conditions are 9.1 percent for maximum stick-deflection rolls, 4.4 percent for 5-g windup turns, and 7.7 percent for 4-g rolling pullouts.

Allen, Michael J.↗

Loads Model Development and Analysis for the F/A-18 Active Aeroelastic Wing Airplane

The Active Aeroelastic Wing airplane was successfully flight-tested in March 2005. During phase 1 of the two-phase program, an onboard excitation system provided independent control surface movements that were used to develop a loads model for the wing structure and wing control surfaces. The resulting loads model, which was used to develop the control laws for phase 2, is described. The loads model was developed from flight data through the use of a multiple linear regression technique. The loads model input consisted of aircraft states and control surface positions, in addition to nonlinear inputs that were calculated from flight-measured parameters. The loads model output for each wing consisted of wing-root bending moment and torque, wing-fold bending moment and torque, inboard and outboard leading-edge flap hinge moment, trailing-edge flap hinge moment, and aileron hinge moment. The development of the Active Aeroelastic Wing loads model is described, and the ability of the model to predict loads during phase 2 research maneuvers is demonstrated. Results show a good match to phase 2 flight data for all loads except inboard and outboard leading-edge flap hinge moments at certain flight conditions. The average load prediction errors for all loads at all flight conditions are 9.1 percent for maximum stick-deflection rolls, 4.4 percent for 5-g windup turns, and 7.7 percent for 4-g rolling pullouts.

Allen, Michael J.↗

Polarized Bidirectional Reflectance of Optically Thick Sparse Particulate Layers: an Efficient Numerically Exact Radiative-Transfer Solution

We describe a simple yet efficient numerical algorithm for computing polarized bidirectional reflectance of an optically thick (semi-infinite), macroscopically flat layer composed of statistically isotropic and mirror symmetric random particles. The spatial distribution of the particles is assumed to be sparse, random, and statistically uniform. The 44 Stokes reflection matrix is calculated by iterating the Ambartsumian's vector nonlinear integral equation. The result is a numerically exact solution of the vector radiative transfer equation and as such fully satisfies the energy conservation law and the fundamental reciprocity relation. Since this technique bypasses the computation of the internal radiation field, it is very fast and highly accurate. The FORTRAN implementation of the technique is publicly available on the World Wide Web at http://www.giss.nasa.gov/staff/ mmishchenko/brf. It can be combined with several existing computer programs providing the requisite single-scattering properties of spherical or morphologically complex particles and applied to a wide range of optical characterization problems. Benchmark results obtained with this program can be used for testing alternative solvers of the vector radiative transfer equation.

radiative transfer↗

Flight Performance and Stability of Space Launch System Core Stage Thrust Vector Control

The Space Launch System (SLS) Core Stage (CS) Thrust Vector Control (TVC) system is comprised of eight mechanical feedback Shuttle heritage Type III TVC actuators and four RS-25 engines, each attached to a Shuttle heritage gimbal block/bearing. The Core Stage TVC shares vehicle control authority with the SLS 5-segment Solid Rocket Boosters (SRBs) during boost phase flight, and is the sole means of vehicle flight control during in exoatmospheric flight following SRB separation. TVC responses during Green Run Hot Fire (GRHF) testing revealed that the TVC did not meet its performance specifications. Step and frequency responses exhibited unexpected departures from prior laboratory data and modeled behavior. Post-test analysis determined that the characteristics of the structure and gimbal friction are significantly influenced by the thrust-loaded conditions, and the command avionics exhibited a small but important gain nonlinearity. Using the available test data, the design team augmented the flight control TVC models to bound the observed results and include the additional fidelity needed for vehicle flight control analysis so as to build sufficient rationale for flight certification. Prior to the Green Run tests, “simplex” linear models typically used for flight control analysis did not include gimbal friction and other nonlinearities owing to long-standing assumptions that these effects were negligible in the Shuttle Orbiter TVC system. Following the Green Run findings, simulation analysis of the flight dynamics in the time and frequency domain revealed the propensity for a flight control limit cycle oscillation (LCO) if friction and structural compliances fell near the edges of test-predicted bounds. While the “most probable” models did not predict an in-flight LCO, the SLS Program conservatively proceeded with a system-wide evaluation and ultimate acceptance of the possibility for a small amplitude, low-frequency TVC LCO in flight. A final validation of the extensive test and modeling effort occurred when the first flight of SLS successfully demonstrated the fully integrated performance of the vehicle’s TVC system This paper is the final installment in a seven-paper series surveying the design, engineering, test validation, and flight performance of the Core Stage Thrust Vector Control system. In this paper, the development of flight rationale in light of the TVC responses observed in Green Run is discussed, along with a review of the flight telemetry illustrating the correlation of the preflight predictions with the observed performance.

John H. Wall↗

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗

SCI Identification (SCIDNT) program user's guide

The computer program Linear SCIDNT which evaluates rotorcraft stability and control coefficients from flight or wind tunnel test data is described. It implements the maximum likelihood method to maximize the likelihood function of the parameters based on measured input/output time histories. Linear SCIDNT may be applied to systems modeled by linear constant-coefficient differential equations. This restriction in scope allows the application of several analytical results which simplify the computation and improve its efficiency over the general nonlinear case.

Source record↗

Derivatives of buckling loads and vibration frequencies with respect to stiffness and initial strain parameters

A uniform variational approach to sensitivity analysis of vibration frequencies and bifurcation loads of nonlinear structures is developed. Two methods of calculating the sensitivities of bifurcation buckling loads and vibration frequencies of nonlinear structures, with respect to stiffness and initial strain parameters, are presented. A direct method requires calculation of derivatives of the prebuckling state with respect to these parameters. An adjoint method bypasses the need for these derivatives by using instead the strain field associated with the second-order postbuckling state. An operator notation is used and the derivation is based on the principle of virtual work. The derivative computations are easily implemented in structural analysis programs. This is demonstrated by examples using a general purpose, finite element program and a shell-of-revolution program.

Haftka, Raphael T.↗

Resource Letter GrW-1: Gravitational Waves

The phenomenon of gravitational radiation was one of the first predictions of Einstein's general theory of relativity. Progress in understanding this radiation theoretically was slow at first, owing to the difficulty of the nonlinear field equations and the subtleties of their physical effects. The experimental side of this subject also has taken a long time to develop, with efforts at detection severely challenged by the extreme weakness of the waves impinging on the Earth. However, as the 21st century begins, observations of the gravitational waves from astrophysical sources such as black holes, neutron stars, and stellar collapse are expected to open a new window on the universe. Vigorous experimental programs centered on ground-based detectors are being carried out worldwide, and a space-based detector is in the planning stages. On the theoretical side, much effort is being expended to produce robust models of the astrophysical sources and accurate calculations of the waveforms they produce. In this Resource Letter, a set of basic references will be presented first, to provide a general introduction to and overview of the literature in this field. The focus then will shift to highlighting key resources in more specialized areas at the forefront of current research.

White, Nicholas E.↗

Pressure Testing of a Minimum Gauge PRSEUS Panel

Advanced aircraft configurations that have been developed to increase fuel efficiency require advanced, novel structural concepts capable of handling the unique load conditions that arise. One such concept is the Pultruded Rod Stitched Efficient Unitized Structure (PRSEUS) developed by the Boeing Company. The PRSEUS concept is being investigated by NASA s Environmentally Responsible Aviation (ERA) Program for use in a hybrid-wing body (HWB) aircraft. This paper summarizes the analysis and test of a PRSEUS panel subjected to internal pressure, the first such pressure test for this structural concept. The pressure panel used minimum gauge skin, with stringer and frame configurations consistent with previous PRSEUS tests. Analysis indicated that for the minimum gauge skin panel, the stringer locations exhibit fairly linear response, but the skin bays between the stringers exhibit nonlinear response. Excellent agreement was seen between nonlinear analysis and test results in the critical portion at the center of the panel. The pristine panel was capable of withstanding the required 18.4 psi pressure load condition without exhibiting any damage. The impacted panel was capable of withstanding a pressure load in excess of 28 psi before initial failure occurred at the center stringer, and the panel was capable of sustaining increased pressure load after the initial failure. This successful PRSEUS panel pressure panel test was a critical step in the building block approach for enabling the use of this advanced structural concept on future aircraft, such as the HWB.

Lovejoy, Andrew J.↗

Solving the Dynamics-Aware Economic Dispatch Problem with the Koopman Operator

The dynamics-aware economic dispatch (DED) problem embeds low-level generator dynamics and operational constraints to enable near real-time scheduling of generation units in a power network. DED produces a more dynamic supervisory control policy than traditional economic dispatch (T-ED) that reduces overall generation costs. However, in contrast to T-ED, DED is a nonlinear, non-convex optimization problem that is computationally prohibitive to solve. We introduce a machine learning-based operator-theoretic approach for solving the DED problem efficiently. Specifically, we develop a novel discrete-time Koopman Operator (KO) formulation that embeds domain information into the structure of the KO to learn high-fidelity approximations of the generator dynamics. Using the KO approximation, the DED problem can be reformulated as a computationally tractable linear program (abbreviated DED-KO). We demonstrate the high solution quality and computational-time savings of the DED-KO model over the original DED formulation on a 9-bus test system.

King, Ethan↗

Optimal feedback control for nonlinear control moment gyro systems

A study was made, leading to development and analysis of a technique allowing the control moment gyro (CMG) momentum vectors to avoid adverse orientations and thus to provide optimal control for the Apollo Telescope Mount. The technique was formulated in terms of a feedback scheme using a check factor, K, formed as the product of the dot products of the three CMG momentum unit vectors. The parameter K approaches unity for all adverse orientations and was used to initiate a single correction maneuver to rotate the CMG momenta away from a pending adverse orientation and to reduce K to an acceptable value. A digital computer program for ATM simulations was modified to incorporate this correction maneuver and performance was evaluated in simulated cases. Antiparallelism was used as the primary example for the evaluation of the feedback optimization technique.

Kyle, R. G.↗

Computer studies of baroclinic flow

Programs necessary for computing the transition curve on the regime diagram for the atmospheric general circulation experiment (AGOE) were completed and used to determine the regime diagram for the rotating annulus and some axisymmetric flows for one possible AGOE configuration. The effect of geometrical constraints on the size of eddies developing from a basic state is being examined. In AGOE, the geometric constraint should be the width of the shear zone or the baroclinic zone. Linear and nonlinear models are to be used to examine both barotropic and baroclinic flows. The results should help explain the scale selection mechanism of baroclinic eddies in the atmosphere experimental models such as AGOE, and the multiple vortex phenomenon in tornadoes.

Gall, R.↗

Application of traditional CFD methods to nonlinear computational aeroacoustics problems

This paper describes an implementation of a high order finite difference technique and its application to the category 2 problems of the ICASE/LaRC Workshop on Computational Aeroacoustics (CAA). Essentially, a popular Computational Fluid Dynamics (CFD) approach (central differencing, Runge-Kutta time integration and artificial dissipation) is modified to handle aeroacoustic problems. The changes include increasing the order of the spatial differencing to sixth order and modifying the artificial dissipation so that it does not significantly contaminate the wave solution. All of the results were obtained from the CM5 located at the Numerical Aerodynamic Simulation Laboratory. lt was coded in CMFortran (very similar to HPF), using programming techniques developed for communication intensive large stencils, and ran very efficiently.

Chyczewski, Thomas S.↗