ALBERT

Description: The objectives of this program are as follows: modelling of guided waves in fiber-reinforced plates in terms of different modes; and analysis of scattering by transverse cracks using modal representation. A hybrid numerical method combining the finite element representation of a region around the crack with the modal representation in the exterior region will be used in this program. Modes will be obtained using the through-the-thickness discretization of the displacement field.

Description: The Multispectral Atmospheric Mapping Sensor (MAMS) is an airborne scanner which collects imagery in eight visible and three IR bands, with ground resolution of 100 m at a nominal aircraft altitude of 20 km. The visible channels provide continuous spectral coverage from 0.42 to 1.05 microns, and the IR channels measure upwelling and reflected radiation at combinations of 3.7, 6.5, 11.1, and 12.5 microns. These include channels nearly identical to the Landsat TM channels 1-4, and similar to those on the VISSR Atmospheric Sounder (VAS) and AVHRR instruments. The high spatial resolution permits investigation of surface thermal and atmospheric water vapor structure at scales not available from current stabilities. Significant modifications have been made since engineering flights in 1985 to improve the radiometric performance. The increased data quality and recent multidisciplinary applications of this data are presented.

Description: In summer 1988 an intensive field experiment was conducted in the vicinity of the Pisgah lava flow in the Mojave Desert. Physical properties such as microtopography, composition, soil moisture, and dielectric constant were measured at five sites representing surfaces with rms heights varying from less than one centimeter to tens of centimeters. In addition, polarimetric radar images at P-band, L-band, and C-band were acquired at three different incidence angles with the NASA/JPL airborne imaging radar polarimeter. Using trihedral corner reflectors deployed in the area prior to imaging, the L- and C-band images were calibrated to provide sigma0 values for each resolution element in the scene. The results of inferring geophysical parameters such as rms surface height and correlation length of the measured surfaces by fitting the observed signatures with those predicted by the small-perturbation model are presented: for smoother surfaces, the rms height values inferred are in good agreement with in situ measurements.

Description: Radar modeling of mangal forest stands in the Sundarbans area of southern Bangladesh is discussed. The modeling uses radar system parameters with forest data on tree height, spacing, biomass, species combinations, and water content (including slightly conductive water) in leaves and trunks of the mangal. For Sundri and Gewa tropical mangal forests, six model components are proposed; the models are required to explain the contributions of various combinations of forest species in the attenuation and scattering of mangal-vegetated nonflooded or flooded surfaces. Statistical data of simulated images have been compared with those of SIR-B images both to refine the modeling procedures and to characterize the model output appropriately. The possibility of delineation of flooded or nonflooded boundaries is discussed.

Description: Surface-air temperature differences are commonly used in a bulk resistance equation for estimating sensible heat flux (H), which is inserted in the one-dimensional energy balance equation to solve for the latent heat flux (LE) as a residual. Serious discrepancies between estimated and measured LE have been observed for partial-canopy-cover conditions, which are mainly attributed to inappropriate estimates of H. To improve the estimates of H over sparse canopies, one- and two-layer resistance models that account for some of the factors causing poor agreement are developed. The utility of the two models is tested with remotely sensed and micrometeorological data for a furrowed cotton field with 20 percent cover and a dry soil surface. It is found that the one-layer model performs better than the two-layer model when a theoretical bluff-body correction for heat transfer is used instead of an empirical adjustment; otherwise, the two-layer model is better.

Description: Experimental treatments in a Douglas-fir forest in NE New Mexico were carried out to determine whether differences in forest canopy chemistry could be detected using data from the Airborne Imaging Spectrometer (AIS-2). Experimental treatments consisted of nitrogen fertilizer additions, sawdust additions, and control plots. After AIS-2 data were collected, the digital number of a given pixel was extracted from each channel, yielding 128 values that were used to form a spectrum. Four spectra were extracted from each treatment plot. Multiple stepwise linear regressions between first and second difference transformations of AIS-2 spectra and the canopy characteristics of biomass, nitrogen concentration, and nitrogen content were performed. The results showed a coefficient of multiple determination of 0.71 between first-difference AIS-2 spectra and measured nitrogen concentration in foliage, indicating that it may be possible to predict nitrogen concentration in Douglas fir using AIS-2 spectra.

Description: The information content in polarimetric SAR images is examined, and the polarimetric image variables containing the information that is important to the classification of terrain features in the images are determined. It is concluded that accurate classification can be done when just over half of the image variables are retained. A reduction in image data dimensionality gives storage savings, and can lead to the improvement of classifier performance. In addition, it is shown that a simplified radar system with only phase-calibrated CO-POL or SINGLE TX channels can give classification performance which approaches that of a fully polarimetric radar.

Description: The spatial and phenological distribution of vegetation was examined in a remote sensing geobotanical study of the East Africa Rift region. Six normalized difference vegetation index (NDVI) AVHRR (advanced very high resolution radiometer) scenes (February, June, and September of both 1984 and 1987) were used as a measure of vegetation presence. NOAA ETOPO-5 elevation and geographic coordinate data were coregistered with the multitemporal NDVI images. Univariate and multivariate statistics indicate that the NDVI values are significantly associated with elevation, latitude, and longitude. This supports the concept of quantifiable, regional gradients that affect large-scale geobotanical studies. A quadratic regression line was fitted to the NDVI, elevation, latitude, and longitudinal data. In this way a regional trend, as affected by elevation, was determined. The deviation of the actual data from this regional trend was displayed as an image, and shows the local climatic and geological influences.

Description: The statistical characteristics of image speckle are reviewed. Existing segmentation techniques that have been used for speckle filtering, edge detection, and texture extraction are sumamrized. The relative effectiveness of each technique is briefly discussed.

Description: Theory and experiments have shown that passive microwave radiometers can be used to measure soil moisture. However, the presence of a vegetative cover alters the measurement that might be obtained under bare conditions. Deterministically accounting for the effect of vegetation and developing algorithms for extracting soil moisture from observations of a vegetable-soil complex present significant obstacles to the practical use of this approach. The presence of a vegetation canopy reduces the sensitivity of passive microwave instruments to soil moisture variations. The reduction in sensitivity, as compared to a bare-soil relationship, increases as microwave frequency increases, implying that the longest wavelength sensors should provide the most information. Sensitivity also decreases as the amount of vegetative wet biomass increases for a given type of vegetation.

Description: Spectral bidirectional reflectances were measured over three natural soil sites using a specially designed radiometer called the Parabola. Two of the sites were bare soils, and the third had a sparse cover of desert scrub. The reflectances were strongly non-Lambertian for all three surfaces, but with markedly different patterns. The measured data were fitted with a quasi-physical reflectance model in which the surface backscattering and forwardscattering are separately formulated. A soil reflectance characterization was obtained by assessing the contributions of the forward, backward, and Lambertian components. This three-parameter characterization produced a satisfactory fit to the measured reflectances and appears promising as a basis for soils categorization.

Description: A radiative-transfer model for the bidirectional reflectance-distribution function (BRDF) shows that snow is moderately anisotropic in the NIR wavelengths. Although the directional-hemispherical albedo of snow decreases as the grains become larger, the forward scattering also increases, with the result that the illumination and viewing geometry must be considered when interpreting physical properties of the surface layer of the snow pack from remote sensing data. Measurements of the BRDF and the transmittance for a variety of snow conditions were made throughout the winter and spring seasons with an SE-590 spectroradiometer, for wavelengths from 0.38 to 1.11 microns. Coincident with these, the surface grain properties were analyzed by stereological methods. A sphere with the same surface-to-volume ratio as the ice grains is used as the equivalent sphere.

