ALBERT

Description: The development of an integrated approach to the modeling of forest dynamics encompassing submodels of forest growth and succession, soil processes and radiation interactions, is reported. Remote sensing technology is a key element of this study in that it provides data for developing, initializing, updating, and validating the models. The objectives are reviewed, the data collected and models in use are discussed, and a framework for studying interactions between the forest growth, soil process and energy interaction components, is described. Remote sensing technology used in the study includes optical and microwave field, aircraft and satellite borne instruments. The types of data collected during intensive field and aircraft campaigns included bidirectional reflectance, thermal emittance and multifrequency, multipolarization synthetic aperture radar backscatter. Synthetic imagery of derived products such as forest biomass and NDVI (Normalized Difference Vegetative Index), and collections of ground data are being assembled in a georeferenced data base. These data are used to drive or test multidiscipline simulations of forested ecosystems. Enhancements to the modeling environment permit considerable flexibility in configuring simulations and selecting results for reporting and graphical display.

Description: The quantitative interpretation of satellite observations requires the use of mathematical tools to extract the desired information on terrestrial environments from the radiation data collected in space. A whole range of approaches can be pursued, from the development of models capable of explaining the nature of the physical signal being measured and of characterizing the state of the system under observation, to the empirical correlations between the variables of interest and the space measurements. The premises and implications of these approaches are outlined, paying special attention to the mathematical and numerical requirements. The role and specific applications of empirical bidirectional reflectance models is also discussed, even though these models do not contribute to the understanding of the theory of radiation transfer or to the assessment of the variables of interest. The advantages and drawbacks of these various approaches and the research priorities for the next few years are discussed in the context of the planned availability of new sensors.

Description: Management of crop residues, the portion of a crop left in the field after harvest, is an important conservation practice for minimizing soil erosion and for improving water quality. Quantification of crop residue cover is required to evaluate the effectiveness of conservation tillage practices. Methods are needed to quantify residue cover that are rapid, accurate, and objective. The fluorescence of crop residue was found to be a broadband phenomenon with emission maxima at 420 to 495 nm for excitations of 350 to 420 nm. Soils had low intensity broadband emissions over the 400 to 690 nm region for excitations of 300 to 600 nm. The range of relative fluorescence intensities for the crop residues was much greater than the fluorescence observed of the soils. As the crop residues decompose their blue fluorescence values approach the fluorescence of the soil. Fluorescence techniques are concluded to be less ambiguous and better suited for discriminating crop residues and soils than reflectance methods. If properly implemented, fluorescence techniques can be used to quantify, not only crop residue cover, but also photosynthetic efficiency in the field.

Description: The importance of the measurement of wind fields is discussed. Wind regime data can be used to infer the amount and type of wind induced (aerolian) transport of sand and dust, or to establish global circulation models, for example on other planets. Since local measurements are costly and often impossible, it is desired to infer such data from remotely sensed information. A potential mechanism for remotely inferring the wind regime by using synthetic aperture radar data to describe the roughness of the surface is described. A project to estimate the practicality of using such a mechanism is described. An experiment that extends the mechanism to vegetated sites, where the goal is to measure potential for erosion, is reported.

Description: Surface reflectance is required to quantitatively investigate molecular absorption and particle scattering properties of materials on the Earth's surface. Atmospheric aerosol optical depth, surface pressure and water vapor are required to constrain a radiative transfer code for the inversion of measured spectral radiance to apparent surface reflectance. A suite of algorithms using nonlinear least squares fitting techniques are described that directly estimate these atmospheric parameters from spectral radiance measured by the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). The derived atmospheric parameters are used to constrain a radiative transfer code for the inversion of the imaging spectrometer radiance to apparent reflectance. The derived apparent reflectance is validated with respect to in situ measurement on the same target.

Description: Investigations designed to study land surface hydrologic-atmospheric interactions, showing the potential of L band passive microwave radiometry for measuring surface soil moisture over large areas, are discussed. Satisfying the data needs of these investigations requires the ability to map large areas rapidly. With aircraft systems this means a need for more beam positions over a wider swath on each flightline. For satellite systems the essential problem is resolution. Both of these needs are currently being addressed through the development and verification of Electronically Scanned Thinned Array Radiometer (ESTAR) technology. The ESTAR L band radiometer was evaluated for soil moisture mapping applications in two studies. The first was conducted over the semiarid rangeland Walnut Gulch watershed located in south eastern Arizona (U.S.). The second was performed in the subhumid Little Washita watershed in south west Oklahoma (U.S.). Both tests showed that the ESTAR is capable of providing soil moisture with the same level of accuracy as existing systems.

Description: A weather resistant automatic scanning Sun photometer system is assessed and demonstrated as practical for measurements of aerosol concentrations and properties at remote sites. Interfaced with a transmitter using the Geostationary Data Collection System (GDCS), the data are processed in near real time. The processing allows a time dependence of the aerosols and water vapor and an ongoing assessment of the health and calibration of the instruments. The system's automatic data acquisition, transmission, and processing offer immediate application to atmospheric monitoring and modeling on a regional to global scale and validation of satellite retrievals. It is estimated that under normal circumstances the retrieved aerosol optical thickness has a network wide accuracy of +/- 0.02 from 340 nm to 1020 nm, water vapor +/- 0.2 cm and size distribution from 0.1 to 3 micrometers.

Description: As part of a global program to validate the ocean surface sensors on board ERS-1, a joint experiment on the Grand Banks of Newfoundland was carried out in Nov. 1991. The principal objective was to provide a field validation of ERS-1 Synthetic Aperture Radar (SAR) measurement of ocean surface structure. The NASA-P3 aircraft measurements made during this experiment provide independent measurements of the ocean surface along the validation swath. The Radar Ocean Wave Spectrometer (ROWS) is a radar sensor designed to measure direction of the long wave components using spectral analysis of the tilt induced radar backscatter modulation. This technique greatly differs from SAR and thus, provides a unique set of measurements for use in evaluating SAR performance. Also, an altimeter channel in the ROWS gives simultaneous information on the surface wave height and radar mean square slope parameter. The sets of geophysical parameters (wind speed, significant wave height, directional spectrum) are used to study the SAR's ability to accurately measure ocean gravity waves. The known distortion imposed on the true directional spectrum by the SAR imaging mechanism is discussed in light of the direct comparisons between ERS-1 SAR, airborne Canadian Center for Remote Sensing (CCRS) SAR, and ROWS spectra and the use of the nonlinear ocean SAR transform.

Description: Terrain slopes, which can be measured with Synthetic Aperture Radar (SAR) interferometry either from a height map or from the interferometric phase gradient, were used to calculate the local incidence angle and the correct pixel area. Both are required for correct thematic interpretation of SAR data. The interferometric correlation depends on the pixel area projected on a plane perpendicular to the look vector and requires correction for slope effects. Methods for normalization of the backscatter and interferometric correlation for ERS-1 SAR are presented.

Description: The relationship between the gradient of the interferometric phase and the terrain slope, which, it is thought, would allow a derivation of the terrain slopes without phase unwrapping, is presented. A linear relationship between the interferometric phase gradient and the terrain slopes was found. A quantitative error analysis showed that only very small errors are introduced by these approximations for orbital Synthetic Aperture Radar (SAR) geometries. An example of a slope map for repeat pass interferometry from ERS-1 SAR data is given. A number of direct and indirect applications of the terrain slope are indicated: erosion and avalanche hazard studies, radiometric calibration of SAR data, and normalization of the interferometric correlation coefficient.

Description: Research on the use of active microwaves in remote sensing, presented during plenary and poster sessions, is summarized. The main highlights are: calibration techniques are well understood; innovative modeling approaches have been developed which increase active microwave applications (segmentation prior to model inversion, use of ERS-1 scatterometer, simulations); polarization angle and frequency diversity improves characterization of ice sheets, vegetation, and determination of soil moisture (X band sensor study); SAR (Synthetic Aperture Radar) interferometry potential is emerging; use of multiple sensors/extended spectral signatures is important (increase emphasis).

Description: Crop residues, the portion of the crop left in the field after harvest, can be an important management factor in controlling soil erosion. Methods to quantify residue cover are needed that are rapid, accurate, and objective. Scenes with known amounts of crop residue were illuminated with long wave ultraviolet (UV) radiation and fluorescence images were recorded with an intensified video camera fitted with a 453 to 488 nm band pass filter. A light colored soil and a dark colored soil were used as background for the weathered soybean stems. Residue cover was determined by counting the proportion of the pixels in the image with fluorescence values greater than a threshold. Soil pixels had the lowest gray levels in the images. The values of the soybean residue pixels spanned nearly the full range of the 8-bit video data. Classification accuracies typically were within 3(absolute units) of measured cover values. Video imaging can provide an intuitive understanding of the fraction of the soil covered by residue.

Description: Global study of land surface properties uses AVHRR channels 1 and 2, but channel 3 may be of interest, although its use requires preprocessing. It consists of both a reflective part and an emissive part, the former can be derived from T3, T4 and T5. Since the water vapor affects channel 3, its content is retrieved from the channel 4 and 5 using the split window technique. A formula of reflective part retrieval at 3.75 micrometers is tested in the case of sunglint observations where the emissivities of channels 4 and 5 can be set to the unity. The formula is adapted and validated to land surface using the FIFE-87 data set. Preliminary applications of the reflectance at 3.75 micrometers to the studies of surface properties retrieval, aerosol retrieval over land, and desertic aerosol retrieval, are addressed.

Description: The polarization of the sunlight scattered by atmospheric aerosols or cloud droplets and reflected from ground surfaces or plant canopies may convey much information when used for remote sensing purposes. The typical polarization features of aerosols, cloud droplets, and plant canopies, as observed by ground based and airborne sensors, are investigated, looking especially for those invariant properties amenable to description by simple models when possible. The question of polarization measurements from space is addressed. The interest of such measurements for remote sensing purposes is investigated, and their feasibility is tested by using results obtained during field campaigns of the airborne POLDER instrument, a radiometer designed to measure the directionality and polarization of the sunlight scattered by the ground atmosphere system.

Description: Papers focused on land surface, atmospheric, and ocean properties are reported. Specific comments pertaining to polarization, models and inversion, and measurements, are given. Recommendations are: continued research into the application potential of the BRDF (Bidirectional Reflectance Distribution Function) and polarization properties of ground surface and atmospheric targets; three dimensional models, which account for the statistical behavior of remotely sensed data, should be extended and inverted in order to support analysis of data potentially covering rolling terrain such that pixels represent heterogeneous mixtures of surface cover types and project ground footprints with sizes between 10 to 6 km, the ground pixel sizes of planned future sensors; available reflectance models should be further validated by means of multi dimensional (directional, spectral, temporal) field data and existing models should be intercompared in more depth to evaluate their performance and limitations; existing methods for model inversion should be validated in more depth in order to quantify the practical limitations and the expected accuracy of the parameters retrieved and new approaches should be developed based upon apriori knowledge of plant canopy development and spectral BRDF properties; there is a need to establish a protocol of validation and intercomparison of the indices and compositing techniques which have been proposed during these last years.

