Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic strength”

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 613 records · Page 34

Analysis and Characterization of Damage and Failure Utilizing a Generalized Composite Material Model Suitable for Use in Impact Problems

A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in state-of-the art composite impact models is under development. In particular, a next generation composite impact material model, jointly developed by the FAA and NASA, is being implemented into the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage, and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters (such as modulus and strength). The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in the various coordinate directions. Due to the fact that the plasticity and damage models are uncoupled, test procedures and methods to both characterize the damage model and to covert the material stress-strain curves from the true (damaged) stress space to the effective (undamaged) stress space have been developed. A methodology has been developed to input the experimentally determined composite failure surface in a tabulated manner. An analytical approach is then utilized to track how close the current stress state is to the failure surface.

Finite Element Method↗

Incorporation of Failure Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. In the damage model, a semi-coupled approach is employed where the overall damage in a particular coordinate direction is assumed to be a multiplicative combination of the damage in that direction resulting from the applied loads in various coordinate directions. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

Impact↗

Implementation of a Tabulated Failure Model Into a Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased use in the aerospace and automotive communities. The aerospace community has identified several key capabilities which are currently lacking in the available material models in commercial transient dynamic finite element codes. To attempt to improve the predictive capability of composite impact simulations, a next generation material model is being developed for incorporation within the commercial transient dynamic finite element code LS-DYNA. The material model, which incorporates plasticity, damage and failure, utilizes experimentally based tabulated input to define the evolution of plasticity and damage and the initiation of failure as opposed to specifying discrete input parameters such as modulus and strength. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. For the damage model, a strain equivalent formulation is used to allow for the uncoupling of the deformation and damage analyses. For the failure model, a tabulated approach is utilized in which a stress or strain based invariant is defined as a function of the location of the current stress state in stress space to define the initiation of failure. Failure surfaces can be defined with any arbitrary shape, unlike traditional failure models where the mathematical functions used to define the failure surface impose a specific shape on the failure surface. In the current paper, the complete development of the failure model is described and the generation of a tabulated failure surface for a representative composite material is discussed.

polymer matrix composites↗

The structure and circulation of the deep Venus atmosphere

A simple model for the structure of a nonrotating Hadley regime in an atmosphere with large thermal inertia is developed. The radiative fluxes are estimated by using a linearization about the radiative equilibrium state, and the dynamical fluxes are estimated by using scaling analysis. The requirement that differential heating by these fluxes be in balance in both the meridional and vertical directions leads to two equations for the mean static stability and meridional temperature contrast. The solution depends on two parameters: the strength of the radiative heating, as measured by the static stability of the radiative equilibrium state; and the ratio of the time it takes an external gravity wave to traverse the atmosphere to the time it would take the atmosphere to cool off radiatively. It is shown that it is not necessary to invoke convection to explain the approximate adiabatic lapse rate in the Venus atmosphere, but a greenhouse effect is necessary to explain the high surface temperatures.

Stone, P. H.↗

Nonlinear dynamic theory for photorefractive phase hologram formation

A nonlinear dynamic theory is developed for the formation of photorefractive volume phase holograms. A feedback mechanism existing between the photogenerated field and free-electron density, treated explicitly, yields the growth and saturation of the space-charge field in a time scale characterized by the coupling strength between them. The expression for the field reduces in the short-time limit to previous theories and approaches in the long-time limit the internal or photovoltaic field. Additionally, the phase of the space charge field is shown to be time-dependent.

Kim, D. M.↗

Report on the final panel discussion on computational aeroacoustics

Some important conclusions about future prospects for aeroacoustics in general, and for computational aeroacoustics in particular, that were reached in the course of the Final Panel Discussion of the Workshop on Computational Aeroacoustics held from 6 to 9 April 1992 by ICASE and NASA Langley Research Center are summarized by the panel chairman. Aeroacoustics must now be involved in interactions with computational fluid dynamics (as applied not only to deterministic flows but also to the statistical characteristics of turbulence), while additionally incorporating rigorous comparisons with experiment. The new Computational Aeroacoustics will press forward in two parallel ways. In one of them, CFD will be used to determine aeroacoustic source strengths, the associated radiation being derived by the Acoustic Analogy approach in one of its forms. In the other, a direct Computational Aeroacoustics will apply CFD techniques over a region extending beyond the flow field so as to include at least the beginnings of the acoustic far field. There are some particularly important areas of study, including rotor noise, boundary-layer noise, and the noise of supersonic jets, where it is strongly recommended that use of both methods is continued. On the other hand, important problems of the diffraction of radiation from aeroacoustic sources around complicated aircraft shapes will require the use of comprehensively Computational Aeroacoustics, while Acoustic Analogy methods seem better suited to estimating subsonic jet noise. The study of model problems to allow comparisons with experiment will be valuable in both lines of attack.

Lighthill, James↗

Computation of rotor aerodynamic loads in forward flight using a full-span free wake analysis

