ALBERT

Description: The latest technology for structural analysis in relation to the design tasks that lie ahead for the space shuttle is reviewed. For shell-of-revolution structures, the analysis can be formulated as a one-dimensional problem which is readily solved by using finite-difference or numerical-integration techniques. For more general asymmetric shells, a two-dimensional formulation is required. However, the governing equations are readily formulated and are amenable to solution by finite-difference techniques. For a completely general structural arrangement, such as structural frameworks, recourse is usually made to discretized formulations using finite elements. Of course, the finite-element programs could be used for shell structures, but at a loss in accuracy and increase in computer time compared with the special purpose programs.

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: The design, manufacture, and testing of an isolation system placed under a bench which simulates a space platform during testing in a hard vacuum are described. A low natural frequency isolation system is used to decouple ground vibrations from the bench. The materials used are vacuum compatible and do not introduce unwanted matter which would contaminate the payload optics and/or the testing environment. The system accommodates payloads of varying weights and envelopes.

Description: The potential cost and performance advantages of welding was understood but ignored by solar panel manufacturers in the U.S. Although NASA, DOD and COMSAT have supported welding development efforts, soldering remains the only U.S. space qualified method for interconnecting solar cells. The reason is that no U.S. satellite prime contractor found it necessary, due to mission requirements, to abandon the space proven soldering process. It appears that the proposed NASA space station program will provide an array requirement, a 10 year operation in a low Earth orbital environment, that mandates welding. The status of welding technology in the U.S. is assessed.

Description: Researchers define damage-tolerant structural systems as those systems which not only have adeqate intact strength to withstand initial failure but also adequate residual strength to minimize the possibility of, and hence the consequences of, further failure. The incorporation of damage tolerance cannot be done in total isolation of the function being required of the system and the costs associated with obtaining improved damage tolerance. The approach, therefore, is to formulate multiple-objective, multi-level decision support problems (DSP), the solutions of which represent a compromise between higher costs and higher damage tolerance. Mulitple-objective decision support problems are easily solved in the linear domain. These formulations, however, include both linear and nonlinear constraints and goals, which in the past, have not been considered due to the resulting complexity. Here, researchers: (1) present a complete discussion and description of decision support problems; (2) identify what further research needs to be done in order to obtain information that is required but not known for solving problems using these models; and (3) identify what needs to be done to implement this prototype method in practice.

Description: Ocean-going vehicles and platforms are among the largest structures in the world and are subjected to relatively harsh conditions of motions and loads. Some of them, such as semi-submersible platforms, are a relatively new type of structure and hence there is no formal, well evolved and established structural design code as there is for more traditional structures. More recently, efforts have also been made to develop a design method of this type for ships and other ocean structures. One of the many advantages of a rationally based design method is versatility; it can be used for structures that have widely differing purposes, measures of merit, shapes and sizes. The purpose is to describe a rationally based design method that has been developed within the field of ocean structures, in order that persons dealing with other types of structure can judge whether and to what extent its various features may be useful for those other types. Also, even though some features may not be applicable they might stimulate some useful ideas.

Description: The use of microcomputers in the design of a cable catenary large space antenna system is discussed. The development of a system design capability, data base utilization, systems integration, program structure and logic, and integrated graphics output are discussed.

Description: An exceptionally elegant method for structural optimization with constraints on the static response presented by Shield and Prager is discussed. Their derivation of the optimality condition was facilitated by a reformulation of the structural elasticity equations in terms of what was then a new variational principle, the principle of stationary mutual potential energy. Their optimality condition relates the design variable to an appropriately defined mutual strain energy. An alternative but related approach, based upon the principle of stationary mutual complementary energy, presented by N. C. Haung, is also discussed. The simplicity of these principles lies in the facts that the energy functionals are stationary at the solution to the field equations and that their stationary value is proportional to the quantity to be optimized.

Description: The effects of space environments on damping materials and damping designs on flexible structures were investigated. The following items were examined: damping of flexible spacecraft appendages; composite loss factor (n sub s) vs. time in high vacuum for damped test beams and damping of flexible structures. The STEP experiments show inherent damping of flexible structures in space effective possible damping design configurations for space structures, effects of passively damped components on the system loss factor of flexible structures and the effect of space environment on properties of damping materials.

Description: Interplanetary dust particles (micrometeoroids) are expected to form well-defined craters upon impacting exposed material in space. Studying the frequency and features of these craters will provide data on the mass-flux distribution of micrometeoroids and, to a lesser extent, on the velocity magnitude and direction. This experiment will study impact craters produced by micrometeoroids on selected materials (metals and glasses in the form of thick targets) to obtain valuable technological and scientific data. Specifically, the studies will focus on determining micrometeoroid composition and mass-flux distribution. Analyses will also be made on the distribution of impact velocity vectors.

Description: The finite element method is extended to a broad class of practical nonlinear problems, treating both theory and applications from a general and unifying point of view. The thermomechanical principles of continuous media and the properties of the finite element method are outlined, and are brought together to produce discrete physical models of nonlinear continua. The mathematical properties of the models are analyzed, and the numerical solution of the equations governing the discrete models is examined. The application of the models to nonlinear problems in finite elasticity, viscoelasticity, heat conduction, and thermoviscoelasticity is discussed. Other specific topics include the topological properties of finite element models, applications to linear and nonlinear boundary value problems, convergence, continuum thermodynamics, finite elasticity, solutions to nonlinear partial differential equations, and discrete models of the nonlinear thermomechanical behavior of dissipative media.