Description: Remotely sensed data are being used to estimate foliar chemical content. This paper reviews how stepwise multiple regression and deconvolution have been used to extract chemical information from foliar spectra, and concludes that both methods are useful, but neither is ideal. It is recommended that the focus of research be modeling in the long term and experimentation in the short term. Long-term research should increase our understanding of the interaction between radiation and foliar chemistry so that the focus of research can move from leaf model to canopy model to field experiment. Short-term research should aim to design experiments in which remotely sensed data are used to generate unambiguous and accurate estimates of foliar chemical content.

Description: Surface energy balance components were evaluated by combining satellite-based spectral data with on-site measurements of solar irradiance, air temperature, wind speed, and vapor pressure. Maps of latent heat flux density and net radiant flux density were produced using Landsat TM data for three dates. The TM-based estimates differed from Bowen-ratio and aircraft-based estimates by less than 12 percent over mature fields of cotton, wheat, and alfalfa.

Description: Multipolarization P-, L-, and C-band airborne SAR data sets were acquired over a coastal zone and a forested wetland of southern Louisiana. The data sets were used with field-collected surface-parameter data in order to determine the value of SAR systems in assessing and mapping coastal-zone surface features. The coastal-zone surface features in this study are sediments, sediment distribution, and the formation of new isles and banks. Results of the data analysis indicate that the P-band radar with 68-cm wavelength is capable of detecting the submerged sediment if the area is very shallow (i.e., a water depth of less than one meter). The penetration capability of P-band radar is also demonstrated in the forested wetland area. The composition and condition of the ground surface can be detected, as well as the standing water beneath dense tree leaves.

Description: This paper considers complex transcendental eigenvalue problems where one is interested in pairs of eigenvalues that are restricted to take real values only. Such eigenvalue problems arise in dynamic stability analysis of nonconservative physical systems, i.e., flutter analysis of aeroelastic systems. Some available solution methods are discussed and a new method is presented. Two computational approaches are described for analytical evaluation of the sensitivities of these eigenvalues when they are dependent on other parameters. The algorithms presented are illustrated through examples.

Description: Multiband radiometric data from an airborne imaging thermal scanner are being studied for use in finding buried gravel deposits. The techniques are based on measuring relative differences in the thermal properties between gravel-laden targets and the surrounding gravelless background. These properties are determined from modeling the spectral radiant emittance recorded over both types of surfaces in conjunction with ground measurements of the most significant heat flows above and below the surface. Thermodynamic properties of sampled materials from control sites are determined, and diurnal and annual subsurface heat waves are recorded. Thermal models that account for heat exchange at the surface, as well as varying levels of soil moisture, humidity, and vegetation, are needed for adaptation and modification to simulate the physical and radiative environments of this region.

Description: The advanced solid-state array spectroradiometer (ASAS) is an airborne, off-nadir pointing imaging spectroradiometer used to acquire bidirectional radiance data for terrestrial targets. As its platform aircraft flies over a target the sensor can image the target through a sequence of at least seven fore-to-aft view directions ranging up to 45 deg on either side of nadir. ASAS acquires data for 29 spectral bands in the visible and near-infrared portions of the spectrum with a resolution of 15 nm. ASAS data were recently acquired for a prairie ecosystem and a northern forest ecosystem. The data demonstrate the combined effects of reflectance anisotropy and increased atmospheric path length on off-nadir observations. One result of these effects is a variation in vegetation indices as a function of view direction. Normalized-difference-vegetation-indices for prairie grass, coniferous, and deciduous canopies varied up to 14 percent, 23 percent, and 6 percent, respectively, relative to nadir as a function of view zenith angle along the solar principal plane.

Description: During the First International Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE), 80 gigabytes of image data were generated from a variety of satellite and airborne sensors in a multidisciplinary attempt to study energy and mass exchange between the land surface and the atmosphere. To make these data readily available to researchers with a range of image data handling experience and capabilities, unique image-processing software was designed to perform a variety of nonstandard image-processing manipulations and to derive a set of standard-format image products. The nonconventional features of the software include: (1) adding new layers of geographic coordinates, and solar and viewing conditions to existing data; (2) providing image polygon extraction and calibration of data to at-sensor radiances; and, (3) generating standard-format derived image products that can be easily incorporated into radiometric or climatology models. The derived image products consist of easily handled ASCII descriptor files, byte image data files, and additional per-pixel integer data files (e.g., geographic coordinates, and sun and viewing conditions). Details of the solutions to the image-processing problems, the conventions adopted for handling a variety of satellite and aircraft image data, and the applicability of the output products to quantitative modeling are presented. They should be of general interest to future experiment and data-handling design considerations.

Description: In order to evaluate the importance of surface thermal flux and emissivity variations on surface and boundary layer processes, a technique that uses thermal data from an airborne multispectral scanner to determine the surface skin temperature and thermal emissivity over a regional area has been developed. These values are used to estimate the total flux density emanating from the surface and at the top of the atmosphere. Data from the multispectral atmospheric mapping sensor (MAMS) collected during the First ISLSCP Field Experiment (FIFE) are used to develop the technique, and to show the time and space variability of the flux values. The ground truth data available during FIFE provide a unique resource to evaluate this technique.

Description: A common problem in acquiring ground-truth data for use in microwave interaction modeling is the capture of surface roughness data that are both sampled at distances comparable to a fraction of the wavelength and extensive enough to represent the surface statistics in at least one resolution cell of the microwave remote sensor used. A technique has been developed for acquiring the necessary photogrammetric data using twin 70-mm film cameras mounted on a helicopter boom. The apparatus is described, and the accuracy with which ground surface roughness can be characterized using this device is estimated. Both standard and cross-correlation methods were used for data reduction. Stereogrammetry is compared with a completely automated image-matching technique. Dense disparity images were generated from the helicopter stereo pairs. Using interior orientation parameters supplied by the camera manufacturers, and assuming that exterior orientation parameters remain constant between control target and test field photography, an extremely dense DEM (digital elevation model) for a test field has been derived. Results are compared, and accuracy estimates are presented.

Description: Algorithms to correct for atmospheric-scattering effects in high-spatial resolution land-surface images require the ability to perform rapid and accurate computations of the top-of-atmosphere diffuse radiance field for arbitrarily general surface reflectance distributions (which may be both heterogeneous and non-Lambertian) and atmospheric models. Using three-dimensional radiative transfer (3DRT) theory algorithms are being developed. The methodology used to perform the 3DRT calculations is described. It is shown how these calculations are used to perform atmospheric corrections, and the sensitivity of the retrieved surface reflectances to atmospheric structural parameters is illustrated.

Description: A model for the brightness temperature of a canopy which is based on the physical characteristics of the vegetation is being developed. This model uses the Peake model to relate the scattering coefficient of the vegetation to its radiometric brightness temperature and a discrete scattering model to computer the scattering coefficient. A distorted Born approximation is used to provide a tractable solution to the scattering problem. To verify the theory a comparison has been made with experiments. Radiometric measurements of brightness temperature were made at L- and C-bands over a crop of mature soybeans, along with ground truth measurements. All model input parameters were measured; there were no free parameters to adjust to obtain agreement with the measurements. Results showed that agreement between theory and measurements was reasonable for horizontal polarization at both L- and C-bands with dry underlying soil.

