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 541 records · Page 30

Analysis of crack-induced-craze in polymers

In this paper, the viscoelastic boundary element method is used to estimate the opening displacement and the envelope stress on the surface of an isolated crack-induced-craze system. To predict the propagation history of both the crack and the craze in a polymer sheet, the material properties of the glassy polymers are represented by a generalized linear viscoelastic model. Results are compared with the theoretical micromechanics predictions. Good agreements are obtained.

Sun, B. N.↗

Linear and nonlinear mechanical properties of a series of epoxy resins

The linear viscoelastic properties have been measured for a series of bisphenol-A-based epoxy resins cured with the diamine DDS. The linear viscoelastic master curves were constructed via time-temperature superposition of frequency dependent G-prime and G-double-prime isotherms. The G-double-prime master curves exhibited two sub-Tg transitions. Superposition of isotherms in the glass-to-rubber transition (i.e., alpha) and the beta transition at -60 C was achieved by simple horizontal shifts in the log frequency axis; however, in the region between alpha and beta, superposition could not be effected by simple horizontal shifts along the log frequency axis. The different temperature dependency of the alpha and beta relaxation mechanisms causes a complex response of G-double-prime in the so called alpha-prime region. A novel numerical procedure has been developed to extract the complete relaxation spectra and its temperature dependence from the G-prime and G-double-prime isothermal data in the alpha-prime region.

Curliss, D. B.↗

A nonlinear least-squares inverse analysis of strike-slip faulting with application to the San Andreas fault

A nonlinear weighted least-squares analysis was performed for a synthetic elastic layer over a viscoelastic half-space model of strike-slip faulting. Also, an inversion of strain rate data was attempted for the locked portions of the San Andreas fault in California. Based on an eigenvector analysis of synthetic data, it is found that the only parameter which can be resolved is the average shear modulus of the elastic layer and viscoelastic half-space. The other parameters were obtained by performing a suite of inversions for the fault. The inversions on data from the northern San Andreas resulted in predicted parameter ranges similar to those produced by inversions on data from the whole fault.

Williams, Charles A.↗

Compression creep of graphite/epoxy laminates monitored using moire interferometry

Uniaxial compressive creep loadings were applied for ten hours to various graphite/epoxy laminates with a hole. All tests were performed at room temperatures. In-plane displacements were monitored using moire interferometry. Moire fringe patterns were reduced to in-plane strain fields numerically. The fiber-dominated laminates tested exhibited little or no viscoelastic behavior during the ten hour period. Conversely, the matrix-dominated laminates tested did exhibit viscoelastic behavior during the ten hour test time.

Tuttle, M. E.↗

Finite-element models of continental extension

Numerical models of the initial deformation of extending continental lithosphere, computed to investigate the control of preexisting thermal and mechanical heterogeneities on the style of deformation, are presented. The finite element method is used to calculate deformation with a viscoelastic-plastic model for the lithosphere. Comparisons of the results of analytic models and finite-element models using this method show that good results may be obtained by the numerical technique, even with elements containing both brittle and viscoelastic sampling points. It is shown that the gross style of initial extensional deformation is controlled by the depth and width of the initial heterogeneity which localizes deformation.

Lynch, H. David↗

Analytical and experimental study of the vibration of bonded beams with a lap joint

A theoretical model to study the flexural vibration of a bonded lap joint system is described in this paper. First, equations of motion at the joint region are derived using a differential element approach. The transverse displacements of the upper and lower beam are considered to be different. The adhesive is assumed to be linearly viscoelastic and the widely used Kelvin-Voight model is used to represent the viscoelastic behavior of the adhesive. The shear force at the interface between the adhesive and the beam is obtained from the simple bending motion equations of the two beams. The resulting equations of motion are combined with the equations of transverse vibration of the beams in the unjointed regions. These are later solved as a boundary value problem to obtain the eigenvalues and eigenvectors of the system. The model can be used to predict the natural frequencies, modal damping ratios, and mode shapes of the system for free vibration. Good agreement between numerical and experimental results was obtained for a system of graphite epoxy beams lap-jointed by an epoxy adhesive.

Rao, M. D.↗

