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: The uncertainties inherent in the dynamics of flexible spacecraft make robustness questions very important when designing vibration suppression systems for these vehicles. This often leads to the use of sensors and actuators which are collocated on the structure, so as to avoid the potentially destabilizing problem of unknown modal phase differences between noncollocated sensors and actuators. The closed-loop performance attainable is, of course, reduced if only collocated pairs are considered. The object of the present study is to investigate whether the special properties of the transmission zeros of such structures can be used to quantify these performance limitations. The approach is based on recent results concerning the transmission zeros of flexible structures.

Description: A generalized master/slave technique and experimental results for coordinated control of two arms rigidly grasping an object is described. An interactive program has been developed to allow a user the flexibility to select appropriate control modes for a given experiment. This interface allows for control gain adjustments. The results of several experiments performed on this system to demonstrate its capabilities such as transporting an object with or without induced internal forces and movement of a constrained object are offered. The system is further developed to achieve a so-called shared control mode in which an operator specifies the free motion trajectory for a point on the object of manipulation via a joystick while the autonomous control system is used for coordination and control of the arms.

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: Algorithmic complexity is discussed as a computational counterpart to the second law of thermodynamics. It is shown that algorithmic complexity, which is a measure of randomness, sets limits on the thermodynamic cost of computations and casts a new light on the limitations of Maxwell's demon. Algorithmic complexity can also be used to define distance between binary strings.

Description: Two procedures for the evaluation of the performance of expert systems are illustrated: one procedure evaluates predictive accuracy; the other procedure is complementary in that it uncovers the factors that contribute to predictive accuracy. Using these procedures, it is argued that expert systems should be more accurate than human experts in two senses. One sense is that expert systems must be more accurate to be cost-effective. Previous research is reviewed and original results are presented which show that simple statistical models typically perform better than human experts for the task of combining evidence from a given set of information sources. The results also suggest the second sense in which expert systems should be more accurate than human experts. They reveal that expert systems should share factors that contribute to human accuracy, but not factors that detract from human accuracy. Thus the thesis is that one should both require and expect systems to be more accurate than humans.

Description: The decoder error probability for Reed-Solomon codes (more generally, linear maximum distance separable codes) is examined. McEliece and Swanson offered an upper bound on P sub E (u), the decoder error probability given that u symbol errors occur. This upper bound is slightly greater than Q, the probability that a completely random error pattern will cause decoder error. By using a combinatoric technique, the principle of inclusion and exclusion, an exact formula for P sub E (u) is derived. The P sub E (u) for the (255,223) Reed-Solomon Code used by NASA, and for the (31,15) Reed-Solomon code (JTIDS code), are calculated using the exact formula, and the P sub E (u)'s are observed to approach the Q's of the codes rapidly as u gets larger. An upper bound for the expression is derived, and is shown to decrease nearly exponentially as u increases. This proves analytically that P sub E (u) indeed approaches Q as u becomes large, and some laws of large numbers come into play.

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: Simple state-space formulas are derived for all controllers solving the standard H(infinity) problem of finding, for a given number gamma greater than 0, all controllers such that the H(infinity) norm of the closed-loop transfer function is (strictly) less than gamma. It is known that a controller exists if and only if the unique stabilizing solutions to two algebraic Riccati equations are positive definite and the spectral radius of their product is less than gamma squared. Under these conditions, a parameterization of all controllers solving the problem is given as a linear fractional transformation (LFT) on a contractive stable free parameter. The state dimension of the coefficient matrix for the LFT, constructed using the two Riccati solutions, equals that of the plant and has a separation structure reminiscent of classical LQG (i.e., H2) theory. A standard H2 solution is developed in parallel.

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: Recursive forward dynamics algorithms are developed and presented for an arbitrary number of robot arms moving a commonly held object. The multiarm forward dynamics problem is to find the angular accelerations at the joints and the contact forces that the arms impart to the task object. The problem also involves finding the acceleration of this object. The multiarm forward dynamics solutions provide a thorough physical and mathematical understanding of the way several arms behave in response to a set of applied joint moments. Such an understanding simplifies and guides the subsequent control design and experimentation process. The forward dynamics algorithms also provide the necessary analytical foundation for conducting analysis and simulation studies. The multiarm algorithms are based on the filtering and smoothing approach recently advanced for single-arm dynamics, and they can be built up modularly from the single-arm algorithms. The algorithms compute recursively the joint angle accelerations, the contact forces, and the task-object accelerations. Algorithms are also developed to evaluate in closed form the linear transformations from the active joint moments to the joint angle accelerations, to the task object accelerations, and to the task-object contact forces. A possible computing architecture is presented as a precursor to a more complete investigation of the computational performance of the dynamics algorithms.

