ALBERT

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: A constitutive theory is presented for representing the anisotropic viscoplastic behavior of high-temperature alloys that possess directional properties resulting from controlled grain growth or solidification. The theory is an extension of a viscoplastic model that has been applied in structural analyses involving isotropic metals. Anisotropy is introduced through the definition of a vector field that identifies a preferential (solidification) direction at each material point. Following the development of a full multiaxial theory, application is made to homogeneously stressed elements in pure shear and to a uniaxially stressed rectangular block in plane stress with the stress direction oriented at an arbitrary angle with the material direction. It is shown that an additional material parameter introduced to characterize the degree of anisotropy can be determined on the basis of simple creep tests.

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: The force input mobility of an infinite elastic circular cylindrical shell filled with fluid is derived by using the spectral equations of motion. Mobilities are evaluated and their physical interpretations are discussed for a steel shell of thickness h/a = 0.05 filled with water and vibrating in the n = 0, 1 and 2 circumferential modes. The results are subsequently used to analyze the related situations of wave transmission through a radial ring constraint and the far field vibrational energy distributions between the contained fluid and the shell wall for line and point driving forces.

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: The boundary collocation method was used to generate Mode 1 stress intensity and crack mouth opening displacement coefficients for externally radially cracked ring segments subjected to three point radial loading. Numerical results were obtained for ring segment outer-to-inner radius ratios (R sub o/R sub i) ranging from 1.10 to 2.50 and crack length to segment width ratios (a/W) ranging from 0.1 to 0.8. Stress intensity and crack mouth displacement coefficients were found to depend on the ratios R sub o/R sub i and a/W as well as the included angle between the directions of the reaction forces. Previously announced in STAR as N83-35413

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: Two methods for identifying the mass, damping and stiffness matrices of a linear vibrating system are presented. Both methods require the measurement of acceleration, velocity and displacement at various locations of the system. In the first method, the response of the system subjected to known forces is used while the second method employs the free vibration data. The unknown parameters are recovered through the standard least squares procedure. Numerical results are presented for several examples.

Description: The problem of potential instability caused by actuator dynamics is examined. Two ways of overcoming it are discussed: rate feedback with compensation and position feedback with tuning filters.

Description: (Previously announced in STAR as N82-24501)

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: Problems encountered by aerospace designers in attempting to optimize whole aircraft are discussed, along with possible solutions. Large scale optimization, as opposed to component-by-component optimization, is hindered by computational costs, software inflexibility, concentration on a single, rather than trade-off, design methodology and the incompatibility of large-scale optimization with single program, single computer methods. The software problem can be approached by placing the full analysis outside of the optimization loop. Full analysis is then performed only periodically. Problem-dependent software can be removed from the generic code using a systems programming technique, and then embody the definitions of design variables, objective function and design constraints. Trade-off algorithms can be used at the design points to obtain quantitative answers. Finally, decomposing the large-scale problem into independent subproblems allows systematic optimization of the problems by an organization of people and machines.

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 simplified buckling analysis is presented for a family of periodic lattice structures such as those proposed for large space structures. A transcendental 6 x 6 matrix of eigenvalues is shown to be sufficient for modeling buckling behavior because member stiffnesses are based on an exact solution of the beam-column equation. Exact stiffnesses are derived for a curved member, thus allowing modeling of imperfect lattice structures. Comparisons of predictions of the lattice model with those available from shell and beam theory underscore the inaccuracies introduced by treating the lattice structure as a continuum. Sample calculations are provided for an isogrid cylinder and a three element double-laced truss.

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: Truss beams with members having viscous damping are modeled with a Timoshenko beam. Procedures for deriving the equivalent bending rigidity, transverse shear rigidity, and damping are presented. Explicit expressions for these equivalent beam properties are obtained for a specific truss beam. The beam model thus established is then used to investigate the effect of damping in free vibration. Finally, the beam is employed in the estimation of structural parameters in a simply-supported truss beam using a random search algorithm.

Description: A general approach for distributed parameter modeling of complex dynamical systems is described. The method consists of dividing the system in parts which can be modeled by simple partial differential equations and coupling the equations thus obtained by applying Hamilton's variational formalism to the entire system. The modeling of a large, offset-fed, wrap-rib antenna is presented to illustrate the approach. Although such models are perhaps not as precise as finite element models, they can be useful for initial physical insight and parametric design.

Description: The hybrid experimental-numerical stress-analysis technique, which saw limited applications during the 1950's, has been resurrected with the vastly improved numerical techniques of the 1970's. By inputing the experimental results as initial and boundary conditions, modern computer codes are executed in its generation and application modes to yield results which are unobtainable when only one of the two techniques is used. The hybrid technique thus exemplifies the complementary role of the experimental and numerical techniques.

