SEARCH · Engineering Papers

Results for “contact model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

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 NTRS1979

Two-dimensional random surface model for asperity-contact in elastohydrodynamic lubrication

Relations for the asperity-contact time function during elastohydrodynamic lubrication of a ball bearing are presented. The analysis is based on a two-dimensional random surface model, and actual profile traces of the bearing surfaces are used as statistical sample records. The results of the analysis show that transition from 90 percent contact to 1 percent contact occurs within a dimensionless film thickness range of approximately four to five. This thickness ratio is several times large than reported in the literature where one-dimensional random surface models were used. It is shown that low pass filtering of the statistical records will bring agreement between the present results and those in the literature.

Coy, J. J.↗

NASA NTRS1979

Two-dimensional random surface model for asperity-contact in elastohydrodynamic lubrication

Relations for the asperity-contact time fraction during elastohydrodynamic lubrication of a ball bearing are presented. The analysis is based on a two-dimensional random surface model, and actual profile traces of the bearing surfaces are used as statistical sample records. The results of the analysis show that transition from 90 percent contact to 1 percent contact occurs within a dimensionless film thickness range of approximately four to five. This thickness ratio is several times larger than reported in the literature where one-dimensional random surface models were used. It is shown that low pass filtering of the statistical records will bring agreement between the present results and those in the literature.

Coy, J. J.↗

NASA NTRS1983

Subsurface stress evaluations under rolling/sliding contacts

A computer model has been developed for evaluating the subsurface stresses incurred within rolling/sliding (elastic) contacts. The model involves first defining the stress tensor at any point (x, y, or z) beneath the surface in terms of the surface stresses. The stress tensors are analyzed to determine the maximum shear stresses and stress reversals. As a result of computations with the model, several observations were made. For example, the maximum reversing shear stresses are on the plane of the orthogonal shear stress. Further, the magnitude of these stresses is not altered by friction. However, under very high friction (typical of dry contact) surface stresses can dominate over subsurface stresses.

Kannel, J. W.↗

NASA NTRS1978

Statistical model for asperity-contact time fraction in elastohydrodynamic lubrication

Relations for the asperity contact time fraction during elastohydrodynamic (EHD) lubrication of a typical ball bearing are presented. The analysis is based on a two-dimensional random surface model, and actual profile traces of the bearing surfaces were used as statistical sample records. The results of the analysis show that transition from 90 percent contact to 1 percent contact occurs within a dimensionless film thickness range of approximately 4 to 5. This thickness ratio is several times larger than reported in the literature where one-dimensional random surface models were used.

Sidik, S. M.↗

NASA NTRS2018

Comparison of Performance Test Results to CFD and Structural Models of Non-Contacting Finger Seals

Performance tests of non-contacting finger seal designs were conducted at 300, 700, 922 K (70, 800 and 1200 F) at pressure differentials up to 517 kPa (75 psid) and surface speeds up to 366 m/s (1200 ft/s). Room temperature, static analysis of the seal was performed. A simplified CFD model was developed to examine pressure loads within the seal. Results from the CFD model were used as input to a finite element analysis model of a six-finger segment of the non-contacting finger seal. Examination of predicted deflections of individual components of the seal gives insight into the seal behavior. Wear patterns from testing verify the pattern of radial deflection. The models are used to predict maximum pressure differential capability of the seal and compared to experimental results. The CFD model slightly under-predicts the measured leakage flow factor, but has the same trend as the measured flow factor versus pressure differential.

Seals:Brush Seals↗

NASA NTRS1990

Graphics modelling of non-contact thickness measuring robotics work cell

A system was developed for measuring, in real time, the thickness of a sprayable insulation during its application. The system was graphically modelled, off-line, using a state-of-the-art graphics workstation and associated software. This model was to contain a 3D color model of a workcell containing a robot and an air bearing turntable. A communication link was established between the graphics workstations and the robot's controller. Sequences of robot motion generated by the computer simulation are transmitted to the robot for execution.

Warren, Charles W.↗

DOE OSTI2022

Elevated temperature contact creep and friction of nickel-based superalloys using machine learning assisted finite element analysis