Description: During the spring of 1988, the NASA/JPL multifrequency, multipolarization SAR (synthetic aperture radar) flew in series of calibration experiments over the Goldstone dry lake bed. An array of calibration devices was deployed, including dihedral and trihedral corner reflectors, polarimetric active radar calibrators, passive receivers, and CW (continuous wave) tone generators. The aim of the campaign was to calibrate both amplitude and phase of the resulting SAR images, over long and short time scales. The results of the analysis of the calibration data collected in the spring of 1988 are presented. Trihedral corner reflector signatures and certain image background measures are used to externally calibrate relative amplitude and phase between polarizations at a given frequency, and to calibrate across frequencies. Assessments are made of the calibration accuracy as a function of image frame position for each frequency, and the stability of the radar calibration over long and short time-scales.

Description: An unsupervised clustering algorithm is used to segment multifrequency polarimetric radar data from the NASA/JPL airborne SAR (synthetic aperture radar). Twenty-two parameters are evaluated for their discriminatory capability for each pixel of an image. A clustering analysis is then performed using different subsets of these parameters. This analysis relies on data taken as part of an intensive field experiment during the summer of 1988 in the vicinity of the Pisgah lava flow in the Mojave Desert in southern California. As part of the experiment, extensive ground truth was acquired, including dielectric constant and topography measurements. Segmentation results show good agreement with these measurements.

Description: Delamination growth caused by local buckling of a delaminated group of plies was investigated. Delamination growth was assumed to be governed by the strain energy release rates, G(1), G(2) and G(3). The strain energy release rates were calculated using a geometrically nonlinear three-dimensional finite element analysis. The program is described and several checks of the analysis are discussed. Based on a limited parametric study, the following conclusions were reached: (1) the problem is definitely mixed mode (in some cases G(1) is larger than G(2), for other cases the opposite is true); (2) in general, there is a large gradient in the strain energy release rates along the delamination front; (3) the locations of maximum G(1) and G(2) depend on the delamination shape and the applied strain; (4) the mode 3 component was negligible for all cases considered; and (5) the analysis predicted that parts of the delamination would overlap. The results presented did not impose contact constraints to prevent overlapping. Further work is needed to determine the effects of allowing the overlapping.

Description: The propagation of elastic waves of the type induced by impact are studied for a special class of cylinders. In this class are thin finite cylinders, for which Mindlin's equations are appropriate, on which point masses are attached. An approximation is derived valid for frequencies below cutoff of the lowest shear mode. An eigenfunction expansion is used to compute the transient response of the coupled system.

Description: The method of boundary elements is adapted to the dynamics of elastic axisymmetric solids. Three scalar potentials are used, each satisfying the Helmholtz equation and each represented by its own source distribution. The sources are distributed over a surface enclosing the solid boundary. Stress or displacement conditions on the boundary are realized by third-order derivatives of the potentials. Discrete ring elements on the source enclosure combine their influences at discrete ring elements on the solid boundary.

Description: An experiment to study the hypothesis that winter cold fronts are more important that violent summer storms in moving sediment in midlatitude coastal river deltas. The experiment is based on Calibrated Airborne Multispectral Scanner (CAMS) data. A new approach for the georeferencing of aircraft imagery has been developed for use in this experiment. The software for the approach uses the aircraft inertial navigation data and is able to automatically map a three-channel 1600 line file through a 90 deg rotation in less than 14 minutes. Results from tests of the technique show rms errors of less than 2/3 of a pixel.

Description: While the standard mode acceleration formulation in structural dynamics has often been interpreted to suggest that the reason for improved convergence obtainable is that the dynamic correction factor is divided by the modal frequencies-squared, an alternative formulation is presented which clearly indicates that the only difference between mode acceleration and mode displacement data recovery is the addition of a static correction term. Attention is given to the advantages in numerical implementation associated with this alternative, as well as to an illustrative example.

Description: A key step in the application of formal automated design techniques to structures under transient loading is the calculation of sensitivities of response quantities to the design parameters. This paper considers response quantities to the design parameters. This paper considers structures with general forms of damping acted on by general transient loading and addresses issues of computational errors and computational efficiency. The equations of motion are reduced using the traditional basis of vibration modes and then integrated using a highly accurate, explicit integration technique. A critical point constraint formulation is used to place constraints on the magnitude of each response quantity as a function of time. Three different techniques for calculating sensitivities of the critical point constraints are presented. The first two are based on the straightforward application of the forward and central difference operators, respectively. The third is based on explicit differentiation of the equations of motion. Condition errors, finite difference truncation errors, and modal convergence errors for the three techniques are compared by applying them to a simple five-span-beam problem. Sensitivity results are presented for two different transient loading conditions and for both damped and undamped cases.

Description: A new and simple method of finite-element grid improvement is presented. The objective is to improve the accuracy of the analysis. The procedure is based on a minimization of the trace of the stiffness matrix. For a broad class of problems this minimization is seen to be equivalent to minimizing the potential energy. The method is illustrated with the classical tapered bar problem examined earlier by Prager and Masur. Identical results are obtained.

Description: Space photography has been successfully used to extend the space remote sensing data base for environmental monitoring by a decade. In this study of Lake Chad, space photographs were digitized and registered to a topographic base map before water classifications were performed. From 1966 to 1985 over a 21,000-square-kilometer decrease in lake surface area was observed.

Description: Space Shuttle astronauts have used handheld cameras to take about 30,000 photographs of the earth as seen from orbit. These pictures provide valuable, true-color depictions of many geologically significant areas. While the photographs have areal coverages and resolutions similar to the more familiar Landsat MSS and TM images, they differ from the latter in having a wide variety of solar illumination angles and look angles. Astronaut photographs can be used as very small scale aerial photographs for geologic mapping and planning logistical support for field work. Astronaut photography offers unique opportunities, because of the intelligence and training of the on-orbit observer, for documenting dynamic geologic activity such as volcanic eruptions, dust storms, etc. Astronauts have photographed more than 3 dozen volcanic eruption plumes, some of which were not reported otherwise. The stereographic capability of astronaut photography also permits three-dimensional interpretation of geologic landforms which is commonly useful in analysis of structural geology. Astronauts have also photographed about 20 known impact craters as part of project to discover presently unknown examples in Africa, South America, and Australia.

Description: Assuming a single-mode transverse displacement, a simple formula is derived for the transverse load-displacement relationship of a plate under in-plane compression. The formula is used to derive a simple analytical expression for the nonlinear dynamic response of postbuckled plates under sinusoidal or random excitation. The highly nonlinear motion of snap-through can be easily interpreted using the single-mode formula. Experimental results are obtained with buckled and cylindrical aluminum panels using discrete frequency and broadband excitation of mechanical and acoustic forces. Some important effects of the snap-through motion on the dynamic response of the postbuckled plates are described. Static tests were used to identify the deformation shape during snap-through.

Description: An experimental investigation of composite to composite bonded joints was undertaken to study the effect of bond thickness on debond growth rate under cyclic loading and critical strain energy release rate under static loading. Double cantilever beam specimens of graphite/epoxy adherends bonded with EC 3445 were tested under mode I loading. A different behavior of fracture and fatigue strength was observed with variation of bondline thickness.

Description: The environmental, geologic, meteorologic, and oceanographic phenomena documented by earth photography during the Space Shuttle STS-29 mission are reviewed. A map of the nadir point positions of earth-viewing photographs from the mission is given and color photographs of various regions are presented. The mission photographs include atmospheric dust and smoke over parts of Africa and Asia, Sahelian water sites, center pivot irrigation fields in the Middle East, urban smog over Mexico City, isolated burning in the Bolivian Amazon, and various ocean features and cloud formations.

