Search NASA⌕ Search

SEARCH · Search NASA

Results for “viscoelastic modeling”

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 37 records · Page 2

Time-domain approach for the transient responses in stratified viscoelastic Earth models

We have developed the numerical algorithm for the computation of transient viscoelastic responses in the time domain for a radially stratified Earth model. Stratifications in both the elastic parameters and the viscosity profile have been considered. The particular viscosity profile employed has a viscosity maximum with a constrast of O(100) in the mid lower mantle. The distribution of relaxation times reveals the presence of a continuous spectrum situated between O(100) and O(exp 4) years. The principal mode is embedded within this continuous spectrum. From this initial-value approach we have found that for the low degree harmonics the non-modal contributions are comparable to the modal contributions. For this viscosity model the differences between the time-domain and normal-mode results are found to decrease strongly with increasing angular order. These calculations also show that a time-dependent effective relaxation time can be defined, which can be bounded by the relaxation times of the principal modes.

Hanyk, L.↗

Rheological effects on friction in elastohydrodynamic lubrication

An analytical and experimental investigation is presented of the friction in a rolling and sliding elastohydrodynamic lubricated contact. The rheological behavior of the lubricant is described in terms of two viscoelastic models. These models represent the separate effects of non-Newtonian behavior and the transient response of the fluid. A unified description of the non-Newtonian shear rate dependence of the viscosity is presented as a new hyperbolic liquid model. The transient response of viscosity, following the rapid pressure rise encountered in the contact, is described by a compressional viscoelastic model of the volume response of a liquid to an applied pressure step. The resulting momentum and energy equations are solved by an iterative numerical technique, and a friction coefficient is calculated. The experimental study was performed, with two synthetic paraffinic lubricants, to verify the friction predictions of the analysis. The values of friction coefficient from theory and experiment are in close agreement.

Trachman, E. G.↗

Nonlinear viscoelastic characterization of thin polyethylene film

In order to understand the state of stress and strain in a typical balloon fabricated from thin polyethylene film, experiment data in the literature reviewed. It was determined that the film behaves as a nonlinear viscoelasticity material and should be characterized accordingly. A simple uniaxial, nonlinear viscoelastic model was developed for predicting stress given a certain strain history. The simple model showed good qualitative agreement with results of constant rate, uniaxial accurately predicting stresses for cyclic strain histories typical of balloon flights. A program was outlined which will result in the development of a more complex nonlinear viscoelastic model.

Wilbeck, J. S.↗

A Nonlinear Viscoelastic-Viscoplastic Constitutive Model for Epoxy Polymers

The objective of this paper is to develop a nonlinear viscoelastic–viscoplastic constitutive model for epoxy polymers. The classic nonlinear viscoelasticity model is reformulated to yield a closed-form incremental constitutive relation, which relates the stress increments to the viscoelastic strain increments. A viscoplasticity model, which consists of the Drucker–Prager yield function, nonlinear isotropic and kinematic hardening laws, and the Perzyna viscosity function, is subsequently developed. The nonlinear viscoelasticity and the vis-coplasticity models are then implemented in a radial return algorithm. A closed-form algorithmic tangent operator is derived to facilitate model implementation. The nonlinear viscoelasticity model is validated by reproducing the test data for polymethyl methacrylate(PMMA). The present constitutive model’s capabilities are demonstrated through modeling viscoelastic–viscoplastic PMMA loaded at different strain rates. The mechanics of structure genome-based micromechanics approach, embedded with the present constitutive model, is used to homogenize a unidirectional fiber-reinforced composite with a PMMA matrix, subjected to uniaxial and shear loading at different strain rates. The viscoelastic characterization of PMT-F7 epoxy is evaluated by comparing the predicted creep–recovery responses by the present constitutive model with the experimental ones.

Liang Zhang↗

Development of a solid propellant viscoelastic dynamic model

The results of a one year study to develop a dynamic response model for the Space Shuttle Solid Rocket Motor (SRM) propellant are presented. An extensive literature survey was conducted, from which it was concluded that the only significant variables affecting the dynamic response of the SRM propellant are temperature and frequency. Based on this study, and experimental data on propellants related to the SRM propellant, a dynamic constitutive model was developed in the form of a simple power law with temperature incorporated in the form of a modified power law. A computer program was generated which performs a least-squares curve-fit of laboratory data to determine the model parameters and it calculates dynamic moduli at any desired temperature and frequency. Additional studies investigated dynamic scaling laws and the extent of coupling between the SRM propellant and motor cases. It was found, in agreement with other investigations, that the propellant provides all of the mass and damping characteristics whereas the case provides all of the stiffness.

