Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic fields”

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 505 records · Page 28

Mercury Crustal and Lithospheric Properties from High-resolution Gravity Field Models

We have analyzed all the available tracking data from the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) mission to determine high-resolution gravity field models. We have applied a technique based on line-of-sight accelerations that results in models with improved correlations with topography in areas where the spacecraft altitude was low. We have used these models in an analysis of admittance to determine properties of the crust and lithosphere. We focus on four areas where correlations between gravity and topography are high. These locations represent different areas on Mercury: the high-Mg region, the Strindberg crater plus some lobate scarps, heavily cratered terrain, and smooth plains (the northern rise). We find crustal densities around 2600 kg m-3 for the first three areas, and a higher density for the northern rise, indicative of volcanic materials. Elastic thickness is generally low and varies between 11 and 30 km.

Sander Johannes Goossens↗

A Review of the Proposed K (sub Isi) Offset-Secant Method for Size-Independent Linear-Elastic Toughness Evaluation

The proposed size-independent linear-elastic fracture toughness, K (sub Isi), for potential inclusion in ASTM E399 targets a consistent 0.5 millimeters crack extension for all specimen sizes through an offset secant that is a function of the specimen ligament length. The K (sub Isi) method also includes an increase in allowable deformation, and the removal of the P (sub max)/P (sub Q) criterion. A finite element study of the K (sub Isi) test method confirms the viability of the increased deformation limit, but has also revealed a few areas of concern. Findings: 1. The deformation limit, b (sub o) greater than or equal to 1.1 times (K (sub I) divided by delta (sub ys) squared) maintains a K-dominant crack tip field with limited plastic contribution to the fracture energy; 2. The three dimensional effects on compliance and the shape of the force versus CMOD (Crack-Mouth Opening Displacement) trace are significant compared to a plane strain assumption; 3. The non-linearity in the force versus CMOD trace at deformations higher than the current limit of 2.5 times (K (sub I) divided by delta (sub ys) squared) is sufficient to introduce error or even "false calls" regarding crack extension when using a constant offset secant line. This issue is more significant for specimens with W (width) greater than or equal to 2 inches; 4. A non-linear plasticity correction factor in the offset secant may improve the viability of the method at deformations between 2.5 times (K (sub I) divided by delta (sub ys) squared) and 1.1 times (K (sub I) divided by delta (sub ys) squared).

James, Mark↗

Nonlinear Analysis of Bonded Composite Single-LAP Joints

This study presents a semi-analytical solution method to analyze the geometrically nonlinear response of bonded composite single-lap joints with tapered adherend edges under uniaxial tension. The solution method provides the transverse shear and normal stresses in the adhesive and in-plane stress resultants and bending moments in the adherends. The method utilizes the principle of virtual work in conjunction with von Karman s nonlinear plate theory to model the adherends and the shear lag model to represent the kinematics of the thin adhesive layer between the adherends. Furthermore, the method accounts for the bilinear elastic material behavior of the adhesive while maintaining a linear stress-strain relationship in the adherends. In order to account for the stiffness changes due to thickness variation of the adherends along the tapered edges, their in-plane and bending stiffness matrices are varied as a function of thickness along the tapered region. The combination of these complexities results in a system of nonlinear governing equilibrium equations. This approach represents a computationally efficient alternative to finite element method. Comparisons are made with corresponding results obtained from finite-element analysis. The results confirm the validity of the solution method. The numerical results present the effects of taper angle, adherend overlap length, and the bilinear adhesive material on the stress fields in the adherends, as well as the adhesive, of a single-lap joint

Oterkus, E.↗

Validation of Methods to Predict Vibration of a Panel in the Near Field of a Hot Supersonic Rocket Plume