Description: Description of an efficient digital computer method for the determination of the propagation of elastic stresses and deformations in certain geometrically nonlinear structures subjected to high impact loading. The finite-element matrix displacement approach utilizing curved quadrilateral shell elements in conjunction with a nodewise predictor-corrector method employing Runge-Kutta extrapolation techniques has been adopted for the present solution. The related computer program written in FORTRAN V for the UNIVAC 1108 computer has proved to be effective for the solution of a range of practical problems including rectangular and cylindrical panels. Numerical results are presented for a relevant structure, the cell container, and the negative electrode of an impact-resistant battery subjected to high impact, simulating its free landing on a planetary surface.

Description: The effect of slight meridional curvature on the buckling of shells subjected to an external constant directional lateral pressure is examined. A class of nearly cylindrical shells with constant meridional curvature, a central radius-to-thickness ratio of 500, and central rise-to-shell-length ratios varying from plus 5% to minus 5% is considered. A brief study is presented of the dependence of the critical lateral pressure on the meridional curvature and on the length and meridional edge restraint.

Description: The results are discussed of a study of the effect of slots on the buckling of cylinders with circular holes of systematically varied circularity. The results obtained suggest that the substantial increase in the buckling loads of cylindrical shells resulting from increases in slot length is due to a relief of the bending stress field.

Description: The paper defines the experimental phase of an objective to obtain the mechanical characteristics and coefficients required by the generalized form of Hooke's law for nylon-polyurethane-coated fabric. Test specimens were cylindrical fabric sleeves and were loaded in axial tension by an Instron, in hoop tension by pressurizing, and in shear by a torquing fixture. An extensive amount of strain data is included for a wide combination of the three membrane loads. The tests indicate highly nonlinear stress-strain characteristics of the fabric and a strong dependency on all three membrane loads.

Description: Concepts of sonic-boom pressure loading of building structures and the associated responses are reviewed, and results of pertinent theoretical and experimental research programs are summarized. The significance of sonic-boom load time histories, including waveshape effects, are illustrated with the aid of simple structural elements such as beams and plates. Also included are discussions of the significance of such other phenomena as three-dimensional loading effects, air cavity coupling, multimodal responses, and structural nonlinearities. Measured deflection, acceleration, and strain data from laboratory models and full-scale building tests are summarized, and these data are compared, where possible, with predicted values. Damage complaint and claim experience due both to controlled and uncontrolled supersonic flights over communities are summarized with particular reference to residential, commercial, and historic buildings. Sonic-boom-induced building responses are compared with those from other impulsive loadings due to natural and cultural events and from laboratory simulation tests.

Description: Study of the failure of glass windows, which is probably the most severe problem confronting the SST program outside of the psychological and sociological aspects. The structural configuration most likely to suffer a window failure due to sonic boom excitation is one representable by a single large room, a single large window, and an open door. It is suggested that if the open door were replaced with another window the resultant stresses in both windows would be substantially reduced.

Description: A geometrically nonlinear finite element analysis of cohesive failure in typical joints is presented. Cracked-lap-shear joints were chosen for analysis. Results obtained from linear and nonlinear analysis show that nonlinear effects, due to large rotations, significantly affect the calculated mode 1, crack opening, and mode 2, inplane shear, strain-energy-release rates. The ratio of the mode 1 to mode 2 strain-energy-release rates (G1/G2) was found to be strongly affected by the adhesive modulus and the adherend thickness. The ratios between 0.2 and 0.8 can be obtained by varying adherend thickness and using either a single or double cracked-lap-shear specimen configuration. Debond growth rate data, together with the analysis, indicate that mode 1 strain-energy-release rate governs debond growth. Results from the present analysis agree well with experimentally measured joint opening displacements. Previously announced in STAR as N83-13497

Description: Studies have been conducted to develop efficient techniques to simulate crack extension and to examine various local and global fracture criteria. Of the considered criteria, the crack-tip-opening angle (CTOA) or displacement (CTOD) at a specified distance from the crack tip was shown to be most suited for modeling stable crack growth and instability during the fracture process. The results obtained in a number of studies show the necessity for studying different crack configurations when assessing the validity of any fracture criteria. One of the objectives of the present investigation is related to a critical evaluation of the CTOD growth criterion using an elastic-plastic finite element analysis under monotonic loading to failure. The analysis was found to predict three stages of crack growth behavior under monotonic loading to failure. Calculated CTOD values agreed well with experimental values for crack growth initiation.

Description: Attenuation of the load transmitted to a vehicle occupant by the structure, either by modifying structural assembly, changing geometry of its elements, or adding specific load limiting devices to help dissipate kinetic energy was studied. General aviation aircraft have fuselage subfloors of a built up structure which are generally very stiff perpendicular to the floor of the cabin. The subfloor structure, designed to crush at an appropriate force level, can be used to advantage in a crash by dissipating energy through plastic buckling of the floor beams and frames. Simple closed form solutions to predict the mean crushing force levels of subfloor designs is useful in engineering practice. The crushing process of thin walled, plate formed, open structures with particular emphasis on L and cruciform shapes was analyzed. Lower and upper bound solutions for the mean crushing strength of cruciforms are obtained by considering modes of deformation which account for both bending and extensional deformation. It is the importance of extensional deformation to the energy absorption process, representing at least one third of the dissipated energy is shown. Previously announced in STAR as N82-33765

