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: Microgravity offers a unique environment for studying polymer diffusion and polymer polymerization reactions. The absence of convection currents, which are the major mode of mixing at the molecular level on Earth, are eliminated or reduced in the microgravity environment. More importantly, the prediction of unique copolymer composition development in microgravity allows controlled formation of new compositions of matter. The absence of mixing at the molecular level should produce unique short block copolymers available for the first time for comonomer compositions which normally lead to random or long block copolymer under good mixing. The investigation of fundamental polymer diffusion and polymer polymerization processes in microgravity is proposed. This effort will involve fundamental studies of monomer and polymer diffusion; their effects on initiation, propagation, and especially termination kinetics rate constant; and the accurate evaluation of copolymerization reactivity ratios in microgravity. The experimental design is presented for these studies along with an evaluation technique for in situ monitoring of polymer diffusion and polymerization kinetics.

Description: The analysis of light scattered from an ensemble of particles has long been a preferred method for characterizing their physical properties. Instrumentation to perform the measurements which forms the basis for such analysis is available in many forms based upon a variety of different experimental techniques. A system is presented which is singularly applicable for making many types of measurements in a microgravity environment. The commercial version of this device, the DAWN-F, has been used in many labs throughout the world to perform analyses of particular importance for both research and production. Light scattering theory is reviewed and the structure and function of the system is described.

Description: Avalanche photodiodes (APD) are excellent small area, solid state detectors for photon counting. Performance possibilities include: photon detection efficiency in excess of 50 percent; wavelength response from 400 to 1000 nm; count rate to 10 (exp 7) counts per sec; afterpulsing at negligible levels; timing resolution better than 1 ns. Unfortunately, these performance levels are not simultaneously available in a single detector amplifier configuration. By considering theoretical performance predictions and previous and new measurements of APD performance, the anticipated performance of a range of proposed APD-based photon counting modules is derived.

Description: Traditional optical systems for photon correlation spectroscopy and laser anemometry have relied upon physically large and fairly expensive lasers, bulk-optics such as lenses of a few inches diameter, large mechanical mounts and carefully selected, fragile and bulky photon counting photomultiplier detectors. In some cases, experimental fluid dynamics at a desired position in a flow, perhaps deep inside complex machinery, is physically impossible or very difficult. Similar problems exist with photon correlation spectroscopy, e.g., remote and heterodyne experiments. Various optical and electro optical components were investigated and characterized with the aim of replacing existing photon correlation laser spectroscopy and anemometry techniques in miniaturized form, and with significant cost reduction. Very recently, a range of miniature, modular light scattering systems were constructed from little solid state optical and electro optical components, and experimentally verified measurement performance comparable to standard lab photon correlation spectroscopy and laser anemometry equipment.

Description: Images from an airborne, scanning radiometer operating at a frequency of 98 GHz have been analyzed. The millimeter-wave images were obtained in 1985-1986 using the JPL millimeter-wave imaging sensor. The goal of this study was to enhance the information content of these images and make their interpretation easier. A visual-interpretative approach was used for information extraction from the images. This included application of nonlinear transform techniques for noise reduction and for color, contrast, and edge enhancement. Results of using the techniques on selected millimeter-wave images are discussed.

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: It is shown how a commercial time interval counter can be used to measure the relative stability of two signals that are offset in frequency and mixed down to a beat note of about 1 Hz. To avoid the dead-time problem, the counter is set up to read the time interval between each beat note upcrossing and the next pulse of a 10 Hz reference pulse train. The actual upcrossing times are recovered by a simple algorithm whose outputs can be used for computing residuals and Allan variance. A noise floor-test yielded a delta f-f Allan deviation of 1.3 times 10 to the minus 9th power/tau relative to the beat frequency.

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 correction factor to the number density measured by the Forward Scattering Spectrometer Probe (FSSP) which compensates for dead time and coincidence errors was determined by calculating the probabilities of, and the average number of particles in, the six possible types of dead time and coincidence events. These probabilities and averages were calculated by means of a probabilistic model based on Poisson statistics. A Monte Carlo computer simulation of the FSSP operation was also carried out and the number density correction factor was compared with the Monte Carlo data. For an actual number density of 2000/cu cm, it was found that the measured number density was of the order of 300/cu cm.

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: Measurement of times of flight of sound waves can be used to determine temperatures in a gas. This paper describes a system, based on this principle, that is capable of giving the temperature profile in a nonisothermal gas volume, for example, prevalent in a large furnace. The apparatus is simple, rugged, accurate, and capable of being automated for process control applications. It is basically an acoustic waveguide where the outside temperature profile is transferred to a chosen gas contained inside the guide.

