Engineering Papers⌕ Search

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.

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

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↗
NASA NTRS1995

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

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

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↗
NASA NTRS2015

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.↗
NASA NTRSDate not supplied

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↗
NASA NTRS1980

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

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.↗
DOE OSTI2025

The effect of ion pairing on speciation and transport in ion exchange membranes at varying hydration levels: A four-state model

Understanding ion pairing in ion exchange membranes (IEMs) is essential for advancing IEM applications in energy and environmental technologies. Here, this study introduces a four-state molecular dynamics model to quantify speciation and transport within Nafion-117, specifically examining the role of ion pairing in monovalent and divalent counterions (NaCl, Na 2 SO 4 , and MgSO 4 ). By analyzing radial distribution functions (RDFs) and molecular snapshots, we distinguish ion pairing modes and classify counterions into four states: condensed counterion, condensed ion pair, free ion pair, and free counterion. A key finding is that while divalent counterions (e. g., Mg 2+ ) maintain stable speciation across hydration levels, monovalent counterions (e.g., Na + ) show notable speciation shifts with hydration. Both monovalent and divalent counterions are not diffusive when condensed onto the polymer (sorbed to membrane functional groups). In contrast, free counterions are diffusive across all hydration levels. To evaluate the overall diffusivity of counterions, four-state fractions and diffusivities are computed, each contributing to counterion transport. The condensed/free ion speciation for multivalent sulfate salts aligns with previous revisions to the Donnan-Manning framework that include ion pairing, thereby validating its relevance to established membrane theories. The four-state model's diffusivity results support several current ion exchange assumptions, including that the condensed counterions are immobile, while uncondensed counterions are mobile. The four-state model offers insights into contact ion pairing within IEMs, highlighting its potential even when undetected in aqueous solution experiments. This work advances the theoretical understanding of counterion speciation in IEMs while identifying model limitations that suggest avenues for refinement, such as distinguishing water-mediated ion pairs between fully hydrated ions.

Ion exchange membranes↗
NASA NTRS1973

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.↗
DOE OSTI2021

Tool-workpiece stick-slip conditions and their effects on torque and heat generation rate in the friction stir welding

Friction stir welding (FSW) has found increased applications in automotive and aerospace industries due to its advantages of solid-state bonding, no fusion and melting, and versatility in various working conditions and material combinations. However, the relationship among processing parameters, material properties, and bonding extent and fidelity remains largely empirical, primarily because of the lack of the mechanistic understanding of the tool-workpiece frictional behavior that affects our subsequent understanding of microstructural evolution and interface bonding formation. While the tool-workpiece stick-slip condition is believed to dictate the resulting torque and heat generation rate during the welding process, it remains rare and elusive to conduct a quantitative experimental measurement of such interfacial field. On the other hand, numerical simulations based on Computational Fluid Dynamics (CFD) rely on ad hoc assumptions of interfacial pressure and shear-stress conditions, but predictions can only be validated via the medium- and far-range temperature field which is known to be insensitive to the interfacial frictional behavior. This work first presents a comparison among two CFD-based simulation methodologies and the Coupled Eulerian Lagrangian (CEL) model in finite element method, the last of which uses the Coulomb friction so that the stick-slip is naturally developed. Based on the Hill-Bower similarity relationship in the contact analysis, an analytical model is developed here to prove why a constant stick-slip fraction will be developed in the steady state, to correlate the stick-slip fraction to processing parameters such as the tool spin rate, and further to derive dimensionless functions for torque and heat-generation-rate predictions. Pros and cons of various numerical approaches in predicting stick-slip are discussed, and our analytical model has been found to agree well with our numerical simulation and literature experimental results. These analyses provide the critical strain-rate and temperature fields that are needed for the bonding analysis in our future work.

36 MATERIALS SCIENCE↗
NASA NTRS2001

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

TESS Lens-Bezel Assembly Modal Testing

The Transiting Exoplanet Survey Satellite (TESS) program, led by the Kavli Institute for Astrophysics and Space Research at the Massachusetts Institute of Technology (MIT) will be the first-ever spaceborne all-sky transit survey. MIT Lincoln Laboratory is responsible for the cameras, including the lens assemblies, detector assemblies, lens hoods, and camera mounts. TESS is scheduled to be launched in August of 2017 with the primary goal to detect small planets with bright host starts in the solar neighborhood, so that detailed characterizations of the planets and their atmospheres can be performed. The TESS payload consists of four identical cameras and a data handling unit. Each camera consists of a lens assembly with seven optical elements and a detector assembly with four charge-coupled devices (CCDs) including their associated electronics. The optical prescription requires that several of the lenses are in close proximity to a neighboring element. A finite element model (FEM) was developed to estimate the relative deflections between each lens-bezel assembly under launch loads to predict that there are adequate clearances preventing the lenses from making contact. Modal tests using non-contact response measurements were conducted to experimentally estimate the modal parameters of the lens-bezel assembly, and used to validate the initial FEM assumptions. Key Words Non-contact measurements, modal analysis, model validation

Modal Testing↗
NASA NTRSDate not supplied

