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 NTRS1980

Nonlinear transient analysis by energy minimization: A theoretical basis for the ACTION computer code

The formulation basis for establishing the static or dynamic equilibrium configurations of finite element models of structures which may behave in the nonlinear range are provided. With both geometric and time independent material nonlinearities included, the development is restricted to simple one and two dimensional finite elements which are regarded as being the basic elements for modeling full aircraft-like structures under crash conditions. Representations of a rigid link and an impenetrable contact plane are added to the deformation model so that any number of nodes of the finite element model may be connected by a rigid link or may contact the plane. Equilibrium configurations are derived as the stationary conditions of a potential function of the generalized nodal variables of the model. Minimization of the nonlinear potential function is achieved by using the best current variable metric update formula for use in unconstrained minimization. Powell's conjugate gradient algorithm, which offers very low storage requirements at some slight increase in the total number of calculations, is the other alternative algorithm to be used for extremely large scale problems.

Kamat, M. P.↗
NASA NTRS1991

Effects of bolt-hole contact on bearing-bypass damage-onset strength

A combined experimental and analytical study was conducted to investigate the effects of bolt-hole contact on the bearing bypass strength of a graphite-epoxy laminate. Tests were conducted on specimens consisting of 16-ply quasi-isotropic T300/5208 laminates with a centrally located hole. Bearing loads were applied through a clearance-fit steel bolt. Damage onset strength and damage mode were determined for each test case. A finite element procedure was used to calculate the bolt-hole stresses and bolt contact for each test case. A finite element procedure was used to calculate the bolt-hole stresses and bolt contact for each measured damage-onset strength. For the tension bearing-bypass cases tested, the bolt contact half-angle was approximately 60 degrees at damage onset. For compression, the contact angle was 20 degrees as the bypass load increased. A corresponding decrease in the bearing damage onset strength was attributed to the decrease in contact angle which made the bearing loads more severe. Hole boundary stresses were also computed by superimposing stresses for separate bearing and bypass loading. Stresses at the specimen net section were accurately approximated by the superposition procedure. However, the peak bearing stresses had large errors because the bolt contact angles were not represented correctly. For compression, peak bearing stress errors of nearly 50 percent were calculated.

Crews, John H., Jr.↗
DOE OSTI2021

Observations of modal coupling due to bolted joints in an experimental benchmark structure

The goal of this paper is to present a set of measurements from a benchmark structure containing two bolted joints to support future efforts to predict the damping due to the joints and to model nonlinear coupling between the first two elastic modes. Bolted joints introduce nonlinearities in structures, typically causing a softening in the natural frequency and an increase in damping because of frictional slip between the contact interfaces within the joint. These nonlinearities pose significant challenges when characterizing the response of the structure under a large range of load amplitudes, especially when the modal responses become coupled, causing the effective damping and natural frequency to not only depend on the excitation amplitude of the targeted mode, but also the relative amplitudes of other modes. In this work, two nominally identical benchmark structures, known in some prior works as the S4 beam, are tested to characterize their nonlinear properties for the first two elastic modes. Detailed surface measurements are presented and validated through finite element analysis and reveal distinct contact interactions between the two sets of beams. The free-free test structures are excited with an impact hammer and the transient response is analyzed to extract the damping and frequency backbone curves. A range of impact amplitudes and drive points are used to isolate a single mode or to excite both modes simultaneously. Differences in the nonlinear response correlate with the relative strength of the modes that are excited, allowing one to characterize mode coupling. Each of the beams shows different nonlinear properties for each mode, which is attributed to the different contact pressure distributions between the parts, although the mode coupling relationship is found to be consistent between the two. Lastly, the test data key finding are presented in this paper and the supporting data is available on a public repository for interested researchers.

42 ENGINEERING↗
NASA NTRS1985

Transient and steady state viscoelastic rolling contact

Based on moving total Lagrangian coordinates, a so-called traveling Hughes type contact strategy is developed. Employing the modified contact scheme in conjunction with a traveling finite element strategy, an overall solution methodology is developed to handle transient and steady viscoelastic rolling contact. To verify the scheme, the results of both experimental and analytical benchmarking is presented. The experimental benchmarking includes the handling of rolling tires up to their upper bound behavior, namely the standing wave response.

Padovan, J.↗
NASA NTRS1992

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.↗
NASA NTRS1989

Stress analysis method for clearance-fit joints with bearing-bypass loads

