Search NASA⌕ Search

SEARCH · Search NASA

Results for “Surface elasticity”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Elastic-Plastic Fracture Mechanics Analysis of Critical Flaw Size in ARES I-X Flange-to-Skin Welds

NASA's Ares 1 Upper Stage Simulator (USS) is being fabricated from welded A516 steel. In order to insure the structural integrity of these welds it is of interest to calculate the critical initial flaw size (CIFS) to establish rational inspection requirements. The CIFS is in turn dependent on the critical final flaw size (CFS), as well as fatigue flaw growth resulting from transportation, handling and service-induced loading. These calculations were made using linear elastic fracture mechanics (LEFM), which are thought to be conservative because they are based on a lower bound, so called elastic, fracture toughness determined from tests that displayed significant plasticity. Nevertheless, there was still concern that the yield magnitude stresses generated in the flange-to-skin weld by the combination of axial stresses due to axial forces, fit-up stresses, and weld residual stresses, could give rise to significant flaw-tip plasticity, which might render the LEFM results to be non-conservative. The objective of the present study was to employ Elastic Plastic Fracture Mechanics (EPFM) to determine CFS values, and then compare these values to CFS values evaluated using LEFM. CFS values were calculated for twelve cases involving surface and embedded flaws, EPFM analyses with and without plastic shakedown of the stresses, LEFM analyses, and various welding residual stress distributions. For the cases examined, the computed CFS values based on elastic analyses were the smallest in all instances where the failures were predicted to be controlled by the fracture toughness. However, in certain cases, the CFS values predicted by the elastic-plastic analyses were smaller than those predicted by the elastic analyses; in these cases the failure criteria were determined by a breakdown in stress intensity factor validity limits for deep flaws (a greater than 0.90t), rather than by the fracture toughness. Plastic relaxation of stresses accompanying shakedown always increases the calculated CFS values compared to the CFS values determined without shakedown. Thus, it is conservative to ignore shakedown effects.

Chell, G. Graham↗
NASA NTRS1985

Fundamental tribological properties of ceramics

When a ceramic is brought into contact with itself, another ceramic, or a metal, strong bond forces can develop between the materials. Adhesion between a ceramic and itself or another solid are discussed from a theoretical consideration of the nature of the surfaces and experimentally by relating bond forces to the interface resulting from solid state contact. Elastic, plastic, and fracture behavior of ceramics in solid-state contact are discussed as they relate to friction and wear. The contact load necessary to initiate fracture in ceramics is shown to be appreciably reduced with tangential motion. Both friction and wear of ceramics are anisotropic and relate to crystal structure as with metals. Both free energy of oxide formation and the d valence bond character of metals are related to the friction and wear characteristics for metals in contact with ceramics. Lubrication is found to increase the critical load necessary to initiate fracture of ceramics with sliding or rubbing contact.

Buckley, D. H.↗
NASA NTRS1985

Fundamental tribological properties of ceramics

When a ceramic is brought into contact with itself, another ceramic, or a metal, strong bond forces can develop between the materials. Adhesion between a ceramic and itself or another solid are discussed from a theoretical consideration of the nature of the surfaces and experimentally by relating bond forces to the interface resulting from solid state contact. Elastic, plastic, and fracture behavior of ceramics in solid-state contact are discussed as they relate to friction and wear. The contact load necessary to initiate fracture in ceramics is shown to be appreciably reduced with tangential motion. Both friction and wear of ceramics are anisotropic and relate to crystal structure as with metals. Both free energy of oxide formation and the d valence bond character of metals are related to the friction and wear characteristics for metals in contact with ceramics. Lubrication is found to increase the critical load necessary to initiate fracture of ceramics with sliding or rubbing contact.

Buckley, D. H.↗
NASA NTRS1987

Coupled NASTRAN/boundary element formulation for acoustic scattering

A coupled finite element/boundary element capability is described for calculating the sound pressure field scattered by an arbitrary submerged 3-D elastic structure. Structural and fluid impedances are calculated with no approximation other than discretization. The surface fluid pressures and normal velocities are first calculated by coupling a NASTRAN finite element model of the structure with a discretized form of the Helmholtz surface integral equation for the exterior field. Far field pressures are then evaluated from the surface solution using the Helmholtz exterior integral equation. The overall approach is illustrated and validated using a known analytic solution for scattering from submerged spherical shells.

Everstine, Gordon C.↗
NASA NTRSDate not supplied

