ALBERT

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

Description: Through the Surveyor 3 and 7, and Apollo 11-17 missions a knowledge of the mechanical properties of Lunar regolith were gained. These properties, including material cohesion, friction, in-situ density, grain-size distribution and shape, and porosity, were determined by indirect means of trenching, penetration, and vane shear testing. Several of these properties were shown to be significantly different from those of terrestrial soils, such as an interlocking cohesion and tensile strength formed in the absence of moisture and particle cementation. To characterize the strength and deformation properties of Lunar regolith experiments have been conducted on a lunar soil simulant at various initial densities, fabric arrangements, and composition. These experiments included conventional triaxial compression and extension, direct tension, and combined tension-shear. Experiments have been conducted at low levels of effective confining stress. External conditions such as membrane induced confining stresses, end platten friction and material self weight have been shown to have a dramatic effect on the strength properties at low levels of confining stress. The solution has been to treat these external conditions and the specimen as a full-fledged boundary value problem rather than the idealized elemental cube of mechanics. Centrifuge modeling allows for the study of Lunar soil-structure interaction problems. In recent years centrifuge modeling has become an important tool for modeling processes that are dominated by gravity and for verifying analysis procedures and studying deformation and failure modes. Centrifuge modeling is well established for terrestrial enginering and applies equally as well to Lunar engineering. A brief review of the experiments is presented in graphic and outline form.

Description: Nonthermal radio emissions from earth, Jupiter, Saturn, and Uranus are reviewed. The dominant source of emission at each planet appears to be AKR-like auroral emission in the X-mode. O-mode emissions are substantially responsible. There is a remarkably constant scaling factor relating the total solar wind input power into each planetary system and the AKR-like auroral emissions.

Description: This article presents the results of a fit of a model of the Martian satellite orbits to earthbased and spacecraft-based observations. An assessment of the orbit accuracies is given and the orbits are compared with those obtained by previous investigators.

Description: The petrology and chronology of early lunar crust is examined using the least equivocal of the available petrographic and age data on lunar rock samples, and the possible processes which produced the lunar crust are discussed. The results suggest that the lunar anorthositic crust was formed by about 120 Ma after the primary accretion of the moon at 4.56 Ga. At least some members of the diverse Mg-suites of rocks, such as norites, troctolites, and dunites, crystallized within a very few 100s of Ma after 4.56 Ga. A trace-element-rich material (KREEP) was formed by about 4.3 Ga ago, and this residue was subsequently reworked in melting and impact processes such that most samples which contain it have ages around 3.9-4.0 Ga. The findings also suggest that the onset of ferrous mare basalt volcanism began about 4.33 Ga, much earlier than was once assumed, and was still in process before the end of the most intense period of bombardment (3.9-4.0 Ga ago).

Description: Recent observational and theoretical investigations of the cosmic microwave background radiation (CMBR) are reviewed. Particular attention is given to spectral distortions and CMBR temperature anisotropies at large, intermediate, and small angular scales. The implications of the observations for inflationary cosmological models with curvature fluctuation are explored, and it is shown that the limits determined for intermediate-scale CMBR anisotropy almost rule out a baryon-dominated cosmology.

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

Description: The findings made by the Voyager 2 Neptune encounter are reviewed. Data on the bowshock, magnetic field, magnetosphere, rings, plasma sheet, aurora, moons, and dust of Neptune are discussed. Findings made concerning Triton are summarized.

Description: The composition of lunar regolith and its attendant properties are discussed. Tables are provided listing lunar minerals, the abundance of plagioclase feldspar, pyroxene, olivine, and ilmenite in lunar materials, typical compositions of common lunar minerals, and cumulative grain-size distribution for a large number of lunar soils. Also provided are charts on the chemistry of breccias, the chemistry of lunar glass, and the comparative chemistry of surface soils for the Apollo sites. Lunar agglutinates, constructional particles made of lithic, mineral, and glass fragments welded together by a glassy matrix containing extremely fine-grained metallic iron and formed by micrometeoric impacts at the lunar surface, are discussed. Crystalline, igneous rock fragments, breccias, and lunar glass are examined. Volatiles implanted in lunar materials and regolith maturity are also addressed.

Description: Theoretical models of planetary-ring dynamics are discussed in a detailed analytical review and illustrated with graphs and diagrams. The streamline concept is introduced, and the phenomena associated with the transport of angular momentum are described. Particular attention is then given to (1) broad rings like those of Saturn (shepherding, density-wave excitation, gaps, bending-wave excitation, multiringlet structures, inner-edge shepherding, and the possibility of polar rings around Neptune), (2) narrow rings like those of Uranus (shepherding, ring shapes, and a self-gravity model of rigid precession), and (3) ring arcs like those seen in stellar-occultation observations of Neptune.

