Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “creep deformation”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Characterization of Metal Matrix Composites

Experimental methods were developed, adapted, and applied to the characterization of a metal matrix composite system, namely, silicon carbide/aluminim (SCS-2/6061 Al), and its constituents. The silicon carbide fiber was characterized by determining its modulus, strength, and coefficient of thermal expansion. The aluminum matrix was characterized thermomechanically up to 399 C (750 F) at two strain rates. The unidirectional SiC/Al composite was characterized mechanically under longitudinal, transverse, and in-plane shear loading up to 399 C (750 F). Isothermal and non-isothermal creep behavior was also measured. The applicability of a proposed set of multifactor thermoviscoplastic nonlinear constitutive relations and a computer code was investigated. Agreement between predictions and experimental results was shown in a few cases. The elastoplastic thermomechanical behavior of the composite was also described by a number of new analytical models developed or adapted for the material system studied. These models include the rule of mixtures, composite cylinder model with various thermoelastoplastic analyses and a model based on average field theory. In most cases satisfactory agreement was demonstrated between analytical predictions and experimental results for the cases of stress-strain behavior and thermal deformation behavior at different temperatures. In addition, some models yielded detailed three-dimensional stress distributions in the constituents within the composite.

Daniel, I. M.↗
DOE OSTI2023

A constitutive framework for rocks undergoing solid dissolution

Here, we formulate a time-dependent damage theory for rocks subjected to mechanical deformation and solid dissolution. The constitutive description is inspired by the transition state theory, which states that the rate of dissolution is a function of the reactive surface area measured through the crack density in the volume. We use a gradient-enhanced damage framework in which damage depends on the deformation of the material as well as on the amount of solid mass dissolved over time. The gradient-enhanced formulation is characterized by a three-field variational formulation with the solid displacement, nonlocal equivalent strain, and nonlocal rate of solid dissolution as the basic state variables. Traditionally, time-independent damage theories have only allowed damage to increase with increasing external load. In the proposed framework, the degree of damage may increase due to solid dissolution even when the external load is held fixed. In this way, solid dissolution is viewed as a process that is responsible for bringing about rate-dependent effects such as creep and stress-relaxation, which are two common features of geomaterial behavior.

42 ENGINEERING↗
NASA NTRS1972

Fracture Behavior in Nylon 6 Fibers

Electron paramagnetic resonance (EPR) techniques are used to determine the number of free radicals produced during deformation leading to fracture of nylon 6 fibers. A reaction rate molecular model is proposed to explain some of the deformation and bond rupture behavior leading to fracture. High-strength polymer fibers are assumed to consist of a sandwich structure of disordered and ordered regions along the fiber axis. In the disordered or critical flaw regions, tie chains connecting the ordered or crystalline block regions are assumed to have a statistical distribution in length. These chains are, therefore, subjected to different stresses. The effective length distribution was determined by EPR. The probability of bond rupture was assumed to be controlled by reaction-rate theory with a stress-aided activation energy and behavior of various loadings determined by numerical techniques. The model is successfully correlated with experimental stress, strain, and bond rupture results for creep, constant rate loadings, cyclic stress, stress relaxation and step strain tests at room temperature.

Lloyd, B. A.↗
NASA NTRS2002

Three-Dimensional Stress Fields and Slip Systems for Single Crystal Superalloy Notched Specimens

Single crystal superalloys have become increasingly popular for turbine blade and vane applications due to their high strength, and creep and fatigue resistance at elevated temperatures. The crystallographic orientation of a single crystal material greatly affects its material properties, including elastic modulus, shear modulus, and ductility. These directional properties, along with the type of loading and temperature, dictate an anisotropic response in the yield strength, creep resistance, creep rupture ductility, fatigue resistance, etc. A significant amount of research has been conducted to determine the material properties in the <001> orientation, yet the material properties deviating from the <001> orientation have not been assessed for all cases. Based on the desired application and design criteria, a crystal orientation is selected to yield the maximum properties. Currently, single crystal manufacturing is able to control the primary crystallographic orientation within 15 of the target orientation, which is an acceptable deviation to meet both performance and cost guidelines; the secondary orientation is rarely specified. A common experiment is the standard load-controlled tensile test, in which specimens with different orientations can be loaded to observe the material response. The deformation behavior of single-crystal materials under tension and compression is known to be a function of not only material orientation, but also of varying microdeformation (i.e. dislocation) mechanisms. The underlying dislocation motion causes deformation via slip, and affects the activation of specific slip systems based on load and orientation. The slip can be analyzed by observing the visible traces left on the surface of the specimen from the slip activity within the single crystal material. The goal of this thesis was to predict the slip systems activated in three-dimensional stress fields of a notched tensile specimen, as a function of crystal orientation, using finite element analysis without addressing microstructural deformation mechanisms that govern their activation. Out of three orientations tested, the specimen with a [110] load orientation and a [001] growth direction had the lowest maximum resolved shear stress; this specimen orientation appears to be the best design candidate for a tensile application.

