Grain coarsening maps: A new tool to predict microfabric evolution of polymineralic rocks

Polymineralic rocks undergo grain coarsening with increasing temperature in both static and deformational environments, as long as no mineral reactions occur. The grain coarsening in such rocks is complex because the different phases influence each other, and it is this interaction that controls the rate of grain coarsening of the entire aggregate. We present a mathematical approach to investigate coupled grain coarsening using a set of microstructural parameters, including grain size and volume fraction of both second phases and matrix mineral in combination with temperature information. Based on samples from polymineralic carbonate mylonites that were deformed at different temperatures, we demonstrate how the mathematical relation can be calibrated for this natural system. Using such data sets for other lithologies, grain coarsening maps can be generated, which allow the prediction of microstructural evolution in polymineralic rocks. Such predictions are crucial for all subdisciplines in the earth sciences that require fundamental knowledge about microstructural changes and rheology of an orogen at different depths, such as structural geology, geophysics, geodynamics, and metamorphic petrology.

Deformation mechanisms in second-phase affected microstructures and their energy balance

Based on the relationship Zener parameter (Z=second-phase size/second-phase volume fraction) vs. calcite grain size (dg), second-phase controlled aggregates and microstructures that are weakly affected by second-phases are discriminated. The latter are characterized by large but constant grain sizes, high calcite grain boundary fractions and crystallographic preferred orientations (CPO), while calcite grain size and calcite grain boundary fraction decrease continuously and CPO weakens with decreasing Z in second-phase controlled microstructures. These observations suggest that second-phase controlled microstructures predominantly deform via granular flow because pinning of calcite grain boundaries reduces the efficiency of dynamic recrystallization favoring mass transfer processes and grain boundary sliding. In contrast, the balance of grain size reduction and growth by dynamic recrystallization maintains a steady state grain size in microstructures that are only weakly affected by second-phases promoting a predominance of dislocation creep. With increasing temperature, the relationship between Z and dg persists but the calcite grain size increases continuously. Based on microstructures, the energy of each modifying process is calculated and its relative contribution is compared with energies of the competing processes (surface energy, dragging energy, dynamic recrystallization energy). The steady state microstructures result from a temperature-dependent energy minimization procedure of the system.

A new technique to automatically quantify microstructures of fine grained carbonate mylonites: two-step etching combined with SEM imaging and image analysis

A two-step etching technique for fine-grained calcite mylonites using 0.37% hydrochloric and 0.1% acetic acid produces a topographic relief which reflects the grain boundary geometry. With this technique, calcite grain boundaries become more intensely dissolved than their grain interiors but second phase minerals like dolomite, quartz, feldspars, apatite, hematite and pyrite are not affected by the acid and therefore form topographic peaks. Based on digital backscatter electron images and element distribution maps acquired on a scanning electron microscope, the geometry of calcite and the second phase minerals can be automatically quantified using image analysis software. For research on fine-grained carbonate rocks (e.g. dolomite calcite mixtures), this low-cost approach is an attractive alternative to the generation of manual grain boundary maps based on photographs from ultra-thin sections or orientation contrast images.

Extensional orogenic collapse captured by strike-slip tectonics:

The late Paleozoic collision between Gondwana and Laurussia resulted in the polyphase deformation and magmatism that characterizes the Iberian Massif of the Variscan orogen. In the Central Iberian Zone, initial con- tinental thickening (D1; folding and thrusting) was followed by extensional orogenic collapse (D2) responsible for the exhumation of high-grade rocks coeval to the emplacement of granitoids. This study presents a tectonometamorphic analysis of the Trancoso-Pinhel region (Central Iberian Zone) to ex- plain the processes in place during the transition froman extension-dominated state (D2) to a compression-dom- inated one (D3).Wereveal the existence of low-dipping D2 extensional structures later affected by several pulses of subhorizontal shortening, each of them typified by upright folds and strike-slip shearing (D3, D4 and D5, as identified by superimposition of structures). The D2 Pinhel extensional shear zone separates a low-grade domain from an underlying high-grade domain, and it contributed to the thermal reequilibration of the orogen by facil- itating heat advection from lower parts of the crust, crustal thinning, decompression melting, and magma intru- sion. Progressive lessening of the gravitational disequilibrium carried out by this D2 shear zone led to a switch from subhorizontal extension to compression and the eventual cessation and capture of the Pinhel shear zone by strike-slip tectonics during renewed crustal shortening. High-grade domains of the Pinhel shear zone were folded together with low-grade domains to define the current upright folded structure of the Trancoso-Pinhel re- gion, the D3 Tamames-Marofa-Sátão synform. Newdating of syn-orogenic granitoids (SHRIMP U\\Pb zircon dat- ing) intruding the Pinhel shear zone, together with the already published ages of early extensional fabrics constrain the functioning of this shear zone to ca. 331–311 Ma, with maximum tectonomagmatic activity at ca. 321–317 Ma. The capture and apparent cessation of movement of the Pinhel shear zone occurred at ca. 317– 311 Ma.