Description: Traditionally, the remote sensing community has relied totally on spectral knowledge to extract vegetation characteristics. However, there are other knowledge bases (KB's) that can be used to significantly improve the accuracy and robustness of inference techniques. Using AI (artificial intelligence) techniques a KB system (VEG) was developed that integrates input spectral measurements with diverse KB's. These KB's consist of data sets of directional reflectance measurements, knowledge from literature, and knowledge from experts which are combined into an intelligent and efficient system for making vegetation inferences. VEG accepts spectral data of an unknown target as input, determines the best techniques for inferring the desired vegetation characteristic(s), applies the techniques to the target data, and provides a rigorous estimate of the accuracy of the inference. VEG was developed to: infer spectral hemispherical reflectance from any combination of nadir and/or off-nadir view angles; infer percent ground cover from any combination of nadir and/or off-nadir view angles; infer unknown view angle(s) from known view angle(s) (known as view angle extension); and discriminate between user defined vegetation classes using spectral and directional reflectance relationships developed from an automated learning algorithm. The errors for these techniques were generally very good ranging between 2 to 15% (proportional root mean square). The system is designed to aid scientists in developing, testing, and applying new inference techniques using directional reflectance data.

Description: Most earth surfaces, particularly those supporting natural vegetation ecosystems, constitute structurally and spectrally complex surfaces that are distinctly non-Lambertian reflectors. Obtaining meaningful measurements of the directional radiances of landscapes and obtaining estimates of the complete bidirectional reflectance distribution functions of ground targets with complex and variable landscape and radiometric features are challenging tasks. Reasons for the increased interest in directional radiance measurements are presented, and the issues that must be addressed when trying to acquire directional radiances for vegetated land surfaces from different types of remote sensing platforms are discussed. Priority research emphases are suggested, concerning field measurements of directional surface radiances and reflectances for future research. Primarily, emphasis must be given to the acquisition of more complete and directly associated radiometric and biometric parameter data sets that will empower the exploitation of the 'angular dimension' in remote sensing of vegetation through enabling the further development and rigorous validation of state of the art plant canopy models.

Description: Synthetic aperture radar (SAR) images of the Greenland ice sheet collected by an airborne system clearly reveal the four melting facies of this sheet defined 30 years ago from snow stratigraphy studies by glaciologists. In particular, the radar echoes from the percolation facies have radiometric and polarimetric characteristics that are unique among terrestrial surfaces, but that resemble the exotic radar echoes recorded from the icy Galilean satellites. There, the radar signals interact with subsurface, massive ice features created in the cold, dry snow by seasonal melting and refreezing events. The subsurface features act as efficient reflectors of the incident radiation most likely via internal reflections. In the soaked-snow facies, the radar reflectivity is much lower because radar signals are attenuated by the wetter snow before they can interact with subsurface structures. Inversion algorithms to derive geophysical information from the SAR data are developed in both cases to estimate snow wetness in the soaked-snow facies and the mass of ice water retained in the percolation facies.

Description: An operational stratospheric correction scheme used after the Mount Pinatubo (Phillipines) eruption (Jun. 1991) is presented. The stratospheric aerosol distribution is assumed to be only variable with latitude. Each 9 days the latitudinal distribution of the optical thickness is computed by inverting radiances observed in the NOAA AVHRR channel 1 (0.63 micrometers) and channel 2 (0.83 micrometers) over the Pacific Ocean. This radiance data set is used to check the validity of model used for inversion by checking consistency of the optical thickness deduced from each channel as well as optical thickness deduced from different scattering angles. Using the optical thickness profile previously computed and radiative transfer code assuming Lambertian boundary condition, each pixel of channel 1 and 2 are corrected prior to computation of NDVI (Normalized Difference Vegetation Index). Comparison between corrected, non corrected, and years prior to Pinatubo eruption (1989 to 1990) NDVI composite, shows the necessity and the accuracy of the operational correction scheme.

Description: Aspects of aerosol studies and remote sensing are reviewed. Aerosol scatters solar radiation before it reaches the surface and scatters and absorbs it again after it is reflected from the surface and before it reaches the satellite sensor. The effect is spectrally and spatially dependent. Therefore atmospheric aerosol (dust, smoke and air pollution particles) has a significant effect on remote sensing. Correction for the aerosol effect was never achieved on an operational basis though several case studies were demonstrated. Correction can be done in a direct way by deriving the aerosol loading from the image itself and correcting for it using the appropriate radiative transfer model or by an indirect way, by defining remote sensing functions that are less dependent on the aerosol loading. To some degree this was already achieved in global remote sensing of vegetation where a composite of several days of NDVI (Normalized Difference Vegetation Index) measurements, choosing the maximal value, was used instead of a single cloud screened value. The Atmospheric Resistant Vegetation Index (ARVI) introduced recently for the NASA Earth Observing System EOS-MODIS is the most appropriate example of indirect correction, where the index is defined in such a way that the atmospheric effect in the blue spectral channel cancels to a large degree the atmospheric in the red channel in computations of a vegetation index. Atmospheric corrections can also use aerosol climatology and ground based instrumentation.

Description: A session dedicated to high spectral resolution in the solar spectrum, covering topics of calibration, atmospheric correction, geology/pedology, inland water, and vegetation, is reported. The session showed a high degree of diversity in the topics and the approaches used. It was highlighted that high spectral resolution data could provide atmospherically corrected ground level calibrated reflectance values. Important advances were shown in the use of radiative transfer models applied either on water bodies or vegetation. Several studies highlighted the high degree of redundancy contained in high spectral resolution data.

Description: Spectral absorption-coefficients (cross-sections) kappa(sub nu) (/cm/atm) have been measured in the 7.62, 8.97, and 12.3 micrometer bands of HCFC-22 (CHClF2) and the 10.6 micrometer bands of SF6 employing a high-resolution Fourier-transform spectrometer. Temperature and total pressure have been varied to simulate conditions corresponding to tropospheric and stratospheric layers in the atmosphere. The kappa(sub nu) are compared with values measured by us previously using a tunable diode laser spectrometer and with the appropriate entries in HITRAN and GEISA, two of the databases known to the atmospheric scientist. The measured absolute intensities of the bands are compared with previously published values.

Description: The analysis of balloon envelopes by the finite element (FE) method is plagued by convergence problems. A pratical FE analysis approach is based on the fact that in thin shells with non-zero Gaussian curvature the membrane solution component is essentially decoupled from the bending solution component. A proxy-problem is solved by using a small artificial bending stiffness that assures convergence without significantly affecting the membrane solution component. This approach has been previously validated on slightly overpressurized balloon envelopes. Extensions of this approach to more difficult problems in the structural analysis of balloon envelopes are presented. The convergence forcing modelling measures are discussed. Implications of the findings of the analysis results to future balloon designs are also discussed.

Description: Microwave radar and radiometer measurements of grasslands indicate a substantial reduction in sensor sensitivity to soil moisture in the presence of a thatch layer. When this layer is wet it masks changes in the underlying soil, making the canopy appear warm in the case of passive sensors (radiometer) and decreasing backscatter in the active case (scatterometer). A model for a grass canopy with thatch will be presented in this paper to explain this behavior and to compare with observations. The canopy model consists of three layers: grass, thatch, and the underlying soil. The grass blades are modeled by elongated elliptical discs and the thatch is modeled as a collection of disk shaped water droplets (i.e., the dry matter is neglected). The ground is homogeneous and flat. The distorted Born approximation is used to compute the radar cross section of this three layer canopy and the emissivity is computed from the radar cross section using the Peake formulation for the passive problem. Results are computed at L-band (1.4 GHz) and C-band (4.75 GHz) using canopy parameters (i.e., plant geometry, soil moisture, plant moisture, etc.) representative of Konza Prairie grasslands. The results are compared to C-band scatterometer measurements and L-band radiometer measurements at these grasslands.

Description: There have been many significant improvements in the public access to the Space Shuttle Earth Observations Photography Database. New information is provided for the user community on the recently released videodisc of this database. Topics covered included the following: earlier attempts; our first laser videodisc in 1992; the new laser videodisc in 1994; and electronic database access.

Description: The Challenge Awards are designed to provide a unique perspective to students gifted in the arts and humanities from which to understand scientific endeavor by giving students an opportunity to participate in an ongoing research project. In the graduate program, seven students who had participated in previous Challenge Awards programs were selected to help develop the tools for Earth observations for the astronauts on the Space Radar Laboratory (SRL) missions. The goal of the Challenge Awards program was to prepare a training manual for the astronauts on the SRL missions. This paper describes the observations to be made by the astronauts on the SRL missions. The emphasis is on the dynamic seasonal features of the Earth's surface and atmosphere which justify the need for more than one flight of the SRL. Complete notebooks of the sites, global seasonal patterns, examples of radar and the Measurement of Air Pollution from Satellites data, and shuttle photographs have been given to each of the SRL crews.

Description: MacSigma0 is an interactive tool for the Macintosh which allows you to display and make computations from radar data collected by the following sensors: the JPL AIRSAR, ERS-1, JERS-1, and Magellan. The JPL AIRSAR system is a multi-polarimetric airborne synthetic aperture radar developed and operated by the Jet Propulsion Laboratory. It includes the single-frequency L-band sensor mounted on the NASA CV990 aircraft and its replacement, the multi-frequency P-, L-, and C-band sensors mounted on the NASA DC-8. MacSigma0 works with data in the standard JPL AIRSAR output product format, the compressed Stokes matrix format. ERS-1 and JERS-1 are single-frequency, single-polarization spaceborne synthetic aperture radars launched by the European Space Agency and NASDA respectively. To be usable by MacSigma0, The data must have been processed at the Alaska SAR Facility and must be in the "low-resolution" format. Magellan is a spacecraft mission to map the surface of Venus with imaging radar. The project is managed by the Jet Propulsion Laboratory. The spacecraft carries a single-frequency, single-polarization synthetic aperture radar. MacSigma0 works with framelets of the standard MIDR CD-ROM data products. MacSigma0 provides four basic functions: synthesis of images (if necessary), statistical analysis of selected areas, analysis of corner reflectors as a calibration measure (if appropriate and possible), and informative mouse tracking. For instance, the JPL AIRSAR data can be used to synthesize a variety of images such as a total power image. The total power image displays the sum of the polarized and unpolarized components of the backscatter for each pixel. Other images which can be synthesized are HH, HV, VV, RL, RR, HHVV*, HHHV*, HVVV*, HHVV* phase and correlation coefficient images. For the complex and phase images, phase is displayed using color and magnitude is displayed using intensity. MacSigma0 can also be used to compute statistics from within a selected area. The statistics computed depend on the image type. For JPL AIRSAR data, the HH, HV, VV, HHVV* phase, and correlation coefficient means and standard deviation measures are calculated. The mean, relative standard deviation, minimum, and maximum values are calculated for all other data types. A histogram of the selected area is also calculated and displayed. The selected area can be rectangular, linear, or polygonal in shape. The user is allowed to select multiple rectangular areas, but not multiple linear or polygonal areas. The statistics and histogram are displayed to the user and can either be printed or saved as a text file. MacSigma0 can also be used to analyze corner reflectors as a measure of the calibration for JPL AIRSAR, ERS-1, and JERS-1 data types. It computes a theoretical radar cross section and the actual radar cross section for a selected trihedral corner reflector. The theoretical cross section, measured cross section, their ratio in dBs, and other information are displayed to the user and can be saved into a text file. For ERS-1, JERS-1, and Magellan data, MacSigma0 simultaneously displays pixel location in data coordinates and in latitude, longitude coordinates. It also displays sigma0, the incidence angle (for Magellan data), the original pixel value (for Magellan data), and the noise power value (for ERS-1 and JERS-1 data). Grey scale computed images can be saved in a byte format (a headerless format which saves the image as a string of byte values) or a PICT format (a standard format readable by other image processing programs for the Macintosh). Images can also be printed. MacSigma0 is written in C-language for use on Macintosh series computers. The minimum configuration requirements for MacSigma0 are System 6.0, Finder 6.1, 1Mb of RAM, and at least a 4-bit color or grey-scale graphics display. MacSigma0 is also System 7 compatible. To compile the source code, Apple's Macintosh Programmers Workbench (MPW) 3.2 and the MPW C language compiler version 3.2 are required. The source code will not compile with a later version of the compiler; however, the compiled application which will run under the minimum hardware configuration is provided on the distribution medium. In addition, the distribution media includes an executable which runs significantly faster but requires a 68881 compatible math coprocessor and a 68020 compatible CPU. Since JPL AIRSAR data files can be very large, it is often desirable to reduce the size of a data file before transferring it to the Macintosh for use in MacSigma0. A small FORTRAN program which can be used for this purpose is included on the distribution media. MacSigma0 will print statistics on any output device which supports QuickDraw, and it will print images on any device which supports QuickDraw or PostScript. The standard distribution medium for MacSigma0 is a set of five 1.4Mb Macintosh format diskettes. This program was developed in 1992 and is a copyrighted work with all copyright vested in NASA. Version 4.2 of MacSigma0 was released in 1993.

Description: VICAR (Video Image Communication and Retrieval) is a general purpose image processing software system that has been under continuous development since the late 1960's. Originally intended for data from the NASA Jet Propulsion Laboratory's unmanned planetary spacecraft, VICAR is now used for a variety of other applications including biomedical image processing, cartography, earth resources, and geological exploration. The development of this newest version of VICAR emphasized a standardized, easily-understood user interface, a shield between the user and the host operating system, and a comprehensive array of image processing capabilities. Structurally, VICAR can be divided into roughly two parts; a suite of applications programs and an executive which serves as the interfaces between the applications, the operating system, and the user. There are several hundred applications programs ranging in function from interactive image editing, data compression/decompression, and map projection, to blemish, noise, and artifact removal, mosaic generation, and pattern recognition and location. An information management system designed specifically for handling image related data can merge image data with other types of data files. The user accesses these programs through the VICAR executive, which consists of a supervisor and a run-time library. From the viewpoint of the user and the applications programs, the executive is an environment that is independent of the operating system. VICAR does not replace the host computer's operating system; instead, it overlays the host resources. The core of the executive is the VICAR Supervisor, which is based on NASA Goddard Space Flight Center's Transportable Applications Executive (TAE). Various modifications and extensions have been made to optimize TAE for image processing applications, resulting in a user friendly environment. The rest of the executive consists of the VICAR Run-Time Library, which provides a set of subroutines (image I/O, label I/O, parameter I/O, etc.) to facilitate image processing and provide the fastest I/O possible while maintaining a wide variety of capabilities. The run-time library also includes the Virtual Raster Display Interface (VRDI) which allows display oriented applications programs to be written for a variety of display devices using a set of common routines. (A display device can be any frame-buffer type device which is attached to the host computer and has memory planes for the display and manipulation of images. A display device may have any number of separate 8-bit image memory planes (IMPs), a graphics overlay plane, pseudo-color capabilities, hardware zoom and pan, and other features). The VRDI supports the following display devices: VICOM (Gould/Deanza) IP8500, RAMTEK RM-9465, ADAGE (Ikonas) IK3000 and the International Imaging Systems IVAS. VRDI's purpose is to provide a uniform operating environment not only for an application programmer, but for the user as well. The programmer is able to write programs without being concerned with the specifics of the device for which the application is intended. The VICAR Interactive Display Subsystem (VIDS) is a collection of utilities for easy interactive display and manipulation of images on a display device. VIDS has characteristics of both the executive and an application program, and offers a wide menu of image manipulation options. VIDS uses the VRDI to communicate with display devices. The first step in using VIDS to analyze and enhance an image (one simple example of VICAR's numerous capabilities) is to examine the histogram of the image. The histogram is a plot of frequency of occurrence for each pixel value (0 - 255) loaded in the image plane. If, for example, the histogram shows that there are no pixel values below 64 or above 192, the histogram can be "stretched" so that the value of 64 is mapped to zero and 192 is mapped to 255. Now the user can use the full dynamic range of the display device to display the data and better see its contents. Another example of a VIDS procedure is the JMOVIE command, which allows the user to run animations interactively on the display device. JMOVIE uses the concept of "frames", which are the individual frames which comprise the animation to be viewed. The user loads images into the frames after the size and number of frames has been selected. VICAR's source languages are primarily FORTRAN and C, with some VAX Assembler and array processor code. The VICAR run-time library is designed to work equally easily from either FORTRAN or C. The program was implemented on a DEC VAX series computer operating under VMS 4.7. The virtual memory required is 1.5MB. Approximately 180,000 blocks of storage are needed for the saveset. VICAR (version 2.3A/3G/13H) is a copyrighted work with all copyright vested in NASA and is available by license for a period of ten (10) years to approved licensees. This program was developed in 1989.

Description: The Interactive Image Display Program (IMDISP) is an interactive image display utility for the IBM Personal Computer (PC, XT and AT) and compatibles. Until recently, efforts to utilize small computer systems for display and analysis of scientific data have been hampered by the lack of sufficient data storage capacity to accomodate large image arrays. Most planetary images, for example, require nearly a megabyte of storage. The recent development of the "CDROM" (Compact Disk Read-Only Memory) storage technology makes possible the storage of up to 680 megabytes of data on a single 4.72-inch disk. IMDISP was developed for use with the CDROM storage system which is currently being evaluated by the Planetary Data System. The latest disks to be produced by the Planetary Data System are a set of three disks containing all of the images of Uranus acquired by the Voyager spacecraft. The images are in both compressed and uncompressed format. IMDISP can read the uncompressed images directly, but special software is provided to decompress the compressed images, which can not be processed directly. IMDISP can also display images stored on floppy or hard disks. A digital image is a picture converted to numerical form so that it can be stored and used in a computer. The image is divided into a matrix of small regions called picture elements, or pixels. The rows and columns of pixels are called "lines" and "samples", respectively. Each pixel has a numerical value, or DN (data number) value, quantifying the darkness or brightness of the image at that spot. In total, each pixel has an address (line number, sample number) and a DN value, which is all that the computer needs for processing. DISPLAY commands allow the IMDISP user to display all or part of an image at various positions on the display screen. The user may also zoom in and out from a point on the image defined by the cursor, and may pan around the image. To enable more or all of the original image to be displayed on the screen at once, the image can be "subsampled." For example, if the image were subsampled by a factor of 2, every other pixel from every other line would be displayed, starting from the upper left corner of the image. Any positive integer may be used for subsampling. The user may produce a histogram of an image file, which is a graph showing the number of pixels per DN value, or per range of DN values, for the entire image. IMDISP can also plot the DN value versus pixels along a line between two points on the image. The user can "stretch" or increase the contrast of an image by specifying low and high DN values; all pixels with values lower than the specified "low" will then become black, and all pixels higher than the specified "high" value will become white. Pixels between the low and high values will be evenly shaded between black and white. IMDISP is written in a modular form to make it easy to change it to work with different display devices or on other computers. The code can also be adapted for use in other application programs. There are device dependent image display modules, general image display subroutines, image I/O routines, and image label and command line parsing routines. The IMDISP system is written in C-language (94%) and Assembler (6%). It was implemented on an IBM PC with the MS DOS 3.21 operating system. IMDISP has a memory requirement of about 142k bytes. IMDISP was developed in 1989 and is a copyrighted work with all copyright vested in NASA. Additional planetary images can be obtained from the National Space Science Data Center at (301) 286-6695.

Description: PC-SEAPAK is a user-interactive satellite data analysis software package specifically developed for oceanographic research. The program is used to process and interpret data obtained from the Nimbus-7/Coastal Zone Color Scanner (CZCS), and the NOAA Advanced Very High Resolution Radiometer (AVHRR). PC-SEAPAK is a set of independent microcomputer-based image analysis programs that provide the user with a flexible, user-friendly, standardized interface, and facilitates relatively low-cost analysis of oceanographic satellite data. Version 4.0 includes 114 programs. PC-SEAPAK programs are organized into categories which include CZCS and AVHRR level-1 ingest, level-2 analyses, statistical analyses, data extraction, remapping to standard projections, graphics manipulation, image board memory manipulation, hardcopy output support and general utilities. Most programs allow user interaction through menu and command modes and also by the use of a mouse. Most programs also provide for ASCII file generation for further analysis in spreadsheets, graphics packages, etc. The CZCS scanning radiometer aboard the NIMBUS-7 satellite was designed to measure the concentration of photosynthetic pigments and their degradation products in the ocean. AVHRR data is used to compute sea surface temperatures and is supported for the NOAA 6, 7, 8, 9, 10, 11, and 12 satellites. The CZCS operated from November 1978 to June 1986. CZCS data may be obtained free of charge from the CZCS archive at NASA/Goddard Space Flight Center. AVHRR data may be purchased through NOAA's Satellite Data Service Division. Ordering information is included in the PC-SEAPAK documentation. Although PC-SEAPAK was developed on a COMPAQ Deskpro 386/20, it can be run on most 386-compatible computers with an AT bus, EGA controller, Intel 80387 coprocessor, and MS-DOS 3.3 or higher. A Matrox MVP-AT image board with appropriate monitor and cables is also required. Note that the authors have received some reports of incompatibilities between the MVP-AT image board and ZENITH computers. Also, the MVP-AT image board is not necessarily compatible with 486-based systems; users of 486-based systems should consult with Matrox about compatibility concerns. Other PC-SEAPAK requirements include a Microsoft mouse (serial version), 2Mb RAM, and 100Mb hard disk space. For data ingest and backup, 9-track tape, 8mm tape and optical disks are supported and recommended. PC-SEAPAK has been under development since 1988. Version 4.0 was updated in 1992, and is distributed without source code. It is available only as a set of 36 1.2Mb 5.25 inch IBM MS-DOS format diskettes. PC-SEAPAK is a copyrighted product with all copyright vested in the National Aeronautics and Space Administration. Phar Lap's DOS_Extender run-time version is integrated into several of the programs; therefore, the PC-SEAPAK programs may not be duplicated. Three of the distribution diskettes contain DOS_Extender files. One of the distribution diskettes contains Media Cybernetics' HALO88 font files, also licensed by NASA for dissemination but not duplication. IBM is a registered trademark of International Business Machines. MS-DOS is a registered trademark of Microsoft Corporation. HALO88 is a registered trademark of Media Cybernetics, but the product was discontinued in 1991.

Description: The Land Analysis System (LAS) is an image analysis system designed to manipulate and analyze digital data in raster format and provide the user with a wide spectrum of functions and statistical tools for analysis. LAS offers these features under VMS with optional image display capabilities for IVAS and other display devices as well as the X-Windows environment. LAS provides a flexible framework for algorithm development as well as for the processing and analysis of image data. Users may choose between mouse-driven commands or the traditional command line input mode. LAS functions include supervised and unsupervised image classification, film product generation, geometric registration, image repair, radiometric correction and image statistical analysis. Data files accepted by LAS include formats such as Multi-Spectral Scanner (MSS), Thematic Mapper (TM) and Advanced Very High Resolution Radiometer (AVHRR). The enhanced geometric registration package now includes both image to image and map to map transformations. The over 200 LAS functions fall into image processing scenario categories which include: arithmetic and logical functions, data transformations, fourier transforms, geometric registration, hard copy output, image restoration, intensity transformation, multispectral and statistical analysis, file transfer, tape profiling and file management among others. Internal improvements to the LAS code have eliminated the VAX VMS dependencies and improved overall system performance. The maximum LAS image size has been increased to 20,000 lines by 20,000 samples with a maximum of 256 bands per image. The catalog management system used in earlier versions of LAS has been replaced by a more streamlined and maintenance-free method of file management. This system is not dependent on VAX/VMS and relies on file naming conventions alone to allow the use of identical LAS file names on different operating systems. While the LAS code has been improved, the original capabilities of the system have been preserved. These include maintaining associated image history, session logging, and batch, asynchronous and interactive mode of operation. The LAS application programs are integrated under version 4.1 of an interface called the Transportable Applications Executive (TAE). TAE 4.1 has four modes of user interaction: menu, direct command, tutor (or help), and dynamic tutor. In addition TAE 4.1 allows the operation of LAS functions using mouse-driven commands under the TAE-Facelift environment provided with TAE 4.1. These modes of operation allow users, from the beginner to the expert, to exercise specific application options. LAS is written in C-language and FORTRAN 77 for use with DEC VAX computers running VMS with approximately 16Mb of physical memory. This program runs under TAE 4.1. Since TAE 4.1 is not a current version of TAE, TAE 4.1 is included within the LAS distribution. Approximately 130,000 blocks (65Mb) of disk storage space are necessary to store the source code and files generated by the installation procedure for LAS and 44,000 blocks (22Mb) of disk storage space are necessary for TAE 4.1 installation. The only other dependencies for LAS are the subroutine libraries for the specific display device(s) that will be used with LAS/DMS (e.g. X-Windows and/or IVAS). The standard distribution medium for LAS is a set of two 9~track 6250 BPI magnetic tapes in DEC VAX BACKUP format. It is also available on a set of two TK50 tape cartridges in DEC VAX BACKUP format. This program was developed in 1986 and last updated in 1992.

Description: Spaceborne Synthetic Aperture Radar (SAR) images are useful for planetary mapping and Earth sciences investigations. However, swath widths rarely exceed 100 Kilometers, and images must be patched together to create a mosaic in order to analyze larger areas. The primary function of this program is to generate large digital mosaics of SAR imagery without manually marked tiepoints. MOSK can produce multiframe mosaics by combining images in the along-track, adjacent cross-track swaths, or ascending and descending passes. Geocoded map registered images, such as the ones produced by MAPJTC (NPO-17718), are required as input. The output is a geocoded mosaic on a standard map grid which permits easy registration with other geocoded data sets. Mosaicking of geocoded SAR imagery involves three steps. First, a match point is selected at the center of the overlapping area, then an image patch around the match point is extracted from both images and cross-correlation is done on this area. Then, images with their refined match points are merged together to form a mosaic. To handle the large data volume of overlapping intermediate stages, large mosaics are divided into equal size quadrants with each quadrant cut from an intermediate mosaic. The full mosaic can then be assembled from the individual quadrants. Finally, radiometric disparities at the image seams are smoothed by a "feathering" technique. The automatic mosaic system generates output with minimal operator interaction. However, manual tiepointing is required in cases of a large registration error or two images with smooth surfaces such as ocean images. MOSK is implemented on a DEC VAX 11/785 running VMS 4.5. Most subroutines are in FORTRAN, but three are in MAXL and one is in APAL. The program requires 1 Mb of memory and a Floating Point Systems AP-5210 array processor. The system memory usage is approximately 1000 pages and the requirement of page file size is 2000 blocks. MOSK was developed in 1988.

