Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Strain calculation”

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 685 records · Page 38

Elevated Temperature Deformation of Fe-39.8Al and Fe-15.6Mn-39.4Al

The elevated temperature compressive properties of binary Fe-39.8 at % Al and Fe-15.6Mn-39.4Al have been measured between 1000 and 1300 K at strain rates between 10(exp 7) and 10(exp 3)/ s. Although the Mn addition to iron aluminide did not change the basic deformation characteristics, the Mn-modified alloy was slightly weaker. In the regime where deformation of FeAl occurs by a high stress exponent mechanism (n = 6), strength increases as the grain size decreases at least for diameters between approx. 200 and approx. 10 microns. Due to the limitation in the grain size-flow stress-temperature-strain rate database, the influence of further reductions of the grain size on strength is uncertain. Based on the appearance of subgrains in deformed iron aluminide, the comparison of grain diameters to expected subgrain sizes, and the grain size exponent and stress exponent calculated from deformation experiments, it is believed that grain size strengthening is the result of an artificial limitation on subgrain size as proposed by Sherby, Klundt and Miller.

Whittenberger, J. Daniel↗

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗

Predicting Unreinforced Fabric Mechanical Behavior with Recurrent Neural Networks

Unreinforced woven fabrics are widely employed in various high-performance applications, including parachute deployment systems, airbags, and ballistic armor. The analysis of such materials is inherently complex due to the multiscale structure of these materials, and the dependence of macroscale behavior on changes that occur at lower scales. Previously, NASA’s Multiscale Analysis Tool (NASMAT) showed its capability in predicting unreinforced fabric behavior at the macroscale by capturing finite rotations that occur at the mesoscale. Though effective, the tool can face high computational cost for large, complex problems, motivating the need for the development of a surrogate model that can capture the same behavior. A recurrent neural network (RNN) was developed and trained on virtual NASMAT data to mimic the physics-based solutions while improving the computational runtime. The architecture of the RNN to best simulate the fabric behavior was carefully crafted based on heuristic knowledge of predicting physics-based temporal data, manual hyperparameter case studies, and Hyperband optimization.. The resultant model was able to predict a variety of stress-strain curves for fabrics with different mesoscale geometries, and was further validated by comparing to experimental data for the K706 style Kevlar plain-weave fabric, demonstrating the ability of the model to effectively capture the geometric changes in the fabric without explicitly calculating them, as is done in NASMAT. Furthermore, the tool showed its ability to improve on the runtime by a factor of 10 for fabric solutions compared to the multiscale tool, which would further enable the simulation of complex loading scenarios on unreinforced fabrics.

Fabric↗

Nucleation and growth of rolling contact failure of 440C bearing steel