Mesh Deformation Boundary Conditions for Three-Dimensional Ablation Solvers

A method for determining boundary node displacements for three-dimensional ablation solvers is developed. The approach is applicable to a wide range of mesh motion algorithms used to determine internal node displacements. The method leverages radial basis functions to analytically define the geometry to which the boundary nodes should be constrained. Constrained optimization, elasticity, and spring analogy solvers are utilized to redistribute the nodes on the evolved geometry surface while maintaining mesh quality. The method is implemented in the CHarring Ablator Response (CHAR) code, and a realistic demonstration case is presented on the Boundary Layer Transition (BOLT) flight experiment configuration. The method is shown to be robust for highly-curved geometries with large domain deformations.

Adam J Amar↗
NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗
NASA NTRS2017

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗
NASA NTRS1974

The use of COD and plastic instability in crack propagation and arrest in shells

The initiation, growth, and possible arrest of fracture in cylindrical shells containing initial defects are dealt with. For those defects which may be approximated by a part-through semi-elliptic surface crack which is sufficiently shallow so that part of the net ligament in the plane of the crack is still elastic, the existing flat plate solution is modified to take into account the shell curvature effect as well as the effect of the thickness and the small scale plastic deformations. The problem of large defects is then considered under the assumptions that the defect may be approximated by a relatively deep meridional part-through surface crack and the net ligament through the shell wall is fully yielded. The results given are based on an 8th order bending theory of shallow shells using a conventional plastic strip model to account for the plastic deformations around the crack border.

Erdogan, F.↗
NASA NTRS1995

Resistance to Rolling in the Adhesive Contact of Two Elastic Spheres

For the stability of agglomerates of micron sized particles it is of considerable importance to study the effects of tangential forces on the contact of two particles. If the particles can slide or roll easily over each other, fractal structures of these agglomerates will not be stable. We use the description of contact forces by Johnson, Kendall and Roberts, along with arguments based on the atomic structure of the surfaces in contact, in order to calculate the resistance to rolling in such a contact. It is shown that the contact reacts elastically to torque forces up to a critical bending angle. Beyond that, irreversible rolling occurs. In the elastic regime, the moment opposing the attempt to roll is proportional to the bending angle and to the pull-off force P(sub c). Young's modulus of the involved materials has hardly any influence on the results. We show that agglomerates of sub-micron sized particles will in general be quite rigid and even long chains of particles cannot be bent easily. For very small particles, the contact will rather break than allow for rolling. We further discuss dynamic properties such as the possibility of vibrations in this degree of freedom and the typical amount of rolling during a collision of two particles.

Dominik, C.↗
NASA NTRS1995

Resistance to Rolling in the Adhesive Contact of Two Elastic Spheres

For the stability of agglomerates of micron sized particles it is of considerable importance to study the effects of tangential forces on the contact of two particles. If the particles can slide or roll easily over each other, fractal structures of these agglomerates will not be stable. We use the description of contact forces by Johnson, Kendall and Roberts, along with arguments based on the atomic structure of the surfaces in contact, in order to calculate the resistance to rolling in such a contact. It is shown that the contact reacts elastically to torque forces up to a critical bending angle. Beyond that, irreversible rolling occurs. In the elastic regime, the moment opposing the attempt to roll is proportional to the bending angle and to the pull-off force P(sub c). Young's modulus of the involved materials has hardly any influence on the results. We show that agglomerates of sub-micron sized particles will in general be quite rigid and even long chains of particles cannot be bent easily. For very small particles, the contact will rather break than allow for rolling. We further discuss dynamic properties such as the possibility of vibrations in this degree of freedom and the typical amount of rolling during a collision of two particles.

Dominik, C.↗
NASA NTRS2016

Tangle-Free Mesh Motion for Ablation Simulations