Description: A status review is presented for a generic blackboard-based distributed problem-solving environment in which multiple-agent cooperation can be effected. This environment is organized into a shared information panel, a chairman control panel, and a metaplanning panel. Each panel contains a number of embedded AI techniques that facilitate its operation and that provide heuristics for solving the underlying team-agent decision problem. The status of these panels and heuristics is described along with a number of robustness considerations. The techniques for each of the three panels and for four sets of paradigm-related advances are described, along with selected results from classroom teaching experiments and from three applications.

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: The linear optimal control problem of an nth-order time-invariant dynamic system with a quadratic performance functional is usually solved by the Hamilton-Jacobi approach. This leads to the solution of the differential matrix Riccati equation with a terminal condition. The bulk of the computation for the optimal control problem is related to the solution of this equation. There are various algorithms in the literature for solving the matrix Riccati equation. However, computational complexities and storage requirements as a function of numbers of state variables, control variables, and sensors are not available for all these algorithms. In this work, the computational complexities and storage requirements for some of these algorithms are given. These expressions show the immensity of the computational requirements of the algorithms in solving the Riccati equation for large-order systems such as the control of highly flexible space structures. The expressions are also needed to compute the speedup and efficiency of any implementation of these algorithms on concurrent machines.

Description: The authors propose a model, or filter, for debiasing opinions from multiple experts and combining them into a single consistent estimate of some variable of interest. A distinguishing feature of the approach consists of making the calibration of experts an integral part of filtering. This enables the filter to learn from previous experience with the experts. The theoretical development takes a Bayesian perspective, using de Finetti's (1964) notion of exchangeability. Experimental results with a preliminary computer implementation of the filter show that its estimates are better than those from comparable filters that do not involve calibration.

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 approximate inversion of an internally unknown linear system, given by its impulse response sequence, by an inverse system having a finite impulse response, is considered. The recursive least-squares procedure is shown to have an exact initialization, based on the triangular Toeplitz structure of the matrix involved. The proposed approach also suggests solutions to the problem of system identification and compensation.

Description: A methodology for robust eigenstructure assignment for multivariable feedback systems is presented. The algorithm is based upon a pole placement technique using projections onto subspaces of admissible eigenvectors. New ideas are introduced to generate target (desired) sets of unitary eigenvectors and determine optimal feasible eigenvectors in a least-square sense. Useful connections are established between the pole-placement by independent modal space control and the method introduced in this paper. A multicriterion optimization algorithm is also presented, which takes efficient advantage of the present eigenstructure assignment method. These developments show significant improvement over an earlier version of this algorithm in both computational cost and accuracy. This optimization process appears to be numerically robust and suitable for high-dimensional multicriterion optimizations; it is especially attractive for computer-aided design of control systems.

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: This paper corrects an unclear statement of the conventional Kalman filter implementation as presented by M. H. Verhaegen and P. van Dooren in Numerical aspects of different Kalman filter implementations, IEEE Trans. Automat. Contr., v. AC-31, no. 10, pp. 907-917, 1986. It is shown that the habitual, incorrect implementation of the Kalman filter has been the major cause of its insensitivity to the so-called loss-of-symmetry phenomenon.

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: Problems of implementing rule-based expert systems using fuzzy sets are considered. A fuzzy logic software development shell is used that allows inclusion of both crisp and fuzzy rules indecision making and process control problems. Results are given that compare this type of expert system to a human expert in some specific applications. Advantages and disadvantages of such systems are discussed.

Description: A large scale domain-independent spatial data management expert system that serves as a front-end to databases containing spatial data is described. This system is unique for two reasons. First, it uses spatial search techniques to generate a list of all the primary keys that fall within a user's spatial constraints prior to invoking the database management system, thus substantially decreasing the amount of time required to answer a user's query. Second, a domain-independent query expert system uses a domain-specific rule base to preprocess the user's English query, effectively mapping a broad class of queries into a smaller subset that can be handled by a commercial natural language processing system. The methods used by the spatial search module and the query expert system are explained, and the system architecture for the spatial data management expert system is described. The system is applied to data from the International Ultraviolet Explorer (IUE) satellite, and results are given.

