Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “viscoelasticity”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Nonlinear effects in thermal stress analysis of a solid propellant rocket motor

Direct characterization procedures were used to determine the relaxation modulus as a function of time, temperature, and state of strain. Using the quasi-elastic method of linearviscoelasticity, these properties were employed in a finite element computer code to analyze a thick-walled, nonlinear viscoelastic cylinder in the state of plane strain bonded to a thin (but stiff) elastic casing and subjected to slow thermal cooling. The viscoelastic solution is then expressed as a sequence of elastic finite element solutions. The strain-dependent character of the relaxation modulus is included by replacing the single relaxation curve used in the linear viscoelastic theory by a family of relaxation functions obtained at various strain levels. These functions may be regarded as a collection of stress histories or responses to specific loads (in this case, step strains) with which the cooldown solution is made to agree by iterations on the modulus and strain level.

Francis, E. C.↗
NASA NTRS1982

A time dependent theory of crazing behavior in polymers

The development of crazing is not only a function of stress, but also a function of time. Under a simple state of tension, a craze opening displacement is closely associated with the viscoelastic behavior of the original bulk polymer medium in which individual crazes initiate and develop. Within each craze region, molecular orientation takes place when conditions permit, and a new phase of rearranged molecules governs its local behavior. Based upon a time-dependent viscoelastic two-dimensional model, using a computer program the craze opening displacement field has been calculated, time-dependent craze length was also computed by taking into consideration the molecular orientation mechanism and large deformations in the craze region. Examples are given for simple viscoelastic media with simplified stress distributions. It is interesting to find out that the occurrence of crazing may be interpreted in terms of the stability or instability of the constitutive behavior of the bulk polymer.

Chern, S. S.↗
NASA NTRS1986

Prediction of the long-term creep compliance of general composite laminates

An accelerated viscoelastic characterization procedure for use with polymer-based composite materials is presented which employs short term test data obtained using unidirectional specimens to predict the long term viscoelastic behavior of general composite laminates. This procedure is here illustrated using the Schapery (1966, 1969) nonlinear theory as the required viscoelastic constitutive model, as well as classical lamination theory for the lamination scheme. The technique is applied to T300/5208 graphite/epoxy.

Tuttle, M. E.↗
NASA NTRS1988

Compression creep of filamentary composites

Axial and transverse strain fields induced in composite laminates subjected to compressive creep loading were compared for several types of laminate layups. Unidirectional graphite/epoxy as well as multi-directional graphite/epoxy and graphite/PEEK layups were studied. Specimens with and without holes were tested. The specimens were subjected to compressive creep loading for a 10-hour period. In-plane displacements were measured using moire interferometry. A computer based data reduction scheme was developed which reduces the whole-field displacement fields obtained using moire to whole-field strain contour maps. Only slight viscoelastic response was observed in matrix-dominated laminates, except for one test in which catastrophic specimen failure occurred after a 16-hour period. In this case the specimen response was a complex combination of both viscoelastic and fracture mechanisms. No viscoelastic effects were observed for fiber-dominated laminates over the 10-hour creep time used. The experimental results for specimens with holes were compared with results obtained using a finite-element analysis. The comparison between experiment and theory was generally good. Overall strain distributions were very well predicted. The finite element analysis typically predicted slightly higher strain values at the edge of the hole, and slightly lower strain values at positions removed from the hole, than were observed experimentally. It is hypothesized that these discrepancies are due to nonlinear material behavior at the hole edge, which were not accounted for during the finite-element analysis.

Graesser, D. L.↗
NASA NTRS1989

Calculating Dynamic Shear Moduli Of Polymers

Improved nonlinear mathematical model developed to fit experimental data on relaxation of stresses in viscoelastic materials. Developed specifically to study shear moduli of rubbery solid rocket propellants, also used to characterize polymers in general, other viscoelastic materials, and composites of viscoelastic materials. Facilitates and improves accuracy of analysis and numerical simulation of mechanical behavior of structural components (for example, tires) made of such materials.

Moore, Carleton J.↗
NASA NTRS1991

Composite propellant technology research: Mechanical property characterization