Description: MAPJTC was designed to rectify and transform the standard image output of the digital Synthetic Aperture Radar (SAR) correlator into a geocoded map registered image without operator interaction or manual tiepointing. This is accomplished by modeling the distortions and predicting the pixel displacements based on platform and radar parameters. The map projection implemented in MAPJTC is the Universal Transverse Mercator. Since the re-sampling operation is independent of the transformation data generation, other cartographic projections can be implemented with few software modifications. MAPJTC makes a precise determination of the geodetic location of an arbitrary pixel within the image frame based on the simultaneous solution of a set of earth model equations, SAR Doppler equations, and SAR range equations that identify the slant range from the sensor to the target at a specific image pixel. Based on a table of geodetic coordinates of the image pixels, the image is then mapped onto the desired cartographic projection by applying the appropriate transformation equations. Typically, mapping involves a two-dimensional re-sampling and is very computationally intensive. MAPJTC reduces the procedure to two one-dimensional passes, which saves computer time. Geocoding transforms the rectified image into a grid defined by a specific map projection. (The image is rotated and rectified to match the map projection.) Again, the two dimensional re-sampling process can be separated into two one-dimensional re-sampling processes. Optionally, MAPJTC can correct terrain-induced distortions in SAR imagery when a digital elevation map is available. MAPJTC was developed on a DEC VAX 11/785 under VMS 4.5. The program is written in FORTRAN (84%), APAL (2%), and MAXL (14%). It requires 6Mb of memory and a Floating Point Systems AP-5210 Array Processor equipped with 1Mb of memory. MAPJTC can run interactively or as a batch job. MAPJTC was developed in 1987.

Description: The Spectral Analysis Manager (SPAM) was developed to allow easy qualitative analysis of multi-dimensional imaging spectrometer data. Imaging spectrometers provide sufficient spectral sampling to define unique spectral signatures on a per pixel basis. Thus direct material identification becomes possible for geologic studies. SPAM provides a variety of capabilities for carrying out interactive analysis of the massive and complex datasets associated with multispectral remote sensing observations. In addition to normal image processing functions, SPAM provides multiple levels of on-line help, a flexible command interpretation, graceful error recovery, and a program structure which can be implemented in a variety of environments. SPAM was designed to be visually oriented and user friendly with the liberal employment of graphics for rapid and efficient exploratory analysis of imaging spectrometry data. SPAM provides functions to enable arithmetic manipulations of the data, such as normalization, linear mixing, band ratio discrimination, and low-pass filtering. SPAM can be used to examine the spectra of an individual pixel or the average spectra over a number of pixels. SPAM also supports image segmentation, fast spectral signature matching, spectral library usage, mixture analysis, and feature extraction. High speed spectral signature matching is performed by using a binary spectral encoding algorithm to separate and identify mineral components present in the scene. The same binary encoding allows automatic spectral clustering. Spectral data may be entered from a digitizing tablet, stored in a user library, compared to the master library containing mineral standards, and then displayed as a timesequence spectral movie. The output plots, histograms, and stretched histograms produced by SPAM can be sent to a lineprinter, stored as separate RGB disk files, or sent to a Quick Color Recorder. SPAM is written in C for interactive execution and is available for two different machine environments. There is a DEC VAX/VMS version with a central memory requirement of approximately 242K of 8 bit bytes and a machine independent UNIX 4.2 version. The display device currently supported is the Raster Technologies display processor. Other 512 x 512 resolution color display devices, such as De Anza, may be added with minor code modifications. This program was developed in 1986.

Description: The Spectral Analysis Manager (SPAM) was developed to allow easy qualitative analysis of multi-dimensional imaging spectrometer data. Imaging spectrometers provide sufficient spectral sampling to define unique spectral signatures on a per pixel basis. Thus direct material identification becomes possible for geologic studies. SPAM provides a variety of capabilities for carrying out interactive analysis of the massive and complex datasets associated with multispectral remote sensing observations. In addition to normal image processing functions, SPAM provides multiple levels of on-line help, a flexible command interpretation, graceful error recovery, and a program structure which can be implemented in a variety of environments. SPAM was designed to be visually oriented and user friendly with the liberal employment of graphics for rapid and efficient exploratory analysis of imaging spectrometry data. SPAM provides functions to enable arithmetic manipulations of the data, such as normalization, linear mixing, band ratio discrimination, and low-pass filtering. SPAM can be used to examine the spectra of an individual pixel or the average spectra over a number of pixels. SPAM also supports image segmentation, fast spectral signature matching, spectral library usage, mixture analysis, and feature extraction. High speed spectral signature matching is performed by using a binary spectral encoding algorithm to separate and identify mineral components present in the scene. The same binary encoding allows automatic spectral clustering. Spectral data may be entered from a digitizing tablet, stored in a user library, compared to the master library containing mineral standards, and then displayed as a timesequence spectral movie. The output plots, histograms, and stretched histograms produced by SPAM can be sent to a lineprinter, stored as separate RGB disk files, or sent to a Quick Color Recorder. SPAM is written in C for interactive execution and is available for two different machine environments. There is a DEC VAX/VMS version with a central memory requirement of approximately 242K of 8 bit bytes and a machine independent UNIX 4.2 version. The display device currently supported is the Raster Technologies display processor. Other 512 x 512 resolution color display devices, such as De Anza, may be added with minor code modifications. This program was developed in 1986.

Description: The Composite Plate Buckling Analysis Program (COMPPAP) was written to help engineers determine buckling loads of orthotropic (or isotropic) irregularly shaped plates without requiring hand calculations from design curves or extensive finite element modeling. COMPPAP is a one element finite element program that utilizes high-order displacement functions. The high order of the displacement functions enables the user to produce results more accurate than traditional h-finite elements. This program uses these high-order displacement functions to perform a plane stress analysis of a general plate followed by a buckling calculation based on the stresses found in the plane stress solution. The current version assumes a flat plate (constant thickness) subject to a constant edge load (normal or shear) on one or more edges. COMPPAP uses the power method to find the eigenvalues of the buckling problem. The power method provides an efficient solution when only one eigenvalue is desired. Once the eigenvalue is found, the eigenvector, which corresponds to the plate buckling mode shape, results as a by-product. A positive feature of the power method is that the dominant eigenvalue is the first found, which is this case is the plate buckling load. The reported eigenvalue expresses a load factor to induce plate buckling. COMPPAP is written in ANSI FORTRAN 77. Two machine versions are available from COSMIC: a PC version (MSC-22428), which is for IBM PC 386 series and higher computers and compatibles running MS-DOS; and a UNIX version (MSC-22286). The distribution medium for both machine versions includes source code for both single and double precision versions of COMPPAP. The PC version includes source code which has been optimized for implementation within DOS memory constraints as well as sample executables for both the single and double precision versions of COMPPAP. The double precision versions of COMPPAP have been successfully implemented on an IBM PC 386 compatible running MS-DOS, a Sun4 series computer running SunOS, an HP-9000 series computer running HP-UX, and a CRAY X-MP series computer running UNICOS. COMPPAP requires 1Mb of RAM and the BLAS and LINPACK math libraries, which are included on the distribution medium. The COMPPAP documentation provides instructions for using the commercial post-processing package PATRAN for graphical interpretation of COMPPAP output. The UNIX version includes two electronic versions of the documentation: one in LaTex format and one in PostScript format. The standard distribution medium for the PC version (MSC-22428) is a 5.25 inch 1.2Mb MS-DOS format diskette. The standard distribution medium for the UNIX version (MSC-22286) is a .25 inch streaming magnetic tape cartridge (Sun QIC-24) in UNIX tar format. For the UNIX version, alternate distribution media and formats are available upon request. COMPPAP was developed in 1992.

Description: The Composite Plate Buckling Analysis Program (COMPPAP) was written to help engineers determine buckling loads of orthotropic (or isotropic) irregularly shaped plates without requiring hand calculations from design curves or extensive finite element modeling. COMPPAP is a one element finite element program that utilizes high-order displacement functions. The high order of the displacement functions enables the user to produce results more accurate than traditional h-finite elements. This program uses these high-order displacement functions to perform a plane stress analysis of a general plate followed by a buckling calculation based on the stresses found in the plane stress solution. The current version assumes a flat plate (constant thickness) subject to a constant edge load (normal or shear) on one or more edges. COMPPAP uses the power method to find the eigenvalues of the buckling problem. The power method provides an efficient solution when only one eigenvalue is desired. Once the eigenvalue is found, the eigenvector, which corresponds to the plate buckling mode shape, results as a by-product. A positive feature of the power method is that the dominant eigenvalue is the first found, which is this case is the plate buckling load. The reported eigenvalue expresses a load factor to induce plate buckling. COMPPAP is written in ANSI FORTRAN 77. Two machine versions are available from COSMIC: a PC version (MSC-22428), which is for IBM PC 386 series and higher computers and compatibles running MS-DOS; and a UNIX version (MSC-22286). The distribution medium for both machine versions includes source code for both single and double precision versions of COMPPAP. The PC version includes source code which has been optimized for implementation within DOS memory constraints as well as sample executables for both the single and double precision versions of COMPPAP. The double precision versions of COMPPAP have been successfully implemented on an IBM PC 386 compatible running MS-DOS, a Sun4 series computer running SunOS, an HP-9000 series computer running HP-UX, and a CRAY X-MP series computer running UNICOS. COMPPAP requires 1Mb of RAM and the BLAS and LINPACK math libraries, which are included on the distribution medium. The COMPPAP documentation provides instructions for using the commercial post-processing package PATRAN for graphical interpretation of COMPPAP output. The UNIX version includes two electronic versions of the documentation: one in LaTex format and one in PostScript format. The standard distribution medium for the PC version (MSC-22428) is a 5.25 inch 1.2Mb MS-DOS format diskette. The standard distribution medium for the UNIX version (MSC-22286) is a .25 inch streaming magnetic tape cartridge (Sun QIC-24) in UNIX tar format. For the UNIX version, alternate distribution media and formats are available upon request. COMPPAP was developed in 1992.

Description: Structural flaws and cracks may grow under fatigue inducing loads and, upon reaching a critical size, cause structural failure to occur. The growth of these flaws and cracks may occur at load levels well below the ultimate load bearing capability of the structure. The Fatigue Crack Growth Computer Program, NASA/FLAGRO, was developed as an aid in predicting the growth of pre-existing flaws and cracks in structural components of space systems. The earlier version of the program, FLAGRO4, was the primary analysis tool used by Rockwell International and the Shuttle subcontractors for fracture control analysis on the Space Shuttle. NASA/FLAGRO is an enhanced version of the program and incorporates state-of-the-art improvements in both fracture mechanics and computer technology. NASA/FLAGRO provides the fracture mechanics analyst with a computerized method of evaluating the "safe crack growth life" capabilities of structural components. NASA/FLAGRO could also be used to evaluate the damage tolerance aspects of a given structural design. The propagation of an existing crack is governed by the stress field in the vicinity of the crack tip. The stress intensity factor is defined in terms of the relationship between the stress field magnitude and the crack size. The propagation of the crack becomes catastrophic when the local stress intensity factor reaches the fracture toughness of the material. NASA/FLAGRO predicts crack growth using a two-dimensional model which predicts growth independently in two directions based on the calculation of stress intensity factors. The analyst can choose to use either a crack growth rate equation or a nonlinear interpolation routine based on tabular data. The growth rate equation is a modified Forman equation which can be converted to a Paris or Walker equation by substituting different values into the exponent. This equation provides accuracy and versatility and can be fit to data using standard least squares methods. Stress-intensity factor numerical values can be computed for making comparisons or checks of solutions. NASA/FLAGRO can check for failure of a part-through crack in the mode of a through crack when net ligament yielding occurs. NASA/FLAGRO has a number of special subroutines and files which provide enhanced capabilities and easy entry of data. These include crack case solutions, cyclic load spectrums, nondestructive examination initial flaw sizes, table interpolation, and material properties. The materials properties files are divided into two types, a user defined file and a fixed file. Data is entered and stored in the user defined file during program execution, while the fixed file contains already coded-in property value data for many different materials. Prompted input from CRT terminals consists of initial crack definition (which can be defined automatically), rate solution type, flaw type and geometry, material properties (if they are not in the built-in tables of material data), load spectrum data (if not included in the loads spectrum file), and design limit stress levels. NASA/FLAGRO output includes an echo of the input with any error or warning messages, the final crack size, whether or not critical crack size has been reached for the specified stress level, and a life history profile of the crack propagation. NASA/FLAGRO is modularly designed to facilitate revisions and operation on minicomputers. The program was implemented on a DEC VAX 11/780 with the VMS operating system. NASA/FLAGRO is written in FORTRAN77 and has a memory requirement of 1.4 MB. The program was developed in 1986.

Description: The Science and Technology Laboratory Applications Software (ELAS) was originally designed to analyze and process digital imagery data, specifically remotely-sensed scanner data. This capability includes the processing of Landsat multispectral data; aircraft-acquired scanner data; digitized topographic data; and numerous other ancillary data, such as soil types and rainfall information, that can be stored in digitized form. ELAS has the subsequent capability to geographically reference this data to dozens of standard, as well as user created projections. As an integrated image processing system, ELAS offers the user of remotely-sensed data a wide range of capabilities in the areas of land cover analysis and general purpose image analysis. ELAS is designed for flexible use and operation and includes its own FORTRAN operating subsystem and an expandable set of FORTRAN application modules. Because all of ELAS resides in one "logical" FORTRAN program, data inputs and outputs, directives, and module switching are convenient for the user. There are over 230 modules presently available to aid the user in performing a wide range of land cover analyses and manipulation. The file management modules enable the user to allocate, define, access, and specify usage for all types of files (ELAS files, subfiles, external files etc.). Various other modules convert specific types of satellite, aircraft, and vector-polygon data into files that can be used by other ELAS modules. The user also has many module options which aid in displaying image data, such as magnification/reduction of the display; true color display; and several memory functions. Additional modules allow for the building and manipulation of polygonal areas of the image data. Finally, there are modules which allow the user to select and classify the image data. An important feature of the ELAS subsystem is that its structure allows new applications modules to be easily integrated in the future. ELAS has as a standard the flexibility to process data elements exceeding 8 bits in length, including floating point (noninteger) elements and 16 or 32 bit integers. Thus it is able to analyze and process "non-standard" nonimage data. The VAX (ERL-10017) and Concurrent (ERL-10013) versions of ELAS 9.0 are written in FORTRAN and ASSEMBLER for DEC VAX series computers running VMS and Concurrent computers running MTM. The Sun (SSC-00019), Masscomp (SSC-00020), and Silicon Graphics (SSC-00021) versions of ELAS 9.0 are written in FORTRAN 77 and C-LANGUAGE for Sun4 series computers running SunOS, Masscomp computers running UNIX, and Silicon Graphics IRIS computers running IRIX. The Concurrent version requires at least 15 bit addressing and a direct memory access channel. The VAX and Concurrent versions of ELAS both require floating-point hardware, at least 1Mb of RAM, and approximately 70Mb of disk space. Both versions also require a COMTAL display device in order to display images. For the Sun, Masscomp, and Silicon Graphics versions of ELAS, the disk storage required is approximately 115Mb, and a minimum of 8Mb of RAM is required for execution. The Sun version of ELAS requires either the X-Window System Version 11 Revision 4 or Sun OpenWindows Version 2. The Masscomp version requires a GA1000 display device and the associated "gp" library. The Silicon Graphics version requires Silicon Graphics' GL library. ELAS display functions will not work with a monochrome monitor. The standard distribution medium for the VAX version (ERL~10017) is a set of two 9-track 1600 BPI magnetic tapes in DEC VAX BACKUP format. This version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. The standard distribution medium for the Concurrent version (ERL-10013) is a set of two 9-track 1600 BPI magnetic tapes in Concurrent BACKUP format. The standard distribution medium for the Sun version (SSC-00019) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Masscomp version, (SSC-00020) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Silicon Graphics version (SSC-00021) is a .25 inch streaming magnetic IRIS tape cartridge in UNIX tar format. Version 9.0 was released in 1991. Sun4, SunOS, and Open Windows are trademarks of Sun Microsystems, Inc. MIT X Window System is licensed by Massachusetts Institute of Technology.