Description: The development of a unique noncontact temperature measurement device utilizing rotating analyzer ellipsometry is described. The technique circumvents the necessity of spectral emissivity estimation by direct measurement concomitant with radiance brightness. Simultaneous determinations of dielectric constants and refractive indices allow changes in the physical and chemical state of a heated surface to be monitored. The results of optical property measurements at 633 nm as functions of temperature between 1000 and 2500 K for eight transition metals including Hf, Ir, Mo, Nb, Pd, Pt, Ta, and V are presented together with preliminary results of oxidation studies on iridium.

Description: A design concept developed for a polarimeter on the vector magnetograph of the SAMEX satellite that would be very sensitive to solar vector magnetic fields is described. A description of the Poincare sphere is presented, along with the instrument scientific requirements, to provide an understanding of how the polarimeter design has been selected. It is shown that the design goal of a polarimetric sensitivity of 0.0001 can be achieved in the linear measurements using a hybrid analyzer. It is also noted that the systematic errors that produce linear crosstalk into the circular measurement will require the use of the redundant polarization measurements for in-flight calibrations and corrections of the data.

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

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

Description: The present analysis method for hot-wire data in supersonic turbulence takes sound field effects into account and yields greater accuracy in its treatment of flow variable fluctuations than existing methods despite requiring only a moderately accurate estimate of static pressure fluctuations. The method demonstrates the way in which neglecting pressure fluctuations will affect hot-wire data analysis, as well as indicating the probable direction the errors will take.

Description: Displays of multifrequency passive microwave data from the Special Sensor Microwave/Imager (SSM/I) flying on the Defense Meteorological Satellite Program (DMSP) spacecraft are presented. Observed brightness temperatures at 85.5 GHz (vertical and horizontal polarizations) and 37 GHz (vertical polarization) are respectively used to 'drive' the red, green, and blue 'guns' of a color monitor. The resultant false-color images can be used to distinguish land from water, highlight precipitation processes and structure over both land and water, and detail variations in other surfaces such as deserts, snow cover, and sea ice. The observations at 85.5 GHz also add a previously unavailable frequency to the problem of rainfall estimation from space. Examples of mesoscale squall lines, tropical and extra-tropical storms, and larger-scale land and atmospheric features as 'viewed' by the SSM/I are shown.

Description: A postdispersion system for astronomical observations with Fourier transform spectrometers in the thermal infrared has been developed which improves the sensitivity of radiation noise limited observations by reducing the spectral range incident on the detector. Special attention is given to the first-generation blocked impurity band detector. Planetary, solar, and stellar observations are reported.

Description: The measurement by a three-dimensional laser Doppler velocimeter of a turbulent flow has been numerically simulated. Errors associated with the probe volume geometry and the coincidence time window concept are revealed. One type of error occurs for high system data rates when multiple particles lead to system realizations. Another error occurs associated with a geometric bias discovered in the present study. This three-dimensional ldv geometric bias exists even for single-particle realizations and regardless of the system data rate. A technique for the elimination of the geometric bias is presented.

Description: A computer-controlled tunable diode laser spectrometer and spectral analysis software are described. The three-channel system records simultaneously the transmission of a subject gas, a temperature-stabilized etalon, and a calibration gas. The software routines are applied to diode laser spectra of HNO3 and NO2 to illustrate the procedures adopted for conversion of raw spectral data to useful transmission and harmonic spectra. Extraction of line positions, absorption intensities, collisional broadening coefficients, and gas concentrations from recorded spectra is also described.

Description: Optical interconnects are being considered for control signal distribution in phased array antennas. A packaged hybrid GaAs optical controller with a 1:16 demultiplexed output that is suitable for this application is described. The controller, which was fabricated using enhancement/depletion mode MESFET technology, operates at demultiplexer-limited input data rates up to 305 Mb/s and requires less than 200 microW optical input power.

Description: A technique has been developed to support the study of the effects of cosmic rays on integrated circuits. The system is designed to determine the particle distribution across the surface of an integrated circuit accurately while the circuit is bombarded by a particle beam. The system uses photomultiplier tubes, an octal discriminator, a computer-controlled NIM quad counter, and an IBM PC. It provides real-time operator feedback for fast beam tuning and monitors momentary fluctuations in the particle beam. The hardware, software, and system performance are described.

