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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

An Anisotropic Yield and Damage Material Model to Improve the Contact Pressure Analysis in a Biomass Shredding System

Size reduction systems used in biomass processing break biomass into smaller pieces by utilizing the kinetic energy from the sharp rotating blades. Abrasive and/or erosive wear caused by biomass comminution results in blade wear of the sharp edged cutters, deteriorating the process efficiency. Here, this study aims to optimize the blade design and improve the system efficiency by attempting to understand the interactions between the blades and biomass particles. Since real-time monitoring of these interactions is impractical during operation, mechanical simulations offer a viable alternative for investigating the shredding process. Yet, the irregular geometry and complex mechanical properties of biomass—such as the anisotropic nature of woodchips and their nonlinear fracture behavior—pose significant challenges for accurately simulating contact pressure. In this work an anisotropic yield material model, along with a damage initiation and evolution function, is applied to the woodchip particle to study the contact pressure on shredder blade, offering a scientific basis for improved blade design and process efficiency. This approach can be extended to other biomass processing systems with similar anisotropic feedstocks, making it a valuable tool for advancing sustainable biomass utilization.

09 - BIOMASS FUELS↗
NASA NTRS1990

Simulation of crack growth and crack closure under large cyclic plasticity

The propagation of a crack during cyclic loading with large plasticity was simulated using a finite element method. The numerical procedure included the determination of crack closing and opening by examining the contact condition of the crack surfaces in the plasticity iteration. The method was applied to the growth of cracks in Alloy 718 single edge notch specimens at 538 C subject to strain cycling under both nominally elastic and elastoplastic conditions. The results show, among other things, that the crack closure and opening stresses become more compressive as the loading becomes more plastic, and that, under a given control condition, the closure and opening stresses become less compressive as the crack propagates.

Kim, K. S.↗
NASA NTRS1993

A study of subsurface crack initiation produced by rolling contact fatigue

Results of subsurface crack initiation studies produced by pure rolling contact fatigue in 7075-T6 aluminum alloy are presented in this article. Microstructural changes and subsequent crack initiation below the contacting surface in cylindrical test specimens subjected to repeated rolling contact are illustrated. The rolling conditions are simulated in a three-dimensional elastic-plastic finite element model in order to estimate the plastic strains and residual stresses in the test material. The numerically estimated distribution of plastic strains in the model correlate well with the extent of microstructural changes observed in the test specimen. Results also indicate that a combination of plastic strains and low values of residual stresses is conducive to subsurface crack initiation and growth.

Kumar, Arun M.↗
NASA NTRS2007

Investigation of Low-Cycle Bending Fatigue of AISI 9310 Steel Spur Gears

An investigation of the low-cycle bending fatigue of spur gears made from AISI 9310 gear steel was completed. Tests were conducted using the single-tooth bending method to achieve crack initiation and propagation. Tests were conducted on spur gears in a fatigue test machine using a dedicated gear test fixture. Test loads were applied at the highest point of single tooth contact. Gear bending stresses for a given testing load were calculated using a linear-elastic finite element model. Test data were accumulated from 1/4 cycle to several thousand cycles depending on the test stress level. The relationship of stress and cycles for crack initiation was found to be semilogarithmic. The relationship of stress and cycles for crack propagation was found to be linear. For the range of loads investigated, the crack propagation phase is related to the level of load being applied. Very high loads have comparable crack initiation and propagation times whereas lower loads can have a much smaller number of cycles for crack propagation cycles as compared to crack initiation.

Handschuh, Robert F.↗
NASA NTRS2007

Investigation of Low-Cycle Bending Fatigue of AISI 9310 Steel Spur Gears

