ALBERT

Description: This thesis presents the on-going development of methodology for a probabilistic material strength degradation model. The probabilistic model, in the form of a postulated randomized multifactor equation, provides for quantification of uncertainty in the lifetime material strength of aerospace propulsion system components subjected to a number of diverse random effects. This model is embodied in the computer program entitled PROMISS, which can include up to eighteen different effects. Presently, the model includes four effects that typically reduce lifetime strength: high temperature, mechanical fatigue, creep, and thermal fatigue. Statistical analysis was conducted on experimental Inconel 718 data obtained from the open literature. This analysis provided regression parameters for use as the model's empirical material constants, thus calibrating the model specifically for Inconel 718. Model calibration was carried out for four variables, namely, high temperature, mechanical fatigue, creep, and thermal fatigue. Methodology to estimate standard deviations of these material constants for input into the probabilistic material strength model was developed. Using the current version of PROMISS, entitled PROMISS93, a sensitivity study for the combined effects of mechanical fatigue, creep, and thermal fatigue was performed. Results, in the form of cumulative distribution functions, illustrated the sensitivity of lifetime strength to any current value of an effect. In addition, verification studies comparing a combination of mechanical fatigue and high temperature effects by model to the combination by experiment were conducted. Thus, for Inconel 718, the basic model assumption of independence between effects was evaluated. Results from this limited verification study strongly supported this assumption.

Description: An anisotropic viscoplasticity model is developed as an extension of the well known Bodner model. The extension is made by replacing the effective stress of the isotropic Bodner model by one involving invariants for transverse isotropy. The resulting model retains the simplicity of Bodner's in the ease with which the material constants are determined experimentally. It allows a representation of strong initial anisotropy yet is based on the scalar state variable under the assertion that induced anisotropy is negligible relative to the strong initial anisotropy. Temperature dependence is taken as in the original Bodner theory. Account is made of fiber volume fraction through nonlinear rules of mixture applied to the stress history and anisotropy parameters. Focus is on the theoretical development of the model, however, application to a W/Cu composite is in progress and will be reported as a sequel to this report.

Description: This work formulates a method for the modeling of material damping characteristics in distributed parameter models which may be easily applied to models such as rod, plate, and beam equations. The general linear boundary value vibration equation is modified to incorporate hysteresis effects represented by complex stiffness using the transfer function approach proposed by Golla and Hughes. The governing characteristic equations are decoupled through separation of variables yielding solutions similar to those of undamped classical theory, allowing solution of the steady state as well as transient response. Example problems and solutions are provided demonstrating the similarity of the solutions to those of the classical theories and transient responses of nonviscous systems.

Description: We consider the modeling and active damping of an elastic beam using distributed actuators and sensors. The piezoelectric ceramic material (PZT) is used to build the actuator. The sensor is made of the piezoelectric polymer polyvinylidene fluoride (PVDF). These materials are glued on both sides of the beam. For the simple clamped beam, the closed loop controller has been shown to be able to extract energy from the beam. The shape of the actuator and its influence on the closed loop system performance are discussed. It is shown that it is possible to suppress the selected mode by choosing the appropriate actuator layout. It is also shown that by properly installing the sensor and determining the sensor shape we can further extract and manipulate the sensor signal for our control need.

Description: The equations of motion for a two-segment deploying telescopic beam are derived through application of Lagrange's equation. The outer tube of the beam is fixed at one end and the inner tube slides freely relative to the fixed segment. The resulting nonlinear, non-autonomous set of equations is linearized and simplified to the standard Euler-Bernoulli partial differential equations for an elastic beam by freezing the deployment process at various stages of deployment, and examining the small amplitude and natural modes of vibration of the resulting configuration. Application of the natural boundary conditions and compatibility of motion relations for the two segments in their common region of overlap leads to a transcendental characteristic equation in the frequency parameter Beta(L). Numerical solution of the equation for the characteristic roots determines the modal frequencies, and the corresponding mode shapes are obtained from the general solution of the Euler-Bernoulli equation tailored to the natural boundary conditions. Sample results of modal frequencies and shapes are presented for various stages of deployment and discussed. It is shown that for all intermediate stages of deployment (between 0 and 100 percent) the spectral distribution is drastically altered by the appearance of regions of very closely spaced modal frequencies. The sources of this modal agglomeration are explored.

Description: An experimental investigation of the post-impact fatigue response of integrally woven carbon/epoxy composites was conducted. Five different through-the-thickness (TTT) reinforcing fibers were used in an experimental textile process that produced an integrally woven (0/90/0/90/0/90/0/90/0)(sub T) ply layup with 21K AS4 carbon tow fiber. The resin was Hercules 3501-6, and the five TTT reinforcing fibers were Kevlar, Toray carbon, AS4 carbon, glass, and IM6 carbon. The purpose of this investigation was to study the post-impact fatigue response of these material systems and to identify the optimum TTT fiber. Samples were impacted with one half inch diameter aluminum balls with an average velocity of 543 ft/sec. Post-impact static compression and constant amplitude tension-compression fatigue tests were conducted. Fatigue tests were conducted with a loading ratio of R=-5, and frequency of 4 Hz. Damage growth was monitored using x-radiographic and sectioning techniques and by examining the stress-strain response (across the impact site) throughout the fatigue tests. The static compressive stress versus far-field strain response was nearly linear for all material groups. All the samples had a transverse shear failure mode. The average compressive modulus (from far-field strain) was about 10 Msi. The average post-impact static compressive strength was about 35.5 Ksi. The IM6 carbon sample had a strength of over 40 Ksi, more than 16 percent stronger than average. There was considerable scatter in the S-N data. However, the IM6 carbon samples clearly had the best fatigue response. The response of the other materials, while worse than IM6 carbon, could not be ranked definitively. The initial damage zones caused by the impact loading and damage growth from fatigue loading were similar for all five TTT reinforcing materials. The initial damage zones were circular and consisted of delaminations, matrix cracks and ply cracks. Post-impact fatigue loading caused delamination growth, ply cracking and fiber bundle failures, typically 45 deg from impact load direction. During the initial 97 percent of fatigue life, delaminations, ply cracks and fiber bundle failures primarily grew at and near the impact site. During the final 3 percent of life, damage grew rapidly transverse to the loading direction as a through-the-thickness transverse shear failure. The stress-strain response was typically linear during the initial 50 percent of life, and stiffness dropped about 20 percent during this period. During the next 47 percent of life, stiffness dropped about 34 percent, and the stress-strain response was no longer linear. The stiffness decreased about 23 percent during the final 3 percent of life. These trends were typical of all the materials tested. Therefore, by monitoring stiffness loss, fatigue failure could be accurately anticipated.

Description: Minimizing the weight of structural components of the space station launched into orbit in a space shuttle can save cost, reduce the number of space shuttle missions, and facilitate on-orbit fabrication. Traditional manual design of such components, although feasible, cannot represent a minimum weight condition. At NASA Lewis Research Center, a design capability called CometBoards (Comparative Evaluation Test Bed of Optimization and Analysis Routines for the Design of Structures) has been developed especially for the design optimization of such flight components. Two components of the space station - a spacer structure and a support system - illustrate the capability of CometBoards. These components are designed for loads and behavior constraints that arise from a variety of flight accelerations and maneuvers. The optimization process using CometBoards reduced the weights of the components by one third from those obtained with traditional manual design. This paper presents a brief overview of the design code CometBoards and a description of the space station components, their design environments, behavior limitations, and attributes of their optimum designs.

Description: To assist analysis of riveted lap joints, stress intensity factors are determined for surface and corner cracks emanating from a wedge-loaded hole by using a 3-D weight function method in conjunction with a 3-D finite element method. A stress intensity factor equation for surface cracks is also developed to provide a closed-form solution. The equation covers commonly-encountered geometrical ranges and retains high accuracy over the entire range.

Description: The Technology Benefit Estimator (T/BEST) system is a formal method to assess advanced technologies and quantify the benefit contributions for prioritization. T/BEST may be used to provide guidelines to identify and prioritize high payoff research areas, help manage research and limited resources, show the link between advanced concepts and the bottom line, i.e., accrued benefit and value, and to communicate credibly the benefits of research. The T/BEST software computer program is specifically designed to estimating benefits, and benefit sensitivities, of introducing new technologies into existing propulsion systems. Key engine cycle, structural, fluid, mission and cost analysis modules are used to provide a framework for interfacing with advanced technologies. An open-ended, modular approach is used to allow for modification and addition of both key and advanced technology modules. T/BEST has a hierarchical framework that yields varying levels of benefit estimation accuracy that are dependent on the degree of input detail available. This hierarchical feature permits rapid estimation of technology benefits even when the technology is at the conceptual stage. As knowledge of the technology details increases the accuracy of the benefit analysis increases. Included in T/BEST's framework are correlations developed from a statistical data base that is relied upon if there is insufficient information given in a particular area, e.g., fuel capacity or aircraft landing weight. Statistical predictions are not required if these data are specified in the mission requirements. The engine cycle, structural fluid, cost, noise, and emissions analyses interact with the default or user material and component libraries to yield estimates of specific global benefits: range, speed, thrust, capacity, component life, noise, emissions, specific fuel consumption, component and engine weights, pre-certification test, mission performance engine cost, direct operating cost, life cycle cost, manufacturing cost, development cost, risk, and development time. Currently, T/BEST operates on stand-alone or networked workstations, and uses a UNIX shell or script to control the operation of interfaced FORTRAN based analyses. T/BEST's interface structure works equally well with non-FORTRAN or mixed software analysis. This interface structure is designed to maintain the integrity of the expert's analyses by interfacing with expert's existing input and output files. Parameter input and output data (e.g., number of blades, hub diameters, etc.) are passed via T/BEST's neutral file, while copious data (e.g., finite element models, profiles, etc.) are passed via file pointers that point to the expert's analyses output files. In order to make the communications between the T/BEST's neutral file and attached analyses codes simple, only two software commands, PUT and GET, are required. This simplicity permits easy access to all input and output variables contained within the neutral file. Both public domain and proprietary analyses codes may be attached with a minimal amount of effort, while maintaining full data and analysis integrity, and security. T/BESt's sotware framework, status, beginner-to-expert operation, interface architecture, analysis module addition, and key analysis modules are discussed. Representative examples of T/BEST benefit analyses are shown.