Description: The feasibility of using a viscoplastic model developed by Robinson to perform a nonlinear structural analysis was investigated. The paper presents analytical solutions for three classical problems: (1) a pressurized thick-walled cylinder, (2) a thin rotating disk, and (3) a pressurized thick-walled sphere. The analytical expressions derived for the stress and the strain rates for these components are general in nature as they consider both the mechanical and thermal loadings to be time-dependent. A computer program VISTAN (VIscoplastic STructural ANalyzer) was developed to obtain the stress and strain distributions. The finite element solutions for these problems are also presented. The numerical results pertaining to isothermal loading conditions are provided. The results obtained demonstrate the feasibility of using the viscoplastic model developed by Robinson to perform nonlinear structural analyses.

Description: Remote sensing instruments for monitoring global changes are examined. The use of the earth observing system, a set of instrument platforms in polar, sun-synchoronous orbit that provide coverage of the entire globe, is discussed. The radar and imaging spectrometers utilized to obtain surface measurements are described. Atmospheric data is collected by the atmospheric IR sounder, the tropospheric emission spectrometer, and the stratospheric wind IR limb sounder. Consideration is given to the operation of the microwave limb sounder, the active cavity radiometer, and the TDRSS.

Description: A combined experimental and analytical study of bearing and bypass loading on single-fastener specimens of a 16-ply, quasi-isotropic T300/5208 graphite/epoxy laminate with a centrally located hole is reported. The specimens were loaded in either tension or compression, and onset damage, ultimate strengths, and corresponding failure modes were determined. The tension data showed the expected linear interaction for combined bearing-bypass loading with damage developing in the net-section tension mode. However, the bearing-onset strengths showed an unexpected interaction of the bearing and compressive bypass loads in which the latter reduced the bearing-onset strength. A linear finite element analysis showed that bearing-bypass loading had a marked influence on the bolt-hole contact which in turn had a significant effect on local stresses.

Description: Multispectral imaging systems provide much of the basic data used by the land and ocean civilian remote-sensing community. There are numerous multispectral imaging systems which have been and are being developed. A common way to compare the radiometric performance of these systems is to examine their noise-equivalent change in reflectance, NE Delta-rho. The NE Delta-rho of a system is the reflectance difference that is equal to the noise in the recorded signal. A comparison is made of the noise equivalent change in reflectance of seven different multispectral imaging systems (AVHRR, AVIRIS, ETM, HIRIS, MODIS-N, SPOT-1, HRV, and TM) for a set of three atmospheric conditions (continental aerosol with 23-km visibility, continental aerosol with 5-km visibility, and a Rayleigh atmosphere), five values of ground reflectance (0.01, 0.10, 0.25, 0.50, and 1.00), a nadir viewing angle, and a solar zenith angle of 45 deg.

Description: Analysis of SAR data of the Greenland ice sheet in summer and winter suggest the use of SAR to monitor the temporal hydrology of ice sheets. Comparisons of each SAR data set with summer Landsat TM imagery show an areal-positive correlation with summer SAR data and a negative correlation with winter SAR data. It is proposed that the summer SAR data are most sensitive to the variable concentrations of free water in the surface snow and that the winter SAR data indicate variations in snow grain size.

Description: The NASA Ames Ecosystem Science and Technology Branch and the U.S. Army Medical Research Institute of Infectious Diseases are conducting research to detect Rift Valley fever (RVF) vector habitats in eastern Africa using active and passive remote-sensing. The normalized difference vegetation index (NDVI) calculated from Landsat TM and SPOT data is used to characterize the vegetation common to the Aedes mosquito. Relationships have been found between the highest NDVI and the 'dambo' habitat areas near Riuru, Kenya on both wet and dry data. High NDVI values, when combined with the vegetation classifications, are clearly related to the areas of vector habitats. SAR data have been proposed for use during the rainy season when optical systems are of minimal use and the short frequency and duration of the optimum RVF mosquito habitat conditions necessitate rapid evaluation of the vegetation/moisture conditions; only then can disease potential be stemmed and eradication efforts initiated.

Description: The possibility of using detailed spatial soil moisture maps as input to a runoff model was investigated. The water balance of a small drainage basin was simulated using a simple storage model. Aircraft microwave measurements of soil moisture were used to construct two-dimensional maps of the spatial distribution of the soil moisture. Data from overflights on different dates provided the temporal changes resulting from soil drainage and evapotranspiration. The study site and data collection are described, and the soil measurement data are given. The model selection is discussed, and the simulation results are summarized. It is concluded that a time series of soil moisture is a valuable new type of data for verifying model performance and for updating and correcting simulated streamflow.

Description: Multiwavelength, multipolarization radar image data can be processed to allow the clear depiction of most information in color composite displays. The first principal component, depicted as variations in image intensity, captures the dominant variation in the data and reduces noise without degradation of spatial detail. Other principal components, depicted primarily as variations in image chromaticity, are derived from spatially filtered data so that image colors can be attributed to significant wavelength and polarization contrasts, and not to speckle. Processing paths can be designed to emphasize wavelength contrasts, polarization contrasts, or overall discrimination of surface features.

Description: An analysis is made of time-dependent covariance of the greenness vegetation index with mapped terrain variables over the Konza Prarie (Kansas) during the 1987 growing season. The analysis was part of an ongoing project to establish appopriate ground-sampling and data-integration strategies for satellite-based monitoring of land surface climate conditions. Greenness images for six dates between May and October were derived from atmospherically corrected thematic mapper (TM) data and coregistered with maps of woody vegetation, fire, and soils. Local variance in greenness peaked in mid-June, falling rapidly until mid-August, and declining gradually thereafter. Greenness images exhibited positive autocorrelation up to distances of 180-210 m, but the dominant scale of pattern occurred at a block size of 60 m by 60 m throughout the growing season. 40-44 percent of total scene variance in July and August was accounted for by the effects of woody vegetation (8.9 percent of the area), prairie burning, and soil type. The effect of these terrain variables was fairly consistent between June and late August and was manifested as additional high-frequency spatial variation in imagery from that period.

Description: The Alaska SAR (synthetic-aperture radar) Facility (ASF) will be capable of receiving, processing, archiving, and producing a variety of SAR image products from three satellite-borne SARs: E-ERS-1 (ESA), J-ERS-1 (NASDA) and Radarsat (Canada). Crucial to the success of the ASF is the Alaska SAR processor (ASP), which will be capable of processing over 200 100-km x 100-km (Seasat-like) frames per day from the raw SAR data, at a ground resolution of about 30 m x 30 m. The processed imagery is of high geometric and radiometric accuracy, and is geolocated to within 500 m. Special-purpose hardware has been designed to execute a SAR processing algorithm to achieve this performance. This hardware is currently undergoing acceptance testing for delivery to the University of Alaska. Particular attention has been devoted to making the operations semi-automated and to providing a friendly operator interface via a computer workstation. The operations and control of the Alaska SAR processor are described.

Description: The problem of estimating areal averages from point measurement has been extensively studied by mining engineers and hydrologists. Its application to satellite measurements has recently been introduced. The semivariaogram has been used in many geostatistical applications to estimate spatial structures of observed properties, such as mineral distributions. An examination is made of snow variations in Colorado from daily snow data collected in 11 SNOTEL stations. The associated semivariogram is estimated. The objective is to estimate the spatial structure of the snow field so that the point data can be used for comparison with, and validation for, satellite measurements.

Description: Spaceborne observations have become increasingly accepted for studies of the hydrological cycle and as a source of information that can be used for management and assessment of water resources over large areas. As data records are extended over long periods of time, using sensors that improve their reliability and quality, the utility of spaceborne observations will become increasingly evident. The Earth Observing System (EOS) concept being studied by NASA holds potential for providing improved means and procedures for hydrological studies, particularly on a global basis, using spaceborne sensors and a comprehensive information system. Continuing emphasis needs to be placed on research exploring the interpretation of remotely sensed observations through the assimilation and use of these data in models of hydrological processes.

