Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural modeling”

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 829 records · Page 46

System/observer/controller identification toolbox

System Identification is the process of constructing a mathematical model from input and output data for a system under testing, and characterizing the system uncertainties and measurement noises. The mathematical model structure can take various forms depending upon the intended use. The SYSTEM/OBSERVER/CONTROLLER IDENTIFICATION TOOLBOX (SOCIT) is a collection of functions, written in MATLAB language and expressed in M-files, that implements a variety of modern system identification techniques. For an open loop system, the central features of the SOCIT are functions for identification of a system model and its corresponding forward and backward observers directly from input and output data. The system and observers are represented by a discrete model. The identified model and observers may be used for controller design of linear systems as well as identification of modal parameters such as dampings, frequencies, and mode shapes. For a closed-loop system, an observer and its corresponding controller gain directly from input and output data.

Juang, Jer-Nan↗

Real-Time Adaptive Drag Minimization Wind Tunnel Investigation of a Flexible Wing with Variable Camber Continuous Trailing Edge Flap System

This paper reports the results of a recently completed real-time adaptive drag minimization wind tunnel investigation of a highly flexible wing wind tunnel model equipped with the Variable Camber Continuous Trailing Flap (VCCTEF) technology at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel investigation is funded by NASA SBIR Phase II contract with Scientific Systems Company, Inc. (SSCI) and University of Washington (UW) as a subcontractor. The wind tunnel model is a sub-scale Common Research Model (CRM) wing constructed of foam core and fiberglass skin and is aeroelastically scaled to achieve a wing tip deflection of 10% of the wing semi-span which represents a typical wing tip deflection for a modern transport such as Boeing 787. The jig-shape twist of the CRM wing is optimized using a CART3D aero-structural model to achieve the minimum induced drag for the design cruise lift coefficient of 0.5. The wing is equipped with two chord wise cambered segments for each of the six span wise flap sections for a total of 12 individual flap segments that comprise the VCCTEF system. Each of the 12 flap segments is actively controlled by an electric servo-actuator. The real-time adaptive drag optimization strategy includes an on-board aerodynamic model identification, a model excitation, and a real-time drag optimization. The on-board aerodynamic model is constructed parametrically as a function of the angle of attack and flap positions to model the lift and drag coefficients of the wing. The lift coefficient models include a linear model and a second-order model. The drag coefficient models include a quadratic model and a higher-order up to 6th-order model to accurately model the drag coefficient at high angles of attack. The onboard aerodynamic model identification includes a recursive least-squares (RLS) algorithm and a batch least-squares (BLS) algorithm designed to estimate the model parameters. The model excitation method is designed to sample the input set that comprises the angle of attack and the flap positions. Three model excitation methods are developed: random excitation method, sweep method, and iterative angle-of-attack seeking method. The real-time drag optimization includes a generic algorithm developed by SSCI and several optimization methods developed by NASA which include a second-order gradient Newton-Raphson optimization method, an iterative gradient optimization method, a pseudo-inverse optimization method, an analytical optimization method, and an iterative refinement optimization method. The first wind tunnel test entry took place in September 2017. This test revealed major hardware issues and required further redesign of the flap servo mechanisms. The second test entry took place in April 2018. However, the test was not successful due to the issues with the onboard aerodynamic model identification RLS algorithm which incorrectly identified model parameters. This test also provides an experimental comparison study between the VCCTEF and a variable camber discrete trailing edge flap (VCDTEF) without the elastomer transition mechanisms. The experimental result confirms the benefit of the VCCTEF which produces lower drag by 5% than the VCDTEF. The third and final test entry took place in June 2018 after the issues with the RLS algorithm have been identified and corrected. Additional improvements were implemented. These include the BLS algorithm, the iterative angle-of-attack seeking method, the iterative gradient optimization method, and the pseudo-inverse optimization method. The test objectives were successfully demonstrated as the real-time drag optimization identifies several optimal solutions at off-design lift coefficients. The iterative gradient optimization method is found to achieve up to 4.7% drag reduction for the off-design lift coefficient of 0.7. The pseudo-inverse optimization method which does not require the drag coefficient model is found to be quite effective in reducing drag. Up to 9.4% drag reduction for the off-design lift coefficient of 0.7 is achieved with the pseudo-inverse optimization method. The wind tunnel investigation demonstrates the potential of real-time drag optimization technology. Several new capabilities are developed that could enable future adaptive wing technologies for flexible wings equipped with drag control devices such as the VCCTEF.

