ALBERT

Description: The rotational stiffness of hinge joints, and the gap of the joints applied in large deployable trusses, have been experimentally shown to have a significant role in determining such structures' dynamic behavior; an analytical validation of these results is presented for the case where linear rotation springs are used to model the hinge joints employed in a simple beam in trusses. The results obtained indicate that the natural frequencies of these structures depend not only on joint stiffness but also on joint location. Such gap parameters as gap size, stiffness, position, and excitation-force levels, are discussed with a view to a deeper understanding of their effects on a space interferometry system's simulated dynamic responses.

Description: Microfracture (fiber/matrix fracture, interphase debonding and interply delamination) in high temperature metal matrix composites (HTMMC), subjected to thermal loading, is computationally simulated. Both unidirectional and crossply SiC/Ti15 composites are evaluated for microfracture driven by thermal loads, using multicell finite element models. Results indicate that under thermal loads alone, microfracture propagation is not as sensitive as it is under mechanical loads.

Description: This paper presents results from a finite element micromechanics analysis of thermally induced stresses in composites at cryogenic temperatures typical of spacecraft operating environments. The influence of microstructural geometry, constituent and interphase properties, and laminate orientation were investigated. Stress field results indicated that significant matrix stresses occur in composites exposed to typical spacecraft thermal excursions; these stresses varied with laminate orientation and circumferential position around the fiber. The major difference in the predicted response of unidirectional and multidirectional laminates was the presence of tensile radial stresses, at the fiber/matrix interface, in multidirectional laminates with off-axis ply angles greater than 15 deg. The predicted damage initiation temperatures and modes were in good agreement with experimental data for both low (207 GPa) and high (517 GPa) modulus carbon fiber/epoxy composites.

Description: The structural integrity of proposed high speed aircraft can be seriously affected by the extremely high surface temperatures and large temperature gradients throughout the vehicle's structure. Variations in the structure's elastic characteristics as a result of thermal effects can be observed by changes in vibration frequency, damping, and mode shape. Analysis codes that predict these changes must be correlated and verified with experimental data. The experimental modal test techniques and procedures used to conduct uniform, nonuniform, and transient thermoelastic vibration tests are presented. Experimental setup and elevated temperature instrumentation considerations are also discussed. Modal data for a 12 by 50 inch aluminum plate heated to a temperature of 475 F are presented. These data show the effect of heat on the plate's modal characteristics. The results indicated that frequency decreased, damping increased, and mode shape remained unchanged as the temperature of the plate was increased.

Description: This paper presents two approaches that will improve the accuracy of explicit approximations of frequency responses. They are important components of structural design optimization and system identification. The first is a product form approximation that improves approximated frequency responses near resonant conditions. This approach is applicable to both design optimization and system identification problems. The second is aimed at improving the computed sensitivity coefficients that are used to solve system identification problems.

Description: An experimental technique was used to measure structural intensity through an aircraft fuselage with an excitation load applied near one of the wing attachment locations. The fuselage was relatively large, requiring several measurement locations to analyze the intensity flow through the whole of the structure. For the measurement of structural intensity, the use of a transducer array was necessary at every location of interest. A trade-off was therefore required between the number of measurement transducers, the mounting of these transducers, and the accuracy of the measurements. Using four accelerometers mounted on a bakelite platform, structural intensity vectors were measured at locations distributed throughout the fuselage. The results of these measurements, together with a discussion on the suitability of the approach for measuring structural intensity on a real structure, are presented in this paper.

Description: The concept of pressure proof testing of fuselage structures with fatigue cracks to insure structural integrity was evaluated from a fracture mechanics viewpoint. A generic analytical and experimental investigation was conducted on uniaxially loaded flat panels with crack configurations and stress levels typical of longitudinal lap splice joints in commercial transport aircraft fuselages. The results revealed that the remaining fatigue life after a proof cycle was longer than that without the proof cycle because of crack growth retardation due to increased crack closure. However, based on a crack length that is slightly less than the critical value at the maximum proof stress, the minimum assured life or proof test interval must be no more than 550 pressure cycles for a 1.33 proof factor and 1530 pressure cycles for a 1.5 proof factor to prevent in-flight failures.

Description: The experiment is aimed at controlling the boundary layer transition location and the plate vibration when excited by a flow and an upstream sound source. Sound has been found to affect the flow at the leading edge and the response of a flexible plate in a boundary layer. Because the sound induces early transition, the panel vibration is acoustically coupled to the turbulent boundary layer by the upstream radiation. Localized surface heating at the leading edge delays the transition location downstream of the flexible plate. The response of the plate excited by a turbulent boundary layer (without sound) shows that the plate is forced to vibrate at different frequencies and with different amplitudes as the flow velocity changes indicating that the plate is driven by the convective waves of the boundary layer. The acoustic disturbances induced by the upstream sound dominate the response of the plate when the boundary layer is either turbulent or laminar. Active vibration control was used to reduce the sound induced displacement amplitude of the plate.

Description: A spectral analysis method is presented for the extraction of frequency dependent strain energy release rates and stress intensity factors. The approach is based on the crack closure technique and is formulated directly in the frequency domain. Because of this, it is computationally more efficient than an alternative time domain approach for systems exhibiting peak behavior. The method is demonstrated for a center cracked panel subject to static in-plane and random acoustic loading.

Description: A finite element approach is developed for beam type pipes undergoing large deflections subjected to random loadings. The influence of fluid velocity on the random response is investigated. The root-mean-square (rms) deflections and frequencies for different sound spectrum level values are determined for pipes with both ends either simply supported or clamped. The required number of modes to achieve accurate rms deflections is studied. The prediction of fatigue life is then based on the maximum rms stress. This analytical investigation will help to broaden the basic understanding of the role of fluid flow within structures subjected to random excitations.

