Measurement of Opening Mode Facesheet/Core Interfacial Fracture Toughness of Sandwich Structure using J-Integral
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Typical damage modes in light honeycomb sandwich structures include face sheet/core disbonding and core fracture, both of which can pose a threat to the structural integrity of a component. These damage modes are of particular interest to aviation certification authorities since several in-service occurrences, such as rudder structural failure and other control surface malfunctions, have been attributed to face sheet/core disbonding. Extensive studies have shown that face sheet/core disbonding and core fracture can lead to damage propagation caused by internal pressure changes in the core. The increasing use of composite sandwich construction in aircraft applications makes it vitally important to understand the effect of ground-air-ground (GAG) cycles and conditions such as maneuver and gust loads on face sheet/core disbonding. The objective of the present study was to use a fracture mechanics based approach developed earlier to evaluate the loading at the disbond front caused by ground-air-ground pressurization and in-plane loading. A honeycomb sandwich panel containing a circular disbond at one face sheet/core interface was modeled with three-dimensional (3D) solid finite elements. The disbond was modeled as a discrete discontinuity and the strain energy release rate along the disbond front was computed using the Virtual Crack Closure Technique (VCCT). Special attention was paid to the pressure-deformation coupling which can decrease the pressure load within the disbonded sandwich section significantly when the structure is highly deformed. The commercial finite element analysis software, Abaqus/Standard, was used for the analyses. The recursive pressure-deformation coupling problem was solved by representing the entrapped air in the honeycomb cells as filled cavities in Abaqus/Standard. The results show that disbond size, face sheet thickness and core thickness are important parameters that determine crack tip loading at the disbond front. Further, the pressure-deformation coupling was found to have an important load decreasing effect [6]. In this paper, a detailed problem description is provided first. Second, the analysis methodology is presented. The fracture mechanics approach used is described and the specifics of the finite element model, including the fluid-filled cavities, are introduced. Third, the initial model verification and validation are discussed. Fourth, the findings from a closely related earlier study [6] are summarized. These findings provided the basis for the current investigation. Fifth, an aircraft ascent scenario from 0 to 12192 m (0 to 40000 ft) is considered and the resulting crack tip loading at the disbond front is determined. In-plane loading to simulate maneuvers and gust conditions are also considered. Sixth, the results are shown for a curved panel, which was used to simulate potential fuselage applications. Finally, a brief summary of observations is presented and recommendations for improvement are provided.
Instabilities associated with fluid handling and operation in liquid rocket propulsion systems and test facilities usually manifest themselves as structural vibrations or some form of structural damage. While the source of the instability is directly related to the performance of a component such as a turbopump, valve or a flow control element, the associated pressure fluctuations as they propagate through the system have the potential to amplify and resonate with natural modes of the structural elements and components of the system. In this paper, the authors have developed an innovative multi-level approach that involves analysis at the component and systems level. The primary source of the unsteadiness is modeled with a high-fidelity hybrid RANS/LES based CFD methodology that has been previously used to study instabilities in feed systems. This high fidelity approach is used to quantify the instability and understand the physics associated with the instability. System response to the driving instability is determined through a transfer matrix approach wherein the incoming and outgoing pressure and velocity fluctuations are related through a transfer (or transmission) matrix. The coefficients of the transfer matrix for each component (i.e. valve, pipe, orifice etc.) are individually derived from the flow physics associated with the component. A demonstration case representing a test loop/test facility comprised of a network of elements is constructed with the transfer matrix approach and the amplification of modes analyzed as the instability propagates through the test loop.
The radial Alfven velocity, V(A), in the magnetosphere was investigated using magnetic field and plasma data from the ISEE-1 spacecraft. The harmonic frequency structure calculated from these data was compared with the observed resonant wave structure in the frequency range 5-100 mHz and over the region L = 3-9. It was found that the frequency of these resonances depend on the A(V), which is a function of the magnetic field strength and plasma density. The toroidal mode resonance frequency structure was modeled for the first four harmonics, with excellent agreement between the modeled and observed harmonics for a passs which does not show a plasmapause gradient. In the presence of a steep plasmapause gradient, the model was found to predict a correspondingly steep increase in the resonance frequencies at the plasmapause, by about 5 mHz for the first harmonic and by 20 mHz for the fourth harmonic. Special consideration is given to the effect of heavy ion mass loading on the modeling.
To assess the dynamic characteristics of a deployable space truss, a finite element model of the Scientific Applications Space Platform (SASP) truss has been formulated. The model incorporates all additional degrees of freedom associated with the pin-jointed members. Comparison of results with SPAR models of the truss show that the joints of the deployable truss significantly affect the vibrational modes of the structure only if the truss is relatively short.
Two very different, highly flexible, space structures proposed during the last decade which involve analytical solution of certain partial differential equations are discussed. The first structure is an array, or hanging chain of aluminum beads which would serve as a convenient communicator if, under the influence of the Earth's gravity gradient, it assumes a local vertical orientation as it circles the Earth. Two passive schemes have been proposed for damping the rigid modes: twist the wire at the ends to provide non-zero moment of inertia about the vertical, thereby inducing relative motion of the two tips during rigid pitch or roll, and thus exercising a damper; and introduce weak lossy springs between the end sections and the main section; thereby providing linear coupling between the springs and all the in-plane (pitch) modes. The second structure is a complete ring of satellites cabled together at synchronous altitude. By a slight increase in altitude, the cable is in tension. This configuration is, however, unstable, and an active feedback control scheme is required to stabilize it. A possible scheme involves local cable length adjustment based on measurement of local altitude and shape variations and their rates.
The Galileo scan platform is controlled in two degrees of freedom. A clock (Spin Bearing) actuator controls the relative position between the rotor and stator, and a cone actuator controls the position between the stator and the platform. Instruments on the platform are required to point to 140 micro-rad accuracy and 50 micro-rad per second stability. The system identification objectives were to identify dominant structural resonance frequencies, mode shapes, and damping ratio which exist in the transfer function between the clock actuator and the gyro sensor; and position the notch filter to limit undesirable actuator torque output to ensure stability and performance.
Autonomous modal identification automates the calculation of natural vibration frequencies, damping, and mode shapes of a structure from experimental data. This technology complements damage detection techniques that use continuous or periodic monitoring of vibration characteristics. The approach shown in the paper incorporates the Eigensystem Realization Algorithm (ERA) as a data analysis engine and an autonomous supervisor to condense multiple estimates of modal parameters using ERA's Consistent-Mode Indicator and correlation of mode shapes. The procedure was applied to free-decay responses of a Space Shuttle tail rudder and successfully identified the seven modes of the structure below 250 Hz. The final modal parameters are a condensed set of results for 87 individual ERA cases requiring approximately five minutes of CPU time on a DEC Alpha computer.
Modern aircraft employ a significant fraction of their weight in composite materials to reduce weight and improve performance. Aircraft aeroservoelastic models are typically characterized by significant levels of model parameter uncertainty due to the composite manufacturing process. Small modeling errors in the finite element model will eventually induce errors in the structural flexibility and mass, thus propagating into unpredictable errors in the unsteady aerodynamics and the control law design. One of the primary objectives of Multi Utility Technology Test-bed (MUTT) aircraft is the flight demonstration of active flutter suppression, and therefore in this study, the identification of the primary and secondary modes for the structural model tuning based on the flutter analysis of MUTT aircraft. The ground vibration test-validated structural dynamic finite element model of the MUTT aircraft is created in this study. The structural dynamic finite element model of MUTT aircraft is improved using the in-house Multi-disciplinary Design, Analysis, and Optimization tool. In this study, two different weight configurations of MUTT aircraft have been improved simultaneously in a single model tuning procedure.
Small modeling errors in the finite element model will eventually induce errors in the structural flexibility and mass, thus propagating into unpredictable errors in the unsteady aerodynamics and the control law design. One of the primary objectives of Multi Utility Technology Test Bed, X-56A, aircraft is the flight demonstration of active flutter suppression, and therefore in this study, the identification of the primary and secondary modes for the structural model tuning based on the flutter analysis of X-56A. The ground vibration test validated structural dynamic finite element model of the X-56A is created in this study. The structural dynamic finite element model of the X-56A is improved using a model tuning tool. In this study, two different weight configurations of the X-56A have been improved in a single optimization run.
Solar radiation pressure on the vibrating antenna structure, temperature gradients induced by solar heating, and stabilizing gravity-gradient torques were considered when the linear regulator theory was used to obtain orientation and shape control of a hoop/column antenna system being considered for the land mobile satellite system. A finite element model of the antenna system which includes all six rigid modes and seven flexible modes was used. Results show that the environmental disturbances affect only the rigid modes of the structure. The effect of solar radiation pressure interacting with the vibrating system is smaller than the solar radiation pressure acting on the thermally deformed structure. To reduce control effort in order to maintain shape and orientation, thermal deformations must be minimized in the preliminary design of the system.
Solar radiation pressure on the vibrating antenna structure, temperature gradients induced by solar heating, and stabilizing gravity-gradient torques were considered when the linear regulator theory was used to obtain orientation and shape control of a hoop/column antenna system being considered for the land mobile satellite system. A finite element model of the antenna system which includes all six rigid modes and seven flexible modes was used. Results show that the environmental disturbances affect only the rigid modes of the structure. The effect of solar radiation pressure interacting with the vibrating system is smaller than the solar radiation pressure acting on the thermally deformed structure. To reduce control effect in order to maintain shape and orientation, thermal deformations must be minimized in the preliminary design of the system.
The effects of interaction between a landing gear and a flexible airplane structure on the behavior of the landing gear and the loads in the structure have been studied by treating the equations of motion of the airplane and the landing gear as a coupled system. The landing gear is considered to have nonlinear characteristics typical of conventional gears, namely, velocity-squared damping, polytropic air-compression springing, and exponential tire force-deflection characteristics. For the case where only two modes of the structure are considered, an equivalent three-mass system is derived for representing the airplane and landing-gear combination, which may be used to simulate the effects of structural flexibility in jig drop tests of landing gears. As examples to illustrate the effects of interaction, numerical calculations, based on the structural properties of two large airplanes having considerably different mass and flexibility characteristics, are presented.
Our research focused on Kelvin waves in the tropical lower stratosphere, and resulted in three papers published or submitted to the Journal of Geophysical Research. The first of these, published in 2002, used temperature data from the Microwave Limb Sounder to examine the amplitude, frequency, phase, and spatial structure of leading modes of Kelvin waves. The second and third, submitted late in 2002 and currently in revision, described the response to Kelvin waves of various trace constituents measured by MLS and CLAES (methane, nitrous oxide, CFC-12, and ozone in the second paper; water vapor in the third paper). Water vapor is a special case because the vertical structure induced by Kelvin waves is convolved with water vapor's seasonally varying vertical profile induced by seasonal variations in temperature at the tropical tropopause. Forward modeling indicated that the vertical resolution of MLS was indeed adequate to capture this complicated structure, yet it was not visible in the MLS data, though the Kelvin wave signature was clear on certain UARS levels. The effects of Kelvin waves on the tropical tropopause and on stratosphere- troposphere exchange cannot be quantified from UARS data because of poor vertical resolution and sensitivity in that region. It is recommended that this analysis be repeated using data from the new MLS and HIRDLS instruments aboard Aura, and that priority be given to fine-scale retrievals of temperature, water vapor, and ozone in the tropical tropopause region.
A number of valuable conclusions can be drawn from this study. First, knockdown factors for a specific fluid are not constant but instead are dependent on the mode shape, although the largest this variability gets is about 10% for LOX, the densest fluid. The factors decrease the most for lower frequency shapes and less for higher ones. It follows, therefore, that mode number mismatch between air and fluid operation becomes not only possible, but common, as a knockdown factor for a particular mode shape may be higher than for another mode shape. Since this is a function of added mass, the mismatch is more prevalent for higher density fluids, but it initiates even for very low density ones. Another important conclusion reached is that it appears that the basic mode shapes of a structure do not change if it is fully symmetric, which includes its geometry and boundary conditions. There is some indication of small changes in the relative magnitudes within the mode shape. This conclusion is evident in the results from the cantilever rectangular plate and the inducer, which are not symmetric, and the fixed-fixed plate and the annular disk, which are. For non-symmetric structures, though, the mode shapes almost universally change for dense fluids, as shown by the very low MAC calculations. For the inducer in particular, the changes follow a trend of reduced parabolic and sine wavelengths with increasing density. It is critical to recognize the change in mode shape for several reasons. First, model updating with modal test becomes problematic if the shapes change. Second, design to avoid resonance is highly critical on the mode shape for modes other than the primary ones, as resonance is only a factor when the excitation shape matches the mode shape. Finally, application of the modal superposition method of forced response analysis is dependent on the use of accurate mode shapes. A more-refined assessment of the "knockdown" factor values and ranges than any previously reported in the literature for a realistic engineering structure is also presented in this paper. This data is of tremendous benefit for preliminary analysis and design, where a quick estimate is necessary. These results are important not just for rocket engine turbomachinery, but for water pumps and turbines, propellers, and any other structure operating in a heavy fluid with dynamic excitation. The clear avenue for future work for this endeavor is to expand the analytical techniques discussed in the literature to develop analytical expressions and justification for the mode shape changes and associated frequency knockdowns. These expressions must be able to accurately predict the functional relationship to the shapes, which will enable accurate tracing of the mode number from vacuum analysis (or testing in air) to analysis and operation in the intended fluid environment.
Subspace and Lanczos iterations have been developed, well documented, and widely accepted as efficient methods for obtaining p-lowest eigen-pair solutions of large-scale, practical engineering problems. The focus of this paper is to incorporate recent developments in vectorized sparse technologies in conjunction with Subspace and Lanczos iterative algorithms for computational enhancements. Numerical performance, in terms of accuracy and efficiency of the proposed sparse strategies for Subspace and Lanczos algorithm, is demonstrated by solving for the lowest frequencies and mode shapes of structural problems on the IBM-R6000/590 and SunSparc 20 workstations.
Coupled waveguides structure proposed as leaky- wave antenna, extending previous coupled-mode theory
A Turbine-Based Combined Cycle (TBCC) dynamic simulation model has been developed to demonstrate all modes of operation, including mode transition, for a turbine-based combined cycle propulsion system. The High Mach Transient Engine Cycle Code (HiTECC) is a highly integrated tool comprised of modules for modeling each of the TBCC systems whose interactions and controllability affect the TBCC propulsion system thrust and operability during its modes of operation. By structuring the simulation modeling tools around the major TBCC functional modes of operation (Dry Turbojet, Afterburning Turbojet, Transition, and Dual Mode Scramjet) the TBCC mode transition and all necessary intermediate events over its entire mission may be developed, modeled, and validated. The reported work details the use of the completed model to simulate a TBCC propulsion system as it accelerates from Mach 2.5, through mode transition, to Mach 7. The completion of this model and its subsequent use to simulate TBCC mode transition significantly extends the state-of-the-art for all TBCC modes of operation by providing a numerical simulation of the systems, interactions, and transient responses affecting the ability of the propulsion system to transition from turbine-based to ramjet/scramjet-based propulsion while maintaining constant thrust.