Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic instability”

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.

At least 73 records · Page 4

Phonon thermodynamics and elastic behavior of GaAs at high temperatures and pressures

The phonons of wurtzite and zinc blende GaAs were calculated at simultaneously elevated temperature and pressure, and elastic constants were calculated as functions of pressure. Pressure caused instabilities of shorter-wavelength transverse acoustic modes in both wurtzite and zinc blende GaAs, causing them to fall to zero at 18 and 20GPa, respectively. Here, the Born stability criteria, which depend on elastic constants and only long wavelength phonons, therefore overestimated the pressure needed to induce instability at 0 K. At elevated temperatures, explicit anharmonicity pushes the onset of instability to higher pressures in both wurtzite and zincblende GaAs. Phonon linewidth and densities of states data showed that the quasiharmonic approximation failed to account for temperature-induced phonon frequency shifts, and the quasiharmonic approximation became less reliable at elevated pressure. In general, the number of three-phonon processes increased with pressure, thereby increasing the temperature-driven broadening of phonon spectral lineshapes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Wrinkling phenomenon in structures. II - Wrinkling criteria

The wrinkling phenomenon in elasticity is mathematically defined and described as a special type of instability associated with the loss of hyperbolicity of the governing equations when some of the characteristic speeds of elastic wave propagation become imaginary. The developed criteria of wrinkling are applied to such structural elements as membranes, strings, soft shells, soft rods, and laminated material. The cause of wrinkling in these elements is due to compression in longitudinal directions that exceeds the through-the-thickness shear modulus. For elements with internal fluid flow the critical compression is reduced due to the destabilizing effect of the flow. The most effective mathematical analysis technique is a combination of tensor analysis and index notations.

Zak, M.↗

Iterated Stretching of Viscoelastic Jets

We examine, with asymptotic analysis and numerical simulation, the iterated stretching dynamics of FENE and Oldroyd-B jets of initial radius r(sub 0), shear viscosity nu, Weissenberg number We, retardation number S, and capillary number Ca. The usual Rayleigh instability stretches the local uniaxial extensional flow region near a minimum in jet radius into a primary filament of radius [Ca(1 - S)/ We](sup 1/2)r(sub 0) between two beads. The strain-rate within the filament remains constant while its radius (elastic stress) decreases (increases) exponentially in time with a long elastic relaxation time 3We(r(sup 2, sub 0)/nu). Instabilities convected from the bead relieve the tension at the necks during this slow elastic drainage and trigger a filament recoil. Secondary filaments then form at the necks from the resulting stretching. This iterated stretching is predicted to occur successively to generate high-generation filaments of radius r(sub n), (r(sub n)/r(sub 0)) = square root of 2[r(sub n-1)/r(sub 0)](sup 3/2) until finite-extensibility effects set in.

Chang, Hsueh-Chia↗

A Multi-Parameter Approach for Calculating Crack Instability

An accurate fracture control analysis of spacecraft pressure systems, boosters, rocket hardware and other critical low-cycle fatigue cases where the fracture toughness highly impacts cycles to failure requires accurate knowledge of the material fracture toughness. However, applicability of the measured fracture toughness values using standard specimens and transferability of the values to crack instability analysis of the realistically complex structures is refutable. The commonly used single parameter Linear Elastic Fracture Mechanics (LEFM) approach which relies on the key assumption that the fracture toughness is a material property would result in inaccurate crack instability predictions. In the past years extensive studies have been conducted to improve the single parameter (K-controlled) LEFM by introducing parameters accounting for the geometry or in-plane constraint effects]. Despite the importance of the thickness (out-of-plane constraint) effects in fracture control problems, the literature is mainly limited to some empirical equations for scaling the fracture toughness data] and only few theoretically based developments can be found. In aerospace hardware where the structure might have only one life cycle and weight reduction is crucial, reducing the design margin of safety by decreasing the uncertainty involved in fracture toughness evaluations would result in lighter hardware. In such conditions LEFM would not suffice and an elastic-plastic analysis would be vital. Multi-parameter elastic plastic crack tip field quantifying developments combined with statistical methods] have been shown to have the potential to be used as a powerful tool for tackling such problems. However, these approaches have not been comprehensively scrutinized using experimental tests. Therefore, in this paper a multi-parameter elastic-plastic approach has been used to study the crack instability problem and the transferability issue by considering the effects of geometrical constraints as well as the thickness. The feasibility of the approach has been examined using a wide range of specimen geometries and thicknesses manufactured from 7075-T7351 aluminum alloy.