Description: A three-channel transportable radiometer that operates at 20.6, 31.65, and 90.0 GHz has been developed. The design of the radiometer is described and the results of various experimental applications of the radiometer are presented. The experiments include a comparison of brightness temperature measurements at various frequencies with the water vapor absorption models of Waters (1976) and Lieve (1989). Measurements of absorption ratios for cloud attenuation are given and the radiometer measurements are compared with oxygen absorption calculated using the Rasenkranz (1988) model.

Description: This paper describes a program of airborne radiometric imaging at 90 GHz and 140 GHz. Using high sensitivity (below 1 K) and high angular resolution (0.5-1.0 degree), high quality images have been made. The following measurements are discussed: cloud and fog penetration at 90 GHz, discrimination between agricultural and urban areas, discrimination between different vegetation types, detection of vehicles on roads, detection and classification of airports and airplanes, ship detection and quantitative oil spill sensing. The application of information enhancement techniques with automatic and real time application aspects is also described, and results of applied techniques for contrast and contour enhancement are shown.

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: Mercuric iodide X-ray detectors have been undergoing tests in a prototype scanning electron microscope system being developed for unmanned space flight. The detector program addresses the issues of geometric configuration in the SEM, compact packaging that includes separate thermoelectric coolers for the detector and FET, X-ray transparent hermetic encapsulation and electrical contacts, and a clean vacuum environment.

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 recent development of Scanning Tunneling Microscopy technology allows the application of electron tunneling to position detectors for the first time. The vacuum tunnel junction is one of the most sensitive position detection mechanisms available. It is also compact, simple, and requires little power. A prototype accelerometer based on electron tunneling, and other sensor applications of this promising new technology are described.

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

Description: The 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: New methods have been developed to implement position sensors based on electron tunneling. The electron tunneling methods enable position to be detected with sub-Angstrom resolution using a compact mechanical structure and simple electronic control elements. A prototype accelerometer is the first sensor based on these principles; it shows reliable operation with a noise-limited sensitivity and a bandwidth of approximately 3 kHz. Based on these results, it is expected that tunnel sensors optimized for many applications will enable an entirely new class of sensors to be developed.

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 Particle Measuring Systems (PMS) Forwared Scattering Spectrometer Probe (FSSP-100) which was flown on the ER-2 during the Airborne Antarctic Ozone Experiment for the measurement of particles in polar stratospheric clouds has been evaluated and calibrated. The sample volume of the probe per 10-s sampling period increases from 257 cu cm for 1-micron particles to 412 for 15-micron particles, but there is substantial uncertainty in this value. Limitations in the measurements from this instrument and possible corrections are discussed. The uncertainty in the total particle mass measured by the probe may be as large as + or - 100 percent. Recommendations are given for the processing of data from the FSSP used in this project.

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: A recent study by Rocketdyne for NASA identified laser anemometry, using a compact optical head, as a feasible diagnostic instrument for the Space Shuttle Main Engine (SSME) Model Verification experiments. Physical Research, Inc. (PRI) is presently under contract from NASA Lewis to develop and deliver such a laser anemometer system. For this application, it is desired to place the laser at a remote distance from the engine, and use single mode polarization preserving fiber optics for the transmission of the laser light to and from the measurement head. Other requirements are given. Analytical and experimental tools are being used to develop the technologies required for the laser anemometer. These include finite element analysis of the optical head and vibration tests for various optical and mechanical components. Design of the optical head and the fiber optic connectors are driven by the temperature and vibration requirements for the measurement environment. Results of the finite element analysis and the vibration tests of the components are included. Conceptual design of the fiber optic launcher and the optical probe has also been complete. Detailed design of the probe as well as the fabrication and assembly of the components is in progress.

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: Spacecraft range measurements have provided the most accurate tests, to date, of some relativistic gravitational parameters, even though the measurements were made with ranging systems having error budgets of about 10 meters. Technology is now available to allow an improvement of two orders of magnitude in the accuracy of spacecraft ranging. The largest gains in accuracy result from the replacement of unstable analog components with high speed digital circuits having precisely known delays and phase shifts.

