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 829 records · Page 46

Analysis of crack propagation in roller bearings using the boundary integral equation method - A mixed-mode loading problem

Crack propagation in a rotating inner raceway of a high-speed roller bearing is analyzed using the boundary integral method. The model consists of an edge plate under plane strain condition upon which varying Hertzian stress fields are superimposed. A multidomain boundary integral equation using quadratic elements was written to determine the stress intensity factors KI and KII at the crack tip for various roller positions. The multidomain formulation allows the two faces of the crack to be modeled in two different subregions, making it possible to analyze crack closure when the roller is positioned on or close to the crack line. KI and KII stress intensity factors along any direction were computed. These calculations permit determination of crack growth direction along which the average KI times the alternating KI is maximum.

Ghosn, L. J.↗

Global crustal response model

A computer program was developed for calculating the radial displacement due to the body tide as a function of spatial position and time. The positions of the moon and sun were evaluated by means of the Hill-Brown and Newcomb theories, respectively. Green's functions were utilized to simulate the elastic Earth responses due to point loads, a convolution program has been developed to evaluate the elastic Earth responses due to ocean loading by a given ocean tide constituent. These responses include radial and horizontal displacement, tilt and gravity acceleration, and strain tensor components.

Source record↗

Specimen size and geometry effects on fracture toughness of Al2O3 measured with short rod and short bar chevron-notch specimens

Plane strain fracture toughness measurements were made on Al2O3 using short rod and short bar chevron notch specimens previously calibrated by the authors for their dimensionless stress intensity factor coefficients. The measured toughness varied systematically with variations in specimen size, proportions, and chevron notch angle apparently due to their influence on the amount of crack extension to maximum load (the measurement point). The toughness variations are explained in terms of a suspected rising R curve for the material tested, along with a discussion of an unavoidable imprecision in the calculation of K sub Ic for materials with rising R curves when tested with chevron notch specimens.

Shannon, J. L., Jr.↗

A time-dependent ice sheet model - Preliminary results

A numerical model of ice sheet flow is developed, and preliminary results are described. This model includes vertical resolution of temperature, stress, and strain rate which represents a considerable improvement over previous vertically averaged ice sheet models. The model follows the flow of ice along a flow line within an ice sheet drainage basin. Longitudinal stresses and basal sliding are included. Basal sliding is dependent on the base shear stress and a specified distribution of basal water pressure. The numerical methods used to solve the coupled set of stress and velocity equations for the static and time-evolutionary cases are discussed. A steady state profile simulating an ice stream is calculated for a particular set of input parameters, and changes in the profile are examined for different choices of parameters. Preliminary studies of response behavior are completed using a simplified ice sheet geometry with a fixed terminus or grounding line. The results of these studies illustrate ice sheet thinning in response to a lowered sea level or to a reduction in the extent of ice rises (or pinning points) within ice shelves.

Bindschadler, R. A.↗

A real time neural net estimator of fatigue life

A neural net architecture is proposed to estimate, in real-time, the fatigue life of mechanical components, as part of the Intelligent Control System for Reusable Rocket Engines. Arbitrary component loading values were used as input to train a two hidden-layer feedforward neural net to estimate component fatigue damage. The ability of the net to learn, based on a local strain approach, the mapping between load sequence and fatigue damage has been demonstrated for a uniaxial specimen. Because of its demonstrated performance, the neural computation may be extended to complex cases where the loads are biaxial or triaxial, and the geometry of the component is complex (e.g., turbopump blades). The generality of the approach is such that load/damage mappings can be directly extracted from experimental data without requiring any knowledge of the stress/strain profile of the component. In addition, the parallel network architecture allows real-time life calculations even for high frequency vibrations. Owing to its distributed nature, the neural implementation will be robust and reliable, enabling its use in hostile environments such as rocket engines. This neural net estimator of fatigue life is seen as the enabling technology to achieve component life prognosis, and therefore would be an important part of life extending control for reusable rocket engines.

Troudet, T.↗

