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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Tidal dissipation in a viscoelastic planet

Tidal dissipation is examined using Maxwell standard liner solid (SLS), and Kelvin-Voigt models, and viscosity parameters are derived from the models that yield the amount of dissipation previously calculated for a moon model with QW = 100 in a hypothetical orbit closer to the earth. The relevance of these models is then assessed for simulating planetary tidal responses. Viscosities of 10 exp 14 and 10 ex 18 Pa s for the Kelvin-Voigt and Maxwell rheologies, respectively, are needed to match the dissipation rate calculated using the Q approach with a quality factor = 100. The SLS model requires a short time viscosity of 3 x 10 exp 17 Pa s to match the Q = 100 dissipation rate independent of the model's relaxation strength. Since Q = 100 is considered a representative value for the interiors of terrestrial planets, it is proposed that derived viscosities should characterize planetary materials. However, it is shown that neither the Kelvin-Voigt nor the SLS models simulate the behavior of real planetary materials on long time scales. The Maxwell model, by contrast, behaves realistically on both long and short time scales. The inferred Maxwell viscosity, corresponding to the time scale of days, is several times smaller than the longer time scale (greater than or equal to 10 exp 14 years) viscosity of the earth's mantle.

Ross, M.↗

NASA NTRS1984

Constraints on the origin of viscoelastic bodies

While the oceans undoubtly contribute to energy dissipation, it seems implausible that such a small amount (0.02%) of the Earth's mass could play such a dominant role in lunar orbital evolution. An alternate hypothesis assumes that most of the dissipation is associated with solid-body tides. This mechanism is also capable of varying significantly over time for viscoelastic bodies. In such an analysis of orbital evolution, the bodies are modeled as incompressible, Kelvin-Voight solids. While this material is simplistic, it does provide an analytical characterization of solid-body dissipation and the resulting tidal moment. The orbital history of a two body system is traced back in time to develop constraints on the satellite's origin. These constraints suggest the Moon appeared in a significantly inclined orbit at a distance greater than ten Earth radii. This distance is outside of the Roche limit and would appear to alleviate some difficulties associated with capture hypotheses. Small eccentricity at close approach suggests an origin by accretion.

Vanarsdale, W. E.↗

NASA NTRS2003

Crustal Deformation in Southcentral Alaska: The 1964 Prince William Sound Earthquake Subduction Zone

This article, for Advances in Geophysics, is a summary of crustal deformation studies in southcentral Alaska. In 1964, southcentral Alaska was struck by the largest earthquake (moment magnitude 9.2) occurring in historical times in North America and the second largest earthquake occurring in the world during the past century. Conventional and space-based geodetic measurements have revealed a complex temporal-spatial pattern of crustal movement. Numerical models suggest that ongoing convergence between the North America and Pacific Plates, viscoelastic rebound, aseismic creep along the tectonic plate interface, and variable plate coupling all play important roles in controlling both the surface and subsurface movements. The geodetic data sets include tide-gauge observations that in some cases provide records back to the decades preceding the earthquake, leveling data that span a few decades around the earthquake, VLBI data from the late 1980s, and GPS data since the mid-1990s. Geologic data provide additional estimates of vertical movements and a chronology of large seismic events. Some of the important features that are revealed by the ensemble of studies that are reviewed in this paper include: (1) Crustal uplift in the region that subsided by up 2 m at the time of the earthquake is as much as 1 m since the earthquake. In the Turnagain Arm and Kenai Peninsula regions of southcentral Alaska, uplift rates in the immediate aftermath of the earthquake reached 150 mm/yr , but this rapid uplift decayed rapidly after the first few years following the earthquake. (2) At some other locales, notably those away the middle of the coseismic rupture zone, postseismic uplift rates were initially slower but the rates decay over a longer time interval. At Kodiak Island, for example, the uplift rates have been decreasing at a rate of about 7mm/yr per decade. At yet other locations, the uplift rates have shown little time dependence so far, but are thought not to be sustainable throughout the several hundred year recurrence time for great earthquake. The nearly 10 mm/yr uplift rate at Seldovia on the Kenai Peninsula is an example.