Crustal deformation along the San Andreas, California

The goal is to achieve a better understanding of the regional and local deformation and crustal straining processes in western North America, particularly the effects of the San Andreas and nearby faults on the spatial and temporal crustal deformation behavior. Construction of theoretical models based on the mechanics of coupled elastic plate, viscoelastic foundation and large scale crack mechanics provide a rational basis for the interpretation of seismic and aseismic anomalies and expedite efforts in forecasting the stability of plate boundary deformation. Special focus is placed on the three dimensional time dependent surface deformation due to localized slippage in a elastic layer coupled to a visco-elastic substrate. The numerical analysis is based on a 3-D boundary element technique. Extension to visco-elastic coupling demands the derivation of 3-D time dependent Green's function. This method was applied to analyze the viscoelastic surface displacements due to a dislocated embedded patch. Surface uplift as a function of time and position are obtained. Comparisons between surface uplift for long and short dislocated patches are made.

Li, Victor C.↗

Graphite/epoxy Composite Laminates with Co-cured Interlaminar Damping Layers

Damped composite laminates were fabricated by co-curing viscoelastic damping film with graphite/epoxy prepreg plies. The dynamic response of the damped plates was measured using an impulse response technique and compared with the response of similar undamped laminates. Modal damping was computed from the frequency response data. Micrographs of the damped laminates showed that the damping layers retained their integrity during the fabrication process. The layers significantly increased the damping in the composite laminates. The use of the constrained viscoelastic film as an integral part of composite structures appears to be a feasible approach to passive vibration control. Composite plates manufactured with co-cured damping layers may have commercial applications in cases where light weight, strength, and vibration and noise reduction are important considerations.

Pereira, J. Michael↗

Numerical simulation of the non-Newtonian mixing layer

This work is a continuing effort to advance our understanding of the effects of polymer additives on the structures of the mixing layer. In anticipation of full nonlinear simulations of the non-Newtonian mixing layer, we examined in a first stage the linear stability of the non-Newtonian mixing layer. The results of this study show that, for a fluid described by the Oldroyd-B model, viscoelasticity reduces the instability of the inviscid mixing layer in a special limit where the ratio (We/Re) is of order 1 where We is the Weissenberg number, a measure of the elasticity of the flow, and Re is the Reynolds number. In the present study, we pursue this project with numerical simulations of the non-Newtonian mixing layer. Our primary objective is to determine the effects of viscoelasticity on the roll-up structure. We also examine the origin of the numerical instabilities usually encountered in the simulations of non-Newtonian fluids.

Azaiez, Jalel↗

Aspects of modelling the tectonics of large volcanoes on the terrestrial planets

Analytic solutions for the responses of planetary lithospheres to volcanic loads have been used to model faulting and infer elastic plate thicknesses. Predictions of the distribution of faulting around volcanic loads, based on the application of Anderson's criteria for faulting to the results of the models, do not agree well with observations. Such models do not give the stress state in the load itself, but only suggest a state of horizontal compressive stress there. Further, these models have considered only the effect of an instantaneously emplaced load. They do not address the time evolution of stresses, nor do they consider the effect of a load which grows. A finite element approach allows us to assign elements to the load itself, and thus permits calculation of the stress state and stress history within the edifice. The effects of episodic load growth can also be treated. When these effects are included, models give much better agreement with observations. We use the finite element code TECTON to construct axisymmetric models of volcanoes resting on an elastic lithospheric plate overlying a viscoelastic asthenosphere. We have implemented time-dependent material properties in order to simulate incremental volcano growth. The viscoelastic layer was taken to extend to a sufficient depth so that a rigid lower boundary has no significant influence on the results. The code first calculates elastic deformations and stresses and then determines the time-dependent viscous deformations and stresses. Time in the model scales as the Maxwell time tau(m) in the asthenosphere. We consider a volcano 25 km in height and 200 km in radius on an elastic lithosphere 40 km thick (parameters approximately appropriate to Ascraeus Mons). The volcano consists of three load increments applied at intervals of 1000 tau(m). Contours of maximum deviatoric stress in the fully-grown edifice at the conclusion of flexure (t = 3000 tau(m)) are shown.

