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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Delamination growth in composite materials

The Double Cantilever Beam (DCB) and the End Notched Flexure (ENF) specimens are employed to characterize MODE I and MODE II interlaminar fracture resistance of graphite/epoxy (CYCOM 982) and graphite/PEEK (APC2) composites. Sizing of test specimen geometries to achieve crack growth in the linear elastic regime is presented. Data reduction schemes based upon beam theory are derived for the ENF specimen and include the effects of shear deformation and friction between crack surfaces on compliance, C, and strain energy release rate, G sub II. Finite element (FE) analyses of the ENF geometry including the contact problem with friction are presented to assess the accuracy of beam theory expressions for C and G sub II. Virtual crack closure techniques verify that the ENF specimen is a pure Mode II test. Beam theory expressions are shown to be conservative by 20 to 40 percent for typical unidirectional test specimen geometries. A FE parametric study investigating the influence of delamination length and depth, span, thickness and material properties on G sub II is presented. Mode I and II interlaminar fracture test results are presented. Important experimental parameters are isolated, such as precracking techniques, rate effects, and nonlinear load-deflection response. It is found that subcritical crack growth and inelastic materials behavior, responsible for the observed nonlinearities, are highly rate-dependent phenomena with high rates generally leading to linear elastic response.

Gillespie, J. W., Jr.↗
NASA NTRS1995

Dynamic sensitivity analysis of frictional contact/impact response of axisymmetric composite structures

A computational procedure is presented for evaluating the sensitivity coefficients of the dynamic frictional contact/impact response of axisymmetric composite structures. The structures are assumed to consist of an arbitrary number of perfectly bonded homogeneous anisotropic layers. The material of each layer is assumed to be hyperelastic, and the effect of geometric nonlinearity is included. The sensitivity coefficients measure the sensitivity of the response to variations in different material, lamination and geometric parameters of the structure. A displacement finite element model is used for the discretization. The normal contact conditions are incorporated into the formulation by using a perturbed Lagrangian approach with the fundamental unknowns consisting of the nodal displacements, and the Lagrange multipliers associated with the contact conditions. The Lagrange multipliers are allowed to be discontinuous at interelement boundaries. Tangential contact conditions are incorporated by using a penalty method in conjunction with the classical Coulomb's friction model. Temporal integration is performed by using Newmark method. The Newton-Raphson iterative scheme is used for the solution of the resulting nonlinear algebraic equations, and for the determination of the contact region, contact conditions (sliding or sticking), and the contact pressures. The sensitivity coefficients are evaluated by using a direct differentiation approach. Numerical results are presented from the frictional contact/impact response of a composite spherical cap impacting on a rigid plate.

Karaoglan, Levent↗
NASA NTRS1989

Stress analysis method for clearance-fit joints with bearing-bypass loads

A simple direct stress analysis method is proposed for a laminate with a clearance-fit fastener subjected to combined bearing and bypass loads in tension or compression, including bearing and bypass loads on bolt-hole contact and local stresses. The approach uses a linear-elastic finite element analysis with an inverse formulation; conditions along the bolt-hole interface are specified by constraint equations that limit nodal displacements to a circular arc corresponding to the bolt diameter. The application of the method to the analysis of the effects of bearing-bypass loading on bolt-hole contact angles and local stresses, with the plate material properties representing a quasi-isotropic T300/5208 graphite/epoxy laminate, is discussed.

Naik, R. A.↗
NASA NTRS1990

Memory switches based on metal oxide thin films

MnO.sub.2-x thin films (12) exhibit irreversible memory switching (28) with an OFF/ON resistance ratio of at least about 10.sup.3 and the tailorability of ON state (20) resistance. Such films are potentially extremely useful as a connection element in a variety of microelectronic circuits and arrays (24). Such films provide a pre-tailored, finite, non-volatile resistive element at a desired place in an electric circuit, which can be electrically turned OFF (22) or disconnected as desired, by application of an electrical pulse. Microswitch structures (10) constitute the thin film element, contacted by a pair of separate electrodes (16a, 16b) and have a finite, pre-selected ON resistance which is ideally suited, for example, as a programmable binary synaptic connection for electronic implementation of neural network architectures. The MnO.sub.2-x microswitch is non-volatile, patternable, insensitive to ultraviolet light, and adherent to a variety of insulating substrates (14), such as glass and silicon dioxide-coated silicon substrates.

