Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics

The way in which basal tractions, associated with mantle convection, couples with the lithosphere is a fundamental problem in geodynamics. A successful lithosphere-mantle coupling model for the Earth will satisfy observations of plate motions, intraplate stresses, and the plate boundary zone deformation. We solve the depth integrated three-dimensional force balance equations in a global finite element model that takes into account effects of both topography and shallow lithosphere structure as well as tractions originating from deeper mantle convection. The contribution from topography and lithosphere structure is estimated by calculating gravitational potential energy differences. The basal tractions are derived from a fully dynamic flow model with both radial and lateral viscosity variations. We simultaneously fit stresses and plate motions in order to delineate a best-fit lithosphere-mantle coupling model. We use both the World Stress Map and the Global Strain Rate Model to constrain the models. We find that a strongly coupled model with a stiff lithosphere and 3-4 orders of lateral viscosity variations in the lithosphere are best able to match the observational constraints. Our predicted deviatoric stresses, which are dominated by contribution from mantle tractions, range between 20-70 MPa. The best-fitting coupled models predict strain rates that are consistent with observations. That is, the intraplate areas are nearly rigid whereas plate boundaries and some other continental deformation zones display high strain rates. Comparison of mantle tractions and surface velocities indicate that in most areas tractions are driving, although in a few regions, including western North America, tractions are resistive. Citation: Ghosh, A., W. E. Holt, and L. M. Wen (2013), Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics.

Influence of system size on spatiotemporal dynamics of a model for plastic instability: Projecting low-dimensional and extensive chaos

This work is a continuation of our efforts to quantify the irregular scalar stress signals from the Ananthakrishna model for the Portevin-Le Chatelier instability observed under constant strain rate deformation conditions. Stress related to the spatial average of the dislocation activity is a dynamical variable that also determines the time evolution of dislocation densities. We carry out detailed investigations on the nature of spatiotemporal patterns of the model realized in the form of different types of dislocation bands seen in the entire instability domain and establish their connection to the nature of stress serrations. We then characterize the spatiotemporal dynamics of the model equations by computing the Lyapunov dimension as a function of the drive parameter. The latter scales with the system size only for low strain rates, where isolated dislocation bands are seen, and at high strain rates, where fully propagating bands are seen. At intermediate applied strain rates corresponding to the partially propagating bands, the Lyapunov dimension exhibits two distinct slopes, one for small system sizes and another for large. This feature is rationalized by demonstrating that the spatiotemporal patterns for small system sizes are altered from the partially propagating band types to isolated burst type. This in turn allows us to reconfirm that low-dimensional chaos is projected from the stress signals as long as there is a one-to-one correspondence between the bursts of dislocation bands and the stress drops. We then show that the stress signals in the regime of partially to fully propagative bands have features of extensive chaos by calculating the correlation dimension density. We also show that the correlation dimension density also depends on the system size. A number of issues related to the system size dependence of the Lyapunov dimension density and the correlation dimension density are discussed.

Evolution of microstructure during hot deformation of pearlitic steel

Hot deformation of pearlitic steel was carried out to examine the overall deformation response to microstructural evolution. To understand the mechanisms operative during hot deformation, compression tests were carried out at various temperatures in the range 400(-)600 degrees C and strain rates in the range 0.001-10 s(-1). The flow curves were analyzed to examine the occurrence of dynamic recrystallization. The evolution of microstructure in hot deformed samples is analysed using EBSD.

Projecting low and extensive dimensional chaos from spatio-temporal dynamics

We review the spatio-temporal dynamical features of the Ananthakrishna model for the Portevin-Le Chatelier effect, a kind of plastic instability observed under constant strain rate deformation conditions. We then establish a qualitative correspondence between the spatio-temporal structures that evolve continuously in the instability domain and the nature of the irregularity of the scalar stress signal. Rest of the study is on quantifying the dynamical information contained in the stress signals about the spatio-temporal dynamics of the model. We show that at low applied strain rates, there is a one-to-one correspondence with the randomly nucleated isolated bursts of mobile dislocation density and the stress drops. We then show that the model equations are spatio-temporally chaotic by demonstrating the number of positive Lyapunov exponents and Lyapunov dimension scale with the system size at low and high strain rates. Using a modified algorithm for calculating correlation dimension density, we show that the stress-strain signals at low applied strain rates corresponding to spatially uncorrelated dislocation bands exhibit features of low dimensional chaos. This is made quantitative by demonstrating that the model equations can be approximately reduced to space independent model equations for the average dislocation densities, which is known to be low-dimensionally chaotic. However, the scaling regime for the correlation dimension shrinks with increasing applied strain rate due to increasing propensity for propagation of the dislocation bands. The stress signals in the partially propagating to fully propagating bands turn to have features of extensive chaos.