Within a multi-fastener joint, fastener holes may be subjected to the combined effects of bearing loads and loads that bypass the hole to be reacted elsewhere in the joint. The analysis of a joint subjected to search combined bearing and bypass loads is complicated by the usual clearance between the hole and the fastener. A simple analysis method for such clearance-fit joints subjected to bearing-bypass loading has been developed in the present study. It uses an inverse formulation with a linear elastic finite-element analysis. Conditions along the bolt-hole contact arc are specified by displacement constraint equations. The present method is simple to apply and can be implemented with most general purpose finite-element programs since it does not use complicated iterative-incremental procedures. The method was used to study the effects of bearing-bypass loading on bolt-hole contact angles and local stresses. In this study, a rigid, frictionless bolt was used with a plate having the properties of a quasi-isotropic graphite/epoxy laminate. Results showed that the contact angle as well as the peak stresses around the hole and their locations were strongly influenced by the ratio of bearing and bypass loads. For single contact, tension and compression bearing-bypass loading had opposite effects on the contact angle. For some compressive bearing-bypass loads, the hole tended to close on the fastener leading to dual contact. It was shown that dual contact reduces the stress concentration at the fastener and would, therefore, increase joint strength in compression. The results illustrate the general importance of accounting for bolt-hole clearance and contact to accurately compute local bolt-hole stresses for combined bearings and bypass loading.

Naik, R. A.↗
NASA NTRS2001

Design, Generation and Tooth Contact Analysis (TCA) of Asymmetric Face Gear Drive With Modified Geometry

A new type of face gear drive for application in transmissions, particularly in helicopters, has been developed. The new geometry differs from the existing geometry by application of asymmetric profiles and double-crowned pinion of the face gear mesh. The paper describes the computerized design, simulation of meshing and contact, and stress analysis by finite element method. Special purpose computer codes have been developed to conduct the analysis. The analysis of this new type of face gear is illustrated with a numerical example.

Litvin, Faydor L.↗
NASA NTRS2011

On Calculation Methods and Results for Straight Cylindrical Roller Bearing Deflection, Stiffness, and Stress

The purpose of this study was to assess some calculation methods for quantifying the relationships of bearing geometry, material properties, load, deflection, stiffness, and stress. The scope of the work was limited to two-dimensional modeling of straight cylindrical roller bearings. Preparations for studies of dynamic response of bearings with damaged surfaces motivated this work. Studies were selected to exercise and build confidence in the numerical tools. Three calculation methods were used in this work. Two of the methods were numerical solutions of the Hertz contact approach. The third method used was a combined finite element surface integral method. Example calculations were done for a single roller loaded between an inner and outer raceway for code verification. Next, a bearing with 13 rollers and all-steel construction was used as an example to do additional code verification, including an assessment of the leading order of accuracy of the finite element and surface integral method. Results from that study show that the method is at least first-order accurate. Those results also show that the contact grid refinement has a more significant influence on precision as compared to the finite element grid refinement. To explore the influence of material properties, the 13-roller bearing was modeled as made from Nitinol 60, a material with very different properties from steel and showing some potential for bearing applications. The codes were exercised to compare contact areas and stress levels for steel and Nitinol 60 bearings operating at equivalent power density. As a step toward modeling the dynamic response of bearings having surface damage, static analyses were completed to simulate a bearing with a spall or similar damage.

Krantz, Timothy L.↗
NASA NTRS1985

Dynamic large deflection response of composite laminates subjected to impact

Impact responses of composite laminates with and without initial stresses are investigated using the finite element method. A nine-node isoparametric quadrilateral element based on the Mindlin plate theory and the von Karman large deflection assumptions is developed. An experimentally established contact law which accounts for the permanent indentation is incorporated into the finite element program to evaluate the impact force. In the time integration, the Newmark constant acceleration algorithm is used in conjunction with successive iterations within each time step. Numerical results, including the contact force histories, deflections and strains in the plate, are presented.

Chen, J. K.↗
NASA NTRS1994

Manual for automatic generation of finite element models of spiral bevel gears in mesh

The goal of this research is to develop computer programs that generate finite element models suitable for doing 3D contact analysis of faced milled spiral bevel gears in mesh. A pinion tooth and a gear tooth are created and put in mesh. There are two programs: Points.f and Pat.f to perform the analysis. Points.f is based on the equation of meshing for spiral bevel gears. It uses machine tool settings to solve for an N x M mesh of points on the four surfaces, pinion concave and convex, and gear concave and convex. Points.f creates the file POINTS.OUT, an ASCI file containing N x M points for each surface. (N is the number of node points along the length of the tooth, and M is nodes along the height.) Pat.f reads POINTS.OUT and creates the file tl.out. Tl.out is a series of PATRAN input commands. In addition to the mesh density on the tooth face, additional user specified variables are the number of finite elements through the thickness, and the number of finite elements along the tooth full fillet. A full fillet is assumed to exist for both the pinion and gear.