Problems involving mesh motion-which should not be mistakenly associated with moving mesh methods, a class of adaptive mesh redistribution techniques-are of critical importance in numerical simulations of the thermal response of melting and ablative materials. Ablation is the process by which material vaporizes or otherwise erodes due to strong heating. Accurate modeling of such materials is of the utmost importance in design of passive thermal protection systems ("heatshields") for spacecraft, the layer of the vehicle that ensures survival of crew and craft during re-entry. In an explicit mesh motion approach, a complete thermal solve is first performed. Afterwards, the thermal response is used to determine surface recession rates. These values are then used to generate boundary conditions for an a posteriori correction designed to update the location of the mesh nodes. Most often, linear elastic or biharmonic equations are used to model this material response, traditionally in a finite element framework so that complex geometries can be simulated. A simple scheme for moving the boundary nodes involves receding along the surface normals. However, for all but the simplest problem geometries, evolution in time following such a scheme will eventually bring the mesh to intersect and "tangle" with itself, inducing failure. This presentation demonstrates a comprehensive and sophisticated scheme that analyzes the local geometry of each node with help from user-provided clues to eliminate the tangle and enable simulations on a wide-class of difficult problem geometries. The method developed is demonstrated for linear elastic equations but is general enough that it may be adapted to other modeling equations. The presentation will explicate the inner workings of the tangle-free mesh motion algorithm for both two and three-dimensional meshes. It will show abstract examples of the method's success, including a verification problem that demonstrates its accuracy and correctness. The focus of the presentation will be on the algorithm; specifics on how the techniques may be used in spacecraft design will be not discussed.

Droba, Justin↗
NASA NTRS1996

Static test induced loads verification beyond elastic limit

Increasing demands for reliable and least-cost high performance aerostructures are pressing design analyses, materials, and manufacturing processes to new and narrowly experienced performance and verification technologies. This study assessed the adequacy of current experimental verification of the traditional binding ultimate safety factor which covers rare events in which no statistical design data exist. Because large, high-performance structures are inherently very flexible, boundary rotations and deflections under externally applied loads approaching fracture may distort their transmission and unknowingly accept submarginal structures or prematurely fracturing reliable ones. A technique was developed, using measured strains from back-to-back surface mounted gauges, to analyze, define, and monitor induced moments and plane forces through progressive material changes from total-elastic to total inelastic zones within the structural element cross section. Deviations from specified test loads are identified by the consecutively changing ratios of moment-to-axial load.

Verderaime, V.↗
NASA NTRS1996

Static test induced loads verification beyond elastic limit

Increasing demands for reliable and least-cost high-performance aerostructures are pressing design analyses, materials, and manufacturing processes to new and narrowly experienced performance and verification technologies. This study assessed the adequacy of current experimental verification of the traditional binding ultimate safety factor which covers rare events in which no statistical design data exist. Because large high-performance structures are inherently very flexible, boundary rotations and deflections under externally applied loads approaching fracture may distort their transmission and unknowingly accept submarginal structures or prematurely fracturing reliable ones. A technique was developed, using measured strains from back-to-back surface mounted gauges, to analyze, define, and monitor induced moments and plane forces through progressive material changes from total-elastic to total-inelastic zones within the structural element cross section. Deviations from specified test loads are identified by the consecutively changing ratios of moment-to-axial load.

Verderaime, V.↗
NASA NTRS1991

Advanced development of the boundary element method for elastic and inelastic thermal stress analysis

The focus of this dissertation is on advanced development of the boundary element method for elastic and inelastic thermal stress analysis. New formulations for the treatment of body forces and nonlinear effects are derived. These formulations, which are based on particular integral theory, eliminate the need for volume integrals or extra surface integrals to account for these effects. The formulations are presented for axisymmetric, two and three dimensional analysis. Also in this dissertation, two dimensional and axisymmetric formulations for elastic and inelastic, inhomogeneous stress analysis are introduced. The derivatives account for inhomogeneities due to spatially dependent material parameters, and thermally induced inhomogeneities. The nonlinear formulation of the present work are based on an incremental initial stress approach. Two inelastic solutions algorithms are implemented: an iterative; and a variable stiffness type approach. The Von Mises yield criterion with variable hardening and the associated flow rules are adopted in these algorithms. All formulations are implemented in a general purpose, multi-region computer code with the capability of local definition of boundary conditions. Quadratic, isoparametric shape functions are used to model the geometry and field variables of the boundary (and domain) of the problem. The multi-region implementation permits a body to be modeled in substructured parts, thus dramatically reducing the cost of analysis. Furthermore, it allows a body consisting of regions of different (homogeneous) material to be studied. To test the program, results obtained for simple test cases are checked against their analytic solutions. Thereafter, a range of problems of practical interest are analyzed. In addition to displacement and traction loads, problems with body forces due to self-weight, centrifugal, and thermal loads are considered.

Henry, Donald P., Jr.↗
NASA NTRS1975

A three-dimensional viscoelastic model of a strike slip fault

