ALBERT

Description: The problem of the hypersonic double ellipse in rarefied flow is treated by a particle method using the collision model first described by McDonald (1988). In the approach used here, the computational overhead is reduced by using simple cubic cells. The problem of the definition of complex geometries is addressed by developing an algorithm to define the relation of a body surface to the network of cells.

Description: Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D viscous hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The code has been applied to Problem 8.2 of this workshop for the case of thermochemical nonequilibrium flow through a nozzle. Chemical reaction rates are defined with the model of Park (1987). Thermal nonequilibrium is modeled using a two-temperature approximation in which the vibrational energies of all molecules are assumed to be in equilibrium at a single temperature which is generally different from the translational-rotational temperature. Two grids were used to define the flow for the original problem, with a stagnation temperature of 6500 K. A third case with a stagnation temperature of 10,000 K is also presented. The solution domain includes the converging nozzle, subsonic flow domain in which the gas is substantially in thermochemical equilibrium and the diverging nozzle, hypersonic flow domain in which the gas is substantially in thermochemical nonequilibrium.

Description: Solutions have been computed and results are presented for Problem 1, the case of Mach 9 transitional flow past a 7 deg half-angle cone at zero incidence. The solutions were computed using a code developed for the integration of the parabolized Navier-Stokes equations. The algorithm employed in the code is based on a Roe-type flux-difference-splitting scheme applied following a finite-volume approach. The basic algorithm has been modified to make it implicit and second-order accurate in the crossflow directions. Results are presented in terms of surface pressure and heat transfer as well as boundary layer profiles of pitot pressure, Mach number, and tangential velocity. The case was recalculated several times in an effort to determine sensitivities to such parameters as grid density, wall temperature, turbulence model parameters, as well as freestream expansion. Comparisons with the experimental data are presented and discussed.

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

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

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

Description: Airloads measured on a two-bladed helicopter rotor in flight, from the Tip Aerodynamic and Acoustic Test, are compared with calculations from a comprehensive helicopter analysis (CAMRAD/JA), and the pressures compared with calculations from a full-potential rotor code (FPR). The flight test results cover an advance ratio range from 0.19 to 0.38. The lowest speed case is characterized by the presence of significant blade-vortex interactions. Good correlation of peak-to-peak vortex-induced loads and the corresponding pressures is obtained. The results of the correlation for this two-bladed rotor are substantially similar to the results for three- and four-bladed rotors, concerning the tip vortex core size for best correlation, calculation of the peak-to-peak loads on the retreating side, and calculation of vortex-induced loads on inboard radial stations.

Description: A new CFD potential code, FPX (eXtended Full-Potential), has been developed for application to both helicopters and tilt-rotors. The code solves the unsteady, three-dimensional full potential equation and is an extension of the rotor code, FPR. Both entropy and viscosity corrections are included to enhance the physical modeling capabilities. A number of efficiency related modifications have yielded a factor of two speed-up in the code. An axial flow capability has been added to treat tilt-rotor in forward flight (cruise mode). In order to employ streamwise periodicity and accurately solve for the propagation of acoustic signals in the tip region, an H-H topology has been added to the basic O-H grid system. Computations are performed for the XV-15 Standard and ATB blades at high-speed conditions. Comparisons are made for the blade aerodynamics and the induced fuselage cabin pressure for a range of Mach numbers. Grid generation, wake treatment, and far-field wall treatment are identified as problem areas with recommendations for future research.

Description: A preliminary test/theory correlation evaluation is conducted for wake measurement test results obtained by LDV for a B360 helicopter rotor, at conditions critical to the understanding of wake-rollup and blade-vortex interaction phenomena. The LDV data were complemented by acoustic, blade pressure, rotor performance, and blade/control load measurements.

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

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

Description: Computations were made for those test cases of Problem 3 which were designated as laminar flows, viz., test cases 3.1, 3.2, 3.4, and 3.5. These test cases corresponded to flows over a flat plate and a compression ramp at high Mach number and at high Reynolds number. The computations over the compression ramps indicate a substantial streamwise extent of separation. Based on previous experience with separated laminar flows at high Mach numbers which indicated a substantial effect with spatial grid refinement, a series of computations with different grid sizes were performed. Also, for the flat plate, comparisons of the results for two different algorithms were made.