Description: The methods for predicting the strength of composite laminates with fastener holes are reviewed for the cases of unloaded as well as loaded holes. Three categories of methods are discussed: (1) the hole-boundary stress (HBS) methods which compare the peak stresses and the laminate strength; (2) the linear elastic fracture (LEFM) methods which assume cracks emanating from the hole and predict failure when the stress intensity factor for a crack equals the laminate fracture toughness; and (3) two different two-parameter (TP) methods: the average-stress (ASTP) method, which predicts failure when the average stress over a distance from the hole equals the laminate strength; and the point-stress (PSTP) Whitney-Nuismer (1974) approach, which predicts failure when the stress near the hole equals the material strength. The PSTP is relatively accurate and is by far the most widely used of all the prediction methods. However, for large ranges of variables and cases, the current strength predictions may not always be accurate.

Description: Optimality criteria methods take advantage of some concepts as those of statically determinate or indeterminate structures, and certain variational principles of structural dynamics, to develop efficient algorithms for the sizing of structures that are subjected to stiffness-related constraints. Some of the methods and iterative strategies developed over the last decade for calculations of the Lagrange multipliers in stressand displacement-limited problems, as well as for satisfying the appropriate optimality criterion, are discussed. The application of these methods are illustrated by solving problems with stress and displacement constraints.

Description: The stiffened composite structural panel analysis and sizing code designated 'PASCO' encompasses both the generality required for the exploitation of composite materials' design flexibility and an accurate buckling analysis for the detection of complex buckling modes. PASCO can accordingly design for buckling, frequency, material strength, and panel stiffness requirements. Attention is given to an additional thermal loading design capability. Design studies illustrate the importance of the multiple load condition capability when thermal loads are present.

Description: A detailed comparison is presented of the predicted eigenfrequencies of twisted rotating plates as obtained by using two different shape functions. Primarily, rotating twisted plates of two different aspect ratios and two different thickness ratios are considered. The effects of rotation are included by using a 'stress smoothing' technique when calculating the augmented stiffness matrix. In addition, the effects of Coriolis acceleration, contributions from membrane behavior, setting angle and sweep angle are considered. The effects of geometric nonlinearity are briefly discussed. Finally, results of a brief study of cambered plates are presented.

Description: Delamination growth in compressively loaded composite laminates was studied analytically and experimentally. The configuration used was a laminate with an across-the-width delamination. An approximate super-position stress analysis was developed to quantify the effects of various geometric, material, and load parameters on mode 1 and mode 2 strain energy release rates G sub 1 and G sub 2, respectively. Calculated values of G sub 1 and G sub 2 were then compared with measured cyclic delamination growth rates to determine the relative importance of G sub 1 and G sub 2. High growth rates were observed only when G sub 1 was large. However, slow growth was observed even when G sub 1 was negligibly small. This growth was apparently due to a large value of G sub 2.

Description: A special three-dimensional element based on the total Lagrangian description of the motion of a layered anisotropic composite medium is developed, validated and employed to analyze laminated anisotropic composite shells. The element contains the following features: geometric nonlinearity, dynamic (transient) behavior and arbitrary lamination scheme and lamina properties. Numerical results of nonlinear bending, natural vibration, and transient response are presented to illustrate the capabilities of the element.

Description: Solutions are given for several plane elastostatic problems of plates having a V-notch on one edge, and subjected to a variety of boundary conditions. The effect of the magnitude of the V-notch angle and specimen geometry on stress intensity factors KI and KII are obtained for unloaded notch surfaces. There is less than one per cent difference in opening model stress intensity factor in going from a zero degree notch angle to a 30 degree notch angle. Notch opening displacements at the plate edge were measured experimentally, and the results obtained were in excellent agreement with the computed results.

Description: Analysis of the panel flutter characteristics of a candidate thermal protection system (TPS) for the space shuttle, using piston theory aerodynamics and Lagrange equations. The results show the TPS candidate panel array to be deep in the 'no-flutter' region during launch and, therefore, safe from panel flutter.

Description: Consideration of Rayleigh wave effects in a homogeneous isotropic linearly elastic half-space subject to an impulsive uniform disk pressure loading. An approximate formula is obtained for the Rayleigh wave effects. It is shown that the Rayleigh waves near the center of loading arise from the portion of the dilatational and shear waves moving toward the axis, after they originate at the edge of the load disk. A study is made of the vertical displacement due to Rayleigh waves at points on the axis near the surface of the elastic half-space.

Description: Consideration of coupled stress waves generated by an impulsive load applied at one end of a semiinfinite plate. For the field equations governing the one-dimensional coupled waves a hyperbolic system of equations is obtained in which a strong coupling in the second derivatives exists. The method of characteristics described by Chou and Mortimer (1967) is extended to cover the case of strong coupling, and a study is made of the transient stress waves in a semiinfinite plate subjected to an initial step input. Coupled discontinuity fronts are found to propagate at different velocities. The normal plate stress and the bending moment at different time regimes are illustrated by graphs.

Description: The optical response to loading of a thin metallic film deposited on a low-modulus structural substrate is studied theoretically and experimentally. Two types of optical properties called total and central-image transmittance (or reflectance) are shown to be related to the mechanical state of the substrate. Empirical optical-mechanical relationships are proposed between these optical properties and the substrate strain field of a general plane-stress problem. A technique based on wrinkle and microfracture patterns is described for determining principal directions of strain. Experimental results for uniaxially loaded specimens show that it is possible to obtain a nearly linear relationship between transmittance and strain for certain materials combinations.