Ramesham, Rajeshuni↗
NASA NTRS1995

Thermal behavior spiral bevel gears

An experimental and analytical study of the thermal behavior of spiral bevel gears is presented. Experimental data were taken using thermocoupled test hardware and an infrared microscope. Many operational parameters were varied to investigate their effects on the thermal behavior. The data taken were also used to validate the boundary conditions applied to the analytical model. A finite element-based solution sequence was developed. The three-dimensional model was developed based on the manufacturing process for these gears. Contact between the meshing gears was found using tooth contact analysis to describe the location, curvatures, orientations, and surface velocities. This information was then used in a three-dimensional Hertzian contact analysis to predict contact ellipse size and maximum pressure. From these results, an estimate of the heat flux magnitude and the location on the finite element model was made. The finite element model used time-averaged boundary conditions to permit the solution to attain steady state in a computationally efficient manner.Then time- and position-varying boundary conditions were applied to the model to analyze the cyclic heating and cooling due to the gears meshing and transferring heat to the surroundings, respectively. The model was run in this mode until the temperature behavior stabilized. The transient flash temperature on the surface was therefore described. The analysis can be used to predict the overall expected thermal behavior of spiral bevel gears. The experimental and analytical results were compared for this study and also with a limited number of other studies. The experimental and analytical results attained in the current study were basically within 10% of each other for the cases compared. The experimental comparison was for bulk thermocouple locations and data taken with an infrared microscope. The results of a limited number of other studies were compared with those obtained herein and predicted the same basic behavior.

Handschuh, Robert F.↗
DOE OSTI2020

Quantifying the resistance to dislocation glide in single phase FeCrAl alloy

The resistance to dislocation glide associated with slip systems {110}<111> and {112}<111> in single phase FeCrAl alloy is measured via micromechanical testing in a scanning electronic microscopy at room temperature. Two important factors, the shape and orientation of a pillar, are discussed with respect to the glide resistance and stress-strain response. Maximizing Schmid factor of one specific slip system while minimizing others is recommended in order to diminish obvious dislocations-induced hardening during in-situ testing. Apparent Schmid factor analysis is conducted to select grains with preferred orientations. Two types of pillars with conventional cylindrical shape or dog-bone shape are tested under compression to estimate the resistance to dislocation glide and evaluate the effect of pillar shape on the compression stress-strain response. One dog-bone pillar is tested under tension to check the tension-compression isotropy of dislocation slip. We discover that the shape of a pillar to a smaller extent affects the measured resistance but strongly influences the stress strain response. Cylindrical pillars exhibit apparent hardening associated with early yielding due to stress or strain concentration at contact region, while dog-bone pillars show an obvious yielding and continuous shearing without hardening. The resistance is 220 MPa for slip system {110}<111> and 230 MPa for slip system {112}<111>. Finite element analysis is performed to account for the influence of pillar shape and contact condition on mechanical response.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Shape optimization including finite element grid adaptation

The prediction of optimal shape design for structures depends on having a sufficient level of precision in the computation of structural response. These requirements become critical in situations where the region to be designed includes stress concentrations or unilateral contact surfaces, for example. In the approach to shape optimization discussed here, a means to obtain grid adaptation is incorporated into the finite element procedures. This facility makes it possible to maintain a level of quality in the computational estimate of response that is surely adequate for the shape design problem.

Kikuchi, N.↗
NASA NTRS1991

A new procedure for calculating contact stresses in gear teeth

A numerical procedure for evaluating and monitoring contact stresses in meshing gear teeth is discussed. The procedure is intended to extend the range of applicability and to improve the accuracy of gear contact stress analysis. The procedure is based upon fundamental solution from the theory of elasticity. It is an iterative numerical procedure. The method is believed to have distinct advantages over the classical Hertz method, the finite-element method, and over existing approaches with the boundary element method. Unlike many classical contact stress analyses, friction effects and sliding are included. Slipping and sticking in the contact region are studied. Several examples are discussed. The results are in agreement with classical results. Applications are presented for spur gears.

Somprakit, Paisan↗
NASA NTRSDate not supplied

Developing Procedures to Implement Geometric Imperfections Beyond Right Circular Cylindrical Shells in Finite Element Method Models