Proof for the existence of a single Poisson's ratio function in isotropic linear viscoelastic materials is presented. An in-depth discussion is given of three dimensional viscoelastic material properties and their relationships to linear isotropic and orthotropic viscoelastic materials. A discussion of the alternate invariant definition as used by Abaqus and how it relates to the form used by Dr. S. Peng is presented.

Bower, Mark V.↗
NASA NTRS1993

Predicting Failure Behaviors Of Thermoset Polymers

Useful predictions made from measurements on test-tube quantities. Mathematical model for viscoelastic failure behavior of epoxy or other thermoset polymer reduces considerably time and expense necessary to evaluate new polymer. One quickly obtains first-order estimates of viscoelastic properties, and decides more easily whether material rejected or good candidate for more thorough evaluation. Only properties required by model are equilibrium pressure vs. volume vs. temperature (PVT) response and linear viscoelastic shear and bulk moduli.

Bowles, Kenneth J.↗
NASA NTRS1996

The Extensional Rheology of Non-Newtonian Materials

The evolution of the transient extensional stresses in dilute and semi-dilute viscoelastic polymer solutions are measured with a filament stretching rheometer of a design similar to that first introduced by Sridhar, et al. The solutions are polystyrene-based (PS) Boger fluids that are stretched at constant strain rates ranging from 0.6 less than or equal to epsilon(0) less than or equal to 4s(exp -1) and to Hencky strains of epsilon greater than 4. The test fluids all strain harden and Trouton ratios exceeding 1000 are obtained at high strains. The experimental data strain hardens at lower strain levels than predicted by bead-spring FENE models. In addition to measuring the transient tensile stress growth, we also monitor the decay of the tensile viscoelastic stress difference in the fluid column following cessation of uniaxial elongation as a function of the total imposed Hencky strain and the strain rate. The extensional stresses initially decay very rapidly upon cessation of uniaxial elongation followed by a slower viscoelastic relaxation, and deviate significantly from FENE relaxation predictions. The relaxation at long times t is greater than or equal to 5 s, is compromised by gravitational draining leading to non-uniform filament profiles. For the most elastic fluids, partial decohension of the fluid filament from the endplates of the rheometer is observed in tests conducted at high strain rates. This elastic instability is initiated near the rigid endplate fixtures of the device and it results in the progressive breakup of the fluid column into individual threads or 'fibrils' with a regular azimuthal spacing. These fibrils elongate and bifurcate as the fluid sample is elongated further. Flow visualization experiments using a modified stretching device show that the instability develops as a consequence of an axisymmetry-breaking meniscus instability in the nonhomogeneous region of highly deformed fluid near the rigid endplate.

Spiegelberg, Stephen H.↗
NASA NTRS1999

Frequency-Dependent Viscosity of Xenon Near the Critical Point

We used a novel, overdamped oscillator aboard the Space Shuttle to measure the viscosity eta of xenon near its critical density rho(sub c), and temperature T(sub c). In microgravity, useful data were obtained within 0.1 mK of T(sub c), corresponding to a reduced temperature t = (T -T(sub c))/T(sub c) = 3 x 10(exp -7). The data extend two decades closer to T(sub c) than the best ground measurements, and they directly reveal the expected power-law behavior eta proportional to t(sup -(nu)z(sub eta)). Here nu is the correlation length exponent, and our result for the small viscosity exponent is z(sub eta) = 0.0690 +/- 0.0006. (All uncertainties are one standard uncertainty.) Our value for z(sub eta) depends only weakly on the form of the viscosity crossover function, and it agrees with the value 0.067 +/- 0.002 obtained from a recent two-loop perturbation expansion. The measurements spanned the frequency range 2 Hz less than or equal to f less than or equal to 12 Hz and revealed viscoelasticity when t less than or equal to 10(exp -1), further from T(sub c) than predicted. The viscoelasticity scales as Af(tau), where tau is the fluctuation-decay time. The fitted value of the viscoelastic time-scale parameter A is 2.0 +/- 0.3 times the result of a one-loop perturbation calculation. Near T(sub c), the xenon's calculated time constant for thermal diffusion exceeded days. Nevertheless, the viscosity results were independent of the xenon's temperature history, indicating that the density was kept near rho(sub c), by judicious choices of the temperature vs. time program. Deliberately bad choices led to large density inhomogeneities. At t greater than 10(exp -5), the xenon approached equilibrium much faster than expected, suggesting that convection driven by microgravity and by electric fields slowly stirred the sample.