Magnan, Shannon M.↗
DOE OSTI2023

Review: Inelastic Constitutive Modeling: Polycrystalline Materials

This article provides a literature review that details the development of inelastic constitutive modeling as it relates to polycrystalline materials. This review distinguishes between inelastic constitutive models that account for nonlinear behavior at the microstructural level, time-independent classic plasticity models, and time-dependent unified models. Particular emphasis is placed on understanding the underlying theoretical framework for unified viscoplasticity models where creep and classical plasticity behavior are considered the result of applied boundary conditions instead of separable rates representing distinct physical mechanisms. This article establishes a clear understanding of the advantages of the unified approach to improve material modeling. This review also discusses recent topics in constitutive modeling that offer new techniques that bridge the gap between the microstructure and the continuum.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗
NASA NTRS1988

Life prediction modeling based on strainrange partitioning

Strainrange partitioning (SRP) is an integrated low-cycle-fatigue life predicting system. It was created specifically for calculating cyclic crack initiation life under severe high-temperature fatigue conditions. The key feature of the SRP system is its recognition of the interacting mechanisms of cyclic inelastic deformation that govern cyclic life at high temperatures. The SRP system bridges the gap between the mechanistic level of understanding that breeds new and better materials and the phenomenological level wherein workable engineering life prediction methods are in great demand. The system was recently expanded to address engineering fatigue problems in the low-strain, long-life, nominally elastic regime. This breakthrough, along with other advances in material behavior and testing technology, has permitted the system to also encompass low-strain thermomechanical loading conditions. Other important refinements of the originally proposed method include procedures for dealing with life-reducing effects of multiaxial loading, ratcheting, mean stresses, nonrepetitive (cumulative loading) loading, and environmental and long-time exposure. Procedure were also developed for partitioning creep and plastic strain and for estimating strainrange versus life relations from tensile and creep rupture properties. Each of the important engineering features of the SRP system are discussed and examples shown of how they help toward predicting high-temperature fatigue life under practical, although complex, loading conditions.

Halford, Gary R.↗
DOE OSTI2021

In situ observation of the deformation and fracture of an alumina-alumina ceramic-matrix composite at elevated temperature using x-ray computed tomography

The deformation and fracture of a 0/90° Nextel 720 mullite fibre / alumina matrix ceramic-matrix composite material has been observed by in situ synchrotron x-ray computed micro-tomography at ambient temperature and at 1100 °C. A three-point loading configuration was used and samples were loaded monotonically until final fracture. Here, the flexural strength was found to be unaffected by temperature. 3D visualisation of the microstructure showed fracture occurred by propagation of internal cracks at stresses close to the flexural strength. Shear fractures propagated across fibre bundles that were oriented perpendicular to the flexural stress at both room temperature and 1100 °C. At the higher temperature, an additional failure mode of delamination was observed, that may arise from relaxation of thermal residual stresses and creep.

36 MATERIALS SCIENCE↗
NASA NTRS2023

Results from InSight Robotic Arm Activities