Description: A lumped-parameter model of a rectangular plate is developed by assuming fundamental mode solutions and using Hamilton's Principle and the Euler equations to set up the differential equation of motion for the system. The plate theory used may be described as the dynamic analogue of the von Karman large-deflection theory. Four sets of symmetrical boundary conditions are considered with the restriction of uniform pressure dynamic loads. The model takes the form of a mass on a cubic-hardening spring with each term defined by algebraic expressions of the plate parameters. The results for some specific problems are compared with two previous solutions. This method is less accurate but simpler to develop and apply.

Description: A generalized eigenvalue algorithm is presented herein along with the complete listing of the associated computer program, which may be conveniently utilized for the efficient solution of certain broad classes of eigenvalue problems. Extensive applications of the procedure are envisaged in the analysis of many important engineering problems, such as stability and natural frequency analysis of practical discrete structural systems, idealized by the finite element technique. The procedure based on the Sturm sequence method is accurate and fast, possessing several significant advantages over other known methods of such analysis. Numerical results are also presented for two representative structural engineering problems.

Description: An experimental and analytical investigation of the effect of a circular hole on the buckling of thin cylindrical shells under axial compression was carried out. The experimental results were obtained from tests performed on seamless electroformed copper shells and Mylar shells with a lap joint seam. These results indicated that the character of the shell buckling was dependent on a parameter which is proportional to the hole radius divided by the square root of the product of the shell radius and thickness. For small values of this parameter, there was no apparent effect of the hole on the buckling load. For slightly larger values of the parameter, the shells still buckled into a general collapse configuration, but the buckling loads were sharply reduced as the parameter increased. For still larger values of the parameter, the buckling loads were further reduced, and the shells buckled into a stable local buckling configuration.

Description: A boron-aluminum structural component, representative of Space Shuttle fuselage structure, has been fabricated and tested as part of a program to demonstrate boron-aluminum feasibility for Shuttle applications. The 4- by 5-foot component consists of eleven 24-ply Z section stringers and a +45 degree 8-ply skin. It was tested at room temperature and at 500 F under combined compressive and shear loadings. The components for this program were fabricated in order to encounter and solve scale-up problems, and to obtain test data representative of production hardware. The data from this program are expected to allow more accurate prediction of weights and costs of boron-aluminum structures for the Space Shuttle.

Description: An explicit formulation is developed to obtain the stored energy function, based on the Valanis-Landel separable symmetric stored energy function. This formulation is applicable to any multiaxial stress state, from simple tensile data alone, for those cases in which the stored energy function is a separable function of the stretch ratios. By direct comparison, it is also shown that styrene butadiene rubber, at least over a specific range, follows this postulated separable form.

Description: This paper gives some examples of projects underway. The projects are grouped according to the primary parameter being studied: loads, configuration, environment, or material. Significant topics include the effect of truncated loads, the effect of proof-tests, interference fasteners, accelerated and real-time tests for supersonic aircraft materials, crack-growth rates in built-up structures, and a model method by which material can be selected for a particular set of operating conditions.

Description: This paper deals with elastic stress intensity factors associated with stress fields near partially debonded interfaces in finite, bi-material plates. Bi-material stress intensity factors K1 and K2 are not associated, respectively, with symmetric and skew-symmetric loading and are functions of an arbitrary length parameter. Any fracture condition involving bi-material stress intensity factors must be independent of the arbitrary length parameter and must be consistent with experiment. We define a new parameter called the interfacial stress intensity factor which satisfies the former condition. A truncated Williams power series stress function for stresses in a cracked bi-material plate, as modified by Sih and Rice, is used with a modified boundary collocation procedure to generate interfacial stress intensity factors for finite bi-material plates.

Description: The results of an analytical study are presented on the use of boron polyimide filamentary composite material for the purpose of increasing the flutter speed of a simple titanium full depth sandwich wing structure designed for strength. The results clearly demonstrate that selective reinforcement of wing surfaces, using judiciously placed filamentary composites, promises sizable mass savings in the design of advanced aircraft structures.

Description: The problem of transient axisymmetric vibrations of thin circular elastic plates due to sonic boom excitation is investigated. The equation of motion for a solid circular plate is solved by applying the modified finite Hankel transform and the Laplace transform, and the numerical results are obtained with the help of the digital computer. From the analysis of the data, obtained for the dynamic deflections of the plates for the boom duration, it is concluded that for a normal flight the boom duration has a significant effect on the vibrations of plates as compared to the overpressure of the boom.

Description: Aluminum alloys of types 1100, 3003, 5052, 6061, and 2024 were used to study the sensitivity of the acousto-elastic constant to changes in the microstructure. Results show that there is a strong relationship between the acousto-elastic constants and the yield strength and hardness. This relationship depends on whether the alloy is strain hardened or precipitation hardened. In strain hardened alloys, the constants increase as the amount of solid solution is decreased, while the behavior is the opposite in precipitation hardened alloys.