Description: The Science and Technology Laboratory Applications Software (ELAS) was originally designed to analyze and process digital imagery data, specifically remotely-sensed scanner data. This capability includes the processing of Landsat multispectral data; aircraft-acquired scanner data; digitized topographic data; and numerous other ancillary data, such as soil types and rainfall information, that can be stored in digitized form. ELAS has the subsequent capability to geographically reference this data to dozens of standard, as well as user created projections. As an integrated image processing system, ELAS offers the user of remotely-sensed data a wide range of capabilities in the areas of land cover analysis and general purpose image analysis. ELAS is designed for flexible use and operation and includes its own FORTRAN operating subsystem and an expandable set of FORTRAN application modules. Because all of ELAS resides in one "logical" FORTRAN program, data inputs and outputs, directives, and module switching are convenient for the user. There are over 230 modules presently available to aid the user in performing a wide range of land cover analyses and manipulation. The file management modules enable the user to allocate, define, access, and specify usage for all types of files (ELAS files, subfiles, external files etc.). Various other modules convert specific types of satellite, aircraft, and vector-polygon data into files that can be used by other ELAS modules. The user also has many module options which aid in displaying image data, such as magnification/reduction of the display; true color display; and several memory functions. Additional modules allow for the building and manipulation of polygonal areas of the image data. Finally, there are modules which allow the user to select and classify the image data. An important feature of the ELAS subsystem is that its structure allows new applications modules to be easily integrated in the future. ELAS has as a standard the flexibility to process data elements exceeding 8 bits in length, including floating point (noninteger) elements and 16 or 32 bit integers. Thus it is able to analyze and process "non-standard" nonimage data. The VAX (ERL-10017) and Concurrent (ERL-10013) versions of ELAS 9.0 are written in FORTRAN and ASSEMBLER for DEC VAX series computers running VMS and Concurrent computers running MTM. The Sun (SSC-00019), Masscomp (SSC-00020), and Silicon Graphics (SSC-00021) versions of ELAS 9.0 are written in FORTRAN 77 and C-LANGUAGE for Sun4 series computers running SunOS, Masscomp computers running UNIX, and Silicon Graphics IRIS computers running IRIX. The Concurrent version requires at least 15 bit addressing and a direct memory access channel. The VAX and Concurrent versions of ELAS both require floating-point hardware, at least 1Mb of RAM, and approximately 70Mb of disk space. Both versions also require a COMTAL display device in order to display images. For the Sun, Masscomp, and Silicon Graphics versions of ELAS, the disk storage required is approximately 115Mb, and a minimum of 8Mb of RAM is required for execution. The Sun version of ELAS requires either the X-Window System Version 11 Revision 4 or Sun OpenWindows Version 2. The Masscomp version requires a GA1000 display device and the associated "gp" library. The Silicon Graphics version requires Silicon Graphics' GL library. ELAS display functions will not work with a monochrome monitor. The standard distribution medium for the VAX version (ERL~10017) is a set of two 9-track 1600 BPI magnetic tapes in DEC VAX BACKUP format. This version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. The standard distribution medium for the Concurrent version (ERL-10013) is a set of two 9-track 1600 BPI magnetic tapes in Concurrent BACKUP format. The standard distribution medium for the Sun version (SSC-00019) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Masscomp version, (SSC-00020) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Silicon Graphics version (SSC-00021) is a .25 inch streaming magnetic IRIS tape cartridge in UNIX tar format. Version 9.0 was released in 1991. Sun4, SunOS, and Open Windows are trademarks of Sun Microsystems, Inc. MIT X Window System is licensed by Massachusetts Institute of Technology.

Description: IMAGEP is a FORTRAN computer algorithm containing various image processing, analysis, and enhancement functions. It is a keyboard-driven program organized into nine subroutines. Within the subroutines are other routines, also, selected via keyboard. Some of the functions performed by IMAGEP include digitization, storage and retrieval of images; image enhancement by contrast expansion, addition and subtraction, magnification, inversion, and bit shifting; display and movement of cursor; display of grey level histogram of image; and display of the variation of grey level intensity as a function of image position. This algorithm has possible scientific, industrial, and biomedical applications in material flaw studies, steel and ore analysis, and pathology, respectively. IMAGEP is written in VAX FORTRAN for DEC VAX series computers running VMS. The program requires the use of a Grinnell 274 image processor which can be obtained from Mark McCloud Associates, Campbell, CA. An object library of the required GMR series software is included on the distribution media. IMAGEP requires 1Mb of RAM for execution. The standard distribution medium for this program is a 1600 BPI 9~track magnetic tape in VAX FILES-11 format. It is also available on a TK50 tape cartridge in VAX FILES-11 format. This program was developed in 1991. DEC, VAX, VMS, and TK50 are trademarks of Digital Equipment Corporation.

Description: Delaminations near the outer surface of a laminate are susceptible to local buckling and buckling-induced delamination propagation when the laminate is subjected to transverse impact loading. This results in a loss of stiffness and strength. TRBUCKL is an unique dynamic delamination buckling and delamination propagation analysis capability that can be incorporated into the structural analysis program, NASTRAN. This capability will aid engineers in the design of structures incorporating composite laminates. The capability consists of: (1) a modification of the direct time integration solution sequence which provides a new analysis algorithm that can be used to predict delamination buckling in a laminate subjected to dynamic loading; and (2) a new method of modeling the composite laminate using plate bending elements and multipoint constraints. The capability now exists to predict the time at which the onset of dynamic delamination buckling occurs, the dynamic buckling mode shape, and the dynamic delamination strain energy release rate. A procedure file for NASTRAN, TRBUCKL predicts both impact induced buckling in composite laminates with initial delaminations and the strain energy release rate due to extension of the delamination. In addition, the file is useful in calculating the dynamic delamination strain energy release rate for a composite laminate under impact loading. This procedure simplifies the simulation of progressive crack extension. TRBUCKL has been incorporated into COSMIC NASTRAN. TRBUCKL is a DMAP Alter for NASTRAN. It is intended for use only with the COSMIC NASTRAN Direct Transient Analysis (RF 9) solution sequence. The program is available as a listing only. TRBUCKL was developed in 1987.

Description: Predictions of fatigue crack growth behavior can be made with the Fatigue Crack Growth Structural Analysis (FASTRAN II) computer program. As cyclic loads are applied to a selected crack configuration with an initial crack size, FASTRAN II predicts crack growth as a function of cyclic load history until either a desired crack size is reached or failure occurs. FASTRAN II is based on plasticity-induced crack-closure behavior of cracks in metallic materials and accounts for load-interaction effects, such as retardation and acceleration, under variable-amplitude loading. The closure model is based on the Dugdale model with modifications to allow plastically deformed material to be left along the crack surfaces as the crack grows. Plane stress and plane strain conditions, as well as conditions between these two, can be simulated in FASTRAN II by using a constraint factor on tensile yielding at the crack front to approximately account for three-dimensional stress states. FASTRAN II contains seventeen predefined crack configurations (standard laboratory fatigue crack growth rate specimens and many common crack configurations found in structures); and the user can define one additional crack configuration. The baseline crack growth rate properties (effective stress-intensity factor against crack growth rate) may be given in either equation or tabular form. For three-dimensional crack configurations, such as surface cracks or corner cracks at holes or notches, the fatigue crack growth rate properties may be different in the crack depth and crack length directions. Final failure of the cracked structure can be modelled with fracture toughness properties using either linear-elastic fracture mechanics (brittle materials), a two-parameter fracture criterion (brittle to ductile materials), or plastic collapse (extremely ductile materials). The crack configurations in FASTRAN II can be subjected to either constant-amplitude, variable-amplitude or spectrum loading. The applied loads may be either tensile or compressive. Several standardized aircraft flight-load histories, such as TWIST, Mini-TWIST, FALSTAFF, Inverted FALSTAFF, Felix and Gaussian, are included as options. FASTRAN II also includes two other methods that will help the user input spectrum load histories. The two methods are: (1) a list of stress points, and (2) a flight-by-flight history of stress points. Examples are provided in the user manual. Developed as a research program, FASTRAN II has successfully predicted crack growth in many metallic materials under various aircraft spectrum loading. A computer program DKEFF which is a part of the FASTRAN II package was also developed to analyze crack growth rate data from laboratory specimens to obtain the effective stress-intensity factor against crack growth rate relations used in FASTRAN II. FASTRAN II is written in standard FORTRAN 77. It has been successfully compiled and implemented on Sun4 series computers running SunOS and on IBM PC compatibles running MS-DOS using the Lahey F77L FORTRAN compiler. Sample input and output data are included with the FASTRAN II package. The UNIX version requires 660K of RAM for execution. The standard distribution medium for the UNIX version (LAR-14865) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. It is also available on a 3.5 inch diskette in UNIX tar format. The standard distribution medium for the MS-DOS version (LAR-14944) is a 5.25 inch 360K MS-DOS format diskette. The contents of the diskette are compressed using the PKWARE archiving tools. The utility to unarchive the files, PKUNZIP.EXE, is included. The program was developed in 1984 and revised in 1992. Sun4 and SunOS are trademarks of Sun Microsystems, Inc. IBM PC is a trademark of International Business Machines Corp. MS-DOS is a trademark of Microsoft, Inc. F77L is a trademark of the Lahey Computer Systems, Inc. UNIX is a registered trademark of AT&T Bell Laboratories. PKWARE and PKUNZIP are trademarks of PKWare, Inc.

Description: The Panel Analysis and Sizing Code (PASCO) was developed for the buckling and vibration analysis and sizing of prismatic structures having an arbitrary cross section. PASCO is primarily intended for analyzing and sizing stiffened panels made of laminated orthotropic materials and is of particular value in analyzing and sizing filamentary composite structures. When used in the analysis mode, PASCO calculates laminate stiffnesses, lamina stress and strains (including the effects of temperature and panel bending), buckling loads, vibration frequencies, and overall panel stiffness. When used in the sizing mode, PASCO adjusts sizing variables to provide a low-mass panel design that carries a set of specified loadings without exceeding buckling or material strength allowables and that meets other design requirements such as upper and lower bounds on sizing variables, upper and lower bounds on overall bending, extensional and shear stiffnesses, and lower bounds on vibration frequencies. Although emphasis in PASCO is placed on flat panels having several identical bays, the only restriction on configuration modeling is that the structure is assumed to be prismatic. In addition, it is assumed that loads and temperatures do not vary along the length of a panel. Because of their wide application in aerospace structures, stiffened panels are readily handled by PASCO. The panel cross section may be composed of an arbitrary assemblage of thin, flat, rectangular plate elements that are connected together along their longitudinal edges. Each plate element consists of a balanced symmetric laminate of any number of layers of orthotropic material. Any group of element widths, layer thicknesses, and layer orientation angles may be selected as sizing variables. Substructuring is available to increase the efficiency of the analysis and to simplify the modeling of complicated structures. The Macintosh version of PASCO includes an interactive, graphic preprocessor called MacPASCO. The main objective of MacPASCO is to make the use of PASCO faster, simpler, and less error-prone. By using a graphical user interface (GUI), MacPASCO simplifies the specification of panel geometry and reduces user input errors, thus making the modeling and analysis of panel designs more efficient. The user draws the initial structural geometry on the computer screen, then uses a combination of graphic and text inputs to: refine the structural geometry, specify information required for analysis such as panel load conditions, and define design variables and constraints for minimum-mass optimization. Composite panel design is an ideal application because the graphical user interface can: serve as a visual aid, eliminate the tedious aspects of text-based input, and eliminate many sources of input errors. The current version of MacPASCO does not implement all the modeling features of PASCO, but has been found to be sufficient for many users. Many difficulties common to text-based inputs are avoided because MacPASCO uses a GUI. First, the graphic displays eliminate syntax errors, like misplaced commas and incorrect command names, because there is no text-based syntax. Second, graphic displays allow the user to see the geometry as it is created and immediately detect and correct any errors. Third, MacPASCO's drawing tools have been designed to avoid modeling errors. Fourth, the graphic displays make revisions to existing structural designs much easier and less error-prone by eliminating the need for the user to conceptualize the text input as geometry. The user can work directly with the geometry displayed on the screen. Finally, MacPASCO automatically generates the correct PASCO input from the geometry displayed on the screen. This input file can be used with any machine version of PASCO to actually perform the analysis and sizing and to output results. The DEC VAX version of PASCO is written in FORTRAN IV for batch execution and has been implemented on a DEC VAX series computer. The Macintosh version of PASCO was developed for Macintosh II series computers with at least 2Mb of RAM running MPW Pascal 3.0 and Language Systems FORTRAN 2.0 under the MPW programming environment. It includes MPW compatible makefiles for compiling the source code. The Macintosh version uses input files compatible with versions of PASCO running on different platforms. MacPASCO is written in Macintosh Programmers Workbench 3.0, MPW Pascal 3.0, and MacAPP 2.0. The Pascal source code is included on the distribution diskette. MacAPP is a development library which is not included. MacPASCO requires a Mac Plus, SE/30, or MacII, IIx, IIcx, IIci, or IIfx running System 6.0 or greater. MacPASCO is System 7.0 compatible. A minimum of 2Mb of RAM is required for execution. The Macintosh version of PASCO is distributed on four 3.5 inch 800K Macintosh format diskettes. The DEC VAX version is distributed on a 9-track 1600 BPI magnetic tape. The PASCO program was developed in 1981, adapted to the DEC VAX in 1983 and to the Macintosh in 1991. MacPASCO was released in 1992.

Description: PLAN2D is a FORTRAN computer program for the plastic analysis of planar rigid frame structures. Given a structure and loading pattern as input, PLAN2D calculates the ultimate load that the structure can sustain before collapse. Element moments and plastic hinge rotations are calculated for the ultimate load. The location of hinges required for a collapse mechanism to form are also determined. The program proceeds in an iterative series of linear elastic analyses. After each iteration the resulting elastic moments in each member are compared to the reserve plastic moment capacity of that member. The member or members that have moments closest to their reserve capacity will determine the minimum load factor and the site where the next hinge is to be inserted. Next, hinges are inserted and the structural stiffness matrix is reformulated. This cycle is repeated until the structure becomes unstable. At this point the ultimate collapse load is calculated by accumulating the minimum load factor from each previous iteration and multiplying them by the original input loads. PLAN2D is based on the program STAN, originally written by Dr. E.L. Wilson at U.C. Berkeley. PLAN2D has several limitations: 1) Although PLAN2D will detect unloading of hinges it does not contain the capability to remove hinges; 2) PLAN2D does not allow the user to input different positive and negative moment capacities and 3) PLAN2D does not consider the interaction between axial and plastic moment capacity. Axial yielding and buckling is ignored as is the reduction in moment capacity due to axial load. PLAN2D is written in FORTRAN and is machine independent. It has been tested on an IBM PC and a DEC MicroVAX. The program was developed in 1988.

Description: Predictions of fatigue crack growth behavior can be made with the Fatigue Crack Growth Structural Analysis (FASTRAN II) computer program. As cyclic loads are applied to a selected crack configuration with an initial crack size, FASTRAN II predicts crack growth as a function of cyclic load history until either a desired crack size is reached or failure occurs. FASTRAN II is based on plasticity-induced crack-closure behavior of cracks in metallic materials and accounts for load-interaction effects, such as retardation and acceleration, under variable-amplitude loading. The closure model is based on the Dugdale model with modifications to allow plastically deformed material to be left along the crack surfaces as the crack grows. Plane stress and plane strain conditions, as well as conditions between these two, can be simulated in FASTRAN II by using a constraint factor on tensile yielding at the crack front to approximately account for three-dimensional stress states. FASTRAN II contains seventeen predefined crack configurations (standard laboratory fatigue crack growth rate specimens and many common crack configurations found in structures); and the user can define one additional crack configuration. The baseline crack growth rate properties (effective stress-intensity factor against crack growth rate) may be given in either equation or tabular form. For three-dimensional crack configurations, such as surface cracks or corner cracks at holes or notches, the fatigue crack growth rate properties may be different in the crack depth and crack length directions. Final failure of the cracked structure can be modelled with fracture toughness properties using either linear-elastic fracture mechanics (brittle materials), a two-parameter fracture criterion (brittle to ductile materials), or plastic collapse (extremely ductile materials). The crack configurations in FASTRAN II can be subjected to either constant-amplitude, variable-amplitude or spectrum loading. The applied loads may be either tensile or compressive. Several standardized aircraft flight-load histories, such as TWIST, Mini-TWIST, FALSTAFF, Inverted FALSTAFF, Felix and Gaussian, are included as options. FASTRAN II also includes two other methods that will help the user input spectrum load histories. The two methods are: (1) a list of stress points, and (2) a flight-by-flight history of stress points. Examples are provided in the user manual. Developed as a research program, FASTRAN II has successfully predicted crack growth in many metallic materials under various aircraft spectrum loading. A computer program DKEFF which is a part of the FASTRAN II package was also developed to analyze crack growth rate data from laboratory specimens to obtain the effective stress-intensity factor against crack growth rate relations used in FASTRAN II. FASTRAN II is written in standard FORTRAN 77. It has been successfully compiled and implemented on Sun4 series computers running SunOS and on IBM PC compatibles running MS-DOS using the Lahey F77L FORTRAN compiler. Sample input and output data are included with the FASTRAN II package. The UNIX version requires 660K of RAM for execution. The standard distribution medium for the UNIX version (LAR-14865) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. It is also available on a 3.5 inch diskette in UNIX tar format. The standard distribution medium for the MS-DOS version (LAR-14944) is a 5.25 inch 360K MS-DOS format diskette. The contents of the diskette are compressed using the PKWARE archiving tools. The utility to unarchive the files, PKUNZIP.EXE, is included. The program was developed in 1984 and revised in 1992. Sun4 and SunOS are trademarks of Sun Microsystems, Inc. IBM PC is a trademark of International Business Machines Corp. MS-DOS is a trademark of Microsoft, Inc. F77L is a trademark of the Lahey Computer Systems, Inc. UNIX is a registered trademark of AT&T Bell Laboratories. PKWARE and PKUNZIP are trademarks of PKWare, Inc.

Description: The SNIP program is a FORTRAN computer code that generates NASTRAN structural model thermal loads when given SINDA (or similar thermal model) temperature results. SNIP correlates thermal nodes to structural elements to interface SINDA finite difference thermal models with NASTRAN finite element structural models. Node-to-element correlation includes determining which SINDA nodes should be related to each NASTRAN element and calculating a weighing factor for temperatures associated with each element-related thermal node. SNIP provides structural model thermal loads that accurately reflect thermal model results while reducing the time required to interface thermal and structural models as compared to other methods. SNIP uses thermal model geometry to search the three-dimensional space around each structural element for the nearest thermal nodes. Thermal model geometry is the combination of standard thermal model temperature results from SINDA and structural model geometry from NASTRAN. Thermal and structural models must both be defined in the same, single Cartesian coordinate system. The thermal nodes located nearest each element are used to determine element temperature for thermal distortion and stress analysis. The program shapes the three-dimensional search region while the user controls the size. With these region specifications, the numerical coding of thermal nodes, and the structural element numbers; the code can provide for the separation of substructures during correlation. The input to SNIP contains a file of thermal model temperature results and a physical location of each thermal node in three-dimensional space, combined in a SNIP-unique format. The input also contains a standard NASTRAN input deck for a model made up of plate, shell, beam, and bar elements. SNIP supports the CTRIA, CQUAD, CBAR, and CBEAM elements of NASTRAN. The user adjusts the input parameters in the source code which control the node-to-element correlation. The program outputs NASTRAN element temperature load cards for each element and NASTRAN case control cards for each temperature load set. SNIP also outputs a list of elements that contains the numbers of the SINDA nodes related to each NASTRAN element and the weight that is given to each node in temperature calculations. SNIP is written in ANSI standard FORTRAN 77. The PC version requires a PC FORTRAN compiler and has compiled successfully using Lahey FORTRAN v. 3.0. A core memory of 300k is recommended. The program was developed in 1987.