Zanganeh, M.↗

INSTABILITY ANALYSIS OF CYLINDRICAL SHELLS UNDER HYDROSTATIC PRESSURE

To determine the elastic buckling pressure of simply supported cylindrical shells subjected to lateral end axial hydrostatic forces, various versions of linear bending theories have been employed in the past. For certain shell dimensions, however the expressions commonly used may yield substantially differing results. In what follows, recent work on this problem by A. E. Armenakas and the writer is briefly reviewed. This work consisted primarily in employing a general bending theory of circular cylindrical shells under the influence of initial stress, developed earlier by the same authors, to re-examine the problem mentioned, and compare the results with those of previous investigations. The outcome was the establishment of a simple but accurate expression for the buckling pressure applicable to a wide range of shell dimensions.

CYLINDRICAL SHELL↗

DESAP 1: A structural design program with stress and displacement constraints. Volume 1: Theoretical and user's manual

A finite element program is presented for computer-automated, minimum weight design of elastic structures with constraints on stresses (including local instability criteria) and displacements. Volume 1 of the report contains the theoretical and user's manual of the program. Sample problems and the listing of the program are included in Volumes 2 and 3. The element subroutines are organized so as to facilitate additions and changes by the user. As a result, a relatively minor programming effort would be required to make DESAP 1 into a special purpose program to handle the user's specific design requirements and failure criteria.

Kiusalaas, J.↗

DESAP 2: A structural design program with stress and buckling constraints. Volume 1: Theoretical and user's manual

DESAP 2 is described as a finite element program for computer-automated, minimum weight design of elastic structures with constraints on stresses (including local instability criteria) and buckling loads. No limits are placed on the number of load conditions for stress-constrained design, but only one of these load conditions can be chosen as the potential buckling load. A substantial portion of DESAP 2, particularly the analysis of the prebuckling state, is derived from the SOLID SAP finite element program. The stress-constrained design is based on the classical stress ratio method, which drives the design towards a fully stressed state. The constraints on the buckling load are handled by solving the appropriate optimality criterion by successive iterations. During each iteration, the element sizes determined by the stress ratio method are used as the minimum size constraints. The element subroutines are organized in a manner that permits the user to make additions and changes with a minimal programming effort. Consequently, DESAP 2 can readily be changed into a special-purpose program to handle the user's specific design requirements and failure criteria.

Kiusalaas, J.↗

Thermoelastic effects in frictionally heated line contact of slablike bodies

Many sealing configurations incorporate the geometry of line contact on bluff, slablike bodies. When sliding speed is sufficiently high along or across the line of contact, instabilities may arise from interactions of thermal expansion, elastic deformation and frictional heating. These lead to concentrated contact with elevated temperatures and pressures. Previous studies have been largely restricted to two-dimensional models of such contact. The present study shows how those earlier results must be modified to apply to the more realistic geometry, and show also that the principal features of the interactions are the same in both geometries.

Heckmann, S. R.↗

A Comparison of Single-Cycle Versus Multiple-Cycle Proof Testing Strategies

Single-cycle and multiple-cycle proof testing (SCPT and MCPT) strategies for reusable aerospace propulsion system components are critically evaluated and compared from a rigorous elastic-plastic fracture mechanics perspective. Earlier MCPT studies are briefly reviewed. New J-integral estimation methods for semi-elliptical surface cracks and cracks at notches are derived and validated. Engineering methods are developed to characterize crack growth rates during elastic-plastic fatigue crack growth (FCG) and the tear-fatigue interaction near instability. Surface crack growth experiments are conducted with Inconel 718 to characterize tearing resistance, FCG under small-scale yielding and elastic-plastic conditions, and crack growth during simulated MCPT. Fractography and acoustic emission studies provide additional insight. The relative merits of SCPT and MCPT are directly compared using a probabilistic analysis linked with an elastic-plastic crack growth computer code. The conditional probability of failure in service is computed for a population of components that have survived a previous proof test, based on an assumed distribution of initial crack depths. Parameter studies investigate the influence of proof factor, tearing resistance, crack shape, initial crack depth distribution, and notches on the MCPT vs. SCPT comparison. The parameter studies provide a rational basis to formulate conclusions about the relative advantages and disadvantages of SCPT and MCPT. Practical engineering guidelines are proposed to help select the optimum proof test protocol in a given application.