Description: Existing NASA supported scientific data bases are usually developed, managed and populated in a tedious, error prone and self-limiting way in terms of what can be described in a relational Data Base Management System (DBMS). The next generation Earth remote sensing platforms, i.e., Earth Observation System, (EOS), will be capable of generating data at a rate of over 300 Mbs per second from a suite of instruments designed for different applications. What is needed is an innovative approach that creates object-oriented databases that segment, characterize, catalog and are manageable in a domain-specific context and whose contents are available interactively and in near-real-time to the user community. Described here is work in progress that utilizes an artificial neural net approach to characterize satellite imagery of undefined objects into high-level data objects. The characterized data is then dynamically allocated to an object-oriented data base where it can be reviewed and assessed by a user. The definition, development, and evolution of the overall data system model are steps in the creation of an application-driven knowledge-based scientific information system.

Description: This paper describes a novel data acquisition system for use with wind-tunnel probe-based measurements, which incorporates a degree of specific fluid dynamics knowledge into a simple expert system-like control program. The concept was developed with a rudimentary expert system coupled to a probe positioning mechanism operating in a small-scale research wind tunnel. The software consisted of two basic elements, a general-purpose data acquisition system and the rulebased control element to take and analyze data and supplying decisions as to where to measure, how many data points to take, and when to stop. The system was validated in an experiment involving a vortical flow field, showing that it was possible to increase the resolution of the experiment or, alternatively, reduce the total number of data points required, to achieve parity with the results of most conventional data acquisition approaches.

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: The spatially constrained clustering (SCC) iterative parallel region-growing technique is applied to image analysis. The SCC algorithm is implemented on the massively parallel processor at NASA Goddard. Most previous region-growing approaches have the drawback that the segmentation produced depends on the order in which portions of the image are processed. The ideal solution to this problem (merging only the single most similar pair of spatially adjacent regions in the image in each iteration) becomes impractical except for very small images, even on a massively parallel computer. The SCC algorithm overcomes these problems by performing, in parallel, the best merge within each of a set of local, possibly overlapping, subimages. A region-splitting stage is also incorporated into the algorithm, but experiments show that region splitting generally does not improve segmentation results. The SCC algorithm has been tested on various imagery data, and test results for a Landsat TM image are summarized.

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

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

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

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

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

Description: 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: A suite of telerobotic tasks has been compiled and assessed for the purpose of selecting viable tasks for near and far term laboratory demonstrations. The primary intent of developing the task is to provide some technical guidelines, with supporting data, for focusing laboratory demonstrations toward application domains that address a wide array of potential telerobot tasks and required technologies. This wide application would then result in a rich technology development environment to meet the broad task requirements of a system such as the Flight Telerobot Servicer. The methodology and results of the telerobot task assessment are described, including a ranking of the final select suite of major tasks. The presented along with guidelines for both interpreting the task ranking results and setting programmatic objectives based on these results. Detailed data about the task candidates and their respective levels of complexity, task primitive actions, and the actual relative measures of task worth as associated with key tradeoff variables such as cost, available research resources, technology availability, and importance to the user community are also presented.

Description: The results of using a new technique for shaping inputs to a model of the space shuttle Remote Manipulator System (RMS) are presented. The shapes inputs move the system to the same location that was originally commanded, however, the oscillations of the machine are considerably reduced. An overview of the new shaping method is presented. A description of RMS model is provided. The problem of slow joint servo rates on the RMS is accommodated with an extension of the shaping method. The results and sample data are also presented for both joint and three-dimensional cartesian motions. The results demonstrate that the new shaping method performs well on large, telerobotic systems which exhibit significant structural vibration. The new method is shown to also result in considerable energy savings during operations of the RMS manipulator.

Description: Four separate projects recently completed or in progress at the MIT Man-Machine Systems Laboratory are summarized. They are: a decision aid for retrieving a tumbling satellite in space; kinematic control and graphic display of redundant teleoperators; real time terrain/object generation: a quad-tree approach; and two dimensional control for three dimensional obstacle avoidance.