Bibel, G. D.↗
DOE OSTI2026

C‐Cracking of Brittle‐Material Spheres From Eccentric Hertzian Contact

Contact of brittle-material spheres can occur while they are manufactured, handled, or used in operation as roller elements or in deformably processed composite material. Statistically, the directional contact between two spheres will nearly always be oblique or eccentric in an unconfined and mobile sphere population. Although analyses of oblique contact exist, the portrayal of its maximum (tensile) first principal stress ( S 1 ) field is lacking. This information is needed for improved judgment of the prospect of crack initiation. Given these, the S 1 due to eccentric contact and subsequent crack initiation was analyzed using finite element analysis (FEA) and corroborative demonstration of c-crack creation. The FEA shows that an asymmetric S 1 field is created about the contact patch and is caused by superimposed shear intrinsic to eccentric contact. A crack will initiate if the maximum S 1 exceeds the material's tensile strength, and its field will cause the crack to propagate and arrest to form a C-shaped crack, or c-crack. C-crack initiation trends with lower applied forces, with greater amount of contact eccentricity, higher coefficient of friction, and smaller spherical radius. Finally, these metrics deserve recognition, especially when c-cracked brittle-material spheres are used in a system whose functionality or reliability is potentially compromised by the c-cracks.

Hertzian contact↗
NASA NTRS2000

Simulation of Hypervelocity Impact on Aluminum-Nextel-Kevlar Orbital Debris Shields

An improved hybrid particle-finite element method has been developed for hypervelocity impact simulation. The method combines the general contact-impact capabilities of particle codes with the true Lagrangian kinematics of large strain finite element formulations. Unlike some alternative schemes which couple Lagrangian finite element models with smooth particle hydrodynamics, the present formulation makes no use of slidelines or penalty forces. The method has been implemented in a parallel, three dimensional computer code. Simulations of three dimensional orbital debris impact problems using this parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials. at an impact velocity of eleven kilometers per second.

Fahrenthold, Eric P.↗
DOE OSTI2023

Application of the van der Pauw method for electrical conductivity measurements at high temperatures using an insulating compressing ring

High quality data on the high temperature electrical properties of ceramics, particularly oxides, is of great value for material selection, design, and modeling for a broad range of emerging applications. Utilizing the mismatch in the coefficient of thermal expansion between two materials, a purely mechanical method for establishing electrical contact in the van der Pauw geometry to measure the bulk resistivity of ceramic disks at high temperatures is presented. Measurements of a reference material, 20 mol. % Gd-doped cerium oxide, are presented up to 1000 °C. The viability of electrical measurements up to a maximum temperature of 1600 °C is also considered. Measurements are performed using multiple techniques and compared to literature values finding excellent agreement. Finally, the approach described in this work enables the van der Pauw method to be applied to many ceramic materials over a wide range of temperatures and environments.

47 OTHER INSTRUMENTATION↗
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↗
DOE OSTI2020

HOSS: an implementation of the combined finite-discrete element method

Nearly thirty years since its inception, the combined finite-discrete element method (FDEM) has made remarkable strides in becoming a mainstream analysis tool within the field of Computational Mechanics. FDEM was developed to effectively “bridge the gap” between two disparate Computational Mechanics approaches known as the finite and discrete element meth-ods. At Los Alamos National Laboratory (LANL) researchers developed the Hybrid Optimization Software Suite (HOSS) as a hybrid multi-physics platform, based on FDEM, for the simulation of solid material behavior complemented with the latest technological enhancements for full fluid–solid interaction. Furthermore, in HOSS, several newly developed FDEM algorithms have been implemented that yield more accurate material deformation formulations, inter-particle interaction solvers, and fracture and fragmentation solutions. Additionally, an explicit computational fluid dynamics solver and a novel fluid–solid interaction algorithms have been fully integrated (as opposed to coupled) into the HOSS’ solid mechanical solver, allowing for the study of an even wider range of problems. Advancements such as this are leading HOSS to become a tool of choice for multi-physics problems. Finally, HOSS has been successfully applied by a myriad of researchers for analysis in rock mechanics, oil and gas industries, engineering application (structural, mechanical and biomedical engineering), mining, blast loading, high velocity impact, as well as seismic and acoustic analysis. This paper intends to summarize the latest development and application efforts for HOSS.

58 GEOSCIENCES↗
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.↗