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

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

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

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

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

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

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

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

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

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

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

Description: 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: 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.

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.

Description: An analytical method is being developed to determine the signature of an acoustic emission waveform from a growing crack and the results of this analysis are compared to experimentally obtained values. Within the assumptions of linear elastic fracture mechanics, a two dimensional model is developed to examine a semi-infinite crack that, after propagating with a constant velocity, suddenly stops. The analytical model employs an integral equation method for the analysis of problems of dynamic fracture mechanics. The experimental procedure uses an interferometric apparatus that makes very localized absolute measurements with very high fidelity and without acoustically loading the specimen.

Description: A computer program was developed to take a NASTRAN output file from a normal modes analysis and calculate the modal strain energies of selected elements. The FORTRAN program can determine the modal strain energies for CROD, CBAR, CELAS, CTRMEM, CQDMEM2, and CSHEAR elements. Modal strain energies are useful in estimating damping in structures.

Description: Many finite element preprocessors describe finite element model geometry with points, lines, surfaces and volumes. One method for describing these basic geometric entities is by use of parametric cubics which are useful for representing complex shapes. The lines, surfaces and volumes may be discretized for follow on finite element analysis. The ability to limit or selectively recover results from the finite element model is extremely important to the analyst. Equally important is the ability to easily apply boundary conditions. Although graphical preprocessors have made these tasks easier, model complexity may not lend itself to easily identify a group of grid points desired for data recovery or application of constraints. A methodology is presented which makes use of the assignment of grid point locations in parametric coordinates. The parametric coordinates provide a convenient ordering of the grid point locations and a method for retrieving the grid point ID's from the parent geometry. The selected grid points may then be used for the generation of the appropriate set and constraint cards.

Description: The DMAP coding was automated to such an extent by using the device of bubble vectors, that it is useable for analyses in its present form. This feasibility study demonstrates that the Ritz Method is so compelling as to warrant coding its modules in FORTRAN and organizing the resulting coding into a new Rigid Format. Even though this Ritz technique was developed for unsymmetric stiffness matrices, it offers advantages to problems with symmetric stiffnesses. If used for the symmetric case the solution would be simplified to one set of modes, because the adjoint would be the same as the primary. Its advantage in either type of symmetry over a classical eigenvalue modal expansion is that information density per Ritz mode is far richer than per eigenvalue mode; thus far fewer modes would be needed for the same accuracy and every mode would actively participate in the response. Considerable economy can be realized in adapting Ritz vectors for modal solutions. This new Ritz capability now makes NASTRAN even more powerful than before.

Description: A monolithic spatial light modulator chip was developed consisting of a large number of micrometer-scale mirror cells which can be rotated through an angle by application of an electrostatic field. The field is generated by electronics integral to the chip. The chip has application in photoreceptor based non-impact printing technologies. Chips containing over 16000 cells were fabricated, and were tested to several billions of cycles. Finite Element Analysis (FEA) of the device was used to model both the electrical and mechanical characteristics.

Description: The identification of power flow paths in dynamically loaded structures is an important, but currently unavailable, capability for the finite element analyst. For this reason, methods for calculating power flows and mechanical intensities in finite element models are developed here. Formulations for calculating input and output powers, power flows, mechanical intensities, and power dissipations for beam, plate, and solid element types are derived. NASTRAN is used to calculate the required velocity, force, and stress results of an analysis, which a post-processor then uses to calculate power flow quantities. The SDRC I-deas Supertab module is used to view the final results. Test models include a simple truss and a beam-stiffened cantilever plate. Both test cases showed reasonable power flow fields over low to medium frequencies, with accurate power balances. Future work will include testing with more complex models, developing an interactive graphics program to view easily and efficiently the analysis results, applying shape optimization methods to the problem with power flow variables as design constraints, and adding the power flow capability to NASTRAN.

