Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite element contact”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Finite Element Contact Analysis and Redesign of Rotorcraft Planetary Gear Sets With Asymmetric Profiles

This article explores the use of asymmetric gear teeth for rotorcraft main rotor planetary gearsets. Two baseline designs having symmetric gear teeth were selected for redesign studies. Finite element-based contact analyses were used to compare designs. Bending stress predictions for the two existing baseline designs were compared to strain gage measurements to confirm the suitability of the modeling. The same modeling method is then used to assess design changes to improve the bending stress and contact pressure magnitudes using asymmetric tooth profiles. For the first redesign study, the selected gearset is mildly asymmetric, with a 4° difference in pressure angles between the two sides of a tooth. For the second redesign study, the selected gearset is strongly asymmetric with a 18° difference in pressure angles. The results demonstrate that the asymmetric tooth approach provides additional design flexibility toward balancing the contact stresses and bending stresses of the sun, planet, and ring members.

gear↗
NASA NTRS1986

A constitutive law for finite element contact problems with unclassical friction

Techniques for modeling complex, unclassical contact-friction problems arising in solid and structural mechanics are discussed. A constitutive modeling concept is employed whereby analytic relations between increments of contact surface stress (i.e., traction) and contact surface deformation (i.e., relative displacement) are developed. Because of the incremental form of these relations, they are valid for arbitrary load-deformation histories. The motivation for the development of such a constitutive law is that more realistic friction idealizations can be implemented in finite element analysis software in a consistent, straightforward manner. Of particular interest is modeling of two-body (i.e., unlubricated) metal-metal, ceramic-ceramic, and metal-ceramic contact. Interfaces involving ceramics are of engineering importance and are being considered for advanced turbine engines in which higher temperature materials offer potential for higher engine fuel efficiency.

Plesha, M. E.↗
NASA NTRS1980

Finite elements for contact problems in two-dimensional elastodynamics

A finite element approach for contact problems in two dimensional elastodynamics was proposed. Sticking, sliding, and frictional contact were taken into account. The method consisted of a modification of the shape functions, in the contact region, in order to involve the nodes of the contacting body. The formulation was symmetric (both bodies were contactors and targets), in order to avoid interpenetration. Compatibility over the interfaces was satisfied. The method was applied to the impact of a block on a rigid target. It is shown that the formulation can be applied to fluid structure interaction, and to problems involving material nonlinearity.

Zimmermann, T. K.↗
NASA NTRS1989

Mixed finite element formulation for frictionless contact problems

Simple mixed finite element models and a computational procedure are presented for the solution of frictionless contact problems. The analytical formulation is based on a form of Reissner's large-rotation theory of the structure with the effect of transverse shear deformation included. The contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the internal forces (or stress resultants), the generalized displacements, and the Lagrange multipliers associated with the contact conditions. Numerical examples of axisymmetric deformations of a hemispherical shell and planar deformations of a circular ring are presented. Both structures are pressed against a rigid plate. Detailed information about the response of the structures is presented. The numerical studies demonstrate the high accuracy of the mixed models and the effectiveness of the computational procedure based on combining the geometrically nonlinear terms and the contact conditions in one iteration loop.

Noor, Ahmed K.↗
NASA NTRS1992

Contact stresses in meshing spur gear teeth: Use of an incremental finite element procedure

Contact stresses in meshing spur gear teeth are examined. The analysis is based upon an incremental finite element procedure that simultaneously determines the stresses in the contact region between the meshing teeth. The teeth themselves are modeled by two dimensional plain strain elements. Friction effects are included, with the friction forces assumed to obey Coulomb's law. The analysis assumes that the displacements are small and that the tooth materials are linearly elastic. The analysis procedure is validated by comparing its results with those for the classical two contacting semicylinders obtained from the Hertz method. Agreement is excellent.

Hsieh, Chih-Ming↗
NASA NTRS1976

On the explicit finite element formulation of the dynamic contact problem of hyperelastic membranes

Contact-impact problems involving finite deformation axisymmetric membranes are solved by the finite element method with explicit time integration. The formulation of the membrane element and the contact constraint conditions are discussed. The hyperelastic, compressible Blatz and Ko material is used to model the material properties of the membrane. Two example problems are presented.