Description: The variability in initial crack size and fatigue crack growth is incorporated in a probabilistic model that is used to predict the fatigue lives for unstiffened aluminum alloy panels containing multiple site damage (MSD). The uncertainty of the damage in the MSD panel is represented by a distribution of fatigue crack lengths that are analytically derived from equivalent initial flaw sizes. The variability in fatigue crack growth rate is characterized by stochastic descriptions of crack growth parameters for a modified Paris crack growth law. A Monte-Carlo simulation explicitly describes the MSD panel by randomly selecting values from the stochastic variables and then grows the MSD cracks with a deterministic fatigue model until the panel fails. Different simulations investigate the influences of the fatigue variability on the distributions of remaining fatigue lives. Six cases that consider fixed and variable conditions of initial crack size and fatigue crack growth rate are examined. The crack size distribution exhibited a dominant effect on the remaining fatigue life distribution, and the variable crack growth rate exhibited a lesser effect on the distribution. In addition, the probabilistic model predicted that only a small percentage of the life remains after a lead crack develops in the MSD panel.

Description: Recent experimental fatigue crack growth studies have concluded an apparent anomalous behavior of short cracks. To investigate the reasons for this unexpected behavior, the present paper focuses on identifying the crack length circumstances under which the requirements for a single parameter (K(sub I) or delta K(sub I) if cyclic loading is considered) characterization are violated. Furthermore, an additional quantity, the T stress, as introduced by Rice, and the related biaxiality ratio, B, are calculated for several crack lengths and two configurations, the single-edge-cracked and the centrally-cracked specimen. It is postulated that a two-parameter characterization by K and T (or B) is needed for the adequate description of the stress and strain field around a short crack. To further verify the validity of this postulate, the influence of the third term of the Williams series on the stress, strain and displacement fields around the crack tip and in particular on the B parameter is also examined. It is found that the biaxiality ratio would be more negative if the third term effects are included in both geometries. The study is conducted using the finite element method with linearly elastic material and isoparametric elements and axial (mode I) loading. Moreover, it is clearly shown that it is not proper to postulate the crack size limits for 'short crack' behavior as a normalized ratio with the specimen width, a/w; it should instead be stated as an absolute, or normalized with respect to a small characteristic dimension such as the grain size. Finally, implications regarding the prediction of cyclic (fatigue) growth of short cracks are discussed.

Description: A methodology for simulating the growth of long through cracks in the skin of pressurized aircraft fuselage structures is described. Crack trajectories are allowed to be arbitrary and are computed as part of the simulation. The interaction between the mechanical loads acting on the superstructure and the local structural response near the crack tips is accounted for by employing a hierarchical modeling strategy. The structural response for each cracked configuration is obtained using a geometrically nonlinear shell finite element analysis procedure. Four stress intensity factors, two for membrane behavior and two for bending using Kirchhoff plate theory, are computed using an extension of the modified crack closure integral method. Crack trajectories are determined by applying the maximum tangential stress criterion. Crack growth results in localized mesh deletion, and the deletion regions are remeshed automatically using a newly developed all-quadrilateral meshing algorithm. The effectiveness of the methodology and its applicability to performing practical analyses of realistic structures is demonstrated by simulating curvilinear crack growth in a fuselage panel that is representative of a typical narrow-body aircraft. The predicted crack trajectory and fatigue life compare well with measurements of these same quantities from a full-scale pressurized panel test.

Description: The crack arrest capability of a tear strap in a pressurized precracked fuselage was studied through instrumented axial rupture tests of small scale models of an idealized fuselage. Upon pressurization, rapid crack propagation initiated at an axial through crack along the stringer and immediately kinked due to the mixed modes 1 and 2 state caused by the one-sided opening of the crack flap. The diagonally running crack further turned at the tear straps. Dynamic finite element analysis of the rupturing cylinder showed that the crack kinked and also ran straight in the presence of a mixed mode state according to a modified two-parameter crack kinking criterion.

Description: Fatigue crack growth associated with cyclic (secondary) plastic flow near a crack front is modelled using an incremental formulation. A new description of threshold behaviour under small load cycles is included. Quasi-static crack extension under high load excursions is described using an incremental formulation of the R-(crack growth resistance)- curve concept. The integration of the equations is discussed. For constant amplitude load cycles the results will be compared with existing crack growth laws. It will be shown that the model also properly describes interaction effects of fatigue crack growth and quasi-static crack extension. To evaluate the more general applicability the model is included in the NASGRO computer code for damage tolerance analysis. For this purpose the NASGRO program was provided with the CORPUS and the STRIP-YIELD models for computation of the crack opening load levels. The implementation is discussed and recent results of the verification are presented.

Description: The bending effects of unsymmetrically bonded composite repairs on cracked aluminum panels were quantified using a plate linear finite element model. Stress intensity factors and strain energy release rates were obtained from the model twice, once with out-of-plane displacement suppressed and another time without these restrictions. Several configurations were examined, crack growth stability was identified, and the effect of a debond was considered. The maximum stress intensity factor was also analyzed. Previous work by other authors was found to underpredict the bending effect.

Description: This paper presents recent developments in advanced analysis methods for the computation of stress site damage. The method of solution is based on the p-version of the finite element method. Its implementation was designed to permit extraction of linear stress intensity factors using a superconvergent extraction method (known as the contour integral method) and evaluation of the J-integral following an elastic-plastic analysis. Coarse meshes are adequate for obtaining accurate results supported by p-convergence data. The elastic-plastic analysis is based on the deformation theory of plasticity and the von Mises yield criterion. The model problem consists of an aluminum plate with six equally spaced holes and a crack emanating from each hole. The cracks are of different sizes. The panel is subjected to a remote tensile load. Experimental results are available for the panel. The plasticity analysis provided the same limit load as the experimentally determined load. The results of elastic-plastic analysis were compared with the results of linear elastic analysis in an effort to evaluate how plastic zone sizes influence the crack growth rates. The onset of net-section yielding was determined also. The results show that crack growth rate is accelerated by the presence of adjacent damage, and the critical crack size is shorter when the effects of plasticity are taken into consideration. This work also addresses the effects of alternative stress-strain laws: The elastic-ideally-plastic material model is compared against the Ramberg-Osgood model.

Description: Computational efficiency in structural optimization can be enhanced if the intensive computations associated with the calculation of the sensitivities, that is, gradients of the behavior constraints, are reduced. Approximation to gradients of the behavior constraints that can be generated with small amount of numerical calculations is proposed. Structural optimization with these approximate sensitivities produced correct optimum solution. Approximate gradients performed well for different nonlinear programming methods, such as the sequence of unconstrained minimization technique, method of feasible directions, sequence of quadratic programming, and sequence of linear programming. Structural optimization with approximate gradients can reduce by one third the CPU time that would otherwise be required to solve the problem with explicit closed-form gradients. The proposed gradient approximation shows potential to reduce intensive computation that has been associated with traditional structural optimization.

Description: Finite element analysis is regularly used during the engineering cycle of mechanical systems to predict the response to static, thermal, and dynamic loads. The finite element model (FEM) used to represent the system is often correlated with physical test results to determine the validity of analytical results provided. Results from dynamic testing provide one means for performing this correlation. One of the most common methods of measuring accuracy is by classical modal testing, whereby vibratory mode shapes are compared to mode shapes provided by finite element analysis. The degree of correlation between the test and analytical mode shapes can be shown mathematically using the cross orthogonality check. A great deal of time and effort can be exhausted in generating the set of test acquired mode shapes needed for the cross orthogonality check. In most situations response data from vibration tests are digitally processed to generate the mode shapes from a combination of modal parameters, forcing functions, and recorded response data. An alternate method is proposed in which the same correlation of analytical and test acquired mode shapes can be achieved without conducting the modal survey. Instead a procedure is detailed in which a minimum of test information, specifically the acceleration response data from a random vibration test, is used to generate a set of equivalent local accelerations to be applied to the reduced analytical model at discrete points corresponding to the test measurement locations. The static solution of the analytical model then produces a set of deformations that once normalized can be used to represent the test acquired mode shapes in the cross orthogonality relation. The method proposed has been shown to provide accurate results for both a simple analytical model as well as a complex space flight structure.

Description: This paper develops a parallelizable multilevel multiple constrained nonlinear equation solver. The substructuring process is automated to yield appropriately balanced partitioning of each succeeding level. Due to the generality of the procedure,_sequential, as well as partially and fully parallel environments can be handled. This includes both single and multiprocessor assignment per individual partition. Several benchmark examples are presented. These illustrate the robustness of the procedure as well as its capability to yield significant reductions in memory utilization and calculational effort due both to updating and inversion.

Description: In the course of previous composite structures test programs, the need for and the feasibility of developing analyses for scale-up effects has been demonstrated. The analysis techniques for scale-up effects fall into two categories. The first category pertains to developing analysis methods independently for a single, unique failure mode in composites, and using this compendium of analysis methods together with a global structural model to identify and predict the response and failure mode of full-scale built-up structures. The second category of scale-up effects pertains to similitude in structural validation testing. In this latter category, dimensional analysis is used to develop scale-up laws that enable extrapolation of sub-scale component test data to full-scale structures. This viewgraph presentation describes the approach taken and some developments accomplished in the first category of analysis for scale-up effects. Layup dependence of composite material properties severely limits the use of the dimensional analysis approach and these limitations are illustrated by examples.

Description: A model which can be used to predict the two-dimensional nonlinear behavior of bridged cracks in orthotropic strips is presented. The results obtained using a singular integral equation formulation which incorporates the anisotropy rigorously show that, although the effects of anisotropy are significant, the nondimensional quantities employed by Cox and Marshall can generate nearly universal results (R-curves, for example) for different levels of relative anisotropy. The role of composite constituent properties in the behavior of bridged cracks is clarified.

Description: This report presents the results of Phase 1 research to develop a methodology for performing large-scale Multi-disciplinary Stochastic Optimization (MSO) for the design of aerospace systems ranging from aeropropulsion components to complete aircraft configurations. The current research recognizes that such design optimization problems are computationally expensive, and require the use of either massively parallel or multiple-processor computers. The methodology also recognizes that many operational and performance parameters are uncertain, and that uncertainty must be considered explicitly to achieve optimum performance and cost. The objective of this Phase 1 research was to initialize the development of an MSO methodology that is portable to a wide variety of hardware platforms, while achieving efficient, large-scale parallelism when multiple processors are available. The first effort in the project was a literature review of available computer hardware, as well as review of portable, parallel programming environments. The first effort was to implement the MSO methodology for a problem using the portable parallel programming language, Parallel Virtual Machine (PVM). The third and final effort was to demonstrate the example on a variety of computers, including a distributed-memory multiprocessor, a distributed-memory network of workstations, and a single-processor workstation. Results indicate the MSO methodology can be well-applied towards large-scale aerospace design problems. Nearly perfect linear speedup was demonstrated for computation of optimization sensitivity coefficients on both a 128-node distributed-memory multiprocessor (the Intel iPSC/860) and a network of workstations (speedups of almost 19 times achieved for 20 workstations). Very high parallel efficiencies (75 percent for 31 processors and 60 percent for 50 processors) were also achieved for computation of aerodynamic influence coefficients on the Intel. Finally, the multi-level parallelization strategy that will be needed for large-scale MSO problems was demonstrated to be highly efficient. The same parallel code instructions were used on both platforms, demonstrating portability. There are many applications for which MSO can be applied, including NASA's High-Speed-Civil Transport, and advanced propulsion systems. The use of MSO will reduce design and development time and testing costs dramatically.