Drag Optimization↗

Aerodynamic Parameters of High Performance Aircraft Estimated from Wind Tunnel and Flight Test Data

A concept of system identification applied to high performance aircraft is introduced followed by a discussion on the identification methodology. Special emphasis is given to model postulation using time invariant and time dependent aerodynamic parameters, model structure determination and parameter estimation using ordinary least squares an mixed estimation methods, At the same time problems of data collinearity detection and its assessment are discussed. These parts of methodology are demonstrated in examples using flight data of the X-29A and X-31A aircraft. In the third example wind tunnel oscillatory data of the F-16XL model are used. A strong dependence of these data on frequency led to the development of models with unsteady aerodynamic terms in the form of indicial functions. The paper is completed by concluding remarks.

Klein, Vladislav↗

Aerodynamic Parameters of High Performance Aircraft Estimated from Wind Tunnel and Flight Test Data

A concept of system identification applied to high performance aircraft is introduced followed by a discussion on the identification methodology. Special emphasis is given to model postulation using time invariant and time dependent aerodynamic parameters, model structure determination and parameter estimation using ordinary least squares and mixed estimation methods. At the same time problems of data collinearity detection and its assessment are discussed. These parts of methodology are demonstrated in examples using flight data of the X-29A and X-31A aircraft. In the third example wind tunnel oscillatory data of the F-16XL model are used. A strong dependence of these data on frequency led to the development of models with unsteady aerodynamic terms in the form of indicial functions. The paper is completed by concluding remarks.

Klein, Vladislav↗

Hypersonic Inflatable Aerodynamic Decelerator Ground Test Development

Hypersonic Inflatable Aerodynamic Decelerator (HIAD) technology readiness levels have been incrementally matured by NASA over the last thirteen years, with most recent support from NASA's Space Technology Mission Directorate (STMD) Game Changing Development Program (GCDP). Recently STMD GCDP has authorized funding and support through fiscal year 2015 (FY15) for continued HIAD ground developments which support a Mars Entry, Descent, and Landing (EDL) study. The Mars study will assess the viability of various EDL architectures to enable a Mars human architecture pathfinder mission planned for mid-2020. At its conclusion in November 2014, NASA's first HIAD ground development effort had demonstrated success with fabricating a 50 W/cm2 modular thermal protection system, a 400 C capable inflatable structure, a 10-meter scale aeroshell manufacturing capability, together with calibrated thermal and structural models. Despite the unquestionable success of the first HIAD ground development effort, it was recognized that additional investment was needed in order to realize the full potential of the HIAD technology capability to enable future flight opportunities. The second HIAD ground development effort will focus on extending performance capability in key technology areas that include thermal protection system, lifting-body structures, inflation systems, flight control, stage transitions, and 15-meter aeroshell scalability. This paper presents an overview of the accomplishments under the baseline HIAD development effort and current plans for a follow-on development effort focused on extending those critical technologies needed to enable a Mars Pathfinder mission.

Del Corso, Jospeh A.↗

Simplified analysis techniques to support the determination of Shuttle aerodynamics

The use of deterministic and regression methods to supplement the maximum likelihood estimation programs for the determination of the stability and control parameters of the Shuttle vehicle is discussed. Examples taken from STS-2 and STS-3 flight data show the use of deterministic methods to verify the effectiveness of the pitch and roll Reaction Control System (RCS), to verify the model structure used in the maximum likelihood estimator and to indicate regions of the descent trajectory where unmodeled moments dominate the vehicle motion and the application of maximum likelihood estimation would lead to incorrect parameter values. The use of regression methods to determine selected parameters when the more sophisticated maximum likelihood programs are not required, to determine the form of the vehicle mathematical model and to select the data time intervals for most effective parameter identification are also discussed.

