ALBERT

Description: Two hybrid methods of deformation analysis combining moire interferometry and finite element techniques are introduced. The hybrid methods are applied to determine the strain fields in the test section of a composite cruciform specimen in order to evaluate the uniformity of the shear strain field. Both methods show that there are significant shear strain gradients in the specimen test section. One of the hybrid methods - the full-field analysis using the random-node mesh generator - detects the presence of material property nonuniformity which is assumed uniform in the boundary-node method.

Description: Results are presented of a preliminary investigation of instrumented fasteners for use as sensors to measure the shear loads transmitted by individual fasteners installed in double splice joints. Calibration and load verification tests were conducted for instrumented fasteners installed at three fastener torque levels. Calibration test results show that the shear strains obtained from the instrumented fasteners vary linearly with the applied load and that the instrumented fasteners can be effectively used to measure shear loads transmitted by individual fasteners installed in double splice joints. The load distribution between individual fasteners is found to be dependent on the location of the fastener in the joint and the fastener torque level.

Description: Recent advances in structural design for control are reviewed. Attention is given to adaptive structures, passive damping, ground testing, and system identification. Directions for future research are suggested.

Description: The benefits of structural dynamics modeling methods in aerospace structures are reviewed. Four major issues in structural dynamics modeling are discussed which encompass most of its subdisciplines: reduced order modeling, constraints in problems with large motions, computational strategies, and fundamental methods. Directions for future research in these areas are addressed.

Description: The results from a study aimed at improving the dynamic and aerodynamic characteristics of composite rotor blades through the use of extension-twist elastic coupling are presented. A set of extension-twist-coupled composite tubular spars, representative of the primary load carrying structure within a helicopter rotor blade, was manufactured using four plies of woven graphite/epoxy cloth 'prepreg.' These spars were non-circular in cross section design and were therefore subject to warping deformations. Three cross-sectional geometries were developed: square, D-shape, and flattened ellipse. Results from free-free vibration tests of the spars were compared with results from normal modes and frequency analyses of companion shell-finite-element models developed in MSC/NASTRAN. Five global or 'non-shell' modes were identified within the 0-2000 Hz range for each spar. The frequencies and associated mode shapes for the D-shape spar were correlated with analytical results, showing agreement within 13.8 percent. Frequencies corresponding to the five global mode shapes for the square spar agreed within 9.5 percent of the analytical results. Five global modes were similarly identified for the elliptical spar and agreed within 4.9 percent of the respective analytical results.

Description: This paper presents an extension of the Coordinate Modal Assurance Criterion (COMAC) for spatial comparison of mode shapes. As is the case for the original COMAC, this enhanced COMAC helps identify degrees of freedom differences between test and analysis modes, but also overcomes some of the limitations of the original COMAC. The original COMAC cannot identify errors caused by scaling or polarity mistakes in the test data. Such mistakes are easily made during modal tests and not always easily detected since they do not affect the subjective quality of a frequency response function.

Description: To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Description: This paper presents an alternative to the correlation of individual components of a mode shape vectors by directly examining the sensitivity of the cross-orthogonality between test and analytical mode shapes. If the test and analysis mode shapes are identical, the diagonal elements of the cross-orthogonality will be identical to the test orthogonality matrix, so the cross-orthogonality matrix provides a concise measure of the 'closeness' between test and analysis mode shapes. There are two major advantages to the cross-orthogonality correlation approach. The first is that a direct correlation of this matrix will more directly meet the goal of the correlation effort (measured by cross-orthogonality). Secondly, and more importantly, the correlation of cross-orthogonality greatly reduces the amount of data that needs to be handled when compared to the correlation of mode shapes.

Description: NASA-Lewis is actively involved in the general effort to research, develop, test, and evaluate advanced theoretical, analytical, experimental, and probabilistic analysis concepts required for life prediction of liquid rocket engines at the subcomponent, component, and engine system levels. The models developed are oriented toward use in advanced health monitoring systems of space propulsion systems. It is planned to demonstrate the methodology considering a representative set of three components such as a main injector element, a combustion chamber liner, and a turbopump blade. This paper describes the initial development and application of this method to a specific location in the main injector element of the SSME. Further enhancements and various elements of the framework will be completed as the work proceeds in subsequent years.