Nickel-based superalloys with superior thermochemical, mechanical, and tribological properties are highly utilized for critical components in several high temperature applications such as gas turbines and nuclear reactors. Inconel 617, in particular, is considered as one of the main candidate superalloys for tribo-components in very-high-temperature gas-cooled nuclear reactors. Recent findings indicate that this alloy grows unique surface oxide especially in a high-temperature helium environment with distinctive wear, friction, and contact properties. Here, this study investigates the high temperature contact area evolution and frictional behavior of Inconel 617 using finite element simulation and provides predictive models for the contact and friction performance at different normal loads, dwell times, and temperatures. High temperature helium-aged Inconel 617 top surface properties (up to 600 °C) are utilized along with a single asperity-based deformable elastic-plastic contact model under combined normal and tangential loading. Machine learning is used to assist the finite element results and to predict friction coefficient as well as contact area evolution. While a small difference is observed in the instantaneous friction coefficient (no dwell time) for all temperatures, friction coefficient increases considerably with dwell time. This shows that the effect of contact creep for longer dwell times significantly dominants the effect of high temperature variation in basic mechanical parameters such as modulus and yield strength. It is found that increasing temperature and dwell times lead to the friction coefficient increase, yet the dominance of dwell time effects decreases at higher temperatures and loads. While the analysis is presented for Inconel 617, the methodology is easy to be generalized and can be applied to other HT alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Big Hill Geomechanical Analysis Using M-D Viscoplastic Material Model

The integrity of wellbores at the interbed between the caprock and salt is a serious concern in the Big Hill site. For the remediation and life extension of wellbores, more accurate predictions from the global model are needed. The Big Hill global model is improved using the M-D viscoplastic contact surface model and the mesh containing the interbed layer with contact surfaces between the salt and caprock layers, and fault blocks in overburden and caprock layers. The model calibration has been performed based on the cavern volumetric closures obtained from the Caveman calculations. The results agree well from 1991 to the early 2000s. The difference starts to widen after that, it might be because of frequent fluid movement and raw water injection. Therefore, the predictions from this improved model could be used to examine the structural integrity of caverns in Big Hill salt dome.

42 ENGINEERING↗

DOE OSTI2020

A hybrid calibration approach to Hertz-type contact parameters for discrete element models

This study aims at providing a hybrid calibration framework to estimate Hertz-type contact parameters (particle-scale shear modulus and Poisson ratio) for both two-dimensional and three-dimensional discrete element modelling (DEM). On the basis of statistically isotropic granular packings, a set of analytical formulae between macroscopic material parameters (Young modulus and Poisson ratio) and particle-scale Hertz-type contact parameters for granular systems are derived under small-strain isotropic stress conditions. However, the derived analytical solutions are only estimated values for general models. By viewing each DEM modelling as an implicit mathematical function taking the particle-level parameters as independent variables and employing the derived analytical solutions as the initial input parameters, an automatic iterative scheme is proposed to obtain the calibrated parameters with higher accuracies. Considering highly nonlinear features and discontinuities of the macro-micro relationship in Hertz-based discrete element models, the adaptive moment estimation algorithm is adopted in this study because of its capacity of dealing with noise gradients of cost functions. Here, the proposed method is validated with several numerical cases including randomly distributed monodisperse and polydisperse packings. Noticeable improvements in terms of calibration efficiency and accuracy have been made.

Constitutive law↗

DOE OSTI2023

Separation of volatile organic contaminants from water using a direct-contact dehumidifier: An experimental study and modeling

A packed-bed, direct-contact dehumidifier (PB-DCD) has been developed to selectively condense water from a gaseous mixture of organic contaminants. In this process, hot air is humidified by mixing with wastewater vapor in a static mixer and is dehumidified by subcooled fresh water passing through the PB-DCD in a counter-current configuration. The goal of this study is to understand the relationship between various operating parameters and how they affect separation of isopropyl alcohol (IPA) from the gaseous mixture with the purpose of expanding the investigation to other contaminants in the future. The process has been successfully used to separate IPA from wastewater. Different parameters including the concentration of IPA in water, feed temperature, air-to-vapor ratio, and cooling-water-to-vapor mass ratio, as well as residence time were examined to determine their influence on contaminant separation and clean water recovery rate. The Buckingham Pi theorem was applied to reduce the number of studying variables by generating non-dimensional groups. Experiments were performed to determine the relationship between the non-dimensional parameters and model mass transfer in the system. Mass transfer coefficients for multicomponent gaseous mixture condensation in a packed-bed direct contact condenser column are calculated and a modified correlation of mass transfer is developed. The parameters of the mass transfer model are obtained from operating conditions of the PB-DCD, stream conditions and thermodynamic properties. A modified Sherwood correlation is developed to predict IPA separation from a gaseous mixture consisting of air, water, and IPA. The modified Sherwood correlation is applied to successfully predict the Sherwood number with a mean absolute error equal to 6.6 %. According to the experimental results, the maximum IPA separation achieved was 76.3 %, which corresponds to a water recovery rate of 58.7 %.