This paper describes the measurement and analysis of surface fluctuating pressure level (FPL) data and vibration data from a plume impingement aero-acoustic and vibration (PIAAV) test to validate NASA s physics-based modeling methods for prediction of panel vibration in the near field of a hot supersonic rocket plume. For this test - reported more fully in a companion paper by Osterholt & Knox at 26th Aerospace Testing Seminar, 2011 - the flexible panel was located 2.4 nozzle diameters from the plume centerline and 4.3 nozzle diameters downstream from the nozzle exit. The FPL loading is analyzed in terms of its auto spectrum, its cross spectrum, its spatial correlation parameters and its statistical properties. The panel vibration data is used to estimate the in-situ damping under plume FPL loading conditions and to validate both finite element analysis (FEA) and statistical energy analysis (SEA) methods for prediction of panel response. An assessment is also made of the effects of non-linearity in the panel elasticity.

Bremner, P. G.↗

Experimental flaw detection by scattering of plate waves

The surface displacements of elastic waves scattered from a surface-breaking crack are measured with a capacitive transducer on a 2D plate, and comparisons with theoretical predictions are reported. The noncontacting detector has a wide and flat frequency response which measures the response of a glass plate to the propagation of Rayleigh-Lamb waves. A theoretical simulation of the surface displacements is generated by multiplying experimental traces with a cosine window, and Green's function is used to compute the source function. Power spectra are also recorded from the experiment to compare the results with a finite-element model of the crack, and good agreement is reported for results in both the time and frequency domains. The crack can be detected in the far field when the spacing of the receiver locations is small, and the transducer/probe appears to be a practical flaw-detection instrument.

Spetzler, Hartmut↗

Interlaminar shear test method development for long term durability testing of composites

The high speed civil transport is a commercial aircraft that is expected to carry 300 passengers at Mach 2.4 over a range of more than 6000 nautical miles. With the existing commercial structural material technology (i.e., aluminum) the performance characteristics of the high speed civil transport would not be realized. Therefore there has been a concerted effort in the development of light weight materials capable of withstanding elevated temperatures for long duration. Thermoplastic composite materials are such candidate materials and the understanding of how these materials perform over the long term under harsh environments is essential to safe and effective design. The matrix dominated properties of thermoplastic composites are most affected by both time and temperature. There is currently an effort to perform short term testing to predict long term behavior of in-plane mechanical properties E22 (transverse modulus of elasticity) and G12 (shear modulus). Out-of-plane properties such as E33, G13, and G23 are inherently more difficult to characterize. This is especially true for the out-of-plane shear modulus G23 and hence there is no existing acceptable standard test method. Since G23 is the most matrix dominated property, it is essential that a test method be developed. A shear test methodology is developed to do just that. The test method, called the double notched specimen, along with the previously developed shear gage was tested at room temperature. Mechanical testing confirmed the attributes of the methodology. A finite element parametric study was conducted for specimen optimization. Moire interferometry, a high sensitivity laser optical method, was used for full-field analysis of the specimen. From this work, material parameters will be determined and thus enable the prediction of long term material behavior of laminates subjected to general loading states.

Ifju, Peter G.↗

Out of Plane Stresses Using a Beam-Based Elastic Foundation Bonded Joint Model

Efficient, semi-analytical models of bonded joints are utilized for initial joint design and sizing due to their convenience and speed. Many of these models are based on beam and plate theories, which are most applicable when in-plane stresses/strains dominate. While this can be a sensible assumption for obtaining the deformations of the joint, out of plane (transverse) stresses are often of critical importance for failure in bonded joints composed of fiber reinforced composite materials. Since through-thickness strains are neglected in these semi-analytical theories, the present study shows a method to back-calculate the out of plane stresses based on stress balance equations, which use the in-plane stresses and displacements that have already been obtained. The method was implemented within the Joint Element Designer (JED) software, and verification and validation examples have been shown. As is typical for these types of models, corner stress singularities are not captured by the JED model, but the general trends of the stress and deformation fields agreed with a higher fidelity finite element model. As such, this method for obtaining out of plane stresses, when coupled with the JED code, is very useful for the complete stress analysis of bonded composite joints.

composites↗

Elastic Distribution of Microshutters, Measurements Obtainable on James Web Space Telescope