Thermomechanical analysis of strain localization in a ductile detachment zone

At the latitude of the Thor-Odin dome (British Columbia) the Columbia River Detachment defines the eastern margin of the Shuswap metamorphic core complex and localizes in a 1 km thick muscovite-bearing quartzite mylonite. We present a combined Ar-40/Ar-39, (micro) structural, and oxygen isotope study of the deformation history in the detachment and evaluate the spatial and temporal relationships between microstructure formation and localization of strain. High-precision Ar-40/Ar-39 geochronology from different levels in the mylonite delineates a pattern of increasingly younger (49.0 to 47.9 Ma) deformation ages in deeper levels of the mylonitic footwall. The correlation of Ar-40/Ar-39 ages with decreasing deformation temperatures (similar to 550 degrees-400 degrees C) in the top 200 m of the mylonite indicates that deformation migrated downward from the contact with the hanging wall. Strain localization was diachronous in progressively deeper levels of the footwall and was likely controlled by fluid-assisted strain hardening due to advective heat removal and contemporaneous reaction weakening due to dissolution-reprecipitation of white mica. The observed constant high-stress microstructures across the entire detachment indicate that flow stress was buffered by the interplay of strain rate and temperature, where high strain rates at elevated temperature produced the same microstructure as lower strain rates under decreasing temperature conditions. The combined data suggest that the complex interplay among temporally nonuniform rates of footwall exhumation, heat advection, and embrittlement by meteoric fluids strongly determines the thermomechanical behavior of extensional detachments.

Kinematics of oblique collision and ramping inferred from microstructures and strain in middle crustal rocks, central Southern Alps, New Zealand

The hanging wall of the Alpine Fault near Franz Josef Glacier has been exhumed during the past similar to2-3 m.y. providing a sample of the ductilely deformed middle crust of a modem obliquely convergent orogen. Presently exposed rocks of the Pacific Plate are inferred to have undergone several phases of ductile deformation as they moved westward above a mid-crustal detachment. Initially they were transpressed across the outboard part of the orogen, resulting in oblate fabrics with a down-dip stretch. Later, they encountered the Alpine Fault, experiencing an oblique-slip backshearing on vertical planes. This escalator-like deformation tilted and thinned the incoming crust onto that crustal-scale oblique ramp. This style of hanging wall deformation may affect only the most rapidly uplifting, central part of the Southern Alps because of the low flexural rigidity of the crust in that region and its displacement over a relatively sharp ramp-angle at depth. A 3D transpressive flow affected mylonites locally near the fault, but their shear direction remained parallel to plate motion, ruling out ductile 'extrusion' as an important process in this orogen. Outside the mylonite zone, late Cenozoic shortening is inferred to be modest (30-40%), as measured from deformation of younger biotite grains. Oblique collision is dominated by translation on the Alpine Fault, and rocks migrate rapidly through the deforming zone, preventing the accumulation of large finite strains. Transpression may play a minor role in oblique collision. (C) 2001 Elsevier Science Ltd. All rights reserved.

Topology optimization for designing strain-gauge load cells

Load cells are used extensively in engineering fields. This paper describes a novel structural optimization method for single- and multi-axis load cell structures. First, we briefly explain the topology optimization method that uses the solid isotropic material with penalization (SIMP) method. Next, we clarify the mechanical requirements and design specifications of the single- and multi-axis load cell structures, which are formulated as an objective function. In the case of multi-axis load cell structures, a methodology based on singular value decomposition is used. The sensitivities of the objective function with respect to the design variables are then formulated. On the basis of these formulations, an optimization algorithm is constructed using finite element methods and the method of moving asymptotes (MMA). Finally, we examine the characteristics of the optimization formulations and the resultant optimal configurations. We confirm the usefulness of our proposed methodology for the optimization of single- and multi-axis load cell structures.

IGFR-I expression and structural analysis of the hard palatine mucosa in an ethanol-drinking rat strain (UChA and UChB)

