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.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Determination of near-surface plastic deformation in sliding contacts

It is pointed out that substantial plastic deformation occurs on and near the contact surfaces, when two solid bodies slide against each other without lubrication. It has been found that this deformation plays an important role in the tribological behavior of the sliding contact. The present investigation has the objective to develop an analytical model to predict the near-surface plastic deformation resulting from a single pass of one metallic surface over another. A finite element viscoplasticity program was written relating velocities to forces in a two-dimensional domain. The program was employed in the study of plastic deformation during a single pass of a hardened tool steel slider over a copper rub specimen. It was found that essentially the only material set in motion by the slider was directly under the contact zone. The agreement between values obtained in the analysis and experimental data is reasonably good.

Kennedy, F. E., Jr.↗

DOE OSTI2022

Removal of Non-Specifically Bound Proteins Using Rayleigh Waves Generated on ST-Quartz Substrates

Label-free biosensors are plagued by the issue of non-specific protein binding which negatively affects sensing parameters such as sensitivity, selectivity, and limit-of-detection. In the current work, we explore the possibility of using the Rayleigh waves in ST-Quartz devices to efficiently remove non-specifically bound proteins via acoustic streaming. A coupled-field finite element (FE) fluid structure interaction (FSI) model of a surface acoustic wave (SAW) device based on ST-Quartz substrate in contact with a liquid loading was first used to predict trends in forces related to SAW-induced acoustic streaming. Based on model predictions, it is found that the computed SAW body force is sufficient to overcome adhesive forces between particles and a surface while lift and drag forces prevent reattachment for a range of SAW frequencies. We further performed experiments to validate the model predictions and observe that the excitation of Rayleigh SAWs removed non-specifically bound (NSB) antigens and antibodies from sensing and non-sensing regions, while rinsing and blocking agents were ineffective. An amplified RF signal applied to the device input disrupted the specific interactions between antigens and their capture antibody as well. ST-quartz allows propagation of Rayleigh and leaky SH-SAW waves in orthogonal directions. Thus, the results reported here could allow integration of three important biosensor functions on a single chip, i.e., removal of non-specific binding, mixing, and sensing in the liquid phase.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1984

Numerical analysis of some problems related to the mechanics of pneumatic tires: Finite deformation/rolling contact of a viscoelastic cylinder and finite deformation of cord-reinforced rubber composites

The formulation and numerical analysis of several problems related to the behavior of pneumatic tires are considered. These problems include the general rolling contact problem of a rubber-like viscoelastic cylinder undergoing finite deformations and the finite deformation of cord-reinforced rubber composites. New finite element models are developed for these problems. Numerical results obtained for several representative cases are presented.

Oden, J. T.↗

NASA NTRS1987

Thermal-stress-free fasteners for joining orthotropic materials

Hot structures fabricated from orthotropic materials are an attractive design option for future high speed vehicles. Joining subassemblies of these materials with standard cylindrical fasteners can lead to loose joints or highly stressed joints due to thermal stress. A method has been developed to eliminate thermal stresses and maintain a tight joint by shaping the fastener and mating hole. This method allows both materials (fastener and structure), with different coefficients of thermal expansion (CTEs) in each of the three material directions, to expand freely with temperature yet remain in contact. For the assumptions made in the analysis, the joint will remain snug, yet free of thermal stress at any temperature. Finite element analysis was used to verify several thermal-stress-free fasteners and to show that conical fasteners, which are thermal-stress-free for isotropic materials, can reduce thermal stresses for transversely isotropic materials compared to a cylindrical fastener. Equations for thermal-stress-free shapes are presented and typical fastener shapes are shown.

Blosser, M. L.↗

NASA NTRS1987

Thermal-stress-free fasteners for joining orthotropic materials