Suit, W. T.↗

A Model-Model and Data-Model Comparison for the Early Eocene Hydrological Cycle

A range of proxy observations have recently provided constraints on how Earth's hydrological cycle responded to early Eocene climatic changes. However, comparisons of proxy data to general circulation model (GCM) simulated hydrology are limited and inter-model variability remains poorly characterised. In this work, we undertake an intercomparison of GCM-derived precipitation and P - E distributions within the extended EoMIP ensemble (Eocene Modelling Intercomparison Project; Lunt et al., 2012), which includes previously published early Eocene simulations performed using five GCMs differing in boundary conditions, model structure, and precipitation-relevant parameterisation schemes. We show that an intensified hydrological cycle, manifested in enhanced global precipitation and evaporation rates, is simulated for all Eocene simulations relative to the preindustrial conditions. This is primarily due to elevated atmospheric paleo-CO2, resulting in elevated temperatures, although the effects of differences in paleogeography and ice sheets are also important in some models. For a given CO2 level, globally averaged precipitation rates vary widely between models, largely arising from different simulated surface air temperatures. Models with a similar global sensitivity of precipitation rate to temperature (dP=dT ) display different regional precipitation responses for a given temperature change. Regions that are particularly sensitive to model choice include the South Pacific, tropical Africa, and the Peri-Tethys, which may represent targets for future proxy acquisition. A comparison of early and middle Eocene leaf-fossil-derived precipitation estimates with the GCM output illustrates that GCMs generally underestimate precipitation rates at high latitudes, although a possible seasonal bias of the proxies cannot be excluded. Models which warm these regions, either via elevated CO2 or by varying poorly constrained model parameter values, are most successful in simulating a match with geologic data. Further data from low-latitude regions and better constraints on early Eocene CO2 are now required to discriminate between these model simulations given the large error bars on paleoprecipitation estimates. Given the clear differences between simulated precipitation distributions within the ensemble, our results suggest that paleohydrological data offer an independent means by which to evaluate model skill for warm climates.

Boundary conditions↗

Hysteresis Modeling in Magnetostrictive Materials Via Preisach Operators

A phenomenological characterization of hysteresis in magnetostrictive materials is presented. Such hysteresis is due to both the driving magnetic fields and stress relations within the material and is significant throughout, most of the drive range of magnetostrictive transducers. An accurate characterization of the hysteresis and material nonlinearities is necessary, to fully utilize the actuator/sensor capabilities of the magnetostrictive materials. Such a characterization is made here in the context of generalized Preisach operators. This yields a framework amenable to proving the well-posedness of structural models that incorporate the magnetostrictive transducers. It also provides a natural setting in which to develop practical approximation techniques. An example illustrating this framework in the context of a Timoshenko beam model is presented.

Smith, R. C.↗

Geophysical data analysis and visualization using the Grid Analysis and Display System

Several problems posed by the rapidly growing volume of geophysical data are described, and a selected set of existing solutions to these problems is outlined. A recently developed desktop software tool called the Grid Analysis and Display System (GrADS) is presented. The GrADS' user interface is a natural extension of the standard procedures scientists apply to their geophysical data analysis problems. The basic GrADS operations have defaults that naturally map to data analysis actions, and there is a programmable interface for customizing data access and manipulation. The fundamental concept of the GrADS' dimension environment, which defines both the space in which the geophysical data reside and the 'slice' of data which is being analyzed at a given time, is expressed The GrADS' data storage and access model is described. An argument is made in favor of describable data formats rather than standard data formats. The manner in which GrADS users may perform operations on their data and display the results is also described. It is argued that two-dimensional graphics provides a powerful quantitative data analysis tool whose value is underestimated in the current development environment which emphasizes three dimensional structure modeling.

Doty, Brian E.↗

Formulation of full state feedback for infinite order structural systems