42 ENGINEERING↗

NASA NTRS1981

Aeroelastic instability stoppers for wind tunnel models

A mechanism for constraining models or sections thereof, was wind tunnel tested, deployed at the onset of aeroelastic instability, to forestall destructive vibrations in the model is described. The mechanism includes a pair of arms pivoted to the tunnel wall and straddling the model. Rollers on the ends of the arms contact the model, and are pulled together against the model by a spring stretched between the arms. An actuator mechanism swings the arms into place and back as desired.

Doggett, R. V., Jr.↗

NASA NTRS1996

An Idealized Direct-Contact Biomass Pyrolysis Reactor Model

A numerical study is performed in order to assess the performance of biomass pyrolysis reactors which utilize direct particle-wall thermal conduction heating. An idealized reactor configuration consisting of a flat-plate turbulent boundary layer flow with particle convection along the heated wall and incorporating particle re-entrainment is considered.

biomass boundary layer modeling pyrolysis tar vort↗

DOE OSTI2021

Effect of workpiece curvature on the tool influence function during hemispherical sub-aperture tool glass polishing

The influence of workpiece curvature on the tool influence function spot during polishing of fused silica glass with cerium oxide slurry, while using a rotating hemispherical pad-foam tool for a wide variety of process conditions (tool displacement, inclination angle, and rotation rate), has been investigated. (Workpiece curvature ranged from 500 mm radius concave to 43 mm radius convex.) The TIF spot decreases in diameter and increases in the peak removal rate on more convex workpieces. In contrast, the TIF spot increases both in diameter and peak removal rate on more concave workpieces. For the range of workpiece curvatures investigated, both the spot size and the peak removal rate changed significantly, as much as 2 times. An elastic sphere–sphere contact mechanics model, which utilizes both a modified displacement (that leads to a change in the applied load) as well as a mismatch factor (that influences the pressure distribution shape), has been developed. Furthermore, the model was validated using both offline load-displacement measurements and finite-element analysis simulations. The model quantitatively describes the measured change in the relative contact diameter and relative pressure distribution, as well as semiquantitively describes the change in the relative volumetric removal rate on a large variety of TIF spots. The change in the volumetric removal rate for convex workpieces is a result of the balance between a decreasing spot size (reducing removal) and an increasing peak pressure (increasing removal), which usually results in relatively small changes in volumetric removal. In the case of concave workpieces, the volumetric removal rate change is also governed by a similar balance, but the spot size increase contribution dominates, resulting in a significant increase in volumetric removal rate. Understanding these trends can enable methods to add greater determinism during the fabrication of freeform optics by adjusting polishing parameters (such as dwell time) while the tool translates along a workpiece surface with different local curvatures.

47 OTHER INSTRUMENTATION↗

NASA NTRS2001

Coupling of Armchair and Zigzag Tubes to a Free Electron Metal

The effect of nanotube chirality is of prime importance in determining its electronic properties. We address the issue of how chirality affects the coupling of a nanotube to metal contacts. We model coupling of armchair and zigzag nanotubes to metal contacts, in both the side- and end-contacted geometries. In the side-contacted geometry, the coupling of armchair and metallic-zigzag nanotubes to a free electron metal are significantly different. Namely, it is possible to drive a larger current through a metallic-zigzag nanotube. The predicted difference holds good when both (a) the entire circumference and (b) only a finite sector of the nanotube makes contact to the metal electrode. It might be possible to observe the predicted difference between armchair and zigzag nanotubes using gold contacts.

Anantram, M. P.↗

DOE OSTI2020

Microscale mechanical modeling of deformable geomaterials with dynamic contacts based on the numerical manifold method

Abstract Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.

58 GEOSCIENCES↗

NASA NTRS1980

Solutions of contact problems by the assumed stress hybrid model

A method was developed for contact problems which may be either frictional or frictionless and may involve extensive sliding between deformable bodies. It was based on an assumed stress hybrid approach and on an incremental variational principle for which the Euler's equations of the functional include the equilibrium and compatibility conditions at the contact surface. The tractions at an assumed contact surface were introduced as Lagrangian multipliers in the formulation. It was concluded from the results of several example solutions that the extensive sliding contact between deformable bodies can be solved by the present method.

Kubomura, K.↗

DOE OSTI2020

BISON Methods and Contact Algorithms

This paper reviews five critical algorithms for BISON’s simulation of nuclear fuel: the burnup model, the coolant channel model, fuel cracking models, mechanical contact, and thermal contact.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