A real time neural net estimator of fatigue life

A neural network architecture is proposed to estimate, in real-time, the fatigue life of mechanical components, as part of the intelligent Control System for Reusable Rocket Engines. Arbitrary component loading values were used as input to train a two hidden-layer feedforward neural net to estimate component fatigue damage. The ability of the net to learn, based on a local strain approach, the mapping between load sequence and fatigue damage has been demonstrated for a uniaxial specimen. Because of its demonstrated performance, the neural computation may be extended to complex cases where the loads are biaxial or triaxial, and the geometry of the component is complex (e.g., turbopumps blades). The generality of the approach is such that load/damage mappings can be directly extracted from experimental data without requiring any knowledge of the stress/strain profile of the component. In addition, the parallel network architecture allows real-time life calculations even for high-frequency vibrations. Owing to its distributed nature, the neural implementation will be robust and reliable, enabling its use in hostile environments such as rocket engines.

Troudet, T.↗

Strain-energy release rate analysis of plates with postbuckled delaminations

The objectives are to present the derivation of the new virtual crack closure technique, evaluate the accuracy of the technique, and finally to present the results of a limited parametric study of laminates with a postbuckled delamination. Although the new virtual crack closure technique is general, only homogeneous, isotropic laminates were analyzed. This was to eliminate the variation of flexural stiffness with orientation, which occurs even for quasi-isotropic laminates. This made it easier to identify the effect of geometrical parameters on G. The new virtual crack closure technique is derived. Then the specimen configurations are described. Next, the stress analyses is discussed. Finally, the virtual crack closure technique is evaluated and then used to calculate the distribution of G along the delamination front of several laminates with a postbuckled delamination.

Whitcomb, John D.↗

On the interactions between large-scale structure and fine-grained turbulence in a free shear flow. III - A numerical solution

The results of a numerical computation of the interactions between the horizontally periodic monochromatic component of a large-scale coherent structure and the fine-grained turbulence in a mixing layer are presented. In the numerical calculations, the appropriate dependent variable is one which comprises both the mean and the large-scale coherent structure. The dynamical equations obtained for such a total coherent structure quantity are identical to the unsteady equations for the mean quantities in the Reynolds sense, except that the fine-grained turbulent stresses are interpreted as being conditionally averaged. It is shown how three dimensional fine-grained turbulence can be produced indirectly from the two dimensional large-scale structure via the isotropizing process of the approximated pressure-strain correlation.

Gatski, T. B.↗

Higher-order derivative correlations and the alignment of small-scale structures in isotropic numerical turbulence

In a three dimensional simulation higher order derivative correlations, including skewness and flatness factors, are calculated for velocity and passive scalar fields and are compared with structures in the flow. The equations are forced to maintain steady state turbulence and collect statistics. It is found that the scalar derivative flatness increases much faster with Reynolds number than the velocity derivative flatness, and the velocity and mixed derivative skewness do not increase with Reynolds number. Separate exponents are found for the various fourth order velocity derivative correlations, with the vorticity flatness exponent the largest. Three dimensional graphics show strong alignment between the vorticity, rate of strain, and scalar-gradient fields. The vorticity is concentrated in tubes with the scalar gradient and the largest principal rate of strain aligned perpendicular to the tubes. Velocity spectra, in Kolmogorov variables, collapse to a single curve and a short minus 5/3 spectral regime is observed.

Kerr, R. A.↗

A simplified computer solution for the flexibility matrix of contacting teeth for spiral bevel gears

A computer code, FLEXM, was developed to calculate the flexibility matrices of contacting teeth for spiral bevel gears using a simplified analysis based on the elementary beam theory for the deformation of gear and shaft. The simplified theory requires a computer time at least one order of magnitude less than that needed for the complete finite element method analysis reported earlier by H. Chao, and it is much easier to apply for different gear and shaft geometries. Results were obtained for a set of spiral bevel gears. The teeth deflections due to torsion, bending moment, shearing strain and axial force were found to be in the order 10(-5), 10(-6), 10(-7), and 10(-8) respectively. Thus, the torsional deformation was the most predominant factor. In the analysis of dynamic load, response frequencies were found to be larger when the mass or moment of inertia was smaller or the stiffness was larger. The change in damping coefficient had little influence on the resonance frequency, but has a marked influence on the dynamic load at the resonant frequencies.