Analysis of aerospace structures is frequently conducted using nominal dimensions and frequently assumes ideal conditions in loading, contact, constraints, et cetera. Off-nominal dimensions and nonideal conditions, however, are present in all structures. These are the result of widely ranging causes from coefficient of thermal expansion mismatches, manufacturing tooling anomalies, to assembly procedures that inadvertently alter the structure. Specifically, geometric imperfections can have potentially significant influence on the response of a structural test article observed in an experiment versus the response given by a numerical simulation. The Python Tool for Implementing Geometric Imperfections in Reduced Structures (Py_TIGIRS) was previously presented as a set of Python scripts to calculate and implement as-manufactured geometric midsurface and thickness imperfections into finite element method (FEM) shell models of nominally right circular cylinders. By taking advantage of the simple shape of a right circular cylinder, interpolations of the measured data points were able to be performed along directions that aligned to the cylindrical coordinate system axes of the entire structure. By taking advantage of the shell representation of the real structure as opposed to modeling using a continuum representation, the thickness variation was able to be implemented by shell section definitions instead of having to modify the position of multiple nodes in the thickness direction. Py_TIGIRS is a useful tool that established a procedural example on how to implement geometric imperfections in right circular cylindrical shell structures. Three new procedures, each expanded from concepts established in Py_TIGIRS, are proposed for various test-article designs and are intended to broaden the range of structures that can be modeled with measured geometric imperfections in the structural analysis community. Each test-article design introduces new challenges to successfully implement geometric imperfections into a FEM model. The first test-article design consists of a carbon fiber reinforced polymer square plate with a hat-shaped stiffener co-cured on one side. This test-article design was for a novel seven-point bend test that was also previously presented. Manufacturing and cure-cycle imperfections are observed using digital image correlation (DIC) techniques. As thermal expansion coefficient mismatches between the plate and stiffener materials were anticipated, a thermal analysis study with continuum shell and solid elements was conducted to capture the global shape observed prior to testing. The second test-article design is of a similar hat-stiffened plate configuration, but with a side length ratio near 3:1 with elongation in the stiffener direction. The test article was used to characterize the response to uniaxial compressive loading in the direction of the stiffener. Due to differing manufacturing steps, a thermal analysis like the one developed for the seven-point bend configuration was unable to mimic the observed geometric imperfections. Instead, a strategy based on applying deformations directly to the structure during analysis was developed for continuum shell and solid element representation of a stiffened panel.

Geometric imperfections↗
NASA NTRS1988

Advances in contact algorithms and their application to tires

Currently used techniques for tire contact analysis are reviewed. Discussion focuses on the different techniques used in modeling frictional forces and the treatment of contact conditions. A status report is presented on a new computational strategy for the modeling and analysis of tires, including the solution of the contact problem. The key elements of the proposed strategy are: (1) use of semianalytic mixed finite elements in which the shell variables are represented by Fourier series in the circumferential direction and piecewise polynomials in the meridional direction; (2) use of perturbed Lagrangian formulation for the determination of the contact area and pressure; and (3) application of multilevel iterative procedures and reduction techniques to generate the response of the tire. Numerical results are presented to demonstrate the effectiveness of a proposed procedure for generating the tire response associated with different Fourier harmonics.

Noor, Ahmed K.↗
DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2021

Up-Tower Investigation of Main Bearing Cage Slip and Loads

This paper described the operational conditions of an SRB in the main bearing position of a three-point mount wind turbine drivetrain. A commercial wind turbine was instrumented to measure load characteristics on the outer ring of the bearing, which were used to examine the roller load characteristics qualitatively and the cage slip quantitatively in a range of operational states. Further, a simple analytic model estimated the bearing loads and was used as input conditions to a three-dimensional finite element model. The finite element roller load results were compared to the measurements and contact stresses were also described.

17 WIND ENERGY↗
NASA NTRS1980

Thermal elastohydrodynamic lubrication of spur gears

An analysis and computer program called TELSGE were developed to predict the variations of dynamic load, surface temperature, and lubricant film thickness along the contacting path during the engagement of a pair of involute spur gears. The analysis of dynamic load includes the effect of gear inertia, the effect of load sharing of adjacent teeth, and the effect of variable tooth stiffness which are obtained by a finite-element method. Results obtained from TELSGE for the dynamic load distributions along the contacting path for various speeds of a pair of test gears show patterns similar to that observed experimentally. Effects of damping ratio, contact ratio, tip relief, and tooth error on the dynamic load were examined. In addition, two dimensionless charts are included for predicting the maximum equilibrium surface temperature, which can be used to estimate directly the lubricant film thickness based on well established EHD analysis.

