Macro–Micro-Coupled Simulations of Dilute Viscoelastic Fluids
Explore the source record for details and available documents.
Michael Cromer↗
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.
Shape recovery simulation of viscoelastic thin-ply composite coilable booms
Explore the source record for details and available documents.
Milinda Yapa Hamillage↗
Synergistic Effects of Physical Aging and Damage on Long-Term Behavior of Polymer Matrix Composites
The research consisted of two major parts, first modeling and simulation of the combined effects of aging and damage on polymer composites and secondly an experimental phase examining composite response at elevated temperatures, again activating both aging and damage. For the simulation, a damage model for polymeric composite laminates operating at elevated temperatures was developed. Viscoelastic behavior of the material is accounted for via the correspondence principle and a variational approach is adopted to compute the temporal stresses within the laminate. Also, the effect of physical aging on ply level stress and on overall laminate behavior is included. An important feature of the model is that damage evolution predictions for viscoelastic laminates can be made. This allows us to track the mechanical response of the laminate up to large load levels though within the confines of linear viscoelastic constitutive behavior. An experimental investigation of microcracking and physical aging effects in polymer matrix composites was also pursued. The goal of the study was to assess the impact of aging on damage accumulation, in ten-ns of microcracking, and the impact of damage on aging and viscoelastic behavior. The testing was performed both at room and elevated temperatures on [+/- 45/903](sub s) and [02/903](sub s) laminates, both containing a set of 90 deg plies centrally located to facilitate investigation of microcracking. Edge replication and X-ray-radiography were utilized to quantify damage. Sequenced creep tests were performed to characterize viscoelastic and aging parameters. Results indicate that while the aging times studied have limited ]Influence on damage evolution, elevated temperature and viscoelastic effects have a profound effect on the damage mode seen. Some results are counterintuitive, including the lower strain to failure for elevated temperature tests and the catastrophic failure mode observed for the [+/- 45/9O3](sub s), specimens. The fracture toughness for transverse cracks increases with increasing temperature for both systems: transverse cracking was completely absent prior to failure in [+/- 45/903](sub s), and was suppressed for [02/903](sub s). No significant effect of damage on aging or viscoelastic parameters was observed.
Brinson, L. Cate↗
HITEMP Material and Structural Optimization Technology Transfer
The feasibility of adding viscoelasticity and the Generalized Method of Cells (GMC) for micromechanical viscoelastic behavior into the commercial HyperSizer structural analysis and optimization code was investigated. The viscoelasticity methodology was developed in four steps. First, a simplified algorithm was devised to test the iterative time stepping method for simple one-dimensional multiple ply structures. Second, GMC code was made into a callable subroutine and incorporated into the one-dimensional code to test the accuracy and usability of the code. Third, the viscoelastic time-stepping and iterative scheme was incorporated into HyperSizer for homogeneous, isotropic viscoelastic materials. Finally, the GMC was included in a version of HyperSizer. MS Windows executable files implementing each of these steps is delivered with this report, as well as source code. The findings of this research are that both viscoelasticity and GMC are feasible and valuable additions to HyperSizer and that the door is open for more advanced nonlinear capability, such as viscoplasticity.
Collier, Craig S.↗
Postseismic rebound due to creep of the lower lithosphere and asthenosphere
Postseismic surface deformations are attributed to the inelastic flow of the subcrustal regions of the Earth following an earthquake. A multilayer representation of the Earth's rheological properties is used in conjunction with a finite element computational scheme to calculate time dependent displacements and strains subsequent to a strike slip earthquake. The deviatoric stress strain relations for the uppermost layer is assumed elastic. Lower layers are assumed to be, in order of increasing depth, a standard, linear, three element, viscoelastic solid; a linear viscoelastic fluid; and another elastic solid. Physically these layers correspond to the upper lithosphere, lower lithosphere, asthenosphere, and lower mantle, respectively. Elastic dilatational properties are assumed throughout. Appreciable postseismic displacements, possibly approaching meters, for large earthquakes, arise from viscoelastic relaxation following sudden coseismic slip. Furthermore, compared to the simpler case of an elastic lithosphere over a viscoelastic asthenosphere and the near fault postseismic shear, strain is increased, by a factor of two or more in some cases, by the presence of a viscoelastic lower lithosphere. Also, the duration of postseismic straining is increased by the relatively slow relaxation of this layer.
Cohen, S. C.↗
A multilayer model of time dependent deformation following an earthquake on a strike-slip fault
A multilayer model of the Earth to calculate finite element of time dependent deformation and stress following an earthquake on a strike slip fault is discussed. The model involves shear properties of an elastic upper lithosphere, a standard viscoelastic linear solid lower lithosphere, a Maxwell viscoelastic asthenosphere and an elastic mesosphere. Systematic variations of fault and layer depths and comparisons with simpler elastic lithosphere over viscoelastic asthenosphere calculations are analyzed. Both the creep of the lower lithosphere and astenosphere contribute to the postseismic deformation. The magnitude of the deformation is enhanced by a short distance between the bottom of the fault (slip zone) and the top of the creep region but is less sensitive to the thickness of the creeping layer. Postseismic restressing is increased as the lower lithosphere becomes more viscoelastic, but the tendency for the width of the restressed zone to growth with time is retarded.
Cohen, S. C.↗
Thermoviscoelastic characterization and prediction of Kevlar/epoxy composite laminates
The thermoviscoelastic characterization of Kevlar 49/Fiberite 7714A epoxy composite lamina and the development of a numerical procedure to predict the viscoelastic response of any general laminate constructed from the same material were studied. The four orthotropic material properties, S sub 11, S sub 12, S sub 22, and S sub 66, were characterized by 20 minute static creep tests on unidirectional (0) sub 8, (10) sub 8, and (90) sub 16 lamina specimens. The Time-Temperature Superposition-Principle (TTSP) was used successfully to accelerate the characterization process. A nonlinear constitutive model was developed to describe the stress dependent viscoelastic response for each of the material properties. A numerical procedure to predict long term laminate properties from lamina properties (obtained experimentally) was developed. Numerical instabilities and time constraints associated with viscoelastic numerical techniques were discussed and solved. The numerical procedure was incorporated into a user friendly microcomputer program called Viscoelastic Composite Analysis Program (VCAP), which is available for IBM PC type computers. The program was designed for ease of use. The final phase involved testing actual laminates constructed from the characterized material, Kevlar/epoxy, at various temperatures and load level for 4 to 5 weeks. These results were compared with the VCAP program predictions to verify the testing procedure and to check the numerical procedure used in the program. The actual tests and predictions agreed for all test cases which included 1, 2, 3, and 4 fiber direction laminates.
Gramoll, K. C.↗
Composite passive damping struts for large precision structures
In the field of viscoelastic dampers, a new strut design comprises a viscoelastic material sandwiched between multiple layers, some of which layers bear and dampen load force. In one embodiment, the layers are composite plies of opposing orientation. In another embodiment, the strut utilizes a viscoelastic layer sandwiched between V-shaped composite plies. In a third embodiment, a viscoelastic layer is sandwiched between sine-shaped plies. Strut strength is equal to or greater than conventional aluminum struts due to the unique high interlaminar shear ply design.
Dolgin, Benjamin P.↗
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.↗
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.↗
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.↗
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.↗
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.↗
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.↗
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.↗
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.↗
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.↗
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↗