Microasperity model for elastohydrodynamic lubrication of a spinning ball on a flat surface
Dry contact deformation theory and point contact elastohydrodynamic theory for microasperity model of spinning ball on flat surface
Allen, C. W.↗
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Surface temperatures and glassy state investigations in tribology, part 2
Measurements of lubricant shear rheological behavior in the amorphous solid region and near the liquid solid transition are reported. Elastic, plastic and viscous behavior was observed. The maximum yield shear stress (limiting shear stress) is a function of temperature and pressure and is believed to be the property which determines the maximum traction in elastohydrodynamic contacts such as traction drives. A shear rheological model based on primary laboratory data is proposed for concentrated contact lubrication. The model is Maxwell model modified with a limiting shear stress. Three material properties are required: low shear stress viscosity, limiting elastic shear modulus, and the limiting shear stress the material can withstand. All three are functions of temperature and pressure.
Bair, S. S.↗
Model Correlation and Thermal Analysis of xEMU Boot at Lunar South Pole Temperatures
The National Aeronautics and Space Administration (NASA) Artemis program plans to send astronauts to the lunar south pole, a region of the moon that is colder than previous lunar missions and low earth orbit operations. The spacesuit boots that will be used on these missions will be directly impacted by these extremely cold temperatures (down to ~50 K). To assess the performance of the Exploration Extravehicular Mobility Unit (xEMU) lunar boot in these extreme temperatures, the boot was tested at the Jet Propulsion Lab (JPL) in the Cryogenic Ice Transfer, Acquisition Development, and Excavation Laboratory (CITADEL) thermal vacuum (TVAC) chamber. This data was used to correlate thermal models to predict operational performance of the boots on the lunar south pole. This paper documents the data analysis, thermal boot model correlation, and operational predictions conducted using data from the xEMU TVAC test. Data from the test series was used to determine expected thermal conductances within the boot and between the boot and environment. These conductances were used to correlate a Thermal Desktop model of the xEMU boot across 10 different test points which varied external temperature, internal boot ventilation flowrate, and contact pressure. The correlated model was then used to predict operational performance in the lunar south pole. While the predictions indicate promising evidence for performance of the boots at the 100K test point, there is still substantial uncertainty in performance, particularly at the 48K test point. The results of this test series and model correlation stress the importance of improved testing for characterizing the expected thermal resistance between the outside of the boot and the lunar surface.
thermal analysis↗
Model Correlation and Thermal Analysis of xEMU Boot at Lunar South Pole Temperatures
The National Aeronautics and Space Administration (NASA) Artemis program plans to send astronauts to the lunar south pole, a region of the moon that is colder than previous lunar missions and low earth orbit operations. The spacesuit boots that will be used on these missions will be directly impacted by these extremely cold temperatures (down to ~50 K). To assess the performance of the Exploration Extravehicular Mobility Unit (xEMU) lunar boot in these extreme temperatures, the boot was tested at the Jet Propulsion Lab (JPL) in the Cryogenic Ice Transfer, Acquisition Development, and Excavation Laboratory (CITADEL) thermal vacuum (TVAC) chamber. This data was used to correlate thermal models to predict operational performance of the boots on the lunar south pole. This paper documents the data analysis, thermal boot model correlation, and operational predictions conducted using data from the xEMU TVAC test. Data from the test series was used to determine expected thermal conductances within the boot and between the boot and environment. These conductances were used to correlate a Thermal Desktop model of the xEMU boot across 10 different test points which varied external temperature, internal boot ventilation flowrate, and contact pressure. The correlated model was then used to predict operational performance in the lunar south pole. While the predictions indicate promising evidence for performance of the boots at the 100K test point, there is still substantial uncertainty in performance, particularly at the 48K test point. The results of this test series and model correlation stress the importance of improved testing for characterizing the expected thermal resistance between the outside of the boot and the lunar surface.
xEMU↗
Theoretical and experimental study of thermal contact resistance in a vacuum environment Semiannual status report no. 8, Jun. - Dec. 1965
Thermal contact resistance in vacuum environment using model predictions
CONTACT RESISTANCE↗
Results of tests on a Rockwell International space shuttle orbiter (-139 configuration) 0.0175-scale model (no. 29-0) in AEDC tunnel B to determine boundary layer characteristics
Results of wind tunnel tests were conducted to determine boundary layer characteristics on the lower surface of a space shuttle orbiter. Total pressure and temperature profile data at various model stations were obtained using a movable, four-degree-of-freedom probe mechanism and static pressure taps on the model surface. During a typical run, the probe was located over a preselected model location, then driven down through the bondary layer until contact was made with the model surface.
Quan, M.↗
Grain dynamics in zero gravity
The dynamics of granular materials has proved difficult to model, primarily because of the complications arising from inelastic losses, friction, packing, and the effect of many grains being in contact simultaneously. The kinetic model of granular systems is similar to the kinetic theory of gases, except that collisional energy losses are always present in the former and must be treated explicity. Few granular materials on Earth are describable by this limiting model, since gravity tends to collapse the grains into a high density state where Coulombic friction effects are dominant. The planned Space Station offers an unusual opportunity to test the kinetic grain model and to explore its predictions. Without gravity, the regime of low interparticle velocities (where an elastic description of the collision is still valid) can be investigated. This will allow for direct interpretation by dynamical computer simulations as well as by the kinetic theory. The dynamics of spherical grains inside a clear box would be examined. Results would be compared with the predictions of the kinetic theory and computer simulations.
Werner, B. T.↗
Computer Simulation of the Forces Acting on the Polystyrene Probe Submerged into the Succinonitrile Near Phase Transition
Modeling approach to simulate both mesoscale and microscopic forces acting in a typical AFM experiment is presented. At mesoscale level interaction between the cantilever tip and the sample surface is primarily described by the balance of attractive Van der Waals and repulsive forces. The model of cantilever oscillations is applicable to both non-contact and "tapping" AFM. This model can be farther enhanced to describe nanoparticle manipulation by cantilever. At microscopic level tip contamination and details of tip-surface interaction can be simulated using molecular dynamics approach. Integration of mesoscale model with molecular dynamic model is discussed.
Bune, Andris V.↗
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↗
Working toward a three-dimensional fatigue closure model for surface cracks
The first reliable elastic fracture mechanics solutions for a surface crack in a plate were obtained by Newman and Raju. The authors, both from the Mechanics of Materials Branch at NASA-Langley, used a highly detailed finite element solution requiring substantial computational resources. Computers have since become more powerful and available; however, many important related problems remain computationally expensive. The problem of three-dimensional fatigue crack growth taking into account plasticity-induced crack closure is one such problem. It is the goal of this research to provide an efficient method to account for three-dimensional crack closure in fatigue. Newman developed a two-dimensional plasticity-induced crack closure model for center cracked specimens. This model requires iterations to determine both the contact solution at each growth step and the extent of the plastic zone at the crack tip. A three-dimensional version of this model would require obtaining these nonlinear variables all along the crack front. This model must be efficient so that repeated calculations can be performed during crack growth simulations. The highly versatile line spring model (LSM) with contact, fatigue, and plasticity will form the basis of the closure model. There are several required additions to past work to address the three-dimensional crack closure problem. Initially, these additions will include (1) an improved LSM to more accurately obtain the crack opening displacement, stress intensity factors, and elastic T-stress near the ends of the surface crack; (2) a method to determine the extent of the plastic zone all along the crack front; (3) a method to determine the contact zone given a perfectly plastic layer of material on the crack surfaces; (4) a method to determine the magnitude of the compressive contact stress; and (5) a way to implement the degree of constraint along the curved crack front. During the summer ASEE program an enhanced LSM was developed. A method similar to that of 'strip synthesis' first introduced by Fujimoto was used. Briefly, the crack opening displacements of 'slices' of the surface crack in a direction parallel to the plate surface are considered in addition to the standard LSM approach that makes use of springs obtained from slices perpendicular to the plate surface. This enhancement is necessary so that an accurate three-dimensional representation of quantities such as contact zone size, plastic zone size, stress intensity factors, T-stress, and crack opening displacement can be determined. By combining results of previous investigations with the LSM, the problem of three-dimensional crack closure will be addressed. In addition to a closure model, the enhanced LSM can be used for many other problems including interacting surface cracks and fatigue crack growth of a through crack with a curved crack front.
Joseph, Paul F.↗
Accurate upwind methods for the Euler equations
A new class of piecewise linear methods for the numerical solution of the one-dimensional Euler equations of gas dynamics is presented. These methods are uniformly second-order accurate, and can be considered as extensions of Godunov's scheme. With an appropriate definition of monotonicity preservation for the case of linear convection, it can be shown that they preserve monotonicity. Similar to Van Leer's MUSCL scheme, they consist of two key steps: a reconstruction step followed by an upwind step. For the reconstruction step, a monotonicity constraint that preserves uniform second-order accuracy is introduced. Computational efficiency is enhanced by devising a criterion that detects the 'smooth' part of the data where the constraint is redundant. The concept and coding of the constraint are simplified by the use of the median function. A slope steepening technique, which has no effect at smooth regions and can resolve a contact discontinuity in four cells, is described. As for the upwind step, existing and new methods are applied in a manner slightly different from those in the literature. These methods are derived by approximating the Euler equations via linearization and diagonalization. At a 'smooth' interface, Harten, Lax, and Van Leer's one intermediate state model is employed. A modification for this model that can resolve contact discontinuities is presented. Near a discontinuity, either this modified model or a more accurate one, namely, Roe's flux-difference splitting. is used. The current presentation of Roe's method, via the conceptually simple flux-vector splitting, not only establishes a connection between the two splittings, but also leads to an admissibility correction with no conditional statement, and an efficient approximation to Osher's approximate Riemann solver. These reconstruction and upwind steps result in schemes that are uniformly second-order accurate and economical at smooth regions, and yield high resolution at discontinuities.
Huynh, Hung T.↗
Optimization of transistor design including large signal device/circuit interactions at extremely high frequencies (20-100+GHz)
Transistor design for extremely high frequency applications requires consideration of the interaction between the device and the circuit to which it is connected. Traditional analytical transistor models are to approximate at some of these frequencies and may not account for variations of dopants and semiconductor materials (especially some of the newer materials) within the device. Physically based models of device performance are required. These are based on coupled systems of partial differential equations and typically require 20 minutes of Cray computer time for a single AC operating point. A technique is presented to extract parameters from a few partial differential equation solutions for the device to create a nonlinear equivalent circuit model which runs in approximately 1 second of personal computer time. This nonlinear equivalent circuit model accurately replicates the contact current properties of the device as computed by the partial differential solver on which it is based. Using the nonlinear equivalent circuit model of the device, optimization of systems design can be performed based on device/circuit interactions.
Levy, Ralph↗
Subharmonic Oscillations and Chaos in Dynamic Atomic Force Microscopy
The increasing use of dynamic atomic force microscopy (d-AFM) for nanoscale materials characterization calls for a deeper understanding of the cantilever dynamics influencing scan stability, predictability, and image quality. Model development is critical to such understanding. Renormalization of the equations governing d- AFM provides a simple interpretation of cantilever dynamics as a single spring and mass system with frequency dependent cantilever stiffness and damping parameters. The renormalized model is sufficiently robust to predict the experimentally observed splitting of the free-space cantilever resonance into multiple resonances upon cantilever-sample contact. Central to the model is the representation of the cantilever sample interaction force as a polynomial expansion with coefficients F(sub ij) (i,j = 0, 1, 2) that account for the effective interaction stiffness parameter, the cantilever-to-sample energy transfer, and the amplitude of cantilever oscillation. Application of the Melnikov method to the model equation is shown to predict a homoclinic bifurcation of the Smale horseshoe type leading to a cascade of period doublings with increasing drive displacement amplitude culminating in chaos and loss of image quality. The threshold value of the drive displacement amplitude necessary to initiate subharmonic generation depends on the acoustic drive frequency, the effective damping coefficient, and the nonlinearity of the cantilever-sample interaction force. For parameter values leading to displacement amplitudes below threshold for homoclinic bifurcation other bifurcation scenarios can occur, some of which lead to chaos.
Cantrell, John H.↗
Modeling of Response of Thermoplastic Composites Under High Energy Impact
This work investigates the response of thermoplastic composites under impact loading. Panels fabricated with the material system under consideration – carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK) – were subject to impact with an Aluminum projectile, at low and high velocities. At lower velocities, the projectile rebounded from the panel after impact and caused damage e.g. delaminations. Finite element models were developed to model this test using LS-DYNA®, with material model MAT299, a continuum damage model based on the LaRC04 criteria, to capture the dynamic response and damage modes in the T700/LMPAEK panels. Interlaminar delaminations were modeled using the cohesive contact formulation available within LS-DYNA. Previous work by the authors showed that delaminations were consistently underpredicted in size, relative to the test. It was surmised that this was due to effects of strain rate on intralaminar and interlaminar response of the T700/LMPAEK material system. Therefore, the current work is focused on incorporating strain rate effects in modeling the interlaminar response. Rate dependence was incorporated using ‘one-way tie-break contact’ in LS-DYNA, wherein, the Mode-II fracture toughness of the LMPAEK resin, G IIc , was assumed to be an exponential function of the strain rate. In particular, G IIc , was modeled to exponentially decrease at higher strain rates. Models with this form of rate dependence favorably correlated with experimental data regarding delamination extent and rebound velocity, relative to prior models without strain rate effects.
Composite Materials↗
A film thickness analysis for line contacts under pure rolling conditions with a non-Newtonian rheological model
The non-Newtonian constitutive equation of Winer and Bair (1979) is applied to a conventional isothermal film thickness analysis of line contact lubrication of rolling elements providing dimensionless film thickness equations for four lubrication methods. The equations can be applied to high viscosity and rolling speed and low limiting shear stress cases where Newtonian methods are not applicable. This analysis is based on a reasonable range of limiting shear stress which is smaller than the Newtonian values for low viscosity lubricants.
Gecim, B.↗
Plasma contacting - An enabling technology
An experimental study of plasma contacting with an emphasis on the electron collection mode of this process is described. Results illustrating variations in plasma property profiles and potential differences that develop at hollow cathode plasma contactors are presented. A model of the electron collection plasma contacting process that is consistent with experimentally measured results is reviewed. The shortcomings of laboratory results as direct predictors of contactor performance in space and their usefulness in validating numerical models of the contacting process, that can be used to predict such performance, are discussed.
Williams, John D.↗
Analysis of electrical contact occurrences between rolling surface with application to elastohydrodynamic lubrication
A mathematical model enabling the relation of electrical contact occurrences between rolling surfaces to the mean thickness of an elastohydrodynamic lubricant film is discussed. A previous development by other authors is shown to contain an error invalidating their results. Some potentially useful approximations are developed for the case when the lubricant film is large with respect to the surface roughness. Complete solutions that are valid for all film thicknesses are not yet available.
Sidik, S. M.↗
Computer Simulation of the Forces Acting on a Submerged Polystyrene Probe as it Approaches the Succinonitrile Melt-Solid Interface
A Modeling approach to simulate both mesoscale and microscopic forces acting in a typical AFM experiment is presented. A mesoscale level interaction between the cantilever tip and the sample surface is primarily described by the balance of attractive Van der Waals and repulsive forces. Ultimately, the goal is to measure the forces between a particle and the crystal-melt interface. Two modes of AFM operation are considered in this paper - a stationary and a "tapping" one. The continuous mechanics approach to model tip-surface interaction is presented. At microscopic levels, tip contamination and details of tip-surface interaction are modeled using a molecular dynamics approach for the case of polystyrene - succinonitrile contact. Integration of the mesoscale model with a molecular dynamic model is discussed.
Bune, Andris V.↗