Berg, Robert F.↗
NASA NTRS1999

Measuring Greenland Ice Mass Variation With Gravity Recovery and the Climate Experiment Gravity and GPS

The response of the Greenland ice sheet to climate change could significantly alter sea level. The ice sheet was much thicker at the last glacial maximum. To gain insight into the global change process and the future trend, it is important to evaluate the ice mass variation as a function of time and space. The Gravity Recovery and Climate Experiment (GRACE) mission to fly in 2001 for 5 years will measure gravity changes associated with the current ice variation and the solid earth's response to past variations. Our objective is to assess the separability of different change sources, accuracy and resolution in the mass variation determination by the new gravity data and possible Global Positioning System (GPS) bedrock uplift measurements. We use a reference parameter state that follows a dynamic ice model for current mass variation and a variant of the Tushingham and Peltier ICE-3G deglaciation model for historical deglaciation. The current linear trend is also assumed to have started 5 kyr ago. The Earth model is fixed as preliminary reference Earth model (PREM) with four viscoelastic layers. A discrete Bayesian inverse algorithm is developed employing an isotropic Gaussian a priori covariance function over the ice sheet and time. We use data noise predicted by the University of Texas and JPL for major GRACE error sources. A 2 mm/yr uplift uncertainty is assumed for GPS occupation time of 5 years. We then carry out covariance analysis and inverse simulation using GRACE geoid coefficients up to degree 180 in conjunction with a number of GPS uplift rates. Present-day ice mass variation and historical deglaciation are solved simultaneously over 146 grids of roughly 110 km x 110 km and with 6 time increments of 3 kyr each, along with a common starting epoch of the current trend. For present-day ice thickness change, the covariance analysis using GRACE geoid data alone results in a root mean square (RMS) posterior root variance of 2.6 cm/yr, with fairly large a priori uncertainties in the parameters and a Gaussian correlation length of 350 km. Simulated inverse can successfully recover most features in the reference present-day change. The RMS difference between them over the grids is 2.8 cm/yr. The RMS difference becomes 1.1 cm/yr when both are averaged with a half Gaussian wavelength of 150 km. With a fixed Earth model, GRACE alone can separate the geoid signals due to past and current load fairly well. Shown are the reference geoid signatures of direct and elastic effects of the current trend, the viscoelastic effect of the same trend starting from 5 kyr ago, the Post Glacial Rebound (PGR), and the predicted GRACE geoid error. The difference between the reference and inverse modeled total viscoelastic signatures is also shown. Although past and current ice mass variations are allowed the same spatial scale, their geoid signals have different spatial patterns. GPS data can contribute to the ice mass determination as well. Additional information is contained in the original.

Wu, Xiao-Ping↗
NASA NTRS2000

Wormhole Formation in RSRM Nozzle Joint Backfill

The RSRM nozzle uses a barrier of RTV rubber upstream of the nozzle O-ring seals. Post flight inspection of the RSRM nozzle continues to reveal occurrence of "wormholes" into the RTV backfill. The term "wormholes", sometimes called "gas paths", indicates a gas flow path not caused by pre-existing voids, but by a little-understood internal failure mode of the material during motor operation. Fundamental understanding of the mechanics of the RSRM nozzle joints during motor operation, nonlinear viscoelastic characterization of the RTV backfill material, identification of the conditions that predispose the RTV to form wormholes, and screening of candidate replacement materials is being pursued by a joint effort between Thiokol Propulsion, NASA, and the Army Propulsion & Structures Directorate at Redstone Arsenal. The performance of the RTV backfill in the joint is controlled by the joint environment. Joint movement, which applies a tension and shear load on the material, coupled with the introduction of high pressure gas in combination create an environment that exceeds the capability of the material to withstand the wormhole effect. Little data exists to evaluate why the material fails under the modeled joint conditions, so an effort to characterize and evaluate the material under these conditions was undertaken. Viscoelastic property data from characterization testing will anchor structural analysis models. Data over a range of temperatures, environmental pressures, and strain rates was used to develop a nonlinear viscoelastic model to predict material performance, develop criteria for replacement materials, and quantify material properties influencing wormhole growth. Three joint simulation analogs were developed to analyze and validate joint thermal barrier (backfill) material performance. Two exploratory tests focus on detection of wormhole failure under specific motor operating conditions. A "validation" test system provides data to "validate" computer models and predictions. Finally, two candidate replacement materials are being screened and "validated" using the developed test systems.

