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 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1976

A finite-element analysis of fatigue crack closure

The paper is concerned with the application of a two-dimensional, nonlinear, finite-element analysis using an incremental theory of plasticity to predict crack-closure and crack-opening stresses during the crack-growth process under cyclic loading. A two-dimensional finite-element computer program is developed, 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. An efficient technique to account for changing boundary conditions under cyclic loading is also incorporated into the nonlinear analysis program. This program is used to study crack extension and crack closure behavior in a center-cracked panel under constant-amplitude and two-level block loading. The calculated crack-opening stresses are found to be quantitatively consistent with experimental measurements.

Newman, J. C., Jr.↗

NASA NTRS1975

Finite element solution of low bond number sloshing

The dynamics of liquid propellant in a low Bond number environment which are critical to the design of spacecraft systems with respect to orbital propellant transfer and attitude control system were investigated. Digital computer programs were developed for the determination of liquid free surface equilibrium shape, lateral slosh natural vibration mode shapes, and frequencies for a liquid in a container of arbitrary axisymmetric shape with surface tension forces the same order of magnitude as acceleration forces. A finite volume element representation of the liquid was used for the vibration analysis. The liquid free surface equilibrium shapes were computed for several tanks at various contact angles and ullage volumes. A configuration was selected for vibration analysis and lateral slosh mode shapes and natural frequencies were obtained. Results are documented.

Wohlen, R. L.↗

NASA NTRS1987

Significance of foundation-soil separation in dynamic soil-structure interaction

THe dynamic response of flexible surface strip-foundations allowed to uplift is numerically obtained for externally applied forces of a transient time variation. The soil medium is represented by an isotropic, homogeneous and linear half-space. The soil is treated by a time domain boundary element method, while the flexible foundation is treated by the finite element method. It was concluded that intermediate relative stiffness leads to moderate deformations when uplift is permitted. Very flexible footings produce higher deformations in unilateral contact compared to bilateral contact, and thus should be considered in their design. Unilateral contact does not significantly increase deformations for stiff footings subjected to concentrated central loading. However, relatively large deformation differences occur when the loading is eccentric, necessitating consideration of uplift in their design.

Spyrakos, Constantine C.↗

NASA NTRS1990

ZIP3D: An elastic and elastic-plastic finite-element analysis program for cracked bodies

ZIP3D is an elastic and an elastic-plastic finite element program to analyze cracks in three dimensional solids. The program may also be used to analyze uncracked bodies or multi-body problems involving contacting surfaces. For crack problems, the program has several unique features including the calculation of mixed-mode strain energy release rates using the three dimensional virtual crack closure technique, the calculation of the J integral using the equivalent domain integral method, the capability to extend the crack front under monotonic or cyclic loading, and the capability to close or open the crack surfaces during cyclic loading. The theories behind the various aspects of the program are explained briefly. Line-by-line data preparation is presented. Input data and results for an elastic analysis of a surface crack in a plate and for an elastic-plastic analysis of a single-edge-crack-tension specimen are also presented.

Shivakumar, K. N.↗

DOE OSTI2022

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗

DOE OSTI2021

Stratigraphic Identification with Airborne Electromagnetic Methods at the Hanford Site, Washington

Stratigraphic units can influence the fate and transport of subsurface contaminants within groundwater. Units having coarse-grained sediments act as preferential flow pathways, and therefore can accelerate the transport of contaminants to reach human and ecological receptors. At legacy waste sites, detailed knowledge of subsurface stratigraphy can be used for effective monitoring and remediation planning to help minimize risk to human health and the environment. Airborne electromagnetic (AEM) methods can non-invasively provide information on kilometer-scale or larger subsurface stratigraphic features and fill informational gaps in directly sampled data from sparsely located boreholes. In this paper, we present inversion results of a 412 line-km frequency-domain AEM survey to delineate subsurface stratigraphic features at the Hanford Site, located in southeastern Washington State. The inversion was performed using a massively parallel 3D electromagnetic modeling and inversion code, where the modeling is based on solving frequency-domain Maxwell’s equations using an unstructured-mesh finite-element method and the inversion employs a Gauss-Newton optimization scheme. The results are compared to an underlying geologic framework model (GFM), built by interpolating contact depths of stratigraphic units interpreted from site borehole datasets. In areas with good borehole coverage, the inversion results show a good match with the GFM to a depth of about 60 m. Outside of these areas, the inversion results exhibit inconsistencies from the assumptions made to create the GFM, demonstrating that the AEM survey results can be used to improve the understanding of the geological conceptual model.