The development of an advanced computational analysis of unsteady aerodynamic loads on isolated helicopter rotors in forward flight is described. The primary technical focus of the development was the implementation of a freely distorting filamentary wake model composed of curved vortex elements laid out along contours of constant vortex sheet strength in the wake. This model captures the wake generated by the full span of each rotor blade and makes possible a unified treatment of the shed and trailed vorticity in the wake. This wake model was coupled to a modal analysis of the rotor blade dynamics and a vortex lattice treatment of the aerodynamic loads to produce a comprehensive model for rotor performance and air loads in forward flight dubbed RotorCRAFT (Computation of Rotor Aerodynamics in Forward Flight). The technical background on the major components of this analysis are discussed and the correlation of predictions of performance, trim, and unsteady air loads with experimental data from several representative rotor configurations is examined. The primary conclusions of this study are that the RotorCRAFT analysis correlates well with measured loads on a variety of configurations and that application of the full span free wake model is required to capture several important features of the vibratory loading on rotor blades in forward flight.

Quackenbush, Todd R.↗

Evolution of microwave limb sounder ozone and the polar vortex during winter

The evolution of polar ozone observed by the Upper Atmosphere Research Satellite (UARS) Microwave Limb Sounder (MLS) is described for the northern hemisphere (NH) winters of 1991/1992, 1992/1993, and 1993/1994 and the southern hemisphere (SH) winters of 1992 and 1993. Imterannual and interhemispheric variability in polar ozone evolution are closely related to differences in the polar vortex and to the frequency, duration and strength of stratospheric sudden warmings. Ozone in the midstratospheric vortices increases during the winter, with largest increases associated with stratospheric warmings and a much larger increase in the NH than in the SH. A smaller NH increase was observed in 1993/1994, when the middle stratospheric vortex was stronger. During strong stratospheric warmings in the NH, the upper stratospheric vortex may be so much eroded that it presents little barrier to poleward transport; in contrast, the SH vortex remains strong throughout the stratosphere during wintertime warmings, and ozone increases only below the mixing ratio peak, due to enhanced diabatic descent. Ozone mixing ratios decrease rapidly in the lower stratosphere in both SH late winters, as expected from chemical destruction due to enhanced reactive chlorine. The interplay between dynamics and chemistry is more complex in the NH lower stratosphere and interannual variability is greater. Evidence has previously been shown for chemical ozone destruction in the 1991/1992 and 1992/1993 winters. We show here evidence suggesting some chemical destruction in late February and early March 1994. In the NH late winter lower stratosphere the pattern of high-ozone values (typical of the vortex) seen in mid-latitudes is related to the strength of the lower-stratospheric vortex, with the largest areal extent of high ozone outside the vortex in 1994, when the lower stratospheric vortex is relatively weak, and the least extent in 1993 when the lower stratospheric vortex is strongest.

Manney, G. L.↗

Statics and dynamics of Giacobini-Zinner magnetic tail

The data gathered by the International Cometary Explorer during its traversal of Comet Giacobini-Zinner's tail are subjected to stress balance requirements, yielding estimates of unmeasured quantities. It is noted that the comet's tail is embedded in an ionosheath whose static pressure is nearly equal to the solar wind stagnation pressure, leading to a large lobe field strength. A systematic variation is found in the ion temperature across the tail, implying the variation of the pickup velocity of new ions. Axial stress balance yields an expression for the strength of the lobe field which reveals weak variation with axial distance.

Siscoe, G. L.↗

Design considerations for attaining 250-knot test velocities at the aircraft landing dynamics facility

Preliminary design studies are presented which consider the important parameters in providing 250 knot test velocities at the Aircraft Landing Dynamics Facility. Four major components of this facility are: the hydraulic jet catapult, the test carriage structure, the reaction turning bucket, and the wheels. Using the hydraulic-jet catapult characteristics, a target design point was selected and a carriage structure was sized to meet the required strength requirements. The preliminary design results indicate that to attain 250 knot test velocities for a given hydraulic jet catapult system, a carriage mass of 25,424 kg (56,000 lbm.) cannot be exceeded.

Gray, C. E., Jr.↗

Vapor Cavitation in Dynamically Loaded Journal Bearings

High speed motion camera experiments were performed on dynamically loaded journal bearings. The length to diameter ratio of the bearing, the speed of the roller and the tube, the surface material of the roller, and the static and dynamic eccentricity of the bearing were varied. One hundred and thirty-four cases were filmed. The occurrence of vapor cavitation was clearly evident in the films and figures presented. Vapor cavitation was found to occur when the tensile stress applied to the oil exceeded the tensile strength of the oil or the binding of the oil to the surface. The physical situation in which vapor cavitation occurs is during the squeezing and sliding motion within a bearing. Besides being able to accurately capture the vapor cavitation on film, an analysis of the formation and collapse of the cavitation bubbles and characteristics of the bubble content are presented.

Jacobson, B. O.↗

Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations

Parameter studies are conducted using the Euler and potential flow equation models for unsteady and steady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux-vector-splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flows results are made as well as with experimental data where available.