Cohen, Steven C.↗

NASA NTRS1993

State of stress, faulting, and eruption characteristics of large volcanoes on Mars

The formation of a large volcano loads the underlying lithospheric plate and can lead to lithospheric flexure and faulting. In turn, lithospheric stresses affect the stress field beneath and within the volcanic edifice and can influence magma transport. Modeling the interaction of these processes is crucial to an understanding of the history of eruption characteristics and tectonic deformation of large volcanoes. We develop models of time-dependent stress and deformation of the Tharsis volcanoes on Mars. A finite element code is used that simulates viscoelastic flow in the mantle and elastic plate flexural behavior. We calculate stresses and displacements due to a volcano-shaped load emplaced on an elastic plate. Models variously incorporate growth of the volcanic load with time and a detachment between volcano and lithosphere. The models illustrate the manner in which time-dependent stresses induced by lithospheric plate flexure beneath the volcanic load may affect eruption histories, and the derived stress fields can be related to tectonic features on and surrounding martian volcanoes.

Mcgovern, Patrick J.↗

NASA NTRS1997

A General Reversible Hereditary Constitutive Model: Theoretical Developments - Part 1

Using an internal-variable formalism as a starting point, we describe the viscoelastic extension of a previously-developed viscoplasticity formulation of the complete potential structure type. It is mainly motivated by experimental evidence for the presence of rate/time effects in the so-called quasilinear, reversible, material response range. Several possible generalizations are described, in the general format of hereditary-integral representations for non-equilibrium, stress-type, state variables, both for isotropic as well as anisotropic materials. In particular, thorough discussions are given on the important issues of thermodynamic admissibility requirements for such general descriptions, resulting in a set of explicit mathematical constraints on the associated kernel (relaxation and creep compliance) functions. In addition, a number of explicit, integrated forms are derived, under stress and strain control to facilitate the parametric and qualitative response characteristic studies reported here, as well as to help identify critical factors in the actual experimental characterizations from test data that will be reported in Part II.

Saleeb, A. F.↗

NASA NTRS1992

Direct simulation of polymer drag reduction in free shear flows and vortex dipoles

One of the most efficient techniques for drag reduction is the injection of polymers near a wall which can achieve a reduction in drag up to 80 percent. Several experimental observations tend to indicate that polymers modify the turbulence structures within the buffer layer and show that the changes consist of a weakening of the strength of the streamwise vortices. In this paper, we investigate the effects of viscoelasticity on two different types of flows: the vortex dipole impinging walls to model streamwise vortices in a turbulent boundary layer and the mixing layer that represents free shear flows. For this purpose, we examined three different rheological models: the Oldroyd-B model, the Jeffrey's corotational model, and the FENE-P model.

Orlandi, P.↗

NASA NTRS2002

A Thermodynamic Theory of Solid Viscoelasticity: Nonlinear Thermo-viscoelasticity - Part II:

This paper, second in the series of three papers, develops a general, nonlinear, non-isothermal, compressible theory for finite rubber viscoelasticity and specifies it in a form convenient for solving problems important to the rubber, tire, automobile, and air-space industries, among others. Based on the quasi-linear approach of non-equilibrium thermodynamics, a general nonlinear theory of differential type has been developed for arbitrary non-isothermal deformations of viscoelastic solids. In this theory, the constitutive equations were presented as the sum of a rubber elastic (equilibrium) and a liquid type viscoelastic (non-equilibrium) terms. These equations have then been simplified using several modeling and simplicity arguments.

Freed, Alan D.↗

NASA NTRS2003

Coulomb Stress Accumulation along the San Andreas Fault System