The InSight lander carried an Instrument Deployment System (IDS) that included an Instrument Deployment Arm (IDA), scoop, five finger “claw” grapple, forearm-mounted Instrument Deployment Camera (IDC) requiring arm motion to image a target, and landermounted Instrument Context Camera (ICC), designed to image the workspace, and to place the instruments onto the surface. As originally proposed, the IDS included a previously built arm and flight spare black and white cameras and had no science objectives or requirements, or expectation to be used after instrument deployment (90 sols). During project development the detectors were upgraded to color, and it was recognized that the arm could be used to carry out a wide variety of activities that would enable both geology and physical properties investigations. During surface operations for two martian years, the IDA was used during major campaigns to image the surface around the lander, to deploy the instruments, to assist the mole in penetrating beneath the surface, to bury a portion of the seismometer tether, to clean dust from the solar arrays to increase power, and to conduct a surface geology investigation including soil mechanics and physical properties experiments. No other surface mission has engaged in such a sustained and varied campaign of arm and scoop activities directed at such a diverse suite of objectives. Images close to the surface and continuous meteorology measurements provided important constraints on the threshold friction wind speed needed to initiate aeolian saltation and surface creep. The IDA was used extensively for almost 22 months to assist the mole in penetrating into the subsurface. Soil was scraped into piles and dumped onto the seismometer tether six times in an attempt to bury the tether and ∼ 30% was entrained in the wind and dispersed downwind 1-2 m, darkening the surface. Seven solar array cleaning experiments were conducted by dumping scoops of soil from 35 cm above the lander deck during periods of high wind that dispersed the sand onto the panels that kicked dust off of the panels into suspension in the atmosphere, thereby increasing the power by ∼15% during this period. Final IDA activities included an indentation experiment that used the IDA scoop to push on the ground to measure the plastic deformation of the soil that complemented soil mechanics measurements from scoop interactions with the surface, and two experiments in which SEIS measured the tilt from the arm pressing on the ground to derive near surface elastic properties.

Mars↗
NASA NTRS1990

Elevated temperature crack growth in aluminum alloys: Tensile deformation of 2618 and FVS0812 aluminum alloys

Understanding the damage tolerance of aluminum alloys at elevated temperatures is essential for safe applications of advanced materials. The objective of this project is to investigate the time dependent subcritical cracking behavior of powder metallurgy FVS0812 and ingot metallurgy 2618 aluminum alloys at elevated temperatures. The fracture mechanics approach was applied. Sidegrooved compact tension specimens were tested at 175, 250, and 316 C under constant load. Subcritical crack growth occurred in each alloy at applied stress intensity levels (K) of between about 14 and 25 MPa/m, well below K (sub IC). Measured load, crack opening displacement and displacement rate, and crack length and growth rate (da/dt) were analyzed with several continuum fracture parameters including, the C-integral, C (sub t), and K. Elevated temperature growth rate data suggest that K is a controlling parameter during time dependent cracking. For FVS0812, da/dt is highest at 175 C when rates are expressed as a function of K. While crack growth rate is not controlled by C (sub t) at 175 C, da/dt appears to better correlate with C (sub t) at higher temperatures. Creep brittle cracking at intermediate temperatures, and perhaps related to strain aging, is augmented by time dependent transient creep plasticity at higher temperatures. The C (sub t) analysis is, however, complicated by the necessity to measure small differences in the elastic crack growth and creep contributions to the crack opening displacement rate. A microstructural study indicates that 2618 and FVS0812 are likely to be creep brittle materials, consistent with the results obtained from the fracture mechanics study. Time dependent crack growth of 2618 at 175 C is characterized by mixed transgranular and intergranular fracture. Delamination along the ribbon powder particle boundaries occurs in FVS0812 at all temperatures. The fracture mode of FVS0812 changes with temperature. At 175 C, it is characterized as dimpled rupture, and at 316 C as mixed matrix superplastic rupture and matrix-dipersoid debonding.

Leng, Yang↗
NASA NTRS2000

HOTCFGM-2D: A Coupled Higher-Order Theory for Cylindrical Structural Components with Bi-Directionally Components with Bi-Directionally Graded Microstructures