Description: The implementation and experimental validation of a new direct adaptive control scheme on a PUMA 560 industrial robot is described. The testbed facility consists of a Unimation PUMA 560 six-jointed robot and controller, and a DEC MicroVAX II computer which hosts the Robot Control C Library software. The control algorithm is implemented on the MicroVAX which acts as a digital controller for the PUMA robot, and the Unimation controller is effectively bypassed and used merely as an I/O device to interface the MicroVAX to the joint motors. The control algorithm for each robot joint consists of an auxiliary signal generated by a constant-gain Proportional plus Integral plus Derivative (PID) controller, and an adaptive position-velocity (PD) feedback controller with adjustable gains. The adaptive independent joint controllers compensate for the inter-joint couplings and achieve accurate trajectory tracking without the need for the complex dynamic model and parameter values of the robot. Extensive experimental results on PUMA joint control are presented to confirm the feasibility of the proposed scheme, in spite of strong interactions between joint motions. Experimental results validate the capabilities of the proposed control scheme. The control scheme is extremely simple and computationally very fast for concurrent processing with high sampling rates.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: A recursive algorithm is developed for the solution of the simulation dynamics problem for a chain of rigid bodies. Arbitrary joint constraints are permitted, that is, joints may allow translational and/or rotational degrees of freedom. The recursive procedure is shown to be identical to that encountered in a discrete-time optimal control problem. For each relevant quantity in the multibody dynamics problem, there exists an analog in the context of optimal control. The performance index that is minimized in the control problem is identified as Gibbs' function for the chain of bodies.

Description: Space telerobots are recognized to require cooperation with human operators in various ways. Multi-level manual and autonomous control superposition in telerobot task execution is described. The object model, the structured master-slave manipulation system, and the motion understanding system are proposed to realize the concept. The object model offers interfaces for task level and object level human intervention. The structured master-slave manipulation system offers interfaces for motion level human intervention. The motion understanding system maintains the consistency of the knowledge through all the levels which supports the robot autonomy while accepting the human intervention. The superposing execution of the teleoperational task at multi-levels realizes intuitive and robust task execution for wide variety of objects and in changeful environment. The performance of several examples of operating chemical apparatuses is shown.

Description: A classical control perspective is used to characterize performance constraints and evaluate compensation techniques for teleoperation with delay. Use of control concepts such as open and closed loop performance, stability, and bandwidth yield insight to the delay problem. Teleoperator performance constraints are viewed as an open loop time delay lag and as a delay-induced closed loop bandwidth constraint. These constraints are illustrated with a simple analytical tracking example which is corroborated by a real time, 'man-in-the-loop' tracking experiment. The experiment also provides insight to those controller characteristics which are unique to a human operator. Predictive displays and feedforward commands are shown to provide open loop compensation for delay lag. Low pass filtering of telemetry or feedback signals is interpreted as closed loop compensation used to maintain a sufficiently low bandwidth for stability. A new closed loop compensation approach is proposed that uses a reactive (or force feedback) hand controller to restrict system bandwidth by impeding operator inputs.

Description: Hierarchical structures offer numerous advantages over conventional structures for the control of telerobotic systems. A hierarchically organized system can be controlled via undetailed task assignments and can easily adapt to changing circumstances. The distributed and modular structure of these systems also enables fast response needed in most telerobotic applications. On the other hand, most of the hierarchical structures proposed in the literature are based on functional properties of a system. These structures work best for a few given functions of a large class of systems. In telerobotic applications, all functions of a single system needed to be explored. This approach requires a hierarchical organization based on physical properties of a system and such a hierarchical organization is introduced. The decomposition, organization, and control of the hierarchical structure are considered, and a system with two robot arms and a camera is presented.

Description: Kolmogorov's axiomatic principles of the probability theory, are reconsidered in the scope of their applicability to the processes of knowledge acquisition and interpretation. The model of uncertainty generation is modified in order to reflect the reality of engineering problems, particularly in the area of intelligent control. This model implies algorithms of learning which are organized in three groups which reflect the degree of conceptualization of the knowledge the system is dealing with. It is essential that these algorithms are motivated by and consistent with the multiresolutional model of knowledge representation which is reflected in the structure of models and the algorithms of learning.

Description: A neural network for controlling the configuration of frame structure with elastic members is proposed. In the present network, the structure is modeled not by using the relative angles of the members but by using the distances between the joint locations alone. The relationship between the environment and the joints is also defined by their mutual distances. The analog neural network attains the reaching motion of the manipulator as a minimization problem of the energy constructed by the distances between the joints, the target, and the obstacles. The network can generate not only the final but also the transient configurations and the trajectory. This framework with flexibility and parallelism is very suitable for controlling the Space Telerobotic systems with many degrees of freedom.

Description: Two neural learning controller designs for manipulators are considered. The first design is based on a neural inverse-dynamics system. The second is the combination of the first one with a neural adaptive state feedback system. Both types of controllers enable the manipulator to perform any given task very well after a period of training and to do other untrained tasks satisfactorily. The second design also enables the manipulator to compensate for unpredictable perturbations.