Description: Fourier transforms and the Wiener-Hopf technique are used in conjunction with plane elastostatics to examine the singular crack tip stress field in the double cantilever beam (DCB) specimen. In place of the Dirac delta function, a family of functions which duplicates the important features of the concentrated forces without introducing unmanageable mathematical complexities is used as a loading function. With terms of order h-squared/a-squared retained in the series expansion, the dimensionless stress intensity factor is found to be K (h to the 1/2)/P = 12 to the 1/2 (a/h + 0.6728 + 0.0377 h-squared/a-squared), in which P is the magnitude of the concentrated forces per unit thickness, a is the distance from the crack tip to the points of load application, and h is the height of each cantilever beam. The result is similar to that obtained by Gross and Srawley by fitting a line to discrete results from their boundary collocation analysis.

Description: A fracture mechanics approach to the well-known delamination problem in composite materials is presented. Based on the theory of anisotropic laminate elasticity and interlaminar fracture mechanics concepts, the composite delamination problem is formulated and solved. The exact order of the delamination crack-tip stress singularity is determined. Asymptotic stress and displacement fields for an interlaminar crack are obtained. Fracture mechanics parameters such as mixed-mode stress intensity factors, KI, KII, KIII, and the energy release rate, G, for composite delamination problems are defined. To illustrate the fundamental nature of the delamination crack behavior, solutions for edge-delaminated graphite-epoxy composites under uniform axial extension are presented. Effects of fiber orientation, ply thickness, and delamination length on the interlaminar fracture are examined.

Description: Previously cited in issue 13, p. 2112, Accession no. A82-30182

Description: Stress-intensity factors are obtained for point loaded equal length cracks emanating from a circular hole in an infinite plate. A series approach and the Muskhelishvili formulation in the two-dimensional theory of elasticity are used to derive the solution. The applicability of the solution is demonstrated by using it as a Green's function to obtain stress-intensity factors in the case of (1) biaxial tension and pure shear of an infinite plate and (2) tension and pin loading of a plate with cracks emanating from one hole in a row of holes.

Description: Previously cited in issue 13, p. 2111, Accession no. A82-30171

Description: Previously cited in issue 13, p. 2110, Accession no. A82-30166

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: The objectives of the hot section technology (HOST) burner liner cyclic rig test program are basically threefold: (1) to assist in developing predictive tools needed to improve design analyses and procedures for the efficient and accurate prediction of burner liner structural response; (2) to calibrate, evaluate and validate these predictive tools by comparing the predicted results with the experimental data generated in the tests; and (3) to evaluate existing as well as advanced temperature and strain measurement instrumentation, both contact and noncontact, in a simulated engine cycle environment. The data generated will include measurements of the thermal environment (metal surface temperatures) as well as structural (strain) and life (fatigue) responses of simulated burner liners and specimens under controlled boundary and operating conditions. These data will be used to calibrate, compare and validate analytical theories, methodologies and design procedures, as well as improvements in them, for predicting liner temperatures, stress-strain responses and cycles to failure. Comparison of predicted results with experimental data will be used to show where the predictive theories, etc. need improvements. In addition, as the predictive tools, as well as the tests, test methods, and data acquisition and reduction techniques, are developed and validated, a proven, integrated analysis/experiment method will be developed to determine the cyclic life of a simulated burner liner.

Description: The HOST Liner Cyclic Program is utilizing two types of test apparatus, rectangular box rigs and a full annular rig. To date two quartz lamp cyclic box rigs have been tested and a third is to begin testing in late October 1983. The box rigs are used to evaluate 5x8 inch rectangular linear samples. A 21 inch diameter outer liner simulator is also being built up for testing beginning in April 1984. All rigs are atmospheric rigs. The first box rig, a three 6-kVA lamp installation, was operated under adverse conditions to determine feasibility of using quartz lamps for cyclic testing. This work was done in December 1981 and looked promising. The second box rig, again using three 6-kVA lamps, was operated to obtain instrumentation durability information and initial data input to a Finite Element Model. This limited test program was conducted in August 1983. Five test plates were run. Instrumentation consisted of strain gages, thermocouples and thermal paint. The strain gages were found to fail at 1200 F as expected though plates were heated to 1700 F. The third box rig, containing four 6-kVA lamps, is in build up for testing to begin in late October 1983. In addition to 33 percent greater power input, this rig has provision for 400 F backside line cooling air and a viewing port suitable for IR camera viewing. The casing is also water cooled for extended durability.

Description: The objective of this program is to produce a series of new computer codes that permit more accurate and efficient three dimensional inelastic analysis of selected hot section components - combustor liners, turbine blades and turbine vanes. The computer codes embody a progression of mathematical models and are streamlined to take advantage of geometrical features, loading conditions, and forms of material response that distinguish each group of selected components.