Description: The equilibrium equations and the compatibility conditions are fundamental to the analyses of structures. However, anyone who undertakes even a cursory generic study of the compatibility conditions can discover, with little effort, that historically this facet of structural mechanics had not been adequately researched by the profession. Now the compatibility conditions (CC's) have been researched and are understood to a great extent. For finite element discretizations, the CC's are banded and can be divided into three distinct categories: (1) the interface CC's, (2) the cluster or field CC's, and (3) the external CC's. The generation of CC's requires the separating of a local region, then writing the deformation displacement relation (ddr) for the region, and finally, the eliminating of the displacements from the ddr. The procedure to generate all three types of CC's is presented and illustrated through examples of finite element models. The uniqueness of the CC's thus generated is shown. The utilization of CC's has resulted in the novel integrated force method (IFM). The solution that is obtained by the IFM converges with a significantly fewer number of elements, compared to the stiffness and the hybrid methods.

Description: After briefly describing the principles of frozen stress photoelastic and moire interferometric analyses, and the corresponding algorithms for converting optical data from each method into stress intensity factors (SIF), the methods are applied to the determination of crack shapes, SIF determination, crack closure displacement fields, and pre-crack damage mechanisms in typical aircraft component configurations.

Description: The bounds on the equivalent elastic material properties of a composite are presently addressed by a unified energy approach which is valid for both unidirectional and 2D and 3D woven composites. The unit cell considered is assumed to consist, first, of the actual composite arrangement of the fibers and matrix material, and then, of an equivalent pseudohomogeneous material. Equating the strain energies due to the two arrangements yields an estimate of the upper bound for the material equivalent properties; successive increases in the order of displacement field that is assumed in the composite arrangement will successively produce improved upper bound estimates.

Description: The evolution of high-temperature, creep-fatigue, life-prediction methods used for cyclic crack initiation is traced from inception in the late 1940s. The methods reviewed are material models as opposed to structural life prediction models. Material life models are used by both structural durability analysts and by material scientists. The latter use micromechanistic models as guidance to improve a material's crack initiation resistance. Nearly one hundred approaches and their variations have been proposed to date. This proliferation poses a problem in deciding which method is most appropriate for a given application. Approaches have been identified as being combinations of fourteen different classifications. This review is intended to aid both developers and users of high-temperature fatigue life prediction methods by providing a background from which choices can be made.

Description: State-of-the-art nonlinear finite element analysis techniques are evaluated by applying them to a realistic aircraft structural component. A wing panel from the V-22 tiltrotor aircraft is chosen because it is a typical modern aircraft structural component for which there is experimental data for comparison of results. From blueprints and drawings, a very detailed finite element model containing 2284 9-node Assumed Natural-Coordinate Strain elements was generated. A novel solution strategy which accounts for geometric nonlinearity through the use of corotating element reference frames and nonlinear strain-displacement relations is used to analyze this detailed model. Results from linear analyses using the same finite element model are presented in order to illustrate the advantages and costs of the nonlinear analysis as compared with the more traditional linear analysis.

Description: As a means of determining a stress intensity factor solution, the compliance properties of an ARALL-2 laminated-sheet composite were investigated. Fatigue crack growth rate (FCGR) tests were conducted on middle crack tension (MT) specimens fabricated from a layup consisting of three sheets of 2024-T3 aluminum bonded together with unidirectional aramid fibers embedded in epoxy. Excellent fatigue crack growth properties are obtained by the presence of unbroken aramid fibers in the wake of the crack tip. These unbroken fibers act as a bridging mechanism to inhibit further crack growth. To quantify the effect of maximum fatigue load on compliance, a series of FCGR tests were performed. Effective crack lengths were determined to be at least 10 mm shorter than surface measured crack lengths for a 76-mm-wide specimen. The bridging zone was estimated to be at least 5 mm. Compliance and stress intensity factor as functions of effective crack length were determined.

Description: Tapered (0 deg) laminates of S2/CE9000 and S2/SP250 glass/epoxies, and IM6/1827I graphite/epoxy were tested in cyclic tension. The specimens usually showed some initial stable delaminations in the tapered region, but these did not affect the stiffness of the specimens, and loading was continued until the specimens either delaminated unstably, or reached 10(exp 6) to 2 x 10 (exp 7) million cycles with no unstable delamination. The final unstable delamination originated at the junction of the thin and tapered regions. A finite-element model was developed for the tapered laminate with and without the initial stable delaminations observed in the tests. The analysis showed that for both cases the most likely place for an opening (Mode I) delamination to originate is at the junction of the taper and thin regions. For each material type, the models were used to calculate the strain energy release rate, G, associated with delaminations originating at that junction and growing either into the thin region or tapered region. For the materials tested, cyclic G(sub Imax) values from DCB tests were used with the maximum strain energy release rates calculated from the finite-element analysis to predict the onset of unstable delamination at the junction as a function of fatigue cycles. The predictions were compared to experimental values of maximum cyclic load as a function of cycles to unstable delamination from fatigue tests in tapered laminates. For the IM6/1827I and S2/SP250 laminates, the predictions agreed very well with the test data. Predicted values for the S2/CE9000 were conservative compared to the test data.

Description: The present study proposes a mathematical model utilizing the internal state variable concept for predicting the upper bounds of the reduced axial and shear stiffnesses in cross-ply laminates with matrix cracks. The displacement components at the matrix crack surfaces are explicitly expressed in terms of the observable axial and shear strains and the undamaged material properties. The reduced axial and shear stiffnesses are predicted for glass/epoxy and graphite/epoxy laminates. Comparison of the model with other theoretical and experimental studies is also presented to confirm direct applicability of the model to angle-ply laminates with matrix cracks subjected to general in-plane loading.

Description: A simple strength of materials analysis was developed for a double-cantilever beam (DCB) specimen to account for geometric nonlinearity effects due to large deflections and T-tabs. A new DCB data analysis procedure was developed to include the effects of these nonlinearities. The results of the analysis were evaluated by DCB tests performed for materials having a wide range of toughnesses. The materials used in the present study were T300/5208, IM7/8551-7, and AS4/PEEK. Based on the present analysis, for a typical deflection/crack length ratio of 0.3 (for AS4/PEEK), T-tabs and large deflections cause a 15 and 3 percent error, respectively, in the computer Mode I strain energy release rate. Design guidelines for DCB specimen thickness and T-tab height were also developed in order to keep errors due to these nonlinearities within 2 percent. Based on the test results, for both hinged and tabbed specimens, the effects of large deflection on the Mode I fracture toughness (G sub Ic) were almost negligible (less than 1 percent) in the case of T300/5208 and IM7/8551-7; however, AS4/PEEK showed a 2 to 3 percent effect. The effects of T-tabs G sub Ic were more significant for all the materials with T300/5208 showing a 5 percent error, IM7/8551-7 a 15 percent error, and, AS4/PEEK a 20 percent error.

Description: A test rig for testing a thick split cantilever beam for scissoring delamination (mode 3) fracture toughness was developed. A 3-D finite element analysis was conducted on the test specimen to determine the strain energy release rate, G, distribution along the delamination front. The virtual crack closure technique was used to calculate the G components resulting from interlaminar tension, GI, interlaminar sliding shear, GII, and interlaminar tearing shear, GIII. The finite element analysis showed that at the delamination front no GI component existed, but a GII component was present in addition to a GIII component. Furthermore, near the free edges, the GII component was significantly higher than the GIII component. The GII/GIII ratio was found to increase with delamination length but was insensitive to the beam depth. The presence of GII at the delamination front was verified experimentally by examination of the failure surfaces. At the center of the beam, where the failure was in mode 3, there was significant fiber bridging. However, at the edges of the beam where the failure was in mode 3, there was no fiber bridging and mode 2 shear hackles were observed. Therefore, it was concluded that the split cantilever beam configuration does not represent a pure mode 3 test. The experimental work showed that the mode 2 fracture toughness, GIIc, must be less than the mode 3 fracture toughness, GIIIc. Therefore, a conservative approach to characterizing mode 3 delamination is to equate GIIIc to GIIc.

Description: Thermal deformations and stresses were studied in a silicon-carbide/aluminum filamentary composite at temperatures up to 370 C (700 F). Longitudinal and transverse thermal strains were measured with strain gages and a dilatometer. An elastoplastic micromechanical analysis based on a one-dimensional rule-of-mixtures model and an axisymmetric two-material composite cylinder model was performed. It was established that beyond a critical temperature thermal strains become nonlinear with decreasing longitudinal and increasing transverse thermal-expansion coefficients. This behavior was attributed to the plastic stresses in the aluminum matrix above the critical temperature. An elastoplastic analysis of both micromechanical models was performed to determine the stress distributions and thermal deformation in the fiber and matrix of the composite. While only axial stresses can be determined by the rule-of-mixtures model, the complete triaxial state of stress is established by the composite cylinder model. Theoretical predictions for the two thermal-expansion coefficients were in satisfactory agreement with experimental results.

Description: Some of the basic issues of ratchetting behavior that are being addressed by the viscoplastic modeling community are discussed. Some of the shortcomings of existing viscoplastic models are examined in the light of the difficulty involved in using established viscoplastic modeling techniques to predict ratchetting accurately.

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

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

Description: The present volume on impact damage in composite plates presents an extensive compendium of visual and graphic data regarding a variety of material and impactor parameters. The photographs are taken with X-ray and C-scan imaging in conjunction with a dye penetrant to show matrix cracks and delaminations. Impact and static-loading tests are performed on plates of graphite-epoxy, graphite-toughened epoxy, and graphite-PEEK materials. The images are presented to yield specific visual data regarding such parameters as impactor velocity, thickness of the back ply group, impactor nose radius, and the effects of multiple delaminations. The images are grouped in eight subsets that correspond to parameters including plate length, material, and the difference in fiber orientation between adjacent ply groups. This substantial volume represents a systematic effort to study the effects of several material parameters on impact damage.

Description: A thermodynamic foundation using the concept of internal state variables is given for a general theory of viscoplasticity for initially isotropic materials. Three, fundamental, internal, state variables are admitted; they are: a tensorial back stress for kinematic effects, and scalar drag and yield strengths for isotropic effects. All three are considered to evolve phenomenologically according to competitive processes between strain hardening, deformation induced dynamic recovery, and thermally induced static recovery. Within this phenomenological framework, a thermodynamically admissible set of evolution equations is proposed. The theory allows each of the three internal variables to be composed as a sum of independently evolving constituents. The evolution of internal state can also include terms that vary linearly with the external variable rates, whose presence affects the energy dissipation properties of a material.

Description: The bearing matrix formulations proposed by Lim and Singh (1990) are extended to analyze the overall dynamics of a geared rotor system which includes a spur gear pair, shafts, rolling-element bearings, a motor, a load, a casing, and flexible or rigid mounts. For this purpose, discrete vibration models are developed and used to predict vibration transmission through the bearings and to investigate the effects of the bearing, casing, and mount dynamics on the dynamic characteristics of the internal rotating system. Analytical predictions show that the theory is capable of predicting the bearing and mount moment transmissibilities in addition to the force transmissibilities. The predicted flexural vibrations of the casing plate are in good agreement with measurements conducted on an experimental set-up that consisted of a high-precision beam and pinion, and four identical rolling element bearings contained in a flexible casing mounted rigidly on a massive foundation.

Description: It is pointed out that real parametric modal frequency and damping variation for lightly damped systems does not resemble disks on the complex plane. Using complex uncertainty (and therefore disklike) models can introduce conservativeness in the design method. An alternative means of developing suitable complex uncertainty models is presented. It involves treating the uncertainty as perturbations to the system eigenvalues and using a particular linear fractional transformation to cover this uncertainty. This approach is applicable to modes within the bandwidth of control. It can be used in conjunction with the standard perturbation modeling approaches. A simple SISO (single-input single-output) example, motivated by a flexible truss experiment at the Jet Propulsion Laboratory, is discussed in order to illustrate the proposed approach.

Description: Using an approach based on a time domain analysis, a method of equivalent linearization is applied for an estimation of the strain response of complex nonlinear structures having nearly arbitrary complexity. Fatigue lives estimated for a nonlinear beam with random excitation using the approximate method were compared with results obtained using a conventional numerical simulation, yielding nearly identical results.

Description: Most methods of system identification of large flexible structures by far are based on the lumped parameter approach. Because of the considerable computational burden due to the large number of unknown parameters, distributed parameter approach, which greatly decreases the number of unknowns, has being investigated. In this paper a distributed parameter model for the estimation of modal characteristics of NASA Mini-Mast truss has been formulated. Both Bernoulli-Euler beam and Timoshenko beam equations are used to characterize the lateral bending vibrations of the truss. The measurement of the lateral displacement at the tip of the truss is provided to the maximum likelihood estimator. Closed-form solutions of the partial differential equations and closed-form expressions of the sensitivity functions are derived so that the estimation algorithm is highly efficient. The resulting estimates from test data by using Timoshenko beam model are found to be comparable to those derived from finite element analysis.

Description: A fiber pushthrough process was computationally simulated using three-dimensional finite element method. The interface material is replaced by an anisotropic material with greatly reduced shear modulus in order to simulate the fiber pushthrough process using a linear analysis. Such a procedure is easily implemented and is computational very effective. It can be used to predict fiber pushthrough load for a composite system at any temperature. The average interface shear strength obtained from pushthrough load can easily be separated into its two components: one that comes from frictioal stresses and the other that comes from chemical adhesion between fiber and the matrix and mechanical interlocking that develops due to shrinkage of the composite because of phase change during the processing. Step-by-step procedures are described to perform the computational simulation, to establish bounds on interfacial bond strength and to interpret interfacial bond quality.

Description: A method is described for calculating shape sensitivities, within MSC/NASTRAN, in a simple manner without resort to external programs. The method uses natural design variables to define the shape changes in a given structure. Once the shape sensitivities are obtained, the shape optimization process is carried out in a manner similar to property optimization processes. The capability of this method is illustrated by two examples: the shape optimization of a cantilever beam with holes, loaded by a point load at the free end (with the shape of the holes and the thickness of the beam selected as the design variables), and the shape optimization of a connecting rod subjected to several different loading and boundary conditions.

Description: The equilibrium equations and the compatibility conditions are fundamental to the analyses of structures. However, anyone who undertakes even a cursory generic study of the compatibility conditions can discover, with little effort, that historically this facet of structural mechanics had not been adequately researched by the profession. Now the compatibility conditions (CC's) have been researched and are understood to a great extent. For finite element discretizations, the CC's are banded and can be divided into three distinct categories: (1) the interface CC's; (2) the cluster or field CC's; and (3) the external CC's. The generation of CC's requires the separating of a local region, then writing the deformation displacement relation (ddr) for the region, and finally, the eliminating of the displacements from the ddr. The procedure to generate all three types of CC's is presented and illustrated through examples of finite element models. The uniqueness of the CC's thus generated is shown.

Description: A six-degree-of-freedom nonlinear semi-definite model with time varying mesh stiffness has been developed for the dynamic analysis of spur gears. The model includes a spur gear pair, two shafts, two inertias representing load and prime mover, and bearings. As the shaft and bearing dynamics have also been considered in the model, the effect of lateral-torsional vibration coupling on the dynamics of gears can be studied. In the nonlinear model developed several factors such as time varying mesh stiffness and damping, separation of teeth, backlash, single- and double-sided impacts, various gear errors and profile modifications have been considered. The dynamic response to internal excitation has been calculated by using the 'static transmission error method' developed. The software prepared (DYTEM) employs the digital simulation technique for the solution, and is capable of calculating dynamic tooth and mesh forces, dynamic factors for pinion and gear, dynamic transmission error, dynamic bearing forces and torsions of shafts. Numerical examples are given in order to demonstrate the effect of shaft and bearing dynamics on gear dynamics.

Description: Active structural control experiments conducted on a 24-ft pinned-free beam derived feedback compensators on the basis of a traveling-wave approach. A compensator is thus obtained which eliminates resonant behavior by absorbing all impinging power. A causal solution is derived for this noncausal compensator which mimics its behavior in a given frequency range, using the Wiener-Hopf. This optimal Wiener-Hopf compensator's structure-damping performance is found to exceed any obtainable by means of rate feedback. Performance limitations encompassed the discovery of frequencies above which the sensor and actuator were no longer dual and an inadvertent coupling of the control hardware to unmodeled structure torsion modes.

Description: The problem of impact is of considerable interest in laminated composite materials. Although important contributions have been made in the understanding the impact problem through numerical solutions, an analytical solution has not been available for the problem of impact of laminated plates. The present work gives an analytical solution to this problem, based on the usual Fourier series expansion for simply-supported plates, combined with Laplace transform techniques for the impact problem solution.

Description: A micromechanics analysis is used to study the effects of constituent properties on thermally induced stresses in continuous fiber reinforced composites. A finite element formulation is described, and results are presented for unidirectional carbon/epoxy laminates. It is shown that significant stresses develop in composites exposed to thermal excursions typical of spacecraft operating environments and that the fiber thermoelastic properties have a minimal effect on the magnitude of these stresses. The finite element micromechanics analysis is then extended to the study of multidirectional laminates using a simple global/local formulation. Damage initiation predictions are compared with experimental data, and factors controlling the initiation of damage are identified. Ways of improving the durability of composites are discussed.

Description: Finite-element solution technology developed for use in conjunction with advanced viscoplastic models is described. The development of such solution technology is necessary for performing stress/life analyses of engineering structural problems where the complex geometries and loadings make the conventional analytical solutions difficult. The versatility of the solution technology is demonstrated by applying it to viscoplastic models possessing different mathematical structures and encompassing isotropic and anisotropic materials. The computational results qualitatively replicate deformation behavior observed in experiments on prototypical structural components.

Description: A curved beam type of test specimen is evaluated for use in determining the through-the-thickness strength of laminated composites. Two variations of a curved beam specimen configuration (semicircular and elliptical) were tested to failure using static and fatigue loads. The static failure load for the semicircular specimens was found to be highly sensitive to flaw content, with the specimens falling into two distinct groups. This result supports the use of proof testing for structural validation. Static design allowables are derived based on the Weibull distribution. Fatigue data indicates no measured increase in specimen compliance prior to final fracture. All static and fatigue failures at room temperature dry conditions occurred catastrophically. The elliptical specimens demonstrated unusually high failure strengths indicating the presence of phenomena requiring further study. Results are also included for specimens exposed to a wet environment showing a matrix strength degradation due to moisture content. Further testing is underway to evaluate a fatigue methodology for matrix dominated failures based on residual static strength (wearout).

Description: The behavior of elliptical sublaminates created by delaminations in composite plates was investigated for plates subjected to in-plane compressive, shear, and thermal loads. The axes of the ellipse may be arbitrarily oriented with respect to the applied loads. A model was developed which provides the stresses, strains, and displacements of the sublaminate, and the loads applied to the plate at which the sublaminate buckles and at which it grows in size. A series of experiments was conducted on sandwich plates made of Fiberite T300/976 graphite-epoxy laminates bonded to an aluminum honeycomb core. Teflon film in either a circle or an ellipse was embedded in the laminate, simulating the presence of a delamination. Each plate was loaded in compression, and the load-strain history of the sublaminate, the sublaminate buckling load, and the sublaminate growth load were measured. The predictions of the model agreed reasonably well with the data.

Description: Experimental PdCr temperature-compensated resistance static-strain gages are described. The gages are developed in both fine-wire and thin-film forms. It is found that a PdCr wire strain gage coated with a flame-sprayed mixture of alumina and 4 wt pct zirconia demonstrates the smallest variation in and the best repeatability of apparent strain among the existing gages used at temperatures up to 800 C. Results of preliminary tests indicate uncompensated uncoated thin-film gages have potential usefulness at temperatures up to 1000 C.

Description: A nonlinear elastic force-displacement relationship is used to calculate the transient impact force and local deformation at the point of contact between impactor and target. The nonlinear analysis and transfer function capabilities of NASTRAN are used to define a finite element model that behaves globally linearly elastic, and locally nonlinear elastic to model the local contact behavior. Results are presented for two different structures: a uniform cylindrical rod impacted longitudinally; and an orthotropic plate impacted transversely. Calculated impact force and transient structural response of the targets are shown to compare well with results measured in experimental tests.

Description: A special class of nonlinear damping models is investigated in which the damping force is proportional to the product of positive integer or the fractional power of the absolute values of displacement and velocity. For a one-degree-of-freedom system, the classical Krylov-Bogoliubov 'averaging' method is used, whereas for a distributed system, both an ad hoc perturbation technique and the finite difference method are employed to study the effects of nonlinear damping. The results are compared with linear viscous damping models. The amplitude decrement of free vibration for a single mode system with nonlinear models depends not only on the damping ratio but also on the initial amplitude, the time to measure the response, the frequency of the system, and the powers of displacement and velocity. For the distributed system, the action of nonlinear damping is found to reduce the energy of the system and to pass energy to lower modes.

Description: Continuum representations of micromechanical phenomena in structured materials are described, with emphasis on cellular solids. These phenomena are interpreted in light of Cosserat elasticity, a generalized continuum theory which admits degrees of freedom not present in classical elasticity. These are the rotation of points in the material, and a couple per unit area or couple stress. Experimental work in this area is reviewed, and other interpretation schemes are discussed. The applicability of Cosserat elasticity to cellular solids and fibrous composite materials is considered as is the application of related generalized continuum theories. New experimental results are presented for foam materials with negative Poisson's ratios.

Description: A detailed study is made of the effects of variations in lamination and material parameters of thin-walled composite frames on their vibrational characteristics. The structures considered are semicircular thin-walled frames with I and J sections. The flanges and webs of the frames are modeled by using two-dimensional shell and plate finite elements. A mixed formulation is used with the fundamental unknowns consisting of both the generalized displacements and stress resultants in the frame. The frequencies and modes predicted by the two-dimensional finite-element model are compared with those obtained from experiments, as well as with the predictions of a one-dimensional, thin-walled-beam, finite-element model. A detailed study is made of the sensitivity of the vibrational response to variations in the fiber orientation, material properties of the individual layers, and boundary conditions.

Description: A computational methodology is developed for the prediction of passive damping in composite structures. The method involves multiple levels of damping modeling by integrating micromechanics, laminate, and structural damping theories. The effects of temperature and moisture on structural damping are included. The simulation of damping in the structural level is accomplished with finite-element discretization. Applications are performed on graphite/epoxy composite beams, plates, and shells to illustrate the methodology. Additional parametric studies demonstrate the variation of structural modal damping with respect to ply angles, fiber volume ratio, and temperature.

Description: A novel formulation termed the integrated force method (IFM) has been developed in recent years for analyzing structures. In this method all the internal forces are taken as independent variables, and the system equilibrium equations (EEs) are integrated with the global compatibility conditions (CCs) to form the governing set of equations. In IFM the CCs are obtained from the strain formulation of St. Venant, and no choices of redundant load systems have to be made, in constrast to the standard force method (SFM). This property of IFM allows the generation of the governing equation to be automated straightforwardly, as it is in the popular stiffness method (SM). In this report IFM and SM are compared relative to the structure of their respective equations, their conditioning, required solution methods, overall computational requirements, and convergence properties as these factors influence the accuracy of the results. Overall, this new version of the force method produces more accurate results than the stiffness method for comparable computational cost.

Description: This paper presents an approximate method for analyzing the two-dimensional friction contact problem so as to compute the dynamic response of a structure constrained by friction interfaces. The friction force at the joint is formulated based on the Coulomb model. The single-term harmonic balance scheme, together with the receptance approach of decoupling the effect of the friction force on the structure from those of the external forces has been utilized to obtain the steady state response. The computational efficiency and accuracy of the method are demonstrated by comparing the results with long-term time solutions.

Description: A computer program called UWCODA is presented. UWCODA is intended to assist in the design, analysis and optimization of composite plates. UWCODA combines a state-of-the-art global optimization algorithm (Improving Hit and Run) with classical lamination theory. Optimization results are presented for simple loading conditions as well as for complex, biaxial load conditions. The computer code proved to be very effective in the design of composite plates.

Description: The objective was to examine the room temperature fatigue and nonlinear cumulative fatigue damage behavior of the cast nickel-based superalloy, MAR M-247. The fatigue test matrix consisted of single-level, fully reversed fatigue experiments. Two series of tests were performed: one of the two baseline fatigue LCF (Low-Cycle Fatigue) life levels was used in the first loading block, and the HCF (High-Cycle Fatigue) baseline loading level was used in the second block in each series. For each series, duplicate tests were performed at each applied LCF life fraction.

Description: The Probabilistic Fracture Mechanics (PFM) is a promising method for estimating the fatigue life and inspection cycles for mechanical and structural components. The Probability Finite Element Method (PFEM), which is based on second moment analysis, has proved to be a promising, practical approach to handle problems with uncertainties. As the PFEM provides a powerful computational tool to determine first and second moment of random parameters, the second moment reliability method can be easily combined with PFEM to obtain measures of the reliability of the structural system. The method is also being applied to fatigue crack growth. Uncertainties in the material properties of advanced materials such as polycrystalline alloys, ceramics, and composites are commonly observed from experimental tests. This is mainly attributed to intrinsic microcracks, which are randomly distributed as a result of the applied load and the residual stress.

Description: A brief description is provided of the fundamental aspects of a quantification process. Progress since the last structural durability conference in 1989 is summarized. The methodology to date and that to be developed during the life of the program is presented. The uncertain factors are presented. The approach is outlined that is required to achieve component and/or system certification in the shortest possible time for affordable reliability risk. Two new elements appear in a block diagram: (1) uncertainties in human factor, and (2) uncertainties in the computer code. Research to quantify the uncertainties in the human factor was initiated and is discussed.

Description: The Instrumentation and Control Technology Division of the Lewis Research Center has been developing an in-house capability to make one dimensional and two dimensional optical strain measurements on high temperature test specimens. The measurements are based on a two-beam speckle-shift technique. The development of composite materials for use in high temperature applications is generating interest in using the speckle-shift technique to measure strains on small diameter fibers and wires of various compositions. The results of preliminary speckle correlation tests on wire and fiber specimens are covered, and the advanced system currently under development is described.

Description: Viewgraphs on intelligent structures technology are presented. Topics covered include: embedding electronics; electrical and mechanical compatibility; integrated circuit chip packaged for embedding; embedding devices within composite structures; test of embedded circuit in G/E coupon; temperature/humidity/bias test; single-chip microcomputer control experiment; and structural shape determination.

Description: Viewgraphs on robust control for uncertain structures are presented. Topics covered include: robust linear quadratic regulator (RLQR) formulas; mismatched LQR design; RLQR design; interpretations of RLQR design; disturbance rejection; and performance comparisons: RLQR vs. mismatched LQR.

Description: Viewgraphs on the application of controlled structures technology (CST) to adaptive optics are presented. Topics covered include: a typical large optical system (LDR); overview of current optical programs; typical adaptive optics system; actuation approaches for deforming a mirror; control approach comparison; control structure interaction (CSI) control bandwidth limitations; applications of CST concepts to optics; distributed control approaches; quasistatic error correction; low authority control (LAC); and high authority control (HAC) design methodology.

Description: Viewgraphs on active impedance matching of complex structural systems are presented. Topics covered include: traveling wave model; dereverberated mobility model; computation of dereverberated mobility; control problem: optimal impedance matching; H2 optimal solution; statistical energy analysis (SEA) solution; experimental transfer functions; interferometer actuator and sensor locations; active strut configurations; power dual variables; dereverberation of complex structure; dereverberated transfer function; compensators; and relative power flow.

Description: Viewgraphs on inhibiting multiple mode vibration in controlled flexible systems are presented. Topics covered include: input pre-shaping background; developing multiple-mode shapers; Middeck Active Control Experiment (MACE) test article; and tests and results.

Description: Viewgraphs on the Middeck Active Control Experiment (MACE) are presented. Topics covered include: science program objectives and rationale; science requirements; capturing the essential physics; science development approach; development model hardware; development model test plan; and flight hardware and operations.

Description: The Plane Wave Approximation (PWA) is widely used in finite element analysis to implement the loading generated by an underwater shock wave. The method required to implement the PWA in NASTRAN is presented along with example problems. A theoretical background is provided and the limitations of the PWA are discussed.

Description: An identification procedure is proposed to identify damage characteristics (location and size of the damage) from dynamic measurements. This procedure was based on minimization of the mean-square measure of difference between measurement data (natural frequencies and mode shapes) and the corresponding predictions obtained from the computational model. The procedure is tested for simulated damage in the form of stiffness changes in a simple fixed free spring mass system and symmetric cracks in a simply supported Bernoulli Euler beam. It is shown that when all the mode information is used in the identification procedure it is possible to uniquely determine the damage properties. Without knowing the complete set of modal information, a restricted region in the initial data space has been found for realistic and convergent solution from the identification process.

Description: Residual stresses in finished parts have often been regarded as factors contributing to premature part failure and geometric distortions. Currently, residual stresses in welded structures and railroad components are being investigated. High residual stresses formed in welded structures due primarily to the differential contractions of the weld material as it cools and solidifies can have a profound effect on the surface performance of the structure. In railroad wheels, repeated use of the brakes causes high residual stresses in the rims which may lead to wheel failure and possible derailment. The goals of the study were: (1) to develop strategies for using x-ray diffraction to measure residual stress; (2) to subject samples of Inconel 718 to various mechanical and heat treatments and to measure the resulting stress using x-ray diffraction; and (3) to measure residual stresses in ferromagnetic alloys using magnetoacoustics.

Description: A new method of conducting lab vibration tests of spacecraft equipment was developed to more closely simulate the vibration environment experienced when the spacecraft is launched on a rocket. The improved tests are tailored to identify equipment design and workmanship problems without inducing artificial failures that would not have occurred at launch. These new, less destructive types of vibration tests are essential to JPL's protoflight test approach in which lab testing is conducted using the flight equipment, often one of a kind, to save time and money. In conventional vibration tests, only the input vibratory motion is specified; the feedback, or reaction force, between the test item and the vibration machine is ignored. Most test failures occur when the test item goes into resonance, and the reaction force becomes very large. It has long been recognized that the large reaction force is a test artifact which does not occur with the lightweight, flexible mounting structures characteristic of spacecraft and space vehicles. In new vibration tests, both the motion and the force provided to the test item by the vibration machine are controlled, so that the vibration ride experienced by the test item is as in flight.

Description: A finite element approach using NASTRAN is developed for solving time-dependent fluid-structure interaction problems, with emphasis on the transient scattering of acoustic waves from submerged elastic structures. Finite elements are used for modeling both structure and fluid domains to facilitate the graphical display of the wave motion through both media. For the liquid, the use of velocity potential as the fundamental unknown results in a symmetric matrix equation. The approach is illustrated for the problem of transient scattering from a submerged elastic spherical shell subjected to an incident tone burst. The use of an analogy between the equations of elasticity and the wave equation of acoustics, a necessary ingredient to the procedure, is summarized.

Description: The purpose here is to identify the structural parameters (plate thickness, liquid level, beam size, number of beams, tray diameter, etc.) that affect the structural integrity of distillation trays in distillation columns. Once the sensitivity of the trays' dynamic response to these parameters has been established, the designer will be able to use this information to prepare more accurate specifications for the construction of new trays. Information is given on both static and dynamic analysis, modal response, and tray failure details.

Description: The capabilities of the postprocessing program CANDI (Color Animation of Nastran DIsplacements) were expanded to accept results from axisymmetric analysis. An auxiliary program, ANIMATE, was developed to allow color display of CANDI output on the IRIS 4D-series workstations. The user can interactively manipulate the graphics display by three-dimensional rotations, translations, and scaling through the use of the keyboard and/or dials box. The user can also specify what portion of the model is displayed. These developments are limited to the display of complex displacements calculated with the NASHUA/NASTRAN procedure for structural acoustics analysis.

Description: Quite often the dynamicist will be faced with having an electric drive motor as a link in the elastic path of a structure such that the motor's characteristics must be taken into account to properly represent the dynamics of the primary structure. He does not want to model it so accurately that he could get detailed stress and displacements in the motor proper, but just sufficiently to represent its inertia loading and elastic behavior from its mounting bolts to its drive coupling. Described here is how the rotor and stator of such a motor can be adequately modeled as a colinear pair of beams.

Description: Discussed here are some recent finite element modeling experiences using the axisymmetric elements CONEAX, TRAPAX, and TRIAAX, from the COSMIC NASTRAN element library. These experiences were gained in the practical application of these elements to the static analysis of helicopter rotor force measuring systems for two design projects for the NASA Ames Research Center. These design projects were the Rotor Test Apparatus and the Large Rotor Test Apparatus, which are dedicated to basic helicopter research. Here, a genetic axisymmetric model is generated for illustrative purposes. Modeling considerations are discussed, and the advantages and disadvantages of using axisymmetric elements are presented. Asymmetric mechanical and thermal loads are applied to the structure, and single and multi-point constraints are addressed. An example that couples the axisymmetric model to a non-axisymmtric model is demonstrated, complete with DMAP alters. Recommendations for improving the elements and making them easier to use are offered.

Description: A method for computing the mass and center of gravity for basic and combined substructures stored in the NASTRAN Substructure Operating File (SOF) is described. The three step method recovers SOF data blocks for the relevant substructure, processes these data blocks using a specially developed FORTRAN routine, and generates the NASTRAN gridpoint weight generator (GPWG) table for the substructure in a PHASE2 SOF execution using a Direct Matrix Abstraction Program (DMAP) sequence. Verification data for the process is also provided.

Description: In the particular NASTRAN application discussed here, a nut was attached to a stationary structure. The object of the analysis was to determine the vibration characteristics of the whole structure for various configurations; i.e., the evaluation of the mode shapes and frequencies when parts were moved to different mating positions. Details of the analysis are given.

Description: It is demonstrated here that complete checkout of a basic substructure can be done under the special circumstance of a sliding connection with offsets. Stiff bar connectors make this possible so long as the bar coordinates are aligned with the displacement coordinates at the sliding surface.

Description: The 1990 UNIX version of NASTRAN was ported to two new platforms that are not supported by COSMIC: the Sun SPARC workstation and the Apple Macintosh using the A/UX version of UNIX. The experiences of the authers in porting NASTRAN is summarized here. Suggestions for users who might attempt similar ports are given.

Description: The distribution of NASTRAN User Manual information has been difficult because of the delay in printing and difficulty in identification of all the users. This has caused many users not to have the current information for the release of NASTRAN that could be available to them. The User Manual updates have been supplied with the NASTRAN Releases, but distribution within organizations was not coordinated with access to releases. The Executive Control, Case Control, and Bulk Data sections are supplied in machine readable format with the 91 Release of NASTRAN. This information is supplied on the release tapes in ASCII format, and a FORTRAN program to access this information is supplied on the release tapes. This will allow each user to have immediate access to User Manual level documentation with the release. The sections on utilities, plotting, and substructures are expected to be prepared for the 92 Release.

Description: The new 1991 COSMIC/NASTRAN version, compatible with the older versions, tries to remove some old constraints and make it easier to extract information from the plot file. It also includes some useful improvements and new enhancements. New features available in the 1991 version are described. They include a new PLT1 tape with simplified ASCII plot commands and short records, combined hidden and shrunk plot, an x-y-z coordinate system on all structural plots, element offset plot, improved character size control, improved FIND and NOFIND logic, a new NASPLOT post-prosessor to perform screen plotting or generate PostScript files, and a BASIC/NASTPLOT program for PC.

Description: Future proposed NASA missions will require large precision truss type structures that are deployed or assembled in space. To date, space structures that are important to missions success have been ground tested to validate their performance. Evaluation of the performance requirements of future systems has shown that using current and projected design and test approaches, the structure cannot be adequately validated by ground test. A problem exists since it is believed that unless important structure systems can be validated by ground tests, they will never be adopted for future missions. New design or ground test approaches are necessary to enable future missions. The inability of current approaches to validate future structural systems is discussed.

Description: The present FEM reduction method for the generation of test-analysis-models (TAMs) in test-analysis correlation addresses contentions that the current modal TAM is hypersensitive to differences between test model shapes and analysis mode shapes, thereby generating large off-diagonal terms within the orthogonality and cross-orthogonality matrices employed in test-analysis mode-shape correlation. A hybrid TAM methodology is accordingly developed which combines the exact representation of the FEM target modes from the modal TAM with the more accurate TAM representation of the residual modes; the superior residual dynamics representation of the hybrid TAM is demonstrated for a detailed representation of a large space structure.

Description: Results are presented of elastic, elastic-plastic and elastic-plastic-creep analyses of a test-rig component of an actively cooled cowl lip. A cowl lip is part of the leading edge of an engine inlet of proposed hypersonic aircraft and is subject to severe thermal loadings and gradients during flight. Values of stresses calculated by elastic analysis are well above the yield strength of the cowl lip material. Such values are highly unrealistic, and thus elastic stress analyses are inappropriate. The inelastic (elastic-plastic and elastic-plastic-creep) analyses produce more reasonable and acceptable stress and strain distributions in the component. Finally, using the results from these analyses, predictions are made for the cyclic crack initiation life of a cowl lip. A comparison of predicted cyclic lives shows the cyclic life prediction from the elastic-plastic-creep analysis to be the lowest and, hence, most realistic.

Description: A 24-ply composite double cantilever-beam specimen under mode I (opening) loading has been analyzed by a 3D FEM code that calculated along a straight delamination starter for several different specimen materials. An isotropic specimen was found to have a strain-energy release rate distribution which varied along its delamination front due to the boundary-layer effect and another effect associated with the anticlastic curvature of the bent specimen arms. A 0-deg graphite-reinforced epoxy specimen had a nearly-uniform strain-energy release rate distribution which dropped only near the edge, due to the boundary-layer effect, and a +/- 45-deg graphite/epoxy specimen exhibited a pronounced strain-energy release rate variation across the specimen width.

Description: A method is developed for predicting the S-N curve of a composite laminate which is subjected to an arbitrary stress ratio, R (minimum stress/maximum stress). The method is based on the measuring of the S-N behavior of two distinct cases, tension-tension and compression-compression fatigue loadings. Using these parameters, expressions are formulated that estimate the fatigue behavior under any stress ratio loading. Experimental results from the testing of graphite/epoxy laminates, with various structures, are compared with the predictions and show good agreement.

Description: The energy-absorption capability as a function of crushing speed was determined for Thornel 300-Fiberite 934 (Gr-E) and Kevlar-49-Fiberite 934 (K-E) composite material. Circular cross section tube specimens were crushed at speeds ranging from 0.01 to 12 m/sec. Ply orientations of the tube specimens were (0/ +/- theta)2 and (+/- theta)2 where theta = 15, 45, and 75 deg. Based upon the results of these tests, the energy-absorption capability of Gr-E and K-E was determined to be a function of crushing speed. The magnitude of the effects of crushing speed on energy-absorption capability was determined to be a function of the mechanisms that control the crushing process. The effects of crushing speed on the energy-absorption capability is related to whether the mechanical response of the crushing mechanism that controls the crushing process is a function of strain rate. Energy-absorption capability of Gr-E and K-E tubes ranged between 0 and 35 percent and 20 and 45 percent, respectively, depending upon ply orientation.

Description: A Rayleigh-Ritz analysis is presented to determine the load or strain at which delamination buckling will occur for a composite laminate containing a single delamination. This delamination is assumed to be elliptical in shape, with local axes of symmetry which may be at an angle relative to the global plate axes. The laminate is assumed to be sujected to global loading which transforms to a loading with a relatively small shear component when taken about the local axes of the elliptical delamination. Special cases of this type of loading include uniaxial and biaxial compression. Possible bending-stretching coupling behavior of the delaminated region and the effect of Poisson's ratio mismatch between the delaminated and base region are investigated. Analytical predictions for delamination buckling loads are shown to correlate well with experimental results for a number of different geometries.

Description: The fatigue behavior of orthotropic laminates (0, +/-45, 0 deg)2s and (90, +/-45, 90 deg)2s, has been evaluated under alternating tension-compression loading. Even though the first laminate is much stronger than the second, both started to fail by delamination. Visual damage started to show only at the very end of the fatigue life but measurement of the stiffness showed that degradation starts at about 80 percent of the fatigue life. The first laminate failed in compression after delamination between the 0 and the 45 deg laminae, while the second failed in tension after delamination between the +45 and -45 deg laminae. It is shown that the interlaminar fatigue strength of both laminate structures can be correlated to the applied fatigue load.

Description: Two methods for calculating linear frequency domain aerodynamic coefficients from a time marching Full Potential cascade solver are developed and verified. In the first method, the Influence Coefficient, solutions to elemental problems are superposed to obtain the solutions for a cascade in which all blades are vibrating with a constant interblade phase angle. The elemental problem consists of a single blade in the cascade oscillating while the other blades remain stationary. In the second method, the Pulse Response, the response to the transient motion of a blade is used to calculate influence coefficients. This is done by calculating the Fourier Transforms of the blade motion and the response. Both methods are validated by comparison with the Harmonic Oscillation method and give accurate results. The aerodynamic coefficients obtained from these methods are used for frequency domain flutter calculations involving a typical section blade structural model. An eigenvalue problem is solved for each interblade phase angle mode and the eigenvalues are used to determine aeroelastic stability. Flutter calculations are performed for two examples over a range of subsonic Mach numbers.

Description: Information on probabilistic structural analysis methods for space propulsion systems is given in viewgraph form. Information is given on deterministic certification methods, probability of failure, component response analysis, stress responses for 2nd stage turbine blades, Space Shuttle Main Engine (SSME) structural durability, and program plans. .

Description: The use of probabilistic structural analysis methods (PSAM) to predict structural reliability is the subject of an on-going NASA research program. The elements of the new technology developed to date is reported. Applications of the developed software to structural problems are demonstrated for simple validation problems and for large scale application problems. On-going research to support component and system reliability predictions suitable for analytical certification of aerospace structures is briefly reviewed.

Description: The purpose is to develop procedures to extract sensitivity analysis information from COSMIC/NASTRAN and to couple it with a mathematical optimization package. At present, the analysis will be limited to stress, displacement, and frequency constraints with structures modeled with membrane elements, rods, and bar elements. Two types of sensitivity analysis are discussed: an adjoint variable approach which is most effective when the number of active constraints is significantly less than the number of physical variables, and an approach based on a first order approximation of a Taylor series. The latter approach is more effective when the number of independent design variables is significantly less than the number of active constraints.

Description: Post-processing algorithms are given to compute the vibratory elastic-rigid coupling matrices and the modal contributions to the rigid-body mass matrices and to the effective modal inertias and masses. Recomputation of the elastic-rigid coupling matrices for a change in origin is also described. A computational example is included. The algorithms can all be executed by using standard finite-element program eigenvalue analysis output with no changes to existing code or source programs.