ALBERT

Description: Viewgraphs on current analytical gaps for future needs for large space systems are presented. Topics covered include: future spacecraft; common control objectives; accuracy requirements; increasing complexity; promising approaches; and analytical gaps.

Description: Large scale space missions of the near future will depend upon successful multi-launch coordination and construction in the space environment. One of the main challenges is how to accomplish a valid global analysis of a construction project with the intent of improving safety, reducing overall mission cost, and total construction time. These three items are dependent on the interruptability of the project, which is the ability of the project to recover from unplanned interruptions; such as failure of the launch vehicle; sudden, on-orbit, crew illness; or damage from a space debris impact on the partially completed space structure. A new method for addressing and analyzing this type of problem is being developed. The method is called Program Interruptability and Risk Evaluation Technique, or PIRET. PIRET has been developed in order to model and analyze potential interruptability concerns of the construction of the U.S. Space Station Freedom (SSF), although PIRET is applicable to any complex, multi-launch structural assembly. This paper is a progress report on the continuing research of the NASA Center for Space Construction at the University of Colorado, Boulder into this area of space construction interruptability. The paper will define the problem of interruptability, will diagram the PIRET approach to space construction, will share results from a preliminary PIRET analysis of SSF, and will show that PIRET is a useful tool for modelling space construction interruptability.

Description: Large Space Structures do not have much damping, which necessitates the installation of a controller onto the structure. If the controller is improperly designed, the structure may become unstable and be destroyed. Since Large Space Structures are extremely expensive pieces of hardware, new controllers must not be tested first on the structure. They must first be tested in computer simulations. Until now, the usual procedure for simulating controlled Large Space Structures is to compute a reduced order modal representation of the structure and then apply the controller. However, this procedure entails modal truncation error. A new software package which is free from this error is currently under development within the Center for Space Construction. The more accurate finite element representation of the structure is used in the simulation, instead of the less accurate reduced order modal representation. This software also features an efficient matrix storage scheme, which effectively deals with the asymmetric system matrices which occur when control is added to the structure. Also, an integration algorithm was chosen so that the simulation is a reliable indicator of system stability or instability. The software package is fairly general in nature. Linearity of the finite element model and of the controller is the only assumption made. Actuator dynamics, sensor dynamics, noise, and disturbances can be handled by the package. In addition, output feedback of displacement, velocity, and/or acceleration signals can be simulated. Kalman state estimation was also implemented. This software was tested on a finite element model of a real Large Space Structure: The Mini-Mast Truss. Mini-Mast is a testbed at NASA-Langley which is currently under development. A 714 degree of freedom finite element model was computed, and a 19 state controller was designed for it. Torque wheel dynamics were added to the model, and the entire closed loop system was simulated with the software package.

Description: Guided wave modes in cladded or uncladded fiber-reinforced composite plates and tubes have been analyzed using a stiffness method in which the displacement variation through the thickness is approximated by polynomial interpolation functions. This allows for an arbitrary number of laminations and fiber orientations different from lamina to lamina. It is shown that dispersive behavior of guided modes depends significantly on the cladding, number of laminae, and interfaces between the adjacent laminae. A hybrid modeling technique is described in which an inner region containing cracks (or other defects) is discretized by finite elements, and the field in the exterior region is represented in terms of modes that are found using the stiffness method described above. It is found that the reflected and transmitted amplitudes of modes vary significantly with the size of a transverse or longitudinal (delamination) crack and frequency. We have also studied the impact response of a unidirectional fiber-reinforced plate. Received signals at the epicentral and other locations are shown. Strong longitudinal anisotropy of the graphite/epoxy plate causes the signal to be considerably different from that in an isotropic plate.