47 OTHER INSTRUMENTATION↗

NASA NTRS1974

Elastic-plastic analysis of a propagating crack under cyclic loading

Development and application of a two-dimensional finite-element analysis to predict crack-closure and crack-opening stresses during specified histories of cyclic loading. An existing finite-element computer program which accounts for elastic-plastic material behavior under cyclic loading was modified to account for changing boundary conditions - crack growth and intermittent contact of crack surfaces. This program was subsequently used to study the crack-closure behavior under constant-amplitude and simple block-program loading.

Newman, J. C., Jr.↗

DOE OSTI2020

Elevated temperature mechanical properties of Inconel 617 surface oxide using nanoindentation

Inconel 617 is a principal candidate material for helium gas cooled very-high-temperature reactors with outlet temperatures of 700–950 °C. Recent findings showed that this alloy develops unique surface oxide layers especially at high temperature (HT) helium environment with distinctive wear, friction and contact properties. This study investigates the elevated temperature mechanical properties of Inconel 617 top surface layers aged in HT helium environment. Nanoindentation technique is used to obtain load-displacement graphs of the alloy top surface oxide in temperatures ranging from 25 °C up to 600 °C. In addition, using finite element analysis along with an iterative regression technique, a semi-numerical method is developed to further measure and quantify the material parameters and, in particular, time-independent and creep characteristics of the oxide. Although Young's modulus of the oxide is found to be relatively close to the bulk for the tested temperatures, the yield strength and hardness, in comparison to the bulk material, increase significantly as the material is oxidized after aging. The oxide exhibits significant softening as the temperature increases to 600 °C. Unlike the bulk material, diffusion through the grains is found to be the dominant creep mechanism for the oxide. Considerable difference between the mechanical properties of the oxide and the bulk material shows the need for accurate measurements of near surface mechanical properties, if reliable predictive contact and tribological models are sought at elevated temperatures.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Delamination imaging in composites using cross-correlation method by non-contact air-coupled Lamb waves

Abstract This paper presents a cross-correlation damage detection technique using damage scattered waves of all directions cross-correlate with the incident waves for delamination imaging in carbon fiber reinforced polymer (CFRP) composite plate. The presented imaging can not only detect the delamination but also precisely quantify the damage with locations and sizes. First, a CFRP composite sample with interlaminar delamination by inserting the Teflon layer is designed and manufactured in house using the hot-press machine. And a three-dimensional model of composite sample is simulated using the finite element method. Next, the cross-correlation imaging algorithm is introduced in detail, and the cross-correlation imaging proof-of-concept study is conducted with the simulated Lamb waves in the composite sample. Finally, a fully non-contact air-coupled transducer and scanning laser Doppler vibrometer system with a single-mode Lamb wave method is established to actuate and sense the interrogating Lamb waves in the structure. The imaging method is experimentally implemented for the one delamination and two delaminations imaging and quantification.

Instruments & Instrumentation↗

NASA NTRS1974

Application of a substructuring technique to the problem of crack extension and closure

A substructuring technique, originally developed for the efficient reanalysis of structures, is incorporated into the methodology associated with the plastic analysis of structures. An existing finite-element computer program that accounts for elastic-plastic material behavior under cyclic loading was modified to account for changing kinematic constraint conditions - crack growth and intermittent contact of crack surfaces in two dimensional regions. Application of the analysis is presented for a problem of a centercrack panel to demonstrate the efficiency and accuracy of the technique.

Armen, H., Jr.↗

NASA NTRS1982

On elastic-plastic analysis of an overloaded breech ring using NASTRAN

The piece-wise linear analysis option of the NASTRAN code was used to analyze a photoplastic model for a sliding breech mechanism. A two dimensional finite element representation for the breech ring was chosen and the material was made of polycarbonate resin. The aluminum block was regarded as rigid and the width of contact was assumed to remain unchanged during loading. The displacements and stresses in the breech ring were obtained for loading in the elastic as well as plastic range. The maximum tensile stresses before and after complete unloading were obtained and compared with numerical and experimental results.

Chen, P. C. T.↗

NASA NTRS1985

Analysis of impact response of buckled composite laminates

The dynamic response of buckled composite plates impacted by a hard object is studied. A nine-node isoparametric finite element, developed based on the Mindlin plate theory and von Karman large deflection assumptions, is used for both postbuckling and impact analyses. In order to evaluate the contact force, an experimentally established contact law which accounts for the permanent indentation is employed. The static postbuckling problem of laminated plates is first considered. The impact responses include contact force histories, dynamic deflections, and dynamic strains in the plate for various buckling conditions. In addition, free vibration of buckled laminated plates is also solved to determine the range of natural frequencies which are needed to choose a proper time increment for time integration in the impact analysis.