Description: The effect of nonlinearity on the higher harmonic control (HHC) of helicopter vibration is investigated using a nonlinear aeroelastic simulation. A nonlinear solution is proposed which relates the HHC inputs to vibration outputs on the basis of a Volterra functional series. The Volterra series solution is shown to reduce to a vector polynomial equation relating HHC inputs to vibration outputs at any harmonic frequency. The nonlinear transfer relationship is identified from a nonlinear vibration analysis computer program, and the identification model is examined in detail. Improvements to current HHC algorithms are presented, and several Kalman filter divergence problems are quantified.

Description: A general crack opening stress equation is presented which may be used to correlate crack growth rate data for various materials and thicknesses, under constant amplitude loading, once the proper constraint factor has been determined. The constraint factor, alpha, is a constraint on tensile yielding; the material yields when the stress is equal to the product of alpha and sigma. Delta-K (LEFM) is plotted against rate for 2024-T3 aluminum alloy specimens 2.3 mm thick at various stress ratios. Delta-K sub eff was plotted against rate for the same data with alpha = 1.8; the rates correlate well within a factor of two.

Description: The transient response of a 0.4 by 0.6 cm rectangular bore rail accelerator was analyzed by a three dimensional finite element code. The copper rail deflected to a peak value of 0.08 mm in compression and then oscillated at an amplitude of 0.02 mm. Simultaneously the insulating side wall of glass fabric base, epoxy resin laminate (G-10) was compressed to a peak value of 0.13 mm and rebounded to a steady state in extension. Projectile pinch or blowby due to the rail extension or compression, respectively, can be identified by examining the time history of the rail displacement. The effect of blowby was most significant at the side wall characterized by mm size displacement in compression. Dynamic stress calculations indicate that the G-10 supporting material behind the rail is subjected to over 21 MPa at which the G-10 could fail if the laminate was not carefully oriented. Results for a polycarbonate resin (Lexan) side wall show much larger displacements and stresses than for G-10. The tradeoff between the transparency of Lexan and the mechanical strength of G-10 for sidewall material is obvious. Displacement calculations from the modal method are smaller than the results from the direct integration method by almost an order of magnitude, because the high frequency effect is neglected. Previously announced in STAR as N83-35412

Description: The unified constitutive theories for application to typical isotropic cast nickel base supperalloys used for air-cooled turbine blades were evaluated. The specific modeling aspects evaluated were: uniaxial, monotonic, cyclic, creep, relaxation, multiaxial, notch, and thermomechanical behavior. Further development of the constitutive theories to model thermal history effects, refinement of the material test procedures, evaluation of coating effects, and verification of the models in an alternate material will be accomplished in a follow-on for this base program.

Description: A state-of-the-art review of applicable constitutive models with selection of two for detailed comparison with a wide range of experimental tests was conducted. The experimental matrix contained uniaxial and biaxial tensile, creep, stress relaxation, and cyclic fatigue tests at temperatures to 1093 C and strain rates from .0000001 to .001/sec. Some nonisothermal cycles will also be run. The constitutive models will be incorporated into the MARC finite element structural analysis program with a demonstration computation made for advanced turbine blade configuration. In the code development work, particular emphasis is being placed on developing efficient integration algorithms for the highly nonlinear and stiff constitutive equations. Another area of emphasis is the appropriate and efficient methodology for determing constitutive constants from a minimum extent of experimental data.

Description: Static aerothermoelastic design/analysis of axial-flow compressors, modal flutter analysis of axial-flow turbomachines, forced vibration analysis of rotating cyclic structures and modal flutter analysis of advanced turbopropellers with highly swept blades are four new capabilities developed and implemented in NASTRAN Level 17.7. The contents, applicability and usefulness of these capabilities which were developed and documented under the sponsorship of NASA's Lewis Research Center are discussed. Overall flowcharts and selected examples are presented.

Description: COSMIC/NASTRAN was converted to the CRAY computer systems. The CRAY version is currently available and provides users with access to all of the machine independent source code of COSMIC/NASTRAN. Future releases of COSMIC/NASTRAN will be made available on the CRAY soon after they are released by COSMIC.

Description: Computer aided engineering systems are invaluable tools in performing NASTRAN finite element analysis. These techniques are implemented in both the pre-processing and post-processing phases of the NASTRAN analysis. The finite element model development, or pre-processing phase, was automated with a computer aided modeling program called Supertabl, and the review and interpretation of the results of the NASTRAN analysis, or post-processing phase, was automated with a computer aided plotting program called Output Display. An intermediate program, Nasplot, which was developed in-house, has also helped to cut down on the model checkout time and reduce errors in the model. An interface has been established between the finite element computer aided engineering system and the Learjet computer aided design system whereby data can be transferred back and forth between the two. These systems have significantly improved productivity and the ability to perform NASTRAN analysis in response to product development requests.

Description: A methodology that supports forced transient response dynamic solutions when both static and kinetic friction effects are included in a structural system model is described. Modifications that support this type of nonlinear transient response solution are summarized for the transient response dynamics (TRD) NASTRAN module. An overview of specific modifications for the NASTRAN processing subroutines, INITL, TRD1C, and TRD1D, are described with further details regarding inspection of nonlinear input definitions to define the type of nonlinear solution required, along with additional initialization requirements and specific calculation subroutines to successfully solve the transient response problem. The extension of the basic NASTRAN nonlinear methodology is presented through several stages of development to the point where constraint equations and residual flexibility effects are introduced into the finite difference Newmark-Beta recurrsion formulas. Particular emphasis is placed on cost effective solutions for large finite element models such as the Space Shuttle with friction degrees of freedom between the orbiter and payloads mounted in the cargo bay. An alteration to the dynamic finite difference equations of motion is discussed, which allows one to include friction effects at reasonable cost for large structural systems such as the Space Shuttle. Data are presented to indicate the possible impact of transient friction loads to the payload designer for the Space Shuttle. Transient response solution data are also included, which compare solutions without friction forces and those with friction forces for payloads mounted in the Space Shuttle cargo bay. These data indicate that payload components can be sensitive to friction induced loads.