Description: One of the most difficult problems in controlling large systems and structures is compensating for the destructive interaction which can occur between the reduced-order model (ROM) of the plant, which is used by the controller, and the unmodeled dynamics of the plant, often called the residual modes. The problem is more significant in the case of large space structures because their naturally light damping and high performance requirements lead to more frequent, destructive residual mode interaction (RMI). Using the design/compensation technique of residual mode filters (RMF's), effective compensation of RMI can be accomplished in a straightforward manner when using linear controllers. The use of RMF's has been shown to be effective for a variety of large structures, including a space-based laser and infinite dimensional systems. However, the dynamics of space structures is often uncertain and may even change over time due to on-orbit erosion from space debris and corrosive chemicals in the upper atmosphere. In this case, adaptive control can be extremely beneficial in meeting the performance requirements of the structure. Adaptive control for large structures is also based on ROM's and so destructive RMI may occur. Unfortunately, adaptive control is inherently nonlinear, and therefore the known results of RMF's cannot be applied. The purpose is to present the results of new research showing the effects of RMI when using adaptive control and the work which will hopefully lead to RMF compensation of this problem.

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: The objective of this program is to develop technology needed for structural evaluation of alternative space construction concepts. Those concepts are as follows: interactive effects on dynamic performance of various environment and self-generated disturbance; new materials concepts, failure mechanisms, and non-destructive evaluation/failure detection; develop stable control algorithms and design effective combination of hierarchical and adaptive controls; and assess control-structure integrated performance and stability.

Description: Partially constructed/assembled structures in space are complicated enough but their dynamics will also be operating in closed-loop with feedback controllers. The dynamics of such structures are modeled by large-scale finite element models. The model dimension L is extremely large (approximately 10,000) while the numbers of actuators (M) and sensors (P) are small. The model parameters M(sub m) mass matrix, D(sub o) damping matrix, and K(sub o) stiffness matrix, are all symmetric and sparse (banded). Thus simulation of open-loop structure models of very large dimension can be accomplished by special integration techniques for sparse matrices. The problem of simulation of closed-loop control of such structures is complicated by the addition of controllers. Simulation of closed-loop controlled structures is an essential part of the controller design and evaluation process. Current research in the following areas is presented: high-order simulation of actively controlled aerospace structures; low-order controller design and SCI compensation for unmodeled dynamics; prediction of closed-loop stability using asymptotic eigenvalue series; and flexible robot manipulator control experiment.

Description: Within the Center for Space Construction, the SIMSTRUC project's objectives center around the development of simulation tools for the realistic analysis of large space structures. The word 'tools' is the broad sense; it designates mathematical models, finite element/finite difference formulations, computational algorithms, implementations on advanced computer architectures, and visualization capabilities. The results of our activities during the first year within the SIMSTRUC project are reported. On the modeling side, an alternative approach to fluid/thermal/structure interaction analysis that is a departure from the 'loosely coupled' and 'unified' approaches that are being currently practiced are described. The advantages of our approach both in terms of accuracy and computational efficiency were demonstrated. On the computational side, a software architecture for parallel/vector and massively parallel supercomputers that speeds up finite element and finite difference computations by several orders of magnitude is presented. As an example, the simulation of the deployment of a space structure that used to require over six hours of a workstation using a conventional finite element software, now runs on a multiprocessor using a parallel computation strategy in less than three seconds. In order to promote the physical understanding of the simulation behavior, a real-time visualization capability on the Connection Machine, which allows the analyst to watch the graphical animation of the results at the same time these are generated, was also developed. It is believed that by combining efficient analytical formulations with the state-of-the-art high performance computer implementations and superfast visualization capabilities, SIMSTRUC is moving fast towards the real-time simulation of large space structures. The designers as well as the researchers will certainly benefit from this technology.

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: Future spacecraft design will be affected by collisions with man-made debris orbiting the earth. Most of this orbital space debris comes from spent rocket stages. It is projected that the source of future debris will be the result of fragmentation of large objects through hypervelocity collisions. Orbiting spacecraft will have to be protected from hypervelocity debris in orbit. The options are to armor the spacecraft, resulting in increased mass, or actively removing the debris from orbit. An active space debris sweeper is described which will utilize momentum transfer to the debris through laser-induced ablation to alter its orbital parameters to reduce orbital lifetime with eventual entry into the earth's atmosphere where it will burn. The paper describes the concept, estimates the amount of velocity change (Delta V) that can be imparted to an object through laser-induced ablation, and investigates the use of a neutral particle beam for the momentum transfer. The space sweeper concept could also be extended to provide a collision avoidance system for the space station and satellites, or could be used for collision protection during interplanetary travel.

Description: The design and use of the Space Station Freedom are discussed. The contributions to the Station from EAS, Japan, Canada, and the US are described. Consideration is given to the capability of the Station, the internal accommodations for crew and payloads, various applications for the modules, and the planning and operation of the payloads.

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: Procedures for attitude determination based on Wahba's loss function are generalized to include the estimation of parameters other than the attitude, such as sensor biases. Optimization with respect to the attitude is carried out using the q-method, which does not require an a priori estimate of the attitude. Optimization with respect to the other parameters employs an iterative approach, which does require an a priori estimate of these parameters. Conventional state estimation methods require a priori estimates of both the parameters and the attitude, while the algorithm presented in this paper always computes the exact optimal attitude for given values of the parameters. Expressions for the covariance of the attitude and parameter estimates are derived.

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: The status and results to date of the NASA-Lewis/USAF Astronautics study of technology for large spacecraft heat-dissipation by means of liquid-droplet radiation (LDR) are discussed. The LDR concept uses a droplet generator to create billions of 200-micron droplets of a heatsink fluid which will cool through radiation into deep space as they fly toward a dropet collector. This exposure to the space environment entails the maintenance of vapor pressure as low as 10 to the -7th torr; the fluid must also be very stable chemically. While certain oils are good fluids for LDR use at low temperatures, higher-temperature heatsink fluids include Li, Sn, and Ga liquid metals.

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 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 review of unfurlable satellite antennas is presented. Typical application requirements for future space missions are first outlined. Then, U.S. and European mesh and inflatable antenna concepts are described. Precision deployables using rigid panels or petals are not included in the survey. RF modeling and performance analysis of gored or faceted mesh reflector antennas are then reviewed. Finally, both on-ground and in-orbit RF test techniques for large unfurlable antennas are discussed.

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 development and use of microspacecraft are examined. It is proposed that 10-50 microspacecraft per year can be launched; up to 50 microspacecraft can be dispatched with traditional launchers; and 1-3 experiments can be performed on the spacecraft. Various applications for the microspacecraft are discussed and specific examples of proposed missions are presented. Some systems and instruments designed for the microspacecraft 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: 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: 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: For 22 years (from ATS-1 to GOES-H) a single technology has dominated imaging from geosynchronous altitudes. In 1990, with the scheduled launch of GOES-I, a major change will occur which will in turn open the way for the Geostationary Platform. The need for improved observations of severe storms has led NOAA to a decision to replace spinning geostationary spacecraft with a three-axis-stabilized type (non-spinning) vehicle already common among communications spacecraft and demonstrated by INSAT. Also, the current spin-scan imager with sounder channels will be replaced by separate instruments capable of independent aiming. The advantages and challenges of the changes are discussed.

Description: A new, near-optimal feedback control technique is introduced that is shown to provide excellent vibration attenuation for those distributed parameter systems that are often encountered in the areas of aeroservoelasticity and large space systems. The technique relies on a novel solution methodology for the classical optimal control problem. Specifically, the quadratic regulator control problem for a flexible vibrating structure is first cast in a weak functional form that admits an approximate solution. The necessary conditions (first-order) are then solved via a time finite-element method. The procedure produces a low dimensional, algebraic parameterization of the optimal control problem that provides a rigorous basis for a discrete controller with a first-order like hold output. Simulation has shown that the algorithm can successfully control a wide variety of plant forms including multi-input/multi-output systems and systems exhibiting significant nonlinearities. In order to firmly establish the efficacy of the algorithm, a laboratory control experiment was implemented to provide planar (bending) vibration attenuation of a highly flexible beam (with a first clamped-free mode of approximately 0.5 Hz).

Description: Flexibility is important for high speed, high precision operation of lightweight manipulators. Accurate dynamic modeling of flexible robot arms is needed. Previous work has mostly been based on linear elasticity with prescribed rigid body motions (i.e., no effect of flexible motion on rigid body motion). Little or no experimental validation of dynamic models for flexible arms is available. Experimental results are also limited for flexible arm control. Researchers include the effects of prismatic as well as revolute joints. They investigate the effect of full coupling between the rigid and flexible motions, and of axial shortening, and consider the control of flexible arms using only additional sensors.

Description: A flexible aircraft or space structure with active control is typically modeled by a large-order state space system of equations in order to accurately represent the rigid and flexible body modes, unsteady aerodynamic forces, actuator dynamics and gust spectra. The control law of this multi-input/multi-output (MIMO) system is expected to satisfy multiple design requirements on the dynamic loads, responses, actuator deflection and rate limitations, as well as maintain certain stability margins, yet should be simple enough to be implemented on an onboard digital microprocessor. A software package for performing an analog or digital control law synthesis for such a system, using optimal control theory and constrained optimization techniques is described.

Description: Simulation of structural response of multi-flexible-body systems by linearized flexible motion combined with nonlinear rigid motion is discussed. Advantages and applicability of such an approach for accurate simulation with greatly reduced computational costs and turnaround times are described, restricting attention to the control design environment. Requirements for updating the linearized flexibility model to track large angular motions are discussed. Validation of such an approach by comparison with other existing codes is included. Application to a flexible robot manipulator system is described.

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: A simulation of the system made up of the Orbiter, Remote Manipulator System (RMS), and payload grappled by the RMS was completed. The simulation was used to study the stability of this overall system when its attitude is under control of the Orbiter's on-orbit Flight Control System (FCS). The simulation was also used to study the dynamics of the system when the RMS and its associated command software are in active control of the relative Orbiter to payload position and orientation. The simulation models all of the following elements: RMS boom bending (represented by two cubic bending models); RMS boom torsion; RMS joint gearbox compliance (represented by a non-linear wind-up model); flexibility at the RMS to Orbiter interface; flexibility at the RMS to payload interface; joint motor dynamics; joint servo-loop dynamics; RMS on-board computer command logic; data transfer delays between the RMS sensor and the RMS on-board computer and between the RMS on-board computer and RMS joint servos; on-orbit flight control nonlinear control logic; and the Reaction Control System (both Primary and Vernier) jet forces and moments.

Description: A three-dimensional finite element formulation for modeling the transient dynamics of constrained multibody space sructures with truss-like configurations is presented. Convected coordinate systems are used to define rigid-body motion of individual elements in the system. These systems are located at one end of each element and are oriented such that one axis passes through the other end of the element. Deformation of each element, relative to its convected coordinate system, is defined by cubic flexural shape functions as used in finite element methods of structural analysis. The formulation is oriented toward joint dominated structures and places the generalized coordinates at the joint. A transformation matrix is derived to integrate joint degree-of-freedom into the equations of motion of the element. Based on the derivation, a general-purpose code LATDYN (Large Angle Transient DYNamics) was developed. Two examples are presented to illustrate the application of the code. For the spin-up of a flexible beam, results are compared with existing solutions available in the literature. For the deployment of one bay of a deployable space truss (the Minimast), results are verified by the geometric knowledge of the system and converged solution of a successively refined model.

Description: Described here are Boeing software tools used for the development of control laws of flexible structures. The Boeing Company has developed a software tool called Modern Control Software Package (MPAC). MPAC provides the environment necessary for linear model development, analysis, and controller design for large models of flexible structures. There are two features of MPAC which are particularly appropriate for use with large models: (1) numerical accuracy and (2) label-driven nature. With the first feature MPAC uses double precision arithmetic for all numerical operations and relies on EISPAC and LINPACK for the numerical foundation. With the second feature, all MPAC model inputs, outputs, and states are referenced by user-defined labels. This feature allows model modification while maintaining the same state, input, and output names. In addition, there is no need for the user to keep track of a model variable's matrix row and colunm locations. There is a wide range of model manipulation, analysis, and design features within the numerically robust and flexible environment provided by MPAC. Models can be built or modified using either state space or transfer function representations. Existing models can be combined via parallel, series, and feedback connections; and loops of a closed-loop model may be broken for analysis.

Description: The goal is to develop the next generation guidance and control analysis and design tools to enable future missions and to improve productivity and reliability.

Description: A NASA program is about to start which has the objective to advance Controls-Structures Interaction (CSI) technology to a point where it can be used in spacecraft design for future missions. Because of the close interrelationships between the structure, the control hardware, and the analysis/design, a highly interdisciplinary activity is defined in which structures, dynamics, controls, computer and electronics engineers work together on a daily basis and are co-located to a large extent. Methods will be developed which allow the controls and structures analysis and design functions to use the same mathematical models. Hardware tests and applications are emphasized and will require development of concepts and test methods to carry out. Because of a variety of mission application problem classes, several time-phased, focus ground test articles are planned. They will be located at the Langley Researdh Center (LaRC), the Marshall Space Flight Center (MSFC) and at the Jet Propulsion Laboratory (JPL). It is anticipated that the ground tests will be subject to gravity and other environmental effects to the extent that orbital flights tests will be needed for verification of some technology items. The need for orbital flight experiments will be quantified based on ground test results and mission needs. Candidate on-orbit experiments will be defined and preliminary design/definition and cost studies will be carried out for one or more high-priority experiments.

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: Actual design turn-around time has become shorter due to the use of optimization techniques which have been introduced into the design process. It seems that what, how and when to use these optimization techniques may be the key factor for future aircraft engineering operations. Another important aspect of this technique is that complex physical phenomena can be modeled by a simple mathematical equation. The new powerful multilevel methodology reduces time-consuming analysis significantly while maintaining the coupling effects. This simultaneous analysis method stems from the implicit function theorem and system sensitivity derivatives of input variables. Use of the Taylor's series expansion and finite differencing technique for sensitivity derivatives in each discipline makes this approach unique for screening dominant variables from nondominant variables. In this study, the current Computational Fluid Dynamics (CFD) aerodynamic and sensitivity derivative/optimization techniques are applied for a simple cone-type forebody of a high-speed vehicle configuration to understand basic aerodynamic/structure interaction in a hypersonic flight condition.

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: The problem is to develop analysis methods, modeling stategies, and simulation tools to predict with assurance the on-orbit performance and integrity of large complex space structures that cannot be verified on the ground. The problem must incorporate large reliable structural models, multi-body flexible dynamics, multi-tier controller interaction, environmental models including 1g and atmosphere, various on-board disturbances, and linkage to mission-level performance codes. All areas are in serious need of work, but the weakest link is multi-body flexible dynamics.

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: The double lead spiral platen parallel jaw end effector is an extremely powerful, compact, and highly controllable end effector that represents a significant improvement in gripping force and efficiency over the LaRC Puma (LP) end effector. The spiral end effector is very simple in its design and has relatively few parts. The jaw openings are highly predictable and linear, making it an ideal candidate for remote control. The finger speed is within acceptable working limits and can be modified to meet the user needs; for instance, greater finger speed could be obtained by increasing the pitch of the spiral. The force relaxation is comparable to the other tested units. Optimization of the end effector design would involve a compromise of force and speed for a given application.

Description: The Deployable Retrievable Boom which was developed as a part of the joint U.S.-Italian Tethered Satellite System (TSS) is described. The design mission of the boom is to support, deploy, and retrieve an experiment package for the study of the electromagnetic field surrounding the satellite. The mechanism includes a jettisoning provision and deployable harness for the supported payloads connection. The boom is based on a tubular telescopic concept. Particular emphasis is placed on the payload harness connection capability and safety provisions.

Description: A brief listing of the concerns of the working group on spacecraft charging is presented. Brief conclusions for each concern is also given.

Description: The effects of spacecraft charging on spacecraft materials are studied. Spacecraft charging interactions seem to couple environment to system performance through materials. Technology is still developing concerning both environment-driven and operating system-driven interactions. The meeting addressed environment but lacked specific mission requirements, as a result system definition are needed to prioritize interactions.

Description: An overview is presented of potential interactions that can occur on spacecraft operating in space environments. These interactions are discussed in detail. The environment acts on spacecraft in such a way that charging of exterior surfaces occurs. The consequences from this charging then affect system operational performance. Hence, it is the coupling of this exterior charging to system performance that is of concern here. These interactions were first discovered in the spacecraft charging phenomena in which the geomagnetic substorms charged external surfaces to a level that discharges occurred. As a result of the discharge, electronic systems either changed logic state (anomalous switching) or failed. These interactions can occur in all orbits. The type associated with geosynchronous orbits is called passive since the environment provides the charging mechanism. This type can also occur in polar orbits due to auroral charging environments. In low Earth orbits, the thermal plasma alleviates charging environment concerns, but system operations can induce similar effects.

Description: Charging of dielectrics is a bulk, not a surface property. Radiation driven charge stops within the bulk and is not quickly conducted to the surface. Very large electric fields develop in the bulk due to this stopped charge. At space radiation levels, it typically requires hours or days for the internal electric fields to reach steady state. The resulting electric fields are large enough to produce electrical failure within the insulator. This type failure is thought to produce nearly all electric discharge anomalies. Radiation also induces bond breakage, creates reactive radicals, displaces atoms and, in general, severely changes the chemistry of the solid state material. Electric fields can alter this process by reacting with charged species, driving them through the solid. Irradiated polymers often lose as much as a percent of their mass, or more, at exposures typical in space. Very different aging or contaminant emission can be induced by the stopped charge electric fields. These radiation effects are detailed.

Description: The effects of thermal cycling on the thermal and mechanical properties of composite materials that are candidates for space structures are briefly described. The results from a thermal analysis of the orbiting Space Station Freedom is used to define a typical thermal environment and the parameters that cause changes in the thermal history. The interactions of this environment with composite materials are shown and described. The effects of this interaction on the integrity as well as the properties of GR/thermoset, Gr/thermoplastic, Gr/metal and Gr/glass composite materials are discussed. Emphasis is placed on the effects of the interaction that are critical to precision spacecraft. Finally, ground test methodology are briefly discussed.

Description: Component qualification and acceptance temperatures are derived from worst case thermal analyses and analytic uncertainty margin subject to certain specified temperature extremes. Temperature requirements are shown for equipment operation within specification and for survival and turn-on (need not operate within specification, but must not experience any degradation when returned to operational range). Temperature excursions for most equipment are seen to be 20 to 50 C above and below room temperature. Components without active electronics which are mounted outboard, such as solar arrays and antennas, are usually designed to withstand wider temperature excursions, particularly at the cold end. Batteries are tightly controlled at cold temperatures to increase life. Payload components such as extremely accurate clocks for precise navigation are controlled over a relatively narrow temperature range.

Description: Effective contamination control must encompass all aspects of ground and flight from design of the system through the end of mission life. Design systems are needed to minimize sensitivity to contamination, ease of cleaning, and contaminant production. Facilities and procedures are critical to maintaining cleanliness during ground operations. Flight operations should be planned so as to minimize contamination. More data from flights are required to assess the adequacy of designs and operations. Standards and specifications should include contamination control requirements.

Description: In the flight materials exposure data base extensive quantitative data is available from limited exposures in a narrow range of orbital environments. More data is needed in a wider range of environments as well as longer exposure times. Synergistic effects with other environmental factors; polar orbit and higher altitude environments; and real time materials degradation data is needed to understand degradation kinetics and mechanism. Almost no laboratory data exists from high fidelity simulations of the LEO environment. Simulation and test system are under development, and the data base is scanty. Theoretical understanding of hyperthermal atom surface reactions in the LEO environment is not good enough to support development of reliable accelerated test methods. The laser sustained discharge, atom beam sources are the most promising high fidelity simulation-test systems at this time.

Description: The Spacecraft Anomaly Data Base was useful in identifying trends in anomaly occurrence. Trends alone do not provide quantitative testimony to a spacecraft's reliability, but they do indicate areas that command closer study. An in-depth analysis of a specific anomaly can be expensive and difficult without access to the spacecraft. Statistically verified anomaly trends can provide a good reference point to begin anomaly analysis. Many spacecraft experience an increase in anomalies during the period of several days centered on the solar equinox, a period that is also correlated with sun eclipse at geostationary altitude and an increase in major geomagnetic storms. Increase anomaly occurrence can also be seen during the local time interval between midnight and dawn. This local time interval represents a region in Earth's near space that experiences an enhancement in electron plasma density due to a migration from the magnetotail during or following a geomagnetic substorm.

Description: The long-term performance of structural materials in the space environment is a key research activity within NASA. The primary concerns for materials in low Earth orbit (LEO) are atomic oxygen erosion and space debris impact. Atomic oxygen studies have included both laboratory exposures in atomic oxygen facilities and flight exposures using the Shuttle. Characterization of atomic oxygen interaction with materials has included surface recession rates, residual mechanical properties, optical property measurements, and surface analyses to establish chemical changes. The Long Duration Exposure Facility (LDEF) is scheduled to be retrieved in 1989 and is expected to provide a wealth of data on atomic oxygen erosion in space. Hypervelocity impact studies have been conducted to establish damage mechanisms and changes in mechanical properties. Samples from LDEF will be analyzed to determine the severity of space debris impact on coatings, films, and composites. Spacecraft placed in geosynchronous Earth orbit (GEO) will be subjected to high doses of ionizing radiation which for long term exposures will exceed the damage threshold of many polymeric materials. Radiation interaction with polymers can result in chain scission and/or cross-linking. The formation of low molecular weight products in the epoxy plasticize the matrix at elevated temperatures and embrittle the matrix at low temperatures. This affects both the matrix-dominated mechanical properties and the dimensional stability of the composite. Embrittlement of the matrix at low temperatures results in enhanced matrix microcracking during thermal cycling. Matrix microcracking changes the coefficient of thermal expansion (CTE) of composite laminates and produces permanent length changes. Residual stress calculations were performed to estimate the conditions necessary for microcrack development in unirradiated and irradiated composites. The effects of UV and electron exposure on the optical properties of transparent polymer films were also examined to establish the optimum chemical structure for good radiation resistance. Thoughts on approaches to establishing accelerated testing procedures are discussed.

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.