Lunar Lander Landing Modeling and Simulation in Simscape

The landing gear is a critical subsystem of a lander for a successful landing on the moon or a planet. During a landing event, structural loading, energy absorption and toppling stability are important factors for the landing gear development. These aspects can be studied using simulations. The landing gear comprising the primary struct, secondary structs, footpad, and joints, is easily modeled using Simscape Multibody toolbox. The energy absorption mechanisms in the structs is readily modeled using Simulink and Multibody blocks based on the relative motions of the inner and outer cylinder of the structures. During touchdown, the footpads strike the lunar surface, and the contact forces can be modeled including the soil mechanical properties. Slosh dynamics are included using a mechanical pendulum model. Tip-over controllers are developed using Simulink blocks and toolboxes. This model was intensively tested manually then additionally with Monte Carlo simulation. It offers the opportunity for landing stability study, i.e. the max slope a lander can land without tipping over in 3D realistic landing situations.

Wei Wu↗
NASA NTRSDate not supplied

Impact of Thermoplastic Composites: Testing and Modeling

Thermoplastic Composites (TPCs) are being increasingly considered for aerospace applications given their faster manufacturing cycles and lower cost. The NASA High-Rate Composite Aircraft Manufacturing (HiCAM) project aims to evaluate and mature thermoplastic manufacturing technologies to achieve a 4X – 6X increase in the production rate of Next Generation Single Aisle (NGSA) commercial aircraft. Within the purview of the NASA HiCAM project, this work investigates the response of thermoplastic composites under High Energy Dynamic Impact (HEDI) conditions. HEDI tests were conducted on panels fabricated with a carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK) material system. A blunt metal projectile was chosen as the impactor which impacted the T700/LMPAEK panels over a range of impact velocities imparted by a gas gun test setup. As expected, lower velocity impacts caused the projectile to rebound, whereas higher velocity impacts resulted in the projectile penetrating the test panels. For the test cases in which the projectile rebounded, the predominant damage modes recorded with Ultrasonic (UT) scans included interlaminar delaminations which exhibited a “rotating fan” like structure when viewed in the through-thickness direction. Additionally, the dynamic deflection of the center point of the back face of the test panels and the velocity of the projectile was also measured during the test. Finite element models were developed to model a test case in which the projectile rebounded. These models were developed using LS-DYNA® wherein the primary objective was to evaluate the ability of the material model MAT299 to capture the dynamic response and damage modes in the T700/LMPAEK test panels. MAT299 is a Deformation Gradient Decomposition (DGD) based Continuum Damage Mechanics (CDM) material model. The interlaminar delaminations were modeled using the cohesive contact formulation available within LS-DYNA, wherein the mixed-mode fracture is captured via the Benzegaggh-Kenane (B-K) mode-dependent fracture energy interpolation law. The above-mentioned modeling methods have mostly been applied to thermoset composites and adapting them to TPCs involves addressing challenges associated with appropriately representing material behavior. Therefore, the proposed paper would discuss the systematic approach taken to adapt thermoset composites modeling practices for applications to TPCs, while appropriately addressing material behavior. The proposed paper would show that the predicted back face center point deflection correlated reasonably well with experimental data and the peak deflection was predicted to be within 5% of the experimental measurements. The projectile rebound velocity, while predicted to be higher than the experimental measurement, was within reasonable bounds. Additionally, it will be shown that the predicted delamination shapes were in reasonable agreement with experimental data.

Composite Materials↗
DOE OSTI2023

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Concurrent measurement of strain and chemical reaction rates in a calcite grain pack undergoing pressure solution: Evidence for surface-reaction controlled dissolution

Pressure solution is inferred to be a significant contributor to sediment compaction and lithification, especially in carbonate sediments. For a sediment deforming primarily by pressure solution, the compaction rate should be directly related to the rate of calcite dissolution, transport along grain contacts, and calcite reprecipitation. Previous experimental work has shown that there is evidence that deformation in wet calcite grain packs is consistent with control by pressure solution, but considerable ambiguity remains regarding the rate limiting mechanism. We present the results of laboratory compaction experiments designed to directly measure calcite dissolution and precipitation rates (recrystallization rates) concurrently with strain rate to test whether measured rates are consistent with predicted rates both in absolute magnitude and time evolution. Recrystallization rates are measured using trace element chemistry (Sr/Ca, Mg/Ca) and isotopes (87Sr/86Sr) of fluids flowing slowly through a compacting grain pack as it is being triaxially compressed. Imaging techniques are used to characterize the grain contacts and strain effects in the post-experiment grain pack. Our data show that calcite recrystallization rates calculated from all three geochemical parameters are in approximate agreement and that the rates closely track strain rate. The geochemically inferred rates are close to predicted rates in absolute magnitude. Uncertainty in grain contact dimensions makes distinguishing between surface reaction control and diffusion control difficult. Measured reaction rates decrease faster than predicted from standard pressure solution creep flow laws. This inconsistency may indicate that calcite dissolution rates at grain contacts are more complex, and more time-dependent, than suggested by geometric models designed to predict grain contact stresses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