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 91 records · Page 5

Designing for time-dependent material response in spacecraft structures

To study the influence on overall deformations of the time-dependent constitutive properties of fiber-reinforced polymeric matrix composite materials being considered for use in orbiting precision segmented reflectors, simple sandwich beam models are developed. The beam models include layers representing the face sheets, the core, and the adhesive bonding of the face sheets to the core. A three-layer model lumps the adhesive layers with the face sheets or core, while a five-layer model considers the adhesive layers explicitly. The deformation response of the three-layer and five-layer sandwich beam models to a midspan point load is studied. This elementary loading leads to a simple analysis, and it is easy to create this loading in the laboratory. Using the correspondence principle of viscoelasticity, the models representing the elastic behavior of the two beams are transformed into time-dependent models. Representative cases of time-dependent material behavior for the facesheet material, the core material, and the adhesive are used to evaluate the influence of these constituents being time-dependent on the deformations of the beam. As an example of the results presented, if it assumed that, as a worst case, the polymer-dominated shear properties of the core behave as a Maxwell fluid such that under constant shear stress the shear strain increases by a factor of 10 in 20 years, then it is shown that the beam deflection increases by a factor of 1.4 during that time. In addition to quantitative conclusions, several assumptions are discussed which simplify the analyses for use with more complicated material models. Finally, it is shown that the simpler three-layer model suffices in many situations.

Hyer, M. W.↗

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

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

crustal mass changes↗

Postseismic surface deformations due to lithospheric and asthenospheric viscoelasticity

This paper proposes a model for postseismic surface deformations by attributing them to lithospheric and asthenospheric viscoelasticity. The model predicts that the deformations due to lithospheric viscoelasticity depend on the decrease in the effective shear modulus acting long after the lithospheric relaxation compared to that acting immediately following the earthquake. While such deformations are generally smaller than those associated with asthenospheric viscoelasticity, they occur on a shorter time scale and may be in opposite direction to both the motion occurring at the time of the earthquake and that occurring as the asthenospheric relaxation occurs.

Cohen, S. C.↗

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 deformation affects 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 for the Tharsis volcanoes on Mars. By means of a finite element code, we calculate stresses and displacements due to a volcano-shaped load emplaced on an elastic plate overlying a viscoelastic mantle. 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. As a result of flexure there are three regions where stresses become sufficiently large to cause failure by faulting, according to the Mohr-Coulomb criterion: at the surface of the plate just outward of the volcano, near the base of the elastic lithosphere beneath the center of the volcano, and on the upper flanks of the volcano early in its growth history.

Mcgovern, Patrick J.↗

Requirements for energy based constitutive modeling in tire mechanics

The history, requirements, and theoretical basis of a new energy based constitutive model for (rubber) material elasticity, hysteresis, and failure are presented. Energy based elasticity is handled by many constitutive models, both in one dimension and in three dimensions. Conversion of mechanical energy to heat can be modeled with viscoelasticity or as structural hysteresis. We are seeking unification of elasticity, hysteresis, and failure mechanisms such as fatigue and wear. An energy state characterization for failure criteria of (rubber) materials may provide this unification and also help explain the interaction of temperature effects with failure mechanisms which are described as creation of growth of internal crack surface. Improved structural modeling of tires with FEM should result from such a unified constitutive theory. The theory will also guide experimental work and should enable better interpretation of the results of computational stress analyses.

Luchini, John R.↗

Hydrodynamic stability of viscoelastic fluids - Importance of fluid model, overstability, and form of disturbance.

A linearized stability analysis has been applied to a fluid flowing in a gravity field between horizontal planes in Couette flow under conditions such that the temperature of the bottom plane exceeds that of the top. It is shown that, under conditions likely to be encountered with polymer solutions, oscillatory instabilities will not be controlling. Criteria are offered for ascertaining when an analysis based upon a second-order fluid model may be expected to yield physically meaningful results. It is also shown that, for the fluid model considered, critical conditions for stability are not changed when disturbances which vary in the flow direction are substituted for those which are a function of the coordinate transverse to the flow.

Mcintire, L. V.↗

Computer simulation of earthquakes

Two computer simulation models of earthquakes were studied for the dependence of the pattern of events on the model assumptions and input parameters. Both models represent the seismically active region by mechanical blocks which are connected to one another and to a driving plate. The blocks slide on a friction surface. In the first model elastic forces were employed and time independent friction to simulate main shock events. The size, length, and time and place of event occurrence were influenced strongly by the magnitude and degree of homogeniety in the elastic and friction parameters of the fault region. Periodically reoccurring similar events were frequently observed in simulations with near homogeneous parameters along the fault, whereas, seismic gaps were a common feature of simulations employing large variations in the fault parameters. The second model incorporated viscoelastic forces and time-dependent friction to account for aftershock sequences. The periods between aftershock events increased with time and the aftershock region was confined to that which moved in the main event.

