Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Elastoviscoplastic buckling behavior of simply supported columns

The elastoviscoplastic buckling behavior of simply supported columns is investigated. The analysis is based on an incremental formulation that can deal with nonlinear geometrical and constitutive relations. Both viscoelastic and viscoplastic models are discussed, and Bodner-Partom's constitutive equations are used in the viscoplastic case. A numerical example is given to demonstrate the buckling behavior of columns under various temperature and loading levels.

Song, Y.↗
NASA NTRS2005

Two Dimensional Viscoelastic Stress Analysis of a Prototypical JIMO Turbine Wheel

The designers of the Jupiter Icy Moons Orbiter (JIMO) are investigating the potential of nuclear powered-electric propulsion technology to provide deep space propulsion. In one design scenario a closed-Brayton-cycle power converter is used to convert thermal energy from a nuclear reactor to electrical power for the spacecraft utilizing an inert gas as the working fluid to run a turboalternator as described in L.S. Mason, "A Power Conversion for the Jupiter Icy Moons Orbiter," Journal of Propulsion and Power, vol. 20, no. 5, pp. 902-910. A key component in the turboalternator is the radial flow turbine wheel which may be fabricated from a cast superalloy. This turbine wheel is envisioned to run continuously over the life of the mission, which is anticipated to be about ten years. This scenario places unusual material requirements on the turbine wheel. Unlike the case of terrestrial turbine engines, fatigue, associated with start-up and shut-down of the engine, foreign-object damage, and corrosion issues are insignificant and thus creep issues become dominate. The purpose of this paper is to present estimates for creep growth of a prototypical JIMO turbine wheel over a ten year life. Since an actual design and bill of materials does not exist, the results presented in this paper are based on preliminary concepts which are likely to evolve over time. For this reason, as well as computational efficiency, a simplified 2-D, in lieu of a 3-D, viscoelastic, finite element model of a prototypical turbine wheel will be utilized employing material properties for the cast superalloy MAR-M247. The creep data employed in this analysis are based on preliminary data being generated at NASA Glenn Research Center.

Gayda, John↗
NASA NTRS1999

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↗
NASA NTRS1985

Viscoelastic behavior and lifetime (durability) predictions

A method for lifetime or durability predictions for laminated fiber reinforced plastics is given. The procedure is similar to but not the same as the well known time-temperature-superposition principle for polymers. The method is better described as an analytical adaptation of time-stress-super-position methods. The analytical constitutive modeling is based upon a nonlinear viscoelastic constitutive model developed by Schapery. Time dependent failure models are discussed and are related to the constitutive models. Finally, results of an incremental lamination analysis using the constitutive and failure model are compared to experimental results. Favorable results between theory and predictions are presented using data from creep tests of about two months duration.

Brinson, R. F.↗
NASA NTRS1974

A model of convergent plate margins based on the recent tectonics of Shikoku, Japan

A viscoelastic finite element plate tectonic model is applied to displacement data for the island of Shikoku, Japan. The flow properties and geometry of the upper portions of the earth are assumed known from geophysical evidence, and the loading characteristics are determined from the model. The nature of the forces acting on the Philippine Sea plate, particularly in the vicinity of the Nankai trough, is determined. Seismic displacement data related to the 1946 Nankaido earthquake are modeled in terms of a thick elastic plate overlying a fluidlike substratum. The sequence of preseismic and seismic displacements can be explained in terms of two independent processes operating on elastic lithospheric plates: a strain accumulation process caused by vertical downward forces acting on or within the lithosphere in the vicinity of the trench, and a strain release process caused by plate failure along a preexisting zone on weakness. This is a restatement of Reid's elastic rebound theory in terms of elastic lithospheric plates.

Bischke, R. E.↗
NASA NTRS1988

Radiation and temperature effects on the time-dependent response of T300/934 graphite/epoxy

A time-dependent characterization study was performed on T300/934 graphite/epoxy in a simulated space environment. Creep tests on irradiated and nonirradiated graphite/epoxy and bulk resin specimens were carried out at temperatures of 72 and 250 F. Irradiated specimens were exposed to dosages of penetrating electron radiation equal to 30 years exposure at GEO-synchronous orbit. Radiation was shown to have little effect on the creep response of both the composite and bulk resin specimens at 72 F while radiation had a significant effect at 250 F. A healing process was shown to be present in the irradiated specimens where broken bonds in the epoxy due to radiation recombined over time to form cross-links in the 934 resin structure. An analytical micromechanical model was also developed to predict the viscoelastic response of fiber reinforced composite materials. The model was shown to correlate well with experimental results for linearly viscoelastic materials with relatively small creep strains.