Hsu, C. Y.↗

GPS measurements of strain accumulation across the Imperial Valley, California: 1986-1989

The Global Positioning System (GPS) data collected in southern California from 1986 to 1989 indicate considerable strain accumulation across the Imperial Valley. Displacements are computed at 29 stations in and near the valley from 1986 to 1988, and at 11 sites from 1988 to 1989. The earlier measurements indicate 5.9 +/- 1.0 cm/yr right-lateral differential velocity across the valley, although the data are heavily influenced by the 1987 Superstition Hills earthquake sequence. Some measurements, especially the east-trending displacements, are suspects for large errors. The 1988 to 1989 GPS displacements are best modeled by 5.2 +/- 0.9 cm/yr of valley crossing deformation, but rates calculated from conventional geodetic measurements (3.4 to 4.3 cm/yr) fit the data nearly as well. There is evidence from GPS and Very Long Base Interferometry (VLBI) observations that the present slip rate along the southern San Andreas fault is smaller than the long-term geologic estimate, suggesting a lower earthquake potential than is currently assumed. Correspondingly, a higher earthquake potential is indicated for the San Jacinto fault. The Imperial Valley GPS sites form part of a 183 station network in southern California and northern Baja California, which spans a cross-section of the North American-Pacific plate boundary.

Larsen, Shawn↗

Global Positioning System measurements of strain accumulation across the Imperial Valley, California - 1986-1989

The Global Positioning System (GPS) data collected in southern California from 1986 to 1989 indicate considerable strain accumulation across the Imperial Valley. Displacements are computed at 29 stations in and near the valley from 1986 to 1988, and at 11 sites from 1988 to 1989. The earlier measurements indicate 5.9 =/- 1.0 cm/yr right-lateral differential velocity across the valley, although the data are heavily influenced by the 1987 Superstition Hills earthquake sequence. Some measurements, especially the east-trending displacements, are suspects for large errors. The 1988 to 1989 GPS displacements are best modeled by 5.2 =/- 0.9 cm/yr of valley crossing deformation, but rates calculated from conventional geodetic measurements (3.4 to 4.3 cm/yr) fit the data nearly as well. There is evidence from GPS and Very Long Base Interferometry (VLBI) observations that the present slip rate along the southern San Andreas fault is smaller than the long-term geologic estimate, suggesting a lower earthquake potential than is currently assumed. Correspondingly, a higher earthquake potential is indicated for the San Jacinto fault. The Imperial Valley GPS sites form part of a 183 station network in southern California and northern Baja California, which spans a cross-section of the North American-Pacific plate boundary.

Larsen, Shawn↗

Strain-energy release rate analysis of a laminate with a postbuckled delamination

The objectives are to present the derivation of the new virtual crack closure technique, evaluate the accuracy of the technique, and finally to present the results of a limited parametric study of laminates with a postbuckled delamination. Although the new virtual crack closure technique is general, only homogeneous, isotropic laminates were analyzed. This was to eliminate the variation of flexural stiffness with orientation, which occurs even for quasi-isotropic laminates. This made it easier to identify the effect of geometrical parameters on G. The new virtual crack closure technique is derived. Then the specimen configurations are described. Next, the stress analyses is discussed. Finally, the virtual crack closure technique is evaluated and then used to calculate the distribution of G along the delamination front of several laminates with a postbuckled delamination.

Whitcomb, John D.↗

The production of turbulent stress in a shear flow by irrotational fluctuations