Description: A numerical model was developed to calculate the interference pattern at the end of a multimode weakly guiding optical fiber under stress. Whenever an optical fiber is under stress, the modal phase in the interference term of the intensity formula changes. Plots of the simulated output of a stressed fiber are presented. For multimode fibers, very complicated patterns result. Under stress, lobes in the pattern are generated, displaced and power is exchanged among them.

Description: Tests were performed measuring the damage initiation loads and the locations, shapes, and sizes of delaminations in Fiberite T300/976 graphite/epoxy, Fiberite IM7/977-2 graphite-toughened epoxy, and ICI APC-2 graphite-PEEK plates subjected to transverse static loads. The data were compared to the results of the Finn-Springer model, and good agreements were found between the measured and calculated delamination lengths and widths.

Description: Present interest in hypersonic vehicles has resulted in a renewed interest in thermal stress analysis of airframe structures. While there are numerous texts and papers on thermal stress analysis, practical examples and experience on light gage aircraft structures are fairly limited. A research program has been undertaken at General Dynamics to demonstrate the present state of the art, verify methods of analysis, gain experience in their use, and develop engineering judgement in thermal stress analysis. The approach for this project has been to conduct a series of analyses of this sample problem and compare analysis results with test data. This comparison will give an idea of how to use our present methods of thermal stress analysis, and how accurate we can expect them to be.

Description: This paper presents the fuzzy dynamical reliability and failure probability as well as the basic principles and the analytical method of loss assessment for nonlinear seismic steel structures. Also presented is the optimization formulation and a numerical example for double objectives, initial construction cost and expected failure loss, and dynamical reliability constraints. The earthquake ground motion is based on a stationary filtered non-white noise and the fuzzy damage grade is described by damage index.

Description: Effective methods of approximate eigensolution reanalysis of modified nonclassically damped structures are developed. For structures with passive or active discrete damping devices or with damping treatment, the system becomes nonproportionally damped and the computation of its dynamic responses may require the use of complex modes. For larger systems, the computation of complex modes is very expensive. Thus it is desirable to have approximate reanalysis techniques for the efficient evaluation of the effect of design changes. In recent years, the assumed mode reanalysis method was successfully applied to minimum weight design of undamped structures with natural frequency constraints. The accuracy of the assumed mode reanalysis method can be improved dramatically if the global approximation function includes the normal modes of the original system and their derivatives. This approach was demonstrated to be effective even for a system with shape changes. The approach used by Noor et al. for eigensolution reanalysis of undamped structures is extended to treat a nonclassically damped system.

Description: With structural design in mind, a new unified variational model has been developed which represents the mechanics of deformation elasto-plasticity with unilateral contact conditions. For a design problem formulated as maximization of the load carrying capacity of a structure under certain constraints, the unified model allows for a simultaneous analysis and design synthesis for a whole range of mechanical behavior.

Description: The real benefit of structural optimization techniques is in the application of these techniques to large structures such as full vehicles or full aircraft. For these structures, however, the sequential computer's time and memory requirements prohibit the solutions. With the rapid development of parallel computers, parallel processing of large scale structural optimization problems is achievable. In this paper we discuss the parallel processing of structural optimization problems with parallel structural analysis. Two different types of interface between the optimization and analysis routines are developed and tested.