Hallquist, J. O.↗
NASA NTRS2003

Simulating Fatigue Crack Growth in Spiral Bevel Pinion

This project investigates computational modeling of fatigue crack growth in spiral bevel gears. Current work is a continuation of the previous efforts made to use the Boundary Element Method (BEM) to simulate tooth-bending fatigue failure in spiral bevel gears. This report summarizes new results predicting crack trajectory and fatigue life for a spiral bevel pinion using the Finite Element Method (FEM). Predicting crack trajectories is important in determining the failure mode of a gear. Cracks propagating through the rim may result in catastrophic failure, whereas the gear may remain intact if one tooth fails and this may allow for early detection of failure. Being able to predict crack trajectories is insightful for the designer. However, predicting growth of three-dimensional arbitrary cracks is complicated due to the difficulty of creating three-dimensional models, the computing power required, and absence of closed- form solutions of the problem. Another focus of this project was performing three-dimensional contact analysis of a spiral bevel gear set incorporating cracks. These analyses were significant in determining the influence of change of tooth flexibility due to crack growth on the magnitude and location of contact loads. This is an important concern since change in contact loads might lead to differences in SIFs and therefore result in alteration of the crack trajectory. Contact analyses performed in this report showed the expected trend of decreasing tooth loads carried by the cracked tooth with increasing crack length. Decrease in tooth loads lead to differences between SIFs extracted from finite element contact analysis and finite element analysis with Hertz contact loads. This effect became more pronounced as the crack grew.

Ural, Ani↗
NASA NTRS1984

Finite element analysis of steadily moving contact fields

By introducing a moving updated Lagrangian observer, this paper develops traveling finite elements with the capacity to handle the global response resulting from steadily moving contact fields. The generality of the results is such that large deformation kinematics and kinetics as well as the full compliment of inertial fields can be handled. To streamline the handling of nonlinear behavior, an elliptically constrained solution algorithm is also developed. Employing this algorithm, the results of several numerical benchmarking studies are presented which illustrate the capacity of the moving updated Lagrangian formulation as well as the potential effects of nonlinearity.

Padovan, J.↗
NASA NTRS2004

Subsurface Stress Fields in FCC Single Crystal Anisotropic Contacts

Single crystal superalloy turbine blades used in high pressure turbomachinery are subject to conditions of high temperature, triaxial steady and alternating stresses, fretting stresses in the blade attachment and damper contact locations, and exposure to high-pressure hydrogen. The blades are also subjected to extreme variations in temperature during start-up and shutdown transients. The most prevalent high cycle fatigue (HCF) failure modes observed in these blades during operation include crystallographic crack initiation/propagation on octahedral planes, and non-crystallographic initiation with crystallographic growth. Numerous cases of crack initiation and crack propagation at the blade leading edge tip, blade attachment regions, and damper contact locations have been documented. Understanding crack initiation/propagation under mixed-mode loading conditions is critical for establishing a systematic procedure for evaluating HCF life of single crystal turbine blades. This paper presents analytical and numerical techniques for evaluating two and three dimensional subsurface stress fields in anisotropic contacts. The subsurface stress results are required for evaluating contact fatigue life at damper contacts and dovetail attachment regions in single crystal nickel-base superalloy turbine blades. An analytical procedure is presented for evaluating the subsurface stresses in the elastic half-space, based on the adaptation of a stress function method outlined by Lekhnitskii. Numerical results are presented for cylindrical and spherical anisotropic contacts, using finite element analysis (FEA). Effects of crystal orientation on stress response and fatigue life are examined. Obtaining accurate subsurface stress results for anisotropic single crystal contact problems require extremely refined three-dimensional (3-D) finite element grids, especially in the edge of contact region. Obtaining resolved shear stresses (RSS) on the principal slip planes also involves considerable post-processing work. For these reasons it is very advantageous to develop analytical solution schemes for subsurface stresses, whenever possible.

Arakere, Nagaraj K.↗
NASA NTRS1989