Description: A major goal in the control of complex mechanical system such as spacecraft rocket engine's advanced aircraft, and power plants is to achieve high performance with increased reliability, component durability, and maintainability. The current practice of decision and control systems synthesis focuses on improving performance and diagnostic capabilities under constraints that often do not adequately represent the materials degradation. In view of the high performance requirements of the system and availability of improved materials, the lack of appropriate knowledge about the properties of these materials will lead to either less than achievable performance due to overly conservative design, or over-straining of the structure leading to unexpected failures and drastic reduction of the service life. The key idea in this report is that a significant improvement in service life could be achieved by a small reduction in the system dynamic performance. The major task is to characterize the damage generation process, and then utilize this information in a mathematical form to synthesize a control law that would meet the system requirements and simultaneously satisfy the constraints that are imposed by the material and structural properties of the critical components. The concept of damage mitigation is introduced for control of mechanical systems to achieve high performance with a prolonged life span. A model of fatigue damage dynamics is formulated in the continuous-time setting, instead of a cycle-based representation, for direct application to control systems synthesis. An optimal control policy is then formulated via nonlinear programming under specified constraints of the damage rate and accumulated damage. The results of simulation experiments for the transient upthrust of a bipropellant rocket engine are presented to demonstrate efficacy of the damage-mitigating control concept.

Description: This paper will briefly set forth some of the basic tenets of mistuned rotating bladed-disk assemblies. The experience with an existing three stage compressor in a transonic wind tunnel will be documented. The manner in which the theoretical properties manifest themselves in this non-ideal compressor will be described. A description of mistuning behaviors that can and cannot be accurately substantiated will be discussed.

Description: A buckling analysis was performed on a hat-stiffened panel subjected to shear loading. Both local buckling and global buckling were analyzed. The global shear buckling load was found to be several times higher than the local shear buckling load. The classical shear buckling theory for a flat plate was found to be useful in predicting the local shear buckling load of the hat-stiffened panel, and the predicted local shear buckling loads thus obtained compare favorably with the results of finite element analysis.

Description: The formulation and application of element-level, element-independent error indicators is investigated. This research culminates in the development of an error indicator formulation which is derived based on the projection of element deformation onto the intrinsic element displacement modes. The qualifier 'element-level' means that no information from adjacent elements is used for error estimation. This property is ideally suited for obtaining error values and driving adaptive mesh refinements on parallel computers where access to neighboring elements residing on different processors may incur significant overhead. In addition such estimators are insensitive to the presence of physical interfaces and junctures. An error indicator qualifies as 'element-independent' when only visible quantities such as element stiffness and nodal displacements are used to quantify error. Error evaluation at the element level and element independence for the error indicator are highly desired properties for computing error in production-level finite element codes. Four element-level error indicators have been constructed. Two of the indicators are based on variational formulation of the element stiffness and are element-dependent. Their derivations are retained for developmental purposes. The second two indicators mimic and exceed the first two in performance but require no special formulation of the element stiffness mesh refinement which we demonstrate for two dimensional plane stress problems. The parallelizing of substructures and adaptive mesh refinement is discussed and the final error indicator using two-dimensional plane-stress and three-dimensional shell problems is demonstrated.

Description: Thermocryogenic buckling and stress analyses were conducted on a horizontally oriented cryogenic tank using the finite element method. The tank is a finite-length circular cylindrical shell with its two ends capped with hemispherical shells. The tank is subjected to cylindrical strip heating in the region above the liquid-cryogen fill level and to cryogenic cooling below the fill level (i.e., under thermocryogenic loading). The effects of cryogen fill level on the buckling temperature and thermocryogenic stress field were investigated in detail. Both the buckling temperature and stress magnitudes were relatively insensitive to the cryogen fill level. The buckling temperature, however, was quite sensitive to the radius-to-thickness ratio. A mechanical stress analysis of the tank also was conducted when the tank was under: (1) cryogen liquid pressure loading; (2) internal pressure loading; and (3) tank-wall inertia loading. Deformed shapes of the cryogenic tanks under different loading conditions were shown, and high-stress domains were mapped on the tank wall for the strain-gage installations. The accuracies of solutions from different finite element models were compared.

Description: Issues which impact on the interpretation and quantification of proof test benefits are reviewed. The importance of each issue in contributing to the extra quality assurance conferred by proof testing components is discussed, particularly with respect to the application of advanced fracture mechanics concepts to enhance the flaw screening capability of a proof test analysis. Items covered include the role in proof testing of elastic-plastic fracture mechanics, ductile instability analysis, deterministic versus probabilistic analysis, single versus multiple cycle proof testing, and non-destructive examination (NDE). The effects of proof testing on subsequent service life are reviewed, particularly with regard to stress redistribution and changes in fracture behavior resulting from the overload. The importance of proof test conditions are also addressed, covering aspects related to test temperature, simulation of service environments, test media and the application of real-time NDE. The role of each issue in a proof test methodology is assessed with respect to its ability to: promote proof test practice to a state-of-the-art; aid optimization of proof test design; and increase awareness and understanding of outstanding issues.

Description: A parallel adaptive dynamic relaxation (ADR) algorithm has been developed for nonlinear structural analysis. This algorithm has minimal memory requirements, is easily parallelizable and scalable to many processors, and is generally very reliable and efficient for highly nonlinear problems. Performance evaluations on single-processor computers have shown that the ADR algorithm is reliable and highly vectorizable, and that it is competitive with direct solution methods for the highly nonlinear problems considered. The present algorithm is implemented on the 512-processor Intel Touchstone DELTA system at Caltech, and it is designed to minimize the extent and frequency of interprocessor communication. The algorithm has been used to solve for the nonlinear static response of two and three dimensional hyperelastic systems involving contact. Impressive relative speedups have been achieved and demonstrate the high scalability of the ADR algorithm. For the class of problems addressed, the ADR algorithm represents a very promising approach for parallel-vector processing.

Description: The FPCAS2D computer code has been developed for aeroelastic stability analysis of bladed disks such as those in fans, compressors, turbines, propellers, or propfans. The aerodynamic analysis used in this code is based on the unsteady two-dimensional full potential equation which is solved for a cascade of blades. The structural analysis is based on a two degree-of-freedom rigid typical section model for each blade. Detailed explanations of the aerodynamic analysis, the numerical algorithms, and the aeroelastic analysis are not given in this report. This guide can be used to assist in the preparation of the input data required by the FPCAS2D code. A complete description of the input data is provided in this report. In addition, four test cases, including inputs and outputs, are provided.

Description: Large-displacement elastic and elastic-plastic, finite-element stress-strain analyses of an oxygen-tree high-conductivity (OFHC) copper plate specimen were performed using an updated Lagrangian formulation. The plate specimen is intended for low-cost experiments that emulate the most important thermomechanical loading and failure modes of a more complex rocket nozzle. The plate, which is loaded in bending at 593 C, contains a centrally located and internally pressurized channel. The cyclic crack initiation lives were estimated using the results from the analyses and isothermal strain-controlled low-cycle fatigue data for OFHC copper. A comparison of the predicted and experimental cyclic lives showed that an elastic analysis predicts a longer cyclic life than that observed in experiments by a factor greater than 4. The results from elastic-plastic analysis for the plate bend specimen, however, predicted a cyclic life in close agreement with experiment, thus justifying the need for the more rigorous stress-strain analysis.

Description: The formal approach described herein computationally simulates the probable ranges of uncertainties for the human factor in probabilistic assessments of structural reliability. Human factors such as marital status, professional status, home life, job satisfaction, work load, and health are studied by using a multifactor interaction equation (MFIE) model to demonstrate the approach. Parametric studies in conjunction with judgment are used to select reasonable values for the participating factors (primitive variables). Subsequently performed probabilistic sensitivity studies assess the suitability of the MFIE as well as the validity of the whole approach. Results show that uncertainties range from 5 to 30 percent for the most optimistic case, assuming 100 percent for no error (perfect performance).

Description: A mixed method of approximation based on Reissner's variational principle is developed for the linear analysis of interlaminar stresses in laminated composites, with special interest in laminates that contain terminated internal plies (dropped-ply laminates). Two models are derived, one for problems of generalized plane deformation and the other for the axisymmetric response of shells of revolution. A layerwise approach is taken in which the stress field is assumed with an explicit dependence on the thickness coordinate in each layer. The dependence of the stress field on the thickness coordinate is determined such that the three-dimensional equilibrium equations are satisfied by the approximation. The solution domain is reduced to one dimension by integration through the thickness. Continuity of tractions and displacements between layers is imposed. The governing two-point boundary value problem is composed of a system of both differential and algebraic equations (DAE's) and their associated boundary conditions. Careful evaluation of the system of DAE's was required to arrive at a form that allowed application of a one-step finite difference approximation. A two-stage Gauss implicit Runge-Kutta finite difference scheme was used for the solution because of its relatively high degree of accuracy. Patch tests of the two models revealed problems with solution accuracy for the axisymmetric model of a cylindrical shell loaded by internal pressure. Parametric studies of dropped-ply laminate characteristics and their influence on the interlaminar stresses were performed using the generalized plane deformation model. Eccentricity of the middle surface of the laminate through the ply drop-off was found to have a minimal effect on the interlaminar stresses under longitudinal compression, transverse tension, and in-plane shear. A second study found the stiffness change across the ply termination to have a much greater influence on the interlaminar stresses.