Mcgovern, Patrick J.↗

On the design of structural components using materials with time-dependent properties

The application of the elastic-viscoelastic correspondence principle is presented as a design tool for structural design engineers for composite material applications. The classical problem of cantilever beams is used as the illustration problem. Both closed-form and approximate numerical solutions are presented for several different problems. The application of the collocation method is presented as a viable and simple design tool to determine the time-dependent behavior and response of viscoelastic composite beams under load.

Rodriguez, Pedro I.↗

Viscosity estimates for the crust and upper mantle from patterns of lacustrine shoreline deformation in the Eastern Great Basin

The deformed shorelines of Lake Bonneville constitute a classic source of information on lithospheric elastic thickness and upper mantle viscosity. We describe and apply a new model to a recently augmented data set. New data better constrain both the complex spatio-temporal pattern of the lake load and the crustal deformation response to that load. The history of lake level fluctuations has been significantly refined and somewhat modified. This is due to both more radiocarbon dates from within the Bonneville basin and to an improved calibration of the radiocarbon timescale itself. The data which constrain the crustal deformation pattern consist of ages and shoreline elevations from several hundred points which sample three major levels of Lake Bonneville and corresponding elevations from the high stands of three smaller lakes situated to the west of Lake Bonneville. The geometry of the Earth model incorporates an arbitrary number of layers overlying a half-space, and the rheology of each level can accommodate an arbitrary number of Maxwell viscoelastic elements in parallel. The inverse modeling comprises three complementary approaches: for the simplest configurations, we performed a direct search of the parameter space and delineated the irregular boundary of the subspace of acceptable models. For more complex configurations, we constrained the elastic parameters to their seismically determined values and then solved for viscosity versus depth profiles by either expressing the log(viscosity) versus log(depth) profile as a series of specially constructed orhtogonal polynomials, or by allowing each of 8-10 layers (plus the half-space) to have an independently determined viscosity. We found that the data do not strongly support (nor can they conclusively exclude) a more complex rheology than simple Maxwell viscoelasticity. The orthogonal polynomial solution exhibits an essentially monotonic decrease in viscosity with depth.

Bills, Bruce G.↗

Reply to the comment by Y. Ricard, R. Sabadini, and G. Spada

We welcome this opportunity which Drs. Ricard, Sabadini and Spada (RSS) have offered us to elaborate on some points concerning the connection between convective circulations and rotational dynamics. The relationship of the linear viscoelastic normal-mode theory to other approaches, such as viscous and the non-linear viscoelastic treatments, can be found in Moser et al. (1993). We will address the issues raised by RSS by first going over the reference system. This will be followed by a discussion concerning long-term polar wander and relative angular momentum. We will also go over the issues concerning the off-diagonal terms in the total moment of inertia and polar wander.

Moser, Jiri↗

Transient Creep of a Composite Lower Crust: Constitutive Theory - 1

A composite model is proposed to describe the time-dependent response of the Earth's lower crust. The motivation for such it model is twofold: First, new observations of widespread postseismic deformation indicate that the deep continental crust responds viscoelastically, having both long-and short-term decay times. Second, by any number of observationally based rationales, the lower crust is compositionally and structurally heterogeneous over many length scales. For heterogeneities that have much smaller characteristic lengths than the minimum deformation wavelength of interest, the aggregate rheology can be described by composite media theory. For wavelengths of the order of the thickness of the lower crust (approx. = 25-40 km) and larger, composite theory may be applied to heterogeneities that are smaller than about several hundred meters, or equivalent to the vertical extent of a thick lower crustal mylonitic shear zone. The composite media theory developed here is constructed using both Eshelhy-Mori-Tanaka theory for aligned generalized spheroidal inclusions and a generalized self-consistent method. The inclusions and matrix are considered to be Maxwellian viscoelastic: a rheology that is consistent with past homogeneous models of postseismic stress relaxation. The composite theory presented here introduces a transient response to a suddenly imposed stress field which does not appear in homogeneous Maxwell models. Analytic expressions for the amplitude and duration of the transient and for the effective long-and short-term viscosities of the composite are given which describe the sensitivity to inclusion concentration (phi), to shape, and to ratio of inclusion-to-matrix viscosity (R).

