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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: The application of a generalized optimality criteria to framed structures is presented. The optimality conditions, Lagrangian multipliers, resizing algorithm, and scaling procedures are all represented as a function of the objective and constraint functions along with their respective gradients. The optimization of two plane frames under multiple loading conditions subject to stress, displacement, generalized stiffness, and side constraints is presented. These results are compared to those found by optimizing the frames using a nonlinear mathematical programming technique.

Description: Recent developments in multilevel optimization are briefly reviewed. The general nature of the multilevel design task, the use of approximations to develop and solve the analysis design task, the structure of the formal multidiscipline optimization problem, a simple cantilevered beam which demonstrates the concepts of multilevel design and the basic mathematical details of the optimization task and the system level are among the topics discussed.

Description: A two-level design optimization metholology is described. A progress report of its application to Printed Wiring Board (PWB) assembly examples is given. The design of PWB assemblies is a complex task which is generally conducted as a sequential process. Individual PWBs are usually designed first, followed by the composition of the PWBs into an assembly. As a result, optimizing design considerations such as assembly reliability cannot be accomplished. This study showed that a two-level decomposition method can be employed to optimize for reliability at both the PWB- and the assembly-level in a coupled manner. The two-level decomposition method also resolved the mixed-integer nonlinear programming nature of the problem rather easily.

Description: Independent Modal Space Control (IMSC) is a technique that is often used for the control of large order structural systems. The pertinent optimization problem in the simultaneous design and control of structures is a min - min problem that minimizes with respect to the structural design variables, the minimum value of the performance index with respect to the control forces obtained using the IMSC technique. The minimization process requires derivatives of eigenvalues and eigenvectors with respect to the design variables. These derivatives can be computed by a rather involved analytical procedure or a relatively simple finite difference procedure. The computer cost effectiveness of these two procedures for the derivative calculations is examined.

Description: The objective of this presentation is to introduce the attendees to the DYSCO program. The emphasis will be on the features which make it multidisciplinary. DYSCO is a very general and versatile software program which couples and solves dynamic systems. It was initiated in the late 1970's in response to a helicopter analysis requirement. The system development, however, resulted in an executive which was completely separated from any particular area of technology, except that of second order ordinary differential equations. During the course of its development, it was funded by the Army Aviation Applied Technology Directorate, the Air Force Wright Aeronautical Laboratories, and by the Kaman Aerospace Corporation. It is completely written in FORTRAN and is operational on IBM and VAX computers.

Description: A numerical method is presented for design sensitivity analysis, using an iterative-method reanalysis of the structure generated by a small perturbation in the design variable; a forward-difference scheme is then employed to obtain the approximate sensitivity. Algorithms are developed for displacement and stress sensitivity, as well as for eignevalues and eigenvector sensitivity, and the iterative schemes are modified so that the coefficient matrices are constant and therefore decomposed only once.

Description: The considerations and the resultant approach used to implement design sensitivity capability for grids into a large scale, general purpose finite element system (MSC/NASTRAN) are presented. The design variables are grid perturbations with a rather general linking capability. Moreover, shape and sizing variables may be linked together. The design is general enough to facilitate geometric modeling techniques for generating design variable linking schemes in an easy and straightforward manner. Test cases have been run and validated by comparison with the overall finite difference method. The linking of a design sensitivity capability for shape variables in MSC/NASTRAN with an optimizer would give designers a powerful, automated tool to carry out practical optimization design of real life, complicated structures.

Description: Structural optimization has been available to the structural analysis community as a tool for many years. The popular use of displacement method finite-element techniques to analyze linearly elastic structures has resulted in an ability to calculate the weight and constraint gradients inexpensively for numerical optimization of structures. Here, recent experiences in the investigation and use of structural optimization are discussed. In particular, experience with the commercially available ADS/NASOPT code is addressed. An overview of the ADS/NASOPT procedure and how it was implemented is given. Two example problems are also discussed.

Description: ASTROS (Automated Structural Optimization System) is a finite-element-based multidisciplinary structural optimization procedure developed under Air Force sponsorship to perform automated preliminary structural design. The design task is the determination of the structural sizes that provide an optimal structure while satisfying numerous constraints from many disciplines. In addition to its automated design features, ASTROS provides a general transient and frequency response capability, as well as a special feature to perform a transient analysis of a vehicle subjected to a nuclear blast. The motivation for the development of a single multidisciplinary design tool is that such a tool can provide improved structural designs in less time than is currently needed. The role of such a tool is even more apparent as modern materials come into widespread use. Balancing conflicting requirements for the structure's strength and stiffness while exploiting the benefits of material anisotropy is perhaps an impossible task without assistance from an automated design tool. Finally, the use of a single tool can bring the design task into better focus among design team members, thereby improving their insight into the overall task.

Description: Large space structures (LSSs) and other dynamical systems of current interest are often extremely complex assemblies of rigid and flexible bodies subjected to kinematical constraints. A formulation is presented for the governing equations of constrained multibody systems via the application of singular value decomposition (SVD). The resulting equations of motion are shown to be of minimum dimension.

Description: A recently developed finite-element capability for general nonlinear shell analysis, featuring the use of three-dimensional constitutive equations within an efficient resultant-oriented framework, is employed to simulate the postbuckling response of an axially compressed composite cylindrical panel with a circular cutout. The problem is a generic example of modern composite aircraft components for which postbuckling strength (i.e., fail-safety) is desired in the presence of local discontinuities such as holes and cracked stiffeners. While the computational software does a reasonable job of predicting both the buckling load and the qualitative aspects of postbuckling (compared both with experiment and another code) there are some discrepancies due to: (1) uncertainties in the nominal layer material properties, (2) structural sensitivity to initial imperfections, and (3) the neglect of dynamic and local material delamination effects in the numerical model. Corresponding refinements are suggested for the realistic continuation of this type of analysis.

Description: The development of commercial finite element software is addressed. This software provides practical tools that are used in an astonishingly wide range of engineering applications that include critical aspects of the safety evaluation of nuclear power plants or of heavily loaded offshore structures in the hostile environments of the North Sea or the Arctic, major design activities associated with the development of airframes for high strength and minimum weight, thermal analysis of electronic components, and the design of sports equipment. In the more advanced application areas, the effectiveness of the product depends critically on the quality of the mechanics and mechanics related algorithms that are implemented. Algorithmic robustness is of primary concern. Those methods that should be chosen will maximize reliability with minimal understanding on the part of the user. Computational efficiency is also important because there are always limited resources, and hence problems that are too time consuming or costly. Finally, some areas where research work will provide new methods and improvements is discussed.