Description: A development status evaluation is presented for CFD methods applicable to fuselage-integrated scramjet powerplant incorporating hypersonic vehicles; these methods are critically important due to the unavailability of experimental facilities for such elevated Mach number/high-enthalphy conditions. Advancements are required in algorithm robustness and speed, geometric flexibility, and the inclusion of more complete flow physics. The most serious deficiencies lie in turbulence modeling, the lack of complete transition-prediction methods, and combustion modeling.

Description: An overview is given of research activity on the application of computational fluid dynamics (CDF) for hypersonic propulsion systems. After the initial consideration of the highly integrated nature of air-breathing hypersonic engines and airframe, attention is directed toward computations carried out for the components of the engine. A generic inlet configuration is considered in order to demonstrate the highly three dimensional viscous flow behavior occurring within rectangular inlets. Reacting flow computations for simple jet injection as well as for more complex combustion chambers are then discussed in order to show the capability of viscous finite rate chemical reaction computer simulations. Finally, the nozzle flow fields are demonstrated, showing the existence of complex shear layers and shock structure in the exhaust plume. The general issues associated with code validation as well as the specific issue associated with the use of CFD for design are discussed. A prognosis for the success of CFD in the design of future propulsion systems is offered.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: A central artificial-viscosity and an upwind-biased difference method are contrived to solve the Euler equations for flowfields over typical spacecrafts. The spatial discretization is based on either nodal or cell-vertex formulation in the domain extending from free stream to the end of the vehicle. The outer boundary is treated as a bow shock in the first method but is placed in the free stream in the second, which captures both bow and internal shocks using an approximate Riemann solver based on high-order extrapolation to the cell face. These methods were tested for the Shuttle and Hermes orbiters at wind-tunnel conditions and angles of attack ranging from 0 to 60 deg. The artificial-viscosity method incorporated with a shock-fitting procedure shows smeared crossflow and wing-shock positions and required 15 percent more CPU per node than the upwind method. Greater flexibility and robustness is demonstrated by the latter on a fixed grid for all cases considered.

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

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

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

Description: Research requirements to an ultra-high-Reynolds-number liquid helium facility are reviewed. Aerodynamic research areas under consideration include wave vortex hazard reduction, vortex control and diagnostics for maneuvering fighter aircraft, and performance of high-lift devices.

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

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

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

Description: A static aeroelastic analysis capability that calculates flexible air loads for generic configuration wings was developed. It was made possible by integrating a finite element structural analysis code (MSC/NASTRAN) and a panel code of aerodynamic analysis based on linear potential flow theory. The framework already built in MSC/NASTRAN was used, and the aerodynamic influence coefficient matrix was computed externally and inserted in the NASTRAN by means of a DMAP program. It was shown that deformation and flexible air loads of an oblique wing configuration including asymmetric wings can be calculated reliably by this code both in subsonic and supersonic speeds.

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

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

Description: A recent study of total-pressure probes for use in highly turbulent streams is extended herein by developing probe systems that measure time-averaged static or ambient pressure and turbulence intensity. Arrangements of tubular probes of circular and elliptical cross section are described that measure the pressure at orifices on the sides of the probes to obtain different responses to the cross-stream velocity fluctuations. When the measured data are combined to remove the effect of the presence of the probes on the local pressure, the time-averaged static pressure and the cross-stream components of turbulence intensity can be determined. If a system of total pressure tubes, as described in an accompanying paper, is added to the static pressure group to form a single cluster, redundant measurements are obtained that permit accuracy and consistency checks.

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

Description: Flow over slender prolate spheroids at incidence is examined. The incidence angle is chosen high enough to cause streamwise separation of the flow in addition to crossflow separation generally observed at lower incidence angles. The freestream Mach number for the cases investigated here is subsonic, thus precluding the use of parabolized procedures. Laminar, transitional and turbulent flow cases are investigated.