Elasto-plastic finite element analyses of two-dimensional rolling and sliding contact deformation of bearing steel

An elastoplastic finite element analysis of repeated rolling-plus-sliding contact is presented for the elastic-linear-kinematic hardening-plastic properties of hardened bearing steel at room temperature. Contacts in the absence of frictional heating are found to produce peak cyclic strain ranges and residual stresses that are similar to the values obtained for pure rolling. For the cases considered, the energy dissipated as heat by plastic deformation is found to be 3-15 percent of the energy dissipated at the surface by friction.

Kumar, A. M.↗
NASA NTRS2002

Design and Stress Analysis of Low-Noise Adjusted Bearing Contact Spiral Bevel Gears

An integrated computerized approach for design and stress analysis of low-noise spiral bevel gear drives with adjusted bearing contact has been developed. The computation procedure is an iterative process, requiring four separate steps that provide: (a) a parabolic function of transmission errors that is able to reduce the effect of errors of alignment, and (b) reduction of the shift of bearing contact caused by misalignment. Application of finite element analysis permits the contact and bending stresses to be determined and investigate the formation of the bearing contact. The design of finite element models and boundary conditions is automated and does not require an intermediate CAD computer program. A commercially available finite element analysis computer program with contact capability was used to conduct the stress analysis. The theory developed is illustrated with numerical examples.

Litvin, Faydor L.↗
NASA NTRS1994

Comparison of Gap Elements and Contact Algorithm for 3D Contact Analysis of Spiral Bevel Gears

Three dimensional stress analysis of spiral bevel gears in mesh using the finite element method is presented. A finite element model is generated by solving equations that identify tooth surface coordinates. Contact is simulated by the automatic generation of nonpenetration constraints. This method is compared to a finite element contact analysis conducted with gap elements.

Bibel, G. D.↗
NASA NTRS2000

Calcaneal loading during walking and running

PURPOSE: This study of the foot uses experimentally measured kinematic and kinetic data with a numerical model to evaluate in vivo calcaneal stresses during walking and running. METHODS: External ground reaction forces (GRF) and kinematic data were measured during walking and running using cineradiography and force plate measurements. A contact-coupled finite element model of the foot was developed to assess the forces acting on the calcaneus during gait. RESULTS: We found that the calculated force-time profiles of the joint contact, ligament, and Achilles tendon forces varied with the time-history curve of the moment about the ankle joint. The model predicted peak talocalcaneal and calcaneocuboid joint loads of 5.4 and 4.2 body weights (BW) during walking and 11.1 and 7.9 BW during running. The maximum predicted Achilles tendon forces were 3.9 and 7.7 BW for walking and running. CONCLUSIONS: Large magnitude forces and calcaneal stresses are generated late in the stance phase, with maximum loads occurring at approximately 70% of the stance phase during walking and at approximately 60% of the stance phase during running, for the gait velocities analyzed. The trajectories of the principal stresses, during both walking and running, corresponded to each other and qualitatively to the calcaneal trabecular architecture.

NASA Center ARC↗
NASA NTRS2003

Finite Element Modeling of a Cylindrical Contact Using Hertzian Assumptions

The turbine blades in the high-pressure fuel turbopump/alternate turbopump (HPFTP/AT) are subjected to hot gases rapidly flowing around them. This flow excites vibrations in the blades. Naturally, one has to worry about resonance, so a damping device was added to dissipate some energy from the system. The foundation is now laid for a very complex problem. The damper is in contact with the blade, so now there are contact stresses (both normal and tangential) to contend with. Since these stresses can be very high, it is not all that difficult to yield the material. Friction is another non-linearity and the blade is made out of a Nickel-based single-crystal superalloy that is orthotropic. A few approaches exist to solve such a problem and computer models, using contact elements, have been built with friction, plasticity, etc. These models are quite cumbersome and require many hours to solve just one load case and material orientation. A simpler approach is required. Ideally, the model should be simplified so the analysis can be conducted faster. When working with contact problems determining the contact patch and the stresses in the material are the main concerns. Closed-form solutions for non-conforming bodies, developed by Hertz, made out of isotropic materials are readily available. More involved solutions for 3-D cases using different materials are also available. The question is this: can Hertzian1 solutions be applied, or superimposed, to more complicated problems-like those involving anisotropic materials? That is the point of the investigation here. If these results agree with the more complicated computer models, then the analytical solutions can be used in lieu of the numerical solutions that take a very long time to process. As time goes on, the analytical solution will eventually have to include things like friction and plasticity. The models in this report use no contact elements and are essentially an applied load problem using Hertzian assumptions to determine the contact patch dimensions.