Description: In most structural optimization problems the implicit behavior constraints are evaluated for successive modifications in the design. For each trial design, the analysis equations must be solved and the multiple repeated analyses usually involve extensive computational effort. This difficulty motivated several studies on explicit approximations of the structural behavior in terms of the design variables. The latter approach can considerably reduce the amount of computations, but the quality of the approximations might not be sufficient. Many of the approximate behavior models proposed in the past are valid only for relatively small changes in the design variables. The accuracy of the results is often insufficient for large changes in the design. The object of this study is to present efficient and high quality approximations of the structural behavior. It will be shown that the quality of the approximations can greatly be improved by combining scaling of the initial design, using intervening variables, and scaling a set of fictitious loads. Integrating these means, a powerful solution procedure can be introduced. In addition, the errors in satisfying the analysis equations can readily be evaluated. A numerical example illustrates the solution methodology and the effectiveness of the proposed approach.

Description: The weight savings due to usage of composite materials in aircraft structural applications is well known. Significant weight and cost benefits are achievable by developing structurally tailored concepts and efficiently integrating them with suitable material and manufacturing technologies. The proposed paper will describe such an efficient concept for application to primary aircraft structures.

Description: This paper first presents a brief review on the application of optimization and active control of seismic structures along with some of the author's recent work. It then assesses the practicality and future development of seismic structural optimization, and some practical problems associated with active control.

Description: The paper presents the theoretical bases and implementation techniques of sensitivity analyses for efficient structural optimization of large structures, based on finite element static and dynamic analysis methods. The sensitivity analyses have been implemented in conjunction with two methods for optimization, namely, the Mathematical Programming and Optimality Criteria methods. The paper discusses the implementation of the sensitivity analysis method into our in-house software package, AutoDesign.

Description: Ill-conditioned systems arising in analysis and optimization can display a high sensitivity to numerical precision for changes and errors in data input. Such data may be in the form of system parameter input or desired system response. The ill-conditioning referred to generally arises from the lack of sufficient independent data to define a complex system or the weak sensitivity of response to source input parameters. It is shown how small errors in data and assumed fixed and known parameters can lead to highly erroneous results in ill-conditioned linear algebraic equations. A simplified detection and correction of critical input data arising in the coefficient matrix and desired response (i.e., right hand side) is proposed.

Description: An identification procedure proposed by Shen and Taylor to determine the crack characteristics (location and size of the crack) from dynamic measurements is tested. This procedure was based on minimization of either the 'mean-square' measure of difference between measurement data (natural frequencies and mode shapes) and the corresponding predictions obtained from the computational model. The procedure is tested for simulated damage in the form of symmetric cracks in a simply-supported Bernoulli-Euler beam. The sensitivity of the solution of damage identification problems to the values of parameters that characterize damage is discussed. A sensitivity formula is derived.

Description: The issue of the utility of multilevel decomposition and optimization remains controversial. To date, only the structural optimization community has actively developed and promoted multilevel optimization techniques. However, even this community acknowledges that multilevel optimization is ideally suited for a rather limited set of problems. It is warned that decomposition typically requires eliminating local variables by using global variables and that this in turn causes ill-conditioning of the multilevel optimization by adding equality constraints. The purpose is to suggest a new multilevel optimization technique. This technique uses behavior variables, in addition to design variables and constraints, to decompose the problem. The new technique removes the need for equality constraints, simplifies the decomposition of the design problem, simplifies the programming task, and improves the convergence speed of multilevel optimization compared to conventional optimization.

Description: In recent years there have been several hierarchic multilevel optimization algorithms proposed and implemented in design studies. Equality constraints are often imposed between levels in these multilevel optimizations to maintain system and subsystem variable continuity. Equality constraints of this nature will be referred to as coupling equality constraints. In many implementation studies these coupling equality constraints have been handled indirectly. This indirect handling has been accomplished using the coupling equality constraints' explicit functional relations to eliminate design variables (generally at the subsystem level), with the resulting optimization taking place in a reduced design space. In one multilevel optimization study where the coupling equality constraints were handled directly, the researchers encountered numerical difficulties which prevented their multilevel optimization from reaching the same minimum found in conventional single level solutions. The researchers did not explain the exact nature of the numerical difficulties other than to associate them with the direct handling of the coupling equality constraints. The coupling equality constraints are handled directly, by employing the Generalized Reduced Gradient (GRG) method as the optimizer within a multilevel linear decomposition scheme based on the Sobieski hierarchic algorithm. Two engineering design examples are solved using this approach. The results show that the direct handling of coupling equality constraints in a multilevel optimization does not introduce any problems when the GRG method is employed as the internal optimizer. The optimums achieved are comparable to those achieved in single level solutions and in multilevel studies where the equality constraints have been handled indirectly.