Stress accumulation rates along the primary segments of the San Andreas Fault system are computed using a three-dimensional (3-D) elastic half-space model with realistic fault geometry. The model is developed in the Fourier domain by solving for the response of an elastic half-space due to a point vector body force and analytically integrating the force from a locking depth to infinite depth. This approach is then applied to the San Andreas Fault system using published slip rates along 18 major fault strands of the fault zone. GPS-derived horizontal velocity measurements spanning the entire 1700 x 200 km region are then used to solve for apparent locking depth along each primary fault segment. This simple model fits remarkably well (2.43 mm/yr RMS misfit), although some discrepancies occur in the Eastern California Shear Zone. The model also predicts vertical uplift and subsidence rates that are in agreement with independent geologic and geodetic estimates. In addition, shear and normal stresses along the major fault strands are used to compute Coulomb stress accumulation rate. As a result, we find earthquake recurrence intervals along the San Andreas Fault system to be inversely proportional to Coulomb stress accumulation rate, in agreement with typical coseismic stress drops of 1 - 10 MPa. This 3-D deformation model can ultimately be extended to include both time-dependent forcing and viscoelastic response.

Smith, Bridget↗

NASA NTRS1988

Energy dissipation in a rolling aircraft tire

The project is extending an existing finite element tire model to calculate the energy dissipation in a free-rolling aircraft tire and temperature buildup in the tire carcass. The model will provide a means of calculating the influence of tire design on the distribution of tire temperature. Current focus is on energy loss measurements of aircraft tire material. The feasibility of taking test specimens directly from the tire carcass for measurements of viscoelastic properties was demonstrated. The interaction of temperature and frequency effects on material loss properties was studied. The tire model was extended to calculate the cyclic energy change in a tire during rolling under load. Input data representing the 40 by 14 aircraft tire whose material loss properties were measured are being used.

Tielking, John T.↗

NASA NTRS2022

Tidal Insights into Rocky and Icy Bodies: An Introduction and Overview

Solid body tides provide key information on the interior structure, evolution, and origin of the planetary bodies. Our Solar system harbours a very diverse population of planetary bodies, including those composed of rock, ice, gas, or a mixture of all. While a rich arsenal of geophysical methods has been developed over several years to infer knowledge about the interior of the Earth, the inventory of tools to investigate the interiors of other Solar-system bodies remains limited. With seismic data only available for the Earth, the Moon, and Mars, geodetic measurements, including the observation of the tidal response, have become especially valuable and therefore, has played an important role in understanding the interior and history of several Solar system bodies. To use tidal response measurements as a means to obtain constraints on the interior structure of planetary bodies, appropriate understanding of the viscoelastic reaction of the materials from which the planets are formed is needed. Here, we review the fundamental aspects of the tidal modeling and the information on the present-day interior properties and evolution of several planets and moons based on studying their tidal response. We begin with an outline of the theory of viscoelasticity and tidal response. Next, we proceed by discussing the information on the tidal response and the inferred structure of Mercury, Venus, Mars and its moons, the Moon, and the largest satellites of giant planets, obtained from the analysis of the data that has been provided by space missions. We also summarise the upcoming possibilities offered by the currently planned missions.

Tides↗

NASA NTRSDate not supplied

Experimental Characterization of the Dimensional Stability of Deployable Composite Booms During Stowage

Thin-ply carbon fiber reinforced polymer materials offer the opportunity to create long thin-shell booms with high strength to weight ratios. These booms can be rolled up and stowed in small volumes to be later deployed. These features make them ideal for use in gossamer solar sails as well as other deployable space structures. It has been observed that stowing the booms in rolled-up configurations can cause time-dependent deformations. A comprehensive understanding of the detrimental effects of stowing the coiled booms is necessary to implement these structures in spaceflight missions. During stowage inside the spacecraft, the boom is subjected to a constant deformation/strain that causes it to relax over time, measurable by a decrease in stress. Using large deformation bending tests, the complete fold-stow-unfold-recover cycle of thin-laminate candidates for the rollable booms is characterized. Further bending and stowage testing is done at the boom level. The results from this test campaign will be used to evaluate the dimensional stability of deployable booms made from thin-ply composites that exhibit visco-elasto-plastic behavior and calibrate and validate numerical finite element models that include visco-elasto-plastic material models.