Hufferd, W. L.↗

A K-BKZ Formulation for Soft-Tissue Viscoelasticity

A viscoelastic model of the K-BKZ (Kaye 1962; Bernstein et al. 1963) type is developed for isotropic biological tissues, and applied to the fat pad of the human heel. To facilitate this pursuit, a class of elastic solids is introduced through a novel strain-energy function whose elements possess strong ellipticity, and therefore lead to stable material models. The standard fractional-order viscoelastic (FOV) solid is used to arrive at the overall elastic/viscoelastic structure of the model, while the elastic potential via the K-BKZ hypothesis is used to arrive at the tensorial structure of the model. Candidate sets of functions are proposed for the elastic and viscoelastic material functions present in the model, including a regularized fractional derivative that was determined to be the best. The Akaike information criterion (AIC) is advocated for performing multi-model inference, enabling an objective selection of the best material function from within a candidate set.

Freed, Alan D.↗

A viscoelastic higher-order beam finite element

A viscoelastic internal variable constitutive theory is applied to a higher-order elastic beam theory and finite element formulation. The behavior of the viscous material in the beam is approximately modeled as a Maxwell solid. The finite element formulation requires additional sets of nodal variables for each relaxation time constant needed by the Maxwell solid. Recent developments in modeling viscoelastic material behavior with strain variables that are conjugate to the elastic strain measures are combined with advances in modeling through-the-thickness stresses and strains in thick beams. The result is a viscous thick-beam finite element that possesses superior characteristics for transient analysis since its nodal viscous forces are not linearly dependent an the nodal velocities, which is the case when damping matrices are used. Instead, the nodal viscous forces are directly dependent on the material's relaxation spectrum and the history of the nodal variables through a differential form of the constitutive law for a Maxwell solid. The thick beam quasistatic analysis is explored herein as a first step towards developing more complex viscoelastic models for thick plates and shells, and for dynamic analyses. The internal variable constitutive theory is derived directly from the Boltzmann superposition theorem. The mechanical strains and the conjugate internal strains are shown to be related through a system of first-order, ordinary differential equations. The total time-dependent stress is the superposition of its elastic and viscous components. Equations of motion for the solid are derived from the virtual work principle using the total time-dependent stress. Numerical examples for the problems of relaxation, creep, and cyclic creep are carried out for a beam made from an orthotropic Maxwell solid.

Johnson, Arthur R.↗

Postseismic viscoelastic deformation and stress. Part 2: Stress theory and computation; dependence of displacement, strain, and stress on fault parameters

A viscoelastic model for deformation and stress associated with earthquakes is reported. The model consists of a rectangular dislocation (strike slip fault) in a viscoelastic layer (lithosphere) lying over a viscoelastic half space (asthenosphere). The time dependent surface stresses are analyzed. The model predicts that near the fault a significant fraction of the stress that was reduced during the earthquake is recovered by viscoelastic softening of the lithosphere. By contrast, the strain shows very little change near the fault. The model also predicts that the stress changes associated with asthenospheric flow extend over a broader region than those associated with lithospheric relaxation even though the peak value is less. The dependence of the displacements, stresses on fault parameters studied. Peak values of strain and stress drop increase with increasing fault height and decrease with fault depth. Under many circumstances postseismic strains and stresses show an increase with decreasing depth to the lithosphere-asthenosphere boundary. Values of the strain and stress at distant points from the fault increase with fault area but are relatively insensitive to fault depth.

Cohen, S. C.↗

The analytical representation of viscoelastic material properties using optimization techniques