Yancey, Robert N.↗
NASA NTRS1982

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.↗
NASA NTRS1988

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.↗
NASA NTRS2003

Algorithms for the Fractional Calculus: A Selection of Numerical Methods

Many recently developed models in areas like viscoelasticity, electrochemistry, diffusion processes, etc. are formulated in terms of derivatives (and integrals) of fractional (non-integer) order. In this paper we present a collection of numerical algorithms for the solution of the various problems arising in this context. We believe that this will give the engineer the necessary tools required to work with fractional models in an efficient way.

Diethelm, K.↗
NASA NTRS1983

NMR (Nuclear Magnetic Resonance) and macromolecular migration in a melt or in concentrated solutions

The purpose of this paper is to analyze the migration process of long polymer molecules in a melt or in concentrated solutions as it may be observed from the dynamics of the transverse magnetization of nuclear spins linked to these chains. The low frequency viscoelastic relaxation of polymer systems is known to be mainly controlled by the mechanism of dissociation of topological constraints excited on chains and which are called entanglements. This mechanism exhibits a strong dependence upon the chain molecular weight. These topological constraints also govern the diffusion process of polymer chains. So, the accurate description of the diffusion motion of a chain may be a convenient way to characterize disentanglement processes necessarily involved in any model proposed to explain viscoelastic effects.

Addad, J. P. C.↗
NASA NTRS1985

Constraints on the rheological structure of the mantle

Rheological models of the mantle are at present limited to radial symmetry, usually with homogeneous linearly viscous or viscoelastic incompressible layers. While such models are probably overly simple, they readily allow calculation of geophysical effects, such as post glacial rebound and related changes in the Earth's shape. They are also directly applicable to problems of global mantle flow and plate motions.

Oconnell, R. J.↗
NASA NTRSDate not supplied

Finite Element Thermal Model for Ultrasonic Welding of Thermoplastic Composites

Ultrasonic welding, UW, is a fast and energy-efficient technique for joining thermoplastic composites. It involves the use of high-frequency mechanical vibrations and a static welding force to melt and join adherends. Ultrasonic welding is an enabling technology to reduce the cost and complexity of in-space construction because lightweight thermoplastic composite components can be packaged compactly for launch and then efficiently assembled using supervised autonomous robotic technologies on site. However, the temperatures in space present challenges to UW, and it is critical that the efficacy of process parameters selected for manufacturing in space is understood prior to launch. To this end, a three-dimensional finite element model is presented in this technical presentation. The model incorporates equations for effects of viscoelastic heating and heat transfer on the welding process. The proposed model is applied to predict the temperature distribution in single lap shear, SLS, samples composed of AS4/PEEK (TC1200) composite that were welded using a terrestrial machine as part of a comprehensive weldability study. Thermocouple and infrared spot sensor data from the SLS samples provide empirical temperature measurements for calibration and validation of the UW thermal model. Calibration and validation of the model is an important step given the significant uncertainties in material properties such as the loss modulus and necessary assumptions in the physics implementations which allow the model to converge in an acceptable amount of time. The validated thermal model can be used to simulate the process for the space environment. Future validation sample testing is planned in a vacuum chamber. The result will be a model capable of guiding process parameter selection to ensure acceptable weld bonds when manufacturing in space.

Josh Fody↗
NASA NTRS1992

Dynamics and stability of a tethered centrifuge in low earth orbit

The three-dimensional attitude dynamics of a spaceborne tethered centrifuge for artificial gravity experiments in low earth orbit is analyzed using two different methods. First, the tethered centrifuge is modeled as a dumbbell with a straight viscoelastic tether, point tip-masses, and sophisticated environmental models such as nonspherical gravity, thermal perturbations, and a dynamic atmospheric model. The motion of the centrifuge during spin-up, de-spin, and steady-rotation is then simulated. Second, a continuum model of the tether is developed for analyzing the stability of lateral tether oscillations. Results indicate that the maximum fluctuation about the 1-g radial acceleration level is less than 0.001 g; the time required for spin-up and de-spin is less than one orbit; and lateral oscillations are stable for any practical values of the system parameters.