Description: NASA's Aeroassist Flight Experiment (AFE) vehicle will be deployed from the Space Shuttle Orbiter in 1994 to make a data-gethering aeropass through the upper atmosphere before returning to orbit for Shuttle pickup. An axisymmetric, chemically-reacting viscous shock-layer code is presently used to calculate AFE heating rates which automatically accounts for the viscous-inviscid interaction and entropy layer-swallowing effects which are ignored by the conventional boundary-layer methods. Results are presented for the stagnation-point heating of the current AFE baseline trajectory.

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

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

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

Description: An unstructured-grid, finite-volume method has been developed for simulating the inviscid flow over spacecrafts of realistic configuration. The grid generation is accomplished by a new technique on the basis of the advancing-front concept. This simple technique is shown to be equally as powerful for a complex multibody as for a single vehicle. Second- or third-order accuracy is obtained via an innovative interpolation procedure similar to the conventional MUSCL approach. This method has been applied to the Shuttle orbiter and a representative Shuttle launch vehicle consisting of the orbiter, the external tank, and the solid rocket boosters. A comparison is discussed between the present results and other results obtained from structured- and unstructured-grid methods.

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

Description: Finite-difference approximations for steady-state compressible Navier-Stokes equations, whose two spatial dimensions are written in generalized curvilinear coordinates and strong conservation-law form, are presently solved by means of Newton's method in order to obtain a lifting-airfoil flow field under subsonic and transonnic conditions. In addition to ascertaining the computational requirements of an initial guess ensuring convergence and the degree of computational efficiency obtainable via the approximate Newton method's freezing of the Jacobian matrices, attention is given to the need for auxiliary methods assessing the temporal stability of steady-state solutions. It is demonstrated that nonunique solutions of the finite-difference equations are obtainable by Newton's method in conjunction with a continuation method.

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

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

Description: An engineering approach was used to include the nonlinear effects of thickness and camber in an analytical aeroelastic analysis of casecades in supersonic acial flow (supersonic leading-edge locus). A hybrid code using Lighthill's nonlinear piston theory and Lane's linear potential theory was developed to include these nonlinear effects. Lighthill's theory was used to calculate the unsteady pressures on the noninterference surface regions of the airfoils in cascade. Lane's theory was used to calculate the unsteady pressures on the remaining interference surface regions. Two airfoil profiles were investigated (a supersonic throughflow fan design and a NACA 66-206 airfoil with a sharp leading edge). Results show that compared with predictions of Lane's potential theory for flat plates, the inclusion of thickness (with or without camber) may increase or decrease the aeroelastic stability, depending on the airfoil geometry and operating conditions. When thickness effects are included in the aeroelastic analysis, inclusion of camber will influence the predicted stability in proportion to the magnitude of the added camber. The critical interblade phase angle, depending on the airfoil profile and operating conditions, may also be influenced by thickness and camber. Compared with predictions of Lane's linear potential theory, the inclusion of thickness and camber decreased the aerodynamic stiffness and increased the aerodynamic damping at Mach 2 and 2.95 for a cascade of supersonic throughflow fan airfoils oscillating 180 degrees out of phase at a reduced frequency of 0.1.

Description: An efficient and highly accurate algorithm based on a spectral collocation method is developed for numerical solution of the compressible, two-dimensional and axisymmetric boundary layer equations. The numerical method incorporates a fifth-order, fully implicit marching scheme in the streamwise (timelike) dimension and a spectral collocation method based on Chebyshev polynomial expansions in the wall-normal (spacelike) dimension. The spectral collocation algorithm is used to derive the nonsimilar mean velocity and temperature profiles in the boundary layer of a 'fuselage' (cylinder) in a high-speed (Mach 5) flow parallel to its axis. The stability of the flow is shown to be sensitive to the gradual streamwise evolution of the mean flow and it is concluded that the effects of transverse curvature on stability should not be ignored routinely.

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

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