Anderson, W. Kyle↗

Accurate solutions, parameter studies and comparisons for the Euler and potential flow equations

Parameter studies are conducted using the Euler and potential flow equation models for steady and unsteady flows in both two and three dimensions. The Euler code is an implicit, upwind, finite volume code which uses the Van Leer method of flux vector splitting which has been recently extended for use on dynamic meshes and maintain all the properties of the original splitting. The potential flow code is an implicit, finite difference method for solving the transonic small disturbance equations and incorporates both entropy and vorticity corrections into the solution procedures thereby extending its applicability into regimes where shock strength normally precludes its use. Parameter studies resulting in benchmark type calculations include the effects of spatial and temporal refinement, spatial order of accuracy, far field boundary conditions for steady flow, frequency of oscillation, and the use of subiterations at each time step to reduce linearization and factorization errors. Comparisons between Euler and potential flow results are made, as well as with experimental data where available.

Anderson, W. Kyle↗

Free and forced response of nearly periodic multi-span beams and multi-bay trusses

The influence of several effects on mode localization in free-free multispan beams is investigated. Finite-element methods are used to study localization as a function of Timoshenko beam effects, imperfection strength, coupling strength, and imperfection type. The multispan beam study is conducted with random realizations of a forty-span free-free beam. In addition, free and forced response of free-free trusses and frames is investigated. Results indicate that, for the trusses and frames considered, the dynamic behavior can be approximated by continuum modeling approaches for only the lowest few global modes. Also, for reasonable values of modal damping, localized modes play a subtle role in the overall structural response.

Lust, S. D.↗

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗

Improved Understanding of the Modeled QBO Using MLS Observations and MERRA Reanalysis

The Quasi-Biennial Oscillation (QBO) dominates the variability of the tropical stratosphere on interannual time scales. The QBO has been shown to extend its influence into the chemical composition of this region through dynamical mechanisms. We have started our analysis using the realistic QBO internally generated by the Goddard Earth Observing System Version 5 (GEOS-5) general circulation model coupled to a comprehensive stratospheric and tropospheric chemical mechanism forced with observed sea surface temperatures over the past 33 years. We will show targeted comparisons with observations from NASAs Aura satellite Microwave Limb Sounder (MLS) and the Modern Era Retrospective-Analysis for Research and Applications (MERRA) reanalysis to provide insight into the simulation of the primary and secondary circulations associated with the QBO. Using frequency spectrum analysis and multiple linear regression we can illuminate the resulting circulations and deduce the strengths and weaknesses in their modeled representation. Inclusion of the QBO in our simulation improves the representation of the subtropical barriers and overall tropical variability. The QBO impact on tropical upwelling is important to quantify when calculating trends in sub-decadal scale datasets.

MERRA↗

Mixed Layer Sub-Mesoscale Parameterization - Part 1: Derivation and Assessment

Several studies have shown that sub-mesoscales (SM 1km horizontal scale) play an important role in mixed layer dynamics. In particular, high resolution simulations have shown that in the case of strong down-front wind, the re-stratification induced by the SM is of the same order of the de-stratification induced by small scale turbulence, as well as of that induced by the Ekman velocity. These studies have further concluded that it has become necessary to include SM in ocean global circulation models (OGCMs), especially those used in climate studies. The goal of our work is to derive and assess an analytic parameterization of the vertical tracer flux under baroclinic instabilities and wind of arbitrary directions and strength. To achieve this goal, we have divided the problem into two parts: first, in this work we derive and assess a parameterization of the SM vertical flux of an arbitrary tracer for ocean codes that resolve mesoscales, M, but not sub-mesoscales, SM. In Part 2, presented elsewhere, we have used the results of this work to derive a parameterization of SM fluxes for ocean codes that do not resolve either M or SM. To carry out the first part of our work, we solve the SM dynamic equations including the non-linear terms for which we employ a closure developed and assessed in previous work. We present a detailed analysis for down-front and up-front winds with the following results: (a) down-front wind (blowing in the direction of the surface geostrophic velocity) is the most favorable condition for generating vigorous SM eddies; the de-stratifying effect of the mean flow and re-stratifying effect of SM almost cancel each other out,

Canuto, V. M.↗

Design of helicopter rotor blades for optimum dynamic characteristics

The possibilities and limitations of tailoring blade mass and stiffness distributions to give an optimum blade design in terms of weight, inertia, and dynamic characteristics are discussed. The extent that changes in mass of stiffness distribution can be used to place rotor frequencies at desired locations is determined. Theoretical limits to the amount of frequency shift are established. Realistic constraints on blade properties based on weight, mass, moment of inertia, size, strength, and stability are formulated. The extent that the hub loads can be minimized by proper choice of E1 distribution, and the minimum hub loads which can be approximated by a design for a given set of natural frequencies are determined. Aerodynamic couplings that might affect the optimum blade design, and the relative effectiveness of mass and stiffness distribution on the optimization procedure are investigated.

Peters, D. A.↗