Description: Reliability was determined for two degrading dynamic systems subject to random load processes. Damage is caused by loss of components for Daniels systems and crack extension for plates with cracks. The analysis accounted for the coupling between response and current damage state of the system. It is based on mean crossing rates of conditional processes and properties of diffusion models. Simple systems are used to illustrate proposed methods for estimating reliability.

Description: A passive vibration damping technique that is referred to as 'Non-Obstructive Particle Damping (NOPD)' is presented. The NOPD technique consists of making small diameter holes (or cavities) at appropriate locations inside vibrating structures and filling these holes to appropriate levels with particles which yield the maximum damping effectiveness for the desired mode (or modes). Powders, spherical shaped, metallic, non-metallic or liquid particles (or mixtures) with different densities, viscosities and adhesive or cohesive characteristics can be used.

Description: Recent results from sensitivity analysis for strain energy with anisotropic elasticity are applied to thickness and orientational design of laminated membranes. Primarily, the first order gradients of the total elastic energy are used in an optimality criteria based method. This traditional method is shown to give slow convergence with respect to design parameters, although the convergence of strain energy is very good. To get a deeper insight into this rather general characteristic, second order derivatives are included and it is shown how they can be obtained by first order sensitivity analysis. Examples of only thickness design, only orientational design, and combined thickness--orientational design are presented.

Description: The objective is to set up and analyze benchmark problems on multibody dynamics and to verify the predictions of two multibody computer simulation codes. TREETOPS and DISCOS have been used to run three example problems - one degree-of-freedom spring mass dashpot system, an inverted pendulum system, and a triple pendulum. To study the dynamics and control interaction, an inverted planar pendulum with an external body force and a torsional control spring was modeled as a hinge connected two-rigid body system. TREETOPS and DISCOS affected the time history simulation of this problem. System state space variables and their time derivatives from two simulation codes were compared.

Description: The non-linear dynamic analysis of large structures is always very time, effort and CPU consuming. Whenever possible the reduction of the size of the mathematical model involved is of main importance to speed up the computational procedures. Such reduction can be performed for the part of the structure which perform linearly. Most of the time, the classical Guyan reduction process is used. For non-linear dynamic process where the non-linearity is present at interfaces between different structures, Craig-Bampton models can provide a very rich information, and allow easy selection of the relevant modes with respect to the phenomenon driving the non-linearity. The paper presents the employment of Craig-Bampton models combined with Newmark direct integration for solving non-linear friction problems appearing at the interface between the Hubble Space Telescope and its solar arrays during in-orbit maneuvers. Theory, implementation in the FEM code ASKA, and practical results are shown.

Description: Structural optimization problems are mostly solved under constraints from statics, such as stresses, strains, or displacements under static loads. But in the design process, dynamic quantities like eigenfrequencies or accelerations under dynamic loads become more and more important. Therefore, it is obvious that constraints from dynamics must be considered in structural optimization packages. This paper addresses the dynamics branch in MBB-LAGRANGE. It will concentrate on two topics, namely on the different formulations for eigenfrequency constraints and on frequency response constraints. For the latter the necessity of a system reduction is emphasized. The methods implemented in LAGRANGE are presented and examples are given.

Description: This paper addresses the weight minimization of a circular plate-like structure which resulted in a 26 percent weight reduction. The optimization was performed numerically with the COPES/ADS program using the modified method of feasible directions. Design parameters were the inner thickness and outer thickness of the plate with constraints on maximum yield stress and maximum transverse displacement. Also, constraints were specified for the upper and lower bounds of the fundamental frequency and plate thicknesses. The MSC/NASTRAN finite element program was used for the evaluation of response variables. Original and final designs of the plate were tested using an Instron tension-compression machine to compare finite element results to measured strain data. The difference between finite element strain components and measured strain data was within engineering accuracy.