Hot structures fabricated from orthotropic materials are an attractive design option for future high speed vehicles. Joining subassemblies of these materials with standard cylindrical fasteners can lead to loose joints or highly stressed joints due to thermal stress. A method has been developed to eliminate thermal stress and maintain a tight joint by shaping the fastener and mating hole. This method allows both materials (fastener and structure), with different coefficients of thermal expansion (CTE's) in each of the three principal material directions, to expand freely with temperature yet remain in contact. For the assumptions made in the analysis, the joint will remain snug, yet free of thermal stress at any temperature. Finite element analysis was used to verify several thermal-stress-free fasteners and to show that conical fasteners, which are thermal-stress-free for isotropic materials, can reduce thermal stresses for transversely isotropic materials compared to a cylindrical fastener. Equations for thermal-stress-free shapes are presented and typical fastener shapes are shown.

Blosser, Max L.↗

NASA NTRS2011

Ascent Heating Thermal Analysis on Spacecraft Adaptor Fairings

When the Crew Exploration Vehicle (CEV) is launched, the spacecraft adaptor (SA) fairings that cover the CEV service module (SM) are exposed to aero heating. Thermal analysis is performed to compute the fairing temperatures and to investigate whether the temperatures are within the material limits for nominal ascent aeroheating case. The ascent heating is analyzed by using computational fluid dynamics (CFD) and engineering codes at Marshall Space Flight Center. The aeroheating environment data used for this work is known as Thermal Environment 3 (TE3) heating data. One of the major concerns is with the SA fairings covering the CEV SM and the SM/crew launch vehicle (CLV) flange interface. The TE3 heating rate is a function of time, wall temperature, and the spatial locations. The implementation of the TE3 heating rate as boundary conditions in the thermal analysis becomes challenging. The ascent heating thermal analysis on SA fairings and SM/CLV flange interface are performed using two commercial software packages: Cullimore & Ring (C&R) Thermal Desktop (TD) 5.1 and MSC Patran 2007r1 b. TD is the pre-and post-processor for SINDA, which is a finite-difference-based solver. In TD, the geometry is built and meshed, the boundary conditions are defined, and then SINDA is used to compute temperatures. MSC Pthermal is a finite-element- based thermal solver. MSC Patran is the pre- and post-processor for Pthermal. Regarding the boundary conditions, the convection, contact resistance, and heat load can be imposed in different ways in both programs. These two software packages are used to build the thermal model for the same analysis to validate each other and show the differences in the modeling details.

Wang, Xiao Yen↗

DOE OSTI2025

Influence of Initial fabric and water-wetting on particle fracture and force chain evolution in natural silica sand

Particle fracture has a significant influence on the engineering behavior of granular materials. However, the combined influence of the initial fabric and wetting on particle fracture and force chain evolution in granular materials, particularly in sands, remains insufficiently explored in the current literature. Here, in this paper, one-dimensional (1D) confined compression experiments were conducted on dry and wet specimens, and particle-level fracture was captured using three-dimensional (3D) in-situ Synchrotron Micro-Computed Tomography (SMT). A quantitative assessment of particle fracture revealed that wet specimens, which have disturbed fabric, exhibited a higher fracture percentage in comparison to dry specimens. Complementary 3D finite element (FE) simulations were performed using dry and wet (only solid fabric was considered and the influence of pore water is not considered) sand assemblies to assess interparticle contact forces and particle-scale stress distributions within the specimens. Force chain analysis was subsequently conducted, encompassing quantification of interparticle forces, characterization of force network structures, and monitoring the dynamic evolution of force chains under different strain levels. The results show that the specimens with disturbed fabric led to a more dynamic and less persistent force network, more fabric instability, and thus more reorganization of force chain structures. In addition, the frequent rearrangement of the force network in the presence of water (with reduced inter-particle contact friction) likely exacerbates localized stress concentrations, promoting failure in previously unengaged or weakly connected particles. The results reported in this paper offer a new insight into how the initial fabric and wetting cause different fracture behavior. The findings can also pave the way for more in-depth future investigations into the mechanics governing particle fracture in granular assemblies.

Finite element analysis↗

NASA NTRS2023

A Finite Element Approach for Simplified 2D Nonlinear Dynamic Contact/Impact Analysis

In this paper, a simplified numerical approach for finite element dynamic analysis of an inelastic solid structure subjected to solid object impact is presented. The approach approximates the impacting solid as the selected multiple nodes, for which mass of the impactor is distributed. The node-to-segment contact formulation with the penalty constraint technique incorporated is employed to impose contact conditions between the nodes and the surface of the receiver structure. The node-to-segment algorithm is integrated into Newton–Raphson time integration scheme and the Lagrange multiplier technique is applied to enforce the identical displacements for the selected nodes throughout the analysis process. The approach is verified using two-dimensional plane strain models considering elastic-perfectly-plastic material behavior. The results obtained using the proposed approach are in a good agreement with those simulated using a commercial finite element code, ABAQUS dynamic/implicit, in terms of displacements and stress distribution fields. The proposed approach is shown to be computationally superior to general finite element method-based contact/impact analysis without significantly sacrificing the accuracy.

Seungwook Seok↗

DOE OSTI2022

Validation and Parametric Investigations of an Internal Permanent Magnet Motor Using a Lumped Parameter Thermal Model

One of the key challenges for the electric vehicle industry is to develop high-power-density electric motors. Achieving higher power density requires efficient heat removal from inside the motor. In order to improve thermal management, a multiphysics modeling framework that is able to accurately predict the behavior of the motor, while being computationally efficient, is essential. This paper first presents a detailed validation of a lumped parameter thermal network (LPTN) model of an Internal Permanent Magnet synchronous motor within the commercially available MOTOR-CAD modeling environment. The validation is based on temperature comparison with experimental data and with more detailed finite element analysis (FEA). All critical input parameters of the LPTN are considered in detail for each layer of the stator, especially the contact resistances between the impregnation, liner, laminations, and housing. Finally, a sensitivity analysis for each of the critical input parameters is provided. A maximum difference of 4% - for the highest temperature in the slot-winding and the end-winding - was found between the LPTN and the experimental data. Comparing the results from the LPTN and the FEA model, the maximum difference was 2% for the highest temperature in the slot-winding and end-winding. As for the LPTN sensitivity analysis, the thermal parameter with the highest sensitivity was found to be the liner-to-lamination contact resistance.

42 ENGINEERING↗

DOE OSTI2020

Validation and Parametric Investigations Using a Lumped Thermal Parameter Model of an Internal Permanent Magnet Motor

One of the key challenges for the electric vehicle industry is to develop high-power-density electric motors. Achieving higher power density requires efficient heat removal from inside the motor. In order to improve thermal management, a multi-physics modeling framework that is able to accurately predict the behavior of the motor, while being computationally efficient, is essential. This paper first presents a detailed validation of a Lumped Parameter Thermal Network (LPTN) model of an Internal Permanent Magnet (IPM) synchronous motor within the commercially available Motor-CAD® modeling environment. The IPM motor’s stator is studied at steady state, and winding losses are generated by a constant DC current. The validation is based on temperature comparison with experimental data and with more detailed Finite Element Analysis (FEA). All critical input parameters of the LPTN are considered in detail for each layer of the stator, especially the contact resistances between the impregnation, liner, laminations and housing. Finally, a sensitivity analysis for each of the critical input parameters is provided. A maximum difference of 4% — for the highest temperature in the slot windings and the end windings — was found between the LPTN and the experimental data. Comparing the results from the LPTN and the FEA model, the maximum difference was 2% for the highest temperature in the slot windings and end windings. As for the LTPN sensitivity analysis, the thermal parameter with the highest sensitivity was found to be the liner-to-lamination contact resistance. The latter is often ignored in the literature, whereas its impact on temperature rise was found to be more significant than any other contact resistance within the stator.

47 OTHER INSTRUMENTATION↗

NASA NTRS1974

A finite element analysis of fatigue crack closure

Experiments have shown that fatigue cracks close at positive loads during constant-amplitude load cycling. The crack-closure phenomenon is caused by residual plastic deformations remaining in the wake of an advancing crack tip. The present paper is concerned with the application of a two-dimensional, nonlinear, finite-element analysis for predicting crack-closure and crack-opening stresses during cyclic loading. A two-dimensional finite-element computer program, which accounts for both elastic-plastic material behavior and changing boundary conditions associated with crack extension and intermittent contact of the crack surfaces under cyclic loading, has been developed. An efficient technique to account for changing boundary conditions was also incorporated into the nonlinear analysis program. This program was subsequently used to study crack extension and crack closure under constant-amplitude and two-level block loading. The calculated crack-closure and crack-opening stresses were qualitatively consistent with experimental observations.

Newman, J. C., Jr.↗

NASA NTRS2022

Design of Rigid-Flex PCB Robotics Leveraging Validated Finite Element Simulations

The use of rigid-flex printed circuit board (PCB) as primary structure has the potential to reduce the weight and volume of robotic systems. In the case of robotics for interplanetary exploration, these systems can leverage origami-inspired folding for increased mobility options and reduced storage volume. Folding rigid-flex PCB robotics can be constructed with rigid PCB connected by short Nomex fabric hinges coupled with flex PCB ribbon cables that permits enhanced system flexibility and energy dissipation to promote impact survivability. This paper presents a design methodology of rigid-flex PCB systems with an emphasis on impact resistance. The design process considers solder joint adequacy, panel bending, and fracture using a finite element (FE) model. The proposed design methodology is developed using a case study with NASA JPL's Pop-Up Folding Flat Explorer Robot (PUFFER). First, the finite-element (FE) modeling methodology is presented with consideration to both frequency and time-domain modeling applications, which include operational self-contact analysis and high impact scenarios. The time-domain impact modeling methodology utilizes hyperelastic material properties for the Nomex hinges. This modeling method is validated using image correlation of PUFFER drop tests. A flowchart is presented to guide users through a validated Abaqus modeling procedure for highly flexible rigid-flex systems. Next, a case study is presented in which PUFFER is subject to drop heights representative of falls into Lunar pits and then the design is refined for a more optimum impact performance. Finally, the results of the case study are used to inform a generalized design methodology for rigid-flex PCB robotics subject to high-impact loads with the considerations presented.

John Bell↗

NASA NTRS2022

Design of Rigid-Flex PCB Robotics Leveraging Validated Finite Element Simulations

The use of rigid-flex printed circuit board (PCB) as primary structure has the potential to reduce the weight and volume of robotic systems. In the case of robotics for interplanetary exploration, these systems can leverage origamiinspired folding for increased mobility options and reduced storage volume. Folding rigid-flex PCB robotics can be constructed with rigid PCB connected by short Nomex fabric hinges coupled with flex PCB ribbon cables that permits enhanced system flexibility and energy dissipation to promote impact survivability. This paper presents a design methodology of rigidflex PCB systems with an emphasis on impact resistance. The design process considers solder joint adequacy, panel bending, and fracture using a finite element (FE) model. The proposed design methodology is developed using a case study with NASA JPL’s Pop-Up Folding Flat Explorer Robot (PUFFER). First, the finite-element (FE) modeling methodology is presented with consideration to both frequency and time-domain modeling applications, which include operational self-contact analysis and high impact scenarios. The time-domain impact modeling methodology utilizes hyperelastic material properties for the Nomex hinges. This modeling method is validated using image correlation of PUFFER drop tests. A flowchart is presented to guide users through a validated Abaqus modeling procedure for highly flexible rigid-flex systems. Next, a case study is presented in which PUFFER is subject to drop heights representative of falls into Lunar pits and then the design is refined for a more optimum impact performance. Finally, the results of the case study are used to inform a generalized design methodology for rigid-flex PCB robotics subject to high-impact loads with the considerations presented.

de la Croix, Jean-Pierre↗

DOE OSTI2020

Validation and Parametric Investigations Using a Lumped Thermal Parameter Model of an Internal Permanent Magnet Motor: Preprint

One of the key challenges for the electric vehicle industry is to develop high-power-density electric motors. Achieving higher power density requires efficient heat removal from inside the motor. In order to improve thermal management, a multi-physics modeling framework that is able to accurately predict the behavior of the motor, while being computationally efficient, is essential. This paper first presents a detailed validation of a Lumped Parameter Thermal Network (LPTN) model of an Internal Permanent Magnet (IPM) synchronous motor within the commercially available Motor-CAD® modeling environment. The IPM motor’s stator is studied at steady state, and winding losses are generated by a constant DC current. The validation is based on temperature comparison with experimental data and with more detailed Finite Element Analysis (FEA). All critical input parameters of the LPTN are considered in detail for each layer of the stator, especially the contact resistances between the impregnation, liner, laminations and housing. Finally, a sensitivity analysis for each of the critical input parameters is provided. A maximum difference of 4% - for the highest temperature in the slot windings and the end windings - was found between the LPTN and the experimental data. Comparing the results from the LPTN and the FEA model, the maximum difference was 2% for the highest temperature in the slot windings and end windings. As for the LTPN sensitivity analysis, the thermal parameter with the highest sensitivity was found to be the liner-to-lamination contact resistance. The latter is often ignored in the literature, whereas its impact on temperature rise was found to be more significant than any other contact resistance within the stator.

ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗

NASA NTRS1996

Design of a Mode Conversion Ultrasonic Motor for Position Control

The many useful characteristics of ultrasonic motors, such as high holding torques, and high torque at low speeds, have made them the subject of increasing interest. In addition, several of their characteristics make them attractive for aerospace applications: they have a torque to weight ratio, and they require neither gearing mechanisms nor lubrication. Moreover, they create negligible magnetic fields, and conversely, they are not affected by external magnetic fields. Ultrasonic motors based on bolt-tightened structures offer simplicity and high stress capability. They use the inverse piezoelectric effect in the stator to produce vibrational energy, which is transferred to the rotor by friction. We designed a bolt-tightened ultrasonic motor using numerical modelling tools (finite element and electromechanical circuit analyses), creating an equivalent circuit model that takes into account the electromechanical energy conversion in the stator and the contact between the stator and the rotor. Analysis of the circuit gives insight into the behavior of the motor and allows its performance to be calculated. Two prototypes of the motor were built; their transient responses and other quantities, such as starting torque, were measured. In this paper, we discuss the numerical and the experimental results, and demonstrate the usefulness of numerical analysis in designing ultrasonic motors and estimating their performance.

LeLetty, Ronan↗

DOE OSTI2025

Finite element analysis of the impact of beam heating mode in molten salt corrosion experiments employing simultaneous ion irradiation

Finite element analysis was used to investigate the temperature and stress profiles that develop in 316L stainless steel membranes being irradiated using different proton beam conditions in contact with a molten salt environment. It was shown that in addition to a nonuniform irradiation profile, a focused 2 MeV proton beam leads to very strong temperature and stress gradients in the membrane, introducing highly localized driving forces that complicate and even compromise the integrity and reliability of the experimental results of corrosion studies. Here, the use of a focused beam in corrosion studies can create experimental artifacts that may misrepresent the true corrosion behavior. In contrast, the use of a rastered beam is shown to distribute the protons and resulting radiation damage uniformly across the membrane face, and more importantly, results in temperature and stress profiles that are not only very uniform but are of much lower magnitude. The use of a rastered beam during molten salt corrosion experiments is therefore recommended to achieve uniform damage rates, thereby reducing both gradients and magnitudes of the temperature and stress distributions.

43 - PARTICLE ACCELERATORS↗

DOE OSTI2020

Ultrasonic waves from radial mode excitation of a piezoelectric disc on the surface of an elastic solid

In this work, we discuss the ultrasonic wave propagation characteristics in an isotropic elastic solid due to the radial mode excitation of a piezoelectric disc actuator bonded to its surface. Finite element simulations using coupled electromechanical modeling is employed to investigate the wave propagation behavior. We find that the radial mode vibrations on the surface of the elastic solid generates all the three types of ultrasonic waves: longitudinal, shear, and surface waves. The waves in the solid are comprised of a central lobe and multiple side lobes based on the frequency of the radial mode excitation. The central lobe is predominantly composed of longitudinal waves and the side lobes are composed of shear waves. While the longitudinal waves have a strong central lobe and weak side lobes that are not fully developed, shear waves consists of fully developed side lobes. In addition, we observe that the longitudinal waves have fewer side lobes inside the solid compared to the shear waves. A semi-analytical approach is presented to explain the above observations. The interface traction boundary condition between the piezoelectric disc and the elastic solid is approximated as a truncated 'Bessel' excitation over the area of the piezoelectric disc in contact with the solid. A Fourier–Bessel series expansion technique is used to investigate the wave propagation behavior from the above traction boundary condition. The results obtained explain the observations from the finite element simulations with regard to the number and distribution of side lobes pertaining to the longitudinal and shear waves generated from the radial modes of the piezoelectric disc. The proposed semi-analytical approach is general and can be applied to any arbitrary axisymmetric excitations on the surface of an elastic solid.

36 MATERIALS SCIENCE↗