Description: Supersonic flows over a sharp and a flat-faced blunt fin mounted on a flat plate are simulated numerically. Several basic issues involved in the resultant three-dimensional steady flow separation are studied. Using the same number of grid points, different grid spacings are employed to investigate the effects of a grid resolution on the origin of the line of separation. Various shock strengths are used to study the so-called separation and unseparated boundary layer and to establish the existence or absence of secondary separation. The length of separation ahead of the flat-faced blunt fin, bifurcation of a horseshoe vortex, and the accessibility of a closed-type separation are investigated. The usual interpretation of the flow field from previous studies and new interpretations arising from the present simulation are discussed.

Description: The Vaisala-Brunt frequencies of atmospheric waves lie sufficiently near He-filled balloon buoyancy oscillation frequencies to cause a near-resonance condition, especially in the 11.6-24.4 km altitude range. This excitation will hold for such balloons irrespective of size and payload variations due to the independence of balloon buoyant frequencies from size. Vertical perturbations of approximately 5-min period should accordingly be anticipated by balloon flight-control system designers.

Description: Numerical solutions of the thin-layer approximation of the 3D Navier-Stokes equations are presented for flows around an ogive-cylinder body with and without a splitter plate. It is suggested that the presence of the splitter plate prevents the interaction between flows on either side of the symmetry plane. It is concluded that, as a result of the enforced symmetry, the antisymmetric mode of the convective instability near the apex of the body cannot be excited and, therefore, the vortices remain symmetric, staying low and parallel to the upper body surface. The antisymmetric mode of the absolute instability mechanism cannot be initiated, which suppresses the alternate shedding of vortices from the cylindrical portion of the body. High-frequency fluctuations of the shear layer are found to remain virtually unaffected by the presence of the splitter plate.

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

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

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

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

Description: Improved algorithms for the solution of the three-dimensional, time-dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured-grid flow solvers. The spatial discretization involves a flux-split approach that is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time-integration scheme using a multistage Runge-Kutta procedure or an implicit time-integration scheme using a Gauss-Seidel relaxation procedure, which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the Office National d'Etudes et de Recherches Aerospatiales M6 wing to demonstrate applications of the new Euler solvers. The paper presents a description of the Euler solvers along with results and comparisons that assess the capability.

Description: An extension is undertaken of a previous numerical study in order to improve current understanding of nonequilibrium flow effects over slender bodies. Upon extension of the parametric study to encompass smaller nose radii, the downstream influence of equilibrium nonequilibrium flow is found to be much smaller than for the larger nose radii. A comparison of individual stagnation nonequilibrium heating rates demonstrates that a relative comparison of the ratios was not indicative of the actual heating reduction.

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

Description: Hypersonic rarefied flow about the Aeroassist Flight Experiment vehicle has been investigated using a three-dimensional direct simulation Monte Carlo method. Calculations are performed for the transitional flows encountered during the vehicle's atmospheric entry for altitudes of 110 and 100 km with an entry velocity of 9.9 km/s. The simulations are performed using a five-species reacting gas model that account for rotational and vibrational internal energies. The solutions indicate that dissociation is important at altitudes of 110 km and below. Results are presented for surface pressures, convective heating, flowfield structure, and aerodynamic coefficient variations with altitude.

Description: The mean flow structure upstream, around within, and in the near wake of a turbulent junction or horseshoe vortex is reported for an incompressible subsonic flow. Measurements of the primitive variables of velocity and pressure are reported on all surfaces bounding the region of the vortex flow and on three transverse and one streamwise plane within the flowfield itself for comparisons between the measured and any calculated flow variables. Detailed surface flow visualizations and some direct force measurements of surface shear stress are also available. The data is a highly detailed, coherent, self-consistent set offered to the computational fluid mechanics community as a standard test case for the evaluation of the capability of numerical solvers intended for predicting the flowfield in such a complex, separated, three-dimensional turbulent flow. This data base is available for copying to user supplied tapes or for transmission via BITNET, as well as in two National Technical Information Service (NTIS) reports.

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