Description: NPLOT is an interactive computer graphics program for plotting undeformed and deformed NASTRAN finite element models (FEMs). Although there are many commercial codes already available for plotting FEMs, these have limited use due to their cost, speed, and lack of features to view BAR elements. NPLOT was specifically developed to overcome these limitations. On a vector type graphics device the two best ways to show depth are by hidden line plotting or haloed line plotting. A hidden line algorithm generates views of models with all hidden lines removed, and a haloed line algorithm displays views with aft lines broken in order to show depth while keeping the entire model visible. A haloed line algorithm is especially useful for plotting models composed of many line elements and few surface elements. The most important feature of NPLOT is its ability to create both hidden line and haloed line views accurately and much more quickly than with any other existing hidden or haloed line algorithms. NPLOT is also capable of plotting a normal wire frame view to display all lines of a model. NPLOT is able to aid in viewing all elements, but it has special features not generally available for plotting BAR elements. These features include plotting of TRUE LENGTH and NORMALIZED offset vectors and orientation vectors. Standard display operations such as rotation and perspective are possible, but different view planes such as X-Y, Y-Z, and X-Z may also be selected. Another display option is the Z-axis cut which allows a portion of the fore part of the model to be cut away to reveal details of the inside of the model. A zoom function is available to terminals with a locator (graphics cursor, joystick, etc.). The user interface of NPLOT is designed to make the program quick and easy to use. A combination of menus and commands with help menus for detailed information about each command allows experienced users greater speed and efficiency. Once a plot is on the screen the interface becomes command driven, enabling the user to manipulate the display or execute a command without having to return to the menu. NPLOT is also able to plot deformed shapes allowing it to perform post-processing. The program can read displacements, either static displacements or eigenvectors, from a MSC/NASTRAN F06 file or a UAI/NASTRAN PRT file. The displacements are written into a unformatted scratch file where they are available for rapid access when the user wishes to display a deformed shape. All subcases or mode shapes can be read in at once. Then it is easy to enable the deformed shape, to change subcases or mode shapes and to change the scale factor for subsequent plots. NPLOT is written in VAX FORTRAN for DEC VAX series computers running VMS. As distributed, the NPLOT source code makes calls to the DI3000 graphics package from Precision Visuals; however, a set of interface routines is provided to translate the DI3000 calls into Tektronix PLOT10/TCS graphics library calls so that NPLOT can use the standard Tektronix 4010 which many PC terminal emulation software programs support. NPLOT is available in VAX BACKUP format on a 9-track 1600 BPI DEC VAX BACKUP format magnetic tape (standard media) or a TK50 tape cartridge. This program was developed in 1991. DEC, VAX, VMS, and TK50 are trademarks of Digital Equipment Corporation. Tektronix, PLOT10, and TCS are trademarks of Tektronix, Inc. DI3000 is a registered trademark of Precision Visuals, Inc. NASTRAN is a registered trademark of the National Aeronautics and Space Administration. MSC/ is a trademark of MacNeal-Schwendler Corporation. UAI is a trademark of Universal Analytics, Inc.

Description: LATDYN is a computer code for modeling the Large Angle Transient DYNamics of flexible articulating structures and mechanisms involving joints about which members rotate through large angles. LATDYN extends and brings together some of the aspects of Finite Element Structural Analysis, Multi-Body Dynamics, and Control System Analysis; three disciplines that have been historically separate. It combines significant portions of their distinct capabilities into one single analysis tool. The finite element formulation for flexible bodies in LATDYN extends the conventional finite element formulation by using a convected coordinate system for constructing the equation of motion. LATDYN's formulation allows for large displacements and rotations of finite elements subject to the restriction that deformations within each are small. Also, the finite element approach implemented in LATDYN provides a convergent path for checking solutions simply by increasing mesh density. For rigid bodies and joints LATDYN borrows extensively from methodology used in multi-body dynamics where rigid bodies may be defined and connected together through joints (hinges, ball, universal, sliders, etc.). Joints may be modeled either by constraints or by adding joint degrees of freedom. To eliminate error brought about by the separation of structural analysis and control analysis, LATDYN provides symbolic capabilities for modeling control systems which are integrated with the structural dynamic analysis itself. Its command language contains syntactical structures which perform symbolic operations which are also interfaced directly with the finite element structural model, bypassing the modal approximation. Thus, when the dynamic equations representing the structural model are integrated, the equations representing the control system are integrated along with them as a coupled system. This procedure also has the side benefit of enabling a dramatic simplification of the user interface for modeling control systems. Three FORTRAN computer programs, the LATDYN Program, the Preprocessor, and the Postprocessor, make up the collective LATDYN System. The Preprocessor translates user commands into a form which can be used while the LATDYN program provides the computational core. The Postprocessor allows the user to interactively plot and manage a database of LATDYN transient analysis results. It also includes special facilities for modeling control systems and for programming changes to the model which take place during analysis sequence. The documentation includes a Demonstration Problem Manual for the evaluation and verification of results and a Postprocessor guide. Because the program should be viewed as a byproduct of research on technology development, LATDYN's scope is limited. It does not have a wide library of finite elements, and 3-D Graphics are not available. Nevertheless, it does have a measure of "user friendliness". The LATDYN program was developed over a period of several years and was implemented on a CDC NOS/VE & Convex Unix computer. It is written in FORTRAN 77 and has a virtual memory requirement of 1.46 MB. The program was validated on a DEC MICROVAX operating under VMS 5.2.

Description: The increasing number of applications of fiber-reinforced composites in industry demands a detailed understanding of their material properties and behavior. A three-dimensional finite-element computer program called PAFAC (Plastic and Failure Analysis of Composites) has been developed for the elastic-plastic analysis of fiber-reinforced composite materials and structures. The evaluation of stresses and deformations at edges, cut-outs, and joints is essential in understanding the strength and failure for metal-matrix composites since the onset of plastic yielding starts very early in the loading process as compared to the composite's ultimate strength. Such comprehensive analysis can only be achieved by a finite-element program like PAFAC. PAFAC is particularly suited for the analysis of laminated metal-matrix composites. It can model the elastic-plastic behavior of the matrix phase while the fibers remain elastic. Since the PAFAC program uses a three-dimensional element, the program can also model the individual layers of the laminate to account for thickness effects. In PAFAC, the composite is modeled as a continuum reinforced by cylindrical fibers of vanishingly small diameter which occupy a finite volume fraction of the composite. In this way, the essential axial constraint of the phases is retained. Furthermore, the local stress and strain fields are uniform. The PAFAC finite-element solution is obtained using the displacement method. Solution of the nonlinear equilibrium equations is obtained with a Newton-Raphson iteration technique. The elastic-plastic behavior of composites consisting of aligned, continuous elastic filaments and an elastic-plastic matrix is described in terms of the constituent properties, their volume fractions, and mutual constraints between phases indicated by the geometry of the microstructure. The program uses an iterative procedure to determine the overall response of the laminate, then from the overall response determines the stress state in each phase of the composite material. Failure of the fibers or matrix within an element can also be modeled by PAFAC. PAFAC is written in FORTRAN IV for batch execution and has been implemented on a CDC CYBER 170 series computer with a segmented memory requirement of approximately 66K (octal) of 60 bit words. PAFAC was developed in 1982.

Description: The Panel Analysis and Sizing Code (PASCO) was developed for the buckling and vibration analysis and sizing of prismatic structures having an arbitrary cross section. PASCO is primarily intended for analyzing and sizing stiffened panels made of laminated orthotropic materials and is of particular value in analyzing and sizing filamentary composite structures. When used in the analysis mode, PASCO calculates laminate stiffnesses, lamina stress and strains (including the effects of temperature and panel bending), buckling loads, vibration frequencies, and overall panel stiffness. When used in the sizing mode, PASCO adjusts sizing variables to provide a low-mass panel design that carries a set of specified loadings without exceeding buckling or material strength allowables and that meets other design requirements such as upper and lower bounds on sizing variables, upper and lower bounds on overall bending, extensional and shear stiffnesses, and lower bounds on vibration frequencies. Although emphasis in PASCO is placed on flat panels having several identical bays, the only restriction on configuration modeling is that the structure is assumed to be prismatic. In addition, it is assumed that loads and temperatures do not vary along the length of a panel. Because of their wide application in aerospace structures, stiffened panels are readily handled by PASCO. The panel cross section may be composed of an arbitrary assemblage of thin, flat, rectangular plate elements that are connected together along their longitudinal edges. Each plate element consists of a balanced symmetric laminate of any number of layers of orthotropic material. Any group of element widths, layer thicknesses, and layer orientation angles may be selected as sizing variables. Substructuring is available to increase the efficiency of the analysis and to simplify the modeling of complicated structures. The Macintosh version of PASCO includes an interactive, graphic preprocessor called MacPASCO. The main objective of MacPASCO is to make the use of PASCO faster, simpler, and less error-prone. By using a graphical user interface (GUI), MacPASCO simplifies the specification of panel geometry and reduces user input errors, thus making the modeling and analysis of panel designs more efficient. The user draws the initial structural geometry on the computer screen, then uses a combination of graphic and text inputs to: refine the structural geometry, specify information required for analysis such as panel load conditions, and define design variables and constraints for minimum-mass optimization. Composite panel design is an ideal application because the graphical user interface can: serve as a visual aid, eliminate the tedious aspects of text-based input, and eliminate many sources of input errors. The current version of MacPASCO does not implement all the modeling features of PASCO, but has been found to be sufficient for many users. Many difficulties common to text-based inputs are avoided because MacPASCO uses a GUI. First, the graphic displays eliminate syntax errors, like misplaced commas and incorrect command names, because there is no text-based syntax. Second, graphic displays allow the user to see the geometry as it is created and immediately detect and correct any errors. Third, MacPASCO's drawing tools have been designed to avoid modeling errors. Fourth, the graphic displays make revisions to existing structural designs much easier and less error-prone by eliminating the need for the user to conceptualize the text input as geometry. The user can work directly with the geometry displayed on the screen. Finally, MacPASCO automatically generates the correct PASCO input from the geometry displayed on the screen. This input file can be used with any machine version of PASCO to actually perform the analysis and sizing and to output results. The DEC VAX version of PASCO is written in FORTRAN IV for batch execution and has been implemented on a DEC VAX series computer. The Macintosh version of PASCO was developed for Macintosh II series computers with at least 2Mb of RAM running MPW Pascal 3.0 and Language Systems FORTRAN 2.0 under the MPW programming environment. It includes MPW compatible makefiles for compiling the source code. The Macintosh version uses input files compatible with versions of PASCO running on different platforms. MacPASCO is written in Macintosh Programmers Workbench 3.0, MPW Pascal 3.0, and MacAPP 2.0. The Pascal source code is included on the distribution diskette. MacAPP is a development library which is not included. MacPASCO requires a Mac Plus, SE/30, or MacII, IIx, IIcx, IIci, or IIfx running System 6.0 or greater. MacPASCO is System 7.0 compatible. A minimum of 2Mb of RAM is required for execution. The Macintosh version of PASCO is distributed on four 3.5 inch 800K Macintosh format diskettes. The DEC VAX version is distributed on a 9-track 1600 BPI magnetic tape. The PASCO program was developed in 1981, adapted to the DEC VAX in 1983 and to the Macintosh in 1991. MacPASCO was released in 1992.

Description: BUCKO is a computer program developed to predict the buckling load of a rectangular compression-loaded orthotropic plate with a centrally located cutout. The plate is assumed to be a balanced, symmetric laminate of uniform thickness. The cutout shape can be elliptical, circular, rectangular, or square. The BUCKO package includes sample data that demonstrates the essence of the program and its ease of usage. BUCKO uses an approximate one-dimensional formulation of the classical two-dimensional buckling problem following the Kantorovich method. The boundary conditions are considered to be simply supported unloaded edges and either clamped or simply supported loaded edges. The plate is loaded in uniaxial compression by either uniformly displacing or uniformly stressing two opposite edges of the plate. The BUCKO analysis consists of two parts: calculation of the inplane stress distribution prior to buckling, and calculation of the plate axial load and displacement at buckling. User input includes plate planform and cutout geometry, plate membrane and bending stiffnesses, finite difference parameters, boundary condition data, and loading data. Results generated by BUCKO are the prebuckling strain energy, inplane stress resultants, buckling mode shape, critical end shortening, and average axial and transverse strains at buckling. BUCKO is written in FORTRAN V for batch execution and has been implemented on a CDC CYBER 170 series computer operating under NOS with a central memory requirement of approximately 343K of 60 bit words. This program was developed in 1984 and was last updated in 1990.

Description: A thermoviscoplastic finite element method employing the Bodner-Partom constitutve model is used to investigate the response of simplified thermal-structural models to intense local heating. The computational method formulates the problem in rate and advances the solution in time by numerical integration. The thermoviscoplastic response of simplified structures with prescribed temperatures is investigated. With rapid rises of temperature, the nickel alloy structures display initially higher yield stresses due to strain rate effects. As temperatures approach elevated values, yield stress and stiffness degrade rapidly and pronounced plastic deformation occurs.

Description: A three-bay, space, cantilever truss is probabilistically evaluated to describe progressive buckling and truss collapse in view of the numerous uncertainties associated with the structural, material, and load variables that describe the truss. Initially, the truss is deterministically analyzed for member forces, and members in which the axial force exceeds the Euler buckling load are identified. These members are then discretized with several intermediate nodes, and a probabilistic buckling analysis is performed on the truss to obtain its probabilistic buckling loads and the respective mode shapes. Furthermore, sensitivities associated with the uncertainties in the primitive variables are investigated, margin of safety values for the truss are determined, and truss end node displacements are noted. These steps are repeated by sequentially removing buckled members until onset of truss collapse is reached. Results show that this procedure yields an optimum truss configuration for a given loading and for a specified reliability.

Description: A brief mission overview of STS-49 is given, and some of the pictorially outstanding and scientifically interesting photographs obtained during the mission are presented. The Earth observations are described and include the following: the Southwestern Pacific Ocean -- wind and water; the Southwestern Pacific Ocean -- coasts and volcanoes; the US; Cuba and the Bahamas; South America; Africa; the Red Sea and Western Indian Ocean; and the Indian Subcontinent.

Description: The Science and Technology Laboratory Applications Software (ELAS) was originally designed to analyze and process digital imagery data, specifically remotely-sensed scanner data. This capability includes the processing of Landsat multispectral data; aircraft-acquired scanner data; digitized topographic data; and numerous other ancillary data, such as soil types and rainfall information, that can be stored in digitized form. ELAS has the subsequent capability to geographically reference this data to dozens of standard, as well as user created projections. As an integrated image processing system, ELAS offers the user of remotely-sensed data a wide range of capabilities in the areas of land cover analysis and general purpose image analysis. ELAS is designed for flexible use and operation and includes its own FORTRAN operating subsystem and an expandable set of FORTRAN application modules. Because all of ELAS resides in one "logical" FORTRAN program, data inputs and outputs, directives, and module switching are convenient for the user. There are over 230 modules presently available to aid the user in performing a wide range of land cover analyses and manipulation. The file management modules enable the user to allocate, define, access, and specify usage for all types of files (ELAS files, subfiles, external files etc.). Various other modules convert specific types of satellite, aircraft, and vector-polygon data into files that can be used by other ELAS modules. The user also has many module options which aid in displaying image data, such as magnification/reduction of the display; true color display; and several memory functions. Additional modules allow for the building and manipulation of polygonal areas of the image data. Finally, there are modules which allow the user to select and classify the image data. An important feature of the ELAS subsystem is that its structure allows new applications modules to be easily integrated in the future. ELAS has as a standard the flexibility to process data elements exceeding 8 bits in length, including floating point (noninteger) elements and 16 or 32 bit integers. Thus it is able to analyze and process "non-standard" nonimage data. The VAX (ERL-10017) and Concurrent (ERL-10013) versions of ELAS 9.0 are written in FORTRAN and ASSEMBLER for DEC VAX series computers running VMS and Concurrent computers running MTM. The Sun (SSC-00019), Masscomp (SSC-00020), and Silicon Graphics (SSC-00021) versions of ELAS 9.0 are written in FORTRAN 77 and C-LANGUAGE for Sun4 series computers running SunOS, Masscomp computers running UNIX, and Silicon Graphics IRIS computers running IRIX. The Concurrent version requires at least 15 bit addressing and a direct memory access channel. The VAX and Concurrent versions of ELAS both require floating-point hardware, at least 1Mb of RAM, and approximately 70Mb of disk space. Both versions also require a COMTAL display device in order to display images. For the Sun, Masscomp, and Silicon Graphics versions of ELAS, the disk storage required is approximately 115Mb, and a minimum of 8Mb of RAM is required for execution. The Sun version of ELAS requires either the X-Window System Version 11 Revision 4 or Sun OpenWindows Version 2. The Masscomp version requires a GA1000 display device and the associated "gp" library. The Silicon Graphics version requires Silicon Graphics' GL library. ELAS display functions will not work with a monochrome monitor. The standard distribution medium for the VAX version (ERL~10017) is a set of two 9-track 1600 BPI magnetic tapes in DEC VAX BACKUP format. This version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. The standard distribution medium for the Concurrent version (ERL-10013) is a set of two 9-track 1600 BPI magnetic tapes in Concurrent BACKUP format. The standard distribution medium for the Sun version (SSC-00019) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Masscomp version, (SSC-00020) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Silicon Graphics version (SSC-00021) is a .25 inch streaming magnetic IRIS tape cartridge in UNIX tar format. Version 9.0 was released in 1991. Sun4, SunOS, and Open Windows are trademarks of Sun Microsystems, Inc. MIT X Window System is licensed by Massachusetts Institute of Technology.