Description: Any optimization of structures for maximum stability or for maximum dynamic stiffness deals with an eigenvalue problem. The goal of this optimization is to raise the lowest eigenvalue (or eigenvalues) of the problem to its highest (optimal) level at a constant volume of the structure. Likely the lowest eigenvalue may be either inherently multi-modal or it can become multi-modal as a result of the optimization process. The multimodeness introduces some ambiguity to the eigenvalue problem and make the optimization difficult to handle. Thus far, only the simplest cases of multi-modal structures have been effectively optimized using rather elaborate analytical methods. Numerous publications report design of a minimum volume structure with different eigenvalues constraints, in which, however, the modality of the problem is assumed a priori. The method presented here utilizes a multi-modal optimality criteria and allows for inclusion of an arbitrary number of buckling or vibrations modes which might influence the optimization process. The real multi-modality of the problem, that is the number of modes participating in the final optimal design is determined iteratively. Because of a natural use of the FEM technique the method is easy to program and might be helpful in design of large flexible space structures.

Description: The structure-control system optimization problem is formulated with constraints on the closed-loop eigenvalues and the efficiency of the reduced order system. The feasibility of the approach is illustrated by designing the ACOSS-FOUR structure with a reduced order system and improving the efficiency characteristics of the structures-control system.

Description: In this paper an optimal design procedure is introduced to improve the overall performance of nonlinear framed structures. The design methodology presented here is a multiple-objective optimization procedure whose objective functions involve the buckling eigenvalues and eigenvectors of the structure. A constant volume with bounds on the design variables is used in conjunction with an optimality criterion approach. The method provides a general tool for solving complex design problems and generally leads to structures with better limit strength and stability. Many algorithms have been developed to improve the limit strength of structures. In most applications geometrically linear analysis is employed with the consequence that overall strength of the design is overestimated. Directly optimizing the limit load of the structure would require a full nonlinear analysis at each iteration which would be prohibitively expensive. The objective of this paper is to develop an algorithm that can improve the limit-load of geometrically nonlinear framed structures while avoiding the nonlinear analysis. One of the novelties of the new design methodology is its ability to efficiently model and design structures under multiple loading conditions. These loading conditions can be different factored loads or any kind of loads that can be applied to the structure simultaneously or independently. Attention is focused on optimal design of space framed structures. Three-dimensional design problems are more complicated to carry out, but they yield insight into real behavior of the structure and can help avoiding some of the problems that might appear in planar design procedure such as the need for out-of-plane buckling constraint. Although researchers in the field of structural engineering generally agree that optimum design of three-dimension building frames especially in the seismic regions would be beneficial, methods have been slow to emerge. Most of the research in this area has dealt with the optimization of truss and plane frame structures.

Description: In many structures the final stress states are dependent on the sequence of construction or the stress states at various stages of construction are of interest. Such problems can be analyzed using finite element programs that have the capability of adding (birthing) elements to simulate the progress of construction. However, the usual procedure of assembling elements may lead to numerical instabilities or stress states that are unrealistic. Both problems are demonstrated in the analysis of a structure using the program ADINA. A technique which combines application of a preload with element birthing to overcome these problems is described and illustrated.

Description: A technique utilizing finite element analysis, liquid impact kinematics, and momentum theory is described and compared to single-drop impact test data performed on various configurations of coated ceramic material. The method correlates well with test data and is useful in predicting the single-drop impact damage velocity threshold for low-density, coated ceramic materials.

Description: An algorithm is outlined for simulating the contact surface erosion for impact problems. The algorithm dynamically relocates the contact surface as projectile and target materials exceed their failure criterion. Example computations of axisymmetric and oblique impacts are compared with experimental data.