Experimental Characterization↗

NASA NTRS1995

Expendable launch vehicle studies

Analytical support studies of expendable launch vehicles concentrate on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed and linearized. The information from LeRC about the grids, masses, and modes is incorporated into the model. The eigenvalues of the plant are analyzed for several modeling factors: utilizing diagonal mass matrix, uniform beam assumption, inclusion of aerodynamics, and the interaction between the aerodynamics and the flexible bending motion. Preliminary PID, LQR, and LQG control designs with sensor and actuator dynamics for this system and simulations are also conducted. The initial analysis for comparison of PD (proportional-derivative) and full state feedback LQR Linear quadratic regulator) shows that the split weighted LQR controller has better performance than that of the PD. In order to meet both the performance and robustness requirements, the H(sub infinity) robust controller for the expendable launch vehicle is developed. The simulation indicates that both the performance and robustness of the H(sub infinity) controller are better than that for the PID and LQG controllers. The modelling and analysis support studies team has continued development of methodology, using eigensensitivity analysis, to solve three classes of discrete eigenvalue equations. In the first class, the matrix elements are non-linear functions of the eigenvector. All non-linear periodic motion can be cast in this form. Here the eigenvector is comprised of the coefficients of complete basis functions spanning the response space and the eigenvalue is the frequency. The second class of eigenvalue problems studied is the quadratic eigenvalue problem. Solutions for linear viscously damped structures or viscoelastic structures can be reduced to this form. Particular attention is paid to Maxwell and Kelvin models. The third class of problems consists of linear eigenvalue problems in which the elements of the mass and stiffness matrices are stochastic. dynamic structural response for which the parameters are given by probabilistic distribution functions, rather than deterministic values, can be cast in this form. Solutions for several problems in each class will be presented.

Bainum, Peter M.↗

NASA NTRS1988

Advanced composite materials for precision segmented reflectors

The objective in the NASA Precision Segmented Reflector (PSR) project is to develop new composite material concepts for highly stable and durable reflectors with precision surfaces. The project focuses on alternate material concepts such as the development of new low coefficient of thermal expansion resins as matrices for graphite fiber reinforced composites, quartz fiber reinforced epoxies, and graphite reinforced glass. Low residual stress fabrication methods will be developed. When coupon specimens of these new material concepts have demonstrated the required surface accuracies and resistance to thermal distortion and microcracking, reflector panels will be fabricated and tested in simulated space environments. An important part of the program is the analytical modeling of environmental stability of these new composite materials concepts through constitutive equation development, modeling of microdamage in the composite matrix, and prediction of long term stability (including viscoelasticity). These analyses include both closed form and finite element solutions at the micro and macro levels.

Stein, Bland A.↗

NASA NTRS1992

Gravitational spreading of high terrain in Ishtar Terra, Venus

A plane-strain viscoelastic finite element algorithm derived by Melosh and Rafesky (1980, 1983) is used to model the gravitational spreading of the high topography of Ishtar Terra. Brittle failure and relaxation of topographic relief are predicted as functions of time for ranges of crustal thickness and thermal gradient and for observed ranges in topographic elevation and slope. For a crustal thickness greater than 10 km, a thermal gradient of 15 K/km or more, and average values of relief (3 km) and slope (3 deg), the topography relaxes to 25 percent of its original height in less than 10 m.y. Although initial failure, in the form of normal faulting on the highlands and margin and occasional shallow thrusting in the lowlands, is predicted to occur much earlier, values of horizontal surface strain that are likely to be observable (about 1 percent) do not accumulate until significant relaxation of relief begins.

Smrekar, Suzanne E.↗

NASA NTRS1981

Viscoelastic propellant effects on Space Shuttle Dynamics