Description: Images of Neptune obtained by the narrow-angle camera of the Voyager 2 spacecraft reveal large-scale cloud features that persist for several months or longer. The features' periods of rotation about the planetary axis range from 15.8 to 18.4 hours. The atmosphere equatorward of -53 deg rotates with periods longer than the 16.05-hour period deduced from Voyager's planetary radio astronomy experiment (presumably the planet's internal rotation period). The wind speeds computed with respect to this radio period range from 20 meters per second eastward to 325 meters per second westward. Thus, the cloud-top wind speeds are roughly the same for all the planets ranging from Venus to Neptune, even though the solar energy inputs to the atmospheres vary by a factor of 1000.

Description: The theory of turbulent plumes maintained above steady line sources of buoyancy is worked out in detail within the limitations of Taylor's entrainment assumption. It is applied to the structure of a pure plume injected into a stably stratified atmosphere. Volcanic basalt eruptions that develop from long, narrow vents create line source plumes, which rise well above the magmatic fire fountains playing near the ground level. The eruption of Laki in 1783 may provide an example of this style of eruption. Flood basalts are more ancient examples. Evidence of enormous fissure eruptions that occurred in the past on Mars and Venus also exists. Owing to the different properties of the atmospheres on these two planets from those on the earth, heights of line source plumes are expected to vary in the ratios 1:6:0.6 (earth:Mars:Venus). It is very unlikely that the observed increase of sulfur dioxide above the Venusian cloud deck in 1978 could have been due to a line source volcanic eruption, even if it had been a flood basalt eruption.

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

Description: An analytical and numerical study of the stability of polar rings around Neptune is presented. The stability proofs are based on various methods used to study gas disks in galaxies. It is shown that stable polar rings can exist despite energy dissipation by collisions between particles. Also, four equilibrium orientations which pass nearly over the pole of Neptune are found, two of which are stable in the presence of dissipation.

Description: Issues related to the establishment of lunar observatories are briefly addressed. The advantages of placing an observatory in a crater at one of the poles, where permanent darkness exists, are pointed out, and the methods required to emplace and operate such an observatory are considered. Planning for the installation of the first set of observatory instruments is discussed.

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

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

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

Description: The phenomena of the Saturn system discovered by the Voyager missions are addressed. The characteristics of the rings, including grooves, spikes, and warps, are described, and the discovery of bombarded rings and young rings is discussed. The unique and unexplained properties of Enceladus are summarized, and the possibility that some moons of Saturn have been reassembled is addressed. The interaction of moons with the rings is examined, and the possibility that Titan's atmosphere may contain complex molecules that are the precursors of life is considered.

Description: Spectral data from earth observations have indicated the presence of N2 and CH4 on Triton. This paper outlines the use of the 1-D radiative-convective model developed for Titan to calculate the current pressure of N2 and CH4 on Triton. The production of haze material is obtained by scaling down from the Titan value. Results and predictions for the Voyager Triton encounter are as follows: A N2-CH4 atmosphere on Triton is thermodynamically self consistent and would have a surface pressure of approximately 50 millibar; due to the chemically produced haze, Triton has a hot atmosphere with a temperature of approximately 130 K; Triton's troposphere is a region of saturation of the major constituent of the atmosphere, N2.

Description: The process by which ambient gases can be implanted into silicates by shocks was investigated by analyzing the noble-gas content of several experimentally and naturally shocked silicate samples. The retentivity of shock-implanted gas during stepwise heating in the laboratory was defined in terms of two parameters, namely, the activation energy for diffusion and the extraction temperature at which 50 percent of the gas is released, both of which correlate with the shock pressure. The experiments indicate that, with increasing shock pressure, gas implantation occurs through an increasing production of microcracks/defects in the silicate lattice. The degree of annealing of these defects control the degree of diffusive loss of implanted gas.

Description: The long-term modulation of Saturn's nonthermal radio emission in the kilometric wavelength range has been studied based upon data obtained by Voyagers 1 and 2. A comparison of the ballistic and hydrodynamic propagation of solar wind features from the spacecraft to Saturn allows the uncertainty inherent in the projection to be determined. The results confirm the previous suggestion that momentum, ram pressure, and kinetic energy flux are the primary solar wind parameters that drive the nonthermal radio emission. It is suggested that, under certain conditions and for limited periods of time, the magnetic properties and time derivatives of the solar wind have increased importance.

Description: It has been observed from the plasma experiments on the Pioneer Venus Orbiter that the altitude of the upper boundary of the ionosphere decreases in response to increasing solar wind dynamic pressure. However, at pressures above about 2.5 x 10 to the -8th dynes/sq cm, the further decrease in the ionopause height is rather small. Following the model of Cloutier et al. (1969), it is suggested that during high solar wind conditions, when the ionopause is formed at lower altitudes, the solar wind induces vertical and horizontal flows which sweep away the ionospheric plasma that is produced locally by photoionization. As a result, a disturbed photodynamical ionosphere is formed which has the scale height of the ionizable neutral constituent. It is shown that such a photodynamical ionosphere is observed at the subsolar ionopause under these conditions. As a consequence of this interaction, the ionopause altitude is observed to follow the small-scale height of the ionizable species, atomic oxygen, showing only small changes with solar wind pressure.

Description: Asteroids appear in light of telescopic and meteority studies to be the most accessible repositories of early solar system history available. In the cooler regions of the outer asteroid belt, apparently unaffected by severe heating, the C, P, and D populations appear to harbor significant inventories of volatiles; the larger primordial belt population may have had an even greater percentage of volatile-rich, low-albedo asteroids, constituting a potent asteroid for veneering early terrestrial planet atmospheres. The volatile-rich asteroids contain carbon, structurally bound and adsorbed water, as well as remnants of interstellar material predating the solar system.

Description: The development and implementation of the Voyager missions are reviewed. The interplanetary missions preceding Voyager are discussed, focusing on the technological development leading up to the Voyager spacecraft. The main results from Voyager observations of Jupiter, Saturn, and Uranus are outlined. Also, consideration is given to the prospects for observations of Neptune.

Description: Isotopic analysis of nesquehonite recovered from the surface of the LEW 85320 H5 ordinary chondrite shows that the delta C-13 and delta O-18 values of the two generations of bicarbonate (Antarctic and Texas) are different: delta C-13 = + 7.9 per mil and + 4.2 per mil; delta O-18 = + 17.9 per mil and + 12.1 per mil, respectively. Carbon isotopic compositions are consistent with equilibrium formation from atmospheric carbon dioxide at - 2 + or - 4 C (Antarctic) and + 16 + or - 4 C (Texas). Oxygen isotopic data imply that the water required for nesquehonite precipitation was derived from atmospheric water vapor or glacial meltwater which had locally exchanged with silicates, either in the meteorite or in underlying bedrock. Although carbonates with similar delta C-13 values have been identified in the SNC meteorites EETA 79001 and Nakhla, petrographic and temperature constraints argue against their simply being terrestrial weathering products.

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

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

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

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

Description: Results from a Monte Carlo simulation of gamma-ray production under astrophysical conditions indicate that it is possible to produce gamma-ray spectra as hard as that of 2CG013+00 reported by COS-B using a monoenergetic beam of cosmic rays hitting a thin target material. It is suggested that the low mass X-ray binary GX13+1 might provide such a mechanism to produce the observed gamma-rays from the region.

Description: If cyclotron scattering, rather than absorption, is responsible for the line features observed recently in two gamma-ray burst spectra (Murakami et al., 1988), then the second and higher harmonics are due to resonant scattering events that excite the electron to Landau levels above the ground state. Here, relativistic Compton scattering cross sections are used to estimate the expected ratio of third to second harmonics in the presence of Doppler broadening. At the field strength (1.7 TG) required to give first and second harmonics at 19 keV and 38 keV, there should be no detectable third harmonic in the spectrum.

Description: Gamma-ray observations obtained by the SMM gamma-ray spectrometer in the energy range 4.1-6.4 MeV are used to provide limits on the possible radiative decay of neutrinos and axions emitted by SN 1987A. For branching ratio values for the radiative decay modes of less than about 0.0001, the present limits are more stringent than those based upon the photon flux from decaying relic neutrinos. The data are also used to set an axion mass limit.

Description: SMM gamma-ray spectrometer data are examined to look for gamma-ray emission coincident with the about-10-s neutrino burst from SN 1987A. The absence of a detectable signal suggests that the energy radiated into MeV gamma rays by neutrino decay (or any other process) is less than 10 to the -10th of that in supernova neutrinos above 9 MeV. The results are used to set a direct limit on the lifetime of any massive neutrino type generated in the core collapse leading to SN 1987A.

Description: A primary concern in the formation of a multinational Mars-exploration program is that of technology transfer, which would be most acute in the U.S./U.S.S.R. case but would exist in lesser degree among other nations. Another concern is that of management-complexity, which could inflate total costs and substantially counteract the anticipated benefits of spreading program costs among a number of nations. It is presently suggested that these problems can be substantially reduced by having each nation become responsible for one (or more) complete flight system. Each nation's role in mission operations must be clear, and the cleanliness of interfaces among flight systems must carry over into mission operations.

Description: The spectrum of antiprotons from dark matter annihilation are calculated using the Lund Monte Carlo program, and simple analytic expressions for the spectrum and low-energy antiproton/proton ratio are derived. Comparing the results with recent upper limits on low energy antiprotons, it is concluded that the reported 4-13 GeV antiproton flux cannot be accounted for by dark matter annihilation. The new upper limits do not provide useful constraints on dark matter particles. They restrict the annihilation rate and imply that annihilation gamma ray and e(+) fluxes would be far below the fluxes produced by cosmic-ray collisions. It may be possible to look for a dark matter halo annihilation signal at antiprotons energies below 0.5 GeV, where the flux from cosmic-ray collisions is expected to be negligible.

Description: The HEAO 1 A-2 all-sky survey data have been analyzed to investigate the autocorrelation function of the surface brightness fluctuations of the extragalactic 2-10 keV X-ray background on angular scales ranging from 3 to 27 deg. The derived upper limits are compatible with optical estimates of the spatial correlation function of Abell's clusters and set interesting constraints on possible AGN clustering.

Description: The physical theory underlying calculations of nuclear absorption and fragmentation cross sections and EM dissociation cross sections is reviewed; typical results are presented in graphs; and the adequacy of the experimental data bases for each type of cross section are discussed. The emphasis is on computations applicable to the estimation of secondary radiation in shielded spacecraft.

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

Description: Observations of the stratospheric IR emission structure on Saturn are presented. The high-spatial-resolution global images show a variety of new features, including a narrow equatorial belt of enhanced emission at 7.8 micron, a prominent symmetrical north polar hotspot at all three wavelengths, and a midlatitude structure which is asymmetrically brightened at the east limb. The results confirm the polar brightening and reversal in position predicted by recent models for seasonal thermal variations of Saturn's stratosphere.

Description: An overview of the Voyager 2 encounter with Neptune is presented, including a brief discussion of the trajectory, the planned observations, and highlights of the results described in the 11 companion papers. Neptune's blue atmosphere has storm systems reminiscent of those in Jupiter's atmosphere. An optically thin methane ice cloud exists near the 1.5-bar pressure level, and an optically thick cloud exists below 3 bars. Neptune's magnetic field is highly tilted and offset from the planet's center; it rotates with a period of 16.11 hours. Two narrow and two broad rings circle the planet; the outermost of these rings has three optically thicker arc segments. Six new moons were discovered in circular prograde orbits, all well inside Triton's retrograde orbit. Triton has a highly reflective and geologically young surface, a thin nitrogen atmosphere, and at least two active geyser-like plumes.

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

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

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

Description: Recent models of the internal structure of Pluto and Charon, made possible by analysis of the Pluto/Charon mutual events are reviewed. At a mean density of just over 2 g/cu cm and a predicted rock mass fraction of roughly 0.7, the Pluto/Charon system is significantly rockier than the icy satellites of the giant planets, a contrast which may reflect its formation in a CO-rich outer solar nebula rather than a circumplanetary nebula. Pluto and Charon may in fact be so rocky that they lost volatiles early in their history (possibly during a Charon-forming impact event), although this is still an open issue.

Description: Preliminary results from electric field measurements in the environment of Mars using the plasma-wave system on board Phobos 2 are reported. Electron-plasma oscillations observed upstream of the bow shock correspond to a solar-wind density of 2/cu cm. The shock-foot boundary was crossed up to three times on each orbit. The shock ramp was detected at altitudes between 0.45 and 0.75 Mars radii R(M) above the planetary surface. The density increased by about a factor of two at the ramp. The shock position, although variable, seems to be consistent with previous measurements. The downstream magnetosheath contained broadband electric-field noise below the plasma frequency. The boundary of th obstacle, or plasmapause, was crossed at altitudes of the order of 0.28 R(M); the cold plasma density was highly variable within the planetopause and reached the unexpected value of 700/cu cm on the third orbit, at 0.25 R(M) altitude. Bursts of waves with frequencies below the electron cyclotron frequency occur within the planetopause.

Description: The magnetic fields of Mars have been measured from Phobos 2 with high temporal resolution in the tail and down to an 850-km altitude. During four successive highly elliptical orbits, the position of the bow shock as well as that of a transition layer, the 'planetopause', were identified. Subsequent circular orbits at 6000-km altitude provided the first high-resolution data in the planetary tail and indicate that the interplanetary magnetic field mainly controls the magnetic tail. Magnetic turbulence was also detected when the spacecraft crossed the orbit of Phobos, indicating the possible existence of a torus near the orbit of this moon.

Description: The formation of the ice-rock cores of the giant planets by density wave-assisted accretion is outlined. The process could be rapid (100,000-1,000,000 yr) and completed within the probable lifetime of the solar nebula. The mechanism works for both Jupiter and Saturn and does not require a large excess of mass over that believed present in their cores.

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

Description: This paper identifies and documents several well-preserved examples of Martian strike-slip faults and examines their relationships to wrinkle-ridges. The strike-slip faulting predates or overlaps periods of wrinkle-ridge growth southeast of Valles Marineris, and some wrinkle ridges may have nucleated and grown as a result of strike-slip displacements along the echelon fault arrays. Lateral displacements of several km inferred along these arrays may be related to tectonism in Tharsis.

Description: Results are presented on the long-term behavior of the main peak electron density in the Venus ionosphere during the solar cycle 21, based on 104 radio occultation measurements of the vertical electron density profile in the dayside ionosphere of Venus carried out aboard the Pioneer Venus Orbiter spacecraft (along with published data on 11 Venera 9-10 measurements). The equation representing the electron temperature at h = 140 km is presented. The results imply that the electron temperature at h = 140 km decreased by about 25 percent from solar maximum to solar minimum, compared to a decrease of 50-75 percent above 200 km found by Kliore and Mullen (1989).

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

Description: Theoretical models of early Martian atmospheric evolution describe the maintenance of a dense CO2 atmosphere and a warm, wet climate until the end of the heavy-bombardment phase of impacting. However, the presence of very young, earthlike fluvial valleys on the northern flank of Alba Patera conflicts with this scenario. Whereas the widespread ancient Martian valleys generally have morphologies indicative of sapping erosion by the slow outflow of subsurface water, the local Alba valleys were probably formed by surface-runoff processes. Because subsurface water flow might be maintained by hydrothermal energy inputs and because surface-runoff valleys developed late in Martian history, it is not necessary to invoke drastically different planet-wide climatic conditions to explain valley development on Mars. The Alba fluvial valleys can be explained by hydrothermal activity or outflow-channel discharges that locally modified the atmosphere, including precipitation and local overland flow on low-permeability volcanic ash.

Description: On the Mars rover sample-return mission, the rover vehicle will collect and select samples from different locations on the Martian surface to be brought back to earth for laboratory studies. It is anticipated that an in situ energy-dispersive X-ray fluorescence (XRF) spectrometer will be on board the rover. On such a mission, sample selection is of higher priority than in situ quantitative chemical anlaysis. With this in mind, a pattern recognition technique is proposed as a simple, direct, and speedy alternative to detailed chemical analysis of the XRF spectra. The validity and efficacy of the pattern recognition technique are demonstrated by the analyses of laboratory XRF spectra obtained from a series of geological samples, in the form both of standardized pressed pellets and as unprepared rocks. It is found that pattern recognition techniques applied to the raw XRF spectra can provide for the same discrimination among samples as a knowledge of their actual chemical composition.

Description: The present photoelectric lightcurve observations of 44 Nysa at one aspect, and 64 Angelina at both polar aspects, indicate the presence of an opposition spike or brightening of approximately 0.25 mag within several deg of zero-phase angle. Since the three curves are identical within the bounds of observational scatter, the opposition spike cannot be judged anomalous within the asteroid taxonomic class E. The general similarity of the observed phase curves to those of the Uranian satellites and the rings of Saturn supports the status of the spike as an ordinary property of moderate-to-high albedo atmosphereless surfaces.

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

Description: In order to ascertain the features of organic compound-production in planetary atmospheres under the effects of plasmas and shocks, various mixtures of N2, CH4, and H2 modeling the atmosphere of Titan were subjected to discrete sparks, laser-induced plasmas, and UV radiation. The experimental results obtained suggest that UV photolysis from the plasma is an important organic compound synthesis process, as confirmed by the photolysis of gas samples that were exposed to the light but not to the shock waves emitted by the sparks. The thermodynamic equilibrium theory is therefore incomplete in the absence of photolysis.

Description: An account is given of the results of photometric light curve observations in the V band for asteroids, yielding estimates of the mean and maximum reduced magnitudes for each object. On the basis of fits of the H-G magnitude relation for 33 objects, the mean values of the slope parameter are examined for different taxonomic classes; these values are then applied to the analysis of less complete data sets. While in the case of the moderate albedo S and M class asteroids the H-G relation appears to fit the available data, the relation for dark asteroids appears to predict more of an opposition than is characteristically present.

Description: Results are reported from the Low Energy Antiproton Experiment (LEAP), a balloon-borne instrument which was flown in August, 1987. No evidence of antiproton fluxes is found in the kinetic energy range of 120 MeV to 360 MeV, at the top of the atmosphere. The 90-percent is found confidence upper limit on the antiproton/proton ratio in this energy range is 3.5 x 10 to the -5th. In particular, this new experiment places an upper limit on the flux almost an order of magnitude below the reported flux of Buffington et al. (1981).

Description: Direct measurements on cosmic ray protons through iron above about 1 TeV/amu have been performed in a series of balloon-borne experiments with emulsion chambers. The measured energy spectra of protons and helium are power laws with exponents of 2.77 + or - 0.09 and 2.72 + or - 0.11 in the energy range 5 to 500 TeV and 2 to 50 TeV/amu, respectively. The proton spectrum shows no evidence of the steepening near 2 TeV which was reported by other experiments. Helium has a slightly higher intensity compared to extrapolations from lower energy measurements. The heavier elements, carbon to sulfur, show a small tendency for intensity enhancement in the relative abundance obove 10 TeV/amu.

Description: Voyager imaging, radio occultation, and stellar occultation data for the regular structure of Saturn's inner Cassini Division are presently analyzed. The regular optical depth variation observed by the radio occultation experiment scan and the feature noted in Voyager images is the same structure, namely the gravitational wakes of two 10-km radius satellites orbiting within the division. The structure is azimuthally symmetric, and is judged to rule out the possibility that large moonlets may be responsible for the observed structure.

Description: TEM and IR spectroscopy investigations of the interplanetary dust particles (IDPs) collected in the stratosphere have shown the majority of IDPs in the layer-lattice silicate and pyroxene classes to not have been heated to temperatures above 600 C during atmospheric entry. This implies that they arrive at the upper atmosphere with low geocentric encounter velocities, and limits the possible encounter trajectories for these particles to relatively circular prograde orbits. On this basis, it is judged unlikely that these IDPs are from earth-crossing comets or asteroids; collected IDPs dominated by olivine include a larger portion of above-600 C-heated particles, suggesting their capture from more eccentric orbits.

Description: The present examination of the west Deuteronilus Mensae region of Mars notes the changes in fretted terrain across the gradational boundary from uplands to lowlands to include a reduction of canyon wall slopes and depths, so that the fretted terrain north of the gradational boundary appears matted, but not obscured. The two process-classes that may account for the lateral overlap are the erosion of stratified upland terrain, and the deposition of plains materials onto the sloping upland margin and fretted terrain. A variety of plains-emplacement mechanisms is considered.

Description: The compositional contrast between the giant-planet satellites and the significantly rockier Pluto/Charon system is indicative of different formation mechanisms; cosmic abundance calculations, in conjunction with an assumption of the Pluto/Charon system's direct formation from solar nebula condensates, strongly suggest that most of the carbon in the outer solar nebula was in CO form, in keeping with both the inheritance from the dense molecular clouds in the interstellar medium, and/or the Lewis and Prinn (1980) kinetic-inhibition model of solar nebula chemistry. Laboratory studies of carbonaceous chondrites and Comet Halley flyby studies suggest that condensed organic material, rather than elemental carbon, is the most likely candidate for the small percentage of the carbon-bearing solid in the outer solar nebula.

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

Description: An attempt is made to provide a constraint on the combined mass of Janus and Epimetheus from an analysis of Voyager I and Voyager 2 data and ground-based observations obtained during the 1966 and 1980 ring plane crossings. The results of the analysis presented here suggest that the total mass is 2.59 + or - 0.26 x 10 to the 21st g, the mass ratio is 3.61 + or - 0.01, and Janus' density is 0.67 + or - 0.10 g/cu cm. The low density of Janus is attributed to its porosity rather than composition.

Description: A coordinated multinational Mars Surveyor Program involving a series of robotic missions is proposed in order to select worthwhile human landing sites for discoveries, for safe operations, and for testing and proving technologies for making human flight to Mars possible. Some characteristics motivating exploration of the planet are briefly discussed, including the possibility of life, geological features, and meteorological conditions. The necessity for preliminary exploration of Mars by robots prior to human exploration is discussed, and the rationale behind a multinational approach for a Mars Surveyor Program is presented.

Description: Arecibo high-resolution (1.5 to 2 km) radar data of Venus for the area extending from Beta Regio to western Eisila Regio provide strong evidence that the mountains in Beta and Eisila Regiones and plains in and adjacent to Guinevere Planitia are of volcanic origin. Recognized styles of volcanism include large volcanic edifices on the Beta and Eisila rises related to regional structural trends, plains with multiple source vents and a mottled appearance due to the ponding of volcanic flows, and plains with bright features surrounded by extensive quasi-circular radar-dark halos. The high density of volcanic vents in the plains suggests that heat loss by abundant and widely distributed plains volcanism may be more significant than previously recognized. The low density of impact craters greater than 15 km in diameter in this region compared to the average density for the higher northern latitudes suggests that the plains have a younger age.

Description: In February and March 1989 the Phobos 2 spacecraft took 37 television images of Phobos from a distance of 190-1100 km. The data are being used to update the three-dimensional model of Phobos, to provide improved determinations of its density and orbital dynamics, and to study its surface color, composition, and texture. Preliminary findings are presented here which include different integrated photometric behavior in visible and near-infrared bands, observation of a region immediately west of Stickney which is relatively free of large grooves, the prevalence of bright rims on grooves and younger craters, and low bulk density.

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

Description: Requirements on orbital background counting rates are described for the UV and visible detectors on the planned second generation instrument Space Telescope Imaging Spectrograph (STIS) and for the proposed Lyman Far Ultraviolet Spectroscopic Explorer. Expected background rates and shielding requirements are analyzed for detector types in both instruments for low earth orbits, and for high orbits in the case of Lyman.

Description: Gamma-ray background spectra, obtained in the energy range from 45 keV to 10 MeV with the cooled high-resolution Ge spectrometer on the HEAO-3 satellite over a 50-day period in 1979-1980 are reported and analyzed. The instrument design and performance characteristics are reviewed; the data-analysis procedures are outlined; about 130 background lines are listed in a table with energies, count rates, and tentative identifications; and the spectra are presented graphically and discussed in detail.

Description: The highly anisotropic nature of the radiation in the low-earth orbit has been ignored for most spacecraft shielding calculations made to date because the standard environmental models describe the omnidirectional flux only, because the varying attitude of the spacecraft in the environment is assumed to average out the effect and because of the added complexity of the calculation. The Space Station is planned to be stabilized with respect to the velocity vector and local vertical. Thus it will pass through the South Atlantic Anomaly where most of the radiation flux is encountered in much the same attitude on each pass. Any calculation including a complex shielding geometry should thus consider the angular distributuon of the incident radiation. An approximate trapped proton angular distribution is presented which includes both the 'pan caked' distribution relative to the magnetic field direction and the east-west effect which is energy dependent. This distribution is then used with a planar shielding geometry to obtain an estimate of the effect of the anisotropy on radiation dose rates in spacecraft.

Description: The near-earth radiation environment is characterized, summarizing published data on trapped and transient charged particles and their potential effects on hardware systems and the crews of manned missions. Topics examined include the role of the magnetosphere, the five radiation domains, cyclic and sporadic variations in the radiation environment, the potential effect of a high-altitude nuclear explosion, NASA empirical models for predicting trapped proton and electron fluxes, and the South Atlantic anomaly and the estimation of flux-free periods. Consideration is given to solar cosmic rays and heavy ions, Galactic cosmic rays, geomagnetic shielding, secondary radiation, the design of shielding systems, variables affecting dose evaluations, and ionizing-radiation doses. Extensive diagrams, graphs, and tables of numerical data are provided.

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

Description: The nearly ubiquitous presence of a negative Eu anomaly in the mare basalts has been suggested to indicate prior separation and flotation of plagioclase from the basalt source region during its crystallization from a lunar magma ocean (LMO). Are there any mare basalts derived from a mantle source which did not experience prior plagioclase separation? Crystal chemical rationale for REE substitution in pyroxene suggests that the combination of REE size and charge, M2 site characteristics of pyroxene, fO2, magma chemistry, and temperature may account for the negative Eu anomaly in the source region of some types of primitive, low TiO2 mare basalts. This origin for the negative Eu anomaly does not preclude the possibility of the LMO as many mare basalts still require prior plagioclase crystallization and separation and/or hybridization involving a KREEP component.

Description: Pioneer Venus ion composition measurements are used to study the Venus ionosphere during solar minimum. It is suggested that the topside electron density profile at Venus during solar minimum has two distinct regimes. One beween 140 and 180 km is dominated by O2(+) ions which are in photochemical equilibrium. The other regime is above 180 km and is dominated by O(+) ions which are disturbed by the solar wind induced plasma transport. For Pioneer Venus, Mariner 10, and Venera 9 and 10 data, it is found that Venus exhibits a photodynamical type of ionopause during solar minimum.

Description: Voyager 2's IR observations of Neptune encompass thermal emissions and broadband radiometer measurements of reflected solar radiation. Temperature maps were obtained for the planet between 80 deg S and 30 deg N for two atmospheric layers, one in the lower stratosphere and the other in the troposphere. The relatively warm pole and equator, with cooler midlatitudes, are qualitatively similar to Uranus, despite the two planets' very different obliquities and internal heat fluxes. Powerful wavelike longitudinal thermal structure is noted, of which some appears to be associated with the Great Dark Spot; a localized cold region uncorrelated with any visible feature is found in the lower stratosphere.

Description: Detection of very intense short radio bursts from Neptune was possible as early as 30 days before closest approach and at least 22 days after closest approach. The bursts lay at frequencies in the range 100 to 1300 kilohertz, were narrowband and strongly polarized, and presumably originated in southern polar regions of the planet. Episodes of smooth emissions in the frequency range from 20 to 865 kilohertz were detected during an interval of at least 10 days around closest approach. The bursts and the smooth emissions can be described in terms of rotation in a period of 16.11 + or - 0.05 hours. The bursts came at regular intervals throughout the encounter, including episodes both before and after closest approach. The smooth emissions showed a half-cycle phase shift between the five episodes before and after closest approach. This experiment detected the foreshock of Neptune's magnetosphere and the impacts of dust at the times of ring-plane crossings and also near the time of closest approach. Finally, there is no evidence for Neptunian electrostatic discharges.

Description: Based on the success of several two-dimensional (latitude, longitude) linear barotropic instability models at matching some of the observed characteristics of the cloud level, polar region of the Venus atmosphere, a more realistic, linear, three-dimensional (height, latitude and longitude) model has been developed to further test the hypothesis that the observed features can be described by linear instability theory. The approach taken is to vary the model input parameters to see whether it is possible to produce modes that resemble the observations of wave activity and to compare those input parameters with other observations of the mean state. Sensitivity studies show that in addition to a well-documented dependence on the mean zonal wind, the growth and propagation of unstable modes depends on the latitude variation of the mean temperature (and hence static stability). These studies have led to the specification of a model basic state wind and temperature field that produces modes which are matched to observations of spatial structure, preferred wavenumber and phase speed of the polar disturbances. Wavenumber 2 is found to have the shortest growth time and unlike the two-dimensional results, wavenumbers 1-3 share a nearly common period of about 3 days. The derived basic state has a temperature field that is quite similar to Pioneer Venus observations; however, in some regions the model basic state wind field departs from cyclostrophic values based on temperature observations.

Description: A system is developed for making measurements of spatial fluctuations in the Cosmic Microwave Background at the 3-mm wavelength, on an angular scale of .5 to 5 degrees. The system includes a telescope with a Gaussian beam with an FWHM of 20 to 50 arcmin, an SIS coherent receiver operating around 90 GHz, and for balloon flights, a pointing system capable of 1 arc-minute RMS stabilization. Results are reported from ground-based measurements made from the South Pole station during December, 1988.

Description: Observations of the continuum emission from warm dust clouds at 230 GHz, or 1300 microns, enables the determination of the frequency dependence of the optically thin longwave emission. Integrating the emission over the solid angle of the clouds gives a flux independent of the beam size and of the internal temperature structure of the clouds. The frequency resolving power of 64 allows the correction of these fluxes for the contribution of free-free emission from nearby H II regions at a price of reduced sensitivity. These observations were combined with similar observations made by others in the submillimeter and far infrared regimes to determine the continuum spectra of the dust clouds. Mean characteristics were determined for these clouds, fitting the spectra with simple transfer equations. Results are briefly discussed.

Description: Models of circumstellar dust shells around ultracompact (UC) H II regions were constructed that accurately fit the observed IR flux distributions. The models assume spherically symmetric dust shells illuminated by stars whose bolometric luminosity is inferred from the integrated FIR flux densities. Assuming ionization by a single zero age main sequence (ZAMS) star, the relations of Panagia were used to infer the stellar radius and effective temperature for a given luminosity. The grain mixture in the dust shell consists of bare graphite and silicate grains with the optical properties of Draine and Lee and the size distribution of Mathis et al. The computer code of Wolfire et al was used to solve the radiative transfer equations through a spherical dust shell. The model provides monochromatic luminosities, dust temperatures, and opacities through the shell. Aside from the stellar and dust properties, the only other input parameters to the model are the distance to the shell, the form of its density distribution, and its outer radius. Predictions of the model are compared with observations of a typical UC H II region and the run of dust temperature with radius and the optical depth with frequency are discussed.

Description: High resolution far IR observations at 50 and 100 microns were made of the young stellar object (YSO), IRAS 16293-2422. The observations are part of a systematic high resolution study of nearby YSO's. The purpose is to obtain resolution in the far IR comparable to that at other wavelengths. Until recently, the high resolution that has been available in the far IR has been from either IRAS (angular resolution of approx 4 min) or the KAO using standard FIR photometry (approx 35 sec). With scanning techniques, it is possible to obtain 10 sec resolution on bright sources. Such a resolution is necessary to better determine the physical conditions of the YSO, and to compare with model of star formation. In order to better constrain the models for the source, the YSO was observed at both 50 and 100 microns on several flights in 1988 April from the KAO. Estimates are presented of the size both along the major and minor axis of the disk, as well as estimates of the dust temperature and 100 micron opacity for the YSO.

Description: The 3.3 micron aromatic feature has been detected in the diffuse galactic emission with the AROME balloon-borne instrument. The results are presented in the form of an map of the 3.3 micron feature's intensity. The AROME instrument consists in a Cassegrain telescope with wobbling secondary mirrors and a liquid/solid nitrogen cooled photometer. The instrumental output is modified by the impulse response of the system. So the galactic surface brightness was restored in Fourier space by an inverse optimal filtering. The map of the feature's intensity is presented for a region of galactic coordinates. All the known H II giant molecular cloud complexes are visible in the 3.3 micron feature emission showing a good correlation with the infrared dust emission.

Description: Recent models of interstellar extinction have shown the importance of understanding both the UV and IR properties of interstellar dust grains. IRAS data have shown variations in 60 and 100 micron emissions presumably due to the presence of IR cirrus, while recent observations in the UV by Fitzpatrick and Massa have identified components in the UV extinction curve which vary in different star regions. A Draine and Anderson model connects these results by proposing that different size variations in interstellar grains would cause distinct changes in both the IR emission and the UV extinction. In order to test this model it is necessary to make observations in well defined locations away from peculiar extinction regions. In the infrared this means looking away from the galactic plane so as to limit non-local sources of IR radiation. Two open clusters that are out of the galactic plane and which contain a number of late B and early A stars suitable for UV extinction studies, and whose IRAS data show variations in the 60/100 micron ratio were studied. Based on the Drain and Anderson model, variations were expected in their UV extinction curves that correlate with the IR cirrus emission.

Description: If the 11.3 micron emission feature seen in the spectra of many planetary nebulae, H II regions, and reflection nebulae is attributable to polycyclic aromatic hydrocarbons (PAHs), then additional features should be present between 11.3 and 13.0 microns. Moderate resolution spectra of NGC 7027, HD 44179, BD+30 deg 3639, and IRAS 21282+5050 are presented which show evidence for new emission features centered near 12.0 and 12.7 microns. These are consistent with an origin from PAHs and can be used to constrain the molecular structure of the family of PAHs responsible for the infrared features. There is an indication that coronene-like PAHs contribute far more to the emission from NGC 7027 than to the emission from HD 44179. The observed asymmetric profile of the 11.3 micron band in all the spectra is consistent with the slight anharmonicity expected in the C-H out-of-plane bending mode in PAHs. A series of repeating features between 10 and 11 microns in the spectrum of HD 44179 suggests a simple hydride larger than 2 atoms is present in the gas phase in this object.

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

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

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

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

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

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

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

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

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

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

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