Description: A study of the dynamic characteristics of a coupled translational-rotational system is given. The formulation of the problem considers the soil-structure interaction effects by utilizing the impedance functions at the foundation of a structure. Due to the fact that the coefficient matrix in the characteristic equation is frequency dependent in nature, iterations have to be performed to find the nature frequencies of the system. Examples and discussions are presented. Comparisons of the analytical results from various approaches are also given.

Description: Two analytical techniques applicable to large deflection dynamic response calculations for pressure loaded composite sandwich panels are demonstrated. One technique utilizes finite element modeling with a single equivalent layer representing the face sheets and core. The other technique utilizes the modal analysis computer code DEPROP which was recently modified to include transverse shear deformation in a core layer. The example problem consists of a simply supported rectangular sandwich panel. Included are comparisons of linear and nonlinear static response calculations, in addition to dynamic response calculations.

Description: The design of thin shell structures with respect to elastoplastic buckling requires an extended analysis of the influence of initial imperfections. For conservative design, the most critical defect should be assumed with the maximum allowable magnitude. This defect is closely related to the initial postbuckling behavior. An algorithm is given for the quasi-static analysis of the postbuckling behavior of structures that exhibit multiple buckling points. the algorithm based upon an energy criterion allows the computation of the critical perturbation which will be employed for the definition of the critical defect. For computational efficiency, the algorithm uses the reduced basis technique with automatic update of the modal basis. The method is applied to the axisymmetric buckling of cylindrical shells under axial compression, and conclusions are given for future research.

Description: A general approach is required for describing matter of behavior when the failure is likely to involve growth and coalescence of a large number of fractures. Failures of this kind appear frequently in rapid dynamic processes, particularly in the formation of spall fragments. An approach to formulating constitutive relations that accounts for the opening, shear and growth of an ensemble of cracks is discussed. The approach accounts for plastic flow accompanying fragmentation. The resulting constitutive relations were incorporated into a Lagrangian computer program. A theoretical approach to coalescence is described. The simplest formulation uses a linear Liouville equation, with crack growth limited by the mean free path of cracks, assumed constant. This approach allows for an anisotropic distribution of cracks. An alternative approach in which the decrease of the mean free path with increasing crack size is accounted for, but the crack distribution is assumed isotropic is described. A reduction of the governing Liouville equation to an ordinary differential equation of third order is possible, and the result can be used to determine how mean free path decreases with increasing crack size.

Description: The design of a composite panel requires some way of finding the minimum thickness laminate which will withstand the load requirements without failure. The mathematical complexity of this problem dictates the use of nonlinear optimization techniques. Specialized laminate optimization programs were developed which are compact and efficient enough to run on microcomputers. Only stresses at a point and inplane loads and deflections are considered. The programs are simple to use and require no knowledge of optimization. Techniques are developed which find minimum thickness laminates with either ply ratios or ply angles as design variables. A method is presented for finding the optimum orientation for the axis of symmetry of an orthotropic laminate. The orthotropic laminate program uses an approximate failure theory, which speed up computations dramatically.

Description: The evaluation of the total probability of a plastic collapse failure P sub f for a highly redundant structure of random interdependent plastic moments acted on by random interdepedent loads is a difficult and computationally very costly process. The evaluation of reasonable bounds to this probability requires the use of second moment algebra which involves man statistical parameters. A computer program which selects the best strategy for minimizing the interval between upper and lower bounds of P sub f is now in its final stage of development. The relative importance of various uncertainties involved in the computational process on the resulting bounds of P sub f, sensitivity is analyzed. Response sensitivities for both mode and system reliability of an ideal plastic portal frame are shown.

Description: A finite element modeling technique which utilizes a triangular element called TSHEL with 45 degrees of freedom and seven point integration was tested for analysis of thin plate and shell structures. The element formulation is based on the degenerate solid shell concept and the mixed formulation with assumed independent inplane and transverse shear strains. The effectiveness of the present modeling technique which features combined use of elements with kinematic modes and those without kinematic modes to eliminate both locking and spurious kinematic modes at the global structural levels are shown.