Ivins, Erik R.↗

Hubble Space Telescope (HST) Solar Array Damper

This paper describes the design of a solar array damper that will be built into each of two new solar arrays to be installed on the Hubble Space Telescope (HST) during Servicing Mission 3. On this mission, currently scheduled for August, 2000, two "rigid" solar array wings will replace the "flexible" wings currently providing power for HST. Dynamic interaction of these wings with the telescope spacecraft can affect the Pointing Control System. The damper, which is integral to the mast of the solar array, suppresses the fundamental bending modes of the deployed wings at 1.2 Hz (in-plane) and 1.6 Hz (out-of-plane). With the flight version of the damper, modal damping of 2.3% of critical is expected over the temperature range of -4 C to 23 C with a peak damping level of 3.9%. The unique damper design is a combination of a titanium spring and viscoelastic-shear-lap dashpot. The damper was designed using a system finite element model of the solar array wing and measured viscoelastic material properties. Direct complex stiffness (DCS) testing was performed to characterize the frequency- and temperature-dependent behavior of the damping prior to fixed-base modal testing of the wing at NASA/Goddard Space Flight Center (NASA/GSFC).

Maly, J. R.↗

Special Considerations in Selection of Fabric Film Laminates for Use in Inflatable Structures

Inflatable structures are gaining wide support in planetary scientific missions as well as commercial applications. For such applications a new class of fabric/film laminates is being considered for use as a structural gas envelope. The emerging composite materials are a result of recent advances in the manufacturing of lightweight, high strength fibers, fabrics and scrims. The lamination of these load-carrying members with the proper gas barriers results in a wide range of materials suitable for various loading and environmental conditions. Polyester-based woven fabrics laminated to thin homogenous film of polyester are an example of this class. This fabric/film laminate is being considered for the development of a material suitable for building large gas envelopes for use in the NASA ultra long duration balloon program (ULDB). Compared to commercial homogenous films, the material provides relatively high strength to weight ratio as well as better resistance to crack and tear propagation, The mechanical, creep and viscoelastic properties of these fabric film laminates have been studied to form a material model. Preliminary analysis indicates that the material is highly viscoelastic. The mechanical properties of this class of materials will be discussed in some details.

Said, Magdi A.↗

On Lateral Viscosity Contrast in the Mantle and the Rheology of Low-Frequency Geodynamics

Mantle-wide heterogeneity is largely controlled by deeply penetrating thermal convective currents. These thermal currents are likely to produce significant lateral variation in rheology, and this can profoundly influence overall material behaviour. How thermally related lateral viscosity variations impact models of glacio-isostatic and tidal deformation is largely unknown. An important step towards model improvement is to quantify, or bound, the actual viscosity variations that characterize the mantle. Simple scaling of viscosity to shear-wave velocity fluctuations yields map-views of long- wavelength viscosity variation. These give a general quantitative description and aid in estimating the depth dependence of rheological heterogeneity throughout the mantle. The upper mantle is probably characterized by two to four orders of magnitude variation (peak-to-peak). Discrepant time-scales for rebounding Holocene shorelines of Hudson Bay and southern Iceland are consistent with this characterization. Results are given in terms of a local average viscosity ratio, (Delta)eta(bar)(sub i), of volumetric concentration, phi(sub i). For the upper mantle deeper than 340 km the following reasonable limits are estimated for (delta)eta(bar) approx. equal 10(exp -2): 0.01 less than or equal to phi less than or equal to 0.15. A spectrum of ratios (Delta)eta(bar)(sub i) less than 0.1 at concentration level eta(sub i) approx. equal 10(exp -6) - 10(exp -1) in the lower mantle implies a spectrum of shorter time-scale deformational response modes for second-degree spherical harmonic deformations of the Earth. Although highly uncertain, this spectrum of spatial variation allows a purely Maxwellian viscoelastic rheology simultaneously to explain all solid tidal dispersion phenomena and long-term rebound-related mantle viscosity. Composite theory of multiphase viscoelastic media is used to demonstrate this effect.

Ivins, Erik R.↗

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↗