McClung, R. C.↗

A Comparison of Single-Cycle Versus Multiple-Cycle Proof Testing Strategies

Single-cycle and multiple-cycle proof testing (SCPT and MCPT) strategies for reusable aerospace propulsion system components are critically evaluated and compared from a rigorous elastic-plastic fracture mechanics perspective. Earlier MCPT studies are briefly reviewed. New J-integral estimation methods for semi-elliptical surface cracks and cracks at notches are derived and validated. Engineering methods are developed to characterize crack growth rates during elastic-plastic fatigue crack growth (FCG) and the tear-fatigue interaction near instability. Surface crack growth experiments are conducted with Inconel 718 to characterize tearing resistance, FCG under small-scale yielding and elastic-plastic conditions, and crack growth during simulated MCPT. Fractography and acoustic emission studies provide additional insight. The relative merits of SCPT and MCPT are directly compared using a probabilistic analysis linked with an elastic-plastic crack growth computer code. The conditional probability of failure in service is computed for a population of components that have survived a previous proof test, based on an assumed distribution of initial crack depths. Parameter studies investigate the influence of proof factor, tearing resistance, crack shape, initial crack depth distribution, and notches on the MCPT vs. SCPT comparison. The parameter studies provide a rational basis to formulate conclusions about the relative advantages and disadvantages of SCPT and MCPT. Practical engineering guidelines are proposed to help select the optimum proof test protocol in a given application.

McClung, R. C.↗

A Comparison of Single-Cycle Versus Multiple-Cycle Proof Testing Strategies

Single-cycle and multiple-cycle proof testing (SCPT and MCPT) strategies for reusable aerospace propulsion system components are critically evaluated and compared from a rigorous elastic-plastic fracture mechanics perspective. Earlier MCPT studies are briefly reviewed. New J-integral estimation methods for semielliptical surface cracks and cracks at notches are derived and validated. Engineering methods are developed to characterize crack growth rates during elastic-plastic fatigue crack growth (FCG) and the tear-fatigue interaction near instability. Surface crack growth experiments are conducted with Inconel 718 to characterize tearing resistance, FCG under small-scale yielding and elastic-plastic conditions, and crack growth during simulated MCPT. Fractography and acoustic emission studies provide additional insight. The relative merits of SCPT and MCPT are directly compared using a probabilistic analysis linked with an elastic-plastic crack growth computer code. The conditional probability of failure in service is computed for a population of components that have survived a previous proof test, based on an assumed distribution of initial crack depths. Parameter studies investigate the influence of proof factor, tearing resistance, crack shape, initial crack depth distribution, and notches on the MCPT versus SCPT comparison. The parameter studies provide a rational basis to formulate conclusions about the relative advantages and disadvantages of SCPT and MCPT. Practical engineering guidelines are proposed to help select the optimum proof test protocol in a given application.

McClung, R. C.↗

Crack Instability Predictions Using a Multi-Term Approach

Present crack instability analysis for fracture critical flight hardware is normally performed using a single parameter, K(sub C), fracture toughness value obtained from standard ASTM 2D geometry test specimens made from the appropriate material. These specimens do not sufficiently match the boundary conditions and the elastic‐plastic constraint characteristics of the hardware component, and also, the crack instability of most commonly used aircraft and aerospace structural materials have some amount of stable crack growth before fracture which makes the normal use of a K(sub C) single parameter toughness value highly approximate. In the past, extensive studies have been conducted to improve the single parameter (K or J controlled) approaches by introducing parameters accounting for the geometry or in‐plane constraint effects. Using 'J‐integral' and 'A' parameter as a measure of constraint is one of the most accurate elastic‐plastic crack solutions currently available. In this work the feasibility of the J‐A approach for prediction of the crack instability was investigated first by ignoring the effects of stable crack growth i.e. using a critical J and A and second by considering the effects of stable crack growth using the corrected J‐delta a using the 'A' parameter. A broad range of initial crack lengths and a wide range of specimen geometries including C(T), M(T), ESE(T), SE(T), Double Edge Crack (DEC), Three‐Hole‐Tension (THT) and NC (crack from a notch) manufactured from Al7075 were studied. Improvements in crack instability predictions were observed compared to the other methods available in the literature.