The objective of this two-year project was to develop and deliver to the NASA-Glenn Research Center a two-dimensional higher-order theory, and related computer codes, for the analysis and design of cylindrical functionally graded materials/structural components for use in advanced aircraft engines (e.g., combustor linings, rotor disks, heat shields, brisk blades). To satisfy this objective, two-dimensional version of the higher-order theory, HOTCFGM-2D, and four computer codes based on this theory, for the analysis and design of structural components functionally graded in the radial and circumferential directions were developed in the cylindrical coordinate system r-Theta-z. This version of the higher-order theory is a significant generalization of the one-dimensional theory, HOTCFGM-1D, developed during the FY97 for the analysis and design of cylindrical structural components with radially graded microstructures. The generalized theory is applicable to thin multi-phased composite shells/cylinders subjected to steady-state thermomechanical, transient thermal and inertial loading applied uniformly along the axial direction such that the overall deformation is characterized by a constant average axial strain. The reinforcement phases are uniformly distributed in the axial direction, and arbitrarily distributed in the radial and circumferential direction, thereby allowing functional grading of the internal reinforcement in the r-Theta plane. The four computer codes fgmc3dq.cylindrical.f, fgmp3dq.cylindrical.f, fgmgvips3dq.cylindrical.f, and fgmc3dq.cylindrical.transient.f are research-oriented codes for investigating the effect of functionally graded architectures, as well as the properties of the multi-phase reinforcement, in thin shells subjected to thermomechanical and inertial loading, on the internal temperature, stress and (inelastic) strain fields. The reinforcement distribution in the radial and circumferential directions is specified by the user. The thermal and inelastic properties of the individual phases can vary with temperature. The inelastic phases are presently modeled by the power-law creep model generalized to multi-directional loading (within fgmc3dq.cylindrical.f and fgmc3dq.cylindrical.transient.f for steady-state and transient thermal loading, respectively), and incremental plasticity and GVIPS unified viscoplasticity theories (within the steady-state loading versions fgmp3dq.cylindrical.f and fgmgvips3dq.cylindrical.f).

Pindera, Marek-Jerzy↗
DOE OSTI2023

Training material models using gradient descent algorithms

High temperature design requires accurate constitutive models to describe material inelastic deformation and failure behavior. Oftentimes, calibrating accurate models devolves into the problem of fitting the model parameters against experimental test data. Here, we present the pyopmat package, an open source framework for calibrating constitutive models against experiment data subjected to various loading conditions using machine learning techniques. The package calculates the exact gradient of the model response with respect to the parameters using a combination of automatic differentiation and the adjoint method. Given this exact gradient, we compare the performance of several gradient-based optimization techniques in fitting realistic constitutive models against data. Here, we demonstrate the efficiency and accuracy of our package through example problems using both synthetic data, generated using known parameter sets, under monotonic and cyclic loading conditions and also with an example applying the techniques developed here to actual high temperature creep-fatigue test data.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗
DOE OSTI2021

A Time-Dependent Directional Damage Theory for Brittle Rocks Considering the Kinetics of Microcrack Growth

In this work, a novel microcrack damage theory for describing the time-dependent behavior of brittle rocks is proposed. Instead of using scalar or tensorial variables to approximate the effect of microcracks, the concept of directional damage is introduced by upscaling the directional distribution density of microcracks through the noninteracting homogenization scheme. The contribution of an individual crack on global strain is first derived from the crack-opening displacements including sliding, closure, and dilation. The model is then cast in the free energy form, and the range of parameters ensuring compliance with the second law of thermodynamics is derived rigorously. Considering that the subcritical propagation of microcracks is the dominant mechanism behind the time-dependent deformation and failure of brittle rocks, the driving force for directional damage is defined as the corresponding energy release rate, and the effect of time is introduced through a damage evolution law inspired from the kinetics of subcritical crack growth. The model is applied to simulate the behavior of basalt rock under various loading conditions. In addition to predicting the stress–strain–time response, the model offers enhanced resolution in representing the anisotropic characteristics of microcrack-induced damage in brittle rocks and their time-dependent evolutions.

58 GEOSCIENCES↗
DOE OSTI2025

Microstructure and mechanical behavior of a TiC nanoprecipitate strengthened V Alloy