A 'two-body' elasto-plastic finite element model of 2-dimensional rolling and rolling-plus-sliding was developed to treat the effect of surface irregularities. The model consists of a smooth cylinder in contact with a semi-infinite half-space that is either smooth or fitted with one of 0.4 microns deep or 7 microns deep groove, or a 0.4 microns high ridge-like asperity. The model incorporates elastic-linear-kinematic hardening-plastic (ELKP) and non-linear-kinematic hardening-plastic (NLKP) material constitutive relations appropriate for hardened bearing steel and the 440C grade. The calculated contact pressure distribution is Hertzian for smooth body contact, and it displays intense, stationary, pressure spikes superposed on the Hertzian pressure for contact with the grooved and ridged surface. The results obtained for the 0.4 microns deep groove compare well with those reported by Elsharkawy and Hamrock for an EHD lubricated contact. The effect of translating the counterface on the half space as opposed to indenting the half space with the counter face with no translation is studied. The stress and strain values near the surface are found to be similar for the two cases, whereas they are significantly different in the subsurface. It is seen that when tiny shoulders are introduced at the edge of the groove in the finite element model, the incremental plasticity and residual stresses are significantly higher in the vicinity of the right shoulder (rolling direction is from left to right) than at the left shoulder. This may explain the experimental observation that the spall nucleation occurs at the exit end of the artificially planted indents. Pure rolling calculations are compared with rolling + sliding calculations. For a coefficient of friction, mu = 0.1, the effect of friction is found to be small. Efforts were made to identify the material constitutive relations which best describe the deformation characteristics of the bearing steels in the initial few cycles. Elastic-linear-kinematic hardening-plastic (ELKP) material constitutive relations produce less net plastic deformation in the initial stages for a given stress, than seen in experiments. A new set of constitutive relations: non-linear-kinematic hardening-plastic (NLKP) was used. This material model produces more plasticity than the ELKP model and shows promise for treating the net distortions in the early stages. Techniques for performing experimental measurements that can be compared with the finite element calculations were devised. The measurements are being performed on 9mm-diameter, 440C steel cylindrical rolling elements in contact with 12.5 mm-diameter, 52100 steel balls in a 3-ball-rod fatigue test machine operating at 3600 RPM. Artificial, 7 microns deep, indents were inserted on the running track of the cylindrical rolling elements and profilometer measurements of these indents made, before and after the rolling. These preliminary measurements show that the indents are substantially deformed plastically in the process of rolling. The deformations of the groove calculated with the finite element model are comparable to those measured experimentally.

Gupta, V.↗

Variation of refractive index in strained In(x)Ga(1-x)As-GaAs heterostructures

In(x)Ga(1-x)As-GaAs heterostructures and strained-layer superlattices can be used as optical waveguides. For such applications it is important to know explicitly the refractive index variation with mismatch strain and with alloying in the ternary layer. Starting from the Kramers-Kronig integral dispersion relations, a model has been developed from which the refractive index change in the ternary layer of In(x)Ga(1-x)As-GaAs heterojunctions can be calculated. The results are presented and discussed. The expected changes in a superlattice have been qualitatively predicted.

Das, U.↗

Manufacturing stresses and strains in filament wound cylinders

Tests were performed to verify a previously developed model for simulating the manufacturing process of filament wound cylinders. The axial and hoop strains were measured during cure inside a filament wound Fiberite T300/976 graphite-epoxy cylinder. The measured strains were compared to those computed by the model. Good agreements were found between the data and the model, indicating that the model is a useful representation of the process. For the conditions of the test, the manufacturing stresses inside the cylinder were also calculated using the model.

Calius, E. P.↗

Probabilistic nonlinear finite element analysis of composite structures

A probabilistic finite element analysis procedure for laminated composite shells is developed. A total Lagrangian finite element formulation, employing a degenerated three-dimensional laminated composite shell element with the full Green-Lagrange strains and first-order shear deformable kinematics, is used. The first-order second-moment technique for probabilistic finite element analysis of random fields is employed, and results are presented in the form of mean and variance of the structural response. Reliability calculations are made by using the first-order reliability method combined with sensitivity derivatives from the finite element analysis. Both ply-level and micromechanics-level random variables are incorporated, the latter by means of the Aboudi micromechanics model. Two sample problems are solved to verify the accuracy of the procedures developed and to quantify the variability of certain material type/structure combinations. In general, the procedure is quite effective in determining the response statistics and reliability for linear and geometric nonlinear behavior of laminated composite shells.

Engelstad, S. P.↗

Effect of mesh distortion on the accuracy of transverse shear stresses and their sensitivity coefficients in multilayered composites