Spectrographic astronomy measurements in the near-infrared region will be done by functional two-dimensional microshutter arrays that are being fabricated at the NASA Goddard Space Flight Center for the James Webb Space Telescope (JWST). These micro-shutter arrays will represent the first mission-critical MEMS devices to be flown in space. JWST will use microshutter arrays to select focal plane object. 2-D programmable aperture masks of more than 200,000 elements select such space object. The use of silicon wafer material promises high efficiency and high contrast. Microshutter operation temperature is around 35K. Microshutter arrays are fabricated as close-packed silicon nitride membranes with a unit cell size of 105 x 204 micrometers. A layer of magnetic material is deposited onto each shutter. Individual shutters are equipped with a torsion flexure. Reactive ion etching (RIE) releases the shutters so they can open up to 90 degrees using the torsion flexure. Shutter rotation is initiated into a silicon support structure via an external magnetic field. Two electrically independent aluminum electrodes are deposited, one onto each shutter and another onto the support structure side-wall, permitting electrostatic latching and 2-D addressing to hold specific shutters open via external electronics.

Kletetschka, Gunther↗

Estimating lithospheric properties at Atla Regio, Venus

Magellan spehrical harmonic gravity and topography models are used to estimate lithospheric properties at Alta Regio, Venus, a proposed hotspot with dynamic support from mantle plume(s). Global spherical harmonic and local representations of the gravity field share common properties in the Atla region interms of their spectral behavior over a wavelength band from approximately 2100 to approximately 700 km. The estimated free-air admittance spectrum displays a rather featureless long-wavelength portion followed by a sharp rise at wavelengths shorter than about 1000 km. This sharp rise requires significant flexural support of short-wavelength structures. The Bouguer coherence also displays a sharp drop in this wavelength band, indicating a finite flexural rigidity of the lithosphere. A simple model for lithospheric loading from above and below is introduced (D. W. Forsyth, 1985) with four parameters: f, the ratio of bottom loading to top loading; z(sub m), crustal thickness; z(sub l) depth to bottom loading source; and T(sub e) elastic lithosphere thickness. A dual-mode compensation model is introduced in which the shorter wavelengths (lambda approximately less than 1000 km) might be explained best by a predominance of top loading by the large shield volcanoes Maat Mons, Ozza Mons, and Sapas Mons, and the longer wavelengths (lambda approximately greater than 1500 km) might be explained best by a deep depth of compensation, possibly representing bottom loading by a dynamic source. A Monte Carlo inversion technique is introduced to thoroughly search out the four-space of the model parameters and to examine parameter correlation in the solutions. Venus either is a considerabe deficient in heat sources relative to Earth, or the thermal lithosphere is overthickened in response to an earlier episode of significant heat loss from the planet.

Phillips, Roger J.↗

Finite Element Methods for real-time Haptic Feedback of Soft-Tissue Models in Virtual Reality Simulators

We have applied the linear elastic finite element method to compute haptic force feedback and domain deformations of soft tissue models for use in virtual reality simulators. Our results show that, for virtual object models of high-resolution 3D data (>10,000 nodes), haptic real time computations (>500 Hz) are not currently possible using traditional methods. Current research efforts are focused in the following areas: 1) efficient implementation of fully adaptive multi-resolution methods and 2) multi-resolution methods with specialized basis functions to capture the singularity at the haptic interface (point loading). To achieve real time computations, we propose parallel processing of a Jacobi preconditioned conjugate gradient method applied to a reduced system of equations resulting from surface domain decomposition. This can effectively be achieved using reconfigurable computing systems such as field programmable gate arrays (FPGA), thereby providing a flexible solution that allows for new FPGA implementations as improved algorithms become available. The resulting soft tissue simulation system would meet NASA Virtual Glovebox requirements and, at the same time, provide a generalized simulation engine for any immersive environment application, such as biomedical/surgical procedures or interactive scientific applications.

Frank, Andreas O.↗

Models of Active Glacial Isostasy Roofing Warm Subduction: Case of the South Patagonian Ice Field