The V-4Cr-4Ti (V44) alloys have been proposed as the prime candidate structural material for self-cooled liquid Li blanket and other designs for fusion energy applications. However, the applications of the V44 alloy are limited to a narrow operation temperature window, due to reduction in creep strength at or above 700 °C and susceptibility to irradiation hardening and embrittlement when irradiated below 400 °C. Here, in this work, we explore the feasibility of designing a novel V alloy to form a high number density of TiC nanoprecipitates, in order to simultaneously improve creep strength and provide defect sinks to mitigate irradiation hardening. Computational thermodynamics was used to design a new alloy (V44C) to achieve our goal of high TiC nanoprecipitates density within the alloy V44 matrix. To ensure scalability, the new alloy was made through arc-melting and ingot-casting followed by hot forging, cold rolling and heat treatments of homogenization and precipitation aging. The microstructure was characterized by SEM, TEM, XRD and APT, confirming the existence of nanoprecipitates predicted in the thermodynamic calculations. In addition to microstructural evaluation tensile properties at room temperature and 700 °C, and Charpy impact energy at room temperature were measured. The microstructure and mechanical properties were then compared with those from a historic reference V44 alloy. The tensile strength improvement in V44C was rationalized based on particle and solid solution strengthening mechanism. The fracture behavior was discussed based on the fractography results and necking deformation behavior.

Alloys development↗
NASA NTRS2016

Hydrogen Embrittlement

Hydrogen embrittlement (HE) is a process resulting in a decrease in the fracture toughness or ductility of a metal due to the presence of atomic hydrogen. In addition to pure hydrogen gas as a direct source for the absorption of atomic hydrogen, the damaging effect can manifest itself from other hydrogen-containing gas species such as hydrogen sulfide (H2S), hydrogen chloride (HCl), and hydrogen bromide (HBr) environments. It has been known that H2S environment may result in a much more severe condition of embrittlement than pure hydrogen gas (H2) for certain types of alloys at similar conditions of stress and gas pressure. The reduction of fracture loads can occur at levels well below the yield strength of the material. Hydrogen embrittlement is usually manifest in terms of singular sharp cracks, in contrast to the extensive branching observed for stress corrosion cracking. The initial crack openings and the local deformation associated with crack propagation may be so small that they are difficult to detect except in special nondestructive examinations. Cracks due to HE can grow rapidly with little macroscopic evidence of mechanical deformation in materials that are normally quite ductile. This Technical Memorandum presents a comprehensive review of experimental data for the effects of gaseous Hydrogen Environment Embrittlement (HEE) for several types of metallic materials. Common material screening methods are used to rate the hydrogen degradation of mechanical properties that occur while the material is under an applied stress and exposed to gaseous hydrogen as compared to air or helium, under slow strain rates (SSR) testing. Due to the simplicity and accelerated nature of these tests, the results expressed in terms of HEE index are not intended to necessarily represent true hydrogen service environment for long-term exposure, but rather to provide a practical approach for material screening, which is a useful concept to qualitatively evaluate the severity of hydrogen embrittlement. The effects of hydrogen gas on mechanical properties such as tensile strength, ductility, fracture, low and high cycle fatigue, crack growth rate, and creep rupture are analyzed with respect to the general trends established from the HEE index values. It is observed that the severity of HE effects is also influenced by environmental factors such as pressure, temperature, and hydrogen gas purity. The severity of HE effects is also influenced by material factors such as surface finish, heat treatment, and product forms, compositions, grain direction, and crystal orientations.

Jonathan A. Lee↗
DOE OSTI2023

The effect of shear strain and shear localization on fault healing

The seismic cycle of repeated earthquake failure requires that faults regain frictional strength during the interseismic phase, when the fault is locked or undergoing quasi-static creep. Fault healing plays a central role in determining earthquake stress drop, recurrence interval, elastic radiation frequency and other source parameters. In particular, the longer a fault remains quasi-stationary, the stronger it becomes and the larger the potential stress drop can be for the next event. Here, we address the role of shear strain and strain localization on fault healing and healing rate. We performed slide-hold-slide friction experiments on quartz gouge in the double-direct shear configuration for shear strain up to 25 and hold times from 10 to 1000 s. The results show that both healing and healing rate increase nonlinearly with increasing shear strain. Frictional healing scales with volumetric strain within the laboratory fault zone. Using the volumetric strain upon reshear as a proxy for strain localization, we demonstrate that the capacity of a fault to heal is directly proportional to shear bandwidth and degree of strain localization. The more the deformation is localized, the higher are the healing and healing rate, and thus, the fault strength. Our data provide a framework for understanding variations in fault strength over the seismic cycle and the role of brecciation and strain localization on spatiotemporal variations in fault strength.

58 GEOSCIENCES↗
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↗