A systematic investigation of fracture mechanisms in Al-Si based eutectic alloy-Effect of Si modification

The fracture characteristics of Al-Si based eutectic alloy are investigated in the unmodified and modified conditions under compression. The investigations are carried out at different strain rates and temperatures. Fracture of the alloy starts with eutectic Si particle fracture and modification plays an important role in particle fracture. The fraction of fractured particles is found to be always lesser in the modified condition than in the unmodified condition. Particle fracture increases with increase in strain. It is found that the Si particle fracture shows an increase with increase in strain rate and decreases with increase in temperature at 10% strain. Large and elongated particles show a greater tendency for fracture in the unmodified and modified conditions. Particle orientation plays an important role on fracture and the cracks are found to occur almost in a direction normal to the tensile strain imposed upon the particles by the deforming matrix in the unmodified alloy. The modified alloy shows a random distribution of fractured particles and crack orientation. The criteria of fracture based on dislocation pile-up mechanism and fiber loading explain the observed difference in particle fracture characteristics due to modification. The particle fracture for the modified alloy is also discussed in terms of Weibull statistics and the existing models of dispersion hardening. Particle/matrix interface decohesion is observed at higher strain rates and temperatures in the modified alloy. Dendritic rotation of 10 degrees is also observed at higher strain rates, which can increase the amount of particle fracture. (C) 2013 Elsevier B.V. All rights reserved.

Effect of fluid medium on mechanical behavior of carbon nanotube foam

This study reports the constitutive response and energy absorption capabilities of fluid-impregnated carbon nanotube (CNT) foams under compressive loading as a function of fluid viscosity and loading rates. At all strain rates tested, we observe two characteristic regimes: below a critical value, increasing fluid viscosity increases the load bearing and energy absorption capacities; after a critical value of the fluid's viscosity, we observe a rapid decrease in the systems' mechanical performance. For a given fluid viscosity, the load bearing capacity of the structure slightly decreases with strain rate. A phenomenological model, accounting for fluid-CNT interaction, is developed to explain the observed mechanical behavior. (C) 2014 AIP Publishing LLC.

Microstructure Development, Nanomechanical, and Dynamic Compression Properties of Spark Plasma Sintered TiB2-Ti-Based Homogeneous and Bi-layered Composites

The growing threats due to increased use of small-caliber armor piercing projectiles demand the development of new light-weight body armor materials. In this context, TiB2 appears to be a promising ceramic material. However, poor sinterability and low fracture toughness remain two major issues for TiB2. In order to address these issues together, Ti as a sinter-aid is used to develop TiB2-(x wt pct Ti), (x = 10, 20) homogeneous composites and a bi-layered composite (BLC) with each layer having Ti content of 10 and 20 wt pct. The present study uniquely demonstrates the efficacy of two-stage spark plasma sintering route to develop dense TiB2-Ti composites with an excellent combination of nanoscale hardness (similar to 36 GPa) and indentation fracture toughness (similar to 12 MPa m(1/2)). In case of BLC, these properties are not compromised w.r.t. homogeneous composites, suggesting the retention of baseline material properties even in the bi-layer design due to optimal relief of residual stresses. The better indentation toughness of TiB2-(10 wt pct Ti) and TiB2-(20 wt pct Ti) composites can be attributed to the observed crack deflection/arrest, indicating better damage tolerance. Transmission electron microscope investigation reveals the presence of dense dislocation networks and deformation twins in alpha-Ti at the grain boundaries and triple pockets, surrounded by TiB2 grains. The dynamic strength of around 4 GPa has been measured using Split Hopkinson Pressure Bar tests in a reproducible manner at strain rates of the order of 600 s(-1). The damage progression under high strain rate has been investigated by acquiring real time images for the entire test duration using ultra-high speed imaging. An attempt has been made to establish microstructure-property correlation and a simple analysis based on Mohr-Coulomb theory is used to rationalize the measured strength properties.