Description: The failure sequence following crack formation in a chevron-notched four-point bend 1 specimen is examined in a parametric study using the Bluhm slice synthesis model. Premature failure resulting from crack formation forces which exceed those required to propagate a crack beyond alpha (min) is examined together with the critical crack length and critical crack front length. An energy based approach is used to establish factors which forecast the tendency of such premature failure due to crack formation for any selected chevron-notched geometry. A comparative study reveals that, for constant values of alpha (1) and alpha (0), the dimensionless beam compliance and stress intensity factor are essentially independent of specimen width and thickness. The chevron tip position, alpha (0) has its primary effect on the force required to initiate a sharp crack. Small values for alpha (0) maximize the stable region length, however, the premature failure tendency is also high for smaller alpha (0) values. Improvements in premature failure resistance can be realized for larger values of alpha (0) with only a minor reduction in the stable region length. The stable region length is also maximized for larger chevron based positions, alpha (1) but the chance for premature failure is also raised. Smaller base positions improve the premature failure resistance with only minor decreases in the stable region length. Chevron geometries having a good balance of premature failure resistance, stable region length, and crack front length are 0.20 less than or equal to alpha (0) is less than or equal to 0.30 and 0.70 is less than or equal to alpha (1) is less than or equal to 0.80.

Description: The accuracy and robustness of a two-dimensional smoothing methodology is examined for the problem of recovering accurate interlaminar shear stress distributions in laminated composite and sandwich plates. The smoothing methodology is based on a variational formulation which combines discrete least-squares and penalty-constraint functionals in a single variational form. The smoothing analysis utilizes optimal strains computed at discrete locations in a finite element analysis. These discrete strain data are smoothed with a smoothing element discretization, producing superior accuracy strains and their first gradients. The approach enables the resulting smooth strain field to be practically C1-continuous throughout the domain of smoothing, exhibiting superconvergent properties of the smoothed quantity. The continuous strain gradients are also obtained directly from the solution. The recovered strain gradients are subsequently employed in the integration o equilibrium equations to obtain accurate interlaminar shear stresses. The problem is a simply-supported rectangular plate under a doubly sinusoidal load. The problem has an exact analytic solution which serves as a measure of goodness of the recovered interlaminar shear stresses. The method has the versatility of being applicable to the analysis of rather general and complex structures built of distinct components and materials, such as found in aircraft design. For these types of structures, the smoothing is achieved with 'patches', each patch covering the domain in which the smoothed quantity is physically continuous.

Description: This paper describes the mechanical design, analysis, fabrication, testing, and lessons learned by developing a uniquely designed spaceflight-like actuator. The linear proof mass actuator (LPMA) was designed to attach to both a large space structure and a ground test model without modification. Previous designs lacked the power to perform in a terrestrial environment while other designs failed to produce the desired accelerations or frequency range for spaceflight applications. Thus, the design for a unique actuator was conceived and developed at NASA Langley Research Center. The basic design consists of four large mechanical parts (mass, upper housing, lower housing, and center support) and numerous smaller supporting components including an accelerometer, encoder, and four drive motors. Fabrication personnel were included early in the design phase of the LPMA as part of an integrated manufacturing process to alleviate potential difficulties in machining an already challenging design. Operating testing of the LPMA demonstrated that the actuator is capable of various types of load functions.

Description: A component specific modeling software program has been developed for propulsion systems. This expert program is capable of formulating the component geometry as finite element meshes for structural analysis which, in the future, can be spun off as NURB geometry for manufacturing. COSMO currently has geometry recipes for combustors, turbine blades, vanes, and disks. Component geometry recipes for nozzles, inlets, frames, shafts, and ducts are being added. COSMO uses component recipes that work through neutral files with the Technology Benefit Estimator (T/BEST) program which provides the necessary base parameters and loadings. This report contains the users manual for combustors, turbine blades, vanes, and disks.

Description: The object of the research described in this paper was to develop a means of identifying the wave-types propagating between two points in a finite length beam. It is known that different structural wave-types possess different dispersion relations: i.e., that their group speeds and the frequency dependence of their group speeds differ. As a result of those distinct dispersion relationships, different wave-types may be associated with characteristic features when structural responses are examined in the time frequency domain. Previously, the time-frequency character of analytically generated structural responses of both single element and multi-element structures were examined by using the Wigner Distribution (WD) along with filtering techniques that were designed to detect the wave-types present in the responses. In the work to be described here, the measure time-frequency response of finite length beam is examined using the WD and filtering procedures. This paper is organized as follows. First the concept of time-frequency analysis of structural responses is explained. The WD is then introduced along with a description of the implementation of a discrete version. The time-frequency filtering techniques are then presented and explained. The results of applying the WD and the filtering techniques to the analysis of a transient response is then presented.

Description: Response of a structure to a convected pressure field generated by far-field jet-noise or turbulent boundary layer is often of interest in engineering. The structure considered in present study is a periodic beam with equally rotational elastic supports. The convected pressure field is characterized by a broad frequency spectrum and a particular convection velocity. In the first order approximation, all the frequency components in the pressure spectrum are convected with the same velocity. In this paper, both finitely and infinitely long periodic beams under convected harmonic pressure are investigated. The exact harmonic responses are obtained by introducing new expressions of mode shapes for these two types of beams into the normal mode approach. Thus, both locations and values of the response peaks in the propagation bands can be precisely predicted and the important coincidence phenomenon can be investigated in exact terms. This work provides a foundation on which the study on vibration control of structure can be conducted.

Description: In this paper, the beam under general time dependent stationary random excitation is investigated, when exact solution is unavailable. Numerical simulations are carried out to compare its results with those yielded by the conventional linearization techniques. It is found that the modified version of the stochastic linearization technique yields considerably more accurate results for the mean square displacement of the beam than the conventional equivalent linearization technique, especially in the case of large nonlinearity.

Description: Continuous fiber composites have shown tremendous promise in industrial applications. Their microstructures, however, are very complex and in many instances difficult to characterize. In this project, the fracture characteristics of a specially made fiber reinforced composite with different fiber spacing are investigated. The experimental results so far have shown that after an initial transient phase the crack speed reaches a steady phase, i.e., independent of the crack length. Within the steady crack growth phase debonding along the fibers in the bridging zone grows in a self-similar manner. During the steady phase the energy dissipation per cycle is constant. Afterwards, an increase of the energy dissipation is observed that is accompanied by a decrease in crack speed. This latter trend is presumed to be the result of relatively large amounts of energy dissipated in the bulk of the specimen. Using appropriate Green's function and computer simulations, the stress intensity factor at the crack tip is evaluated for various cases of bridging stresses. In this way the effects of specimen size and fiber spacing on the overall fracture behavior of the composite system are analyzed. The steady crack speed and the steady rate of debonding have a similar power dependence on stress level. Dimensional analysis demonstrates that the particular fracture process is not governed by dimensional invariance but on the detailed micromechanisms in the bridging zone.

Description: This paper develops a simplified continuum (continuous wrp to time, stress, etc.) fatigue damage model for use in Life Extending Controls (LEC) studies. The work is based on zero mean stress local strain cyclic damage modeling. New nonlinear explicit equation forms of cyclic damage in terms of stress amplitude are derived to facilitate the continuum modeling. Stress based continuum models are derived. Extension to plastic strain-strain rate models are also presented. Application of these models to LEC applications is considered. Progress toward a nonzero mean stress based continuum model is presented. Also, new nonlinear explicit equation forms in terms of stress amplitude are also derived for this case.

Description: Utilizing a proposed fixture design or some variation thereof, this report presents a verification approach to strength test space flight payload components, electronics boxes, mechanisms, lines, fittings, etc., which traditionally do not lend themselves to classical static loading. The fixture, through use of ordered Euler rotation angles derived herein, can be mounted on existing vibration shakers and can provide an innovative method of applying single axis flight load vectors. The versatile fixture effectively loads protoflight or prototype components in all three axes simultaneously by use of a sinusoidal burst of desired magnitude at less than one-third the first resonant frequency. Cost savings along with improved hardware confidence are shown. The end product is an efficient way to verify experiment hardware for both random vibration and strength.

Description: The Thermal/Structural Tailoring of Engine Blades (T/STAEBL) system is a family of computer programs executed by a control program. The T/STAEBL system performs design optimizations of cooled, hollow turbine blades and vanes. This manual describes the T/STAEBL data block structure and system organization. The approximate analysis and optimization modules are detailed, and a validation test case is provided.

Description: Generic load models are described with multiple levels of progressive sophistication to simulate the composite (combined) load spectra (CLS) that are induced in space propulsion system components, representative of Space Shuttle Main Engines (SSME), such as transfer ducts, turbine blades and liquid oxygen (LOX) posts. These generic (coupled) models combine the deterministic models for composite load dynamic, acoustic, high-pressure and high rotational speed, etc., load simulation using statistically varying coefficients. These coefficients are then determined using advanced probabilistic simulation methods with and without strategically selected experimental data. The entire simulation process is included in a CLS computer code. Applications of the computer code to various components in conjunction with the PSAM (Probabilistic Structural Analysis Method) to perform probabilistic load evaluation and life prediction evaluations are also described to illustrate the effectiveness of the coupled model approach.

Description: This report describes an analytical study of vibration reduction in a four-bladed helicopter rotor using an actively controlled, partial span, trailing edge flap located on the blade. The vibration reduction produced by the actively controlled flap (ACF) is compared with that obtained using individual blade control (IBC), in which the entire blade is oscillated in pitch. For both cases a deterministic feedback controller is implemented to reduce the 4/rev hub loads. For all cases considered, the ACF produced vibration reduction comparable with that obtained using IBC, but consumed only 10-30% of the power required to implement IBC. A careful parametric study is conducted to determine the influence of blade torsional stiffness, spanwise location of the control flap, and hinge moment correction on the vibration reduction characteristics of the ACF. The results clearly demonstrate the feasibility of this new approach to vibration reduction. It should be emphasized than the ACF, used together with a conventional swashplate, is completely decoupled from the primary flight control system and thus it has no influence on the airworthiness of the helicopter. This attribute is potentially a significant advantage when compared to IBC.

Description: Applications are described of high-performance parallel, computation for the analysis of complete jet engines, considering its multi-discipline coupled problem. The coupled problem involves interaction of structures with gas dynamics, heat conduction and heat transfer in aircraft engines. The methodology issues addressed include: consistent discrete formulation of coupled problems with emphasis on coupling phenomena; effect of partitioning strategies, augmentation and temporal solution procedures; sensitivity of response to problem parameters; and methods for interfacing multiscale discretizations in different single fields. The computer implementation issues addressed include: parallel treatment of coupled systems; domain decomposition and mesh partitioning strategies; data representation in object-oriented form and mapping to hardware driven representation, and tradeoff studies between partitioning schemes and fully coupled treatment.

