Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A Method for Thermal Analysis of Spiral Bevel Gears

A modeling method for analyzing the three-dimensional thermal behavior of spiral bevel gears has been developed. The model surfaces are generated through application of differential geometry to the manufacturing process for face-milled spiral bevel gears. Contact on the gear surface is found by combining tooth contact analysis with three-dimensional Hertzian theory. The tooth contact analysis provides the principle curvatures and orientations of the two surfaces. This information is then used directly in the Hertzian analysis to find the contact size and maximum pressure. Heat generation during meshing is determined as a function of the applied load, sliding velocity, and coefficient of friction. Each of these factors change as the point of contact changes during meshing. A nonlinear finite element program was used to conduct the heat transfer analysis. This program permitted the time- and position-varying boundary conditions, found in operation, to be applied to a one-tooth model. An example model and analytical results are presented.

Handschuh, Robert F.↗
NASA NTRS1980

Contact law and impact responses of laminated composites

Static identation tests were performed to determine the law of contact between a steel ball and glass/epoxy and graphite/epoxy laminated composites. For both composites the power law with an index of 1.5 was found to be adequate for the loading curve. Substantial permanent deformations were noted after the unloading. A high order beam finite element was used to compute the dynamic contact force and response of the laminated composite subjected to the impact of an elastic sphere. This program can be used with either the classical Hertzian contact law or the measured contact law. A simple method is introduced for estimating the contact force and contact duration in elastic impacts.

Sun, C. T.↗
DOE OSTI2024

Fabrication and Performance Evaluation of Double-Sided Copper Nanowire Arrays as Thermal and Electrical Interfacial Layers

Reducing contact interface thermal and electrical resistances is in great demand across various industries, particularly in the semiconductor industry. This study introduces an approach using double-sided copper nanowire (Cu NW) arrays on copper sheets as both thermal and electrical interfacial layers, designed to effectively accommodate the topographical inconsistencies between contact surfaces. Further, experimental outcomes reveal a significant reduction in thermal contact resistance (TCR), with a value of 2.5 mm 2 K W –1 , thereby exceeding the efficiency of reported nanostructural thermal interface materials (TIMs). Additionally, when utilized as an electrical interfacial layer, these double-sided Cu NWs arrays dramatically reduced electrical contact resistance (ECR), outperforming traditional conductive grease in applications necessitating separable bonding, though showing comparable performance to costly silver-based conductive epoxies required for permanent, inseparable bonds. The promising results of the double-sided Cu NWs arrays in reducing both TCR and ECR, confirmed by finite element simulation, highlight their substantial potential in advancing TIMs and electrical interconnection applications across various sectors.

42 ENGINEERING↗
DOE OSTI2022

Stress Profile in Coating Layers of TRISO Fuel Particles in Contact with One Another

This work presents a discussion on a series of finite element analyses that assess stress evolution in the coating layers of tristructural isotropic (TRISO) particles in contact with each other while embedded in a matrix. The initial simulations were of applied uniaxial pressure versus matrix elastic modulus. These simulations predicted increasing stress in the silicon carbide coating layers of the TRISO particles with decreasing matrix elastic modulus. The second set of simulations focused on the effects of heating and cooling and the associated dimensional change on the state of stress in the coating layers. The general finding was that there was no significant difference below the coating layer’s deposition temperature. Although, above the deposition temperature, the contacting particles had higher stress compared with those that were separated. The third set of simulations focused on the effects of irradiation, specifically, creep, dimensional change, and swelling. An interface debonding model was introduced since these potential effects have a significant bearing on predicted stresses.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2000

Fretting Stresses in Single Crystal Superalloy Turbine Blade Attachments