Description: The Science and Technology Laboratory Applications Software (ELAS) was originally designed to analyze and process digital imagery data, specifically remotely-sensed scanner data. This capability includes the processing of Landsat multispectral data; aircraft-acquired scanner data; digitized topographic data; and numerous other ancillary data, such as soil types and rainfall information, that can be stored in digitized form. ELAS has the subsequent capability to geographically reference this data to dozens of standard, as well as user created projections. As an integrated image processing system, ELAS offers the user of remotely-sensed data a wide range of capabilities in the areas of land cover analysis and general purpose image analysis. ELAS is designed for flexible use and operation and includes its own FORTRAN operating subsystem and an expandable set of FORTRAN application modules. Because all of ELAS resides in one "logical" FORTRAN program, data inputs and outputs, directives, and module switching are convenient for the user. There are over 230 modules presently available to aid the user in performing a wide range of land cover analyses and manipulation. The file management modules enable the user to allocate, define, access, and specify usage for all types of files (ELAS files, subfiles, external files etc.). Various other modules convert specific types of satellite, aircraft, and vector-polygon data into files that can be used by other ELAS modules. The user also has many module options which aid in displaying image data, such as magnification/reduction of the display; true color display; and several memory functions. Additional modules allow for the building and manipulation of polygonal areas of the image data. Finally, there are modules which allow the user to select and classify the image data. An important feature of the ELAS subsystem is that its structure allows new applications modules to be easily integrated in the future. ELAS has as a standard the flexibility to process data elements exceeding 8 bits in length, including floating point (noninteger) elements and 16 or 32 bit integers. Thus it is able to analyze and process "non-standard" nonimage data. The VAX (ERL-10017) and Concurrent (ERL-10013) versions of ELAS 9.0 are written in FORTRAN and ASSEMBLER for DEC VAX series computers running VMS and Concurrent computers running MTM. The Sun (SSC-00019), Masscomp (SSC-00020), and Silicon Graphics (SSC-00021) versions of ELAS 9.0 are written in FORTRAN 77 and C-LANGUAGE for Sun4 series computers running SunOS, Masscomp computers running UNIX, and Silicon Graphics IRIS computers running IRIX. The Concurrent version requires at least 15 bit addressing and a direct memory access channel. The VAX and Concurrent versions of ELAS both require floating-point hardware, at least 1Mb of RAM, and approximately 70Mb of disk space. Both versions also require a COMTAL display device in order to display images. For the Sun, Masscomp, and Silicon Graphics versions of ELAS, the disk storage required is approximately 115Mb, and a minimum of 8Mb of RAM is required for execution. The Sun version of ELAS requires either the X-Window System Version 11 Revision 4 or Sun OpenWindows Version 2. The Masscomp version requires a GA1000 display device and the associated "gp" library. The Silicon Graphics version requires Silicon Graphics' GL library. ELAS display functions will not work with a monochrome monitor. The standard distribution medium for the VAX version (ERL~10017) is a set of two 9-track 1600 BPI magnetic tapes in DEC VAX BACKUP format. This version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. The standard distribution medium for the Concurrent version (ERL-10013) is a set of two 9-track 1600 BPI magnetic tapes in Concurrent BACKUP format. The standard distribution medium for the Sun version (SSC-00019) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Masscomp version, (SSC-00020) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Silicon Graphics version (SSC-00021) is a .25 inch streaming magnetic IRIS tape cartridge in UNIX tar format. Version 9.0 was released in 1991. Sun4, SunOS, and Open Windows are trademarks of Sun Microsystems, Inc. MIT X Window System is licensed by Massachusetts Institute of Technology.

Description: MacMultiview is an interactive tool for the Macintosh II family which allows one to display and make computations utilizing polarimetric radar data collected by the Jet Propulsion Laboratory's imaging SAR (synthetic aperture radar) polarimeter system. The system includes the single-frequency L-band sensor mounted on the NASA CV990 aircraft and its replacement, the multi-frequency P-, L-, and C-band sensors mounted on the NASA DC-8. MacMultiview provides two basic functions: computation of synthesized polarimetric images and computation of polarization signatures. The radar data can be used to compute a variety of images. The total power image displays the sum of the polarized and unpolarized components of the backscatter for each pixel. The magnitude/phase difference image displays the HH (horizontal transmit and horizontal receive polarization) to VV (vertical transmit and vertical receive polarization) phase difference using color. Magnitude is displayed using intensity. The user may also select any transmit and receive polarization combination from which an image is synthesized. This image displays the backscatter which would have been observed had the sensor been configured using the selected transmit and receive polarizations. MacMultiview can also be used to compute polarization signatures, three dimensional plots of backscatter versus transmit and receive polarizations. The standard co-polarization signatures (transmit and receive polarizations are the same) and cross-polarization signatures (transmit and receive polarizations are orthogonal) can be plotted for any rectangular subset of pixels within a radar data set. In addition, the ratio of co- and cross-polarization signatures computed from different subsets within the same data set can also be computed. Computed images can be saved in a variety of formats: byte format (headerless format which saves the image as a string of byte values), MacMultiview (a byte image preceded by an ASCII header), and PICT2 format (standard format readable by MacMultiview and other image processing programs for the Macintosh). Images can also be printed on PostScript output devices. Polarization signatures can be saved in either a PICT format or as a text file containing PostScript commands and printed on any QuickDraw output device. The associated Stokes matrices can be stored in a text file. MacMultiview is written in C-language for Macintosh II series computers. MacMultiview will only run on Macintosh II series computers with 8-bit video displays (gray shades or color). The program also requires a minimum configuration of System 6.0, Finder 6.1, and 1Mb of RAM. MacMultiview is NOT compatible with System 7.0. It requires 32-Bit QuickDraw. Note: MacMultiview may not be fully compatible with preliminary versions of 32-Bit QuickDraw. Macintosh Programmer's Workshop and Macintosh Programmer's Workshop C (version 3.0) are required for recompiling and relinking. The standard distribution medium for this package is a set of three 800K 3.5 inch diskettes in Macintosh format. This program was developed in 1989 and updated in 1991. MacMultiview is a copyrighted work with all copyright vested in NASA. QuickDraw, Finder, Macintosh, and System 7 are trademarks of Apple Computer, Inc.

Description: The accurate partitioning of available energy into sensible and latent heat flux is crucial to the understanding of surface atmosphere interactions. This issue is more complicated in arid and semi arid regions where the relative contribution to surface fluxes from the soil and vegetation may vary significantly throughout the day and throughout the season. A three component model to estimate sensible heat flux over heterogeneous surfaces is presented. The surface was represented with two adjacent compartments. The first compartment is made up of two components, shrubs and shaded soil, the second of open 'illuminated' soil. Data collected at two different sites in Nevada (U.S.) during the Summers of 1991 and 1992 were used to evaluate model performance. The results show that the present model is sufficiently general to yield satisfactory results for both sites.

Description: Since September 1991, the ERS-1 SAR (Synthetic Aperture Radar) has collected approximately 170 frames of ocean imagery on 28 passes over the western Gulf Stream in support of ESA experiment US8-2c. SAR signatures of the north wall are seen on nearly all passes, with modulation depth varying from 3% to 35% for 100 m samples. Many small and mesoscale circulation features associated with the Stream are evident. The detailed form of the signature varies considerably, however. Narrow bright or dark linear features appear to follow streamlines, and on one occasion a strong dark line was associated with an in situ measurement of a sharp current shear. Similarly, larger spatial scale changes in backscatter over the Stream were associated with in situ measurements of atmospheric stability transitions. Physical explanations for the narrow features are not so obvious. However, the accumulation of surfactants along converging current boundaries or local short wave straining and breaking appear plausible. These preliminary results strongly suggest that a wide swath (approximately 500 km) SAR with at least 100 m resolution would be a useful adjunct to existing satellite AVHRR (Advanced Very High Resolution Radiometer) imagery.

Description: ERS-1 SAR (Synthetic Aperture Radar) data were used to specify Winter physical processes on the surface of the Greenland Sea, and SSM/I (Special Sensor Microwave/Imager) data were used to characterize the regional behavior of the ice cover. Examination of the SSM/I data indicated that the convective water was likely to be confined to small (less than 100 km) domains near the ice edge. Using ERS-1 SAR data crossing the ice edge, it was possible to identify ice edge features that are very similar to modeled plumes; they have a diameter of about 100 m, a spacing of about 300 m, and cover an area about 20 by 90 km. This plume observation is the first such identification of convection. The plumes are topped by ice, and the return-water areas are open, indicating a freshened top layer of the sea and emphasizing the importance of this layer to convection.

Description: There is renewed interest in using optical remote sensing techniques for measuring and monitoring airborne toxins. Many instruments have been used to measure traces of atmospheric molecular species. Some optical remote sensing technology is in use for measuring or monitoring toxic or pollutant gases. The examples in this article provide a number of lessons that may be of use in the new field of optical remote sensing of toxic gases. Topics covered include the following: atmospheric properties, Fourier transform spectrometry, gas-filter correlation, differential absorption laser systems, aircraft missions, laser long-path instruments, differential optical absorption spectrometry, and ozone measurement.

Description: The small-crack effect, where small fatigue cracks grow faster and at lower stress-intensity factors than large cracks, has been found to be significant for many materials and loading conditions. In this paper, plasticity effects and crack-closure modelling of small fatigue cracks are reviewed. A crack-closure model with a cyclic-plastic-zone-corrected effective stress-intensity factor range (related to the cyclic J-integral) and microstructural data on crack-initiation sites were used to calculate small-crack growth rates and fatigue lives for unnotched and notched specimens made of two aluminum alloys. The crack-closure transient from the plastic wake was shown to be the dominant cause of the small-crack effect and plasticity effects on the cyclic-plastic-zone-corrected stress-intensity factor range were negligible except at extremely high stress levels. Small-crack growth rates and fatigue lives under both constant-amplitude and spectrum loading from tests and analyses agreed well.

Description: A study is made of the thermomechanical buckling of flat unstiffened composite panels with central circular cutouts. The panels are subjected to combined temperature changes and applied edge loading (or edge displacements). The analysis is based on a first-order shear deformation plate theory. A mixed formulation is used with the fundamental unknowns consisting of the generalized displacements and the stress resultants of the plate. Both the stability boundary and the sensitivity coefficients are evaluated. The sensitivity coefficients measure the sensitivity of the buckling response to variations in the different lamination and material parameters of the panel. Numerical results are presented showing the effects of the variations in the hole diameter, laminate stacking sequence, fiber orientation, and aspect ratio of the panel on the thermomechanical buckling response and its sensitivity coefficients.

Description: Global and regional temperature variations in the lower troposphere and lower stratosphere are examined for the period 1979-92 from Microwave Sounder Unit (MSU) data obtained by the Television Infrared Observation Satellite (TIROS)-N series of National Oceanic and Atmospheric Administration (NOAA) operational satellites. In the lower troposphere, globally-averaged temperature variations appear to be dominated by tropical El Nino (warm) and La Nina (cool) events and volcanic eruptions. The Pinatubo volcanic eruption in June 1991 appears to have initiated a cooling trend which persisted through the most recent data analyzed (July, 1992), and largely overwhelmed the warming from the 1991-92 El Nino. The cooling has been stronger in the Northern Hemisphere than in the Southern Hemisphere. The temperature trend over the 13.5 year satellite record is small (+0.03 C) compared to the year-to-year variability (0.2-0.4 C), making detection of any global warming signal fruitless to date. However, the future global warming trend, currently predicted to be around 0.3 C/decade, will be much easier to discern should it develop. The lower stratospheric temperature record is dominated by warm episodes from the Pinatubo eruption and the March 1982 eruption of El Chichon volcano.

Description: PATSTAGS translates PATRAN finite model data into STAGS (Structural Analysis of General Shells) input records to be used for engineering analysis. The program reads data from a PATRAN neutral file and writes STAGS input records into a STAGS input file and a UPRESS data file. It is able to support translations of nodal constraints, nodal, element, force and pressure data. PATSTAGS uses three files: the PATRAN neutral file to be translated, a STAGS input file and a STAGS pressure data file. The user provides the names for the neutral file and the desired names of the STAGS files to be created. The pressure data file contains the element live pressure data used in the STAGS subroutine UPRESS. PATSTAGS is written in FORTRAN 77 for DEC VAX series computers running VMS. The main memory requirement for execution is approximately 790K of virtual memory. Output blocks can be modified to output the data in any format desired, allowing the program to be used to translate model data to analysis codes other than STAGSC-1 (HQN-10967). This program is available in DEC VAX BACKUP format on a 9-track magnetic tape or TK50 tape cartridge. Documentation is included in the price of the program. PATSTAGS was developed in 1990. DEC, VAX, TK50 and VMS are trademarks of Digital Equipment Corporation.

Description: The Composite Structure Preliminary Sizing program, COMPSIZE, is an analytical tool which structural designers can use when doing approximate stress analysis to select or verify preliminary sizing choices for composite structural members. It is useful in the beginning stages of design concept definition, when it is helpful to have quick and convenient approximate stress analysis tools available so that a wide variety of structural configurations can be sketched out and checked for feasibility. At this stage of the design process the stress/strain analysis does not need to be particularly accurate because any configurations tentatively defined as feasible will later be analyzed in detail by stress analysis specialists. The emphasis is on fast, user-friendly methods so that rough but technically sound evaluation of a broad variety of conceptual designs can be accomplished. Analysis equations used are, in most cases, widely known basic structural analysis methods. All the equations used in this program assume elastic deformation only. The default material selection is intermediate strength graphite/epoxy laid up in a quasi-isotropic laminate. A general flat laminate analysis subroutine is included for analyzing arbitrary laminates. However, COMPSIZE should be sufficient for most users to presume a quasi-isotropic layup and use the familiar basic structural analysis methods for isotropic materials, after estimating an appropriate elastic modulus. Homogeneous materials can be analyzed as simplified cases. The COMPSIZE program is written in IBM BASICA. The program format is interactive. It was designed on an IBM Personal Computer operating under DOS with a central memory requirement of approximately 128K. It has been implemented on an IBM compatible with GW-BASIC under DOS 3.2. COMPSIZE was developed in 1985.

Description: Several types of algorithms are generally used to process digital imagery such as Landsat data. The most commonly used algorithms perform the task of registration, compression, and classification. Because there are different techniques available for performing registration, compression, and classification, imagery data users need a rationale for selecting a particular approach to meet their particular needs. This collection of registration, compression, and classification algorithms was developed so that different approaches could be evaluated and the best approach for a particular application determined. Routines are included for six registration algorithms, six compression algorithms, and two classification algorithms. The package also includes routines for evaluating the effects of processing on the image data. This collection of routines should be useful to anyone using or developing image processing software. Registration of image data involves the geometrical alteration of the imagery. Registration routines available in the evaluation package include image magnification, mapping functions, partitioning, map overlay, and data interpolation. The compression of image data involves reducing the volume of data needed for a given image. Compression routines available in the package include adaptive differential pulse code modulation, two-dimensional transforms, clustering, vector reduction, and picture segmentation. Classification of image data involves analyzing the uncompressed or compressed image data to produce inventories and maps of areas of similar spectral properties within a scene. The classification routines available include a sequential linear technique and a maximum likelihood technique. The choice of the appropriate evaluation criteria is quite important in evaluating the image processing functions. The user is therefore given a choice of evaluation criteria with which to investigate the available image processing functions. All of the available evaluation criteria basically compare the observed results with the expected results. For the image reconstruction processes of registration and compression, the expected results are usually the original data or some selected characteristics of the original data. For classification processes the expected result is the ground truth of the scene. Thus, the comparison process consists of determining what changes occur in processing, where the changes occur, how much change occurs, and the amplitude of the change. The package includes evaluation routines for performing such comparisons as average uncertainty, average information transfer, chi-square statistics, multidimensional histograms, and computation of contingency matrices. This collection of routines is written in FORTRAN IV for batch execution and has been implemented on an IBM 360 computer with a central memory requirement of approximately 662K of 8 bit bytes. This collection of image processing and evaluation routines was developed in 1979.

Description: The Science and Technology Laboratory Applications Software (ELAS) was originally designed to analyze and process digital imagery data, specifically remotely-sensed scanner data. This capability includes the processing of Landsat multispectral data; aircraft-acquired scanner data; digitized topographic data; and numerous other ancillary data, such as soil types and rainfall information, that can be stored in digitized form. ELAS has the subsequent capability to geographically reference this data to dozens of standard, as well as user created projections. As an integrated image processing system, ELAS offers the user of remotely-sensed data a wide range of capabilities in the areas of land cover analysis and general purpose image analysis. ELAS is designed for flexible use and operation and includes its own FORTRAN operating subsystem and an expandable set of FORTRAN application modules. Because all of ELAS resides in one "logical" FORTRAN program, data inputs and outputs, directives, and module switching are convenient for the user. There are over 230 modules presently available to aid the user in performing a wide range of land cover analyses and manipulation. The file management modules enable the user to allocate, define, access, and specify usage for all types of files (ELAS files, subfiles, external files etc.). Various other modules convert specific types of satellite, aircraft, and vector-polygon data into files that can be used by other ELAS modules. The user also has many module options which aid in displaying image data, such as magnification/reduction of the display; true color display; and several memory functions. Additional modules allow for the building and manipulation of polygonal areas of the image data. Finally, there are modules which allow the user to select and classify the image data. An important feature of the ELAS subsystem is that its structure allows new applications modules to be easily integrated in the future. ELAS has as a standard the flexibility to process data elements exceeding 8 bits in length, including floating point (noninteger) elements and 16 or 32 bit integers. Thus it is able to analyze and process "non-standard" nonimage data. The VAX (ERL-10017) and Concurrent (ERL-10013) versions of ELAS 9.0 are written in FORTRAN and ASSEMBLER for DEC VAX series computers running VMS and Concurrent computers running MTM. The Sun (SSC-00019), Masscomp (SSC-00020), and Silicon Graphics (SSC-00021) versions of ELAS 9.0 are written in FORTRAN 77 and C-LANGUAGE for Sun4 series computers running SunOS, Masscomp computers running UNIX, and Silicon Graphics IRIS computers running IRIX. The Concurrent version requires at least 15 bit addressing and a direct memory access channel. The VAX and Concurrent versions of ELAS both require floating-point hardware, at least 1Mb of RAM, and approximately 70Mb of disk space. Both versions also require a COMTAL display device in order to display images. For the Sun, Masscomp, and Silicon Graphics versions of ELAS, the disk storage required is approximately 115Mb, and a minimum of 8Mb of RAM is required for execution. The Sun version of ELAS requires either the X-Window System Version 11 Revision 4 or Sun OpenWindows Version 2. The Masscomp version requires a GA1000 display device and the associated "gp" library. The Silicon Graphics version requires Silicon Graphics' GL library. ELAS display functions will not work with a monochrome monitor. The standard distribution medium for the VAX version (ERL~10017) is a set of two 9-track 1600 BPI magnetic tapes in DEC VAX BACKUP format. This version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. The standard distribution medium for the Concurrent version (ERL-10013) is a set of two 9-track 1600 BPI magnetic tapes in Concurrent BACKUP format. The standard distribution medium for the Sun version (SSC-00019) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Masscomp version, (SSC-00020) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Silicon Graphics version (SSC-00021) is a .25 inch streaming magnetic IRIS tape cartridge in UNIX tar format. Version 9.0 was released in 1991. Sun4, SunOS, and Open Windows are trademarks of Sun Microsystems, Inc. MIT X Window System is licensed by Massachusetts Institute of Technology.