Stevens, J.↗
NASA NTRS2002

Thermo-Mechanical Analyses of Dynamically Loaded Rubber Cylinders

Thick rubber components are employed by the Army to carry large loads. In tanks, rubber covers road wheels and track systems to protect roadways. It is difficult for design engineers to simulate the details of the hysteretic heating for large strain viscoelastic deformations. In this study, an approximation to the viscoelastic energy dissipated per unit time is investigated for use in estimating mechanically induced viscoelastic heating. Coupled thermo-mechanical simulations of large cyclic deformations of rubber cylinders are presented. The cylinders are first compressed axially and then cyclically loaded about the compressed state. Details of the algorithm and some computational issues are discussed. The coupled analyses are conducted for tall and short rubber cylinders both with and without imbedded metal disks.

Johnson, Arthur R.↗
DOE OSTI2021

Phase-field modeling of rate-dependent fluid-driven fracture initiation and propagation

The rate-dependent behavior associated with deformation and fracturing of materials, such as natural rocks, poses significant challenges for modeling. In addition to the complications of the viscoelastic response, the speed of fracture propagation reflects micromechanical mechanisms in the fracture process zone(FPZ). In order to represent these complicated behaviors, a thermodynamically consistent, rate-dependent fracture model is required. Based on rigorous thermodynamic principles, we derive a rate-dependent phase-field mechanical model coupled with single-phase fluid flow in both the matrix and the fracture. The model is guaranteed to satisfy energy conservation during fracture propagation. Here, the system of equations is solved using the introduced solution procedure and a novel preconditioner that accounts for the complex fluid-structure interaction. The proposed phase-field model is tested against several benchmark problems on solid-fluid coupling, fluid-driven fracture propagation and rate-dependent viscoelastic deformation. The model serves as a strong basis for investigating rate-dependent fracturing experiments and for making predictions of material behaviors under new conditions.

42 ENGINEERING↗
DOE OSTI2021

Flow visualisation in real-size optical injectors of conventional, additised, and renewable gasoline blends

Research on renewable and alternative fuels is crucial for improving the energy and environmental efficiency of modern gasoline internal combustion engines. To highlight the influence of fuel rheological and thermodynamic properties on phase change and atomisation processes, three types of gasoline blends were tested. More specifically, the campaign comprised a reference gasoline, an ethanol/gasoline blend (10% v/v) representative of renewable fuels, and an additised gasoline sample treated with viscoelasticity-inducing agents. High-speed imaging of the transient two-phase flow field arising in the internal geometry and the near-nozzle spray region of gasoline injectors was performed employing Diffuse Backlight Illumination. The metallic body of a commercial injector was modified to fit transparent tips realising two nozzle layouts, namely a two-hole real size model resembling the Engine Combustion Network spray G injector and an enraged replica with an offset hole. Experiments were conducted at realistic operating conditions comprising an injection pressure of 100 bar and ambient pressures in the range of 0.1–6.0 bar to cover the entire range of chamber pressures prevailing in Gasoline Direct Injection engines. The action of viscoelastic additives was verified to have a suppressive effect on in-nozzle cavitation (6% reduction in cavitation extent) , while also enhancing spray atomisation at flash-boing conditions, in a manner resembling the more volatile gasoline/ethanol blends. Finally, persisting liquid ligaments were found to form after the end of injection for the additised sample, owing to the surfactant nature of the additives.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

A constitutive model for glass–ceramic materials