Description: Numerical techniques for calculating the low frequency vibrational resonances of submerged structures are reviewed. Both finite element and boundary element approaches for calculating fully-coupled added mass matrices for use in NASTRAN analysis are described and illustrated. The finite element approach is implemented using existing capability in NASTRAN. The boundary element approach uses the NASHUA structural-acoustics program to compute the added mass matrix. The two procedures are compared to each other for the case of a submerged cylindrical shell with flat end closures. It is concluded that both procedures are capable of computing accurate submerged resonances and that the more elegant boundary element procedure is easier to use but may be more expensive computationally.

Description: The CONTINUE feature in transient analysis as implemented in the standard release of COSMIC/NASTRAN has inherent errors associated with it. As a consequence, the results obtained by a CONTINUEd restart run do not, in general, match the results that would be obtained in a single run without the CONTINUE feature. These inherent errors were eliminated by improvements to the restart logic that were developed by RPK Corporation and that are available on all RPK-supported versions of COSMIC/NASTRAN. These improvements ensure that the results of a CONTINUEd transient analysis run are the same as those of a non-CONTINUEd run. In addition, the CONTINUE feature was extended to transient analysis involving uncoupled modal equations. The improvements and enhancement were illustrated by examples.

Description: A modification to the NASTRAN solution sequence for transient analysis with direct time integration (COSMIC NASTRAN rigid format 9) was developed and incorporated into a DMAP alter. This DMAP alter calculates the buckling stability of a dynamically loaded structure, and is used to predict the onset of structural buckling under stress-wave loading conditions. The modified solution sequence incorporates the linear buckling analysis capability (rigid format 5) of NASTRAN into the existing Transient solution rigid format in such a way as to provide a time dependent eigensolution which is used to assess the buckling stability of the structure as it responds to the impulsive load. As a demonstration of the validity of this modified solution procedure, the dynamic buckling of a prismatic bar subjected to an impulsive longitudinal compression is analyzed and compared to the known theoretical solution. In addition, a dynamic buckling analysis is performed for the analytically less tractable problem of the localized dynamic buckling of an initially flawed composite laminate under transverse impact loading. The addition of this DMAP alter to the transient solution sequence in NASTRAN facilitates the computational prediction of both the time at which the onset of dynamic buckling occurs in an impulsively loaded structure, and the dynamic buckling mode shapes of that structure.

Description: A powerful enhancement to the DMAP alter capability was developed and is available on all RPK-supported versions of COSMIC/NASTRAN. This enhancement involves the addition of two alter control cards, called INSERT and DELETE, to the Executive Control Deck. These cards allow for DMAP alters to be made by referencing DMAP statements by their module names rather than by their statement numbers in the rigid format DMAP sequence. This allows for increased user convenience and flexibility and makes alters more meaningful to the user. In addition, DMAP alter packages employing the alter control cards will be much less susceptible to future changes in rigid format DMAPs than alter packages employing the standard ALTER control cards. The usage of the cards is illustrated by examples.

Description: Major improvements were made to the IBM version of COSMIC/NASTRAN by RPK Corporation under contract to IBM Corporation. These improvements will become part of COSMIC's IBM version and will be available in the second quarter of 1989. The first improvement is the inclusion of code to take advantage of IBM's new Vector Facility (VF) on its 3090 machines. The remaining improvements are modifications that will benefit all users as a result of the extended addressing capability provided by the MVS/XA operating system. These improvements include the availability of an in-memory data base that potentially eliminates the need for I/O to the PRIxx disk files. Another improvement is the elimination of multiple load modules that have to be loaded for every link switch within NASTRAN. The last improvement allows for NASTRAN to execute above the 16 mega-byte line. This improvement allows for NASTRAN to have access to 2 giga-bytes of memory for open core and the in-memory data base.

Description: Changes in both software and hardware are rapidly bringing conceptual engineering tools like finite element analysis into mainstream mechanical design. Systems that integrate all phases of the manufacturing process provide the most cost benefits. The application of programming concepts like object oriented programming allow for the encapsulation of intelligent data within the design geometry. This combined with declining cost in per seat hardware bring new alternatives to the user.