Description: The theory and user instructions for an optimization code based on the method of feasible directions are presented. The code was written for wide distribution and ease of attachment to other simulation software. Although the theory of the method of feasible direction was developed in the 1960's, many considerations are involved in its actual implementation as a computer code. Included in the code are a number of features to improve robustness in optimization. The search direction is obtained by solving a quadratic program using an interior method based on Karmarkar's algorithm. The theory is discussed focusing on the important and often overlooked role played by the various parameters guiding the iterations within the program. Also discussed is a robust approach for handling infeasible starting points. The code was validated by solving a variety of structural optimization test problems that have known solutions obtained by other optimization codes. It has been observed that this code is robust: it has solved a variety of problems from different starting points. However, the code is inefficient in that it takes considerable CPU time as compared with certain other available codes. Further work is required to improve its efficiency while retaining its robustness.

Description: The Thermal/Structural Tailoring of Engine Blades (T/STAEBL) system is a computer code that is able to perform numerical optimizations of cooled jet engine turbine blades and vanes. These optimizations seek an airfoil design of minimum operating cost that satisfies realistic design constraints. This report documents the organization of the T/STAEBL computer program, its design and analysis procedure, its optimization procedure, and provides an overview of the input required to run the program, as well as the computer resources required for its effective use. Additionally, usage of the program is demonstrated through a validation test case.

Description: The structural design and configuration feasibility of the long spacer truss assembly that will be used as part of the Space Station Freedom is the focus of this study. The structural analysis discussed herein is derived from the transient loading events presented in the Space Transportation System Interface Control Document (STS ICD). The transient loading events are liftoff, landing, and emergency landing loads. Quasi-static loading events were neglected in this study since the magnitude of the quasi-static acceleration factors is lower than that of the transient acceleration factors. Structural analysis of the proposed configuration of the long spacer truss with four longerons indicated that negative safety margins are possible. As a result, configuration changes were proposed. The primary configuration change suggested was to increase the number of truss longerons to six. The six-longeron truss appears to be a more promising structure than the four-longeron truss because it offers a positive margin of safety and more volume in its second bay (BAY2). This additional volume can be used for resupply of some of the orbital replacement units (such as a battery box). Note that the design effort on the long spacer truss has not fully begun and that calculations and reports of the negative safety margins are, to date, based on concept only.

Description: The inelastic hysteresis energy applied to the material in a cycle is used as the basis for predicting nonisothermal fatigue life of a wrought cobalt-base superalloy, Haynes 188, from isothermal fatigue data. Damage functions that account for hold-time effects and time-dependent environmental phenomena such as oxidation and hot corrosion are proposed in terms of the inelastic hysteresis energy per cycle. The proposed damage functions are used to predict the bithermal and thermomechanical fatigue lives of Haynes 188 between 316 and 760 C from isothermal fatigue data. Predicted fatigue lives of all but two of the nonisothermal tests are within a factor of 1.5 of the experimentally observed lives.

Description: A shake test was conducted in the 80- by 120-Foot Wind Tunnel at NASA Ames Research Center, using the NASA Ames Rotor Test Apparatus (RTA) and the Sikorsky S-76 rotor hub. The primary objective of this shake test was to determine the modal properties of the RTA, the S-76 rotor hub, and the model support system installed in the wind tunnel. Random excitation was applied at the rotor hub, and vibration responses were measured using accelerometers mounted at various critical locations on the model and the model support system. Transfer functions were computed using the load cell data and the accelerometer responses. The transfer function data were used to compute the system modal parameters with the aid of modal analysis software.

Description: In the aircraft industry, full scale and large component testing is a very necessary, time consuming, and expensive process. It is essential to find ways by which this process can be minimized without loss of reliability. One possible alternative is the use of scaled down models in testing and use of the model test results in order to predict the behavior of the larger system, referred to herein as prototype. This viewgraph presentation provides justifications and motivation for the research study, and it describes the necessary conditions (similarity conditions) for two structural systems to be structurally similar with similar behavioral response. Similarity conditions provide the relationship between a scaled down model and its prototype. Thus, scaled down models can be used to predict the behavior of the prototype by extrapolating their experimental data. Since satisfying all similarity conditions simultaneously is in most cases impractical, distorted models with partial similarity can be employed. Establishment of similarity conditions, based on the direct use of the governing equations, is discussed and their use in the design of models is presented. Examples include the use of models for the analysis of cylindrical bending of orthotropic laminated beam plates, of buckling of symmetric laminated rectangular plates subjected to uniform uniaxial compression and shear, applied individually, and of vibrational response of the same rectangular plates. Extensions and future tasks are also described.

Description: Structural analysis of piping systems, especially dynamic analysis, typically considers the duct structure and the contained fluid column separately. Coupling of these two systems, however, forms a new dynamic system with characteristics not necessarily described by the superposition of the two component system's characteristics. Methods for modeling the two coupled components simultaneously using finite line elements are presented. Techniques for general duct intersections, area or direction changes, long radius bends, hydraulic losses, and hydraulic impedances are discussed. An example problem and results involving time transients are presented. Additionally, a program to enhance post-processing of line element models is discussed.

Description: This report summarizes work that has been performed to date on the strength of a cryotank insulation system using Rohacell foam and TUFI-coated AETB-12 ceramic tiles directly bonded to a simulated graphite-epoxy tank wall. Testing utilized a custom specimen design which consists of a long tensile specimen with eccentric loading to induce curvature similar to the curvature expected due to 'pillowing' of the tank when pressurized. A finite element model was constructed to predict the specific element strains in the test article, and to assist with design of the test specimen to meet the specific goals of curvature and laminate strain. The results indicate that the heat treated 3.25-pcf density Rohacell foam does not provide sufficient strength for the induced stresses due to curvature and stress concentration at the RTV bondline to the TUFI tile. The test was repeated using higher density non-heat treated Rohacell foam (6.9 pcf) without foam failure. The finite element model was shown to predict specimen behavior, and validation of the model was successful. It is pertinent to mention that the analyses described herein accurately predicted the failure of the heat treated foams and based on this analysis method it is expected that the untreated 3.25 pcf Rohacell foam will be successful.

Description: Butt-weld joints are most commonly designed into pressure vessels which then become as reliable as the weakest increment in the weld chain. In practice, weld material properties are determined from tensile test specimen and provided to the stress analyst in the form of a stress versus strain diagram. Variations in properties through the thickness of the weld and along the width of the weld have been suspect but not explored because of inaccessibility and cost. The purpose of this study is to investigate analytical and computational methods used for analysis of welds. The weld specimens are analyzed using classical elastic and plastic theory to provide a basis for modeling the inelastic properties in a finite-element solution. The results of the analysis are compared to experimental data to determine the weld behavior and the accuracy of prediction methods. The weld considered in this study is a multiple-pass aluminum 2219-T87 butt weld with thickness of 1.40 in. The weld specimen is modeled using the finite-element code ABAQUS. The finite-element model is used to produce the stress-strain behavior in the elastic and plastic regimes and to determine Poisson's ratio in the plastic region. The value of Poisson's ratio in the plastic regime is then compared to experimental data. The results of the comparisons are used to explain multipass weld behavior and to make recommendations concerning the analysis and testing of welds.

Description: A force method formulation, the completed Beltrami-Michell formulation (CBMF), has been developed for analyzing boundary value problems in elastic continua. The CBMF is obtained by augmenting the classical Beltrami-Michell formulation with novel boundary compatibility conditions. It can analyze general elastic continua with stress, displacement, or mixed boundary conditions. The CBMF alleviates the limitations of the classical formulation, which can solve stress boundary value problems only. In this report, the CBMF is specialized for plates and shells. All equations of the CBMF, including the boundary compatibility conditions, are derived from the variational formulation of the integrated force method (IFM). These equations are defined only in terms of stresses. Their solution for kinematically stable elastic continua provides stress fields without any reference to displacements. In addition, a stress function formulation for plates and shells is developed by augmenting the classical Airy's formulation with boundary compatibility conditions expressed in terms of the stress function. The versatility of the CBMF and the augmented stress function formulation is demonstrated through analytical solutions of several mixed boundary value problems. The example problems include a composite circular plate and a composite circular cylindrical shell under the simultaneous actions of mechanical and thermal loads.

Description: Personnel in the Quality and Integrity Design Engineering Center (QIDEC) at the University of Utah are working under a two year grant from the FAA to better understand the role of fretting corrosion and fretting fatigue in aircraft rivet hole cracking. The current program follows a one year grant program which was completed in 1993. This paper provides a status report on the results of these grant programs. Recent effort has been focused on developing basic fretting fatigue models which consider variation in the coefficient of friction with time and location within the fretting interface. This is a very important characteristic of the QIDEC model because coefficient of friction varies significantly during the fretting fatigue process. Copies of QIDEC documents discussed in this paper can be obtained by contacting the authors.

Description: Crack propagation tests were conducted to determine crack growth behavior in five helicopter materials for surface cracks between 0.005 to 0.020 inches in depth. Constant amplitude tests were conducted at stress ratios R equals 0.1 and 0.5, and emphasis was placed on near threshold data (i.e., 10-8 to 10-6 inches/cycle). Spectrum tests were conducted using a helicopter spectrum. The test specimen was an unnotched tension specimen, and cracks were initiated from a small EDM notch. An optical/video system was used to monitor crack growth. The material for the test specimens was obtained from helicopter part forgings. Testing was conducted at stresses below yield to reflect actual stresses in helicopter parts.

Description: The focus of this work is the buckling response of symmetrically laminated composite plates having a planform area in the shape of an isosceles trapezoid. The loading is assumed to be inplane and applied perpendicular to the parallel ends of the plate. The tapered edges of the plate are assumed to have simply supported boundary conditions, while the parallel ends are assumed to have either simply supported or clamped boundary conditions. A semi-analytic closed-form solution based on energy principles and the Trefftz stability criterion is derived and solutions are obtained using the Rayleigh-Ritz method. Intrinsic in this solution is a simplified prebuckling analysis which approximates the inplane force resultant distributions by the forms Nx=P/W(x) and Ny=Nxy=0, where P is the applied load and W(x) is the plate width which, for the trapezoidal planform, varies linearly with the lengthwise coordinate x. The out-of-plane displacement is approximated by a double trigonometric series. This analysis is posed in terms of four nondimensional parameters representing orthotropic and anisotropic material properties, and two nondimensional parameters representing geometric properties. For comparison purposes, a number of specific plate geometry, ply orientation, and stacking sequence combinations are investigated using the general purpose finite element code ABAQUS. Comparison of buckling coefficients calculated using the semi-analytical model and the finite element model show agreement within 5 percent, in general, and within 15 percent for the worst cases. In order to verify both the finite element and semi-analytical analyses, buckling loads are measured for graphite/epoxy plates having a wide range of plate geometries and stacking sequences. Test fixtures, instrumentation system, and experimental technique are described. Experimental results for the buckling load, the buckled mode shape, and the prebuckling plate stiffness are presented and show good agreement with the analytical results regarding the buckling load and the prebuckling plate stiffness. However, the experimental results show that for some cases the analysis underpredicts the number of halfwaves in the buckled mode shape. In the context of the definitions of taper ratio and aspect ratio used in this study, it is concluded that the buckling load always increases as taper ratio increases for a given aspect ratio for plates having simply supported boundary conditions on the parallel ends. There are combinations of plate geometry and ply stackling sequences, however, that reverse this trend for plates having clamped boundary conditions on the parallel ends such that an increase in the taper ratio causes a decrease in the buckling load. The clamped boundary conditions on the parallel ends of the plate are shown to increase the buckling load compared to simply supported boundary conditions. Also, anisotropy (the D16 and D26 terms) is shown to decrease the buckling load and skew the buckled mode shape for both the simply supported and clamped boundary conditions.