Description: A study undertaken to determine SIR-C (Shuttle Imaging Radar-C) science requirements for polarimetric calibration is described. A model describing the effect of miscalibration is presented, followed by a description of specific experiments and the resulting requirements levied on the system. The effects of miscalibration on some commonly used measures which can be extracted from SAR (synthetic aperture radar) data are also discussed. It is concluded that polarimetric calibration requirements are strongly application-dependent, and that a compromise must be made to maximize the domain of application of the data without prohibitive cost. Data summarizing SIR-C polarimetric calibration goals are also presented.

Description: Regression and ratio estimators are used to integrate AVHRR-Global Area Coverage (GAC) and Landsat MSS digital data to estimate forest area in the continental United States. Forestlands are enumerated for the 48 contiguous states using five different AVHRR-GAC data sets. Results indicated that the GAC and MSS forest estimates were not highly correlated. Although the ratio of means and linear regression corrections were, on the average, closer to national U.S. Forest Service forest area estimates, these correction procedures did not consistently improve GAC estimates of forest area. GAC forest area estimates tended to be high in densely forested regions such as the northeast and low in sparsely forested areas.

Description: Consideration is given to the possibility of using Landsat TM data to detect differences and changes in the leaf area index (LAI) of closed canopy pine plantations in central Massachusetts. Field measurements of the LAI were obtained using locally-derived allometric relationships between leaf area and trunk diameter. The relationship between the Landsat TM Band 3/Band 4 ratio and regional differences in LAI is examined. A thinning treatment was conducted to demonstrate the ability of TM to detect a known reduction in LAI. It is concluded that, for the detection of moderate changes and differences in the LAI of closed canopy pine plantations at local scales, the TM sensor is better than field measurements involving allometric equations.

Description: Results are presented from field experiments relating spectral reflectance to intercepted photosynthetically active radiation (PAR) and net CO2 exchange in a natural canopy composed of the marsh cordgrass (Spartina alterniflora). Reflectance measurements made by a hand-held radiometer with Landsat TM spectral wavebands are used to compute remote sensing indices such as the normalized difference vegetation index. Consideration is given to the impact of standing dead canopy material on the relationship between intercepted PAR and spectral vegetation indices and the impact of changes in photosynthetic efficiency on the relationship between vegetation indices and CO2 exchange rates. The results suggest that quantitative remote assessment of photosynthesis and net gas exchange in natural vegetation is feasible, especially if the analysis incorporates information on biological responses to environmental variables.

Description: A method to detect plant water stress by remote sensing is proposed using indices of near-IR and mid-IR wavelengths. The ability of the Leaf Water Content Index (LWCI) to determine leaf relative water content (RWC) is tested on species with different leaf morphologies. The way in which the Misture Stress Index (MSI) varies with RWC is studied. On test with several species, it is found that LWCI is equal to RWC, although the reflectances at 1.6 microns for two different RWC must be known to accurately predict unknown RWC. A linear correlation is found between MSI and RWC with each species having a different regression equation. Also, MSI is correlated with log sub 10 Equivalent Water Thickness (EWT) with data for all species falling on the same regression line. It is found that the minimum significant change of RWC that could be detected by appying the linear regression equation of MSI to EWT is 52 percent. Because the natural RWC variation from water stress is about 20 percent for most species, it is concluded that the near-IR and mid-IR reflectances cannot be used to remotely sense water stress.

Description: A finite element formulation is presented for determining the large-amplitude free and steady-state forced vibration response of arbitrarily laminated anisotropic composite rectangular thin plates. The nonlinear stiffness and harmonic force matrices of an arbitrarily laminated composite rectangular plate element are developed for nonlinear free and forced vibration analyses. The linearized updated-mode method with nonlinear time function approximation is employed for the solution of the system nonlinear eigenvalue equations. The amplitude-frequency relations for convergence with gridwork refinement, different boundary conditions, aspect ratios, lamination angles and number of plies are presented. The finite element results are compared with available approximate continuum solutions.

Description: The objectives are to present the derivation of the new virtual crack closure technique, evaluate the accuracy of the technique, and finally to present the results of a limited parametric study of laminates with a postbuckled delamination. Although the new virtual crack closure technique is general, only homogeneous, isotropic laminates were analyzed. This was to eliminate the variation of flexural stiffness with orientation, which occurs even for quasi-isotropic laminates. This made it easier to identify the effect of geometrical parameters on G. The new virtual crack closure technique is derived. Then the specimen configurations are described. Next, the stress analyses is discussed. Finally, the virtual crack closure technique is evaluated and then used to calculate the distribution of G along the delamination front of several laminates with a postbuckled delamination.

Description: A description of the finite element implementation of Robinson's unified viscoplastic model into the General Purpose Finite Element Program (MARC) is presented. To demonstrate its application, the implementation is applied to some uniaxial and multiaxial problems. A comparison of the results for the multiaxial problem of a thick internally pressurized cylinder, obtained using the finite element implementation and an analytical solution, is also presented. The excellent agreement obtained confirms the correct finite element implementation of Robinson's model.

Description: In the dynamic formulation of holonomic and nonholonomic systems based on D'Alembert-Lagrange equation, the forces of constraints are maintained in the dynamic equations by introducing auxiliary variables, called Lagrange multipliers. This approach introduces a set of generalized reaction forces associated with the system generalized coordinates. Different sets of variables can be used as generalized coordinates and accordingly, the generalized reactions associated with these generalized coordinates may not be the actual reaction forces at the joints. In rigid body dynamics, the generalized reaction forces and the actual reaction forces at the joints represent equipollent systems of forces since they produce the same total forces and moments at and about any point on the rigid body. This is not, however, the case in deformable body analyses wherein the generalized reaction forces depend on the system generalized reference and elastic coordinates. In this paper, a method for determining the actual reaction forces at the joints from the generalized reaction forces in deformable multibody systems is presented.

Description: The Boundary Force Method (BFM), a form of indirect boundary element method, is used to analyze composite laminates with cracks. The BFM uses the orthotropic elasticity solution for a concentrated horizontal and vertical force and a moment applied at a point in a cracked, infinite sheet as the fundamental solution. The necessary stress functions for this fundamental solution were formulated using the complex variables theory of orthotropic elasticity. The current method is an improvement over a previous method using only forces and no moment. The improved method was verified by comparing it to accepted solutions for a finite-width, center-crack specimen subjected to uniaxial tension. Four graphite/epoxy laminates were used: (0 + or - 45/90)sub s, (0), (+ or - 45)sub s, and (+ or - 30)sub s. The BFM results agreed well with accepted solutions. Convergence studies showed that with the addition of the moment in the fundamental solution, the number of boundary elements required for a converged solution was significantly reduced. Parametric studies were done for two configurations for which no orthotropic solutions are currently available; a single edge crack and an inclined single edge crack.

Description: The Landsat Thematic Mapper (TM) imagery of the Basin and Range province of southern Nevada was analyzed to identify and map volcanic rock assemblages at three Tertiary calderas. It was found that the longer-wavelength visible and the NIR TM Bands 3, 5, and 7 provide more effective lithologic discrimination than the shorter-wavelength bands, due partly to deeper penetration of the longer-wavelength bands, resulting in more lithologically driven radiances. Shorter-wavelength TM Bands 1 and 2 are affected more by surficial weathering products including desert varnish which may or may not provide an indirect link to lithologic identity. Guidelines for lithologic analysis of volcanic terrains using Landsat TM imagery are outlined.