Zanganeh, Mohammad↗

Acoustic velocities, elasticity, and pressure-induced elastic softening in compressed neodymium

For the first time, the acoustic velocities and elasticity of double hexagonal close-packed neodymium (Nd) were measured at high pressure using ultrasonic interferometry in conjunction with synchrotron X-ray techniques. Here, both the compressional and shear wave velocities increased with an increase in pressure of up to ~5.4 GPa before they exhibited almost constant values upon further compression. This pressure-induced elastic wave velocity softening can be attributed to the structural instability in dhcp-Nd. From the current sound velocities and density data, we obtained new bulk and shear moduli as well as their pressure dependence for polycrystalline dhcp-Nd, yielding $B_0$=33.2(9) GPa, $G_0$=14.1(5) GPa, $\partial$$B/$$\partial$$P$ = 3.6 (4) and $\partial$$=G/$$\partial$$P$ = 1.3 (2) while other elasticity-related/mechanical properties, including Young's modulus, Debye temperature, Poisson's ratio, and Grüneisen parameter, which are important parameters for the application of Nd in extreme environments, were also derived.

36 MATERIALS SCIENCE↗

Aeroelastic stability and control of an oblique wing

Dynamic model tests were staged to compare the aeroelastic stability of a free-flying oblique-winged aircraft and the known divergent instability characterizing swept-forward wing configurations, to test the effect of wing elasticity on aileron control and trim. The aircraft fuselage was clamped in some tests and left free to roll in others. Two rigid wing panels were hinged near the wing center with elastic restraining springs. It is found that the aeroelastic divergence affecting swept-forward wings does not occur with an oblique wing, whose aeroelastic instability sets in at a dynamic pressure exceeding that for static divergence of swept-forward wings, and appears in the form of undamped oscillations. Oblique wing stability is sensitive to the moment of inertia of the fuselage in roll. Roll control and lateral trim are attainable with conventional ailerons, even at speeds at which divergence of a clamped forward wing shows up.

Jones, R. T.↗

Plasticity-mediated deformation instabilities in thin film-compliant substrate systems: direct three-dimensional simulations

Abstract Surface wrinkles driven by mechanical instability commonly form in thin-film structures attached to a compliant substrate. In this study, a recently developed computational approach is employed to simulate the formation and transformation of wrinkles involving plastic yielding of the thin film. The three-dimensional (3D) finite element models contain an embedded imperfection at the film-substrate interface, serving to trigger the bifurcation modes. Successful application of this technique to allow for film plasticity is demonstrated, including the evolution of 3D surface patterns and their correlation with the overall load–displacement response. The simulations reveal that plastic yielding transforms the surface instability patterns into more localized forms. Under uniaxial loading, the sinusoidal elastic wrinkles undergo the wrinkle-to-fold transition. With equi-biaxial loading, the initial square-checkerboard array turns into continuous tall ridges along the 45° directions. In both loading modes, the plasticity-induced instability patterns are only partially relieved upon unloading, leaving permanent features on the surface.

36 MATERIALS SCIENCE↗

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗

Morphological instability in epitaxially strained dislocation-free solid films - Linear stability theory

The morphological instability of a growing epitaxially strained dislocation-free solid film is analyzed. An evolution equation for the film surface is derived in the dilute limit of vacancies based on surface diffusion driven by a stress-dependent chemical potential. From the time-dependent linear stability problem the conditions for which a growing film is unstable are determined. It is found that the instability is driven by the lattice mismatch between the film and the substrate; however, low temperatures as well as elastically stiff substrates are stabilizing influences. The results also reveal that the critical film thickness for instability depends on the growth rate of the film itself. Detailed comparison with experimental observations indicates that the instability described exhibits many of the observed features of the onset of the 'island instability'.

Spencer, B. J.↗