Chen, J. K.↗

NASA NTRS1995

Bicycle ergometer instrumentation to determine muscle and bone forces during exercise

It is hypothesized that bone loss experienced by astronauts in zero gravity conditions may be curtailed by appropriate exercise. According to Wolf's law, bone regenerates when muscles produce stresses by pulling on the bone during daily activity and/or exercise on Earth. to use this theory to prevent or decrease bone loss, one needs to quantify musculoskeletal loads and relate them to bone density changes. In the context of the space program, it is desirable to determine musculoskeletal loads during exercise (using the bicycle ergometer in this case) so that one may make similar measurements on Earth and in space. In this manner, load measurements on Earth may be used as reference to generate similar loads during exercise in space. The work reported in this document entails a musculoskeletal load measurement system that, when complete, will provide forces at muscle insertion points and other contact points, on bone. This data will be used by Dr. Beth A. Todd, who is also a SSF working with Dr. Shackelford, as input to a finite element model of bone sections to determine stress distributions. A bicycle ergometer has been instrumented to measure parameters needed to determine musculoskeletal forces during exercise. A primary feature of the system is its compactness. It uses small/light sensors without line-of-sight requirements. The system developed includes sensors, signal processing, a data acquisition system, and software to collect the data. The sensors used include optical encoders to measure position and orientation of the pedal (foot), accelerometers to determine kinematic parameters of the shank and thigh, load cells to measure pedal forces on the sagittal plane, and EMG probes to measure muscle activity. The signals are processed using anti-aliasing filters and amplifiers. The sensors' output is digitized using 30 channels of a board mounted inside a 486 class PC. A program sets the data acquisition parameters and collects data during a time period specified by the user. The data is put directly into a file on the hard disk in binary form. The 30 channels are sampled at 200 KHz, and each 30 channel scan is done at a rate of 1000 Hz. The instrumented ergometer has been flown in the KC-135 zero-gravity (zero-g) flight to collect information needed to determine musculoskeletal forces under these conditions. Similar information has been collected in 1-g conditions for comparision with the results from the zero-g case. At this time, the sets of data from both experiments are being processed. An existing methodology will be used to determine the kinematic parameters of the shank and thigh using accelerometer and encoder data. This methodology was developed during the fellow's previous NASA/ASEE fellowship and thanks to a Director's Grant. In the future, a methodology to determine the musculoskeletal forces using Newton's Law of Motion and optimization techniques will be developed to determine forces exerted by particular muscles.

Figueroa, Fernando↗

NASA NTRS1999

Fracture Mechanics of Thin, Cracked Plates Under Tension, Bending and Out-of-Plane Shear Loading

Cracks in the skin of aircraft fuselages or other shell structures can be subjected to very complex stress states, resulting in mixed-mode fracture conditions. For example, a crack running along a stringer in a pressurized fuselage will be subject to the usual in-plane tension stresses (Mode-I) along with out-of-plane tearing stresses (Mode-III like). Crack growth and initiation in this case is correlated not only with the tensile or Mode-I stress intensity factor, K(sub I), but depends on a combination of parameters and on the history of crack growth. The stresses at the tip of a crack in a plate or shell are typically described in terms of either the small deflection Kirchhoff plate theory. However, real applications involve large deflections. We show, using the von-Karman theory, that the crack tip stress field derived on the basis of the small deflection theory is still valid for large deflections. We then give examples demonstrating the exact calculation of energy release rates and stress intensity factors for cracked plates loaded to large deflections. The crack tip fields calculated using the plate theories are an approximation to the actual three dimensional fields. Using three dimensional finite element analyses we have explored the relationship between the three dimensional elasticity theory and two dimensional plate theory results. The results show that for out-of-plane shear loading the three dimensional and Kirchhoff theory results coincide at distance greater than h/2 from the crack tip, where h/2 is the plate thickness. Inside this region, the distribution of stresses through the thickness can be very different from the plate theory predictions. We have also explored how the energy release rate varies as a function of crack length to plate thickness using the different theories. This is important in the implementation of fracture prediction methods using finite element analysis. Our experiments show that under certain conditions, during fatigue crack growth, the presence of out-of-plane shear loads induces a great deal of contact and friction on the crack surfaces, dramatically reducing crack growth rate. A series of experiments and a proposed computational approach for accounting for the friction is discussed.