Modern geodetic techniques such as precise Global Positioning System (GPS) and high-resolution space gravity mapping (Gravity Recovery and Climate Experiment, GRACE) make it possible to measure the present-day rate of viscoelastic gravitational Earth response to present and past glacier mass changes. The Andes of Patagonia contain glacial environments of dramatic mass change. These mass load changes occur near a tectonically active boundary between the Antarctic and South American plates. The mechanical strength of the continental side of this boundary is influenced by Neogene ridge subduction and by the subduction of a youthful oceanic slab. A ridge of young volcanos parallels the Pacific coastline. Release of volatiles (such as water) at depth along this ridge creates a unique rheological environment. To assess the influence of this rheological ridge structure on the observational land uplift rate, we apply a two dimensional viscoelastic Earth model. A numerical study is presented which examines the sensitivity of the glacial loading-unloading response to the complex structure at depth related to the subducting slab, the viscous wedge between slab and continental lithosphere, and the increase of elastic thickness from oceanic to continental lithosphere. A key feature revealed by our numerical experiments is a continuum flow wherein the slab subdues the material transport toward oceanic mantle and crust. The restricted flow is sensitive to the details of slab mechanical strength and penetration into the upper mantle. The reduced viscosity within the mantle wedge, however, enhances the load-induced material transport everywhere within the asthenosphere.

crustal mass changes↗

Determining water cloud particle sizes from lidar depolarization measurements and time dependent multiple scattering coefficients

We present a technique of extracting water cloud particle size information from lidar measurements in conjunction with double scattering calculations. In this presentation, we describe the technique and give examples using data taken with the Air Force Phillips Laboratory's (Geophysics Directorate) low altitude Nd:YAG, elastic backscatter lidar. In a related presentation we describe the double scattering lidar model which we developed for this work. The technique uses simultaneous measurements of two orthogonal linear polarization components of lidar returns from water clouds or other media composed of spherical particles. Any depolarization of the incident lidar radiation backscattered by such a media can only be due to multiple scattering. The amount of depolarization is dependent on the extinction coefficient and the single scatter phase matrix, both of which are functions of position in the medium. The phase matrix is dependent on the index of refraction of the particles and the particle size distribution. Our technique is a modification of a procedure presented in Reference 1. There, particle sizes of water clouds are determined from double scattering calculations together with measurements of radiation scattered from volumes outside the lidar receiver's field of view (which can only be multiply scattered radiation). The methodology of our technique is similar but our 'probe' of the scattering phase function (and thus the particle size distribution) is different.

Garner, R. C.↗

Interlaminar Stresses by Refined Beam Theories and the Sinc Method Based on Interpolation of Highest Derivative

Computation of interlaminar stresses from the higher-order shear and normal deformable beam theory and the refined zigzag theory was performed using the Sinc method based on Interpolation of Highest Derivative. The Sinc method based on Interpolation of Highest Derivative was proposed as an efficient method for determining through-the-thickness variations of interlaminar stresses from one- and two-dimensional analysis by integration of the equilibrium equations of three-dimensional elasticity. However, the use of traditional equivalent single layer theories often results in inaccuracies near the boundaries and when the lamina have extremely large differences in material properties. Interlaminar stresses in symmetric cross-ply laminated beams were obtained by solving the higher-order shear and normal deformable beam theory and the refined zigzag theory with the Sinc method based on Interpolation of Highest Derivative. Interlaminar stresses and bending stresses from the present approach were compared with a detailed finite element solution obtained by ABAQUS/Standard. The results illustrate the ease with which the Sinc method based on Interpolation of Highest Derivative can be used to obtain the through-the-thickness distributions of interlaminar stresses from the beam theories. Moreover, the results indicate that the refined zigzag theory is a substantial improvement over the Timoshenko beam theory due to the piecewise continuous displacement field which more accurately represents interlaminar discontinuities in the strain field. The higher-order shear and normal deformable beam theory more accurately captures the interlaminar stresses at the ends of the beam because it allows transverse normal strain. However, the continuous nature of the displacement field requires a large number of monomial terms before the interlaminar stresses are computed as accurately as the refined zigzag theory.

Slemp, Wesley C. H.↗

Drops in Space: Super Oscillations and Surfactant Studies