A study is made of the effect of mesh distortion on the accuracy of transverse shear stresses and their first-order and second-order sensitivity coefficients in multilayered composite panels subjected to mechanical and thermal loads. The panels are discretized by using a two-field degenerate solid element, with the fundamental unknowns consisting of both displacement and strain components, and the displacement components having a linear variation throughout the thickness of the laminate. A two-step computational procedure is used for evaluating the transverse shear stresses. In the first step, the in-plane stresses in the different layers are calculated at the numerical quadrature points for each element. In the second step, the transverse shear stresses are evaluated by using piecewise integration, in the thickness direction, of the three-dimensional equilibrium equations. The same procedure is used for evaluating the sensitivity coefficients of transverse shear stresses. Numerical results are presented showing no noticeable degradation in the accuracy of the in-plane stresses and their sensitivity coefficients with mesh distortion. However, such degradation is observed for the transverse shear stresses and their sensitivity coefficients. The standard of comparison is taken to be the exact solution of the three-dimensional thermoelasticity equations of the panel.

Noor, Ahmed K.↗

Experimental Studies in Helicopter Vertical Climb Performance

Data and analysis from an experimental program to measure vertical climb performance on an eight-foot model rotor are presented. The rotor testing was performed using a unique moving-model facility capable of accurately simulating the flow conditions during axial flight, and was conducted from July 9, 1992 to July 16, 1992 at the Dynamic Model Track, or 'Long Track,' just prior to its demolition in August of 1992. Data collected during this brief test program included force and moment time histories from a sting-mounted strain gauge balance, support carriage velocity, and rotor rpm pulses. In addition, limited video footage (of marginal use) was recorded from smoke flow studies for both simulated vertical climb and descent trajectories. Analytical comparisons with these data include a series of progressively more detailed calculations ranging from simple momentum theory, a prescribed wake method, and a free-wake prediction.

McKillip, Robert M., Jr.↗

Effect of chemistry and turbulence on NO formation in oxygen-natural gas flames

The effects of chemistry and turbulence on NO formation in oxygen-natural turbulent diffusion flames gas flames have been investigated. The chemistry of nitric oxides has been studied numerically in the counterflow configuration. Systematic calculations with the GRI 2.11 mechanism for combustion of methane and NO chemistry were conducted to provide a base case. It was shown that the 'simple' Zeldovich mechanism accounts for more than 75% of N2 consumption in the flame in a range of strain-rates varying between 10 and 1000 s-l. The main shortcomings of this mechanism are: 1) overestimation (15%) of the NO production rate at low strain-rates because it does not capture the reburn due to the hydrocarbon chemistry, and 2) underestimation (25%) of the NO production rate at high strainrates because it ignores NO production through the prompt mechanism. Reburn through the Zeldovich mechanism alone proves to be significant at low strain-rates. A one-step model based on the Zeldovich mechanism and including reburn has been developed. It shows good agreement with the GRI mechanism at low strain-rates but underestimates significantly N2 consumption (about 50%) at high strain-rates. The role of turbulence has been assessed by using an existing 3-D DNS data base of a diffusion flame in decaying turbulence. Two PDF closure models used in practical industrial codes for turbulent NO formation have been tested. A simpler version of the global one-step chemical scheme for NO compared to that developed in this study was used to test the closure assumptions of the PDF models, because the data base could not provide all the necessary ingredients. Despite this simplification, it was possible to demonstrate that the current PDF models for NO overestimate significantly the NO production rate due to the fact that they neglect the correlations between the fluctuations in oxygen concentration and temperature. A single scalar PDF model for temperature that accounts for such correlations based on laminar flame considerations has been developed and showed excellent agreement with the values given by the DNS.

Samaniego, J. -M.↗

Effect of a hypergravity environment on cortical bone elasticity in rats