An analytic approximation to the Green's function for the displacements due to a strike slip point source in an elastic layer over a viscoelastic half-space is developed. This approximate Green's function is useful because it can be analytically integrated over the fault surface. Comparison with a numerical integration of the exact solution integral indicates that the approximation is quite good. The approximate Green's function is integrated analytically to obtain the displacements due to a finite rectangular strike slip fault in an elastic layer over a viscoelastic half-space. Ground displacements and angle changes from a model survey net are computed to illustrate the viscoelastic relaxation which follows a fracture in the elastic region.

Rundle, J. B.↗
NASA NTRSDate not supplied

Simulation of Oxidative Etch Pit Formation and Growth on FiberForm

Erosion of carbon surfaces due to oxidation does not occur uniformly but through the formation of localized etch pits because of active surface sites. These active sites are formed due to the presence of atomic defects on the carbon surface and have much higher reactivity compared to average non-defective sites. Thus, these active sites are the first to react during ablation, resulting in their removal. This causes all the neighboring atoms to be defective, and increasing their reactivity, thus leading to localized carbon removal around these “active” sites. In this manner, these highly reactive defects serve as nucleation sites for the formation and growth of etch pits, with detrimental effects on the structural integrity. In order to understand the influence of these etch pits on the material properties of carbon fiber microstructures, we have developed a new capability within direct simulation Monte Carlo (DSMC) to capture the etch pit formation process. This capability is developed within the DSMC code SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) and can model the material removal in the presence of active sites leading to the formation of etch pits. The focus of the current work will be to study the effect of etch pits on the material properties of FiberForm, the precursor substrate of the PICA Thermal Protection System (TPS) material. The microstructure of virgin FiberForm obtained directly from X-ray microtomography scans is used within SPARTA to obtain the ablated geometries with etch pits. These pitted microstructures are then imported into the Porous Microstructure Analysis (PuMA) software and various material properties such as thermal conductivity, elasticity, and permeability are computed. The variation of these properties because of the complex evolution of the surface topology due to the formation of etch pits is studied and analyzed. Furthermore, the effect of pitting is compared to the case of uniform radial shrinking of fibers, which has been the standard for modelling ablation of carbon structures, and significant differences are observed. Thus, a physically realistic model of material removal through the formation of etch pits will be helpful in predicting the degradation of carbon-based TPS more accurately during oxidation; as well as other mechanisms such as spallation, which involves the removal of chunks of material into the flow due to the growth of etch pits. This will ultimately improve our understanding of the failure modes in these materials due to ablation.

DSMC↗
NASA NTRSDate not supplied

Simulation of Oxidative Etch Pit Formation and Growth on FiberForm

Erosion of carbon surfaces due to oxidation does not occur uniformly but through the formation of localized etch pits because of active surface sites. These active sites are formed due to the presence of atomic defects on the carbon surface and have much higher reactivity compared to average non-defective sites. Thus, these active sites are the first to react during ablation, resulting in their removal. This causes all the neighboring atoms to be defective, and increasing their reactivity, thus leading to localized carbon removal around these “active” sites. In this manner, these highly reactive defects serve as nucleation sites for the formation and growth of etch pits, with detrimental effects on the structural integrity. In order to understand the influence of these etch pits on the material properties of carbon fiber microstructures, we have developed a new capability within direct simulation Monte Carlo (DSMC) to capture the etch pit formation process. This capability is developed within the DSMC code SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) and can model the material removal in the presence of active sites leading to the formation of etch pits. The focus of the current work will be to study the effect of etch pits on the material properties of FiberForm, the precursor substrate of the PICA Thermal Protection System (TPS) material. The microstructure of virgin FiberForm obtained directly from X-ray microtomography scans is used within SPARTA to obtain the ablated geometries with etch pits. These pitted microstructures are then imported into the Porous Microstructure Analysis (PuMA) software and various material properties such as thermal conductivity, elasticity, and permeability are computed. The variation of these properties because of the complex evolution of the surface topology due to the formation of etch pits is studied and analyzed. Furthermore, the effect of pitting is compared to the case of uniform radial shrinking of fibers, which has been the standard for modelling ablation of carbon structures, and significant differences are observed. Thus, a physically realistic model of material removal through the formation of etch pits will be helpful in predicting the degradation of carbon-based TPS more accurately during oxidation; as well as other mechanisms such as spallation, which involves the removal of chunks of material into the flow due to the growth of etch pits. This will ultimately improve our understanding of the failure modes in these materials due to ablation.

DSMC↗