Knudsen, Erik↗
NASA NTRS2004

Subsurface Stress Fields In Single Crystal (Anisotropic) Contacts

Single crystal superalloy turbine blades used in high pressure turbomachinery are subject to conditions of high temperature, triaxial steady and alternating stresses, fretting stresses in the blade attachment and damper contact locations, and exposure to high-pressure hydrogen. The blades are also subjected to extreme variations in temperature during start-up and shutdown transients. The most prevalent HCF failure modes observed in these blades during operation include crystallographic crack initiation/propagation on octahedral planes, and noncrystallographic initiation with crystallographic growth. Numerous cases of crack initiation and crack propagation at the blade leading edge tip, blade attachment regions, and damper contact locations have been documented. Understanding crack initiation/propagation under mixed-mode loading conditions is critical for establishing a systematic procedure for evaluating HCF life of single crystal turbine blades. This paper presents analytical and numerical techniques for evaluating two and three dimensional subsurface stress fields in anisotropic contacts. The subsurface stress results are required for evaluating contact fatigue life at damper contacts and dovetail attachment regions in single crystal nickel-base superalloy turbine blades. An analytical procedure is , presented, for evaluating the subsurface stresses in the elastic half-space, using a complex potential method outlined by Lekhnitskii. Numerical results are presented for cylindrical and spherical anisotropic contacts, using finite element analysis. Effects of crystal orientation on stress response and fatigue life are examined.

Arakere, Nagaraj K.↗
NASA NTRS1991

Elasto-plastic coupled temperature-displacement finite element analysis of two-dimensional rolling-sliding contact with a translating heat source

The present paper describes a transient translating elastoplastic thermomechanical finite element model to study two-dimensional frictional rolling contact. Frictional two-dimensional contact is simulated by repeatedly translating a nonuniform thermomechanical distribution across the surface of an elastoplastic half space. The half space is represented by a two-dimensional finite element mesh with appropriate boundaries. Calculations are for an elastic-perfectly plastic material and the selected thermophysical properties are assumed to be temperature independent. The paper presents temperature variations, stress and plastic strain distributions and deformations. Residual tensile stresses are observed. The magnitude and depth of these stresses depends on (1) the temperature gradients and (2) the magnitudes of the normal and tangential tractions.

Rubin, C. A.↗
NASA NTRS2002

Procedure for Tooth Contact Analysis of a Face Gear Meshing With a Spur Gear Using Finite Element Analysis

A procedure was developed to perform tooth contact analysis between a face gear meshing with a spur pinion using finite element analysis. The face gear surface points from a previous analysis were used to create a connected tooth solid model without gaps or overlaps. The face gear surface points were used to create a five tooth face gear Patran model (with rim) using Patran PCL commands. These commands were saved in a series of session files suitable for Patran input. A four tooth spur gear that meshes with the face gear was designed and constructed with Patran PCL commands. These commands were also saved in a session files suitable for Patran input. The orientation of the spur gear required for meshing with the face gear was determined. The required rotations and translations are described and built into the session file for the spur gear. The Abaqus commands for three-dimensional meshing were determined and verified for a simplified model containing one spur tooth and one face gear tooth. The boundary conditions, loads, and weak spring constraints were determined to make the simplified model work. The load steps and load increments to establish contact and obtain a realistic load was determined for the simplified two tooth model. Contact patterns give some insight into required mesh density. Building the two gears in two different local coordinate systems and rotating the local coordinate systems was verified as an easy way to roll the gearset through mesh. Due to limitation of swap space, disk space and time constraints of the summer period, the larger model was not completed.

Bibel, George↗