An investigation of the low-cycle bending fatigue of spur gears made from AISI 9310 gear steel was completed. Tests were conducted using the single-tooth bending method to achieve crack initiation and propagation. Tests were conducted on spur gears in a fatigue test machine using a dedicated gear test fixture. Test loads were applied at the highest point of single tooth contact. Gear bending stresses for a given testing load were calculated using a linear-elastic finite element model. Test data were accumulated from 1/4 cycle to several thousand cycles depending on the test stress level. The relationship of stress and cycles for crack initiation was found to be semi-logarithmic. The relationship of stress and cycles for crack propagation was found to be linear. For the range of loads investigated, the crack propagation phase is related to the level of load being applied. Very high loads have comparable crack initiation and propagation times whereas lower loads can have a much smaller number of cycles for crack propagation cycles as compared to crack initiation.

Handschuh, Robert F.↗
DOE OSTI2025

Sensitivity analysis of thermal contact conductance modeling to inform MiniFuel irradiation capsule designs

The MiniFuel irradiation platform has been developed by Oak Ridge National Laboratory as a flexible, high-throughput separate effects testing capability within the High Flux Isotope Reactor (HFIR). Finite element thermal models are relied upon to design MiniFuel experiments to achieve a specific time-averaged irradiation temperature for experimental objectives. A previous study identified that uncertainty in the component heat generation rates and thermal contact conductance (TCC) model are the most significant contributors to predicted fuel temperature variance. To address both sources of uncertainty, this work performs sensitivity analysis on the TCC model to identify high-impact, high-uncertainty parameters that contribute to fuel temperature variance. The TCC model is analyzed in increasing detail, first using a standalone Python code, then again after coupling Python to the BISON fuel performance code. Furthermore, the parameters with the largest contributions to fuel temperature variance which can be reduced through design changes are identified as the initial subcapsule gas pressure, contact pressure between the fuel and dish, and the effective surface roughness of the interface. A set of design recommendations for future capsule designs has been established and applied to reduce the previously quantified average fuel temperature uncertainty ranges of ± 40 °C in the HFIR vertical experiment facilities (VXF) and ± 80 °C in the removable beryllium (RB) reflector to approximately ± 32 °C and ± 53 °C, respectively. This equates to a 21 % and 33 % reduction in the uncertainty range of the average fuel temperature for VXF and RB, respectively.

BISON↗
DOE OSTI2025

Predictive dynamic wetting, fluid–structure interaction simulations for braze run-out

Brazing and soldering are metallurgical joining techniques that use a wetting molten metal to create a joint between two faying surfaces. Here, the quality of the brazing process depends strongly on the wetting properties of the molten filler metal, namely the surface tension and contact angle, and the resulting joint can be susceptible to various defects, such as run-out and underfill, if the material properties or joining conditions are not suitable. In this work, we implement a finite element simulation to predict the formation of such defects in braze processes. This model incorporates both fluid–structure interaction through an arbitrary Eulerian–Lagrangian technique and free surface wetting through conformal decomposition finite element modeling. Upon validating our numerical simulations against experimental run-out studies on a silver-Kovar system, we then use the model to predict run-out and underfill in systems with variable surface tension, contact angles, and applied pressure. Finally, we consider variable joint/surface geometries and show how different geometrical configurations can help to mitigate run-out. This work aims to understand how brazing defects arise and validate a coupled wetting and fluid–structure interaction simulation that can be used for other industrial problems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

ECAR-7210 Rev 0 Steady-state Thermo-mechanical Analysis of the MARVEL Fuel Bundle

This document reports the steady-state thermal-hydraulic and static mechanical analyses of the MARVEL micro-reactor core. The objective of this work is to investigate whether fuel-to-fuel contact could occur due to the thermal deformation of the fuel cladding during the normal operation and to assess the resulting impact on the peak cladding temperatures (PCT) if contact occurs. The scope of the work consists of three computational tasks. First, a conjugate heat transfer analysis using computational fluid dynamics (CFD) model was performed to evaluate the PCTs during normal operating conditions. Second, a finite element analysis (FEA) using the cladding temperature field obtained from the CFD analysis was performed to evaluate the thermal deformation of cladding and determine the occurrence of fuel-to-fuel contact. Finally, a CFD analysis of FEA-informed fuel-to-fuel gap (represented by a conservative 0.05 mm uniform gap) was performed to evaluate whether the PCT exceeds the safety criteria.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
NASA NTRS1984