Quadrelli, B. M.↗
NASA NTRS1987

Interseismic deformation for normal fault earthquakes

Comparison between the coseismic vertical deformation associated with normal fault earthquakes and permanent deformation indicated by geologic structure indicates that large deformation must occur during the period between earthquakes. When sufficient geodetic and geologic information is available it is possible to estimate the spatial character of this interseismic deformation. A case in point is the 1983 Borah Peak, Idaho earthquake. This normal fault earthquake produced roughly 5 times as much basin subsidence as it did uplift of the adjacent mountain ranges. In contrast, geophysical and geological observations show that the basin is roughly as deep as the bounding range is high. A similar normal fault event in the Basin and Range, the 1959 Hebgen Lake earthquake, was also accompanied by substantially larger subsidence than uplift. In this case, postseismic geodetic measurements show broad regional uplift with a spatial pattern which is roughly consistent with that proposed for the Borah Peak event. Modeling suggests that postseismic viscoelastic relaxation and strain accumulation in an elastic lithosphere overlying a viscoelastic asthenosphere are possible physical mechanisms to generate interseismic uplift. These mechanisms may be the contemporary expression of those processes responsible for the high elevation of the Basin and Range Province.

Reilinger, Robert↗
NASA NTRS2004

Crustal Deformation and the Seismic Cycle across the Kodiak Islands, Alaska

The Kodiak Islands are located approximately 130 to 250 km from the Alaska-Aleutian Trench where the Pacific plate is underthrusting the North American plate at a rate of about 57 mm/yr. The southern extent of the 1964 Prince William Sound (${M-w}$ = 9.2) earthquake rupture occurred offshore and beneath the eastern portion of the Kodiak Islands. Here we report GPS results (1993-2001) from northern Kodiak Island that span the transition between the 1964 uplift region along the eastern coast and the region of coseismic subsidence further inland. The horizontal velocity vectors range from 22.9 $\pm$ 2.2 mm/yr at N26.3$\deg$W $\pm$ 2.5$\deg$, about 150 km from the trench, to 5.9 $\pm$ 1.3 mm/yr at N65.9$\deg$W $\pm$ 6.6$\deg$, about 190 km from the trench. Near the northeastern coast of Kodiak the velocity vector above the shallow, locked main thrust zone is between the orientation of PCFC-NOAM plate motion (N22$/deg$W) and the trench-normal (N3O$\deg$W). Further west, our geodetic results suggest the accumulation of shear strain that will be released eventually as left-lateral motion on upper plate faults such as the Kodiak Island fault. These results are consistent with the hypothesis that the difference between the Pacific-North American plate motion and the orientation of the down going slab would lead to 4-8 mm/yr of left-lateral slip. Short-term geodetic uplift rates range from 2 - 14 mm/yr, with the maximum uplift located near the axis of maximum subsidence during the 1964 earthquake. We evaluated alternate interseismic models for Kodiak to test the importance of various mechanisms responsible for crustal deformation rates. These models are based on the plate interface slip history inferred from earlier modeling of coseismic and post-seismic geodetic results. The horizontal (trench perpendicular) and vertical deformation rates across Kodiak are consistent with a model that includes the viscoelastic response to : (1) a downgoing Pacific plate interface that is locked at shallow depths,(2) coseismic slip in the 1964 and (3) interseismic creep below the seismogenic zone. The change in orientation of the horizontal velocity vector occurs down-dip from the locked main thrust zone. In southern Kodiak, the coseismic slip in the 1964 earthquake was smaller than in the northern Kodiak region; yet, the horizontal, interseismic velocities as a function of distance from the trench are comparable to those in northern Kodiak. Based on the earthquake history prior to, and following the 1964 earthquake, we hypothesize that the plate interface in southern Kodiak slips in more frequent large earthquakes than in northern Kodiak.

Sauber, Jeanne↗
NASA NTRS2005

Crustal Deformation and the Seismic Cycle Across the Kodiak Islands, Alaska