Zehnder, Alan T.↗

NASA NTRS2016

Post-Buckling Analysis of Curved Honeycomb Sandwich Panels Containing Interfacial Disbonds

A numerical study on the effect of facesheet-core disbonds on the post-buckling response of curved honeycomb sandwich panels is presented herein. This work was conducted as part of the development of a damage tolerance plan for the next-generation Space Launch System heavy lift launch vehicle payload fairing. As such, the study utilized full-scale fairing barrel segments as the structure of interest. The panels were composed of carbon fiber reinforced polymer facesheets and aluminum honeycomb core. The panels were analyzed numerically using the finite element method incorporating geometric nonlinearity. In a predetermined circular region, facesheet and core nodes were detached to simulate a disbond, between the outer mold line facesheet and honeycomb core, induced via low-speed impact. Surface-to-surface contact in the disbonded region was invoked to prevent interpenetration of the facesheet and core elements and obtain realistic stresses in the core. The diameter of this disbonded region was varied and the effect of the size of the disbond on the post-buckling response was observed. Significant changes in the slope of the edge load-deflection response were used to determine the onset of global buckling and corresponding buckling load. Finally, several studies were conducted to determine the sensitivity of the numerical predictions to refinement in the finite element mesh.

elastic buckling↗

NASA NTRS1996

Free-Surface and Contact Line Motion of Liquid in Microgravity

This project involves fundamental studies of the role of nonlinearity in determining the motion of liquid masses under the principal influences of surface tension, viscosity and inertia. Issues to be explored are relevant to aspects of terrestrial processes, as well as being immediately applicable to fluid management in a low-gravity environment. Specific issues include: (1) the mechanic's of liquid masses in large-amplitude motions, (2) the influence of bounding surfaces on the motion, and (3) the ability of such surfaces to control liquid motion by wetting forces, especially when they are augmented by various surface treatments. Mathematical techniques include asymptotic analysis of the governing equations, for problem simplification, and numerical simulation, using both boundary-element and finite-difference methods. The flow problem is divided into an 'outer' or inviscid potential-flow region and one or more inner, or viscous dominated, regions. Relevant to one inner region, the vicinity of the contact line, we discuss time-dependent simulation of slow droplet motion, on a surface of variable wettability, using the lubrication approximation. The simulation uses a disjoining pressure model and reproduces realistic wetting-dewetting behavior.

Schwartz, Leonard W.↗

DOE OSTI2020

Sample shape and boundary dependence of measured transverse thermal properties

In spite of increased interest in thermal Hall measurements for the analysis of insulating quantum materials, there remains large uncertainty in such measurements due to contact misalignment. In this paper, we propose that sample geometry and uncertain boundary conditions provide an additional source of uncertainty in the measurement of $D_{xy}$ or $κ_{xy}$. By running simple simulations in an open source finite-element solver, we demonstrate that measured $D_{xy}$ can be changed by a non-negligible fraction in samples with similar width and length. This geometric corrective factor depends on the distinction between a uniform heat flow and constant temperature boundary couplings to a bath. Sample geometry and boundary conditions can be accounted for through simulation or by using rectangular samples to make thermal Hall measurements more reliable and reproducible. Finally, we detail a contactless optical method for measuring $D_{xy}$ based on an existing photothermal microscope technique. This method is insensitive to the longitudinal diffusivity pollution caused by contact misalignment.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Self-assembled multifunctional nanostructures for surface passivation and photon management in silicon photovoltaics

This work reports the fabrication and characterization of multifunctional, nanostructured passivation layers formed using a self-assembly process that provide both surface passivation and improved light trapping in crystalline silicon photovoltaic (PV) cells. Scalable block copolymer self-assembly and vapor phase infiltration processes are used to form arrays of aluminum oxide nanostructures (Al 2 O 3 ) on crystalline silicon without substrate etching. The Al 2 O 3 nanostructures are characterized using scanning electron microscopy (SEM), transmission electron microscopy (TEM), and spectroscopic ellipsometry. Injection-level dependent photoconductance measurements are used to determine the effective carrier lifetime of the samples to confirm the nanostructures successfully passivate the Si surface. Finite element method simulations and reflectance measurement show that the nanostructures increase the internal rear reflectance of the PV cell by suppressing the parasitic optical losses in the metal contact. An optimized morphology of the structures is identified for their potential use in PV cells as multifunctional materials providing surface passivation, photon management, and carrier transport pathways.

14 SOLAR ENERGY↗