Single crystal nickel base superalloy turbine blades are being utilized in rocket engine turbopumps and turbine engines because of their superior creep, stress rupture, melt resistance and thermomechanical fatigue capabilities over polycrystalline alloys. Currently the most widely used single crystal nickel base turbine blade superalloys are PWA 1480/1493 and PWA 1484. These alloys play an important role in commercial, military and space propulsion systems. High Cycle Fatigue (HCF) induced failures in aircraft gas turbine and rocket engine turbopump blades is a pervasive problem. Blade attachment regions are prone to fretting fatigue failures. Single crystal nickel base superalloy turbine blades are especially prone to fretting damage because the subsurface shear stresses induced by fretting action at the attachment regions can result in crystallographic initiation and crack growth along octahedral planes. Furthermore, crystallographic crack growth on octahedral planes under fretting induced mixed mode loading can be an order of magnitude faster than under pure mode I loading. This paper presents contact stress evaluation in the attachment region for single crystal turbine blades used in the NASA alternate Advanced High Pressure Fuel Turbo Pump (HPFTP/AT) for the Space Shuttle Main Engine (SSME). Single crystal materials have highly orthotropic properties making the position of the crystal lattice relative to the part geometry a significant factor in the overall analysis. Blades and the attachment region are modeled using a large-scale 3D finite element (FE) model capable of accounting for contact friction, material orthotrophy, and variation in primary and secondary crystal orientation. Contact stress analysis in the blade attachment regions is presented as a function of coefficient of friction and primary and secondary crystal orientation, Stress results are used to discuss fretting fatigue failure analysis of SSME blades. Attachment stresses are seen to reach peak values at locations where fretting cracks have been observed. Fretting stresses at the attachment region are seen to vary significantly as a function of crystal orientation. Attempts to adapt techniques used for estimating fatigue life in the airfoil region, for life calculations in the attachment region, are presented. An effective model for predicting crystallographic crack initiation under mixed mode loading is required for life prediction under fretting action.

Arakere, Nagaraj K.↗
NASA NTRS1987

Mechanisms of rolling contact spalling

The results of a study aimed at analyzing the mechanical material interactions responsible for rolling contact spalling of the 440 C steel, high pressure oxygen turbopump bearings are presented. A coupled temperature displacement finite element analysis of the effects of friction heating under the contact is presented. The contact is modelled as a stationary, heat generating, 2 dimensional indent in an elastic perfectly plastic half-space with heat fluxes up to 8.6 x 10000 KW/m sq comparable to those generated in the bearing. Local temperatures in excess of 1000 C are treated. The calculations reveal high levels of residual tension after the contact is unloaded and cools. Efforts to promote Mode 2/Mode 3 fatigue crack growth under cyclic torsion in hardened 440 C steel are described. Spalls produced on 440 C steel by a 3 ball/rod rolling contact testing machine were studied with scanning microscopy. The shapes of the cyclic, stress strain hysteresis loops displayed by hardened 440 C steel in cyclic torsion at room temperature are defined for the plastic strain amplitudes encountered in rolling/sliding contact. Results of these analyses are discussed in detail.

Kumar, A. M.↗
NASA NTRS1994

Modeling Thermal Contact Resistance

One difficulty in using cryocoolers is making good thermal contact between the cooler and the instrument being cooled. The connection is often made through a bolted joint. The temperature drop associated with this joint has been the subject of many experimental and theoretical studies. The low temperature behavior of dry joints have shown some anomalous dependence on the surface condition of the mating parts. There is also some doubts on how well one can extrapolate from the test samples to predicting the performance of a real system. Both finite element and analytic models of a simple contact system have been developed. The model assumes (a) the contact is dry (contact limited to a small portion of the total available area and the spaces in-between the actual contact patches are perfect insulators), (b) contacts are clean (conductivity of the actual contact is the same as the bulk), (c) small temperature gradients (the bulk conductance may be assumed to be temperature independent), (d) the absolute temperature is low (thermal radiation effects are ignored), and (e) the dimensions of the nominal contact area are small compared to the thickness of the bulk material (the contact effects are localized near the contact). The models show that in the limit of actual contact area much less than the nominal area (a much less than A), that the excess temperature drop due to a single point of contact scales as a(exp -1/2). This disturbance only extends a distance approx. A(exp 1/2) into the bulk material. A group of identical contacts will result in an excess temperature drop that scales as n(exp -1/2), where n is the number of contacts and n dot a is constant. This implies that flat rough surfaces will have a lower excess temperature drop than flat polished surfaces.

Kittel, Peter↗
NASA NTRS1983

Analytical and experimental study of a standing torus with normal loads