The Kodiak Islands are located approximately 120 to 250 km from the Alaska-Aleutian Trench - and are within the southern extent of the 1964 Prince William Sound (M(sub w) = 9.2) earthquake rupture zone. Here we report new campaign GPS results (1993-2001) from northern Kodiak Island. The rate and orientation of the horizontal velocities, relative to a fixed North America, range from 25.3 plus or minus 1.4 mm/yr at N32.9 deg. W plus or minus 2.5 to 8.5 plus or minus 1.0 mm/yr at N59.7 deg. W plus or minus 6.5 deg. In addition to the northern Kodiak data, we analyzed data from three southern Kodiak Island stations. The inland stations from both the northern and southern networks indicate a counterclockwise rotation of the velocity vectors. These results are consistent with the hypothesis that the difference between the Pacific-North American plate motion and the orientation of the down going slab would lead to 4-8 mm/yr of left-lateral slip above the unlocked, down-dip portion of the main thrust zone. The northern and southern Kodiak geodetic data are consistent with a model that includes the viscoelastic response to (1) a downgoing Pacific plate interface that is locked at shallow depths, (2) local coseismic slip in the 1964 earthquake, and (3) interseismic creep down dip from the seismogenic zone. Based on the pre-1964 and post-1944 earthquake history, as well as the pattern of interseismic earthquakes across the plate boundary zone, we hypothesize that in southern Kodiak some strain is released in moderate to large earthquakes between the occurrences of great earthquakes like the 1964 event.

Sauber, Jeanne↗
NASA NTRS2005

Satellite Geodetic Constraints On Earthquake Processes: Implications of the 1999 Turkish Earthquakes for Fault Mechanics and Seismic Hazards on the San Andreas Fault

Our principal activities during the initial phase of this project include: 1) Continued monitoring of postseismic deformation for the 1999 Izmit and Duzce, Turkey earthquakes from repeated GPS survey measurements and expansion of the Marmara Continuous GPS Network (MAGNET), 2) Establishing three North Anatolian fault crossing profiles (10 sitedprofile) at locations that experienced major surface-fault earthquakes at different times in the past to examine strain accumulation as a function of time in the earthquake cycle (2004), 3) Repeat observations of selected sites in the fault-crossing profiles (2005), 4) Repeat surveys of the Marmara GPS network to continue to monitor postseismic deformation, 5) Refining block models for the Marmara Sea seismic gap area to better understand earthquake hazards in the Greater Istanbul area, 6) Continuing development of models for afterslip and distributed viscoelastic deformation for the earthquake cycle. We are keeping close contact with MIT colleagues (Brad Hager, and Eric Hetland) who are developing models for S. California and for the earthquake cycle in general (Hetland, 2006). In addition, our Turkish partners at the Marmara Research Center have undertaken repeat, micro-gravity measurements at the MAGNET sites and have provided us estimates of gravity change during the period 2003 - 2005.

Reilinger, Robert↗
NASA NTRS1986

Viscoelastic deformation near active plate boundaries

Model deformations near the active plate boundaries of Western North America using space-based geodetic measurements as constraints are discussed. The first six months of this project were spent gaining familarity with space-based measurements, accessing the Crustal Dynamics Data Information Computer, and building time independent deformation models. The initial goal was to see how well the simplest elastic models can reproduce very long base interferometry (VLBI) baseline data. From the Crustal Dynamics Data Information Service, a total of 18 VLBI baselines are available which have been surveyed on four or more occasions. These data were fed into weighted and unweighted inversions to obtain baseline closure rates. Four of the better quality lines are illustrated. The deformation model assumes that the observed baseline rates result from a combination of rigid plate tectonic motions plus a component resulting from elastic strain build up due to a failure of the plate boundary to slip at the full plate tectonic rate. The elastic deformation resulting from the locked plate boundary is meant to portray interseismic strain accumulation. During and shortly after a large interplate earthquake, these strains are largely released, and points near the fault which were previously retarded suddenly catch up to the positions predicted by rigid plate models. Researchers judge the quality of fit by the sum squares of weighted residuals, termed total variance. The observed baseline closures have a total variance of 99 (cm/y)squared. When the RM2 velocities are assumed to model the data, the total variance increases to 154 (cm/y)squared.

Ward, S. N.↗