Description: The results of a study to determine the accuracy of the three dimensional solid elements available in NASTRAN for predicting displacements is presented. Of particular interest in the study is determining how to effectively use solid elements in analyzing thick optical mirrors, as might exist in a large telescope. Surface deformations due to thermal and gravity loading can be significant contributors to the determination of the overall optical quality of a telescope. The study investigates most of the solid elements currently available in either COSMIC or MSC NASTRAN. Error bounds as a function of mesh refinement and element aspect ratios are addressed. It is shown that the MSC solid elements are, in general, more accurate than their COSMIC NASTRAN counterparts due to the specialized numerical integration used. In addition, the MSC elements appear to be more economical to use on the DEC VAX 11/780 computer.

Description: Some innovative techniques applicable to sensitivity analysis of discretized structural systems are reviewed. These techniques include a finite-difference step-size selection algorithm, a method for derivatives of iterative solutions, a Green's function technique for derivatives of transient response, a simultaneous calculation of temperatures and their derivatives, derivatives with respect to shape, and derivatives of optimum designs with respect to problem parameters. Computerized implementations of sensitivity analysis and applications of sensitivity derivatives are also discussed. Finally, some of the critical needs in the structural sensitivity area are indicated along with Langley plans for dealing with some of these needs.

Description: In developing the constitutive equations an interdisciplinary approach is being pursued. Specifically, both metallurgical and continuum mechanics considerations are recognized in the formulation. Experiments will be utilized to both explore general qualitative features of the material behavior that needs to be modeled and to provide a means of assessing the validity of the equations being developed. The model under development explicitly recognizes crystallographic slip on the individual slip systems. This makes possible direct representation of specific slip system phenomena. The present constitutive formulation takes the anisotropic creep theory and incorporates two state variables into the model to account for the effect of prior inelastic deformation history on the current rate-dependent response of the material.

Description: A brief overview of the status of the fatigue and fracture programs is given. The programs involve the development of appropriate analytic material behavior models for cyclic stress-strain-temperature-time/cyclic crack initiation, and cyclic crack propagation. The underlying thrust of these programs is the development and verification of workable engineering methods for the calculation, in advance of service, of the local cyclic stress-strain response at the critical life governing location in hot section compounds, and the resultant crack initiation and crack growth lifetimes.

Description: For economic reasons, the steam drive for a large compressor was replaced by a large synchronous electric motor. Due to the resulting large increase in mass and because the unit was mounted on a steel frame approximately 18 feet above ground level, it was deemed necessary to determine if a steady state or transient vibration problem existed. There was a definite possibility that a resonant or near resonant condition could be encountered. The ensuing analysis, which led to some structural changes as the analysis proceeded, did not reveal any major steady state vibration problems. However, the analysis did indicate that the system would go through several natural frequencies of the support structure during start-up and shutdown. This led to the development of special start-up and shutdown procedures to minimize the possibility of exciting any of the major structural modes. A coast-down could result in significant support structure and/or equipment damage, especially under certain circumstances. In any event, dynamic field tests verified the major analytical results. The unit has now been operating for over three years without any major vibration problems.

Description: Seismic analysis is a technique which pertains to loading described in terms of boundary accelerations. Earthquake shocks to buildings is the type of excitation which usually comes to mind when one hears the word seismic, but this technique also applied to a broad class of acceleration excitations which are applied at the base of a structure such as vibration shaker testing or shocks to machinery foundations. Four different solution paths are available in NASTRAN for seismic analysis. They are: Direct Seismic Frequency Response, Direct Seismic Transient Response, Modal Seismic Frequency Response, and Modal Seismic Transient Response. This capability, at present, is invoked not as separate rigid formats, but as pre-packaged ALTER packets to existing RIGID Formats 8, 9, 11, and 12. These ALTER packets are included with the delivery of the NASTRAN program and are stored on the computer as a library of callable utilities. The user calls one of these utilities and merges it into the Executive Control Section of the data deck to perform any of the four options are invoked by setting parameter values in the bulk data.

Description: A user's guide to the COSMIC/NASTRAN free field input for the Bulk Data section of the NASTRAN program is proposed. The free field input is designed to be user friendly and the user is not forced out of the computer system due to input errors. It is easy to use, with only a few simple rules to follow. A stand alone version of the COSMIC/NASTRAN free field input is also available. The use of free field input is illustrated by a number of examples.

Description: The annual static loading tests were made on the assembled wing/fuselage structure of the RP-1 sailplane. Simple wing bending and wing bending/torsion tests up to about four G's were made, using sand bag loading on the inverted aircraft. Bending a torsion deflections were measured and compared to those obtained in previous years' tests. It appears that the structure is retaining its original strength and stiffness very well. This aircraft has been disassembled and stored on the Jonsson Engineering Center under ambient, indoor conditions of temperature and humidity.