The estimation of exact displacement and displacement rate feedback kernels from finite dimensional control solutions based on finite element structural models is discussed. These kernels are then transformed to equivalent curvature and curvature rate feedback kernels. These curvature kernels are augmented with single point displacement and rotation feedback to account for rigid body motions. A growing class of sensors known as area-averaging sensors is used to measure the curvature and curvature rate state functions. The output of area-averaging sensors equals the convolution of all structural curvature states with the spatial sensitivity function of the sensors. Transforming the discrete feedback gains into continuous feedback kernels and employing area-averaging sensors make it possible to implement full state feedback for infinite order structural systems.

Miller, David W.↗

Modeling the Proton Radiation Belt With Van Allen Probes Relativistic Electron-Proton Telescope Data

An empirical model of the proton radiation belt is constructed from data taken during 2013-2017 by the Relativistic Electron-Proton Telescopes on the Van Allen Probes satellites. The model intensity is a function of time, kinetic energy in the range 18-600 megaelectronvolts, equatorial pitch angle, and L shell of proton guiding centers. Data are selected, on the basis of energy deposits in each of the nine silicon detectors, to reduce background caused by hard proton energy spectra at low L. Instrument response functions are computed by Monte Carlo integration, using simulated proton paths through a simplified structural model, to account for energy loss in shielding material for protons outside the nominal field of view. Overlap of energy channels, their wide angular response, and changing satellite orientation require the model dependencies on all three independent variables be determined simultaneously. This is done by least squares minimization with a customized steepest descent algorithm. Model uncertainty accounts for statistical data error and systematic error in the simulated instrument response. A proton energy spectrum is also computed from data taken during the 8 January 2014 solar event, to illustrate methods for the simpler case of an isotropic and homogeneous model distribution. Radiation belt and solar proton results are compared to intensities computed with a simplified, on-axis response that can provide a good approximation under limited circumstances.

Selesnick, R. S.↗

Stream Flow Prediction by Remote Sensing and Genetic Programming

A genetic programming (GP)-based, nonlinear modeling structure relates soil moisture with synthetic-aperture-radar (SAR) images to present representative soil moisture estimates at the watershed scale. Surface soil moisture measurement is difficult to obtain over a large area due to a variety of soil permeability values and soil textures. Point measurements can be used on a small-scale area, but it is impossible to acquire such information effectively in large-scale watersheds. This model exhibits the capacity to assimilate SAR images and relevant geoenvironmental parameters to measure soil moisture.

Chang, Ni-Bin↗

Comparison of Model Calculations of Biological Damage from Exposure to Heavy Ions with Measurements

The space environment consists of a varying field of radiation particles including high-energy ions, with spacecraft shielding material providing the major protection to astronauts from harmful exposure. Unlike low-LET gamma or X rays, the presence of shielding does not always reduce the radiation risks for energetic charged-particle exposure. Dose delivered by the charged particle increases sharply at the Bragg peak. However, the Bragg curve does not necessarily represent the biological damage along the particle path since biological effects are influenced by the track structures of both primary and secondary particles. Therefore, the ''biological Bragg curve'' is dependent on the energy and the type of the primary particle and may vary for different biological end points. Measurements of the induction of micronuclei (MN) have made across the Bragg curve in human fibroblasts exposed to energetic silicon and iron ions in vitro at two different energies, 300 MeV/nucleon and 1 GeV/nucleon. Although the data did not reveal an increased yield of MN at the location of the Bragg peak, the increased inhibition of cell progression, which is related to cell death, was found at the Bragg peak location. These results are compared to the calculations of biological damage using a stochastic Monte-Carlo track structure model, Galactic Cosmic Ray Event-based Risk Model (GERM) code (Cucinotta, et al., 2011). The GERM code estimates the basic physical properties along the passage of heavy ions in tissue and shielding materials, by which the experimental set-up can be interpreted. The code can also be used to describe the biophysical events of interest in radiobiology, cancer therapy, and space exploration. The calculation has shown that the severely damaged cells at the Bragg peak are more likely to go through reproductive death, the so called "overkill".

Kim, Myung-Hee Y.↗

Aeroelastic loads and stability investigation of a full-scale hingeless rotor