Description: The objective of this research is to develop an analytical tool capable of economically evaluating the cyclic time dependent plasticity which occurs in hot section engine components in areas of strain concentration resulting from the combination of both mechanical and thermal stresses. The techniques developed must be capable of accommodating large excursions in temperatures with the associated variations in material properties including plasticity and creep. The overall objective of this proposed program is to develop advanced 3-D inelastic structural/stress analysis methods and solution strategies for more accurate and yet more cost effective analysis of combustors, turbine blades, and vanes. The approach will be to develop four different theories, one linear and three higher order with increasing complexities including embedded singularities.

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: The dynamic characteristics of the Solar Array Flight Experiment (SAFE) structure during deployment and retraction are investigated. The SAFE structure consists of a deployable mast with an attached solar blanket designed with accordion type folds to permit packaging in a small volume. The planar form of the blanket geometry during deployment is maintained by a blanket tension/guidewire system. Structurally, the mast is modeled as an Euler beam column with inplane and out of plane bending and finite torsional stiffness. For out of plane motion, the blanket is modeled as a distributed mass uniformly supported by the three guidewires. For inplane motion the blanket displacements are assumed to vary linearly from the mast base to the mast tip. The mathematical model uses a virtual work formulation, required because the axial loading on the mast is nonconservative, combined with assumed beam modes to derive the differential equations of motion. Consideration of the time dependent boundary conditions results in an infinite set of ODE with time dependent coefficients. Finally, correlation of mast tip accelerations to mast base bending moments for specified modal motions are indicated.

Description: Traditional composite lamination theory was used to predict composite graphite-epoxy laminate stiffnesses for comparison to quasi-experimental stiffnesses developed from cylindrical bottle pressure testing. Stiffness sensitivities were examined for variations in constituent materials and geometric properties. The material component interactions examined were the fiber longitudinal and transverse Young's modulus, the fiber shear modulus and primary Poisson's ratio, the resin Young's modulus, shear modulus, and Poisson's ratio. The geometric variation of the helical winding angle was also examined. Two computer programs were written to generate the data used to demonstrate the stiffness variations.

Description: An experiment to measure the elastic deflection of the DSS 14 concrete pedestal under the weight of the antenna was conducted in February 1983 and is compared to a similiar experiment made in 1968. Comparison of the results confirms the decrease in elastic modulus measured on core samples recently taken from the pedestal.

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: An unconditionally stable implicit-explicit method is proposed for the analysis of transient coupled thermal stress waves, and the computer-implementation aspects of the method are discussed. In the method proposed here, the mechanical displacement and temperature are used as independent variables. The resulting coupled finite element matrix equations are symmetric.

Description: Previously cited in issue 13, p. 2112, Accession no. A82-30186

Description: Previously cited in issue 13, p. 2108, Accession no. A82-30124

Description: Previously cited in issue 01, p. 60, Accession no. A83-10900

Description: Previously cited in issue 13, p. 2109, Accession no. A82-30152

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: A fracture mechanics analysis has been developed that describes the progress of delamination damage in composite plates struck by a hard spherical object. The analysis is based on large deflection plate mechanics for circular isotropic plates wherein multiple axisymmetric delaminations grow. Test data show that the analysis predicts the influence of plate thickness, support conditions, and matrix toughness on the onset and propagation of delamination.

Description: A three-dimensional element based on the total Lagrangian description of large deformations of a layered anisotropic composite medium is developed, validated, and used to analyze layered 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: Stiffened panels are widely used in aircraft structures such as wing covers, fuselages, control surfaces, spar webs, bulkheads, and floors. The detailed sizing of minimum-weight stiffened panels involves many considerations. Use of composite materials introduces additional complexities. Many potential modes of failure exist. Analyses for these modes are often not trivial, especially for those involving large out-of-plane displacements. Accurate analyses of all potential failure modes are essential. Numerous practical constraints arise from manufacturing/cost considerations and from damage tolerance, durability, and stiffness requirements. The number of design variables can be large when lamina thicknesses and stacking sequence are being optimized. A significant burden is placed on the sizing code due to the complex analyses, practical constraints, and number of design variables. On the other hand, sizing weight-efficient panels without the aid of an automated procedure is almost out of the question. The sizing code postbuckled Open-Stiffener Optimum Panels (POSTOP) has been developed to aid in the design of minimum-weight panels subject to the considerations mentioned above. Developed for postbuckled composite panels, POSTOP may be used for buckling resistant panels and metallic panels as well. The COPES/CONMIN optimizer is used in POSTOP although other options such as those in the ADS system could be substituted with relative ease. The basic elements of POSTOP are shown. Some of these elements and usage of the program are described.