Description: Forest stand structure and biomass data were collected using conventional forest inventory techniques in tropical, subtropical, and warm temperate forest biomes. The feasibility of detecting tropical forest successional age class and total biomass differences using Landsat-Thematic mapper (TM) data, was evaluated. The Normalized Difference Vegetation Index (NDVI) calculated from Landsat-TM data were not significantly correlated with forest regeneration age classes in the mountain terrain of the Luquillo Experimental Forest, Puerto Rico. The low sun angle and shadows cast on steep north and west facing slopes reduced spectral reflectance values recorded by TM orbital altitude. The NDVI, calculated from low altitude aircraft scanner data, was significatly correlated with forest age classes. However, analysis of variance suggested that NDVI differences were not detectable for successional forests older than approximately 15-20 years. Also, biomass differences in young successional tropical forest were not detectable using the NDVI. The vegetation index does not appear to be a good predictor of stand structure variables (e.g., height, diameter of main stem) or total biomass in uneven age, mixed broadleaf forest. Good correlation between the vegetation index and low biomass in even age pine plantations were achieved for a warm temperate study site. The implications of the study for the use of NDVI for forest structure and biomass estimation are discussed.

Description: Thematic Mapper data were analyzed with respect to its capability for mapping the complex structure and dynamics of suspended matter distribution in the coastal area of the German Bight (North Sea). Three independent pieces of information were found by factor analysis of all seven TM channels: suspended matter concentration, atmospheric scattering, and sea surface temperature. For the required atmospheric correction, the signal-to-noise ratios of Channels 5 and 7 have to be improved by averaging over 25 x 25 pixels, which also makes it possible to monitor the aerosol optical depth and aerosol type over cloud-free water surfaces. Near-surface suspended matter concentrations may be detected with an accuracy of factor less than 2 by using an algorithm derived from radiative transfer model calculation. The patchiness of suspended matter and its relation to underwater topography was analyzed with autocorrelation and cross-correlation.

Description: Reflectances calculated from TM data and corrected for atmospheric effects correspond with in situ measured reflectances in the nadir-viewing mode, and are shown to be related to a glacier's mass balance if measured over a period of years. A reflectance of 0.895 for a test site in the Wrangell Mountains, Alaska, was calculated from TM Band 4 (0.76 - 0.90 micron) data and corrected for atmospheric effects. This value was comparable to the in situ reflectance of 0.90 measured in the same 0.76 - 0.90 micron wavelength region. For the same site, a reflectance value of 0.79 derived from integrating over most (0.40 - 3.0 micron) of the reflective portion of the electromagnetic spectrum was quite different from the integrated reflectance of 0.95 calculated for the spectral range 0.40 - 1.0 micron. This demonstrates the importance of using the full reflective energy spectrum for calculating the albedo of snow, and for obtaining a meaningful computation of a glacier's energy and mass balance change.

Description: This paper's studies were initiated under the NASA program for the purpose of conducting the earth sciences research using the Landsat Thematic Mapper. The goals of the program include studies of the factors influencing the growth, health, condition, and distribution of vegetation on the earth; the processes controlling the evolution of the earth's crust; the earth's water budget and the hydrologic processes that operate at local, regional, and global scales; the physical and chemical interaction between different types of surficial materials; and the interaction between the earth's surface and its atmosphere. Twenty-seven domestic and five foreign investigations were initiated in 1985, with the results from most of them already published (one study was terminated due to the delay in the TDRSS). Twelve of the studies addressed hydrology, snow and ice, coastal processes, and near-shore oceanographic phenomena; seven addressed vegetation, soils, or animal habitat; and twelve addressed geologic subjects.

Description: The fundamental relationship between the morphology of a composite laminate and the resulting free edge effects is explored and related to delamination failures. Cross-ply, angle-ply, and quasi-isotropic laminates are discussed in detail. It is shown that the local mismatch in elastic properties of adjacent layers and the global stacking sequence of a laminate both have a significant influence on the interlaminar stresses and delamination failures.

Description: A method is developed for the optimal design of composite links based on dynamic performance criteria directly related to structural modal damping and dynamic stiffness. An integrated mechanics theory correlates structural composite damping to the parameters of basic composite material systems, laminate parameters, link shape, and modal deformations. The inclusion of modal properties allows the selective minimization of vibrations associated with specific modes. Ply angles and fiber volumes are tailored to obtain optimal combinations of damping and stiffness. Applications to simple composite links indicate wide margins for trade-offs and illustrate the importance of various design variables to the optimal design.

Description: Comparison between the microwave polarized difference temperature (MPDT) derived from 37 GHz band data and the normalized difference vegetation index (NDVI) derived from near-infrared and red bands, from several empirical investigations are summarized. These indicate the complementary character of the two measures in environmental monitoring. Overall the NDVI is more sensitive to green leaf activity, whereas the MPDT appears also to be related to other elements of the above-ground biomass. Monitoring of hydrological phenomena is carried out much more effectively by the MPDT. Further work is needed to explain spectral and temporal variation in MPDT both through modelling and field experiments.

Description: Coincident Scanning Microwave Multi-channel Radiometer 37 GHz and Advanced Very High Resolution Radiometer normalized difference vegetation index satellite data have been compared from drought-affected regions of sub-Saharan Africa and northeastern Brazil for the time period of 1980-1985. Although the two satellite data types can be highly correlated, differences between them were found for the Sahel zone in 1985 and for northeastern Brazil from 1984-1985. These findings suggest that scattering or surface roughness contributions may be greater than previously assumed for the 37 GHz microwave data.

Description: The microwave polarization difference temperature (MPDT) prepared from 37 GHz Scanning Multifrequnecy Microwave Radiometer (SMMR) data from the Nimbus-7 satellite allows a unique vision of hydrological features, especially of humid areas. From the series of monthly images from 1979 to 1985 the major hydrological features of four major river basins of South America are briefly examined. These include rivers of the Amazon, La Plata, Orinoco and Sao Francisco, wetlands in the La Plata and Amazon basins and floods in the La Plata basin. Many hydrological features appear much as they do on conventional maps. An interesting prespective can be derived from the ways in which they differ from conventional maps.

Description: Monthly 37 GHz microwave polarization difference temperatures (MPDT) derived from the Nimbus-7 scanning multichannel microwave radiometer (SMMR) for southern Africa from 1979 to 1985 are compared with rainfall and Advanced Very High Resolution Radiometer (AVHRR) normalized difference vegetation index (NDVI) data. MPDT rose sharply during a drought episode which occurred within the period included in the data. The rise was seen not only in the growing season, but also in the dry season MPDT when no actively photosynthetic, water-containing leaves are present. The results suggest that scattering of the emitted microwave radiation by dead and living vegetation is a more important factor than has previously been recognized.

Description: Simplified radiative transfer models for 37 GHz brightness temperatures are developed and applied to observations taken with the Nimbus-7 scanning multichannel microwave radiometer. It is found that, as the vegetation density increases, polarization difference decreases. The polarization difference also decreases with increasing surface roughness. Empirical relations are established between the satellite observations and annual rainfall over Africa and Australia, primary productivity, and actual evaporation. Color-coded maps of primary productivity and evaporation are presented.

Description: Local buckling can cause large interlaminar stresses along the delamination front, which can lead to delamination growth. This paper examines several methods of calculating strain-energy release rates, which are often used to predict delamination growth. The thin-film plate analysis, which was least expensive, calculated the total strain-energy release rate G(T) quite accurately. However, the stress field along the delamination front is highly mixed-mode and has no fixed ratio of G(I) to G(II). Since plate analysis can only calculate G(T), it would not be useful for accurate predictions of delamination growth if mode mix is important.

