Empirical Model Development for Predicting Shock Response on Composite Materials Subjected to Pyroshock Loading
No abstract available
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.
No abstract available
An evaluation is made of the experimental results obtained in the 720B airliner flight crash testing of an energy-absorbing passenger seat design, which has been derived through modification of a standard three-passenger commercial aviation seat. The highly instrumented 720B aircraft was crashed at the end of a remotely controlled flight on a specially prepared, gravel-covered ground site. In order to limit acceleration upon impact, the seat tested had been designed and fabricated to collapse under a controlled load. The aircraft impact's shock pulse is absorbed by tubes which crush under load as axial columns.
Shock analysis is typically done via SRS (shock response spectrum). Oftentimes shock analysis is not performed because most loads specifications only supply SRS requirements with no representative shock transients that can be used for analysis. SRS analysis is too conservative. There is no structural dynamics involved as all modes are assumed to have peak responses simultaneously and in-phase with each other. It is possible that all other loads and combination of loads show the design to be good, whereas SRS analysis (1000 g’s are involved) results may require redesign. Multiple shock transients can satisfy a SRS, yet most loads documents don’t provide a representative shock transient that could be used for shock analysis. Shock testing is less common than harmonic and random vibration testing. Engineers typically go into shock tests with very little knowledge of what to expect. A pass/fail type of approach with very little engineering done to prepare for the test or to control the type of shock transient is used. Most engineers will not model the shock test fixture due to lack of knowledge about the details of the test setup. Shock analysis and testing is often an open loop type situation with no test setup modeling, no test correlation, and no control over the test shock transient. In this work, if we can create applied test force transient that comes close to matching calculated force transient then a more reasonable acceleration transient will be created at the test reference accelerometer that better matches operational shocks.
Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.
Design improvements and detailed functional analyses are reviewed to trace the development of a pyroactuated release device with segmented thread design from its intermediate design into one that reduces the levels of shock spectra generated during its operation by 50%. Comparisons of shock output and internal load distribution are presented, along with descriptions of mechanical operation for both designs. Results also show the potential areas where design development activity can gain further progress in lowering actuation shock levels.
During the development stage, in order to design/to size the rocket engine components and to reduce the risks, the local dynamic environments as well as dynamic interface loads must be defined. There are two kinds of dynamic environment, i.e. shock transients and steady-state random and sinusoidal vibration environments. Usually, the steady-state random and sinusoidal vibration environments are scalable, but the shock environments are not scalable. In other words, based on similarities only random vibration environments can be defined for a new engine. The methodology covered in this paper provides a way to predict the shock environments and the dynamic loads for new engine systems and new engine components in the early stage of new engine development or engine nozzle modifications.
The PIP-II Project at FNAL requires the assembly of 3 high-beta 650MHz cryomodules at STFC Daresbury (DL) in the UK. These modules must be safely transported from DL to FNAL in the USA. Previous experience with cryomodule transport was leveraged at both labs to design a transport system to protect the cryomodules during transit. Requirements for the system included mitigation of shocks, drops, and vibrations, and acting as a lifting fixture. It is comprised of a tessellated steel frame which encompasses the module with a wire rope isolator arrangement which the module mounts to. The frame was designed to withstand the weight of the 12.5 tonne cryomodule in various load cases. Details of shock and vibration profiles were obtained from MIL-STD-810H and were used to guide the sizing of the isolators. The frame and the isolation system were analysed via FEA using the shock and vibration profiles as an input. The transport system was found to be suitable for the given isolation, frame stiffness, and lifting code requirements. The frame has been fabricated and successfully load tested at FNAL. It will now be road tested with a dummy cryomodule before undergoing a trial run to DL.
An experimental research program providing basic knowledge and establishing a database on the fluctuating pressure loads produced on aerodynamic surfaces beneath three dimensional shock wave/boundary layer interactions is described. Such loads constitute a fundamental problem of critical concern to future supersonic and hypersonic flight vehicles. A turbulent boundary layer on a flat plate is subjected to interactions with swept planar shock waves generated by sharp fins at angle of attack. Fin angles from 10 to 20 deg at freestream Mach numbers of 3 and 4 produce a variety of interaction strengths from weak to very strong. Miniature Kulite pressure transducers flush-mounted in the flat plate are used to measure interaction-induced wall pressure fluctuations. The distributions of properties of the pressure fluctuations, such as their ring levels, amplitude distributions, and power spectra, are also determined. Measurements were made for the first time in the aft regions of these interactions, revealing fluctuating pressure levels as high as 160 dB. These fluctuations are dominated by low frequency (0-5 kHz) signals. The maximum ring levels in the interactions show an increasing trend with increasing interaction strength. On the other hand, the maximum ring levels in the forward portion of the interactions decrease linearly with increasing interaction sweep back. These ring pressure distributions and spectra are correlated with the features of the interaction flowfield. The unsteadiness of the off-surface flowfield is studied using a new, non-intrusive technique based on the shadow graph method. The results indicate that the entire lambda-shock structure generated by the interaction undergoes relatively low-frequency oscillations. Some regions where particularly strong fluctuations are generated were identified. Fluctuating pressure measurements are also made along the line of symmetry of an axisymmetric jet impinging upon a flat plate at an angle. This flow was chosen as a simple analog to the impinging jet region found in the rear portion of the shock wave/boundary layer interactions under study. It is found that a sharp peak in ring pressure level exists at or near the mean stagnation point. It is suggested that the phenomena responsible for this peak may be active in the swept interactions as well, and may cause the extremely high fluctuating pressures observed in the impinging jet region in the present experimental program.
This paper summarizes some of the efforts to advance hypersonic structures technology through the National Aero-Space Plane (NASP) Structures Technology Maturation Program. The ranges of expected structural loads and results from analysis and test activities are described. Topics briefly covered include shock impingement effects on aerothermal loads, actively-cooled structures concepts and test results to date for leading edges and panels, design and fabrication of a carbon-carbon control surface, a program for predicting performance of metal matrix composites, an analysis and sizing procedure for thermal structures, and some recently-developed test fixtures for seals, thermal insulation, and actively-cooled panels.
A method to determine a bound of structural loads for a spacecraft mounted on a launch vehicle is developed. The method utilizes the interface shock spectra and the relative impedance of the spacecraft and launch vehicle. The method is developed for single-degree-of-freedom models and then generalized to multidegree-of-freedom models.
Three-dimensional numerical investigations on the start-up side load physics for a regeneratively cooled, high-aspect-ratio nozzle were performed. The objectives of this study are to identify the three-dimensional side load physics and to compute the associated aerodynamic side load using an anchored computational methodology. The computational methodology is based on an unstructured-grid, pressure-based computational fluid dynamics formulation, and a transient inlet condition based on an engine system simulation. Computations were performed for both the adiabatic and cooled walls in order to understand the effect of boundary conditions. Finite-rate chemistry was used throughout the study so that combustion effect is always included. The results show that three types of shock evolution are responsible for side loads: generation of combustion wave; transitions among free-shock separation, restricted-shock separation, and simultaneous free-shock and restricted shock separations; along with oscillation of shocks across the lip. Wall boundary conditions drastically affect the computed side load physics: the adiabatic nozzle prefers free-shock separation while the cooled nozzle favors restricted-shock separation, resulting in higher peak side load for the cooled nozzle than that of the adiabatic nozzle. By comparing the computed physics with those of test observations, it is concluded that cooled wall is a more realistic boundary condition, and the oscillation of the restricted-shock separation flow pattern across the lip along with its associated tangential shock motion are the dominant side load physics for a regeneratively cooled, high aspect-ratio rocket engine.
The objective of this effort is to develop a computational methodology to capture the startup side load physics and to anchor the computed aerodynamic side loads with the available data from a regeneratively cooled, high-aspect-ratio nozzle, hot-fired at sea level. The computational methodology is based on an unstructured-grid, pressure-based, reacting flow computational fluid dynamics and heat transfer formulation, a transient 5 s inlet history based on an engine system simulation, and a wall temperature distribution to reflect the effect of regenerative cooling. To understand the effect of regenerative wall cooling, two transient computations were performed using the boundary conditions of adiabatic and cooled walls, respectively. The results show that three types of shock evolution are responsible for side loads: generation of combustion wave; transitions among free-shock separation, restricted-shock separation, and simultaneous free-shock and restricted shock separations; along with the pulsation of shocks across the lip, although the combustion wave is commonly eliminated with the sparklers during actual test. The test measured two side load events: a secondary and lower side load, followed by a primary and peak side load. Results from both wall boundary conditions captured the free-shock separation to restricted-shock separation transition with computed side loads matching the measured secondary side load. For the primary side load, the cooled wall transient produced restricted-shock pulsation across the nozzle lip with peak side load matching that of the test, while the adiabatic wall transient captured shock transitions and free-shock pulsation across the lip with computed peak side load 50% lower than that of the measurement. The computed dominant pulsation frequency of the cooled wall nozzle agrees with that of a separate test, while that of the adiabatic wall nozzle is more than 50% lower than that of the measurement. The computed teepee-like formation and the tangential motion of the shocks during lip pulsation also qualitatively agree with those of test observations. Moreover, a third transient computation was performed with a proportionately shortened 1 s sequence, and lower side loads were obtained with the higher ramp rate.
The effectiveness of a fluid spike as a device to protect leading edges of hypersonic atmospheric flight vehicles from high aerothermal loads produced by complex shock-shock interference is studied. The two-dimensional Navier-Stokes equations are solved using an unstructured cell-centered, fully implicit, flux-difference split algorithm. Adaptively generated unstructured meshes are employed. A type IV shock-shock interference for Mach 8 flow on a cylindrical leading edge with and without a small contraflow supersonic jet (fluid spike) placed at two different locations on the body is solved. A typical flow past a blunt body with a type IV shock-shock interference produces very high pressures and heat fluxes on the leading edge. Present results indicate that a fluid spike displaces the bow shock further in front of the body and modifies the shock-shock interference pattern. This leads to reduced peak pressures and heat fluxes on the body.
This paper, which is presented in three parts, is an analytical study of the behavior of landing gear shock struts, with various types of assumptions for the shock-strut characteristics. The effects of tire springing are neglected. The first part compares the behavior of struts with linear and quadratic damping. The second part considers struts with nonlinear spring characteristics and linear or quadratic damping. The third part treats the oleo-pneumatic strut with air-compression springing without damping and with damping proportional to velocity. It is indicated how the damping factor can be determined by experiment.
Both separate and combined wind tunnel and vibration shaker tests were conducted on two structural panels representative of the Shuttle orbiter in the NASA LaRC 8-foot transonic pressure tunnel (TPT) to determine the effects of combined loads on the thermal protection system (TPS). The primary objective of this test was to provide a combined full-scale load environment and realistic time history of the dynamic pressures, Mach numbers (through transonic), and dynamic structural responses of these panels. The panels were selected from orbiter locations where interactive load sources such as aerodynamic shock waves, turbulent boundary layers, strut-induced vorticity, and substrate deformation combined to provide high bonding loads between the TPS and the orbiter structure. The test panels for this program were highly instrumented with static and dynamic pressure gages, accelerometers, deflectometers, strain gages, Schileren and high-speed photography, and special instrumentation necessary to determine TPS/structure interface loads and tile motions. Two test specimens of each orbiter panel were utilized. Both were high-fidelity representations of the selected orbiter location. The first panel (calibration panel) for each test was utilized to calibrate the tunnel and shaker system as well as provide a data-gathering source for applied static and dynamic loads. The second panel (test panel) for each test was subjected to flight time exposures only so that TPS characteristics as a function of flight mission exposure could be determined. These tests were the last in a series that was suce6ssfully completed and necessary to certify the TPS prior to the first launching of the Shuttle.
The existence of hydromagnetic waves in the mass-loaded solar wind upstream from the bow shock of a comet is well-established both for comet Giacobini-Zinner and for comet Halley. Whereas previous reports have been concerned either with the magnetic field observations or with plasma observations, here observations of the magnetic field with the solar wind proton and alpha particle distributions are combined. This allows the three possible modes of propagation for these waves to be separated. The magnetic component is predominantly transverse to the magnetic field and linearly polarized. The flow vector also has a substantial amount of power parallel to the magnetic field. An examination of the pressure variations shows that slow magnetosonic waves are more common than the fast mode.
Understanding the role of microstructural heterogeneities on the shock wave propagation and defect evolution behavior is essential to predicting the dynamic response of metals. Heterogeneities, such as voids, provide challenges to understanding the wave propagation behavior as the shock-void interaction can collapse the void and result in large plastic strains and significant heating (hot spot formation) in the metal. Accurate modeling of this phenomenon requires predicting the void collapse mechanisms and the related heat generation and dissipation mechanisms (hotspot formation) that determine the microstructure evolution. While molecular dynamics (MD) simulations can model the void collapse behavior, the time/length scale capabilities pose a challenge to connect with continuum models or the experimental scales. Here, this study presents the capability of the newly developed quasi-coarse-grained dynamics (QCGD) method that extends the MD simulations to larger system sizes and longer times to model this phenomenon. This study uses QCGD simulations to investigate the mechanisms of shock wave interactions with pre-existing voids in single-crystal Ta microstructures for variations in shock pressures, void size, and loading orientations. For a given orientation, the plasticity contributions and rates of void collapse are observed to vary with shock pressures and void size. The larger void sizes and higher pressures result in increased temperatures (hot spots) and faster void collapse rates and unravel the variations in the plasticity contributions. In addition, QCGD simulations investigate the post-collapse microstructure evolution as a release wave travels through the hot spot region. The simulations reveal that the reduction in temperatures due to heat dissipation initiates the dynamic recrystallization behavior in the hotspot regions.
Mussels use byssal threads to secure themselves to rocks and as shock absorbers during cyclic loading from wave motion. Byssal threads combine high strength and toughness with extensibility of nearly 200%. Researchers attribute tensile properties of byssal threads to their elaborate multi-domain collagenous protein cores. Because the elastic properties have been previously scrutinized, we instead examined byssal thread viscoelastic behaviour, which is crucial for withstanding cyclic loading. By targeting protein domains in the collagenous core via chemical treatments, stress relaxation experiments provided insights on domain contributions and were coupled with in situ small-angle X-ray scattering to investigate relaxation-specific molecular reorganizations. Results show that when silk-like domains in the core were disrupted, the stress relaxation of the threads decreased by nearly 50% and lateral molecular spacing also decreased, suggesting that these domains are essential for energy dissipation and assume a compressed molecular rearrangement when disrupted. A generalized Maxwell model was developed to describe the stress relaxation response. The model predicts that maximal damping (energy dissipation) occurs at around 0.1 Hz which closely resembles the wave frequency along the California coast and implies that these materials may be well adapted to the cyclic loading of the ambient conditions.