Description: Calibration of polarimetric radar systems is a field of research in which great progress has been made over the last few years. POLCAL (Polarimetric Radar Calibration) is a software tool intended to assist in the calibration of Synthetic Aperture Radar (SAR) systems. In particular, POLCAL calibrates Stokes matrix format data produced as the standard product by the NASA/Jet Propulsion Laboratory (JPL) airborne imaging synthetic aperture radar (AIRSAR). POLCAL was designed to be used in conjunction with data collected by the NASA/JPL AIRSAR system. AIRSAR is a multifrequency (6 cm, 24 cm, and 68 cm wavelength), fully polarimetric SAR system which produces 12 x 12 km imagery at 10 m resolution. AIRSTAR was designed as a testbed for NASA's Spaceborne Imaging Radar program. While the images produced after 1991 are thought to be calibrated (phase calibrated, cross-talk removed, channel imbalance removed, and absolutely calibrated), POLCAL can and should still be used to check the accuracy of the calibration and to correct it if necessary. Version 4.0 of POLCAL is an upgrade of POLCAL version 2.0 released to AIRSAR investigators in June, 1990. New options in version 4.0 include automatic absolute calibration of 89/90 data, distributed target analysis, calibration of nearby scenes with calibration parameters from a scene with corner reflectors, altitude or roll angle corrections, and calibration of errors introduced by known topography. Many sources of error can lead to false conclusions about the nature of scatterers on the surface. Errors in the phase relationship between polarization channels result in incorrect synthesis of polarization states. Cross-talk, caused by imperfections in the radar antenna itself, can also lead to error. POLCAL reduces cross-talk and corrects phase calibration without the use of ground calibration equipment. Removing the antenna patterns during SAR processing also forms a very important part of the calibration of SAR data. Errors in the processing altitude or in the aircraft roll angle are possible causes of error in computing the antenna patterns inside the processor. POLCAL uses an altitude error correction algorithm to correctly remove the antenna pattern from the SAR images. POLCAL also uses a topographic calibration algorithm to reduce calibration errors resulting from ground topography. By utilizing the backscatter measurements from either the corner reflectors or a well-known distributed target, POLCAL can correct the residual amplitude offsets in the various polarization channels and correct for the absolute gain of the radar system. POLCAL also gives the user the option of calibrating a scene using the calibration data from a nearby site. This allows precise calibration of all the scenes acquired on a flight line where corner reflectors were present. Construction and positioning of corner reflectors is covered extensively in the program documentation. In an effort to keep the POLCAL code as transportable as possible, the authors eliminated all interactions with a graphics display system. For this reason, it is assumed that users will have their own software for doing the following: (1) synthesize an image using HH or VV polarization, (2) display the synthesized image on any display device, and (3) read the pixel locations of the corner reflectors from the image. The only inputs used by the software (in addition to the input Stokes matrix data file) is a small data file with the corner reflector information. POLCAL is written in FORTRAN 77 for use on Sun series computers running SunOS and DEC VAX computers running VMS. It requires 4Mb of RAM under SunOS and 3.7Mb of RAM under VMS for execution. The standard distribution medium for POLCAL is a .25 inch streaming magnetic tape cartridge in UNIX tar format. It is also available on a 9-track 1600 BPI magnetic tape in DEC VAX FILES-11 format or on a TK50 tape cartridge in DEC VAX FILES-11 format. Other distribution media may be available upon request. Documentation is included in the price of the program. POLCAL 4.0 was released in 1992 and is a copyrighted work with all copyright vested in NASA.

Description: This paper introduces a study on an Electromagnetically Levitated Vibration Isolation System (ELVIS) for isolation control of large-scale vibration. This system features no mechanical contact between the isolation table and the installation floor, using a total of four electromagnetic actuators which generate magnetic levitation force in the vertical and horizontal directions. The configuration of the magnet for the vertical direction is designed to prevent any generation of restoring vibratory force in the horizontal direction. The isolation system is set so that vibration control effects due to small earthquakes can be regulated to below 5(gal) versus horizontal vibration levels of the installation floor of up t 25(gal), and those in the horizontal relative displacement of up to 30 (mm) between the floor and levitated isolation table. In particular, studies on the relative displacement between the installation floor and the levitated isolation table have been made for vibration control in the horizontal direction. In case of small-scale earthquakes (Taft wave scaled: max. 25 gal), the present system has been confirmed to achieve a vibration isolation to a level below 5 gal. The vibration transmission ratio of below 1/10 has been achieved versus continuous micro-vibration (approx. one gal) in the horizontal direction on the installation floor.

Description: The objective of this study is to experimentally determine an empirical model of the vibrational dynamics of the Spacecraft COntrol Laboratory Experiment (SCOLE) facility. The first two flexible modes of this test article are identified using a linear least-square identification procedure and the data utilized for this procedure are obtained by exciting the structure from a quiescent state with torque wheels. The time history data of rate gyro sensors and accelerometers due to excitation and after excitation in terms of free-decay are used in the parameter estimation of the vibrational model. The free-decay portion of the data is analyzed using the Discrete Fourier transform to determine the optimal model order to use in modelling the response. Linear least-square analysis is then used to select the parameters that best fit the output of an Autoregressive (AR) model to the data. The control effectiveness of the torque wheels is then determined using the excitation portion of the test data, again using linear least squares.

Description: A new approach to find homogeneous models for beam-like repeated flexible structures is proposed which conceptually involves two steps. The first step involves the approximation of 3-D non-homogeneous model by a 1-D periodic beam model. The structure is modeled as a 3-D non-homogeneous continuum. The displacement field is approximated by Taylor series expansion. Then, the cross sectional mass and stiffness matrices are obtained by energy equivalence using their additive properties. Due to the repeated nature of the flexible bodies, the mass, and stiffness matrices are also periodic. This procedure is systematic and requires less dynamics detail. The first step involves the homogenization from a 1-D periodic beam model to a 1-D homogeneous beam model. The periodic beam model is homogenized into an equivalent homogeneous beam model using the additive property of compliance along the generic axis. The major departure from previous approaches in literature is using compliance instead of stiffness in homogenization. An obvious justification is that the stiffness is additive at each cross section but not along the generic axis. The homogenized model preserves many properties of the original periodic model.

Description: This report will discuss the design of a liquid hydrogen fuel tank constructed from composite materials. The focus of this report is to recommend a design for a fuel tank which will be able to withstand all static and dynamic forces during manned flight. Areas of study for the design include material selection, material structural analysis, heat transfer, thermal expansion, and liquid hydrogen diffusion. A structural analysis FORTRAN program was developed for analyzing the buckling and yield characteristics of the tank. A thermal analysis Excel spreadsheet was created to determine a specific material thickness which will minimize heat transfer through the wall of the tank. The total mass of the tank was determined by the combination of both structural and thermal analyses. The report concludes with the recommendation of a layered material tank construction. The designed system will include exterior insulation, combination of metal and organize composite matrices and honeycomb.