Description: EDO Corporation, Barnes Engineering Division designed and constructed a high resolution thermal imaging system on contract to Lockheed for use in the SDI Star Lab. This employs a Pt Si CCD array which is sensitive in the spectral range of 3 to 5 microns. Star Lab will be flown in the Shuttle bay and consists basically of a large, reflecting, tracking telescope with associated sensors and electronics. The thermal imaging system is designed to operate in the focal plane of this telescope. The configuration of the system is illustrated. The telescope provides a collimated beam output which is focussed onto the detector array by a silicon objective lens. The detector array subtends a field of view of 1.6 degrees x 1.22 degrees. A beam switching mirror permits bypassing the large telescope to give a field of 4 degrees x 3 degrees. Two 8 position filter wheels are provided, and background radiation is minimized by Narcissus mirrors. The detector is cooled with a Joule-Thompson cryostat fed from a high pressure supply tank. This was selected instead of a more convenient closed-cycle system because of concern with vibration. The latter may couple into the extremely critical Starlab tracking telescope. The electronics produce a digitized video signal for recording. Offset and responsivity correction factors are stored for all pixels and these corrections are made to the digitized output in real time.

Description: A preliminary evaluation was made by ORNL of a two-color ratio pyrometer (TCRP) for temperature control in the Modular Electromagnetic Levitation (MEL) experiment. A discussion was presented by Eric Spjut at the 1987 NASA Non-Contact Temperature Measurement Workshop (NASA Conf. Publ. 2503, pp. 182-213) in which he described the non-linear characteristics of the time response of TCPs. Researchers replicated his model and results and note that the non-linear response behavior is minimized for small temperature steps at high temperatures. They then used the predicted response in a model for a proportional or integral feedback controller and predicted the control characteristics for heating and cooling a 5-mm diameter sphere of niobium at high (1500 to 2750 K) temperatures. The analysis shows that for a slow (25-ms) time response for a commercial RCRP, overshoots of several hundred kelvins will result from a 100-K decrease in the setpoint, and temperature tracking errors of 14 to 45 K will occur for control temperature ramps of 1000K/s. For a fast (greater than 0.1 ms) time response, the overshoot and ramp response errors are largely eliminated.

Description: The development of a noncontact temperature measurement device utilizing rotating analyzer ellipsometry is described. The technique circumvents the necessity of spectral emissivity estimation by direct measurement concomittant with radiance brightness. Using this approach, the optical properties of electromagnetically levitated liquid metals Cu, Ag, Au, Ni, Pd, Pt, and Zr were measured in situ at four wavelengths and up to 600 K superheat in the liquid. The data suggest an increase in the emissivity of the liquid compared with the incandescent solid. The data also show moderate temperature dependence of the spectral emissivity. A few measurements of the optical properties of undercooled liquid metals were also conducted. The data for both solids and liquids show excellent agreement with available values in the literature for the spectral emissivities as well as the optical constants.

Description: Commercially available types of infrared thermal imaging instruments, both viewers (qualitative) and imagers (quantitative) are discussed. The various scanning methods by which thermal images (thermograms) are generated will be reviewed. The performance parameters (figures of merit) that define the quality of performance of infrared radiation thermometers will be introduced. A discussion of how these parameters are extended and adapted to define the performance of thermal imaging instruments will be provided. Finally, the significance of each of the key performance parameters of thermal imaging instruments will be reviewed and procedures currently used for testing to verify performance will be outlined.

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: The optimal placement of discrete actuators and sensors is posed as a combinatorial optimization problem. Two examples for truss structures were used for illustration; the first dealt with the optimal placement of passive dampers along existing truss members, and the second dealt with the optimal placement of a combination of a set of actuators and a set of sensors. Except for the simplest problems, an exact solution by enumeration involves a very large number of function evaluations, and is therefore computationally intractable. By contrast, the simulated annealing heuristic involves far fewer evaluations and is best suited for the class of problems considered. As an optimization tool, the effectiveness of the algorithm is enhanced by introducing a number of rules that incorporate knowledge about the physical behavior of the problem. Some of the suggested rules are necessarily problem dependent.

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 sensor has been designed and tested for precise pointing applications. The device is able to sense extremely small rotary motion and is immune to cross-axis forces. The hardware and design characteristics of the torque sensor are presented. Test data, integrated control methodology, and future applications are included.

Description: The development of a prototype Control Moment Gyroscope (CMG) is discussed. Physical characteristics and the results of functional testing are presented to demonstrate the level of system performance obtained. Particular attention is given to how the man-rated mission requirement influenced the choice of the materials, fabrication, and design details employed. Comparisons are made of the measured system responses against the prediction generated by computer simulation.

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.