Cohen, S. C.↗

Nonlinear viscoelastic characterization of polycarbonate

Uniaxial tensile creep and recovery data from polycarbonate at six temperatures and six stress levels are analyzed for nonlinear viscoelastic constitutive modeling. A theory to account for combined effects of two or more accelerating factors is presented.

Caplan, E. 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.↗

Viscoplastic analysis of an experimental cylindrical thrust chamber liner

A viscoplastic stress-strain analysis of an experimental cylindrical thrust chamber is presented. A viscoelastic constitutive model incorporating a single internal state variable that represents kinematic hardening was employed to investigate whether such a viscoplastic model could predict the experimentally observed behavior of the thrust chamber. Two types of loading cycles were considered: a short cycle of 3.5 sec. duration that corresponded to the experiments, and an extended loading cycle of 485.1 sec. duration that is typical of the Space Shuttle Main Engine (SSME) operating cycle. The analysis qualitatively replicated the deformation behavior of the component as observed in experiments designed to simulate SSME operating conditions. The analysis also showed that the mode and location in the component may depend on the loading cycle. The results indicate that using viscoplastic models for structural analysis can lead to a more realistic life assessment of thrust chambers.

Arya, Vinod K.↗

Lower-mantle viscosity constrained by seismicity around deglaciated regions

It is shown here that seismicity around the margins of deglaciated areas provides a constraint on the viscosity of the lower mantle. Calculations using a spherical, viscoelastic earth model show that the present-day magnitude of the stress fields induced in the lithosphere beneath the Laurentide and Fennoscandian ice sheets is very sensitive to the value of the lower-mantle viscosity. Stress of about 100 bar, sufficient to cause seismicity, can still remain in the lithosphere for lower-mantle viscosities greater than about 10 to the 22nd Pa-s; for lower-mantle viscosities of about 10 to the 21st Pa-s, only a few tens of bars of stress persist in the lithosphere today. This influence of lower-mantle viscosity on the state of stress in the lithosphere also has implications for the migration of stress from earthquakes, and hence for earthquake recurrence times.

Spada, G.↗

Viscoplastic analysis of an experimental cylindrical thrust chamber liner

A viscoplastic stress-strain analysis of an experimental cylindrical thrust chamber is presented. A viscoelastic constitutive model incorporating a single internal state variable that represents kinematic hardening was employed to investigate whether such a viscoplastic model could predict the experimentally observed behavior of the thrust chamber. Two types of loading cycles were considered: a short cycle of 3.5-s duration that corresponded to the experiments, and an extended loading cycle of 485.1 s duration that is typical of the Space Shuttle Main Engine (SSME) operating cycle. The analysis qualitatively replicated the deformation behavior of the component as observed in experiments designed to simulate SSME operating conditions. The analysis also showed that the mode and location of failure in the component may depend on the loading cycle. The results indicate that using viscoplastic models for structural analysis can lead to a more realistic life assessment of thrust chambers.

Arya, Vinod K.↗

Present-day secular variations in the zonal harmonics of earth's geopotential

The mathematical formulation required for predicting secular variation in the geopotential is developed for the case of a spherically symmetric, self-gravitating, viscoelastic earth model and an arbitrary surface load which can include a gravitational self-consistent ocean loading component. The theory is specifically applied to predict the present-day secular variation in the zonal harmonics of the geopotenial arising from the surface mass loading associated with the late Pleistocene glacial cycles. A procedure is outlined in which predictions of the present-day geopotential signal due to the late Pleistocene glacial cycles may be used to derive bounds on the net present-day mass flux from the Antarctic and Greenland ice sheets to the local oceans.

Mitrovica, J. X.↗

A spectral formalism for computing three-dimensional deformations due to surface loads. 2: Present-day glacial isostatic adjustment