Microstructural Features of Hot Deformed Nb-1Zr-0.1C Alloy

Niobium-based alloys are well-established refractory materials; as a result of their high melting temperature and good creep properties, these alloys find their applications in nuclear reactors. The present study deals with a microstructural response of these materials during hot working. The evolution of microstructure and texture during high-temperature deformation has been investigated in the temperature range 1500-1700A degrees C and strain rate range of 0.001-0.1 s(-1). For each deformed sample, the microstructure has been examined in detail. The microstructural features clearly revealed the formation of a substructure and the occurrence of dynamic recrystallization in a proper temperature-strain rate window. At low strain rates, the necklace structure formation was more prominent.

An experimental study on loading rate effect on acoustic emission based b-values related to reinforced concrete fracture

This article reports on analysis of fracture processes in reinforced concrete (RC) beams with acoustic emission (AE) technique. An emphasis was given to study the effect of loading rate on variation in AE based b-values with the development of cracks in RC structures. RC beams of length 3.2 m were tested under load control at a rate of 4 kN/s, 5 kN/s and 6 kN/s and the b-value analysis available in seismology was used to study the fracture process in RC structures. Moreover, the b-value is related to the strain in steel to assess the damage state. It is observed that when the loading rate is higher, quick cracking development lead to rapid fluctuations and drops in the b-values. Also it is observed that concrete behaves relatively more brittle at higher loading rates (or at higher strain rates). The average b-values are lower as a few but larger amplitudes of AE events occur in contrast to more number of low amplitude AE events occur at low loading rates (or at low strain rates). (C) 2014 Elsevier Ltd. All rights reserved.

A critical examination of the paradox of strength and ductility in ultrafine-grained metals

The paradox of strength and ductility is now well established and denotes the difficulty of simultaneously achieving both high strength and high ductility. This paradox was critically examined using a cast Al-7% Si alloy processed by high-pressure torsion (HPT) for up to 10 turns at a temperature of either 298 or 445 K. This processing reduces the grain size to a minimum of similar to 0.4 mu m and also decreases the average size of the Si particles. The results show that samples processed to high numbers of HPT turns exhibit both high strength and high ductility when tested at relatively low strain rates and the strain rate sensitivity under these conditions is similar to 0.14 which suggests that flow occurs by some limited grain boundary sliding and crystallographic slip. The results are also displayed on the traditional diagram for strength and ductility and they demonstrate the potential for achieving high strength and high ductility by increasing the number of turns in HPT.

The contribution of grain boundary sliding in tensile deformation of an ultrafine-grained aluminum alloy having high strength and high ductility

An as-cast Al-7 % Si alloy was processed by high-pressure torsion (HPT) for up to 10 turns at temperatures of 298 or 445 K. The HPT-processed samples had ultrafine-grained structures and they were tested in tension at room temperature at various strain rates in the range from 1.0 x 10(-4) to 1.0 x 10(-2) s(-1). The contributions of grain boundary sliding (GBS) to the total strain were measured directly using atomic force microscopy. Samples simultaneously showing both high strength and high ductility contained the highest fractions of high-angle grain boundaries (HAGB) and exhibited the highest contributions from GBS, whereas samples showing high strength but low ductility gave negligible values for the sliding contributions. It is concluded that high strength and high ductility require both an ultrafine grain size and a high fraction of HAGB.

Compressive flow behavior of Al-Si based alloy: Role of heat treatment