This report presents a technique to model viscoelastic material properties with a function of the form of the Prony series. Generally, the method employed to determine the function constants requires assuming values for the exponential constants of the function and then resolving the remaining constants through linear least-squares techniques. The technique presented here allows all the constants to be analytically determined through optimization techniques. This technique is employed in a computer program named PRONY and makes use of commercially available optimization tool developed by VMA Engineering, Inc. The PRONY program was utilized to compare the technique against previously determined models for solid rocket motor TP-H1148 propellant and V747-75 Viton fluoroelastomer. In both cases, the optimization technique generated functions that modeled the test data with at least an order of magnitude better correlation. This technique has demonstrated the capability to use small or large data sets and to use data sets that have uniformly or nonuniformly spaced data pairs. The reduction of experimental data to accurate mathematical models is a vital part of most scientific and engineering research. This technique of regression through optimization can be applied to other mathematical models that are difficult to fit to experimental data through traditional regression techniques.

Hill, S. A.↗

A General Reversible Hereditary Constitutive Model: Application to a Titanium Alloy - Part 2

Given the mathematical framework and specific viscoelastic model in Part 1 our primary goal in this second part is focused on model characterization and assessment for the specific titanium alloy, TIMETAL 21S. The model is motivated by experimental evidence suggesting the presence of significant rate/time effects in the so-called quasilinear, reversible, material response range. An explanation of the various experiments performed and their corresponding results are also included. Finally, model correlations and predictions are presented for a wide temperature range.

Arnold, S. M.↗

Comparison of observed rheological properties of hard wheat flour dough with predictions of the Giesekus-Leonov, White-Metzner and Phan-Thien Tanner models

The measured rheological behavior of hard wheat flour dough was predicted using three nonlinear differential viscoelastic models. The Phan-Thien Tanner model gave good zero shear viscosity prediction, but overpredicted the shear viscosity at higher shear rates and the transient and extensional properties. The Giesekus-Leonov model gave similar predictions to the Phan-Thien Tanner model, but the extensional viscosity prediction showed extension thickening. Using high values of the mobility factor, extension thinning behavior was observed but the predictions were not satisfactory. The White-Metzner model gave good predictions of the steady shear viscosity and the first normal stress coefficient but it was unable to predict the uniaxial extensional viscosity as it exhibited asymptotic behavior in the tested extensional rates. It also predicted the transient shear properties with moderate accuracy in the transient phase, but very well at higher times, compared to the Phan-Thien Tanner model and the Giesekus-Leonov model. None of the models predicted all observed data consistently well. Overall the White-Metzner model appeared to make the best predictions of all the observed data.

Non-NASA Center↗

Finite-Strain Fractional-Order Viscoelastic (FOV) Material Models and Numerical Methods for Solving Them

Fraction-order viscoelastic (FOV) material models have been proposed and studied in 1D since the 1930's, and were extended into three dimensions in the 1970's under the assumption of infinitesimal straining. It was not until 1997 that Drozdov introduced the first finite-strain FOV constitutive equations. In our presentation, we shall continue in this tradition by extending the standard, FOV, fluid and solid, material models introduced in 1971 by Caputo and Mainardi into 3D constitutive formula applicable for finite-strain analyses. To achieve this, we generalize both the convected and co-rotational derivatives of tensor fields to fractional order. This is accomplished by defining them first as body tensor fields and then mapping them into space as objective Cartesian tensor fields. Constitutive equations are constructed using both variants for fractional rate, and their responses are contrasted in simple shear. After five years of research and development, we now possess a basic suite of numerical tools necessary to study finite-strain FOV constitutive equations and their iterative refinement into a mature collection of material models. Numerical methods still need to be developed for efficiently solving fraction al-order integrals, derivatives, and differential equations in a finite element setting where such constitutive formulae would need to be solved at each Gauss point in each element of a finite model, which can number into the millions in today's analysis.

Freed, Alan D.↗

Transients in Pacific/North American Plate Boundary Deformation: Synthesis and Modeling of GPS and Borehole Strain Observations