Description: Demands for nonlinear time history simulations of large, flexible multibody dynamic systems has created a need for efficient interfaces between finite-element modeling programs and time-history simulations. One such interface, TREEFLX, an interface between NASTRAN and TREETOPS, a nonlinear dynamics and controls time history simulation for multibody structures, is presented and demonstrated via example using the proposed Space Station Mobile Remote Manipulator System (MRMS). The ability to run all three programs (NASTRAN, TREEFLX and TREETOPS), in addition to other programs used for controller design and model reduction (such as DMATLAB and TREESEL, both described), under a UNIX Workstation environment demonstrates the flexibility engineers now have in designing, developing and testing control systems for dynamically complex systems.

Description: Starting with the design objective the operational cycle life of the Swaging Tool was increased. To accomplish this increase in cycle life without increasing the size or weight of the tool would be engineering achievement. However, not only was the operational cycle life increased between 2 to 10 times but simultaneously the size and weight of the Swage Tool was decreased by about 50 percent. This accomplishment now becomes an outstanding engineering achievement. This achievement was only possible because of the computerized Patran, Nastran and Medusa programs.

Description: Patran and Supertab are interactive computer graphics pre- and postprocessors that can be used to generate NASTRAN bulk data decks and to visualize results from a NASTRAN analysis. Both of the programs are in use at the Numerical Structural Mechanics Branch of the David Taylor Research Center (DTRC). Various aspects of Patran and Supertab are discussed including: geometry modeling, finite element mesh generation, bulk data deck creation, results translation and visualization, and the user interface. Some advantages and disadvantages of both programs will be pointed out.

Description: Three areas of improvement in COSMIC/NASTRAN, 1989 release, were incorporated recently that make the analysis program run faster on large problems. Actual log files and actual timings on a few test samples that were run on IBM, CDC, VAX, and CRAY computers were compiled. The speed improvement is proportional to the problem size and number of continuation cards. Vectorizing certain operations in BANDIT, makes BANDIT run twice as fast in some large problems using structural elements with many node points. BANDIT is a built-in NASTRAN processor that optimizes the structural matrix bandwidth. The VAX matrix packing routine BLDPK was modified so that it is now packing a column of a matrix 3 to 9 times faster. The denser and bigger the matrix, the greater is the speed improvement. This improvement makes a host of routines and modules that involve matrix operation run significantly faster, and saves disc space for dense matrices. A UNIX version, converted from 1988 COSMIC/NASTRAN, was tested successfully on a Silicon Graphics computer using the UNIX V Operating System, with Berkeley 4.3 Extensions. The Utility Modules INPUTT5 and OUTPUT5 were expanded to handle table data, as well as matrices. Both INPUTT5 and OUTPUT5 are general input/output modules that read and write FORTRAN files with or without format. More user informative messages are echoed from PARAMR, PARAMD, and SCALAR modules to ensure proper data values and data types being handled. Two new Utility Modules, GINOFILE and DATABASE, were written for the 1989 release. Seven rigid elements are added to COSMIC/NASTRAN. They are: CRROD, CRBAR, CRTRPLT, CRBE1, CRBE2, CRBE3, and CRSPLINE.

Description: A metal-matrix composite (MMC) model was developed which includes the concept of damage evolution. The evolution of damage is assumed to be governed by a Kachanov-type equation. This viscoplastic damage model was implemented in the finite element code, MARC. Both uniaxial (creep) and multiaxial (an internally pressurized thick-walled cylinder) problems were analyzed using this implementation. Some preliminary results are presented which consider monotonic (constant) loadings. The creep curves including damage for four fiber orientations are presented. As expected, the minimum creep occurs when load is applied in a direction parallel to the fibers. The tangential strains at the inner radius of a thick-walled MMC-cylinder for four fiber orientations are shown with damage included. The cylinder exhibits the maximum creep resistance when the fibers are oriented in the circumferential direction, perpendicular to the axis of the cylinder. Time-to-failure for the thick-walled cylinder for the same fiber orientation angles is also shown. As expected, the life of the cylinder can be increased by orientating the fibers in the circumferential direction, perpendicular to the axis of the cylinder. The results, although qualitative, indicate that significant benefits in creep-resistance and service life can be achieved by using MMC materials as structural materials for high-temperature design.

