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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Mechanically fastened composite laminates subjected to combined bearing-bypass and shear loading

Bolts and rivets provide a means of load transfer in the construction of aircraft. However, they give rise to stress concentrations and are often the source and location of static and fatigue failures. Furthermore, fastener holes are prone to cracks during take-off and landing. These cracks present the most common origin of structural failures in aircraft. Therefore, accurate determination of the contact stresses associated with such loaded holes in mechanically fastened joints is essential to reliable strength evaluation and failure prediction. As the laminate is subjected to loading, the contact region, whose extent is not known, develops between the fastener and the hole boundary through this contact region, which consists of slip and no-slip zones due to friction. The presence of the unknown contact stress distribution over the contact region between the pin and the composite laminate, material anisotropy, friction between the pin and the laminate, pin-hole clearance, combined bearing-bypass and shear loading, and finite geometry of the laminate result in a complex non-linear problem. In the case of bearing-bypass loading in compression, this non-linear problem is further complicated by the presence of dual contact regions. Previous research concerning the analysis of mechanical joints subjected to combined bearing-bypass and shear loading is non-existent. In the case of bearing-bypass loading only, except for the study conducted by Naik and Crews (1991), others employed the concept of superposition which is not valid for this non-linear problem. Naik and Crews applied a linear finite element analysis with conditions along the pin-hole contact region specified as displacement constraint equations. The major shortcoming of this method is that the variation of the contract region as a function of the applied load should be known a priori. Also, their analysis is limited to symmetric geometry and material systems, and frictionless boundary conditions. Since the contact stress distribution and the contact region are not known a priori, they did not directly impose the boundary conditions appropriate for modelling the contact and on-contact regions between the fastener and the hole. Furthermore, finite element analysis is not suitable for iterative design calculations for optimizing laminate construction in the presence of fasteners under complex loading conditions. In this study, the solution method developed by Madenci and Ileri (1992a,b) has been extended to determine the contact stresses in mechanical joints under combined bearing-bypass and shear loading, and bearing-bypass loading in compression resulting in dual contact regions.

Madenci, Erdogan↗
NASA NTRS1982

Wave propagation in graphite/epoxy laminates due to impact

The low velocity impact response of graphite-epoxy laminates is investigated theoretically and experimentally. A nine-node isoparametric finite element in conjunction with an empirical contact law was used for the theoretical investigation. Flat laminates subjected to pendulum impact were used for the experimental investigation. Theoretical results are in good agreement with strain gage experimental data. The collective results of the investigation indicate that the theoretical procedure describes the impact response of the laminate up to about 150 in/sec. impact velocity.

Tan, T. M.↗
NASA NTRS1985

Use of statical indentation laws in the impact analysis of laminated composite plates

The low-velocity impact response of graphite/epoxy laminates was investigated theoretically and experimentally. A nine-node isoparametric plate finite element in conjunction with an empirical contact law was used for the theoretical investigation. Theoretical results are in good agreement with strain-gage experimental data. The results of the investigation indicate that the present theoretical procedure describes the impact response of laminate for low-impact velocities.

Tan, T. M.↗
NASA NTRS1987

A fiber-optic technique for the measurement of contact angle in a clearance-fit pin-loaded hole

A fiber-optic technique for measuring contact angle during pin loading of a specimen is proposed. The experimental design and procedures for loading a 49.8-mm-diameter instrumented pin into an quasi-isotropic graphite-epoxy specimen are described. The optical fiber was located just above the surface of the pin outer diameter in order to obtain accurate pin-hole contact-angle measurements at increasing load levels. The movement of the optical fiber through the no-contact, contact, and no-contact regions is discussed; the photodiode output decreased monotonically as the fiber moved from the no-contact to the contact region and then decreased monotonically as the fiber moved from the contact region to the no-contact region. Variations in the contact angle measurements are examined as function of applied load level. The measurements are compared to contact angle values obtained using a finite element analysis and an electrical technique; it is determined that the data correlate well.

Prabhakaran, R.↗
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.↗
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.↗
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.↗
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.↗
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.↗
NASA NTRS1983

A tire contact solution technique

An efficient method for calculating the contact boundary and interfacial pressure distribution was developed. This solution technique utilizes the discrete Fourier transform to establish an influence coefficient matrix for the portion of the pressurized tire surface that may be in the contact region. This matrix is used in a linear algebra algorithm to determine the contact boundary and the array of forces within the boundary that are necessary to hold the tire in equilibrium against a specified contact surface. The algorithm also determines the normal and tangential displacements of those points on the tire surface that are included in the influence coefficient matrix. Displacements within and outside the contact region are calculated. The solution technique is implemented with a finite-element tire model that is based on orthotropic, nonlinear shell of revolution elements which can respond to nonaxisymmetric loads. A sample contact solution is presented.

Tielking, J. T.↗