Finite element analysis is applied to study the large deflection and stress of a thin-walled pressurized torus loaded by normal contact with a plane. Experimentally determined stresses are compared to those resulting from finite element analysis at various meridional and circumferential coordinates of the torus, including the footprint area. Circumferential strains compare favorably while meridional strains are higher in the finite element analysis, probably due to slippage of the boundary at the rim. The STAGS finite element computer program was used. The purpose of the study was to evaluate the various program options for structural loading, for material modeling, for stress analysis, and for grid refinement. The experimental model, a thin-walled rubber tube mounted on a steel cylindrical rim, provided measured results to compare with various analytical trade-offs. It was found that there was almost no difference in predicted deflections or stress distributions between linear and nonlinear material description. However, the difference between linear analysis and that of geometric nonlinearities utilizing incremental loading was marked.

Mack, M. J., Jr.↗
NASA NTRS1973

A thermal, thermoelastic, and wear simulation of a high-energy sliding contact problem

This paper describes an investigation of the sliding contact problem encountered in high-energy disk brakes. The analysis includes a simulation modeling, using the finite element method, of the thermoelastic instabilities that cause transient changes in contact to occur on the friction surface. In order to include the effect of wear of the concentrated contacts on the friction surface, a wear criterion is proposed that results in prediction of wear rates for disk brakes that are quite close to experimentally determined wear rates. The thermal analysis shows that the transient temperature distribution in a disk brake can be determined more accurately by use of this thermomechanical analysis than by a more conventional analysis that assumes constant contact conditions. It is also shown that lower, more desirable, temperatures in disk brakes can be attained by increasing the volume, the thermal conductivity, and especially, the heat capacity of the brake components.

Kennedy, F. E., Jr.↗
NASA NTRS1974

A thermal, thermoelastic, and wear analysis of high-energy disk brakes

A thermomechanical investigation of the sliding contact problem encountered in high-energy disk brakes is described. The analysis includes a modelling, using the finite element method of the thermoelastic instabilities that cause transient changes in contact area to occur on the friction surface. In order to include the effect of wear at the contact surface, a wear criterion is proposed that results in the prediction of wear rates for disk brakes that are quite close to experimentally determined wear rates. The thermal analysis shows that the transient temperature distribution in a disk brake assembly can be determined more accurately by use of this thermomechanical analysis than by a more conventional analysis that assumes constant contact conditions. It also shows that lower, more desirable, temperatures in disk brakes can be attained by increasing the volume, the thermal conductivity, and, especially, the heat capacity of the brake components.

Kennedy, F. E., Jr.↗
DOE OSTI2023

Smoothing tool design and performance during subaperture glass polishing

During subaperture tool grinding and polishing, overlaps of the tool influence function can result in undesirable mid-spatial frequency (MSF) errors in the form of surface ripples, which are often corrected using a smoothing polishing step. Here, in this study, flat multi-layer smoothing polishing tools are designed and tested to simultaneously (1) reduce or remove MSF errors, (2) minimize surface figure degradation, and (3) maximize the material removal rate. A time-dependent convergence model in which spatial material removal varies with a workpiece-tool height mismatch, combined with a finite element mechanical analysis to determine the interface contact pressure distribution, was developed to evaluate various smoothing tool designs as a function of tool material properties, thicknesses, pad textures, and displacements. An improvement in smoothing tool performance is achieved when the gap pressure constant, $\overline{h}$ (which describes the inverse rate at which the pressure drops with a workpiece-tool height mismatch), is minimized for smaller spatial scale length surface features (namely, MSF errors) and maximized for large spatial scale length features (i.e., surface figure). Five specific smoothing tool designs were experimentally evaluated. A two-layer smoothing tool using a thin, grooved IC1000 polyurethane pad (with a high elastic modulus, E pad =360 MPa ), thicker blue foam (with an intermediate modulus, E foam =5.3 MPa ) underlayer, and an optimized displacement (d t =1 mm ) provided the best overall performance (namely, high MSF error convergence, minimal surface figure degradation, and high material removal rate).

47 OTHER INSTRUMENTATION↗
NASA NTRS1981

A method of selecting grid size to account for Hertz deformation in finite element analysis of spur gears

