SEARCH · Engineering Papers
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.
The development of a general mathematical model and solutions of test problems to analyze large nonisothermal elasto-visco-plastic deformatisms of structures is discussed. Geometric and material type nonlinearities of higher order are present in the development of the mathematical model and in the developed solution methodology.
Analysis and synthesis tools are developed to improved the dynamic performance of reconfigurable nonminimum phase, nonstrictly positive real-time variant systems. A novel Spline Varying Optimal (SVO) controller is developed for the kinematic nonlinear system. There are several advantages to using the SVO controller, in which the spline function approximates the system model, observer, and controller gain. They are: The spline function approximation is simply connected, thus the SVO controller is more continuous than traditional gain scheduled controllers when implemented on a time varying plant; ft is easier for real-time implementations in storage and computational effort; where system identification is required, the spline function requires fewer experiments, namely four experiments; and initial startup estimator transients are eliminated. The SVO compensator was evaluated on a high fidelity simulation of the Shuttle Remote Manipulator System. The SVO controller demonstrated significant improvement over the present arm performance: (1) Damping level was improved by a factor of 3; and (2) Peak joint torque was reduced by a factor of 2 following Shuttle thruster firings.
The effect of stitching on mode I and mode II strain energy release rates for debond configurations is studied using an analysis based on plate finite elements and the virtual crack closure technique. The stitches were modeled as discrete nonlinear fastener elements with a compliance determined by experiment. The axial and shear behavior of the stitches was considered with both the compliances and failure loads assumed to be independent. The mode I strain energy release rate, G(sub I), was shown to decrease once the debond had grown beyond the first row of stitches and was reduced to zero for long debonds, however, the mode II strain energy release rate, G(sub II), continued to be of significant magnitude over the range of debond lengths considered.
Three procedures to guide selection of an efficient modal basis in a nonlinear random response analysis are examined. One method is based only on proper orthogonal decomposition, while the other two additionally involve smooth orthogonal decomposition. Acoustic random response problems are employed to assess the performance of the three modal basis selection approaches. A thermally post-buckled beam exhibiting snap-through behavior, a shallowly curved arch in the auto-parametric response regime and a plate structure are used as numerical test articles. The results of a computationally taxing full-order analysis in physical degrees of freedom are taken as the benchmark for comparison with the results from the three reduced-order analyses. For the cases considered, all three methods are shown to produce modal bases resulting in accurate and computationally efficient reduced-order nonlinear simulations.
A parametric study of the effects of test-fixture-induced initial prestress and elastic edge restraints on the prebuckling and buckling responses of a compression-loaded, quasi-isotropic curved panel is presented. The numerical results were obtained by using a geometrically nonlinear finite element analysis code with high-fidelity models. The results presented show that a wide range of prebuckling and buckling behavior can be obtained by varying parameters that represent circumferential loaded-edge restraint and rotational unloaded-edge restraint provided by a test fixture and that represent the mismatch in specimen and test-fixture radii of curvature. For a certain range of parameters, the panels exhibit substantial nonlinear prebuckling deformations that yield buckling loads nearly twice the corresponding buckling load predicted by a traditional linear bifurcation buckling analysis for shallow curved panels. In contrast, the results show another range of parameters exist for which the nonlinear prebuckling deformations either do not exist or are relatively benign, and the panels exhibit buckling loads that are nearly equal to the corresponding linear bifurcation buckling load. Overall, the results should also be of particular interest to scientists, engineers, and designers involved in simulating flight-hardware boundary conditions in structural verification and certification tests, involved in validating structural analysis tools, and interested in tailoring buckling performance.
The goal of this research is to develop the transfer matrix method to treat nonlinear autonomous boundary value problems with multiple branches. The application is the complete nonlinear aeroelastic analysis of multiple-branched rotor blades. Once the development is complete, it can be incorporated into the existing transfer matrix analyses. There are several difficulties to be overcome in reaching this objective. The conventional transfer matrix method is limited in that it is applicable only to linear branch chain-like structures, but consideration of multiple branch modeling is important for bearingless rotors. Also, hingeless and bearingless rotor blade dynamic characteristics (particularly their aeroelasticity problems) are inherently nonlinear. The nonlinear equations of motion and the multiple-branched boundary value problem are treated together using a direct transfer matrix method. First, the formulation is applied to a nonlinear single-branch blade to validate the nonlinear portion of the formulation. The nonlinear system of equations is iteratively solved using a form of Newton-Raphson iteration scheme developed for differential equations of continuous systems. The formulation is then applied to determine the nonlinear steady state trim and aeroelastic stability of a rotor blade in hover with two branches at the root. A comprehensive computer program is developed and is used to obtain numerical results for the (1) free vibration, (2) nonlinearly deformed steady state, (3) free vibration about the nonlinearly deformed steady state, and (4) aeroelastic stability tasks. The numerical results obtained by the present method agree with results from other methods.
An objective of the HPCC Program at NASA Langley has been to promote the use of advanced computing techniques to more rapidly solve the problem of multidisciplinary optimization of a supersonic transport configuration. As a result, a software system has been designed and is being implemented to integrate a set of existing discipline analysis codes, some of them CPU-intensive, into a distributed computational framework for the design of a High Speed Civil Transport (HSCT) configuration. The proposed paper will describe the engineering aspects of integrating these analysis codes and additional interface codes into an automated design system. The objective of the design problem is to optimize the aircraft weight for given mission conditions, range, and payload requirements, subject to aerodynamic, structural, and performance constraints. The design variables include both thicknesses of structural elements and geometric parameters that define the external aircraft shape. An optimization model has been adopted that uses the multidisciplinary analysis results and the derivatives of the solution with respect to the design variables to formulate a linearized model that provides input to the CONMIN optimization code, which outputs new values for the design variables. The analysis process begins by deriving the updated geometries and grids from the baseline geometries and grids using the new values for the design variables. This free-form deformation approach provides internal FEM (finite element method) grids that are consistent with aerodynamic surface grids. The next step involves using the derived FEM and section properties in a weights process to calculate detailed weights and the center of gravity location for specified flight conditions. The weights process computes the as-built weight, weight distribution, and weight sensitivities for given aircraft configurations at various mass cases. Currently, two mass cases are considered: cruise and gross take-off weight (GTOW). Weights information is obtained from correlations of data from three sources: 1) as-built initial structural and non-structural weights from an existing database, 2) theoretical FEM structural weights and sensitivities from Genesis, and 3) empirical as-built weight increments, non-structural weights, and weight sensitivities from FLOPS. For the aeroelastic analysis, a variable-fidelity aerodynamic analysis has been adopted. This approach uses infrequent CPU-intensive non-linear CFD to calculate a non-linear correction relative to a linear aero calculation for the same aerodynamic surface at an angle of attack that results in the same configuration lift. For efficiency, this nonlinear correction is applied after each subsequent linear aero solution during the iterations between the aerodynamic and structural analyses. Convergence is achieved when the vehicle shape being used for the aerodynamic calculations is consistent with the structural deformations caused by the aerodynamic loads. To make the structural analyses more efficient, a linearized structural deformation model has been adopted, in which a single stiffness matrix can be used to solve for the deformations under all the load conditions. Using the converged aerodynamic loads, a final set of structural analyses are performed to determine the stress distributions and the buckling conditions for constraint calculation. Performance constraints are obtained by running FLOPS using drag polars that are computed using results from non-linear corrections to the linear aero code plus several codes to provide drag increments due to skin friction, wave drag, and other miscellaneous drag contributions. The status of the integration effort will be presented in the proposed paper, and results will be provided that illustrate the degree of accuracy in the linearizations that have been employed.
The modified structure method, which is based on Koiter's theory of imperfections, was used to calculate approximate buckling loads of several shells of revolution. The method does not appear to be practical for shells because, in many cases, the prebuckling nonlinearity may be too large to be treated accurately as a small imperfection.
This report documents a validation study on the standard Open Hole Compression (OHC) laminate test specimen conducted as part of the NASA Advanced Composites Project (ACP). Tests were conducted on OHC specimens with hard, quasi-isotropic, and soft layups using digital image correlation and X-Ray computed tomography to capture the structural response and damage evolution. Progressive damage models were constructed for use with the CompDam continuum damage mechanics code following the best practices established during the ACP. Detailed interrogation of the analysis results and comparison with experimental measurements provide a basis for assessing the capability of the modeling approach for OHC. The structural response is found to be captured well, with strengths predicted within 3\% of the experimental values for hard and quasi-isotropic laminates. In the soft laminate, the model predicts failure to be more brittle than the nonlinear, ductile response that was measured. Damage states extracted from the models at the same load level as test measurements are overlaid to show directly the similarities and differences between test and analysis results. Studying the damage evolution predicted by the analysis reveals that the failure process is a competition between fiber damage and delamination/sub-laminate buckling, with fiber damage dominating the collapse in the hard laminate and sub-laminate buckling governing in the soft laminate. Finally, a series of parametric studies varying numerical solution parameters (mesh size, mass scaling) and physical properties (fiber direction compressive strength and toughness) reveal sensitivities and deficiencies of the model. To the authors' knowledge, this study is the first for OHC specimens to include detailed evaluation of damage mode interactions, direct overlay of predicted and measured damage states, and sensitivity of the predicted results to difficult-to-measure fiber direction material properties.
A method of analysis for thin structures that incorporates finite strain, elastic-plastic, strain hardening, time dependent material behavior implemented with respect to a fixed configuration and is consistently valid for finite strains and finite rotations is developed. The theory is formulated systematically in a body fixed system of convected coordinates with materially embedded vectors that deform in common with continuum. Tensors are considered as linear vector functions and use is made of the dyadic representation. The kinematics of a deformable continuum is treated in detail, carefully defining precisely all quantities necessary for the analysis. The finite strain theory developed gives much better predictions and agreement with experiment than does the traditional small strain theory, and at practically no additional cost. This represents a very significant advance in the capability for the reliable prediction of nonlinear transient structural responses, including the reliable prediction of strains large enough to produce ductile metal rupture.
Determining the magnetic field structure, electric currents, and plasma distributions within flux transfer event (FTE)-type flux ropes is critical to the understanding of their origin, evolution, and dynamics. Here the Magnetospheric Multiscale mission's high-resolution magnetic field and plasma measurements are used to identify FTEs in the vicinity of the subsolar magnetopause. The constant-alpha flux rope model is used to identify quasi-force free flux ropes and to infer the size, the core magnetic field strength, the magnetic flux content, and the spacecraft trajectories through these structures. Our statistical analysis determines a mean diameter of 1,700 ± 400 km (~30 ± 9 d(sub i)) and an average magnetic flux content of 100 ± 30 kWb for the quasi-force free FTEs at the Earth's subsolar magnetopause which are smaller than values reported by Cluster at high latitudes. These observed nonlinear size and magnetic flux content distributions of FTEs appear consistent with the plasmoid instability theory, which relies on the merging of neighboring, small-scale FTEs to generate larger structures. The ratio of the perpendicular to parallel components of current density, R(sub J), indicates that our FTEs are magnetically force-free, defined as R(sub J) < 1, in their core regions (<0.6 R(sub flux rope)). Plasma density is shown to be larger in smaller, newly formed FTEs and dropping with increasing FTE size. It is also shown that parallel ion velocity dominates inside FTEs with largest plasma density. Field-aligned flow facilitates the evacuation of plasma inside newly formed FTEs, while their core magnetic field strengthens with increasing FTE size.
A nonlinear aerodynamic prediction technique based on the full potential equation in conservation form has been developed for the treatment of supersonic flows. The method uses the theory of characteristic signal propagation to accurately simulate the flow structure, which includes shock waves and mixed elliptic-hyperbolic crossflow. An implicit approximate factorization scheme is employed to solve the finite-differenced equation. The necessary body-fitted grid system in every marching plane is generated numerically, using an elliptic grid solver. Results are shown for conical and nonconical wing-body combinations and compared with experimental data and Euler calculations. The method demonstrates an enormous savings in execution time and memory requirements over Euler methods.
A method for the estimation of parameters and error analysis in the development of nonlinear modeling for environmental impact assessment studies is presented. The modular computer program can interactively fit different nonlinear models to the same set of data, dynamically changing the error structure associated with observed values. Parameter estimation techniques and sequential estimation algorithms employed in parameter identification and model selection are first discussed. Then, least-square parameter estimation procedures are formulated, utilizing differential or integrated equations, and are used to define a model for association of error with experimentally observed data.
Selected research and technology activities at NASA Dryden Flight Research Center are summarized. These following activities exemplify the Center's varied and productive research efforts: Developing a Requirements Development Guide for an Automatic Ground Collision Avoidance System; Digital Terrain Data Compression and Rendering for Automatic Ground Collision Avoidance Systems; Nonlinear Flutter/Limit Cycle Oscillations Prediction Tool; Nonlinear System Identification Using Orthonormal Bases: Application to Aeroelastic/Aeroservoelastic Systems; Critical Aerodynamic Flow Feature Indicators: Towards Application with the Aerostructures Test Wing; Multidisciplinary Design, Analysis, and Optimization Tool Development Using a Genetic Algorithm; Structural Model Tuning Capability in an Object-Oriented Multidisciplinary Design, Analysis, and Optimization Tool; Extension of Ko Straight-Beam Displacement Theory to the Deformed Shape Predictions of Curved Structures; F-15B with Phoenix Missile and Pylon Assembly--Drag Force Estimation; Mass Property Testing of Phoenix Missile Hypersonic Testbed Hardware; ARMD Hypersonics Project Materials and Structures: Testing of Scramjet Thermal Protection System Concepts; High-Temperature Modal Survey of the Ruddervator Subcomponent Test Article; ARMD Hypersonics Project Materials and Structures: C/SiC Ruddervator Subcomponent Test and Analysis Task; Ground Vibration Testing and Model Correlation of the Phoenix Missile Hypersonic Testbed; Phoenix Missile Hypersonic Testbed: Performance Design and Analysis; Crew Exploration Vehicle Launch Abort System-Pad Abort-1 (PA-1) Flight Test; Testing the Orion (Crew Exploration Vehicle) Launch Abort System-Ascent Abort-1 (AA-1) Flight Test; SOFIA Flight-Test Flutter Prediction Methodology; SOFIA Closed-Door Aerodynamic Analyses; SOFIA Handling Qualities Evaluation for Closed-Door Operations; C-17 Support of IRAC Engine Model Development; Current Capabilities and Future Upgrade Plans of the C-17 Data Rack; Intelligent Data Mining Capabilities as Applied to Integrated Vehicle Health Management; STARS Flight Demonstration No. 2 IP Data Formatter; Space-Based Telemetry and Range Safety (STARS) Flight Demonstration No. 2 Range User Flight Test Results; Aerodynamic Effects of the Quiet Spike(tm) on an F-15B Aircraft; F-15 Intelligent Flight Controls-Increased Destabilization Failure; F-15 Integrated Resilient Aircraft Control (IRAC) Improved Adaptive Controller; Aeroelastic Analysis of the Ikhana/Fire Pod System; Ikhana: Western States Fire Missions Utilizing the Ames Research Center Fire Sensor; Ikhana: Fiber-Optic Wing Shape Sensors; Ikhana: ARTS III; SOFIA Closed-Door Flutter Envelope Flight Testing; F-15B Quiet Spike(TM) Aeroservoelastic Flight Test Data Analysis; and UAVSAR Platform Precision Autopilot Flight Results.
Nonlinear programming algorithms play an important role in structural design optimization. Several such algorithms have been implemented in OpenMDAO framework developed at NASA Glenn Research Center (GRC). OpenMDAO is an open source engineering analysis framework, written in Python, for analyzing and solving Multi-Disciplinary Analysis and Optimization (MDAO) problems. It provides a number of solvers and optimizers, referred to as components and drivers, which users can leverage to build new tools and processes quickly and efficiently. Users may download, use, modify, and distribute the OpenMDAO software at no cost. This paper summarizes the process involved in analyzing and optimizing structural components by utilizing the framework s structural solvers and several gradient based optimizers along with a multi-objective genetic algorithm. For comparison purposes, the same structural components were analyzed and optimized using CometBoards, a NASA GRC developed code. The reliability and efficiency of the OpenMDAO framework was compared and reported in this report.
Nonlinear programming algorithms play an important role in structural design optimization. Several such algorithms have been implemented in OpenMDAO framework developed at NASA Glenn Research Center (GRC). OpenMDAO is an open source engineering analysis framework, written in Python, for analyzing and solving Multi-Disciplinary Analysis and Optimization (MDAO) problems. It provides a number of solvers and optimizers, referred to as components and drivers, which users can leverage to build new tools and processes quickly and efficiently. Users may download, use, modify, and distribute the OpenMDAO software at no cost. This paper summarizes the process involved in analyzing and optimizing structural components by utilizing the framework s structural solvers and several gradient based optimizers along with a multi-objective genetic algorithm. For comparison purposes, the same structural components were analyzed and optimized using CometBoards, a NASA GRC developed code. The reliability and efficiency of the OpenMDAO framework was compared and reported in this report.
This paper is concerned with the high Reynolds number flow over a spanwise periodic array of roughness elements with inter-element spacing of the order of the local boundary-layer thickness. While earlier work by Goldstein, Sescu, Duck and Choudhari (2010) and Goldstein, Sescu, Duck and Choudhari (2011) was mainly concerned with smaller roughness heights that produced relatively weak distortions of the downstream flow, the focus here is on extending the analysis to larger roughness heights and streamwise elongated planform shapes that together produce a qualitatively different, nonlinear behavior of the downstream wakes. The roughness scale flow now has a novel triple-deck structure that is somewhat different from related studies that have previously appeared in the literature. The resulting flow is formally nonlinear in the intermediate wake region, where the streamwise distance is large compared to the roughness dimensions but small compared to the downstream distance from the leading edge, as well as in the far wake region where the streamwise length scale is of the order of the downstream distance from the leading edge. In contrast, the flow perturbations in both of these wake regions were strictly linear in the earlier work by Goldstein et al (2010, 2011). This is an important difference because the nonlinear wake flow in the present case provides an appropriate basic state for studying the secondary instability and eventual breakdown into turbulence.
An analytical method for computing unsteady flow conditions in liquid-filled flow systems is developed. The method which is called the wave plan incorporates distributed parameter and nonlinear effects including the effects of viscous resistance. The wave plan is essentially a solution synthesized from the effects of incremental step pressure pulses. The pressure pulses are generated because of incremental flow-rate changes that originate in a hydraulic system from a variety of sources, including the mechanical motion of the system structure. The pressure pulses propagate throughout the system at sonic velocity and are partially transmitted and reflected at each discontinuity. The velocity change caused by each pressure pulse is obtained from the Joukowski relation. Pressure and velocity time histories at any point in the system are obtained by a timewise summation of the contributions of the incremental pressure pulses passing that point. The analysis is presented in a form general enough to be applied to a variety of hydraulic systems. To illustrate the application of the method to a specific system, the response of a straight hydraulic line to a sinusoidal orifice-area variation of an upstream valve is computed. Both a constant-cross- section line and a tapered line are analyzed in the examples, and various nonlinear effects evaluated. Comparisons are carried out with experimental data obtained for the constant-diameter line and good agreement is shown to exist.