Description: The stress derivative technique for control of keystoning deformation in under-integrated finte elements is based on expansion of the stress in a Taylor series about the element center and retention of additional terms beyond the constant stress term. It has the advantage over other control techniques that keystoning resistance is provided by actual rather than artificial material properties. Application of this technique to the quadrilateral ring elements used for modelling solids of revolution subjected to axisymmetric loads is described. In a cylindrical coordinate system additional terms appear in the formulation which must be dealt with in arriving at a workable keystoning control scheme.

Description: An existing program is currently being adapted to perform finite element analysis by distributing substructures over a network of four Apple IIe microcomputers connected to a shared disk. In this network, one microcomputer controls the entire process while the others perform the analysis on each substructure in parallel. This substructure analysis is used in an iterative, fully stressed, structural resizing procedure. This procedure allows experimentatation with resizing in which all analyses are not completed during a single iteration. This research gives some insight on how to configure multidiscriplinary analysis and optimization procedures for decomposable engineering systems using either high performance engineering workstations or a parallel processor supercomputer. In addition, the operational experience gained facilitates the implementation of analysis programs on these new computers when they become available in an engineering environment.

Description: Described are the experiences gained from solving for the dynamic response of two simple structures on an experimental Multiple Instruction Multiple Data (MIMD) computer called the finite element machine. Introduced are MIMD computing concepts, describing how the concurrent algorithmic techniques implemented and giving results for the two example problems. The results show computational speedups of up to 7.83 using eight of the finite element machine processors and indicate that significant computational speedups are possible for large order structural computations.

Description: The FLEX/32 Multicomputer is a generic environment for cooperating multiple processors. The FLEX/32 supports a number of different processors, making it heterogeneous in terms of the instruction sets it supports, and homogeneous in its ability to provide consistent storage and input/output facilities to its differing processors. These facilities are accessed through standard 32-bit VMEbus connections. The FLEX/32 supports the full UNIX System V Operating System and languages associated with it, plus the extended ConCurrent C and Concurrent FORTRAN 77 languages that allow programming of concurrent software at a high level. Direct programming support at all levels is provided by the environment hardware for concurrent software execution and optimization, including hardware support for shared resource access arbitration, conditional critical region arbitration, and interprocessor messages.

Description: Research aimed at faster, more cost effective parallel machines and algorithms for improving designer productivity with finite element computations is discussed. A set of 8 boards, containing 4 nearest neighbor connected arrays of commercially available floating point chips and substantial memory, are inserted into a commercially available machine. One-tenth Mflop (64 bit operation) processors provide an 89% efficiency when solving the equations arising in a finite element problem for a single variable regular grid of size 40 by 40 by 40. This is approximately 15 to 20 times faster than a much more expensive machine such as a VAX 11/780 used in double precision. The efficiency falls off as faster or more processors are envisaged because communication times become dominant. A novel successive overrelaxation algorithm which uses cyclic reduction in order to permit data transfer and computation to overlap in time is proposed.

Description: Finite displacements of very shallow orthotropic shells and plates at small deformations

Description: Stability statistics in structural mechanics of plates and shells

Description: Contact problem for elastic foundations subjected to bending strains

Description: In the design of complex built-up structures that are made of truss, beam, membrane, shell, and solid, there are five different kinds of design variables: material property, sizing, shape, configuration, and topological variables. Previous research has shown that the improvement in performances obtained by altering the configuration of structural components can be much more significant than those obtained when the geometry is assumed to be fixed. Using the variational approach, a unified design sensitivity was developed for the first three kinds of design variables, and was further extended recently in many structural analysis problems such as nonlinear, structural dynamics, and frequency response analysis. A continuum design sensitivity analysis method is developed for the configuration design variable of built-up structures.

Description: Linear formulation of stress concentration near holes in various type shells