Description: Alloy 718 crack growth experiments were conducted to assess the ability of the selected path-independent (P-I) integrals to describe the elevated temperature crack growth behavior. These tests were performed on single edge notch (SEN) specimens under displacement control with multiple extensometers to monitor the specimen and crack mouth opening displacement (CMOD). The displacements in these tests were sufficiently high to induce bulk cyclic inelastic deformation of the specimen. Under these conditions, the linear elastic fracture mechanics (LEFM) parameter K does not correlate the crack growth data. The experimentally measured displacement gradients at the end of specimen gage length were used as the boundary conditions in elastic-plastic finite element method (FEM) analyses. These analyses were performed with a node release approach using CYANIDE, a GEAE FEM code, which included a gap element which is capable of efficiently simulating crack closure. Excellent correlation was obtained between the experimentally measured and predicted variation of stress and CMOD with crack length and the stress-CMOD loops for Alloy 718 tests conducted at 538 C. This confirmed the accuracy of the FEM crack growth simulation approach. The experimentally measured crack growth rate data correlated well the selected P-I integrals. These investigations have produced significant progress in developing P-I integrals as non-linear fracture mechanics parameters. The results suggest that this methodology has the potential of accurately describing elevated temperature crack growth behavior under the combined influence of thermal cycling and bulk elastic-inelastic deformation states.

Description: Infrared spectroscopy of warm (about 150 to 750 K), dusty astronomical sources has revealed a structured emission spectrum which can be diagnostic of the composition, temperature, and in some cases, even size and shape of the grains giving rise to the observed emission. The identifications of silicate emission in oxygen rich objects and SiC in carbon rich object are two examples of this type of analysis. Cometary spectra at moderate resolution have similarly revealed silicate emission, tying together interstellar and interplanetary dust. However, Goebel has pointed out that some astronomical sources appear to contain a different type of dust which results in a qualitatively different spectral shape in the 8 to 13 micron region. The spectra shown make it appear unlikely that silicon nitride can be identified as the source of the 8 to 13 micron emission in either NGC 6572 or Nova Aql 1982. The similarity between the general wavelength and shape of the 10 micron emission from some silicates and that from the two forms of silicon nitride reported could allow a mix of cosmic grains which include some silicon nitride if only the 8 to 13 micron data are considered.

Description: The x ray generation due to arising of hot dense plasma balls at high-velocity (greater than or equal to 70 km/s) collisions of dust grains in the interstellar medium is considered. Analytical expressions for efficiency of conversions of colliding dust particle kinetic energy into x ray radiation are presented. The observed intensity distribution of the diffuse component of soft cosmic x rays (0.1 to 1 keV) may be partly caused by collisions between the dusty components of high velocity clouds and of the disk of the Galaxy.

Description: Until now it has remained difficult to account for the rather high temperatures seen in many diffuse interstellar clouds. Various heating mechanisms have been considered: photoionization of minor species, ionization of H by cosmic rays, and photoelectric effect on small grains. Yet all these processes are either too weak or efficient under too restricting conditions to balance the observed cooling rates. A major heat source is thus still missing in the thermal balance of the diffuse gas. Using photoionization cross sections measured in the lab, it was shown that in order to balance the observed cooling rates in cold diffuse clouds (T approx. 80 K) the PAHs would have to contain 15 percent of the cosmic abundance of carbon. This value does not contradict the former estimation of 6 percent deduced from the IR emission bands since this latter is to be taken as a lower limit. Further, it was estimated that the contribution to the heating rate due to PAH's in a warm HI cloud, assuming the same PAH abundance as for a cold HI cloud, would represent a significant fraction of the value required to keep the medium in thermal balance. Thus, photoionization of PAHs might well be a major heat source for the cold and warm HI media.