An analytical investigation was conducted to study the influence of various parameters on predicting the aeroelastic loads and stability of a full-scale hingeless rotor in hover and forward flight. The CAMRAD/JA (Comprehensive Analytical Model of Rotorcraft Aerodynamics and Dynamics, Johnson Aeronautics) analysis code is used to obtain the analytical predictions. Data are presented for rotor blade bending and torsional moments as well as inplane damping data obtained for rotor operation in hover at a constant rotor rotational speed of 425 rpm and thrust coefficients between 0.0 and 0.12. Experimental data are presented from a test in the wind tunnel. Validation of the rotor system structural model with experimental rotor blade loads data shows excellent correlation with analytical results. Using this analysis, the influence of different aerodynamic inflow models, the number of generalized blade and body degrees of freedom, and the control-system stiffness at predicted stability levels are shown. Forward flight predictions of the BO-105 rotor system for 1-G thrust conditions at advance ratios of 0.0 to 0.35 are presented. The influence of different aerodynamic inflow models, dynamic inflow models and shaft angle variations on predicted stability levels are shown as a function of advance ratio.

Peterson, Randall L.↗

GRAIL Refinements to Lunar Seismic Structure

Gravity field measurements are perhaps the most numerous of the indirect observations relevant to the Moon's internal structure. Multiple recent missions have mapped the global lunar gravity field, each one improving upon the resolution of the last. The details of the Moon's deepest structure, including the parameters that define the lunar core, however, were still largely unaddressed by pre-GRAIL gravity measurements, which were not high enough resolution to resolve the tidal coefficients at a sufficient degree of accuracy. Current constraints on core size and state arise from other indirect measurements, including lunar laser ranging, magnetic induction studies, and analyses of elemental abundances in depth-derived mare basalts. These inferences vary widely, but when considered together with structure models derived from the seismic data gathered during the Apollo missions, a schematic of the lunar interior containing a partially molten deepest mantle layer overlying molten outer and solid inner core layers was obtained. Seismology provides the most direct constraints on the variables that govern the dynamic properties of the body. However, the GRAIL mission's high-resolution measurements of the lunar gravity field are being used to constrain the interior structure of the Moon using a "crust to core" approach. GRAIL's constraints on crustal thickness, mantle structure, core radius and stratification, and core state (solid vs. molten) therefore complement seismic investigations. This work focuses on expanding our knowledge of the Moon's internal structure using joint gravity and seismic analyses, which will improve constraints on the deep lunar mantle and core.

Weber, Renee C.↗

A comparison of measured and calculated thermal stresses in a hybrid metal matrix composite spar cap element

A hybrid spar of titanium with an integrally brazed composite, consisting of an aluminum matrix reinforced with boron-carbide-coated fibers, was heated in an oven and the resulting thermal stresses were measured. Uniform heating of the spar in an oven resulted in thermal stresses arising from the effects of dissimilar materials and anisotropy of the metal matrix composite. Thermal stresses were calculated from a finite element structural model using anisotropic material properties deduced from constituent properties and rules of mixtures. Comparisons of calculated thermal stresses with measured thermal stresses on the spar are presented. It was shown that failure to account for anisotropy in the metal matrix composite elements would result in large errors in correlating measured and calculated thermal stresses. It was concluded that very strong material characterization efforts are required to predict accurate thermal stresses in anisotropic composite structures.

Jenkins, J. M.↗

The 15-meter diameter mechanically scanned deployable antenna

A preliminary design with structural model data and thermal-performance estimates of a 15-meter mechanically scanned deployable antenna (MSDA) that could be launched onboard a Shuttle Orbiter to provide radiometric brightness temperature maps of the Earth and oceans in selected bands over a frequency range from 1.4 to 11 GHz is provided. The study objectives were met through the design of a unique, integrated, offset feed mast and reflector design that uses the deployable box-truss structure as a building block. The performance of this system is summarized. The all graphite-epoxy, 4.57-meter prototype cube that was completed in 1981 and is proposed for this reflector and feed mast design is presented.

Coyner, J. V.↗