A method of selecting grid size for the finite element analysis of gear tooth deflection is presented. The method is based on a finite element study of two cylinders in line contact, where the criterion for establishing element size was that there be agreement with the classical Hertzian solution for deflection. The results are applied to calculate deflection for the gear specimen used in the NASA spur gear test rig. Comparisons are made between the present results and the results of two other methods of calculation. The results have application in design of gear tooth profile modifications to reduce noise and dynamic loads.

Coy, J. J.↗
NASA NTRS1985

Stress analysis method for a clearance fit bolt under bearing loads

In the present, simple analysis method for bolt hole clearance stresses in mechanically fastened joints, an inverse formulation is used with a finite element analysis. Conditions along the bolt hole contact arc are specified by displacement constraint equations. The method was applied to a single fastener clearance-fit joint with typical clearance values, for the cases of a rigid, frictionless bolt and a quasi-isotropic graphite/epoxy laminate. Results indicate that the contact arc, as well as the peak stresses around the hole (and their locations), were strongly influenced by the clearance. After a minor initial nonlinearity, peak stresses varied linearly with applied load. Clearance levels are found to have only a minor influence on overall joint stiffness.

Crews, J. H., Jr.↗
DOE OSTI2021

Nonlinear Interface Reduction for Time-Domain Analysis of Hurty/Craig-Bampton Superelements with Frictional Contact

Virtual prototyping in engineering design relies today on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.

42 ENGINEERING↗
NASA NTRS2006

Bearing-Load Modeling and Analysis Study for Mechanically Connected Structures

Bearing-load response for a pin-loaded hole is studied within the context of two-dimensional finite element analyses. Pin-loaded-hole configurations are representative of mechanically connected structures, such as a stiffener fastened to a rib of an isogrid panel, that are idealized as part of a larger structural component. Within this context, the larger structural component may be idealized as a two-dimensional shell finite element model to identify load paths and high stress regions. Finite element modeling and analysis aspects of a pin-loaded hole are considered in the present paper including the use of linear and nonlinear springs to simulate the pin-bearing contact condition. Simulating pin-connected structures within a two-dimensional finite element analysis model using nonlinear spring or gap elements provides an effective way for accurate prediction of the local effective stress state and peak forces.

Knight, Norman F., Jr.↗
DOE OSTI2022

Scalable preconditioning for the stabilized contact mechanics problem

We present a family of preconditioning strategies for the contact problem in fractured and faulted porous media. We combine low-order continuous finite elements to simulate the bulk deformation with piecewise constant Lagrange multipliers to impose the frictional contact constraints. This formulation is not uniformly inf-sup stable and requires stabilization. We improve previous work by Franceschini et al. (2020) by introducing a novel jump stabilization technique that requires only local geometrical and mechanical properties. We then design scalable preconditioning strategies that take advantage of the block structure of the Jacobian matrix using a physics-based partitioning of the unknowns by field type, namely displacement and Lagrange multipliers. The key to the success of the proposed preconditioners is a pseudo-Schur complement obtained by eliminating the Lagrange multiplier degrees of freedom, which can then be efficiently solved using an optimal multigrid method. Numerical results, including complex real-world problems, are presented to illustrate theoretical properties, scalability and robustness of the preconditioner. A comparison with other approaches available in the literature is also provided.

58 GEOSCIENCES↗
NASA NTRS1994

A Procedure for 3-D Contact Stress Analysis of Spiral Bevel Gears

Contact stress distribution of spiral bevel gears using nonlinear finite element static analysis is presented. Procedures have been developed to solve the nonlinear equations that identify the gear and pinion surface coordinates based on the kinematics of the cutting process and orientate the pinion and the gear in space to mesh with each other. Contact is simulated by connecting GAP elements along the intersection of a line from each pinion point (parallel to the normal at the contact point) with the gear surface. A three dimensional model with four gear teeth and three pinion teeth is used to determine the contact stresses at two different contact positions in a spiral bevel gearset. A summary of the elliptical contact stress distribution is given. This information will be helpful to helicopter and aircraft transmission designers who need to minimize weight of the transmission and maximize reliability.

Kumar, A.↗