Description: Polycyclic aromatic hydrocarbon (PAH) molecules have recently been proposed as an important and hitherto undetected component of the Interstellar Medium (ISM). The theory was based on an explanation of the Unidentified IR Emission Bands by Leger et al. It has already led to a verified prediction on extended galactic and extragalactic emissions measured by IRAS, or by a recent balloon borne experiment. The physics that rules the motion of such molecules in the ISM was studied, taking into account their coupling with the ambient gas, the radiation field (absorption and emission) and the static magnetic field. This is important for many implications of the PAH theory such as the radio emission by these molecules or the expected polarization of their IR emission. A reflection nebulae is considered where the situation is rather well known. Every day life of a mean PAH molecule in such a region is as follows: every 3 hrs a UV photon is absorbed heating the molecule to a thousand degs; the temperature decay due to cooling by IR emission follows then within a few seconds. A collision with a molecule of gas occurs typically once a week, while an H atom is ejected or captured at the same rate. A typical cooling cycle after a heat impulse is given. The PAH molecules studied as representative of the family has typically 50 atoms, a radius of 4.5 A, is circular and has a molecular mass of M = 300; its permanent dipole moment is 3 Debye.

Description: The emission lines observed in many interstellar IR sources at 3.28, 6.2, 7.7, 8.7, and 11.3 microns are theorized to originate from polycyclic aromatic hydrocarbons (PAHs). These assignments are based on analyses of lab IR spectra of neutral PAHs. However, it is likely that in the interstellar medium that PAHs are ionized, i.e., are positively charged. Besides, as pointed out by Allamandola et al., although the IR emission band spectrum resembles what one might expect from a mixture of PAHs, it does not match in details such as frequency, band profile, or relative intensities predicted from the absorption spectra of any known PAH molecule. One source of more information to test the PAH theory is ab initio molecular orbital theory. It can be used to compute, from first principles, the geometries, vibrational frequencies, and vibrational intensities for model PAH compounds which are difficult to study in the lab. The Gaussian 86 computer program was used to determine the effect of ionization on the infrared absorption spectra of several small PAHs: naphthalene and anthracene. A preliminary report is presented of the results of these calculations.

Description: Several celestial objects, including UV rich regions of planetary and reflection nebulae, stars, H II regions, and extragalactic sources, are characterized by the unidentified infrared emission bands (UIR bands). A few years ago, it was proposed that polycyclic aromatic hydrocarbon species (PAHs) are responsible for most of the UIR bands. This hypothesis is based on a spectrum analysis of the observed features. Comparisons of observed IR spectra with lab absorption spectra of PAHs support the PAH hypothesis. An example spectrum is represented, where the Orion Bar 3.3 micron spectrum is compared with the absorption frequencies of the PAHs Chrysene, Pyrene, and Coronene. The laser excited 3.3 micron emission spectrum is presented from a gas phase PAH (azulen). The infrared fluorescence theory (IRF) is briefly explained, followed by a description of the experimental apparatus, a report of the results, and discussion.

Description: It is well known that a large number of celestial objects exhibit, in the range 3 to 12 micron, a family of emission features called unidentified infrared bands (UIR). They usually appear together and are associated with UV sources. Recently various authors have suggested that these features could be attributed to solid carbonaceous materials. Following this interest, a systematic analysis was performed of various types of amorphous carbon grains and polycyclic aromatic hydrocarbons (PAH), produced in lab. Updating results of Raman measurements performed on several carbonaceous materials, chosen according to their astrophysical interest, are presented. The measurements were made by means of a Jobin-Yvon monochromator HG2S and standard DC electronic. The line at 5145 A of an Ar+ laser was used as excitation source.

Description: The discovery of broad emission features between 3.2 and 3.6 microns were reported in the spectrum of Nova Cen 1986 (V842 Cen) some 300 days following outburst and remaining prominent for several months. The general characteristics of these features are similar to those attributed to polycyclic hydrocarbon (PAH) molecules in other dusty sources, although the relative strengths are different, and these observations provide the first clear evidence for molecular constituents other than graphite particles in the ejecta of novae.