There is considerable interest in determining whether hypergravity can be used as a countermeasure for microgravity-induced bone loss. This study was conducted on 20 immature male rats in order to investigate possible elastic adaptations of cortical bone in rapidly growing rats exposed to chronic hypergravity. Ten rats were continuously centrifuged for 14 days at twice gravitational acceleration (2G) on a 12.75 foot radius centrifuge and 10 rats concurrently acted as stationary controls. The effect of hypergravity on the elastic characteristics of cortical bone was quantified via ultrasonic wave propagation. Propagation velocities of longitudinal and shear waves were measured through cubic cortical specimens from the posterior femoral diaphyses. Density was measured with an Archimedes' technique. The orthotropic elastic properties were calculated and used to compare the difference between groups. Results showed an average increase in both the Young's moduli (Eii, + 2.2%) and shear moduli (Gij, + 4.3%) with a statistically significant increase only in G12 (+15.7%, P = 0.046). The ratio of transverse to axial strain (Poisson's ratio, nuij) demonstrated statistically significant changes in nu12, nu21, nu13, and nu31 (P < 0.05). These findings suggest that although slight elastic changes were incurred via a hypergravity environment, the treatment level or duration in this study do not dramatically perturb the normal elastic behavior of cortical bone and that dramatic biomechanical differences noted in previous studies were due more to structural changes than material elasticity changes. Hypergravity applied post facto to a microgravity environment would offer further illucidation of this method as treatment for a degenerative spaceflight experience.

Non-NASA Center↗

Crystallized N-terminal domain of influenza virus matrix protein M1 and method of determining and using same

The matrix protein, M1, of influenza virus strain A/PR/8/34 has been purified from virions and crystallized. The crystals consist of a stable fragment (18 Kd) of the M1 protein. X-ray diffraction studies indicated that the crystals have a space group of P3.sub.t 21 or P3.sub.2 21. Vm calculations showed that there are two monomers in an asymmetric unit. A crystallized N-terminal domain of M1, wherein the N-terminal domain of M1 is crystallized such that the three dimensional structure of the crystallized N-terminal domain of M1 can be determined to a resolution of about 2.1 .ANG. or better, and wherein the three dimensional structure of the uncrystallized N-terminal domain of M1 cannot be determined to a resolution of about 2.1 .ANG. or better. A method of purifying M1 and a method of crystallizing M1. A method of using the three-dimensional crystal structure of M1 to screen for antiviral, influenza virus treating or preventing compounds. A method of using the three-dimensional crystal structure of M1 to screen for improved binding to or inhibition of influenza virus M1. The use of the three-dimensional crystal structure of the M1 protein of influenza virus in the manufacture of an inhibitor of influenza virus M1. The use of the three-dimensional crystal structure of the M1 protein of influenza virus in the screening of candidates for inhibition of influenza virus M1.

Luo, Ming↗

A kinetic model for the thermal nitridation of SiO2/Si

To explain the observed nitrogen distributions in thermally nitridated SiO2 films, a kinetic model is proposed in which the nitridation process is simulated, using the first-order chemical kinetics and Arrhenius dependence of the diffusion and reaction rates on temperature. The calculations show that initially, as the substrate reacts with diffusing nitrogen, a nitrogen-rich oxynitride forms at the SiO2-Si interface, while at nitridation temperatures above 1000 C, an oxygen-rich oxynitride subsequently forms at the interface, due to reaction of the substrate with an increasingly concentrated oxygen displaced by the slower nitridation of the SiO2. This sequence of events results in a nitrogen distribution in which the peak of the interfacial nitrogen concentration occurs away from the interface. The results are compared with the observed nitrogen distribution. The calculated results have correctly predicted the positions of the interfacial nitrogen peaks at the temperatures of 800, 1000, and 1150 C. To account for the observed width of the interfacial nitrogen distribution, it was found necessary to include in the simulations the effect of interfacial strain.

Vasquez, R. P.↗

Finite element elastic-plastic-creep and cyclic life analysis of a cowl lip

Results are presented of elastic, elastic-plastic, and elastic-plastic-creep analyses of a test-rig component of an actively cooled cowl lip. A cowl lip is part of the leading edge of an engine inlet of proposed hypersonic aircraft and is subject to severe thermal loadings and gradients during flight. Values of stresses calculated by elastic analysis are well above the yield strength of the cowl lip material. Such values are highly unrealistic, and thus elastic stress analyses are inappropriate. The inelastic (elastic-plastic and elastic-plastic-creep) analyses produce more reasonable and acceptable stress and strain distributions in the component. Finally, using the results from these analyses, predictions are made for the cyclic crack initiation life of a cowl lip. A comparison of predicted cyclic lives shows the cyclic life prediction from the elastic-plastic-creep analysis to be the lowest and, hence, most realistic.

Arya, Vinod K.↗

Finite element elastic-plastic-creep and cyclic life analysis of a cowl lip

Results are presented of elastic, elastic-plastic and elastic-plastic-creep analyses of a test-rig component of an actively cooled cowl lip. A cowl lip is part of the leading edge of an engine inlet of proposed hypersonic aircraft and is subject to severe thermal loadings and gradients during flight. Values of stresses calculated by elastic analysis are well above the yield strength of the cowl lip material. Such values are highly unrealistic, and thus elastic stress analyses are inappropriate. The inelastic (elastic-plastic and elastic-plastic-creep) analyses produce more reasonable and acceptable stress and strain distributions in the component. Finally, using the results from these analyses, predictions are made for the cyclic crack initiation life of a cowl lip. A comparison of predicted cyclic lives shows the cyclic life prediction from the elastic-plastic-creep analysis to be the lowest and, hence, most realistic.

Arya, V. K.↗

Mode 1 and Mode 2 Analysis of Graphite/Epoxy Composites Using Double Cantilever Beam and End-Notched Flexure Tests

The critical strain energy release rates associated with debonding of the adhesive bondlines in graphite/epoxy IM6/3501-6 interlaminar fracture specimens were investigated. Two panels were manufactured for this investigation; however, panel two was layed-up incorrectly. As a result, data collected from Panel Two serves no real purpose in this investigation. Double Cantilever Beam (DCB) specimens were used to determine the opening Mode I interlaminar fracture toughness, G1(sub c), of uni-directional fiber re-inforced composites. The five specimens tested from Panel One had an average value of 946.42J/sq m for G1(sub c) with an acceptable coefficient of variation. The critical strain energy release rate, G2(sub c), for initiation of delamination under inplane shear loading was investigated using the End-Notched Flexure (ENF) Test. Four specimens were tested from Panel One and an average value of 584.98J/sq m for G2(sub c) was calculated. Calculations from the DCB and ENF test results for Panel One represent typical values of G1(sub c) and G2(sub c) for the adhesive debonding in the material studied in this investigation.

Hufnagel, Kathleen P.↗

Optimal output fast feedback in two-time scale control of flexible arms

Control of lightweight flexible arms moving along predefined paths can be successfully synthesized on the basis of a two-time scale approach. A model following control can be designed for the reduced order slow subsystem. The fast subsystem is a linear system in which the slow variables act as parameters. The flexible fast variables which model the deflections of the arm along the trajectory can be sensed through strain gage measurements. For full state feedback design the derivatives of the deflections need to be estimated. The main contribution of this work is the design of an output feedback controller which includes a fixed order dynamic compensator, based on a recent convergent numerical algorithm for calculating LQ optimal gains. The design procedure is tested by means of simulation results for the one link flexible arm prototype in the laboratory.

Siciliano, B.↗

Development of a dual strain gage balance system for measuring light loads

A strain-gage-balance (SGB) force-measurement system is described which was designed to meet light-load requirements for an airfoil model tested in the NASA Langley low-turbulence pressure tunnel (LTPT). This system was developed to obtain direct force data needed to verify calculated pressure/force correlations used in previous aerodynamic tests. The three-component force-measurement system was designed so that the SGBs would simultaneously support and measure the following loads: weight of the airfoil (approximately 3.0 lbf); 8.0 lbf of lift; 0.5 lbf of drag; and 16.0 in. lbf of pitching moment. In addition to these design loads, the system was required to withstand 100-percent overload on all three components. The system comprises an airfoil, two SGBs, a thermal flexure, and a mounting plate. The installation of the system in the LTPT is also discussed.

Roberts, Paul W.↗