The study analyzed the effects of chronic alcohol ingestion on the ultrastructure of the lining epithelium of the hard palatine mucosa of rats UChA and UChB (lines with voluntary alcohol consumption) in order to contribute to the understanding of the consequences of alcohol abuse for the morphology of the digestive system. Thirty female adult animals aged 120 days were divided into three experimental groups. (1) Ten UChA rats (genetically low ethanol consumer) with voluntary intake of 10% v/v (5.45 g/kg/day) ethanol solution and water. (2) Ten UChB (genetically high ethanol consumer) rats with voluntary intake of 10% v/v (7.16 g/kg/day) ethanol solution and water. (3) Ten Wistar rats with voluntary ad libitum water intake (control group). Both groups received Nuvital pellets ad libitum. The IGFR-I expression was intense in both experimental groups. The epithelial cells of the alcoholic rats UChA and UChB showed many alterations such as the presence of lipid droplets, altered nuclei, nuclei in corneum layer and disrupted mitochondria. It was concluded that ethanol intake induces ultrastructural lesions in the hard palatine mucosa. (C) 2011 Elsevier Ltd. All rights reserved.

Predicting the mechanical behavior of amorphous polymeric materials under strain through multi-scale simulation

Polymeric materials have become the reference material for high reliability and performance applications. However, their performance in service conditions is difficult to predict, due in large part to their inherent complex morphology, which leads to non-linear and anisotropic behavior, highly dependent on the thermomechanical environment under which it is processed. In this work, a multiscale approach is proposed to investigate the mechanical properties of polymeric-based material under strain. To achieve a better understanding of phenomena occurring at the smaller scales, the coupling of a finite element method (FEM) and molecular dynamics (MD) modeling, in an iterative procedure, was employed, enabling the prediction of the macroscopic constitutive response. As the mechanical response can be related to the local microstructure, which in turn depends on the nano-scale structure, this multiscale approach computes the stress-strain relationship at every analysis point of the macro-structure by detailed modeling of the underlying micro- and meso-scale deformation phenomena. The proposed multiscale approach can enable prediction of properties at the macroscale while taking into consideration phenomena that occur at the mesoscale, thus offering an increased potential accuracy compared to traditional methods.

Fluid-rock-strain feedbacks in extensional shear zones