Description: In order to predict the dynamic response of a flexible structure in a fluid flow, the equations of motion of the structure and the fluid must be solved simultaneously. In this paper, we present several partitioned procedures for time-integrating this focus coupled problem and discuss their merits in terms of accuracy, stability, heterogeneous computing, I/O transfers, subcycling, and parallel processing. All theoretical results are derived for a one-dimensional piston model problem with a compressible flow, because the complete three-dimensional aeroelastic problem is difficult to analyze mathematically. However, the insight gained from the analysis of the coupled piston problem and the conclusions drawn from its numerical investigation are confirmed with the numerical simulation of the two-dimensional transient aeroelastic response of a flexible panel in a transonic nonlinear Euler flow regime.

Description: This report supplements a previous report of the same title submitted in June, 1992. It summarizes additional analytical techniques which have been developed for predicting the response of linear and nonlinear structures to noise excitations generated by large propulsion power plants. The report is divided into nine chapters. The first two deal with incomplete knowledge of boundary conditions of engineering structures. The incomplete knowledge is characterized by a convex set, and its diagnosis is formulated as a multi-hypothesis discrete decision-making algorithm with attendant criteria of adaptive termination.

Description: Parametric studies were performed to assess the effects of various parameters on the free-vibration behavior (natural frequencies) of (+/- theta)(sub 2) angle-ply, fiber composite, thin shell structures in a hygrothermal environment. Knowledge of the natural frequencies of structures is important in considering their response to various kinds of excitation, especially when structures and force systems are complex and when excitations are not periodic. The three dimensional, finite element structural analysis computer code CSTEM was used in the Cray YMP computer environment. The fiber composite shell was assumed to be cylindrical and made from T300 graphite fibers embedded in an intermediate-modulus, high-strength matrix. The following parameters were investigated: the length and the laminate thickness of the shell, the fiber orientation, the fiber volume fraction, the temperature profile through the thickness of the laminate, and laminates with different ply thicknesses. The results indicate that the fiber orientation and the length of the laminated shell had significant effects on the natural frequencies. The fiber volume fraction, the laminate thickness, and the temperature profile through the shell thickness had weak effects on the natural frequencies. Finally, the laminates with different ply thicknesses had an insignificant influence on the behavior of the vibrated laminated shell. Also, a single through-the-thickness, eight-node, three dimensional composite finite element analysis appears to be sufficient for investigating the free-vibration behavior of thin, composite, angle-ply shell structures.

Description: A user's guide for the computer program OPTCOMP is presented in this report. This program provides a capability to optimize the fabrication or service-induced residual stresses in uni-directional metal matrix composites subjected to combined thermo-mechanical axisymmetric loading using compensating or compliant layers at the fiber/matrix interface. The user specifies the architecture and the initial material parameters of the interfacial region, which can be either elastic or elastoplastic, and defines the design variables, together with the objective function, the associated constraints and the loading history through a user-friendly data input interface. The optimization procedure is based on an efficient solution methodology for the elastoplastic response of an arbitrarily layered multiple concentric cylinder model that is coupled to the commercial optimization package DOT. The solution methodology for the arbitrarily layered cylinder is based on the local-global stiffness matrix formulation and Mendelson's iterative technique of successive elastic solutions developed for elastoplastic boundary-value problems. The optimization algorithm employed in DOT is based on the method of feasible directions.

Description: Elastic-plastic finite-element analyses were used to calculate stresses and displacements around a crack in finite-thickness plates for an elastic-perfectly plastic material. Middle- and edge-crack specimens were analyzed under tension and bending loads. Specimens were 1.25 to 20 mm thick with various widths and crack lengths. A global constraint factor alpha(sub g), an averaged normal-stress to flow-stress ratio over the plastic region, was defined to simulate three-dimensional (3D) effects in two-dimensional (2D) models. For crack lengths and uncracked ligament lengths greater than four times the thickness, the global constraint factor was found to be nearly a unique function of a normalized stress-intensity factor (related to plastic-zone size to thickness ratio) from small- to large-scale yielding conditions for various specimen types and thickness. For crack length-to-thickness ratios less than four, the global constraint factor was specimen type, crack length and thickness dependent. Using a 2D strip-yield model and the global constraint factors, plastic-zone sizes and crack-tip displacements agreed reasonably well with the 3D analyses. For a thin sheet aluminum alloy, the critical crack-tip-opening angle during stable tearing was found to be independent of specimen type and crack length for crack length-to-thickness ratios greater than 4.

Description: The National Aeronautics and Space Administration and the Chinese Aeronautical Establishment participated in a Fatigue and Fracture Mechanics Cooperative Program. The program objectives were to identify and characterize crack initiation and growth of small cracks (10 microns to 2 mm long) in commonly used US and PRC aluminum alloys, to improve fracture mechanics analyses of surface- and corner-crack configurations, and to develop improved life-prediction methods. Fatigue and small-crack tests were performed on single-edgenotch tension (SENT) specimens and large-crack tests were conducted on center-crack tension specimens for constant-amplitude (stress ratios of -1, 0, and 0.5) and Mini-TWIST spectrum loading. The plastic replica method was used to monitor the initiation and growth of small fatigue cracks at the semicircular notch. Crack growth results from each laboratory on 7075-T6 bare and LC9cs clad aluminum alloys agreed well and showed that fatigue life was mostly crack propagation from a material defect (inclusion particles or void) or from the cladding layer. Finite-element and weight-function methods were used to determine stress intensity factors for surface and corner cracks in the SENT specimens. Equations were then developed and used in a crack growth and crack-closure model to correlate small- and large-crack data and to make life predictions for various load histories. The cooperative program produced useful experimental data and efficient analysis methods for improving life predictions. The results should ultimately improve aircraft structural reliability and safety.

Description: A new method of design decomposition for structural analysis and optimization is described. For this method, the structure is divided into substructures where each substructure has its structural response described by a structural-response subproblem, and its structural sizing determined from a structural-sizing subproblem. The structural responses of substructures that have rigid body modes when separated from the remainder of the structure are further decomposed into displacements that have no rigid body components, and a set of rigid body modes. The structural-response subproblems are linked together through forces determined within a structural-sizing coordination subproblem which also determines the magnitude of any rigid body displacements. Structural-sizing subproblems having constraints local to the substructures are linked together through penalty terms that are determined by a structural-sizing coordination subproblem. All the substructure structural-response subproblems are totally decoupled from each other, as are all the substructure structural-sizing subproblems, thus there is significant potential for use of parallel solution methods for these subproblems.

Description: A theoretical equation describing the elastic-plastic deformation of a cantilever beam subject to a constant pressure is developed. The theoretical result is compared numerically to the computer program BUCKY for the case of an elastic-perfectly plastic specimen. It is shown that the theoretical and numerical results compare favorably in the plastic range. Comparisons are made to another research code to further validate the BUCKY results. This paper serves as a quality test for the computer program BUCKY developed at NASA Johnson Space Center.

Description: Previous work using photostress to analyze behavior of aluminum welded joints was useful to determine mechanical properties of the weld and parent materials along the centerline of the joint. It was shown that significant differences exist in the stress-strain characteristics at points beginning at the centerline of the weld and extending for a distance of one inch to either side of the weld. Because of the highly variable behavior detected in the previous work, it was decided to extend the work to investigate behavior of joints 1/8, 1/2, and 1.40 inches thick.

Description: Contemporary applications of failure criteria frequently incorporate two-dimensional or simplified three-dimensional methodologies for prediction of stresses. Motivation behind the development of a new multi-dimensional failure criterion is due mainly to the lack of a sufficiently accurate mathematical tool that accounts for the behavior of brittle material with anisotropic properties. Such a criterion should be able to provide a reliable maximum load estimate so that design of the structure is not penalized in terms of excessive weight requirements. The failure criterion developed is represented by a fracture surface in a six-dimensional stress space. The criterion is applied for failure prediction of SR-200 beryllium sheet structures, a non-homogeneous orthotropic material used widely in space applications. Two experiments are used to verify the criterion.

Description: Extensive research activities were performed in the area of multidisciplinary modeling and simulation of aerospace vehicles that are relevant to NASA Dryden Flight Research Facility. The efforts involved theoretical development, computer coding, and debugging of the STARS code. New solution procedures were developed in such areas as structures, CFD, and graphics, among others. Furthermore, systems-oriented codes were developed for rendering the code truly multidisciplinary and rather automated in nature. Also, work was performed in pre- and post-processing of engineering analysis data.

Description: Through the use of equivalent variable analogies, the authors demonstrate how an electrical subsystem can be modeled by an equivalent structural subsystem. This allows the electrical subsystem to be probabilistically analyzed by using available structural reliability computer codes such as NESSUS. With the ability to analyze the electrical subsystem probabilistically, we can evaluate the reliability of systems that include both structural and electrical subsystems. Common examples of such systems are a structural subsystem integrated with a health-monitoring subsystem, and smart structures. Since these systems have electrical subsystems that directly affect the operation of the overall system, probabilistically analyzing them could lead to improved reliability and reduced costs. The direct effect of the electrical subsystem on the structural subsystem is of secondary order and is not considered in the scope of this work.