An unprecedented microgravity observation of maximal shape oscillations of a surfactant-bearing water drop the size of a ping pong ball was observed during a mission of Space Shuttle Columbia as part of the second United States Microgravity Laboratory-USML-2 (STS-73, October 20-November 5, 1995). The observation was precipitated by the action of an intense sound field which produced a deforming force on the drop. When this deforming force was suddenly reduced, the drop executed nearly free and axisymmetric oscillations for several cycles, demonstrating a remarkable amplitude of nonlinear motion. Whether arising from the discussion of modes of oscillation of the atomic nucleus, or the explosion of stars, or how rain forms, the complex processes influencing the motion, fission, and coalescence of drops have fascinated scientists for centuries. Therefore, the axisymmetric oscillations of a maximally deformed liquid drop are noteworthy, not only for their scientific value but also for their aesthetic character. Scientists from Yale University, the Jet Propulsion Laboratory (JPL) and Vanderbilt University conducted liquid drop experiments in microgravity using the acoustic positioning/manipulation environment of the Drop Physics Module (DPM). The Yale/JPL group's objectives were to study the rheological properties of liquid drop surfaces on which are adsorbed surfactant molecules, and to infer surface properties such as surface tension, Gibb's elasticity, and surface dilatational viscosity by using a theory which relies on spherical symmetry to solve the momentum and mass transport equations.

Apfel, Robert E.↗

Effect of Graphene Addition on Shape Memory Behavior of Epoxy Resins

Shape memory polymers (SMPs) and composites are a special class of smart materials known for their ability to change size and shape upon exposure to an external stimulus (e.g. light, heat, pH, or magnetic field). These materials are commonly used for biomedical applications; however, recent attempts have been made towards developing SMPs and composites for use in aircraft and space applications. Implementing SMPs and composites to create a shape change effect in some aircraft structures could potentially reduce drag, decrease fuel consumption, and improve engine performance. This paper discusses the development of suitable materials to use in morphing aircraft structures. Thermally responsive epoxy SMPs and nanocomposites were developed and the shape memory behavior and thermo-mechanical properties were studied. Overall, preliminary results from dynamic mechanical analysis (DMA) showed that thermally actuated shape memory epoxies and nanocomposites possessed Tgs near approximately 168 C. When graphene nanofiller was added, the storage modulus and crosslinking density decreased. On the other hand, the addition of graphene enhanced the recovery behavior of the shape memory nanocomposites. It was assumed that the addition of graphene improved shape memory recovery by reducing the crosslinking density and increasing the elasticity of the nanocomposites.

Williams, Tiffany↗

Shadowgraph Study of Gradient Driven Fluctuations