Description: Determination of moments, forces per unit length, and displacements characterizing state of stress and strain in two-dimensional helicoidal shell of uniform thickness

Description: Stress concentration around curvilinear holes in plates and shell, conformal mapping and other procedures for mathematical problems, and dynamic problems in bodies weakened by holes

Description: Stress concentration in shallow spherical shells determined by method with controlled boundary conditions

Description: Random vibration response data for orbiting geophysical observatory - comparative analysis of flight, acoustic, and vibration tests

Description: Flight performance of structural dampers on solar panels and waveguide antenna of Mariner IV space probe

Description: High temperature effects on compression spring characteristics

Description: Visual observation inspection device for bonded honeycomb sandwich structures using birefringent plastic coatings

Description: Fabrication and development of sandwich structure common bulkheads composed of aluminum alloy skins bonded to glass fiber honeycomb core for Saturn IV vehicles

Description: Design of heavy-gauge bonded honeycomb sandwich

Description: Thick honeycomb core sandwich structures for weight reduction in flight vehicles

Description: Engineering techniques for metal-to-metal and sandwich adhesive bonding

Description: Wind shear and turbulence effects on vertical slender structures

Description: Wind profiles and velocity distributions, and ground roughness effect on boundary layer

Description: Ground wind loading criteria for aerospace vehicle designing

Description: Aerodynamic excitation of tall structures by ground wind using wind tunnel models

Description: Ground wind oscillatory loads and launch vehicle configurations

Description: Elastic stability of slender bar under dynamic load for space rocket applications

Description: Digital computer automation of structural design process with application to launch vehicles and civil engineering structures

Description: Deperming of OGO satellite observatory

Description: High-impact resistant mechanisms for planetary landing of scientific spacecraft payloads

Description: Expandable boom from metal tape and fasteners for spacecraft applications

Description: Wind tunnel apparatus and atmospheric turbulence simulation for structural engineering

Description: Optimization of slosh baffle for large launch vehicle propellant tanks

Description: Linear and nonlinear inviscid supersonic flow in panel flutter research

Description: Magnetic shields of Mumetal for magnetometer sensors of Mariner Mars spacecraft

Description: Magnetic shielding of spacecraft relays and switches

Description: Evaluation of Mumetal shielding for reed switch relays

Description: Effects of geometrical variations in shape of explosive device on direction of projection of metal casing parts in guided missile warheads

Description: Meteoroid flux and puncture models for cislunar space

Description: Concentration of moments at holes in thin plates whose material follows nonlinear law of elasticity

Description: Temperature field and stress state for material nonuniformity caused by macroinclusions or cavities in structural elements of elastic bodies

Description: Perturbation criterion of inhomogeneous elastic body based on energy approach study of crack formation using concepts advanced by Griffith

Description: Propagation of elastic expansion wave from circular hole in cylindrical, anisotropic, inhomogeneous plate

Description: Dynamic stresses in unbound elastic space with cylindrical cavity

Description: Elastic wave propagation along cylindrical cavity filled with conducting fluid

Description: Fourier transform and theory of generalized functions to obtain differential equilibrium equations to determine concentrated bending and tangential forces on shallow shells

Description: Steady stress state and temperature field in thin thermoelastic plate weakened by circular hole

Description: Elimination of stress concentration around holes in laminated objects by use of reinforcement rings in design stage

Description: Stability of infinite plate weakened by hole with sharp notches

Description: Influence of stress concentration on supporting capacity of structural elements made of brittle microscopically nonuniform materials

Description: Nonlinear elastic plates weakened by square holes with reinforced corners

Description: Strain components on middle surface of thin shells determined by variational principles of linear moment theory

Description: Approximate solution for bending of long, narrow anisotropic plate moving in gas at supersonic velocity

Description: Critical loads calculated for cylindrical and conical shells that are simply-supported or have built-in ends

Description: Elastic-plastic stressed states for rapidly revolving disk of arbitrary profile in nonstationary temperature field