Description: Many methods have been developed to aid analysts in identifying component modes which contribute significantly to component responses. These modes, typically targeted for dynamic model correlation via a modal survey, are known as target modes. Most methods used to identify target modes are based on component global dynamic behavior. It is sometimes unclear if these methods identify all modes contributing to responses important to the analyst. These responses are usually those in areas of hardware design concerns. One method used to check the completeness of target mode sets and identify modes contributing significantly to important component responses is mode participation. With this method, the participation of component modes in dynamic responses is quantified. Those modes which have high participation are likely modal survey target modes. Mode participation is most beneficial when it is used with responses from analyses simulating actual flight events. For spacecraft, these responses are generated via a structural dynamic coupled loads analysis. Using MSC/NASTRAN DMAP, a method has been developed for calculating mode participation based on transient coupled loads analysis results. The algorithm has been implemented to be compatible with an existing coupled loads methodology and has been used successfully to develop a set of modal survey target modes.

Description: When using all MSC/NASTRAN eigensolution methods except Lanczos, the analyst can replace the coupled system rigid-body modes calculated within DMAP module READ with mass orthogonalized and normalized rigid-body modes generated from the system stiffness. This option is invoked by defining MSC/NASTRAN r-set degrees of freedom via the SUPORT bulk data card. The newly calculated modes are required if the rigid-body modes calculated by the eigensolver are not 'clean' due to numerical roundoffs in the solution. When performing transient structural dynamic load analysis, the numerical roundoffs can result in inaccurate rigid-body accelerations which affect steady-state responses. Unfortunately, when using the Lanczos method and defining r-set degrees of freedom, the rigid-body modes calculated within DMAP module REIGL are retained. To overcome this limitation and to allow MSC/NASTRAN to handle SUPORT degrees of freedom identically for all eigensolvers, a DMAP Alter has been written which replaces Lanczos-calculated rigid-body modes with stiffness-generated rigid-body modes. The newly generated rigid-body modes are normalized with respect to the system mass and orthogonalized using the Gram-Schmidt technique. This algorithm has been implemented as an enhancement to an existing coupled loads methodology.

Description: A formal procedure for the probabilistic design assessment of a composite structure is described. The uncertainties in all aspects of a composite structure (constituent material properties, fabrication variables, structural geometry, service environments, etc.), which result in the uncertain behavior in the composite structural responses, are included in the assessment. The probabilistic assessment consists of design criteria, modeling of composite structures and uncertainties, simulation methods, and the decision making process. A sample case is presented to illustrate the formal procedure and to demonstrate that composite structural designs can be probabilistically assessed with accuracy and efficiency.

Description: The objective of this study is to develop an anlysis to predict buckling loads of ring-stiffened anisotropic cylinders subject to axial compression, torsion, and internal pressure. This structure is modeled as a branched shell. A nonlinear axisymmetric prebuckling equilibrium state is assumed which is amenable to an exact solution within each branch. Axisymmetric geometric imperfections are included. Buckling displacements are represented by a Fourier series in the circumferential coordination and the finite-element method in the axial coordinate. Application of the Trefftz criterion to the second variation of the total potential energy leads to a nonlinear eigenvalue problem for the buckling load and mode. Results are presented for both unstiffened and ring-stiffened cylinders in the form of buckling interaction diagrams. Imperfections can cause an unexpected buckling mode in the ring web which would not occur for the perfect structure, and pressurization diminishes the benefit of adding rings to the unstiffened shell to increase the buckling load. The implementation of the analysis methodology into a structural sizing algorithm is discussed.

Description: Recently, there has been much interest in adaptive structures that can respond to a varying environment by changing their properties. Piezoelectric materials and shape memory alloys (SMA) are often used as partial thickness actuators to create such adaptivity by applied energy, usually electric curent. These actuators can be used to inducce strains in a structure and reduce stresses in regions of high stress concentration. Two of the present authors show that axisymmetric actuation strains applied troughout the thickness of a plate with a hole can reduce the stress concentration factor (SCF) in an isotropic plate from 3 to 2. However, in most cases actuators are expected to be bonded to or embedded in the plate, so that the actuation strains are applied in the actuators and not directly in the plate. The objective of this note is to show that such partial-thickness actuation cannot be used to reduce the stress concentration factor with axisymmetric actuations strain distribution.

Description: The effects of blade and root-flexure elasticity and dynamic stall on the stability of hingeless rotor blades are investigated. The dynamic stall description is based on the ONERA models of lift, drag, and pitching moment. The structural analysis is based on three blade models that range from a rigid flap-lag model to two elastic flap-lag-torsion models, which differ in representing root-flexure elasticity. The predictions are correlated with the measured lag damping of an experimental isolated three-blade rotor; the correlation covers rotor operations from near-zero-thrust conditions in hover to highly stalled, high-thrust conditions in foward flight. That correlation shows sensitivity of lag-damping predictions to structural refinements in blade and root-flexure modeling. Moreover, this sensitivity increases with increasing control pitch angle and advance ratio. For high-advance-ratio and high-thrust conditions, inclusion of dynamic stall generally improves the correlation.

Description: The Timoshenko beam theory was used to model cracked beams and to calculate the total strain-energy release rate. The root rotations of the beam segments at the crack tip were estimated based on an approximate two-dimensional elasticity solution. By including the strain energy released due to the root rotations of the beams during crack extension, the strain-energy release rate obtained using beam theory agrees very well with the two-dimensional finite element solution. Numerical examples were given for various beam geometries and loading conditions. Comparisons with existing beam models were also given.