La présence de fluide météorique synchrone à l'activité du détachement (Farmin, 2003 ; Mulch et al., 2007 ; Gébelin et al., 2011), implique que les zones de cisaillement sont des systèmes ouverts avec des cellules de convections à l'échelle crustale et un intense gradient géothermique au sein du détachement (Morrison et Anderson, 1998, Gottardi et al., 2011). De plus, les réactions métamorphiques liées à des infiltrations fluides dans les zones de cisaillement extensionnel peuvent influencer les paramètres rhéologiques du système (White and Knipe, 1978), et impliquer la localisation de la déformation dans la croûte. Dans ce manuscrit, deux zones de cisaillement infiltrées par des fluides météoriques sont étudiées, l'une étant largement quartzitique, et l'autre de nature granitique ; les relations entre déformation, fluides, et roches s'appuient sur des approches structurales, microstructurales, chimiques et isotopiques. L'étude du détachement du Columbia river (WA, USA) met en évidence que la déformation mylonitique se développe en un million d'années. La phase de cisaillement principal s'effectue à 365± 30°C d'après les compositions isotopiques en oxygène du quartz et de la muscovite. Ces minéraux atteignent l'équilibre isotopique lors de leur recristallisation dynamique contemporaine à la déformation. La zone de cisaillement enregistre une baisse de température, remplaçant le mécanisme de glissement par dislocation par celui de dissolution- précipitation dans les derniers stades de l'activité du détachement. La dynamique de circulation fluide bascule d'une circulation pervasive à chenalisée, ce qui engendre localement la rupture des équilibres d'échange isotopiques. La zone de cisaillement de Bitterroot (MT, USA) présente une zone mylonitique de 600m d'épaisseur, progressant des protomylonites aux ultramylonites. L'intensité de la localisation de la déformation se reflète directement sur l'hydratation des feldspaths, réaction métamorphique majeure dite de « rock softening ». Une étude sur roche totale indique des transferts de masse latéraux au sein des mylonites, et d'importantes pertes de volume dans les ultramylonites. La composition isotopique en hydrogène des phyllosilicates met en évidence la présence (1) d'une source magmatique/métamorphique originelle, caractérisée par les granodiorites ayant conservé leur foliation magmatique, jusqu'aux protomylonites, et (2) une source météorique qui tamponne les valeurs des phyllosilicates des fabriques mylonitiques jusqu'aux veines de quartz non-déformées. Les compositions isotopiques en oxygène des minéraux illustrent le tamponnement de la composition du fluide météorique par l'encaissant. Ce phénomène cesse lors du processus de chloritisation de la biotite, puisque les valeurs des chlorites sont extrêmement négatives (-10 per mil). La thermométrie isotopique indique une température d'équilibre isotopique de la granodiorite entre 600-500°C, entre 500-300°C dans les mylonites, et entre 300 et 200°C dans les fabriques cassantes (cataclasites et veines de quartz). Basé sur les résultats issus de ce travail, nous proposons un modèle général d'interactions fluide-roches-déformation dans les zones de détachements infiltrées par des fluides météoriques. Les zones de détachements évoluent rapidement (en quelques millions d'années) au travers de la transition fragile-ductile ; celle-ci étant partiellement contrôlée par l'effet thermique des circulations de fluide météoriques. Les systèmes de détachements sont des lieux où la déformation et les circulations fluides sont couplées ; évoluant rapidement vers une localisation de la déformation, et de ce fait, une exhumation efficace. - The presence of meteoric fluids synchronous with the activity of extensional detachment zones (Famin, 2004; Mulch et al., 2007; Gébelin et al., 2011) implies that extensional systems involve fluid convection at a crustal scale, which results in high geothermal gradients within active detachment zones (Morrison and Anderson, 1998, Gottardi et al., 2011). In addition, the metamorphic reactions related to fluid infiltration in extensional shear zones can influence the rheology of the system (White and Knipe, 1978) and ultimately how strain localizes in the crust. In this thesis, two shear zones that were permeated by meteoric fluids are studied, one quartzite-dominated, and the other of granitic composition; the relations between strain, fluid, and evolving rock composition are addressed using structural, microstructural, and chemical/isotopic measurements. The study of the Columbia River detachment that bounds the Kettle core complex (Washington, USA) demonstrates that the mylonitic fabrics in the 100 m thick quartzite- dominated detachment footwall developed within one million years. The main shearing stage occurred at 365 ± 30°C when oxygen isotopes of quartz and muscovite equilibrated owing to coeval deformation and dynamic recrystallization of these minerals. The detachment shear zone records a decrease in temperature, and dislocation creep during detachment shearing gave way to dissolution-precipitation and fracturing in the later stages of detachment activity. Fluid flow switched from pervasive to channelized, leading to isotopic disequilibrium between different minerals. The Bitterroot shear zone detachment (Montana, USA) developed a 600 m thick mylonite zone, with well-developed transitions from protomylonite to ultramylonite. The localization of deformation relates directly to the intensity of feldspar hydration, a major rock- softening metamorphic reaction. Bulk-rock analyses of the mylonitic series indicate lateral mass transfer in the mylonite (no volume change), and significant volume loss in ultramylonite. The hydrogen isotope composition of phyllosilicates shows (1) the presence of an initial magmatic/metamorphic source characterized by the granodiorite in which a magmatic, and gneissic (protomylonite) foliation developed, and (2) a meteoric source that buffers the values of phyllosilicates in mylonite, ultramylonite, cataclasite, and deformed and undeformed quartz veins. The mineral oxygen isotope compositions were buffered by the host-rock compositions until chloritization of biotite started; the chlorite oxygen isotope values are negative (-10 per mil). Isotope thermometry indicates a temperature of isotopic equilibrium of the granodiorite between 600-500°C, between 500-300°C in the mylonite, and between 300 and 200°C for brittle fabrics (cataclasite and quartz veins). Results from this work suggest a general model for fluid-rock-strain feedbacks in detachment systems that are permeated by meteoric fluids. Phyllosilicates have preserved in their hydrogen isotope values evidence for the interaction between rock and meteoric fluids during mylonite development. Fluid flow generates mass transfer along the tectonic anisotropy, and mylonites do not undergo significant volume change, except locally in ultramylonite zones. Hydration of detachment shear zones attends mechanical grain size reduction and enhances strain softening and localization. Self-exhuming detachment shear zones evolve rapidly (a few million years) through the transition from ductile to brittle, which is partly controlled by the thermal effect of circulating surface fluids. Detachment systems are zones in the crust where strain and fluid flow are coupled; these systems. evolve rapidly toward strain localization and therefore efficient exhumation.

A Practical Micro Mechanical Model Based Local Approach Methodology for the Analysis of Ductile Fracture of Welded T-Joints