Description: Progress has been made on the design fabrication and installation of a tension carry-through linkage between the aft parts of the port and starboard wings to resist forward bending deflection on the RP-2 sailplane. In addition, towhook/release mechanisms, landing gear and wheel brake assembly, compound curved sandwich composite fuselage shells and an enclosed trailer to house and transport the complete disassembled aircraft, have all undergone various stages of design, fabrication, installation and testing.

Description: An algorithm is described for synthesizing the mass and stiffness matrices from experimentally derived modal data in a way that preserves the physical significance of the individual mass and stiffness elements. The mass and stiffness matrices are derived for a rollup solar array example, and are then used to define the modal response of a modified array.

Description: The initiation and growth of small cracks (5-500 microns m) from edge notches in 2024-T3 aluminum alloy sheets were studied under constant amplitude loading. Two methods were used to measure crack shape and size. In the first method, striation marker bands were periodically formed along the crack front by interrupting the constant amplitude loading by either an elevated R-ratio load sequence, or by an overload sequence. In the second method, the surface crack length was monitored by taking surface replicas at regular intervals. The marker band techniques did not provide reliable crack length and crack shape information for cracks smaller than 2 mm. The replica technique provided accurate information for surface crack length at all crack lengths, and fracture tests on specimens with small cracks provided crack-shape information. Crack growth rates were plotted against the stress intensity factor ranges. The results exhibited the small crack effect, in that the small cracks grew faster than large cracks at the same stress-intensity factor range. A crack closure model was also used to analyze the growth of small cracks from small (inclusion) defects at the notch surface.

Description: Adhesives are commonly characterized only by their static strength even though they are used in structural joints that are subjected to fatigue loads. This paper reviews the relationship between static and fatigue strength for four different specimen types: single-lap-shear, edge-delamination, double cantilever beam, and cracked-lap-shear. It was found that the ratio of static strength to fatigue strength varied from 2.3 to 4.7, depending on the adhesive and specimen configuration.

Description: The results of a three-dimensional, geometrically nonlinear, finite-element analysis of the bending of cross-ply laminated anisotropic composite plates are presented. Individual laminae are assumed to be homogeneous, orthotropic and linearly elastic. A fully three-dimensional isoparametric finite element with eight nodes (i.e. linear element) and 24 degrees of freedom (three displacement components per node) is used to model the laminated plate. The finite element results of the linear analysis are found to agree very well with the exact solutions of cross-ply laminated rectangular plates under sinusoidal loading. The finite element results of the three-dimensional, geometrically nonlinear analysis are compared with those obtained by using a shear deformable, geometrically nonlinear, plate theory. It is found that the deflections predicted by the shear deformable plate theory are in fair agreement with those predicted by three-dimensional elasticity theory; however stresses were found to be not in good agreement.

Description: The degree of insensitivity of the temperature dependence method to variations in preferred orientations in A508 low carbon steel specimens was determined. A pulsed phase locked loop was used to measure ultrasonic velocity in these specimens in the temperature range of 50 C below and above room temperature. Results are given in tables and indicate that the temperature dependence of ultrasonic longitudinal velocity in four specimens cut at different angles remains unchanged as a function of orientation within + or -2%. The variation of temperature dependence of ultrasonic shear velocity in the A508 steel as a function of applied stresses was also measured. Results show that values at zero stress are equal regardless of whether the polarization direction is parallel or perpendicular to the direction in which the stress is applied.

Description: Pin-loaded holes commonly occur in engineering structures. However, accurate analysis of such holes presents formidable difficulties because of the load-dependent contact of the pin with the plate. Significant progress has recently been achieved in the analysis of holes in isotropic plates. This paper develops a simple and accurate method for the partial contact analysis of pin-loaded holes in composites. The method is based on an inverse formulation that seeks to determine loads in a given contact-separation configuration. A unified approach for all types of fit was used. Continuum solutions were obtained for infinitely large plates of various typical orthotropic properties with holes loaded by smooth rigid pins. These solutions were then compared with those for isotropic plates. The effects of orthotropy and the type of fit were studied through load-contact relationships, distribution of stresses and displacements, and their variation with load. The results are of direct relevance to the analysis and design of pin joints in composite plates.

Description: This paper contains a description of the three-dimensional elasticity equations and the associated finite element model for natural vibrations of laminated rectangular plates. The numerical results for natural frequencies are compared with those obtained by a shear deformable plate theory. A number of cross-ply and angle-ply rectangular plates are analyzed for natural frequencies. For relatively thick plates, the plate element predicts frequencies higher than those predicted by the 3-D element.

Description: Nagtegaal and de Jong (1982) have studied stresses generated by simple finite shear in the case of elastic-plastic and rigid-plastic materials which exhibit anisotropic hardening. They reported that the shear stress is oscillatory in time. It was found that the occurrence of such an 'anomaly' is not restricted to anisotropic plasticity. Similar behavior in finite shear may result even in the case of hypoelasticity and classical isotropic hardening plasticity theory. The present investigation is concerned with the central problem of 'generalizing' with respect to the finite strain case, taking into account the constitutive relations of infinitesimal strain theories of classical plasticity with isotropic or kinematic hardening. The problem of hypoelasticity is also considered. It is shown that current controversies surrounding the choice of stress rate in the finite-strain generalizations of the constitutive relations and the anomalies surrounding kinematic hardening plasticity theory are easily resolvable.