Theoretical model of impact damage in structural ceramics

This paper presents a mechanistically consistent model of impact damage based on elastic failures due to tensile and shear overloading. An elastic axisymmetric finite element model is used to determine the dynamic stresses generated by a single particle impact. Local failures in a finite element are assumed to occur when the primary/secondary principal stresses or the maximum shear stress reach critical tensile or shear stresses, respectively. The succession of failed elements thus models macrocrack growth. Sliding motions of cracks, which closed during unloading, are resisted by friction and the unrecovered deformation represents the 'plastic deformation' reported in the literature. The predicted ring cracks on the contact surface, as well as the cone cracks, median cracks, radial cracks, lateral cracks, and damage-induced porous zones in the interior of hot-pressed silicon nitride plates, matched those observed experimentally. The finite element model also predicted the uplifting of the free surface surrounding the impact site.

Liaw, B. M.↗
DOE OSTI2021

THERMAL/STRUCTURAL ANALYSIS OF THE AXION QUANTUM METROLOGY CAVITY AND ITS COMPONENTS

This research was centered around maximizing the capability to cool dielectric material within a containment unit, or Photonic Band Gap (PBG) cavity, designed for detecting axion dark matter and identifying the unit’s thermal properties. There are multiple types of PBG cavities, but the latest version that axion researchers wish to use has been theorized to contain possible issues related to its thermal properties. Thermal conductivity is an issue with the dielectric material because it is made from alumina which is highly insulative. This is important since the research is being done in a cryogenic environment and the thermal noise affects the quantum bit used for detecting the axion to photon conversion process. Therefore, any improvements to this unit should be justified and implemented but are not entirely limited to thermal contact related aspects of PBG cavities. A prospect of using sapphire in place of alumina also exists, but this is a more expensive, less tested, and more elusive material to justify the creation of a full dielectric structure out of sapphire. Thermal aspects of the cavity were analyzed using finite element method (FEM) and an experiment designed to test different thermal joint materials. FEM was used to check the contraction of the cavity during cooling, the contact quality between the dielectric material and cavity wall, theheat flow rate through the assembly, and helps visualize the cavity’s reaction to different design changes. A simple comparison between thermal conductivity curves justified the usage of sapphire over alumina from a thermal property point of view. Thus, the method of experimentation had an additional dielectric material, sapphire, to test with and compare to alumina. The thermal test identified which material is best to use for a thermal joint but also simultaneously found the conductance of the joint and dielectric material as well as an estimate for what the temperature is inside the larger scale cavity since a thermometer cannot be placed inside the PBG when it cools. Initially, the cavity simulation was tested for structural deformation properties and stress distribution since it is made of copper which contracts heavily in comparison to most other materials. The reduction in volume of the cavity was less than anticipated which gave more room for a possible usage of a modified thermal joint for additional contact area and higher conductivity.

43 PARTICLE ACCELERATORS↗
NASA NTRS1986

Plasticity analysis of wear phenomena as an aid in the development of abradable materials

Plastic deformation accompanying the wear process in fully dense and porous metals was studied in order to develop more abradable gas path seal material systems for turbine engine applications. The finite element method was used to model high-rate viscoplastic deformation in the vicinity of a moving contact. Techniques were developed for studying such deformation in both solid (fully dense) and porous metals. Results of the analysis were verified by comparison with experimental data obtained using a pendulum-type rub test apparatus and a pin-on-disk sliding wear tester. It was found that the near-surface plastic deformation in fully dense solids was different from that for porous materials because of the compressibility of the latter. Both types of materials exhibit deformation features which can be taken advantage of in the design of more abradable material systems.

Kennedy, F. E., Jr.↗
NASA NTRS1989

Bearing-bypass loading in composite joints - Testing and stress analysis