The flow characteristics of a near eutectic Al-Si based cast alloy have been examined in compression at strain rates varying from 3 x 10(-4) to 10(2) s(-1) and at three different temperatures, i.e., room temperature (RT), 100 degrees C and 200 degrees C. The dependence of the flow behavior on heat treatment is studied by testing the alloy in non-heat treated (NHT) and heat treated (HT) conditions. The heat treatment has strong influence on strain rate sensitivity (SRS), strength and work hardening behavior of the alloy. It is observed that the strength of the alloy increases with increase in strain rate and it increases more rapidly above the strain rate of 10(-1) s(-1) in HT condition at all the temperatures, and at 100 degrees C and 200 degrees C in NHT condition. The thermally dependent process taking place in the HT matrix is responsible for the observed greater SRS in HT condition. The alloy in HT condition exhibits a larger work hardening rate than in NHT condition during initial stages of straining. However, the hardening rate decreases more sharply at higher strains in HT condition due to precipitate shearing and higher rate of Si particle fracture. Thermal hardening is observed at 200 degrees C in NHT condition due to precipitate formation, which results in increased SRS at higher temperatures. Thermal softening is observed in HT condition at 200 C due to precipitate coarsening, which leads to a decrease in SRS at higher temperatures. Stress simulations by a finite element method support the experimentally observed particle and matrix fracture behavior. A negative SRS and serrated flow are observed in the lower strain rate regime (3 x 10(-4)-10(-2) s(-1)) at RT and 100 degrees C, in both NHT and HT conditions. The observations show that both dynamic strain aging (DSA) and precipitate shearing play a role in serrated flow. (C) 2015 Elsevier B.V. All rights reserved.

Correlation between band propagation property and the nature of serrations in the Portevin-Le Chatelier effect

We investigate the correlation between the band propagation property and the nature and amplitude of serrations in the Portevin-Le Chatelier effect within the framework of the Ananthakrishna model. Several significant results emerge. First, we find that spatial and temporal correlations continuously increase with strain rate from type C to type A bands. Consequently, the nature of the bands also changes continuously from type C to A bands, and so do the changes in the associated serrations. Second, even the smallest extent of propagation induces small amplitude serrations. The spatial extent of band propagation is directly correlated with the duration of small amplitude serrations, a result that is consistent with recent experiments. This correspondence allows one to estimate the spatial extent of band propagation by just measuring the temporal stretch of small amplitude serrations. Therefore, this should be of practical value when only stress versus strain is recorded. Third, the average stress drop magnitude of the small amplitude serrations induced by the propagating bands remains small and nearly constant with strain rate. As a consequence, the fully propagating type A bands are in a state of criticality. We rationalize the increasing levels of spatial and temporal correlations found with increasing strain rates. Lastly, the model also predicts several band morphologies seen in experiments including the Luders-like propagating band. (C) 2015 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Dynamic compression behavior of reactive spark plasma sintered ultrafine grained (Hf, Zr)B-2-SiC composites

This paper reports the dynamic compression behavior of ultrafine grained (Hf, Zr)B-2-SiC composites, sintered using reactive spark plasma sintering at 1600 degrees C for 10 min. Dynamic strength of similar to 2.3 GPa has been measured using Split Hopkinson Pressure Bar (SHPB) tests in a reproducible manner at strain rates of 800-1300 s(-1). A comparison with competing boride based armor ceramics, in reference to the spectrum of properties evaluated, establishes the potential of (Hf, Zr)B-2-SiC composites for armor applications. (C) 2015 Elsevier Ltd and Techna Group S.r.l. All rights reserved.

Kinematic flow patterns in slow deformation of a dense granular material

The kinematic flow pattern in slow deformation of a model dense granular medium is studied at high resolution using in situ imaging, coupled with particle tracking. The deformation configuration is indentation by a flat punch under macroscopic plane-strain conditions. Using a general analysis method, velocity gradients and deformation fields are obtained from the disordered grain arrangement, enabling flow characteristics to be quantified. The key observations are the formation of a stagnation zone, as in dilute granular flow past obstacles; occurrence of vortices in the flow immediately underneath the punch; and formation of distinct shear bands adjoining the stagnation zone. The transient and steady state stagnation zone geometry, as well as the strength of the vortices and strain rates in the shear bands, are obtained from the experimental data. All of these results are well-reproduced in exact-scale non-smooth contact dynamics simulations. Full 3D numerical particle positions from the simulations allow extraction of flow features that are extremely difficult to obtain from experiments. Three examples of these, namely material free surface evolution, deformation of a grain column below the punch and resolution of velocities inside the primary shear band, are highlighted. The variety of flow features observed in this model problem also illustrates the difficulty involved in formulating a complete micromechanical analytical description of the deformation.