This thesis concentrates on developing a practical local approach methodology based on micro mechanical models for the analysis of ductile fracture of welded joints. Two major problems involved in the local approach, namely the dilational constitutive relation reflecting the softening behaviour of material, and the failure criterion associated with the constitutive equation, have been studied in detail. Firstly, considerable efforts were made on the numerical integration and computer implementation for the non trivial dilational Gurson Tvergaard model. Considering the weaknesses of the widely used Euler forward integration algorithms, a family of generalized mid point algorithms is proposed for the Gurson Tvergaard model. Correspondingly, based on the decomposition of stresses into hydrostatic and deviatoric parts, an explicit seven parameter expression for the consistent tangent moduli of the algorithms is presented. This explicit formula avoids any matrix inversion during numerical iteration and thus greatly facilitates the computer implementation of the algorithms and increase the efficiency of the code. The accuracy of the proposed algorithms and other conventional algorithms has been assessed in a systematic manner in order to highlight the best algorithm for this study. The accurate and efficient performance of present finite element implementation of the proposed algorithms has been demonstrated by various numerical examples. It has been found that the true mid point algorithm (a = 0.5) is the most accurate one when the deviatoric strain increment is radial to the yield surface and it is very important to use the consistent tangent moduli in the Newton iteration procedure. Secondly, an assessment of the consistency of current local failure criteria for ductile fracture, the critical void growth criterion, the constant critical void volume fraction criterion and Thomason's plastic limit load failure criterion, has been made. Significant differences in the predictions of ductility by the three criteria were found. By assuming the void grows spherically and using the void volume fraction from the Gurson Tvergaard model to calculate the current void matrix geometry, Thomason's failure criterion has been modified and a new failure criterion for the Gurson Tvergaard model is presented. Comparison with Koplik and Needleman's finite element results shows that the new failure criterion is fairly accurate indeed. A novel feature of the new failure criterion is that a mechanism for void coalescence is incorporated into the constitutive model. Hence the material failure is a natural result of the development of macroscopic plastic flow and the microscopic internal necking mechanism. By the new failure criterion, the critical void volume fraction is not a material constant and the initial void volume fraction and/or void nucleation parameters essentially control the material failure. This feature is very desirable and makes the numerical calibration of void nucleation parameters(s) possible and physically sound. Thirdly, a local approach methodology based on the above two major contributions has been built up in ABAQUS via the user material subroutine UMAT and applied to welded T joints. By using the void nucleation parameters calibrated from simple smooth and notched specimens, it was found that the fracture behaviour of the welded T joints can be well predicted using present methodology. This application has shown how the damage parameters of both base material and heat affected zone (HAZ) material can be obtained in a step by step manner and how useful and capable the local approach methodology is in the analysis of fracture behaviour and crack development as well as structural integrity assessment of practical problems where non homogeneous materials are involved. Finally, a procedure for the possible engineering application of the present methodology is suggested and discussed.

Structural studies in the Southern Province, South of Sudbury, Ontario

Rocks correlated with the Hough Lake and Quirke Lake Groups of the Huronian Supergroup form part of a northeasterly trending corridor that separates 1750 Ma granitic intrusive rocks of the Chief Lake batholith from the 1850 Ma mafic intrusive rocks of the Sudbury Igneous Complex. This corridor is dissected by two major structural features; the Murray Fault Zone (MFZ) and the Long Lake Fault (LLF). Detailed structural mapping and microstructural analysis indicates that the LLF, which has juxtaposed Huronian rocks of different deformation style and metamorphism grade, was a more significant plane of dislocation than the MFZ. The sense of displacement along the LLF is high angle reverse in which rocks to the southeast have been raised relative to those in the northwest. South of the LLF Huronian rocks underwent ductile defonnation at amphibolite facies conditions. The strain was constrictional, defined by a triaxial strain ellipsoid in which X > Y > z. Calculations of a regional k value were approximately 1.3. Penetrative ductile defonnation resulted in the development of a preferred crystallographic orientation in quartz as well as the elongation of quartz grains to fonn a regional southeast-northwest trending, subvertical lineation. Similar lithologies north of the LLF underwent dominantly brittle deformation under greenschist facies conditions. Deformation north of the LLF is characterized by the thrusting of structural blocks to form angular discordances in bedding orientation which were previously interpreted as folds. Ductile deformation occurred between 1750 and 1238 Ma and is correlated with a regional period of south over north reverse faulting that effected much of the southern Sudbury region. Post dating the reverse faulting event was a period of sedimentation as a conglomerate unit was deposited on vertically bedded Huronian rocks. Rocks in the study area were intruded by both mafic and felsic dykes. The 1238 Ma mafic dykes appear to have been offset during a period of dextral strike slip displacement along the major fault'). Indirect evidence indicates that this event occurred after the thrusting at 950 to 1100 Ma associated with the Grenvillian Orogeny.