Description: A computational methodology is described to probabilistically simulate the stress concentration factors (SCF's) in composite laminates. This new approach consists of coupling probabilistic composite mechanics with probabilistic finite element structural analysis. The composite mechanics is used to probabilistically describe all the uncertainties inherent in composite material properties, whereas the finite element is used to probabilistically describe the uncertainties associated with methods to experimentally evaluate SCF's, such as loads, geometry, and supports. The effectiveness of the methodology is demonstrated by using is to simulate the SCF's in three different composite laminates. Simulated results match experimental data for probability density and for cumulative distribution functions. The sensitivity factors indicate that the SCF's are influenced by local stiffness variables, by load eccentricities, and by initial stress fields.

Description: The linear elastic response is determined for an internally pressurized, long circular cylindrical shell stiffened on the inside by a regular arrangement of identical stringers and identical rings. Periodicity of this configuration permits the analysis of a portion of the shell wall centered over a generic stringer-ring joint; i.e., a unit cell model. The stiffeners are modeled as discrete beams, and the stringer is assumed to have a symmetrical cross section and the ring an asymmetrical section. Asymmetery causes out-of-plane bending and torsion of the ring. Displacements are assumed as truncated double Fourier series plus simple terms in the axial coordinate to account for the closed and pressure vessel effect (a non-periodic effect). The interacting line loads between the stiffeners and the inside shell wall are Lagrange multipliers in the formulation, and they are also assumed as truncated Fourier series. Displacement continuity constraints between the stiffeners and shell along the contact lines are satisfied point-wise. Equilibrium is imposed by the principle of virtual work. A composite material crown panel from the fuselage of a large transport aircraft is the numerical example. The distributions of the interacting line loads, and the out-of-plane bending moment and torque in the ring, are strongly dependent on modeling the deformations due to transverse shear and cross-sectional warping of the ring in torsion. This paper contains the results from the semiannual report on research on 'Pressure Pillowing of an Orthogonally Stiffened Cylindrical Shell'. The results of the new work are illustrated in the included appendix.

Description: In assessing the reliability of a space nuclear propulsion system (SNPS) nozzle, uncertainties associated with the following design parameters were considered: geometry, boundary conditions, material behavior, and thermal and pressure loads. A preliminary assessment of the reliability was performed using NESSUS (Numerical Evaluation of Stochastic Structures Under Stress), a finite-element computer code developed at the NASA Lewis Research Center. The sensitivity of the nozzle reliability to the uncertainties in the random variables was quantified. With respect to the effective stress, preliminary results showed that the nozzle spatial geometry uncertainties have the most significant effect at low probabilities whereas the inner wall temperature has the most significant effect at higher probabilities.

Description: Vehicles, such as the X-15 or the National Aerospace Plane (NASP), traveling at hypersonic speeds through the earth's atmosphere experience aerodynamic heating. The heating can be severe enough that a thermal protection system is required to limit the temperature of the vehicle structure. Although several categories of thermal protection systems are mentioned briefly, the majority of the present paper describes convectively cooled structures for large areas. Convective cooling is a method of limiting structural temperatures by circulating a coolant through the vehicle structure. Efforts to develop convectively cooled structures during the past 30 years, from early engine structures which were intended to be tested on the X-15 to structural panels fabricated and tested under the NASP program, are described. Many of the lessons learned from these research efforts are presented.

Description: The primary effort concerned the development of analytical methods for the accurate prediction of the effect of random loading on a panel with a crack. Of particular concern was the influence of frequency on the stress intensity factor behavior. Many modern structures, such as those found in advanced aircraft, are lightweight and susceptible to critical vibrations, and consequently dynamic response plays a very important role in their analysis. The presence of flaws and cracks can have catastrophic consequences. The stress intensity factor, K, emerges as a very significant parameter that characterizes the crack behavior. In analyzing the dynamic response of panels that contain cracks, the finite element method is used, but because this type of problem is inherently computationally intensive, a number of ways of calculating K more efficiently are explored.

Description: Adaptive composite structures using actuation materials, such as piezoelectric fibers, were assessed probabilistically utilizing intraply hybrid composite mechanics in conjunction with probabilistic composite structural analysis. Uncertainties associated with the actuation material as well as the uncertainties in the regular (traditional) composite material properties were quantified and considered in the assessment. Static and buckling analyses were performed for rectangular panels with various boundary conditions and different control arrangements. The probability density functions of the structural behavior, such as maximum displacement and critical buckling load, were computationally simulated. The results of the assessment indicate that improved design and reliability can be achieved with actuation material.

Description: A new finite element for coupling built-up shell substructures is presented. The present work extends the hybrid variational formulation of the interface element developed by Aminpour and Ransom to permit coupling between two intersecting substructures. Designed for the assembly of independently built-up finite element models, this technique provides a level of modeling flexibility previously unavailable.

Description: The NASA-UVA Light Aerospace Alloy and Structures Technology (LA2ST) Program was initiated in 1986, and continues a high level of activity, with projects being conducted by graduate students and faculty advisors in the Departments of Materials Science and Engineering, and Mechanical and Aerospace Engineering at the University of Virginia. This work is funded by the NASA-Langley Research Center under Grant NAG-1-745. Here, we report on progress achieved between July 1 and December 31, 1993. The objective of the LA2ST Program is to conduct interdisciplinary graduate student research on the performance of next generation, light weight aerospace alloys, composites and thermal gradient structures in collaboration with NASA-Langley researchers. Specific technical objectives are presented for each research project. We generally aim to produce relevant data and basic understanding of material mechanical response, environmental/corrosion behavior, and microstructure; new monolithic and composite alloys; advanced processing methods; new solid and fluid mechanics analyses; measurement and modeling advances; and critically, a pool of educated graduate students for aerospace technologies.

Description: This paper is a summary of two NASA contracts on high temperature fatigue crack propagation in metals. The first evaluated the ability of fairly simple nonlinear fracture parameters to correlate crack propagation. Hastelloy-X specimens were tested under isothermal and thermomechanical cycling at temperatures up to 980 degrees C (1800 degrees F). The most successful correlating parameter was the crack tip opening displacement derived from the J-integral. The second evaluated the ability of several path-independent integrals to correlate crack propagation behavior. Inconel 718 specimens were tested under isothermal, thermomechanical, temperature gradient, and creep conditions at temperatures up to 650 degrees C (1200 degrees F). The integrals formulated by Blackburn and by Kishimoto correlated the data reasonably well under all test conditions.

Description: This document contains presentations and abstracts from the Workshop on Scaling Effects in Composite Materials and Structures jointly sponsored by NASA Langley Research Center, Virginia Tech, and the Institute for Mechanics and Materials at the University of California, San Diego, and held at NASA Langley on November 15-16, 1993. Workshop attendees represented NASA, other government research labs, the aircraft/rotorcraft industry, and academia. The workshop objectives were to assess the state-of-technology in scaling effects in composite materials and to provide guidelines for future research.

Description: A brief discussion of oxygen-extraction technology and equipment developed by the Space Engineering Research Center (SERC) is presented. Lunar and martian oxygen production from in situ materials and the development of oxygen factories are covered briefly.

Description: The boundary element method is utilized in this study to conduct thermal analysis of functionally graded composites, materials in which the internal microstructure or properties are explicitly tailored in order to obtain an optimal response, on the micromechanical (constituent) scale. A unique feature of the boundary element formulations used here is the use of circular shape functions to convert the two-dimensional integrations of the composite fibers to one dimensional integrations. Using the computer code BEST-CMS, the through the thickness temperature profiles are computed for a representative material with varying numbers of fibers and fiber spacing in the thickness direction. The computed temperature profiles are compared to those obtained using an alternate analytical theory which explicitly couples the heterogeneous microstructure to the global analysis. The boundary element results compared favorably to the analytical calculations, with discrepancies that are explainable based on the boundary element formulation. The results serve both to demonstrate the ability of the boundary element method to analyze these types of materials, and to verify the accuracy of the analytical theory.

Description: Computational methods are described for simulating the formation and the material behavior of ice in prevailing transient environments. The methodology developed at the NASA Lewis Research Center was adopted. A three dimensional finite-element heat transfer analyzer was used to predict the thickness of ice formed under prevailing environmental conditions. A multi-factor interaction model for simulating the material behavior of time-variant ice layers is presented. The model, used in conjunction with laminated composite mechanics, updates the material properties of an ice block as its thickness increases with time. A sample case of ice formation in a body of water was used to demonstrate the methodology. The results showed that the formation and the material behavior of ice can be computationally simulated using the available composites technology.

Description: Verification of a dynamic model of a constrained structure requires a modal survey test of the physical structure and subsequent modification of the model to obtain the best agreement possible with test data. Constrained-boundary or fixed-base testing has historically been the most common approach for verifying constrained mathematical models, since the boundary conditions of the test article are designed to match the actual constraints in service. However, there are difficulties involved with fixed-base testing, in some cases making the approach impractical. It is not possible to conduct a truly fixed-base test due to coupling between the test article and the fixture. In addition, it is often difficult to accurately simulate the actual boundary constraints, and the cost of designing and constructing the fixture may be prohibitive. For use when fixed-base testing proves impractical or undesirable, alternate free-boundary test methods have been investigated, including the residual flexibility technique. The residual flexibility approach has been treated analytically in considerable detail and has had limited frequency response measurements for the method. This concern is well-justified for a number of reasons. First, residual flexibilities are very small numbers, typically on the order of 1.0E-6 in/lb for translational diagonal terms, and orders of magnitude smaller for off-diagonal values. This poses difficulty in obtaining accurate and noise-free measurements, especially for points removed from the excitation source. A second difficulty encountered in residual measurements lies in obtaining a clean residual function in the process of subtracting synthesized modal data from a measured response function. Inaccuracies occur since modes are not subtracted exactly, but only to the accuracy of the curve fits for each mode; these errors are compounded with increasing distance from the excitation point. In this paper, the residual flexibility method is applied to a simple structure in both test and analysis. Measured and predicted residual functions are compared, and regions of poor data in the measured curves are described. It is found that for accurate residual measurements, frequency response functions having prominent stiffness lines in the acceleration/force format are needed. The lack of such stiffness lines increases measurement errors. Interface drive point frequency respose functions for shuttle orbiter payloads exhibit dominant stiffness lines, making the residual test approach a good candidate for payload modal tests when constrained tests are inappropriate. Difficulties in extracting a residual flexibility value from noisy test data are discussed. It is shown that use of a weighted second order least-squares curve fit of the measured residual function allows identification of residual flexibility that compares very well with predictions for the simple structure. This approach also provides an estimate of second order residual mass effects.

Description: International technical experts in durability and damage tolerance of metallic airframe structures were assembled to present and discuss recent research findings and the development of advanced design and analysis methods, structural concepts, and advanced materials. The symposium focused on the dissemination of new knowledge and the peer-review of progress on the development of advanced methodologies. Papers were presented on: structural concepts for enhanced durability, damage tolerance, and maintainability; new metallic alloys and processing technology; fatigue crack initiation and small crack effects; fatigue crack growth models; fracture mechanics failure, criteria for ductile materials; structural mechanics methodology for residual strength and life prediction; development of flight load spectra for design and testing; and advanced approaches to resist corrosion and environmentally assisted fatigue.

Description: This document is the final report for the project entitled 'Thermal Structure Analyses for CSM Testbed (COMET),' for the period of May 16, 1992 - August 15, 1994. The project was focused on the investigation and development of finite element analysis capability of the computational structural mechanics (CSM) testbed (COMET) software system in the field of thermal structural responses. The stages of this project consisted of investigating present capabilities, developing new functions, analysis demonstrations, and research topics. The appendices of this report list the detailed documents of major accomplishments and demonstration runstreams for future references.

Description: The Elastic Plastic Fracture Mechanics Methodology has evolved significantly in the last several years. Nevertheless, some of these concepts need to be extended further before the whole methodology can be safely applied to structural parts. Specifically, there is a need to include the effect of constraint in the characterization of material resistance to crack growth and also to extend these methods to the case of 3D defects. As a consequence, this project was started as a 36 month research program with the general objective of developing an elastic plastic fracture mechanics methodology to assess the structural reliability of pressure vessels and other parts of interest to NASA which may contain flaws. The project is divided into three tasks that deal with (1) constraint and thickness effects, (2) three-dimensional cracks, and (3) the Leak-Before-Burst (LBB) criterion. This report period (March 1994 to August 1994) is a continuation of attempts to characterize three dimensional aspects of fracture present in 'two dimensional' or planar configuration specimens (Chapter Two), especially, the determination of, and use of, crack face separation data. Also, included, are a variety of fracture resistance testing results (J(m)R-curve format) and a discussion regarding two materials of NASA interest (6061-T651 Aluminum alloy and 1N718-STA1 nickel-base super alloy) involving a bases for like constraint in terms of ligament dimensions, and their comparison to the resulting J(m)R-curves (Chapter Two).

Description: Silicon based power devices can be used at 200 C. The device measurements made during this program show a predictable shift in device parameters with increasing temperature. No catastrophic or abrupt changes occurred in the parameters over the temperature range. As expected, the most dramatic change was the increase in leakage currents with increasing temperature. At 200 C the leakage current was in the milliAmp range but was still several orders of magnitude lower than the on-state current capabilities of the devices under test. This increase must be considered in the design of circuits using power transistors at elevated temperature. Three circuit topologies have been prototyped using MOSFET's and IGBT's. The circuits were designed using zero current or zero voltage switching techniques to eliminate or minimize hard switching of the power transistors. These circuits have functioned properly over the temperature range. One thousand hour life data have been collected for two power supplies with no failures and no significant change in operating efficiency. While additional reliability testing should be conducted, the feasibility of designing soft switched circuits for operation at 200 C has been successfully demonstrated.

Description: A composite wing with spars and bulkheads is used to demonstrate the effectiveness of probabilistic assessment of smart composite structures to control uncertainties in distortions and stresses. Results show that a smart composite wing can be controlled to minimize distortions and to have specified stress levels in the presence of defects. Structural responses such as changes in angle of attack, vertical displacements, and stress in the control and controlled plies are probabilistically assessed to quantify their respective uncertainties. Sensitivity factors are evaluated to identify those parameters that have the greatest influence on a specific structural response. Results show that smart composite structures can be configured to control both distortions and ply stresses to satisfy specified design requirements.

Description: The accuracy of using a one-dimensional analysis to predict frequencies of elastically-coupled highly-twisted rotor blades is addressed. Degrees of freedom associated with shear deformation are statically condensed from the formulation, so the analysis uses only those degrees of freedom associated with classical beam theory. The effects of cross section deformation (warping) are considered, and are shown to become significant for some types of elastic coupling. Improved results are demonstrated for highly-coupled blade structures through account of warping in a local cross section analysis, without explicit inclusion of these effects in the beam analysis. A convergence study is also provided which investigates the potential for improving efficiency of elastically-coupled beam analysis through implementation of a p-version beam finite element.

Description: Desirable characteristics and benefits of design oriented analysis methods are described and illustrated by presenting a synoptic description of the development and uses of the Equivalent Laminated Plate Solution (ELAPS) computer code. ELAPS is a design oriented structural analysis method which is intended for use in the early design of aircraft wing structures. Model preparation is minimized by using a few large plate segments to model the wing box structure. Computational efficiency is achieved by using a limited number of global displacement functions that encompass all segments over the wing planform. Coupling with other codes is facilitated since the output quantities such as deflections and stresses are calculated as continuous functions over the plate segments. Various aspects of the ELAPS development are discussed including the analytical formulation, verification of results by comparison with finite element analysis results, coupling with other codes, and calculation of sensitivity derivatives. The effectiveness of ELAPS for multidisciplinary design application is illustrated by describing its use in design studies of high speed civil transport wing structures.