Wang, K. L.↗
DOE OSTI2022

Investigation of Cutter–Woodchip Contact Pressure in a New Biomass Comminution System

A Crumbler ® rotary shear system was recently developed to process fuel-grade woodchips into feedstocks having more uniform particle sizes for improved flowability and higher thermochemical output yield than a traditional hammer mill. It represents a significant innovation in the state of the art of biomass comminution equipment. However, the milling unit has experienced significant abrasive wear of the feeding teeth and cutters when processing hard and dusty feedstocks, such as logging residue and corn stover. Here, we present initial simulation results from a finite-element stress analysis of the rotary shear, which is used to investigate the dynamics of interaction between the cutters and woodchips. This is the first step for optimizing the cutter design to improve the system lifetime and processing efficiency, thus reducing downtime and improving overall productivity.

42 ENGINEERING↗
DOE OSTI2025

Feasibility Assessment of Using Electrical Impedance Tomography for Damage Localization in Graphite Microreactor Components

Microreactors have great potential to decrease capital costs and construction timelines, reducing the barriers to implementing advanced nuclear reactor technologies. However, the lower power output of these microreactors introduces economic challenges if the operation and maintenance costs cannot be reduced sufficiently. Many microreactor concepts use graphite materials for in-core moderator and structural purposes, which will require periodic inspection or, ideally, in situ structural health monitoring. This work describes an initial evaluation of the feasibility of leveraging graphite’s semiconducting properties to perform electrical impedance tomography (EIT) for defect localization. First, a study was performed to identify the best methods for bonding electrical contacts to minimize the effects of contact resistance that interfere with impedance measurements of the graphite. After determining the best bonding approach, electrodes were attached to multiple graphite components with varying geometries (e.g., simple block and more representative microreactor hexagonal block). In parallel, finite element analysis approaches were developed and implemented to analyze the impact of defects and contact resistances and to inform an iterative inverse model for reconstructing the conductivity distribution. The results of this study show that defect localization in graphite components is possible using EIT if a sufficient number of electrodes (to improve spatial resolution) can be bonded using a technique with very low contact resistance and those contact resistances remain stable during reactor operation. If the reactor core and vessel design do not allow electrodes to be bonded during reactor operation, then it may be possible to detach/reattach electrodes between reactor operation and outages if a suitable mechanical connection that offers repeatable contact resistance can be identified. Alternatively, future work could focus on identifying contacts that can survive the harsh reactor operating environment so that contacts do not need to be removed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

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

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

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A numerical study of the meso-structure variability in the compaction process of prepreg platelet molded composites

Compaction deformations in a thermoplastic composite system molded from chopped prepreg platelets were analyzed to obtain statistical distributions of local distortions that define the final consolidated meso-structure. This study provides an insight into and allows quantification of the intrinsic geometric variability of the composite meso-structure that controls the variability of the effective composite mechanical properties. A two-phase viscoelastic finite element model was utilized to treat the behavior of molten prepreg platelets, which were assumed to be incompressible and inextensible in the fiber direction during processing. A two-step simulation procedure with discretely modeled platelets and platelet-to-platelet contacts is proposed. The analysis begins with a drop simulation of the uncontrolled material deposition followed by a 3D finite element analysis of the preform compaction and bulk material flow. The proposed model provides a method for generating a composite meso-structure ab initio starting from a statistical sampling of the initial conditions for platelet deposition. The predicted platelet thickness distribution, intra-platelet fiber waviness and out-of-plane platelet undulations were found to be in good agreement with experimental measurements. An estimate of the volume of resin pockets was inferred from the compaction simulation.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Tenth NASTRAN User's Colloquium

The development of the NASTRAN computer program, a general purpose finite element computer code for structural analysis, was discussed. The application and development of NASTRAN is presented in the following topics: improvements and enhancements; developments of pre and postprocessors; interactive review system; the use of harmonic expansions in magnetic field problems; improving a dynamic model with test data using Linwood; solution of axisymmetric fluid structure interaction problems; large displacements and stability analysis of nonlinear propeller structures; prediction of bead area contact load at the tire wheel interface; elastic plastic analysis of an overloaded breech ring; finite element solution of torsion and other 2-D Poisson equations; new capability for elastic aircraft airloads; usage of substructuring analysis in the get away special program; solving symmetric structures with nonsymmetric loads; evaluation and reduction of errors induced by Guyan transformation.

Source record↗