Description: With the aim of applying the Dugdale model to the compact specimen, the equations for stress intensity factor and crack surface displacement are obtained for a pair of pin loads and a uniform stress acting on a segment of the crack surface. The plastic zone size is calculated from the 'finiteness' condition of Dugdale, and the results agree well with collocation results from Terada (1983). The load that causes incipient yielding at the compression point of the compact specimen is calculated from a finite element strip yield analysis for an elastic-perfectly plastic material. The crack tip opening displacement for the compact specimen is calculated by adding the displacement at the tip of the physical crack length due to the pin load and due to the uniform stress. The results are within 1.5 percent of current collocation results.

Description: This work summarizes a comprehensive study made of the free vibrations of twisted, cantilevered plates of rectangular planform. Numerous theoretical and experimental investigations previously made by others have resulted in frequency results which disagree considerably. To clarify the problem a joint industry/government/university research effort was initiated to obtain comprehensive theoretical and experimental results for models having useful ranges of aspect ratios, thickness ratios and twist angles. Theoretical data came from 19 independent computer analyses, including finite element, shell theory and beam theory idealizations. Two independent sets of experimental data were also obtained. The theoretical and experimental results are summarized and compared.

Description: A series of innovative designs and inventions which led to the solution of many aerospace vibration and shock problems through damping techniques is presented. The design of damped airborne structures has presented a need for such creative innovation. The primary concern was to discover what concepts were necessary for good structural damping. Once these concepts are determined and converted into basic principles, the design of hardware follows. The following hardware and techniques were developed in support of aerospace program requirements: shipping containers, alignment cables for precision mechanisms, isolation of small components such as relays and flight instruments, isolation for heavy flight equipment, coupling devices, universal joints, use of wire mesh to replace cable, isolation of 16-dB, 5000 lb horn, and compound damping devices to get better isolation from shock and vibration in a high steady environment.

Description: The vibrations of a cylindrical elastic shell filled with fluid and excited by a monopole acoustic source are studied. The modal mobilities of the shell wall at the source plane are examined for different radial locations of the source. The distribution of vibrational energy between the shell and the fluid field is calculated at various axial locations along the shell and for different radial source locations. An interchange of energy of vibration between the shell and the fluid as the wave field propagates along the shell-fluid system is observed.

Description: The application of substructuring techniques are described for two example problems, a square plate and the static analysis of a frame-wall interaction problem in multistory structures. A DMAP routine for calculating the multipoint constraint forces is presented.

Description: A nuclear reactor vessel is supported by a Z-ring and a box ring girder. The two proposed structural configurations to transmit the loads from the Z-ring and the box ring girder to the foundation are shown. The cantilever concrete ledge transmitting the load from the Z-ring and the box girder via the cavity wall to the foundation is shown, along with the loads being transmitted through one of the six steel columns. Both of these two supporting systems were analyzed by using rigid format 9 of NASTRAN for dynamic loads, and the thermal stresses were analyzed by AXISOL. The six column configuration was modeled by a combination of plate and bar elements, and the concrete cantilever ledge configuration was modeled by plate elements. Both configurations were found structurally satisfactory; however, nonstructural considerations favored the concrete cantilever ledge.

Description: NASTRAN was chosen as the principal analytical tool for structural analysis of the Illinois Center Plaza Hotel Building in Chicago, Illinois. The building is a 40-story, reinforced concrete structure utilizing a monolithic slab-column system. The displacements, member stresses, and foundation loads due to wind load, live load, and dead load were obtained through a series of NASTRAN runs. These analyses and the input technique are described.

Description: The extent to which suitable solutions may be obtained for one physics problem and two engineering type problems is traced. NASTRAN appears to be a practical tool to solve one-group steady-state neutron diffusion equations. Transient diffusion analysis may be performed after new levels that allow time-dependent temperature calculations are developed. NASTRAN piecewise linear anlaysis may be applied to solve those plasticity problems for which a smooth stress-strain curve can be used to describe the nonlinear material behavior. The accuracy decreases when sharp transitions in the stress-strain relations are involved. Improved NASTRAN usefulness will be obtained when nonlinear material capabilities are extended to axisymmetric elements and to include provisions for time-dependent material properties and creep analysis. Rigid formats 3 and 5 proved to be very convenient for the buckling and normal-mode analysis of a nuclear fuel element.

Description: The addition of discrete elements to NASTRAN poses significant interface problems with the level 15.1 assembly modules and geometry modules. Potential problems in designing new modules for higher-order parametric discrete elements are reviewed in both areas. An assembly procedure is suggested that separates grid point degrees of freedom on the basis of admissibility. New geometric input data are described that facilitate the definition of surfaces in parametric space.

Description: The NASTRAN program output is discussed from a structural analysts point of view, and the simple modifications which were made to the program in order to improve it are also described. The convenience of the output for use in original design work is critically appraised and compared with the output from ASTRAL. It is shown that considerable hand calculation is necessary in order to extract useful load distribution data from the available NASTRAN output. For this reason, some effort was directed toward providing additional force output for the NASTRAN shear panel element.

Description: The geometric and kinematic properties of the element are described, along with development of the matrices and vectors which characterize the element. The necessity for considering small deviations from element planeness is discussed, and the approach taken to account for these deviations is described. The improved accuracy over the existing quadrilateral membrane is indicated by a sample calculation for which an analytical solution is available from beam theory. For the same finite element idealization, the errors in maximum displacement and stress were significantly reduced.