Glass–ceramics have received recent attention for use in glass–ceramic to metal hermetic seals. Due to their heterogeneous microstructure, these materials exhibit a number of advantageous responses over conventional glass based seals. Key amongst them is the possibility of a controllable thermal strain response and apparent coefficient of thermal expansion which may be used to minimize thermally induced residual stresses for aforementioned seals. Furthermore, these behaviors result from an inorganic glass matrix and variety of crystalline ceramic phases including silica polymorph(s) that may undergo reversible solid-to-solid transformations with associated inelastic strain. Correspondingly, these materials exhibit complex thermomechanical responses associated with multiple inelastic mechanisms (viscoelasticity and phase transformation). While modeling these behaviors is essential for developing and analyzing the corresponding applications, no such model exists. Therefore, in this work a three-dimensional continuum constitutive model for glass–ceramic materials combining these various inelastic mechanisms is developed via an internal state variable approach. A corresponding fully implicit three dimensional numerical formulation is also proposed and implemented. The model is used to simulate existing experiments and validate the proposed formalism. As an example, the simple seal problem of a glass–ceramic seal inside a concentric metal shell is explored. Finally, the impact of the cooling rates, viscoelastic shift factors, and inelastic strain on final residual stress state are all investigated and the differing contributions highlighted.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Soft sample deformation, damage and induced electromechanical property changes in contact- and tapping-mode atomic force microscopy

Probe-induced soft sample damage in atomic force microscopy (AFM), as well as the resulting alteration of local mechanical and electrical properties of the material are explored, specifically comparing contact-mode and intermittent-contact-mode imaging methods. In our experiments, performed on conductive polymer films, induced changes are present in contact-mode imaging while they are negligible or absent in tapping-mode imaging. Here, to understand this result, a viscoelastic parameter extraction is performed, which suggests that permanent sample deformation can readily occur for tip-sample interactions with a duration on the timescale of contact-mode interactions. Using the extracted viscoelastic parameters, a dynamic AFM simulation is conducted, which suggests that the material responds more elastically with reduced or absent sample damage in tapping-mode AFM, due to the higher rate of mechanical deformation and shorter timescales.

42 ENGINEERING↗
DOE OSTI2022

Application of Multiharmonic QCM-D for Detection of Plasmin at Hydrophobic Surfaces Modified by β-Casein

Plasmin protease plays an important role in many processes in living systems, including milk. Monitoring plasmin activity is important for control of the nutritional quality of milk and other dairy products. We designed a biosensor to detect the proteolytic activity of plasmin, using multiharmonic quartz crystal microbalance with dissipation (QCM-D). The β-casein immobilized on the hydrophobic surface of 1-dodecanethiol on the AT-cut quartz crystal was used to monitor plasmin activity. We demonstrated detection of plasmin in a concentration range of 0.1–20 nM, with the limit of detection about 0.13 ± 0.01 nM. The analysis of viscoelastic properties of the β-casein layer showed rapid changes of shear elasticity modulus, μ, and coefficient of viscosity, η, at plasmin sub-nanomolar concentrations, followed by modest changes at nanomolar concentrations, indicating multilayer architecture β-casein. A comparative analysis of viscoelastic properties of β-casein layers following plasmin and trypsin cleavage showed that the higher effect of trypsin was due to larger potential cleavage sites of β-casein.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Simulation of Deployable Composite Structures Based On Mechanics of Structure Genome

In this paper, a simulation method for analyzing deployable composite structures is presented. With a proper material model, effective plate/shell properties of the composites is obtained based on Mechanics of Structure Genome (MSG), and then implemented into a user subroutine UGENS for global structure simulation in ABAQUS. Column bending test (CBT) and composite boom and hub structure are studied for demonstration. A viscoelastic material model with direct integration implementation is adopted in this paper. CBT simulation shows good agreement with experiments during relaxation, while errors are observed when comparing residual deformation. This simulation can be potentially used as a calibration tool for material properties. After CBT simulation, a demonstrative model with a lenticular boom and the hub is created in ABAQUS. Complete process of flattening, coiling, stowage, deploying and recovery is simulated with the viscoelastic material model. Residual deformation of the boom is analyzed.

Multi-scale modeling↗