Using a spherically symmetric, self-gravitating, linear viscoelastic Earth model, we predict present-day three-dimensional surface deformation rates and baseline evolutions arising as a consequence of the late Pleistocene glacial cycles. In general, we use realistic models for the space-time geometry of the final late Pleistocene deglaciation event and incorporate a gravitationally self-consistent ocean meltwater redistribution. The predictions of horizontal velocity presented differ significantly, in both their amplitude and their spatial variation, from those presented in earlier analysis of others which adopted simplified models of both the late Pleistocene ice history and the Earth rheology. An important characteristic of our predicted velocity fields is that the melting of the Laurentide ice sheet over Canada is capable of contributing appreciably to the adjustment in Europe. The sensitivity of the predictions to variations in mantle rheology is investigated by considering a number of different Earth models, and by computing appropriate Frechet kernels. These calculations suggest that the sensitivity of the deformations to the Earth's rheology is significant and strongly dependent on the location of the site relative to the ancient ice sheet. The effects on the predictions of three-dimensional deformation rates of altering the ice history or adopting approximate models for the ocean meltwater redistribution have also been considered and found to be important (the former especially so). Finally, for a suite of Earth models we provide predictions of the velocity of a number of baselines in North America and Europe. We find that, in general, both radial and tangential motions contribute significantly to baseline length changes, and that these contributions are a strong function of the Earth model. We have, furthermore, found a set of Earth models which, together with the ICE-3G deglaciation chronology, produce predictions of baseline length changes that are consistent with very long baseline interferometry measurements of baselines within Europe.

Mitrovica, J. X.↗

A singularity free approach to post glacial rebound calculations

Calculating the post glacial response of a viscoelastic Earth model using the exponential decay normal mode technique leads to intrinsic singularities if viscosity varies continuously as a function of radius. We develop a numerical implementation of the Complex Real Fourier transform (CRFT) method as an accurate and stable procedure to avoid these singularities. Using CRFT, we investigate the response of a set of Maxwell Earth models to surface loading. We find that the effect of expanding a layered viscosity structure into a continuously varying structure is to destroy the modes associated with the boundary between layers. Horizontal motion is more sensitive than vertical motion to the viscosity structure just below the lithosphere. Horizontal motion is less sensitive to the viscosity of the lower mantle than the vertical motion is. When the viscosity increases at 670 km depth by a factor of about 60, the response of the lower mantle is close to its elastic limit. Any further increase of the viscosity contrast at 670 km depth or further increase of viscosity as a continuous function of depth starting from 670 km depth is unlikely to be resolved.

Fang, Ming↗

Viscoelastic Postseismic Rebound to Strike-Slip Earthquakes in Regions of Oblique Plate Convergence

According to the slip partitioning concept, the trench parallel component of relative plate motion in regions of oblique convergence is accommodated by strike-slip faulting in the overriding continental lithosphere. The pattern of postseismic surface deformation due to viscoelastic flow in the lower crust and asthenosphere following a major earthquake on such a fault is modified from that predicted from the conventual elastic layer over viscoelastic halfspace model by the presence of the subducting slab. The predicted effects, such as a partial suppression of the postseismic velocities by 1 cm/yr or more immediately following a moderate to great earthquake, are potentially detectable using contemporary geodetic techniques.

Cohen, Steven C.↗

Prediction of Long-Term Strength of Thermoplastic Composites Using Time-Temperature Superposition

Accelerated tests for composite failure were investigated. Constant ramp transverse strength tests on thermoplastic composite specimens were conducted at four temperatures from 300 F to 450 F and five duration times from 0.5 sec to 24 hrs. Up to 400 F, the time-temperature-superposition method produces a master curve allowing strength at longer times to be estimated from strength tests conducted over shorter times but at higher temperatures. The shift factors derived from compliance tests applied well to the strength data. To explain why strength behaved similar to compliance, a viscoelastic fracture model was investigated based on the hypothesis that the work of fracture for crack initiation at some critical flaw remains constant with time and temperature. The model, which used compliance as input, was found to fit the strength data only if the critical fracture energy was allowed to vary with stress rate. Fracture tests using double cantilever beam specimens were conducted from 300 F to 450 F over time scales similar to the strength study. The toughness data showed a significant change with loading rate, less variation with temperature, did not form a master curve, and could not be correlated with the fracture model. Since the fracture model did not fit the fracture data, an alternative explanation based on the dilatational strain energy density was proposed. However the usefulness of this model is severely limited because it relies on a critical parameter which varies with loading rate.

Reeder, James R.↗

Tensor Fields for Use in Fractional-Order Viscoelasticity

To be able to construct viscoelastic material models from fractional0order differentegral equations that are applicable for 3D finite-strain analysis requires definitions for fractional derivatives and integrals for symmetric tensor fields, like stress and strain. We define these fields in the body manifold. We then map them ito spatial fields expressed in terms of an Eulerian or Lagrangian reference frame where most analysts prefer to solve boundary problems.

Freed, Alan D.↗