Description: Remote terrestrial sensing (RTS) data is constantly being collected from a variety of space-based and earth-based sensors. The collected data, and especially 'value-added' analyses of the data, are finding growing application for commercial, government, and scientific purposes. The scale of this data collection and analysis is truly enormous; e.g., by 1995, the amount of data available in just one sector, NASA space science, will reach 5 petabytes. Moreover, the amount of data, and the value of analyzing the data, are expected to increase dramatically as new satellites and sensors become available (e.g., NASA's Earth Observing System satellites). Lockheed and other companies are beginning to provide data and analysis commercially. A critical issue for the exploitation of collected data is the dissemination of data and value-added analyses to a diverse and widely distributed customer base. Customers must be able to use their computational environment (eventually the National Information Infrastructure) to obtain timely and complete information, without having to know the details of where the relevant data resides and how it is accessed. Customers must be able to routinely use standard, widely available (and, therefore, low cost) analyses, while also being able to readily create on demand highly customized analyses to make crucial decisions. The diversity of user needs creates a difficult software problem: how can users easily state their needs, while the computational environment assumes the responsibility of finding (or creating) relevant information, and then delivering the results in a form that users understand? A software agent is a self-contained, active software module that contains an explicit representation of its operational knowledge. This explicit representation allows agents to examine their own capabilities in order to modify their goals to meet changing needs and to take advantage of dynamic opportunities. In addition, the explicit representation allows agents to advertize their capabilities and results to other agents, thereby allowing the collection of agents to reuse each others work.

Description: The purpose of this project was to develop a simple model which could be used to study residual stress. The mechanism that results in residual stresses in the welding process starts with the deposition of molten weld metal which heats the immediately adjacent material. After solidification of weld material, normal thermal shrinkage is resisted by the adjacent, cooler material. When the thermal strain exceeds the elastic strain corresponding to the yield point stress, the stress level is limited by this value, which decreases with increasing temperature. Cooling then causes elastic unloading which is restrained by the adjoining material. Permanent plastic strain occurs, and tension is caused in the region immediately adjacent to the weld material. Compression arises in the metal farther from the weld in order to maintain overall static equilibrium. Subsequent repair welds may add to the level of residual stresses. The level of residual stress is related to the onset of fracture during welding. Thus, it is of great importance to be able to predict the level of residual stresses remaining after a weld procedure, and to determine the factors, such as weld speed, temperature, direction, and number of passes, which may affect the magnitude of remaining residual stress. It was hoped to use traditional analytical modeling techniques so that it would be easier to comprehend the effect of these variables on the resulting stress. This approach was chosen in place of finite element methods so as to facilitate the understanding of the physical processes. The accuracy of the results was checked with some existing experimental studies giving residual stress levels found from x-ray diffraction measurements.

Description: Weather forecasting has been called the second oldest profession. To do so accurately and with some consistency requires an ability to understand the processes which create the clouds, drive the winds, and produce the ever changing atmospheric conditions. Measurement of basic parameters such as temperature, water vapor content, pressure, windspeed and wind direction throughout the three dimensional atmosphere form the foundation upon which a modern forecast is created. Doppler radar, and space borne remote sensing have provided forecasters the new tools with which to ply their trade.

Description: All long-duration spacecraft in low-earth-orbit are subject to high speed impacts by meteoroids and pieces of orbital debris. The threat of damage from such impacts is a significant design consideration in the development of long duration earth-orbiting spacecraft. This report presents the results of a study whose objective was to develop an empirical model to predict the magnitude of the various cracking and through-hole creation phenomena accompanying a habitable module penetration. The significance of the work performed is that the model predictions can be fed directly into a survivability analysis to determine whether or not module unzipping would occur under a specific set of impact conditions. The likelihood of module unzipping over a structure's lifetime can also be determined in such an analysis. In addition, effective hole size predictions can be used as part of a survivability analysis to determine the time available for module evacuation prior to the onset of incapacitation due to air loss. Some of the phenomena considered include maximum petal length, maximum tip-to-tip crack distance, depth of petal deformation, number of cracks formed, orientation of the maximum tip-to-tip distance with respect to the inner wall grain direction, and the effective inner wall hole diameter.

Description: It is common practice to use split sleeve coldworking of fastener holes as a means of extending the fatigue life of metal structures. In search of lower manufacturing costs, the aerospace industry is examining the split mandrel (sleeveless) coldworking process as an alternative method of coldworking fastener holes in metal structures. The split mandrel process (SpM) significantly extends the fatigue life of metal structures through the introduction of a residual compressive stress in a manner that is very similar to the split sleeve system (SpSl). Since the split mandrel process is significantly less expensive than the split sleeve process and more adaptable to robotic automation, it will have a notable influence upon other new manufacture of metal structures which require coldworking a significant number of holes, provided the aerospace community recognizes that the resulting residual stress distributions and fatigue life improvement are the same for both processes. Considerable testing has validated the correctness of that conclusion. The findings presented in this paper represent the results of an extensive research and development program, comprising data collected from over 400 specimens fabricated from 2024-T3 and 7075-T651 aluminum alloys in varied configurations, which quantify the benefits (fatigue enhancement and cost savings) of automating a sleeveless coldworking system.

Description: R-curves were predicted for Alclad 2024-T3 and C188-T3 sheet using the results of small-coupon Kahn tear tests in combination with two-dimensional elastic-plastic finite element stress analyses. The predictions were compared to experimental R-curves from 6.3, 16 and 60-inch wide M(T) specimens and good agreement was obtained. The method is an inexpensive alternative to wide panel testing for characterizing the fracture toughness of damage-tolerant sheet alloys. The usefulness of this approach was demonstrated by performing residual strength calculations for a two-bay crack in a representative fuselage structure. C188-T3 was predicted to have a 24 percent higher load carrying capability than 2024-T3 in this application as a result of its superior fracture toughness.

Description: Cold working holes for improved fatigue life of fastener holes are widely used on aircraft. This paper presents methods used by the authors to determine the percent of cold working to be applied and to analyze fatigue crack growth of cold worked fastener holes. An elastic, perfectly-plastic analysis of a thick-walled tube is used to determine the stress field during the cold working process and the residual stress field after the process is completed. The results of the elastic/plastic analysis are used to determine the amount of cold working to apply to a hole. The residual stress field is then used to perform damage tolerance analysis of a crack growing out of a cold worked fastener hole. This analysis method is easily implemented in existing crack growth computer codes so that the cold worked holes can be used to extend the structural life of aircraft. Analytical results are compared to test data where appropriate.

Description: Effective Total Fatigue Life and Crack Growth Scatter Models are proposed. The first of them is based on the power form of the Wohler curve, fatigue scatter dependence on mean life value, cycle stress ratio influence on fatigue scatter, and validated description of the mean stress influence on the mean fatigue life. The second uses in addition are fracture mechanics approach, assumption of initial damage existence, and Paris equation. Simple formulas are derived for configurations of models. A preliminary identification of the parameters of the models is fulfilled on the basis of experimental data. Some new and important results for fatigue and crack growth scatter characteristics are obtained.

Description: Terminating action is a remedial repair which entails the replacement of shear head countersunk rivets with universal head rivets which have a larger shank diameter. The procedure was developed to eliminate the risk of widespread fatigue damage (WFD) in the upper rivet row of a fuselage lap joint. A test and evaluation program has been conducted by Foster-Miller, Inc. (FMI) to evaluate the terminating action repair of the upper rivet row of a commercial aircraft fuselage lap splice. Two full scale fatigue tests were conducted on fuselage panels using the growth of fatigue cracks in the lap joint. The second test was performed to evaluate the effectiveness of the terminating action repair. In both tests, cyclic pressurization loading was applied to the panels while crack propagation was recorded at all rivet locations at regular intervals to generate detailed data on conditions of fatigue crack initiation, ligament link-up, and fuselage fracture. This program demonstrated that the terminating action repair substantially increases the fatigue life of a fuselage panel structure and effectively eliminates the occurrence of cracking in the upper rivet row of the lap joint. While high cycle crack growth was recorded in the middle rivet row during the second test, failure was not imminent when the test was terminated after cycling to well beyond the service life. The program also demonstrated that the initiation, propagation, and linkup of WFD in full-scale fuselage structures can be simulated and quantitatively studied in the laboratory. This paper presents an overview of the testing program and provides a detailed discussion of the data analysis and results. Crack distribution and propagation rates and directions as well as frequency of cracking are presented for both tests. The progression of damage to linkup of adjacent cracks and to eventual overall panel failure is discussed. In addition, an assessment of the effectiveness of the terminating action repair and the occurrence of cracking in the middle rivet row is provided, and conclusions of practical interest are drawn.

Description: The characteristics of widespread fatigue damage (WSFD) in fuselage riveted structure were established by detailed nondestructive and destructive examinations of fatigue damage contained in a full size fuselage test article. The objectives of this work were to establish an experimental data base for validating emerging WSFD analytical prediction methodology and to identify first order effects that contribute to fatigue crack initiation and growth. Detailed examinations were performed on a test panel containing four bays of a riveted lap splice joint. The panel was removed from a full scale fuselage test article after receiving 60,000 full pressurization cycles. The results of in situ examinations document the progression of fuselage skin fatigue crack growth through crack linkup. Detailed tear down examinations and fractography of the lap splice joint region revealed fatigue crack initiation sites, crack morphology and crack linkup geometry. From this large data base, distributions of crack size and locations are presented and discussions of operative damage mechanisms are offered.

Description: The CORPUS (Computation Of Retarded Propagation Under Spectrum loading) crack growth prediction model for variable-amplitude loading, as introduced by De Koning, was based on crack closure. It includes a multiple-overload effect and a transition from plane strain to plane stress. In the modified CORPUS model an underload affected zone (ULZ) is introduced, which is significant for flight-simulation loading in view of the once per flight compressive ground load. The ULZ is associated with reversed plastic deformation induced by the underloads after crack closure has already occurred. Predictions of the crack growth fatigue life are presented for a large variety of flight-simulation test series on 2024-T3 sheet specimens in order to reveal the effects of a number of variables: the design stress level, the gust spectrum severity, the truncation level (clipping), omission of small cycles, and the ground stress level. Tests with different load sequences are also included. The trends of the effects induced by the variables are correctly predicted. The quantitative agreement between the predictions and the test results is also satisfactory.

Description: The minimum weight optimization of structural systems, subject to strength and displacement constraints as well as size side constraints, was investigated by the Simultaneous ANalysis and Design (SAND) approach. As an optimizer, the code NPSOL was used which is based on a sequential quadratic programming (SQP) algorithm. The structures were modeled by the finite element method. The finite element related input to NPSOL was automatically generated from the input decks of such standard FEM/optimization codes as NASTRAN or ASTROS, with the stiffness matrices, at present, extracted from the FEM code ANALYZE. In order to avoid ill-conditioned matrices that can be encountered when the global stiffness equations are used as additional nonlinear equality constraints in the SAND approach (with the displacements as additional variables), the matrix displacement method was applied. In this approach, the element stiffness equations are used as constraints instead of the global stiffness equations, in conjunction with the nodal force equilibrium equations. This approach adds the element forces as variables to the system. Since, for complex structures and the associated large and very sparce matrices, the execution times of the optimization code became excessive due to the large number of required constraint gradient evaluations, the Kreisselmeier-Steinhauser function approach was used to decrease the computational effort by reducing the nonlinear equality constraint system to essentially a single combined constraint equation. As the linear equality and inequality constraints require much less computational effort to evaluate, they were kept in their previous form to limit the complexity of the KS function evaluation. To date, the standard three-bar, ten-bar, and 72-bar trusses have been tested. For the standard SAND approach, correct results were obtained for all three trusses although convergence became slower for the 72-bar truss. When the matrix displacement method was used, correct results were still obtained, but the execution times became excessive due to the large number of constraint gradient evaluations required. Using the KS function, the computational effort dropped, but the optimization seemed to become less robust. The investigation of this phenomenon is continuing. As an alternate approach, the code MINOS for the optimization of sparse matrices can be applied to the problem in lieu of the Kreisselmeier-Steinhauser function. This investigation is underway.

Description: Computational Fluid Dynamics (or CFD) methods are very familiar to the research community. Even the general public has had some exposure to CFD images, primarily through the news media. However, very little attention has been paid to CST--Computational Structures Technology. Yet, no important design can be completed without it. During the first half of this century, researchers only dreamed of designing and building structures on a computer. Today their dreams have become practical realities as computational methods are used in all phases of design, fabrication and testing of engineering systems. Increasingly complex structures can now be built in even shorter periods of time. Over the past four decades, computer technology has been developing, and early finite element methods have grown from small in-house programs to numerous commercial software programs. When coupled with advanced computing systems, they help engineers make dramatic leaps in designing and testing concepts. The goals of CST include: predicting how a structure will behave under actual operating conditions; designing and complementing other experiments conducted on a structure; investigating microstructural damage or chaotic, unpredictable behavior; helping material developers in improving material systems; and being a useful tool in design systems optimization and sensitivity techniques. Applying CST to a structure problem requires five steps: (1) observe the specific problem; (2) develop a computational model for numerical simulation; (3) develop and assemble software and hardware for running the codes; (4) post-process and interpret the results; and (5) use the model to analyze and design the actual structure. Researchers in both industry and academia continue to make significant contributions to advance this technology with improvements in software, collaborative computing environments and supercomputing systems. As these environments and systems evolve, computational structures technology will evolve. By using CST in the design and operation of future structures systems, engineers will have a better understanding of how a system responds and lasts, more cost-effective methods of designing and testing models, and improved productivity. For informational and educational purposes, a videotape is being produced using both static and dynamic images from research institutions, software and hardware companies, private individuals, and historical photographs and drawings. The extensive number of CST resources indicates its widespread use. Applications run the gamut from simpler university-simulated problems to those requiring solutions on supercomputers. In some cases, an image or an animation will be mapped onto the actual structure to show the relevance of the computer model to the structure. Transferring the digital files to videotape presents a number of problems related to maintaining the quality of the original image, while still producing a broadcast quality videotape. Since researchers normally do not create a computer image using traditional composition theories or video production requirements, often the image loses some of its original digital quality and impact when transferred to videotape. Although many CST images are currently available, those that are edited into the final project must meet two important criteria: they must complement the narration, and they must be broadcast quality when recorded on videotape.

Description: Thermal stress analyses are an important aspect in the development of aerospace vehicles at NASA-LaRC. These analyses require knowledge of the temperature distributions within the vehicle structures which consequently necessitates the need for accurate thermal property data. The overall goal of this ongoing research effort is to develop methodologies for the estimation of the thermal property data needed to describe the temperature responses of these complex structures. The research strategy undertaken utilizes a building block approach. The idea here is to first focus on the development of property estimation methodologies for relatively simple conditions, such as isotropic materials at constant temperatures, and then systematically modify the technique for the analysis of more and more complex systems, such as anisotropic multi-component systems. The estimation methodology utilized is a statistically based method which incorporates experimental data and a mathematical model of the system. Several aspects of this overall research effort were investigated during the time of the ASEE summer program. One important aspect involved the calibration of the estimation procedure for the estimation of the thermal properties through the thickness of a standard material. Transient experiments were conducted using a Pyrex standard at various temperatures, and then the thermal properties (thermal conductivity and volumetric heat capacity) were estimated at each temperature. Confidence regions for the estimated values were also determined. These results were then compared to documented values. Another set of experimental tests were conducted on carbon composite samples at different temperatures. Again, the thermal properties were estimated for each temperature, and the results were compared with values obtained using another technique. In both sets of experiments, a 10-15 percent off-set between the estimated values and the previously determined values was found. Another effort was related to the development of the experimental techniques. Initial experiments required a resistance heater placed between two samples. The design was modified such that the heater was placed on the surface of only one sample, as would be necessary in the analysis of built up structures. Experiments using the modified technique were conducted on the composite sample used previously at different temperatures. The results were within 5 percent of those found using two samples. Finally, an initial heat transfer analysis, including conduction, convection and radiation components, was completed on a titanium sandwich structural sample. Experiments utilizing this sample are currently being designed and will be used to first estimate the material's effective thermal conductivity and later to determine the properties associated with each individual heat transfer component.

Description: The previous design philosophies involving safe life, fail-safe and damage tolerance concepts become inadequate for assuring the safety of aging aircraft structures. For example, the failure mechanism for the Aloha Airline accident involved the coalescence of undetected small cracks at the rivet holes causing a section of the fuselage to peel open during flight. Therefore, the fuselage structure should be designed to have sufficient residual strength under worst case crack configurations and in-flight load conditions. Residual strength is interpreted as the maximum load carrying capacity prior to unstable crack growth. Internal pressure and bending moment constitute the two major components of the external loads on the fuselage section during flight. Although the stiffeners in the form of stringers, frames and tear straps sustain part of the external loads, the significant portion of the load is taken up by the skin. In the presence of a large crack in the skin, the crack lips bulge out with considerable yielding; thus, the geometric and material nonlinearities must be included in the analysis for predicting residual strength. Also, these nonlinearities do not permit the decoupling of in-plane and out-of-plane bending deformations. The failure criterion combining the concepts of absorbed specific energy and strain energy density addresses the aforementioned concerns. The critical absorbed specific energy (local toughness) for the material is determined from the global specimen response and deformation geometry based on the uniaxial tensile test data and detailed finite element modeling of the specimen response. The use of the local toughness and stress-strain response at the continuum level eliminates the size effect. With this critical parameter and stress-strain response, the finite element analysis of the component by using STAGS along with the application of this failure criterion provides the stable crack growth calculations for residual strength predictions.

Description: A sandwich construction consists of a low-density core material with high strength face sheets bounded to the top and bottom surfaces. The construction has been widely used in the aerospace and marine industries due to its outstanding characteristics such as noise absorption, weight minimization, heat insulation, and better bending stiffness. In sandwich structures used in high-performance aircraft, the face sheets are often made of fiber-reinforced composite materials and the core is made of honeycomb. The structures may also have variable thickness so as to satisfy aerodynamic requirements. In the stress analysis, the constant-thickness face sheets are usually considered as membrane and the core is assumed to be inextensible but deformable in the thickness direction. The static behavior of variable-thickness, isotropic and homogeneous sandwich beams was successfully studied by employing a constant-thickness theory but allowing stiffnesses to vary in accordance with local thickness variations. It has been recently found in a refined theory that the analyses based on the constant thickness theory locally can lead to significant errors in structural responses if the sandwich beam is thickness-tapered and the cores are deformable in transverse shear. The errors arise mainly from two factors: (1) the transverse shear components of the membrane forces in the face sheets alter the transverse shears carried by the core; and (2) the face-sheet membrane strains arise from transverse shear deformation of the core. In practice the variable thickness may not only exist in core but also in face sheets. The thickness-variations may even be a type of step function. In this case the transverse shear stress in the face sheets and bending stress in the core should be taken into account in the refined theory mentioned. In the present study, energy principles are employed in deriving governing equations for general bending of anisotropic sandwich beams with variable thickness in both face sheets and cores. Solutions to these equations are based on a finite difference scheme. As an example in application, a simply supported thickness-tapered sandwich beam subject to a concentrated load at its center is considered. Let W' be the maximum deflection of the beam in which face sheets are considered as membrane, while W'' is that based on using the modified refined theory. It is found that W' is always larger than W'', however, the magnitude of (W'- W'') appears to be insensitive to the change of the taper of the beam.

Description: Algorithms and techniques for use in the identification and location of large buildings in digitized copies of aerial photographs are developed and tested. The building data would be used in the simulation of objects located in the vicinity of an airport that may be detected by aircraft radar. Two distinct approaches are considered. Most building footprints are rectangular in form. The first approach studied is to search for right-angled corners that characterize rectangular objects and then to connect these corners to complete the building. This problem is difficult because many nonbuilding objects, such as street corners, parking lots, and ballparks often have well defined corners which are often difficult to distinguish from rooftops. Furthermore, rooftops come in a number of shapes, sizes, shadings, and textures which also limit the discrimination task. The strategy used linear sequences of different samples to detect straight edge segments at multiple angles and to determine when these segments meet at approximately right-angles with respect to each other. This technique is effective in locating corners. The test image used has a fairly rectangular block pattern oriented about thirty degrees clockwise from a vertical alignment, and the overall measurement data reflect this. However, this technique does not discriminate between buildings and other objects at an operationally suitable rate. In addition, since multiple paths are tested for each image pixel, this is a time consuming task. The process can be speeded up by preprocessing the image to locate the more optimal sampling paths. The second approach is to rely on a human operator to identify and select the building objects and then to have the computer determine the outline and location of the selected structures. When presented with a copy of a digitized aerial photograph, the operator uses a mouse and cursor to select a target building. After a button on the mouse is pressed, with the cursor fully within the perimeter of the building, the program scans from the position of the cursor to a perimeter position where a shift in grayscale is detected. Once at the perimeter, the process traces along it, around the building, until it eventually returns to the perimeter starting point. Spatial resolution limits cause the perimeter trace to be somewhat course so that a line straightening algorithm is employed. One result is that the building corner positions become more distinctly defined.

Description: This work concerns fracture mechanics modeling of composite delamination problems. In order to predict delamination resistance, an applied stress intensity factor, K, or energy release rate, G, must be compared to a mode-dependent critical value of K or G from experiment. In the interfacial fracture analysis of most applications and some tests, the mode of crack extension is not uniquely defined. It is instead a function of distance from the crack tip due to the oscillating singularity existing at the tip. In this work, a consistent method is presented of extracting crack extension modes in such cases. In particular, use of the virtual crack closure technique (VCCT) to extract modes of crack extension is studied for cases of a crack along the interface between two in-plane orthotropic materials. Modes of crack extension extracted from oscillatory analyses using VCCT are a function of the virtual crack extension length, delta. Most existing efforts to obtain delta-independent modes of crack extension involve changing the analysis in order to eliminate its oscillatory nature. One such method involves changing one or more properties of the layers to make the oscillatory exponent parameter, epsilon, equal zero. Standardized application of this method would require consistent criteria for identifying which properties can be altered without changing the physical aspects of the problem. Another method involves inserting a thin homogeneous layer (typically referred to as a resin interlayer) along the interface and placing the crack within it. The drawbacks of this method are that it requires increased modeling effort and introduces the thickness of the interlayer as an additional length parameter. The approach presented here does not attempt to alter the interfacial fracture analysis to eliminate its oscillatory behavior. Instead, the argument is made that the oscillatory behavior is non-physical and that if its effects were separated from VCCT quantities, then consistent, delta-independent modes of crack extension could be defined. Knowledge of the near-tip fields in a planar orthotropic material interfacial fracture analysis is used to determine the explicit delta dependence of VCCT parameters. Once this delta dependence is determined, energy release rates are defined with this delta dependence factored out. This modified VCCT method is applied to results from two finite element test cases. It is shown that, as predicted, delta-independent modes of crack extension result. The modified VCCT approach shows potential as a consistent method of extracting crack extension modes. It uses the same information from a finite element analysis (i.e., nodal forces and displacements) as the traditional VCCT method does. The A-independent modes extracted using the modified VCCT approach can also be used as guides to test the convergence of finite element solutions.

Description: The optical method of caustics has been successfully extended to enable stress intensity factors as low as 1MPa square root of m to be determined accurately for central fatigue cracks in 2024-T3 aluminium alloy test panels. The feasibility of using this technique to study crack closure, and to determine the effective stress intensity factor range, Delta K(sub eff), has been investigated. Comparisons have been made between the measured values of stress intensity factor, K(sub caus), and corresponding theoretical values, K(sub theo), for a range of fatigue cracks grown under different loading conditions. The values of K(sub caus) and K(sub theo) were in good agreement at maximum stress, where the cracks are fully open, while K(sub caus) exceeded K(sub theo) at minimum stress, due to crack closure. However, the levels of crack closure and values of Delta K(sub eff) obtained could not account for the variations of crack growth rate with loading conditions. It is concluded that the values of Delta K(sub eff), based on caustic measurements in a 1/square root of r stress field well outside the plastic zone, do not fully reflect local conditions which control crack tip behavior.

Description: Fatigue crack growth tests were conducted on an Fe 510 E C-Mn steel and a submerged arc welded joint from the same material under constant, variable, and random loading amplitudes. Paris-Erdogan's crack growth rate law was tested for the evaluation of m and C using the stress intensity factor K, the J-integral, the effective stress intensity factor K(sub eff), and the root mean square stress intensity factor K(sub rms) fracture mechanics concepts. The effect of retardation and residual stresses resulting from welding was also considered. It was found that all concepts gave good life predictions in all cases.

Description: Newman crack-closure model and the relevant crack growth program were applied to the analysis of crack growth under constant amplitude and aircraft spectrum loading on a number of aluminum alloy materials. The analysis was performed for available test data of 2219-T851, 2024-T3, 2024-T351, 7075-T651, 2324-T39, and 7150-T651 aluminum materials. The results showed that the constraint factor is a significant factor in the method. The determination of the constraint factor is discussed. For constant amplitude loading, satisfactory crack growth lives could be predicted. For the above aluminum specimens, the ratio of predicted to experimental lives, Np/Nt, ranged from 0.74 to 1.36. The mean value of Np/Nt was 0.97. For a specified complex spectrum loading, predicted crack growth lives are not in very good agreement with the test data. Further effort is needed to correctly simulate the transition between plane strain and plane stress conditions, existing near the crack tip.

Description: Load separation is the representation of the load in the test records of geometries containing cracks as a multiplication of two separate functions: a crack geometry function and a material deformation function. Load separation is demonstrated in the test records of several two-dimensional geometries such as compact tension geometry, single edge notched bend geometry, and center cracked tension geometry and three-dimensional geometries such as semi-elliptical surface crack. The role of load separation in the evaluation of the fracture parameter J-integral and the associated factor eta for two-dimensional geometries is discussed. The paper also discusses the theoretical basis and the procedure for using load separation as a simplified yet accurate approach for plastic J evaluation in semi-elliptical surface crack which is a three-dimensional geometry. The experimental evaluation of J, and particularly J(sub pl), for three-dimensional geometries is very challenging. A few approaches have been developed in this regard and they are either complex or very approximate. The paper also presents the load separation as a mean to identify the blunting and crack growth regions in the experimental test records of precracked specimens. Finally, load separation as a methodology in elastic-plastic fracture mechanics is presented.

Description: This paper develops improved stochastic models for the description of a large variety of fatigue crack growth phenomena that occur in components of considerable importance to the functionality and reliability of complex engineering structures. In essence, the models are based on the McGill-Markov and Closure-Lognormal stochastic processes. Not only do these models have the capability of predicting the statistical dispersion of crack growth rates, they also, by incorporating the concept of crack closure, have the capability of transferring stochastic crack growth properties measured under ideal laboratory conditions to situations of industrial significance, such as those occurring under adverse loading and/or environmental conditions. The primary data required in order to be in a position to estimate the pertinent parameters of these stochastic models are obtained from a statistically significant number of replicate tests. In this paper, both the theory and the experimental technique are illustrated using a Ti-6Al-4V alloy. Finally, important structural integrity, reliability, availability and maintainability concepts are developed and illustrated.

Description: Some of the environments and loads experienced by the Space Shuttle or future reusable space vehicles are unique, while others are similar to those encountered by commercial and/or military aircraft. Prior to the Space Transportation System (STS) flights, fatigue loads spectra were generated for the Space Shuttle based on anticipated environments and assumptions that were shown not to be applicable to the actual flight environments the vehicle experienced. This resulted in the need to generate a new cycle of fatigue loads spectra, which was based on measured flight data as well as mission profiles, reflecting the various types of service and operations the vehicle and payloads experienced.

Description: The purpose of this view-graph presentation is a computational investigation of the closed-loop output feedback control of a Euler-Bernoulli beam based on finite element approximation. The observer is part of the classical observer plus state feedback control, but it is finite-dimensional. In the theoretical work on the subject it is assumed (and sometimes proved) that increasing the number of finite elements will improve accuracy of the control. In applications, this may be difficult to achieve because of numerical problems. The main difficulty in computing the observer and simulating its work is the presence of high frequency eigenvalues in the finite-element model and poor numerical conditioning of some of the system matrices (e.g. poor observability properties) when the dimension of the approximating system increases. This work dealt with some of these difficulties.

Description: It is known that an effective control system is the key condition for successful implementation of high-performance magnetic servo systems. Major issues to design such control systems are nonlinearity; unmodeled dynamics, such as secondary effects for copper resistance, stray fields, and saturation; and that disturbance rejection for the load effect reacts directly on the servo system without transmission elements. One typical approach to design control systems under these conditions is a special type of nonlinear feedback called gain scheduling. It accommodates linear regulators whose parameters are changed as a function of operating conditions in a preprogrammed way. In this paper, an on-line learning fuzzy control strategy is proposed. To inherit the wealth of linear control design, the relations between linear feedback and fuzzy logic controllers have been established. The exercise of engineering axioms of linear control design is thus transformed into tuning of appropriate fuzzy parameters. Furthermore, fuzzy logic control brings the domain of candidate control laws from linear into nonlinear, and brings new prospects into design of the local controllers. On the other hand, a self-learning scheme is utilized to automatically tune the fuzzy rule base. It is based on network learning infrastructure; statistical approximation to assign credit; animal learning method to update the reinforcement map with a fast learning rate; and temporal difference predictive scheme to optimize the control laws. Different from supervised and statistical unsupervised learning schemes, the proposed method learns on-line from past experience and information from the process and forms a rule base of an FLC system from randomly assigned initial control rules.