Description: Over the years, NASA has been conducting the Advanced Earth-to-Orbit (AETO) Propulsion Technology Program to provide the knowledge, understanding, and design methodology that will allow the development of advanced Earth-to-orbit propulsion systems with high performance, extended service life, automated operations, and diagnostics for in-flight health monitoring. The objective of the Aerothermodynamic Loads Definition Study is to develop methods to more accurately predict the operating environment in AETO propulsion systems, such as the Space Shuttle Main Engine (SSME) powerhead. The approach taken consists of 2 parts: to modify, apply, and disseminate existing computational fluid dynamics tools in response to current needs and to develop new technology that will enable more accurate computation of the time averaged and unsteady aerothermodynamic loads in the SSME powerhead. The software tools are detailed. Significant progress was made in the area of turbomachinery, where there is an overlap between the AETO efforts and research in the aeronautical gas turbine field.

Description: The purpose of the Probabilistic Structural Analysis Method (PSAM) project is to develop structural analysis capabilities for the design analysis of advanced space propulsion system hardware. The boundary element method (BEM) is used as the basis of the Probabilistic Advanced Analysis Methods (PADAM) which is discussed. The probabilistic BEM code (PBEM) is used to obtain the structural response and sensitivity results to a set of random variables. As such, PBEM performs analogous to other structural analysis codes such as finite elements in the PSAM system. For linear problems, unlike the finite element method (FEM), the BEM governing equations are written at the boundary of the body only, thus, the method eliminates the need to model the volume of the body. However, for general body force problems, a direct condensation of the governing equations to the boundary of the body is not possible and therefore volume modeling is generally required.

Description: A major research and technology program in Probabilistic Structural Analysis Methods (PSAM) is currently being sponsored by the NASA Lewis Research Center with Southwest Research Institute as the prime contractor. This program is motivated by the need to accurately predict structural response in an environment where the loadings, the material properties, and even the structure may be considered random. The heart of PSAM is a software package which combines advanced structural analysis codes with a fast probability integration (FPI) algorithm for the efficient calculation of stochastic structural response. The basic idea of PAAM is simple: make an approximate calculation of system response, including calculation of the associated probabilities, with minimal computation time and cost, based on a simplified representation of the geometry, loads, and material. The deterministic solution resulting should give a reasonable and realistic description of performance-limiting system responses, although some error will be inevitable. If the simple model has correctly captured the basic mechanics of the system, however, including the proper functional dependence of stress, frequency, etc. on design parameters, then the response sensitivities calculated may be of significantly higher accuracy.

Description: Two-dimensional optical strain measurements on high temperature test specimens are presented. This two-dimensional capability is implemented through a rotatable sensitive strain axis. Three components of surface strain can be measured automatically, from which the first and second principal strains are calculated. One- and two-dimensional strain measurements at temperatures beyond 750 C with a resolution of 15 microstrain are demonstrated. The system is based on a one-dimensional speckle shift technique. The speckle shift technique makes use of the linear relationship between surface strain and the differential shift of laser speckle patterns in the diffraction plane. Laser speckle is a phase effect that occurs when spatially coherent light interacts with an optically rough surface. Since speckle is generated by any diffusely reflecting surface, no specimen preparation is needed to obtain a good signal. Testing was done at room temperature on a flat specimen of Inconel 600 mounted in a fatigue testing machine. A load cell measured the stress on the specimen before and after acquiring the speckle data. Strain components were measured at 0 C (parallel to the load axis) and at plus or minus 45 C, and plots indicate the calculated values of the first and second principal strains. The measured values of Young's modulus and Poisson's ratio are in good agreement with handbook values. Good linearity of the principal strain moduli at high temperatures indicate precision and stability of the system. However, a systematic error in the high-temperature test setup introduced a scale factor in the slopes of the two-dimensional stress-strain curves. No high temperature effects, however, have been observed to degrade speckle correlation.

Description: The fusion of the probabilistic finite element method (PFEM) and reliability analysis for probabilistic fracture mechanics (PFM) is presented. A comprehensive method for determining the probability of fatigue failure for curved crack growth was developed. The criterion for failure or performance function is stated as: the fatigue life of a component must exceed the service life of the component; otherwise failure will occur. An enriched element that has the near-crack-tip singular strain field embedded in the element is used to formulate the equilibrium equation and solve for the stress intensity factors at the crack-tip. Performance and accuracy of the method is demonstrated on a classical mode 1 fatigue problem.

Description: The probabilistic structural analysis method (PSAM) was developed to analyze the effects of fluctuating loads, variable material properties, and uncertain analytical models especially for high performance structures such as the Space Shuttle Main Engine turbopump blades. Risk is calculated after expensive service experience. However, probabilistic structural analysis provides a rational alternative method to quantify uncertainties in the structural performance and durability. NESSUS (Numerical Evaluation of Stochastic Structures Under Stress) was developed as a probabilistic structural analysis computer code which integrates finite element methods and reliability algorithms, capable to predicting the probability distributions of structural response variables such as stress, displacement, natural frequencies, and buckling loads. This computer code is detailed.

Description: The Probabilistic Structural Analysis Methods (PSAM) project developed at the Southwest Research Institute integrates state-of-the-art structural analysis techniques with probability theory for the design and analysis of complex large-scale engineering structures. An advanced efficient software system (NESSUS) capable of performing complex probabilistic analysis has been developed. NESSUS contains a number of software components to perform probabilistic analysis of structures. These components include: an expert system, a probabilistic finite element code, a probabilistic boundary element code and a fast probability integrator. The NESSUS software system is shown. An expert system is included to capture and utilize PSAM knowledge and experience. NESSUS/EXPERT is an interactive menu-driven expert system that provides information to assist in the use of the probabilistic finite element code NESSUS/FEM and the fast probability integrator (FPI). The expert system menu structure is summarized. The NESSUS system contains a state-of-the-art nonlinear probabilistic finite element code, NESSUS/FEM, to determine the structural response and sensitivities. A broad range of analysis capabilities and an extensive element library is present.

Description: Constitutive and life prediction models are developed and verified for materials typically used in hot gas path components of reusable space propulsion systems over the range of relevant operating environments. The efforts were centered on the development of crack initiation life prediction methods, while the efforts of a counterpart group were centered on the development of cyclic crack propagation life prediction methods. The complexion of the active tasks are presented. A significant new task started this year will incorporate the various material constitutive and life prediction models developed in this program into a comprehensive creep-fatigue damage analysis and life assessment computer code. The program will function as a postprocessor to general structural analysis programs (such as finite element or boundary element codes) using the output of such analyses (stress, strain, and temperature fields as functions of time) as the input to the damage analysis and life assessment code. The code will be designed to execute on engineering/scientific workstations and will feature a windowing, mouse-driven user interface. Current plans call for the code to be finished and made available for use in mid 1991.

Description: In structural analysis the amount of computational time necessary for a solution is proportional to the number of degrees of freedom times the bandwidth squared. In implicit time analysis, this must be done at each discrete point in time. If, in addition, the problem is nonlinear, then this solution must be iterated at each point in time. If the bandwidth is large, the size of the problem that can be analyzed is severely limited. The multi-grid method is a possible algorithm which can make this solution much more computationally efficient. This method has been used for years in computational fluid mechanics. It works on the fact that relaxation is very efficient on the high frequency components of the solution (nearest neighbor interactions) and not very good on low frequency components of the solution (far interactions). The multi-grid method is then to relax the solution on a particular model until the residual stops changing. This indicates that the solution contains the higher frequency components. A coarse model is then generated for the lower frequency components to the solution. The model is then relaxed for the lower frequency components of the solution. These lower frequency components are then interpolated to the fine model. In computational fluid mechanics the equations are usually expressed as finite differences.