A fluid or fluid mixture, subjected to a vertical temperature and/or concentration gradient in a gravitational field, exhibits greatly enhanced light scattering at small angles. This effect is caused by coupling between the vertical velocity fluctuations due to thermal energy and the vertically varying refractive index. Physically, small upward or downward moving regions will be displaced into fluid having a refractive index different from that of the moving region, thus giving rise to the enhanced scattering. The scattered intensity is predicted to vary with scattering wave vector q, as q(sup -4), for sufficiently large q, but the divergence is quenched by gravity at small q. In the absence of gravity, the long wavelength fluctuations responsible for the enhanced scattering are predicted to grow until limited by the sample dimensions. It is thus of interest to measure the mean-squared amplitude of such fluctuations in the microgravity environment for comparison with existing theory and ground based measurements. The relevant wave vectors are extremely small, making traditional low-angle light scattering difficult or impossible because of stray elastically scattered light generated by optical surfaces. An alternative technique is offered by the shadowgraph method, which is normally used to visualize fluid flows, but which can also serve as a quantitative tool to measure fluctuations. A somewhat novel shadowgraph apparatus and the necessary data analysis methods will be described. The apparatus uses a spatially coherent, but temporally incoherent, light source consisting of a super-luminescent diode coupled to a single-mode optical fiber in order to achieve extremely high spatial resolution, while avoiding effects caused by interference of light reflected from the various optical surfaces that are present when using laser sources. Results obtained for a critical mixture of aniline and cyclohexane subjected to a vertical temperature gradient will be presented. The sample was confined between two horizontal parallel sapphire plates with a vertical spacing of 1 mm. The temperatures of the sapphire plates were controlled by independent circulating water loops that used Peltier devices to add or remove heat from the room air as required. For a mixture with a temperature gradient, two effects are involved in generating the vertical refractive index gradient, namely thermal expansion and the Soret effect, which generates a concentration gradient in response to the applied temperature gradient. For the aniline/cyclohexane system, the denser component (aniline) migrates toward the colder surface. Consequently, when heating from above, both effects result in the sample density decreasing with altitude and are stabilizing in the sense that no convective motion occurs regardless of the magnitude of the applied temperature gradient. The Soret effect is strong near a binary liquid critical point, and thus the dominant effect is due to the induced concentration gradient. The results clearly show the divergence at low q and the predicted gravitational quenching. Results obtained for different applied temperature gradients at varying temperature differences from the critical temperature, clearly demonstrate the predicted divergence of the thermal diffusion ratio. Thus, the more closely the critical point is approached, the smaller becomes the temperature gradient required to generate the same signal. Two different methods have been used to generate pure concentration gradients. In the first, a sample cell was filled with a single fluid, ethylene glycol, and a denser miscible fluid, water, was added from below thus establishing a sharp interface to begin the experiment. As time went on the two fluids diffused into each other, and large amplitude fluctuations were clearly observed at low q. The effects of gravitational quenching were also evident. In the second method, the aniline/cyclohexane sample was used, and after applying a vertical temperature gradient for several hours, the top and bottom temperatures were set equal and the thermal gradient died on a time scale of seconds, leaving the Soret induced concentration gradient in place. Again, large-scale fluctuations were observed and died away slowly in amplitude as diffusion destroyed the initial concentration gradient.

Cannell, David↗

Separation of crack extension modes in orthotropic delamination models

In the analysis of an interface crack between dissimilar elastic materials, the mode of crack extension is typically not unique, due to oscillatory behavior of near-tip stresses and displacements. This behavior currently limits the applicability of interfacial fracture mechanics as a means to predict composite delamination. The Virtual Crack Closure Technique (VCCT) is a method used to extract mode 1 and mode 2 energy release rates from numerical fracture solutions. The mode of crack extension extracted from an oscillatory solution using the VCCT is not unique due to the dependence of mode on the virtual crack extension length, Delta. In this work, a method is presented for using the VCCT to extract Delta-independent crack extension modes for the case of an interface crack between two in-plane orthotropic materials. The method does not involve altering the analysis to eliminate its oscillatory behavior. Instead, it is argued that physically reasonable, Delta-independent modes of crack extension can be extracted from oscillatory solutions. Knowledge of near-tip fields is used to determine the explicit Delta dependence of energy release rate parameters. Energy release rates are then defined that are separated from the oscillatory dependence on Delta. A modified VCCT using these energy release rate definitions is applied to results from finite element analyses, showing that Delta-independent modes of crack extension result. The modified technique has potential as a consistent method for extracting crack extension modes from numerical solutions. The Delta-independent modes extracted using this technique can also serve as guides for testing the convergence of finite element models. Direct applications of this work include the analysis of planar composite delamination problems, where plies or debonded laminates are modeled as in-plane orthotropic materials.

Beuth, Jack L.↗

Remote query measurement of pressure, fluid-flow velocity, and humidity using magnetoelastic thick-film sensors

Free-standing magnetoelastic thick-film sensors have a characteristic resonant frequency that can be determined by monitoring the magnetic flux emitted from the sensor in response to a time varying magnetic field. This property allows the sensors to be monitored remotely without the use of direct physical connections, such as wires, enabling measurement of environmental parameters from within sealed, opaque containers. In this work, we report on application of magnetoelastic sensors to measurement of atmospheric pressure, fluid-flow velocity, temperature, and mass load. Mass loading effects are demonstrated by fabrication of a remote query humidity sensor, made by coating the magnetoelastic thick film with a thin layer of solgel deposited Al2O3 that reversibly changes mass in response to humidity. c2000 Elsevier Science S.A. All rights reserved.

NASA Discipline Life Sciences Technologies↗