This is the Final Technical Report on research conducted between 1 June 1997 and 14 September 2001 entitled "Transients in Pacific/North American plate boundary deformation: Synthesis and modeling of GPS and borehole strain observations." As the project title implies, our effort involved a geodetic study of strain transients, i.e., temporal variations in deformation rates, that occur within plate boundary zones and their relationship to earthquakes and plate motions. Important transients occur during and following large earthquakes, and there are also strain transients not apparently associated with earthquakes. A particularly intriguing class of transients, for which there is a modest but growing list of examples, are preseismic anomalies. Such earthquake precursors, if further documented and understood, would have obvious importance for earthquake hazard mitigation. Because the timescales for these diverse transients range over at least 6 orders of magnitude (minutes to years), no single geodetic technique is optimum. We therefore undertook a systematic synthesis of Global Positioning Satellite (GPS) and borehole strainmeter data in three areas in California where there are adequate numbers of both types of instruments (or their equivalent): the San Francisco Bay region (within the Bay Area Regional Deformation network), southern California (within the Southern California Integrated GPS Network), and Parkfield (where a two-color laser system provides a proxy for continuous GPS measurements). An integral component of our study was the elucidation of the physical mechanisms by which such transients occur and propagate. We therefore initiated the development of multiple forward models, using two independent approaches. In the first, we explored the response to specified earthquake slip in viscoelastic models that incorporated failure criteria and the geometry of major faults in California. In the second approach, we examined the dynamical response of a complex rheological medium to the application of a far-field stress imposed by plate motions. The forward models were used both to gain insight into the range of strain transients to be expected under different assumed mechanical conditions and to develop representations for strain fields that allow GPS, borehole, and other strain data to be combined in a self-consistent, yet well-determined, manner. The models also provided a basis for hypothesis testing, by which data from a strain transient well characterized by GPS and borehole observations were utilized to distinguish among competing candidates for the causative physical mechanism and the governing physical characteristics. During the three years of this project, continued to a fourth year through a no-cost extension of the grant, we published 14 papers and presented or co-authored 37 papers at national scientific meetings.

Solomon, Sean C.↗

A nonlinear viscoelastic characterization of graphite/epoxy composites

A nonlinear viscoelastic model based on the procedure of Findley et al., (1948) has been used to accurately represent the creep of several graphite/epoxy composites. Applying this approach to unidirectional 0, 90 and 10 deg off-axis tensile specimens, the viscoelastic response of a lamina can be characterized. The resulting lamina model has been useful for representing the behavior of a lamina in a numerical procedure to predict creep and delayed failures of general laminates. Also, independently, the Findley procedure has been used to characterize the nonlinear viscoelastic behavior of general laminates.

Dillard, D. A.↗

Design and Analysis of Molded Elastomer Seals

Modeling Materials: a) Hyperelastic model used for room temperature and higher, but has limitations at predicting response at different compressive rates and at low temperatures. b) Viscoelastic model is rarely used, and needs more research and application. Test and Characterization of Materials: a) Better material characterization using biaxial tension vs simple compression testing modes is debatable. b) Testing at high and low temperature is challenging.

Yetter, Paul↗

Sea Level Rise in the Samoan Islands Escalated by Viscoelastic Relaxation After the 2009 Samoa-Tonga Earthquake

The Samoan islands are an archipelago hosting a quarter million people mostly residing inthree major islands, Savai'i and Upolu (Samoa), and Tutuila (American Samoa). The islands haveexperienced sea level rise by 2-3 mm/year during the last half century. The rate, however, has dramaticallyincreased following the Mw 8.1 SamoaTonga earthquake doublet (megathrust + normal faulting) inSeptember 2009. Since the earthquake, we found largescale gravity increase (0.5 Gal/year) around theislands and ongoing subsidence (8-16 mm/year) of the islands from our analysis of Gravity Recovery AndClimate Experiment gravity and GPS displacement data. The postseismic horizontal displacement is faster inSamoa, while the postseismic subsidence rate is considerably larger in American Samoa. The analysis oflocal tide gauge records and satellite altimeter data also identified that the relative sea level rise becomesfaster by 7-9 mm/year in American Samoa than Samoa. A simple viscoelastic model with a Maxwellviscosity of 2–3×10(exp 18) Pa s for the asthenosphere explained postseismic deformation at nearby GPS sites aswell as Gravity Recovery And Climate Experiment gravity change. It is found that the constructiveinterference of viscoelastic relaxation from both megathrust and normal faulting has intensified thepostseismic subsidence at American Samoa, causing ~5 times faster sea level rise than the global average.Our model indicates that this trend is likely to continue for decades and result in sea level rise of 30-40 cm,which is independent of and in addition to anticipated climaterelated sea level rise. It will worsen coastalflooding on the islands leading to regular nuisance flooding.

Han, Shin-Chan↗