The program of solid propellant research performed in support of the space shuttle dynamics modeling effort is described. Stiffness, damping, and compressibility of the propellant and the effects of many variables on these properties are discussed. The relationship between the propellant and solid rocket booster dynamics during liftoff and boost flight conditions and the effects of booster vibration and propellant stiffness on free free solid rocket booster modes are described. Coupled modes of the shuttle system and the effect of propellant stiffness on the interfaces of the booster and the external tank are described. A finite shell model of the solid rocket booster was developed.

Bugg, F.↗

NASA NTRS2000

Crustal Deformation Associated with Glacial Fluctuations in the Eastern Chugach Mountains, Alaska

The changes of the solid Earth in south central Alaska in response to two major glacial fluctuations on different temporal and spatial scales have been estimated and we evaluated their influence on the stress state and ongoing tectonic deformation of the region. During the recent (1993-1995) Bering Glacier surge, a large transfer of ice from the Bagley Ice Field to the Bering Glacier terminus region occurred. We estimated the elastic displacement of the solid Earth due to ice mass redistribution from Global Positioning System (GPS) measurements at sites near the surging glacier. We can account for these displacements by transfer of an ice volume of about 14 cubic km from the surge reservoir area to the terminus region. We examined the background seismicity (M(sub L) > 2.5) before, during, and after the surge. We found that the occurrence of small earthquakes (M(sub L) < 4.0) in the surge reservoir region increased during the surge time interval possibly in response to a decrease in ice mass. This suggests that a small decrease in the vertical stress, o,3, could be enough to modulate the occurrence of small, shallow earthquakes in this dominantly thrust fault setting. During this century the southern Alaska coastal glaciers have been undergoing an overall decrease in volume. Based on our compilation of changes in the extent and thickness of the coastal glaciers between the Malaspina and Bering, we calculated surface displacements due to the Earth's viscoelastic response to annual thinning and to the cumulative retreat over the last 100 years. The uplift of the region due to an average annual thinning rate of 1-6 m/yr in the ablation region is 1-12 mm/yr. For our reference model with a viscosity of 5 x 10(exp 19) Pa s for depths between approximately equal 40 and 200 km the total viscoelastic response due to the retreat over the last century may be as much as a couple of meters within the coastal ablation zone near Icy Bay. The maximum decrease in sigma(sub V) between 0 and 10 km was approximately equal 1.0 MPa, which is significant in relation to the stress drops in recent earthquakes (approximately 2 to 10 MPa) but small in relation to the estimated tectonic stress magnitude. Therefore the occurrence of an earthquake such as the St. Elias (1979, M(sub S) = 7.2) may have been advanced in time; however, most of the ongoing stress accumulation would be primarily due to tectonic forces.

Sauber, Jeanne↗

NASA NTRS1974

The response of many-storied structures under seismic actions

A seismic action which is characterized by the acceleration function of the seismic motion is analyzed. The equations of motion in the direction of the seismic action are determined for a completely rigid ground and for a ground with linear deformability characteristics (rotation deformability). The equation for the rotation motion is derived by considering elastic and viscoelastic properties of the foundation ground. The differential equations are written by treating the mechanical model of a many-storied structure as a fixed vertical beam with a uniformly distributed mass as well as with concentrated masses on each floor.

Munteanu, K.↗

NASA NTRS1979

Numerical and laboratory simulation of fault motion and earthquake occurrence

This paper reviews the simulation of earthquake occurrence by numerical and laboratory mechanical block models. Simple linear rheological elements are used with elastic forces driving the main events and viscoelastic forces being important for aftershock and creep occurrence. Friction and its dependence on velocity, stress, and displacement also play a key role in determining how, when, and where fault motion occurs. The discussion of the qualitative behavior of the simulators focuses on the manner in which energy is stored in the system and released by the unstable and stable sliding processes. The numerical results emphasize the statistics of earthquake occurrence and the correlations among source parameters.

Cohen, S. C.↗