SEARCH · Engineering Papers
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.
A system which rapidly determines the elastic stiffness constants of materials in plate shape as well as characterize material flaws that are affecting these constants. Rapid (below a minute) nondestructive evaluation system allow for the determination of material stiffness constants, low noise data acquisition algorithm for measuring spectral data. A real time method of displaying leaky Lamb wave spectral data.
In this paper time and/or rate dependent deformation regions are experimentally mapped out as a function of temperature. It is clearly demonstrated that the concept of a threshold stress (a stress that delineate reversible and irreversible behavior) is valid and necessary at elevated temperatures and corresponds to the classical yield stress at lower temperatures. Also the infinitely slow modulus, (Es) i.e. the elastic modulus of the material if it was loaded at an infinitely slow strain rate, and the "dynamic modulus", modulus, Ed, which represents the modulus of the material if it is loaded at an infinitely fast rate are used to delineate rate dependent from rate independent regions. As demonstrated at elevated temperatures there is a significant difference between the two modulus values, thus indicating both significant time-dependence and rate dependence. In the case of the nickel-based super alloy, ME3, this behavior is also shown to be grain size specific. Consequently, at higher temperatures viscoelastic behavior exist below k (i.e., the threshold stress) and at stresses above k the behavior is viscoplastic. Finally a multi-mechanism, stress partitioned viscoelastic model, capable of being consistently coupled to a viscoplastic model is characterized over the full temperature range investigated for Ti-6-4 and ME3.
C-mode scanning acoustic microscopy (C-SAM) is a nondestructive inspection technique that uses ultrasound to show the internal feature of a specimen. A very high or ultra-high-frequency ultrasound passes through a specimen to produce a visible acoustic microimage (AMI) of its inner features. As ultrasound travels into a specimen, the wave is absorbed, scattered or reflected. The response is highly sensitive to the elastic properties of the materials and is especially sensitive to air gaps. This specific characteristic makes AMI the preferred method for finding "air gaps" such as delamination, cracks, voids, and porosity. C-SAM analysis, which is a type of AMI, was widely used in the past for evaluation of plastic microelectronic circuits, especially for detecting delamination of direct die bonding. With the introduction of the flip-chip die attachment in a package; its use has been expanded to nondestructive characterization of the flip-chip solder bumps and underfill. Figure 1.1 compares visual and C-SAM inspection approaches for defect detection, especially for solder joint interconnections and hidden defects. C-SAM is specifically useful for package features like internal cracks and delamination. C-SAM not only allows for the visualization of the interior features, it has the ability to produce images on layer-by-layer basis. Visual inspection; however, is only superior to C-SAM for the exposed features including solder dewetting, microcracks, and contamination. Ideally, a combination of various inspection techniques - visual, optical and SEM microscopy, C-SAM, and X-ray - need to be performed in order to assure quality at part, package, and system levels. This reports presents evaluations performed on various advanced packages/assemblies, especially the flip-chip die version of ball grid array/column grid array (BGA/CGA) using C-SAM equipment. Both external and internal equipment was used for evaluation. The outside facility provided images of the key features that could be detected using the most advanced C-SAM equipment with a skilled operator. Investigation continued using in-house equipment with its limitations. For comparison, representative X-rays of the assemblies were also gathered to show key defect detection features of these non-destructive techniques. Key images gathered and compared are: Compared the images of 2D X-ray and C-SAM for a plastic LGA assembly showing features that could be detected by either NDE technique. For this specific case, X-ray was a clear winner. Evaluated flip-chip CGA and FCBGA assemblies with and without heat sink by C-SAM. Only the FCCGA package that had no heat sink could be fully analyzed for underfill and bump quality. Cross-sectional microscopy did not revealed peripheral delamination features detected by C-SAM. Analyzed a number of fine pitch PBGA assemblies by C-SAM. Even though the internal features of the package assemblies could be detected, C-SAM was unable to detect solder joint failure at either the package or board level. Twenty times touch ups by solder iron with 700degF tip temperature, each with about 5 second duration, did not induce defects to be detected by C-SAM images. Other techniques need to be considered to induce known defects for characterization. Given NASA's emphasis on the use of microelectronic packages and assemblies and quality assurance on workmanship defect detection, understanding key features of various inspection systems that detect defects in the early stages of package and assembly is critical to developing approaches that will minimize future failures. Additional specific, tailored non-destructive inspection approaches could enable low-risk insertion of these advanced electronic packages having hidden and fine features.
The aerospace research community has been working towards the development of a comprehensive framework for the design of composite structures at the component, subassembly, assembly, and system levels. Such a framework requires accurate, efficient and easy-to-use mathematical and computational models that are in sync with physical reality. The research work in this paper is designed to help an aerospace structural engineer move towards this goal by simulating dynamic impact to composite panels. First, the existing orthotropic visco-elastic visco-plastic material model in LS-DYNA commercial finite element program, MAT_COMPOSITE_TABULATED_PLASTICITY_DAMAGE (MAT_213) is enhanced with the development and use of rate-dependent damage data. Using state-of-the-art experimental procedures, impact validation test data are created. Second, the versatility of the constitutive model is examined by its ability to simulate a low velocity (27 ft/s, 8 m/s) impact test, and two high velocity impact tests (385–417 ft/s, 117–127 m/s) that straddle the ballistic limit of the projectile and the composite panels. Third, these impact events are modeled using three different finite element types – under-integrated and fully-integrated solid elements, and fully-integrated thin shell elements. Insight into accuracy versus numerical efficiency is gained by comparing the performance of solid versus thin shell elements. The numerical results show that accurate predictions of these impact events can be obtained with reasonable effort and care.
Apparatus (10) and methods for combining time reversal and elastic nonlinearity of formation materials for qualtitatively probing for over-pressured regions down hole in advance of a well drilling bit, to determine the distance to the over-pressured region, and for accurately measuring pore pressure downhole in a formation, are described. Classical and reciprocal time reversal methods may be utilized to achieve these measurements.
This document is a PowerPoint presentation basically summarizing the publication found here: https://www.sciencedirect.com/science/article/pii/S0022311522002082
Lightning talk presentation for the ignite off competition consisting of 20 picture-focused slides, as submitted to the national competition hosted by ORNL.
The crustal deformation patterns associated with the earthquake cycle can depend strongly on the rheological properties of subcrustal material. Substantial deviations from the simple patterns for a uniformly elastic earth are expected when viscoelastic flow of subcrustal material is considered. The detailed description of the deformation pattern and in particular the surface displacements, displacement rates, strains, and strain rates depend on the structure and geometry of the material near the seismogenic zone. The origin of some of these differences are resolved by analyzing several different linear viscoelastic models with a common finite element computational technique. The models involve strike-slip faulting and include a thin channel asthenosphere model, a model with a varying thickness lithosphere, and a model with a viscoelastic inclusion below the brittle slip plane. The calculations reveal that the surface deformation pattern is most sensitive to the rheology of the material that lies below the slip plane in a volume whose extent is a few times the fault depth. If this material is viscoelastic, the surface deformation pattern resembles that of an elastic layer lying over a viscoelastic half-space. When the thickness or breath of the viscoelastic material is less than a few times the fault depth, then the surface deformation pattern is altered and geodetic measurements are potentially useful for studying the details of subsurface geometry and structure. Distinguishing among the various models is best accomplished by making geodetic measurements not only near the fault but out to distances equal to several times the fault depth. This is where the model differences are greatest; these differences will be most readily detected shortly after an earthquake when viscoelastic effects are most pronounced.
The crustal deformation patterns associated with the earthquake cycle can depend strongly on the rheological properties of subcrustal material. Substantial deviations from the simple patterns for a uniformly elastic earth are expected when viscoelastic flow of subcrustal material is considered. The detailed description of the deformation pattern and in particular the surface displacements, displacement rates, strains, and strain rates depend on the structure and geometry of the material near the seismogenic zone. The origin of some of these differences are resolved by analyzing several different linear viscoelastic models with a common finite element computational technique. The models involve strike-slip faulting and include a thin channel asthenosphere model, a model with a varying thickness lithosphere, and a model with a viscoelastic inclusion below the brittle slip plane. The calculations reveal that the surface deformation pattern is most sensitive to the rheology of the material that lies below the slip plane in a volume whose extent is a few times the fault depth. If this material is viscoelastic, the surface deformation pattern resembles that of an elastic layer lying over a viscoelastic half-space. When the thickness or breath of the viscoelastic material is less than a few times the fault depth, then the surface deformation pattern is altered and geodetic measurements are potentially useful for studying the details of subsurface geometry and structure. Distinguishing among the various models is best accomplished by making geodetic measurements not only near the fault but out to distances equal to several times the fault depth. This is where the model differences are greatest; these differences will be most readily detected shortly after an earthquake when viscoelastic effects are most pronounced.
and vibration suppression device is composed of statically compressed chains of spherical particles. The device superimposes a combination of dissipative damping and dispersive effects. The dissipative damping resulting from the elastic wave attenuation properties of the bulk material selected for the granular media is independent of particle geometry and periodicity, and can be accordingly designed based on the dissipative (or viscoelastic) properties of the material. For instance, a viscoelastic polymer might be selected where broadband damping is desired. In contrast, the dispersive effects result from the periodic arrangement and geometry of particles composing a linear granular chain. A uniform (monatomic) chain of statically compressed spherical particles will have a low-pass filter effect, with a cutoff frequency tunable as a function of particle mass, elastic modulus, Poisson fs ratio, radius, and static compression. Elastic waves with frequency content above this cutoff frequency will exhibit an exponential decay in amplitude as a function of propagation distance. System design targeting a specific application is conducted using a combination of theoretical, computational, and experimental techniques to appropriately select the particle radii, material (and thus elastic modulus and Poisson fs ratio), and static compression to satisfy estimated requirements derived for shock and/or vibration protection needs under particular operational conditions. The selection of a chain of polymer spheres with an elastic modulus .3 provided the appropriate dispersive filtering effect for that exercise; however, different operational scenarios may require the use of other polymers, metals, ceramics, or a combination thereof, configured as an array of spherical particles. The device is a linear array of spherical particles compressed in a container with a mechanism for attachment to the shock and/or vibration source, and a mechanism for attachment to the article requiring isolation (Figure 1). This configuration is referred to as a single-axis vibration suppressor. This invention also includes further designs for the integration of the single-axis vibration suppressor into a six-degree-of-freedom hexapod "Stewart"mounting configuration (Figure 2). By integrating each singleaxis vibration suppressor into a hexapod formation, a payload will be protected in all six degrees of freedom from shock and/or vibration. Additionally, to further enable the application of this device to multiple operational scenarios, particularly in the case of high loads, the vibration suppressor devices can be used in parallel in any array configuration.
A collision of two rocky objects circling the sun in space, each roughly the size and mass of a large mountain range, was modeled. A fragmentation hydrocode was developed to perform dynamical computations of collisional outcomes. Explosive framentation and fluid dynamics were used and drawn together into a single application. To model a solid, certain material parameters, such as density, elasticity, rigidity, and energies of melting and vaporization were input. These parameters are well-known for a variety of important materials, such as ice, iron, granite, and basalt. Another important parameter used is the distribution of initial flaws within the material.
We have demonstrated electrostatic switching in vertically oriented nanotubes or nanofibers, where a nanoprobe was used as the actuating electrode inside an SEM. When the nanoprobe was manipulated to be in close proximity to a single tube, switching voltages between 10 V - 40 V were observed, depending on the geometrical parameters. The turn-on transitions appeared to be much sharper than the turn-off transitions which were limited by the tube-to-probe contact resistances. In many cases, stiction forces at these dimensions were dominant, since the tube appeared stuck to the probe even after the voltage returned to 0 V, suggesting that such structures are promising for nonvolatile memory applications. The stiction effects, to some extent, can be adjusted by engineering the switch geometry appropriately. Nanoscale mechanical measurements were also conducted on the tubes using a custom-built anoindentor inside an SEM, from which preliminary material parameters, such as the elastic modulus, were extracted. The mechanical measurements also revealed that the tubes appear to be well adhered to the substrate. The material parameters gathered from the mechanical measurements were then used in developing an electrostatic model of the switch using a commercially available finite-element simulator. The calculated pull-in voltages appeared to be in agreement to the experimentally obtained switching voltages to first order.
Background The integration of dissimilar materials poses a significant challenge in engineering, necessitating innovative solutions for robust and reliable joining. Interlocking metasurfaces (ILMs) are a new joining technology comprising arrays of autogenous features patterned across two surfaces that interlock to form robust structural joints. Objective Here, this study elucidates the factors influencing the tensile performance of ILM joints formed between dissimilar materials. Methods We employed parametric optimization to identify optimal unit cell geometries for maximal yield strength based on the hypothesis that the elastic tensile properties of the materials are the primary determinants of tensile performance. Experimental validation was performed by mechanically testing the theorized optimal ILM geometry and a range of ILM geometries to capture the overall behavior trends of joints between two additively manufactured polymers, VeroPureWhite (VW) and RGDA8430-DM (8430). Results Experimental validation of optimized designs revealed that additional factors, e.g. flexural strength and localized plasticity, also strongly influenced the tensile performance of T-slot ILMs joining dissimilar materials. The proposed optimal design remained the best performer. Conclusions This study demonstrates the viability of ILMs as a joining method for dissimilar materials. ILMs can join dissimilar materials with no loss in joint yield strength compared to joints composed solely of the weaker of the two constitutive materials. ILMs demonstrated their potential as a versatile and effective joining technology in diverse engineering applications.
Mechanoluminescent materials, which emit light in response to elastic deformation, are demanded for use as in situ stress sensors. ZnS doped with Mn is known to exhibit one of the lowest reported thresholds for appearance of mechanoluminescence, with repeatable light emission under contact pressure <10 MPa. Here, the physical basis for such behavior remains as yet unclear. Here, reliable microscopic detection of mechanoluminescence of single ZnS:Mn microparticles, in combination with nanoscale structural characterization, provides evidence that the mechanoluminescent properties of these particles result from interplay between a non-centrosymmetric crystal lattice and its defects, viz., dislocations and stacking faults. Statistical analysis of the distributions of mechanoluminescence energy release trajectories reveals two distinct mechanisms of excitation: one attributable to a piezo-phototronic effect and the other due to dislocation motion. At pressures below 8.1 MPa, both mechanisms contribute to mechanoluminescent output, with a dominant contribution from the piezo-phototronic mechanism. In contrast, above 8.1 MPa, dislocation motion is the primary excitation source. For the piezo-phototronic mechanism, we propose a specific model that accounts for elastic ZnS:Mn mechanoluminescence under very low pressure. The charged interfaces in stacking faults lead to the presence of filled traps, which otherwise would be empty in the absence of the built-in electric field. Upon application of external stress, local enhancement of the piezoelectric field at the stacking faults’ interfaces facilitates release of the trapped carriers and subsequent luminescence. This field enhancement explains how <10 MPa pressure produces thousands of photons.
Metallic glasses are a new class of structural materials which exhibit exciting mechanical properties including high strength and elasticity. In terms of fracture toughness, the material class of metallic glasses spans a wide range; Some metallic glasses are extremely brittle and exhibit near ideal brittle behavior whereas others can be exceptional tough with values comparable to the toughest metals out there. Such large range of observed fracture toughness within the material class of metallic glasses is surprising as they have seemingly a very similar atomic structure. Therefore, we developed “artificial microstructures” which allows to decouple the various contribution of sample geometry, imperfection and structure. Specifically, we decouple variations in the alloys’ chemical composition and the atomic structure and quantified the resulting fracture toughness. Atomic structure of a metallic glass can be modified by the fictive temperature. The fictive temperature of a glass is the temperature at which the liquid metallic glass falls out of equilibrium upon colling and forms a glass. Upon further cooling the structure is maintained only thermal oscillations decrease due to a lower absolute temperature. We found that the effect of fictive temperature (same chemistry, different structural stages of the glass) is comparable to the variations of fracture toughness when the chemistry is varied. Hence, it appears that the subtle differences in the glass structure are responsible for the large range of fracture toughness’s observed. Our results reveal that fracture toughness within the material class of metallic glasses varies significantly and we found some example alloys with exceptional high resistance to fracture and others that are almost ideally brittle. Significant influences on the fracture toughness have the structure of the glass, its chemistry, and some imperfections in the structure.
We consider the modeling and active damping of an elastic beam using distributed actuators and sensors. The piezoelectric ceramic material (PZT) is used to build the actuator. The sensor is made of the piezoelectric polymer polyvinylidene fluoride (PVDF). These materials are glued on both sides of the beam. For the simple clamped beam, the closed loop controller has been shown to be able to extract energy from the beam. The shape of the actuator and its influence on the closed loop system performance are discussed. It is shown that it is possible to suppress the selected mode by choosing the appropriate actuator layout. It is also shown that by properly installing the sensor and determining the sensor shape we can further extract and manipulate the sensor signal for our control need.
In isotropic materials, the direction of the energy flux (energy per unit time per unit area) of an ultrasonic plane wave is always along the same direction as the normal to the wave front. In anisotropic materials, however, this is true only along symmetry directions. Along other directions, the energy flux of the wave deviates from the intended direction of propagation. This phenomenon is known as energy flux deviation and is illustrated. The direction of the energy flux is dependent on the elastic coefficients of the material. This effect has been demonstrated in many anisotropic crystalline materials. In transparent quartz crystals, Schlieren photographs have been obtained which allow visualization of the ultrasonic waves and the energy flux deviation. The energy flux deviation in graphite/epoxy (gr/ep) composite materials can be quite large because of their high anisotropy. The flux deviation angle has been calculated for unidirectional gr/ep composites as a function of both fiber orientation and fiber volume content. Experimental measurements have also been made in unidirectional composites. It has been further demonstrated that changes in composite materials which alter the elastic properties such as moisture absorption by the matrix or fiber degradation, can be detected nondestructively by measurements of the energy flux shift. In this research, the effects of nonlinear elasticity on energy flux deviation in unidirectional gr/ep composites were studied. Because of elastic nonlinearity, the angle of the energy flux deviation was shown to be a function of applied stress. This shift in flux deviation was modeled using acoustoelastic theory and the previously measured second and third order elastic stiffness coefficients for T300/5208 gr/ep. Two conditions of applied uniaxial stress were considered. In the first case, the direction of applied uniaxial stress was along the fiber axis (x3) while in the second case it was perpendicular to the fiber axis along the laminate stacking direction (x1).
A non-linear, continuum-based constitutive model is developed for carbon nanotube materials in which bundles of aligned carbon nanotubes have varying amounts of crosslinks between the nanotubes. The model accounts for the non-linear elastic constitutive behavior of the material in terms of strain, and is developed using a thermodynamic energy approach. The model is used to examine the effect of the crosslinking on the overall mechanical properties of variations of the crosslinked carbon nanotube material with varying degrees of crosslinking. It is shown that the presence of the crosslinks has significant effects on the mechanical properties of the carbon nanotube materials. An increase in the transverse shear properties is observed when the nanotubes are crosslinked. However, this increase is accompanied by a decrease in axial mechanical properties of the nanotube material upon crosslinking.