A combined experimental and analytical study of bearing and bypass loading on single-fastener specimens of a 16-ply, quasi-isotropic T300/5208 graphite/epoxy laminate with a centrally located hole is reported. The specimens were loaded in either tension or compression, and onset damage, ultimate strengths, and corresponding failure modes were determined. The tension data showed the expected linear interaction for combined bearing-bypass loading with damage developing in the net-section tension mode. However, the bearing-onset strengths showed an unexpected interaction of the bearing and compressive bypass loads in which the latter reduced the bearing-onset strength. A linear finite element analysis showed that bearing-bypass loading had a marked influence on the bolt-hole contact which in turn had a significant effect on local stresses.

Crews, John H., Jr.↗
NASA NTRS1997

Profiled Roller Stress/Fatigue Life Analysis Methodology and Establishment of an Appropriate Stress/Life Exponent

The objective of this work was to determine the three dimensional volumetric stress field, surface pressure distribution and actual contact area between a 0.50" square roller with different crown profiles and a flat raceway surface using Finite Element Analysis. The 3-dimensional stress field data was used in conjunction with several bearing fatigue life theories to extract appropriate values for stress-life exponents. Also, results of the FEA runs were used to evaluate the laminated roller model presently used for stress and life prediction.

Source record↗
NASA NTRS1999

Delamination and Stitched Failure in Stitched Composite Joints

The effect of stitches on the failure of a single lap joint configuration was determined in a combined experimental and finite element study. The experimental program was conducted to determine debond growth under static monotonic loading. The stitches were shown to delay the initiation of the debond and provide load transfer beyond the load necessary to completely debond the stitched lap joint. The experimentally determined debond length vs. applied load was used as an input parameter in the finite element analysis of both configurations. The strain energy release rates at the debond from were calculated using plate finite elements. Nonlinear fastener elements were used to model the stitches and multipoint constraints were used to model the contact problem. Models of the unstitched configuration showed significant values of modes I and II across the width of the joint and showed that mode III is zero at the centerline but increases near the free edge. Models of the stitched configuration showed that the stitches were effective in reducing mode I to zero, but had less of an effect on modes II and III.

Glaessgen, E. H.↗
NASA NTRS2019

Modeling and simulation challenges for softgoods in space deployable structures

This paper examines modeling and simulation challenges for deployable spacecraft structures with softgood components. Previous studies highlighted the numerical challenges due to the mixture of high stiffness and low stiffness contact in this class of structures. This paper examines the use of the commercial LS-DYNA nonlinear finite element code on a set of relevant benchmark problems.

Mobrem, Mehran↗
NASA NTRS2003

Finite-Element Methods for Real-Time Simulation of Surgery

Two finite-element methods have been developed for mathematical modeling of the time-dependent behaviors of deformable objects and, more specifically, the mechanical responses of soft tissues and organs in contact with surgical tools. These methods may afford the computational efficiency needed to satisfy the requirement to obtain computational results in real time for simulating surgical procedures as described in Simulation System for Training in Laparoscopic Surgery (NPO-21192) on page 31 in this issue of NASA Tech Briefs. Simulation of the behavior of soft tissue in real time is a challenging problem because of the complexity of soft-tissue mechanics. The responses of soft tissues are characterized by nonlinearities and by spatial inhomogeneities and rate and time dependences of material properties. Finite-element methods seem promising for integrating these characteristics of tissues into computational models of organs, but they demand much central-processing-unit (CPU) time and memory, and the demand increases with the number of nodes and degrees of freedom in a given finite-element model. Hence, as finite-element models become more realistic, it becomes more difficult to compute solutions in real time. In both of the present methods, one uses approximate mathematical models trading some accuracy for computational efficiency and thereby increasing the feasibility of attaining real-time up36 NASA Tech Briefs, October 2003 date rates. The first of these methods is based on modal analysis. In this method, one reduces the number of differential equations by selecting only the most significant vibration modes of an object (typically, a suitable number of the lowest-frequency modes) for computing deformations of the object in response to applied forces.

Basdogan, Cagatay↗
NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗
DOE OSTI2024

Adaptively remeshed multiphysical modeling of resistance forge welding with experimental validation of residual stress fields and measurement processes

Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.

36 MATERIALS SCIENCE↗