This paper examines, both theoretically and experimentally, the effect produced by irrotational fluctuations, associated with a nearby turbulent field, in a region where the turbulence is initially very low but where there is a mean shear. Calculations are based on rapid distortion theory and experiments use linearized hot wire anemometers in an open circuit wind tunnel. Turbulent shear stress is observed to grow from zero to significant values in the interaction region. The magnitude and extent of this observed shear stress agree reasonably well with predictions of the analysis, when intermittency effects are included. It is concluded that turbulent stresses can be produced by irrotational fluctuations in a region of mean shear and that this effect can be estimated using rapid distortion theory if the overall strain ratio is not large.

Gartshore, I. S.↗

Finite Element Prediction of Loss Factors for Structures with Frequency-dependent Damping Treatments

A finite element procedure is described for calculating the loss factors for elastic structures to which frequency-dependent viscoelastic damping treatments were applied. The frequency dependence of the viscoelastic damping material is treated by approximating its shear modulus with a second-order polynomial so that the stiffnesses associated with the constant, linear, and quadratic terms can be combined, respectively, with the stiffness, damping, and mass matrices assembled for the rest of the structure. A single complex eigenvalue analysis is then performed in which the eigenvalues are purely imaginary. The loss factor is computed by the modal strain energy (MSE) approach. In the the MSE approach, the loss factor of a composite structure vibrating in one of its natural modes may be visualized as a weighted average of the loss factors of the component parts, with the relative stored energies as weighting constants. The finite element procedure, which can treat very general geometries, is illustrated for the case of a vibrating constrained-layer damped plate.

Everstine, G. C.↗

Acoustic fatigue life prediction for nonlinear structures with multiple resonant modes

An efficient method is presented for estimating the high-cycle fatigue lives of complex nonlinear structures. The procedure is based on an application of the method of equivalent linearization to construct a time-domain representation of the stress or strain response of a nonlinear structure. Fatigue estimates are obtained by processing the time-domain signal, and it is a simple matter to incorporate any damage accumulation model desired in the procedure. The method may also be applied to a wide range of complex nonlinear systems. Comparisons are presented of estimated fatigue lives for a nonlinear beam with random excitation obtained using the approximate method and using a conventional numerical simulation. The two techniques are found to give nearly identical results regardless of the level of nonlinearity. The number of calculations required in the approximate method is nearly the same as in a linear system and is not significantly affected by the type and level of nonlinearity in the structure.

Sun, J. Q.↗

Characterization of the tip field of a discrete dislocation pileup for the development of physically based micromechanics

It is shown, on the basis of calculations by Eshelby et al. (1951), Armstrong et al. (1966), and Chou and Li (1969), that a single parameter, such as the force on the leading dislocation (F), the crack extension force, or the stress intensity factor, is capable of characterizing uniquely the entire tip field of a discrete dislocation pileup, including the positions of mobile dislocations behind the locked leading dislocation at the tip. Conversely, the position of the i-th mobile dislocation X(i) is related to the value of F and is capable of characterizing the entire stress, strain, and displacement fields at the tip of a discrete dislocation pileup. If the interactions between dislocations are linear elastic, the measured positions of the mobile dislocations can be used to determine the value of F, which can then be used as a quantitative measure of the strength of a dislocation barrier resisting the propagation of a microslip or the nucleation of a microfracture.

Gao, Q.↗

Development of procedures for calculating stiffness and damping properties of elastomers in engineering applications. Part 4: Testing of elastomers under a rotating load

A test rig designed to measure stiffness and damping of elastomer cartridges under a rotating load excitation is described. The test rig employs rotating unbalance in a rotor which runs to 60,000 RPM as the excitation mechanism. A variable resonant mass is supported on elastomer elements and the dynamic characteristics are determined from measurements of input and output acceleration. Five different cartridges are considered: three of these are buttons cartridges having buttons located in pairs, with 120 between each pair. Two of the cartridges consist of 360 elastomer rings with rectangular cross-sections. Dynamic stiffness and damping are measured for each cartridge and compared with predictions at different frequencies and different strains.

Darlow, M. S.↗