Description: The deployment, assembly and mission-oriented maneuvering of space structures in orbit will trigger large motions of flexible, truss-type structures. In addition, the presence of on-board controls both for attitude stabilization and specified vibration tolerance requirements may further complicate the dynamic behavior of the orbiting structures. Because of safety and cost considerations, the dynamic response of the combined structural and control systems must be predicted reliably. This need can only be met through the development of reliable and efficient simulation capabilities, since there is general agreement that on-orbit experiments should be limited because of cost, time and facility constraints. The long-term objective of this research effort is to develop a next-generation computer simulator for the dynamics and control of large space structures. The simulator will be based on integrating four research thrusts: a new multibody dynamics formulation methodology, modeling capabilities in long/slender truss-beam components with realistic joints, efficient computational procedures that can be implemented either in sequential or concurrent computers, and prototype simulation modules that can be easily processed into a modern large-scale engineering software system such as the NASA/Computational Structural Mechanics (CSM) testbed.

Description: The 3-D inelastic analysis method is a focused program with the objective to develop computationally effective analysis methods and attendant computer codes for three-dimensional, nonlinear time and temperature dependent problems present in the hot section of turbojet engine structures. Development of these methods was a major part of the Hot Section Technology (HOST) program over the past five years at Lewis Research Center.

Description: A summary of the status of this five-year project which is now in its third year of research and development is presented. The goal of the project is the development of several methodologies for probabilistic structural modeling. Probabilistic structural modeling consists of stochastic models of material properties, part geometries, boundary conditions, as well as loading conditions. The current presentation focuses on one methodology - coupling of an advanced finite element structural analysis code with probabilistic modeling strategies. The essential algorithm developments for combining the finite element and probabilistic analysis methods are reported. The validity of the resulting probabilistic structural analysis method is confirmed through a series of test problems with exact results based on Monte Carlo simulations. Additionally, the applicability of the method to a Space Propulsion System (a turbine blade) is demonstrated for static stresses.

Description: This project supports the conversion of codes in Computational Structural Mechanics (CSM) to a parallel form which will efficiently exploit the computational power available from multiprocessors. The work is a part of a comprehensive, FORTRAN-based system to form a basis for a parallel version of the NICE/SPAR combination which will form the CSM Testbed. The software is macro-based and rests on the force methodology developed by the principal investigator in connection with an early scientific multiprocessor. Machine independence is an important characteristic of the system so that retargeting it to the Flex/32, or any other multiprocessor on which NICE/SPAR might be imnplemented, is well supported. The principal investigator has experience in producing parallel software for both full and sparse systems of linear equations using the force macros. Other researchers have used the Force in finite element programs. It has been possible to rapidly develop software which performs at maximum efficiency on a multiprocessor. The inherent machine independence of the system also means that the parallelization will not be limited to a specific multiprocessor.

Description: A new numerical algorithm for the solution of large-order eigenproblems typically encountered in linear elastic finite element systems is presented. The architecture of parallel processing is utilized in the algorithm to achieve increased speed and efficiency of calculations. The algorithm is based on the frontal technique for the solution of linear simultaneous equations and the modified subspace eigenanalysis method for the solution of the eigenproblem. Assembly, elimination and back-substitution of degrees of freedom are performed concurrently, using a number of fronts. All fronts converge to and diverge from a predefined global front during elimination and back-substitution, respectively. In the meantime, reduction of the stiffness and mass matrices required by the modified subspace method can be completed during the convergence/divergence cycle and an estimate of the required eigenpairs obtained. Successive cycles of convergence and divergence are repeated until the desired accuracy of calculations is achieved. The advantages of this new algorithm in parallel computer architecture are discussed.

Description: The overall objective of this research is to develop efficient methods for the solution of linear and nonlinear systems of equations on parallel and supercomputers, and to apply these methods to the solution of problems in structural analysis. Attention has been given so far only to linear equations. The methods considered for the solution of the stiffness equation Kx=f have been Choleski factorization and the conjugate gradient iteration with SSOR and Incomplete Choleski preconditioning. More detail on these methods will be given on subsequent slides. These methods have been used to solve for the static displacements for the mast and panel focus problems in conjunction with the CSM testbed system based on NICE/SPAR.

Description: The background for the CSM Testbed Architecture and a description of the three architecture related tasks under the NASA/Langley CSM contract (NAS1-18444) are given. For each task the objectives, subtasks chosen to achieve the objectives, and the accomplishments are presented. For Task 2, Near-Term Enhancements to the CSM Testbed Architecture, the primary objectives are to quickly correct: (1) inefficiency in the GAL data manager, and (2) deficiency in the CLIP-interpreted command language, CLAMP. The corresponding modifications should preserve upward-compatibility to a reasonable extent, while at the same time increasing the flexibility and extensibility of the CSM Testbed. Researchers have increased the efficiency of GAL by a factor of 2+ and have made several improvements in CLIP. For Task 5, Matrix Algebra Methods and Utilities, the goal is to investigate the current capabilities of the CSM Testbed and the required improvements to the CSM Testbed for performing the matrix algebraic operations required for the analysis of present and future CSM focal problems using advanced methods. Researchers have made an extensive study of the matrix algebra functions in SPAR that includes documenting these routines and supplying inline comments.

Description: The PISCES project started in 1984 under the sponsorship of the NASA Computational Structural Mechanics (CSM) program. A PISCES 1 programming environment and parallel FORTRAN were implemented in 1984 for the DEC VAX (using UNIX processes to simulate parallel processes). This system was used for experimentation with parallel programs for scientific applications and AI (dynamic scene analysis) applications. PISCES 1 was ported to a network of Apollo workstations by N. Fitzgerald.

Description: In January 1987, SPARTA received a Phase I SBIR award from the NASA Lewis Research Center to investigate the feasibility of a finite element solver implemented on multiple VLSI processors. The transputer was chosen as the processor for the feasibility study since it combined low cost with high performance and was specifically designed to directly link with other transputers to form networks of multiple processors. A brief description of transputers, a summary of the SBIR feasibility study, and a discussion of issues concerning a large scale transputer based finite element solver (TBFES) and the performance levels which can be expected are discussed.

Description: Damage assessment of structural assemblies is treated as an identification problem. A brief review of identification methods is first presented with particular focus on the output error approach. The use of numerical optimization methods in identifying the location and extent of damage in structures is studied. The influence of damage on eigenmode shapes and static displacements is explored as a means of formulating a measure of damage in the structure. Preliminary results obtained in this study are presented and special attention is directed at the shortcomings associated with the nonlinear programming approach to solving the optimization problem.

Description: Within the past several years, the Jet Propulsion Laboratory has designed and built major ground antenna structures in Spain, Australia, and California. One of the antennas at each location is a 70 meter-diameter structure that is a retrofit of the existing 64 meter antenna. The 64 meter existing antennas were first stripped back to a 34 meter interior and then completely new construction with deeper trusses was added to extend the interior to 70 meters. The 70 meter project included the rare opportunity to collect field data to compare with predictions of the finite-element analytical models. The new quadripod design was tested for its lower mode natural frequencies and the main reflector was measured by theodolite to determine deflections of subsets of the backup-structure deformations under load. The emphasis here is to examine measurement results and possibly provide some appreciation of the relationship of predictions made from the design model to actual measurements.