Random vibration of the functionally graded laminates in thermal environments

This work deals with the random free vibration of functionally graded laminates with general boundary conditions and subjected to a temperature change, taking into account the randomness in a number of independent input variables such as Young's modulus, Poisson's ratio and thermal expansion coefficient of each constituent material. Based on third-order shear deformation theory, the mixed-type formulation and a semi-analytical approach are employed to derive the standard eigenvalue problem in terms of deflection, mid-plane rotations and stress function. A mean-centered first-order perturbation technique is adopted to obtain the second-order statistics of vibration frequencies. A detailed parametric study is conducted, and extensive numerical results are presented in both tabular and graphical forms for laminated plates that contain functionally graded material which is made of aluminum and zirconia, showing the effects of scattering in thermo-clastic material constants, temperature change, edge support condition, side-to-thickness ratio, and plate aspect ratio on the stochastic characteristics of natural frequencies. (c) 2005 Elsevier B.V. All rights reserved.

Effect of rolling on dissipation in fault gouges

Sliding and rolling are two outstanding deformation modes in granular media. The first one induces frictional dissipation whereas the latter one involves deformation with negligible resistance. Using numerical simulations on two-dimensional shear cells, we investigate the effect of the grain rotation on the energy dissipation and the strength of granular materials under quasistatic shear deformation. Rolling and sliding are quantified in terms of the so-called Cosserat rotations. The observed spontaneous formation of vorticity cells and clusters of rotating bearings may provide an explanation for the long standing heat flow paradox of earthquake dynamics.

Vibration of symmetrically laminated thick super elliptical plates

This paper reports a free vibration analysis of thick plates with rounded corners subject to a free, simply-supported or clamped boundary condition. The plate perimeter is defined by a super elliptic function with a power defining the shape ranging from an ellipse to a rectangle. To incorporate transverse shear deformation, the Reddy third-order plate theory is employed. The energy integrals incorporating shear deformation and rotary inertia are formulated and the p-Ritz procedures are used to derive the governing eigenvalue equation. Numerical examples for plates with different shapes and boundary conditions are solved and their frequency parameters, where possible, are compared with known results. Parametric studies are carried out to show the sensitivities of frequency parameters by varying the geometry, fibre stacking sequence, and boundary condition. (C) 1999 Academic Press.

First principles calculation of electron-phonon and alloy scattering in strained SiGe

First-principles electronic structure methods are used to predict the mobility of n-type carrier scattering in strained SiGe. We consider the effects of strain on the electron-phonon deformation potentials and the alloy scattering parameters. We calculate the electron-phonon matrix elements and fit them up to second order in strain. We find, as expected, that the main effect of strain on mobility comes from the breaking of the degeneracy of the six Δ and L valleys, and the choice of transport direction. The non-linear effects on the electron-phonon coupling of the Δ valley due to shear strain are found to reduce the mobility of Si-like SiGe by 50% per % strain. We find increases in mobility between 2 and 11 times that of unstrained SiGe for certain fixed Ge compositions, which should enhance the thermoelectric figure of merit in the same order, and could be important for piezoresistive applications.

Evolution of a polymineralic mantle shear zone and the role of second phases in the localization of deformation

The influence of second phases (e.g., pyroxenes) on olivine grain size was studied by quantitative microfabric analyses of samples of the Hilti massif mantle shear zone (Semail ophiolite, Oman). The microstructures range from porphyroclastic tectonites to ultramylonites, from outside to the center of the shear zone. Starting at conditions of ridge-related flow, they formed under continuous cooling leading to progressive strain localization. The dependence of the average olivine grain size on the second-phase content can be split into a second-phase controlled and a dynamic recrystallization-controlled field. In the former, the olivine grain size is related to the ratio between the second-phase grain size and volume fraction (Zener parameter). In the latter, dynamic recrystallization manifested by a balance between grain growth and grain size reduction processes yields a stable olivine grain size. In both fields the average olivine and second-phase grain size decreases with decreasing temperature. Combining the microstructural information with deformation mechanism maps suggests that the porphyroclastic tectonites (similar to 1100 degrees C) and mylonites (similar to 800 degrees C) formed under the predominance of dislocation creep. Since olivine-rich layers are intercalated with layer parallel, polymineralic bands in the mylonites, nearly equiviscous conditions can be assumed. In the ultramylonites, diffusion creep represents the major deformation mechanism in the polymineralic layers. It is this switch in deformation mechanism from dislocation creep to diffusion creep that forces strain to localize in the fine-grained polymineralic domains at low temperatures (<similar to 700 degrees C), underlining the role of the second phases on strain localization in cooling mantle rocks.

TECTONICS OF THE HELVETIC NAPPE SYSTEM (W SWITZERLAND): CONTROL OF STRAIN LOCALIZATION AND BASEMENT-COVER DEFORMATION. Insights from models based on continuum mechanics

The Helvetic nappe system in Western Switzerland is a stack of fold nappes and thrust sheets em-placed at low grade metamorphism. Fold nappes and thrust sheets are also some of the most common features in orogens. Fold nappes are kilometer scaled recumbent folds which feature a weakly deformed normal limb and an intensely deformed overturned limb. Thrust sheets on the other hand are characterized by the absence of overturned limb and can be defined as almost rigid blocks of crust that are displaced sub-horizontally over up to several tens of kilometers. The Morcles and Doldenhom nappe are classic examples of fold nappes and constitute the so-called infra-Helvetic complex in Western and Central Switzerland, respectively. This complex is overridden by thrust sheets such as the Diablerets and Wildhörn nappes in Western Switzerland. One of the most famous example of thrust sheets worldwide is the Glariis thrust sheet in Central Switzerland which features over 35 kilometers of thrusting which are accommodated by a ~1 m thick shear zone. Since the works of the early Alpine geologist such as Heim and Lugeon, the knowledge of these nappes has been steadily refined and today the geometry and kinematics of the Helvetic nappe system is generally agreed upon. However, despite the extensive knowledge we have today of the kinematics of fold nappes and thrust sheets, the mechanical process leading to the emplacement of these nappe is still poorly understood. For a long time geologist were facing the so-called 'mechanical paradox' which arises from the fact that a block of rock several kilometers high and tens of kilometers long (i.e. nappe) would break internally rather than start moving on a low angle plane. Several solutions were proposed to solve this apparent paradox. Certainly the most successful is the theory of critical wedges (e.g. Chappie 1978; Dahlen, 1984). In this theory the orogen is considered as a whole and this change of scale allows thrust sheet like structures to form while being consistent with mechanics. However this theoiy is intricately linked to brittle rheology and fold nappes, which are inherently ductile structures, cannot be created in these models. When considering the problem of nappe emplacement from the perspective of ductile rheology the problem of strain localization arises. The aim of this thesis was to develop and apply models based on continuum mechanics and integrating heat transfer to understand the emplacement of nappes. Models were solved either analytically or numerically. In the first two papers of this thesis we derived a simple model which describes channel flow in a homogeneous material with temperature dependent viscosity. We applied this model to the Morcles fold nappe and to several kilometer-scale shear zones worldwide. In the last paper we zoomed out and studied the tectonics of (i) ductile and (ii) visco-elasto-plastic and temperature dependent wedges. In this last paper we focused on the relationship between basement and cover deformation. We demonstrated that during the compression of a ductile passive margin both fold nappes and thrust sheets can develop and that these apparently different structures constitute two end-members of a single structure (i.e. nappe). The transition from fold nappe to thrust sheet is to first order controlled by the deformation of the basement. -- Le système des nappes helvétiques en Suisse occidentale est un empilement de nappes de plis et de nappes de charriage qui se sont mis en place à faible grade métamorphique. Les nappes de plis et les nappes de charriage sont parmi les objets géologiques les plus communs dans les orogènes. Les nappes de plis sont des plis couchés d'échelle kilométrique caractérisés par un flanc normal faiblement défor-mé, au contraire de leur flanc inverse, intensément déformé. Les nappes de charriage, à l'inverse se caractérisent par l'absence d'un flanc inverse bien défini. Elles peuvent être définies comme des blocs de croûte terrestre qui se déplacent de manière presque rigide qui sont déplacés sub-horizontalement jusqu'à plusieurs dizaines de kilomètres. La nappe de Mordes et la nappe du Doldenhorn sont des exemples classiques de nappes de plis et constitue le complexe infra-helvétique en Suisse occidentale et centrale, respectivement. Ce complexe repose sous des nappes de charriages telles les nappes des Diablerets et du Widlhörn en Suisse occidentale. La nappe du Glariis en Suisse centrale se distingue par un déplacement de plus de 35 kilomètres qui s'est effectué à la faveur d'une zone de cisaillement basale épaisse de seulement 1 mètre. Aujourd'hui la géométrie et la cinématique des nappes alpines fait l'objet d'un consensus général. Malgré cela, les processus mécaniques par lesquels ces nappes se sont mises en place restent mal compris. Pendant toute la première moitié du vingtième siècle les géologues les géologues ont été confrontés au «paradoxe mécanique». Celui-ci survient du fait qu'un bloc de roche haut de plusieurs kilomètres et long de plusieurs dizaines de kilomètres (i.e., une nappe) se fracturera de l'intérieur plutôt que de se déplacer sur une surface frictionnelle. Plusieurs solutions ont été proposées pour contourner cet apparent paradoxe. La solution la plus populaire est la théorie des prismes d'accrétion critiques (par exemple Chappie, 1978 ; Dahlen, 1984). Dans le cadre de cette théorie l'orogène est considéré dans son ensemble et ce simple changement d'échelle solutionne le paradoxe mécanique (la fracturation interne de l'orogène correspond aux nappes). Cette théorie est étroitement lié à la rhéologie cassante et par conséquent des nappes de plis ne peuvent pas créer au sein d'un prisme critique. Le but de cette thèse était de développer et d'appliquer des modèles basés sur la théorie de la méca-nique des milieux continus et sur les transferts de chaleur pour comprendre l'emplacement des nappes. Ces modèles ont été solutionnés de manière analytique ou numérique. Dans les deux premiers articles présentés dans ce mémoire nous avons dérivé un modèle d'écoulement dans un chenal d'un matériel homogène dont la viscosité dépend de la température. Nous avons appliqué ce modèle à la nappe de Mordes et à plusieurs zone de cisaillement d'échelle kilométrique provenant de différents orogènes a travers le monde. Dans le dernier article nous avons considéré le problème à l'échelle de l'orogène et avons étudié la tectonique de prismes (i) ductiles, et (ii) visco-élasto-plastiques en considérant les transferts de chaleur. Nous avons démontré que durant la compression d'une marge passive ductile, a la fois des nappes de plis et des nappes de charriages peuvent se développer. Nous avons aussi démontré que nappes de plis et de charriages sont deux cas extrêmes d'une même structure (i.e. nappe) La transition entre le développement d'une nappe de pli ou d'une nappe de charriage est contrôlé au premier ordre par la déformation du socle. -- Le système des nappes helvétiques en Suisse occidentale est un emblement de nappes de plis et de nappes de chaînage qui se sont mis en place à faible grade métamoiphique. Les nappes de plis et les nappes de charriage sont parmi les objets géologiques les plus communs dans les orogènes. Les nappes de plis sont des plis couchés d'échelle kilométrique caractérisés par un flanc normal faiblement déformé, au contraire de leur flanc inverse, intensément déformé. Les nappes de charriage, à l'inverse se caractérisent par l'absence d'un flanc inverse bien défini. Elles peuvent être définies comme des blocs de croûte terrestre qui se déplacent de manière presque rigide qui sont déplacés sub-horizontalement jusqu'à plusieurs dizaines de kilomètres. La nappe de Morcles and la nappe du Doldenhorn sont des exemples classiques de nappes de plis et constitue le complexe infra-helvétique en Suisse occidentale et centrale, respectivement. Ce complexe repose sous des nappes de charriages telles les nappes des Diablerets et du Widlhörn en Suisse occidentale. La nappe du Glarüs en Suisse centrale est certainement l'exemple de nappe de charriage le plus célèbre au monde. Elle se distingue par un déplacement de plus de 35 kilomètres qui s'est effectué à la faveur d'une zone de cisaillement basale épaisse de seulement 1 mètre. La géométrie et la cinématique des nappes alpines fait l'objet d'un consensus général parmi les géologues. Au contraire les processus physiques par lesquels ces nappes sont mises en place reste mal compris. Les sédiments qui forment les nappes alpines se sont déposés à l'ère secondaire et à l'ère tertiaire sur le socle de la marge européenne qui a été étiré durant l'ouverture de l'océan Téthys. Lors de la fermeture de la Téthys, qui donnera naissance aux Alpes, le socle et les sédiments de la marge européenne ont été déformés pour former les nappes alpines. Le but de cette thèse était de développer et d'appliquer des modèles basés sur la théorie de la mécanique des milieux continus et sur les transferts de chaleur pour comprendre l'emplacement des nappes. Ces modèles ont été solutionnés de manière analytique ou numérique. Dans les deux premiers articles présentés dans ce mémoire nous nous sommes intéressés à la localisation de la déformation à l'échelle d'une nappe. Nous avons appliqué le modèle développé à la nappe de Morcles et à plusieurs zones de cisaillement provenant de différents orogènes à travers le monde. Dans le dernier article nous avons étudié la relation entre la déformation du socle et la défonnation des sédiments. Nous avons démontré que nappe de plis et nappes de charriages constituent les cas extrêmes d'un continuum. La transition entre nappe de pli et nappe de charriage est intrinsèquement lié à la déformation du socle sur lequel les sédiments reposent.

A two-phase composite in simple shear: Effective mechanical anisotropy development and localization potential

We present a combined shape and mechanical anisotropy evolution model for a two-phase inclusion-bearing rock subject to large deformation. A single elliptical inclusion embedded in a homogeneous but anisotropic matrix is used to represent a simplified shape evolution enforced on all inclusions. The mechanical anisotropy develops due to the alignment of elongated inclusions. The effective anisotropy is quantified using the differential effective medium (DEM) approach. The model can be run for any deformation path and an arbitrary viscosity ratio between the inclusion and host phase. We focus on the case of simple shear and weak inclusions. The shape evolution of the representative inclusion is largely insensitive to the anisotropy development and to parameter variations in the studied range. An initial hardening stage is observed up to a shear strain of gamma = 1 irrespective of the inclusion fraction. The hardening is followed by a softening stage related to the developing anisotropy and its progressive rotation toward the shear direction. The traction needed to maintain a constant shear rate exhibits a fivefold drop at gamma = 5 in the limiting case of an inviscid inclusion. Numerical simulations show that our analytical model provides a good approximation to the actual evolution of a two-phase inclusion-host composite. However, the inclusions develop complex sigmoidal shapes resulting in the formation of an S-C fabric. We attribute the observed drop in the effective normal viscosity to this structural development. We study the localization potential in a rock column bearing varying fraction of inclusions. In the inviscid inclusion case, a strain jump from gamma = 3 to gamma = 100 is observed for a change of the inclusion fraction from 20% to 33%.

Small quantity but large effect — How minor phases control strain localization in upper mantle shear zones

Low viscosity domains such as localized shear zones exert an important control on the geodynamics of the uppermost mantle. Grain size reduction and subsequent strain localization related to a switch from dislocation to diffusion creep is one mechanism to form low viscosity domains. To sustain strain localization, the grain size of mantle minerals needs to be kept small over geological timescales. One way to keep olivine grain sizes small is by pinning of mobile grain boundaries during grain growth by other minerals (second phases). Detailed microstructural studies based on natural samples from three shear zones formed at different geodynamic settings, allowed the derivation of the olivine grain-size dependence on the second-phase content. The polymineralic olivine grain-size evolution with increasing strain is similar in the three shear zones. If the second phases are to pin the mobile olivine grain boundary the phases need to be well mixed before grain growth. We suggest that melt-rock and metamorphic reactions are crucial for the initial phase mixing in mantle rocks. With ongoing deformation and increasing strain, grain boundary sliding combined with mass transfer processes and nucleation of grains promotes phase mixing resulting in fine-grained polymineralic mixtures that deform by diffusion creep. Strain localization due to the presence of volumetrically minor minerals in polymineralic mantle rocks is only important at high strain deformation (ultramylonites) at low temperatures (<~800°C). At smaller strain and stress conditions and/or higher temperatures other parameters like overall energy available to deform a given rock volume, the inheritance of mechanical anisotropies or the presence of water or melts needs to be considered to explain strain localization in the upper mantle.

The variability of peridotite composition across a mantle shear zone (Lanzo massif, Italy): interplay of melt focusing and deformation

In this paper we present new data on the spatial variability of peridotite composition across a kilometer-scale mantle shear zone within the Lanzo massif (Western Alps, Italy). The shear zone separates the central from the northern part of the massif. Plagioclase peridotite shows gradually increasing deformation towards the shear zone, from porphyroclastic to mylonitic textures in the central body, while the northern body is composed of porphyroclastic rocks. The peridotite displays a large range of compositions, from fertile peridotite to refractory harzburgite and dunite. Deformed peridotites (proto-mylonite and mylonites) tend to be compositionally more homogeneous and fertile than weakly deformed peridotites. The composition of most plagioclase peridotites show rather high and constant (Ce/Yb) (N) ratios, and Yb (N) that cannot be explained by any simple melting model. Instead, refertilization modeling, consisting of melt increments from spinel peridotite sources, particularly with E-MORB melt, reasonably reproduces the plagioclase peridotite whole rock composition. Combined with constraints from Ce-Nb and Ce-Th systematics, we speculate that peridotites such as those from Lanzo record pervasive refertilization processes in the thermal boundary layer. In this scenario, mantle shear zones might act as important areas of melt focusing in the upper mantle that separates the thermal boundary layer from the conductively cooled mantle.

Microfabrics in carbonate mylonites along a large-scale shear zone (Helvetic Alps)

Carbonate mylonites with varying proportions of second-phase minerals were collected at positions of increasing metamorphic grade along the basal thrust of the Morcles nappe (Helvetic nappes, Switzerland). Variations of temperature, stress, and strain rate, changes in chemistry of solid and fluid phases, and differing degrees of strain localization and annealing were tracked by measuring the shapes, mean sizes, and size distributions of both matrix and second-phase grains, as well as crystal preferred orientation (CPO) of the matrix. Field structures suggest that strain rate was constant along the fault. The mean and distribution of the calcite grain sizes were affected most profoundly by temperature: Increased temperature, presumably accompanied by decreased stress, correlated with larger mean sizes and wider size distributions. At a given location, the matrix grains in mylonites with more second-phase particles are, on average, smaller, have narrower size distributions, and have more elongate shapes. For example, mylonites with 50 vol.% of second phases have matrix grain sizes half that of pure mylonites. Changes in calcite chemistry and the presence of synkinematic fluids seemed to influence microfabric only weakly. Temporal variations in conditions, such as exhumation-induced cooling, apparently provoke changes in temperature, stress, and strain rate along the nappe. These changes result in further strain localization during retrograde conditions and cause the grain size to be reduced by an additional 50%. The matrix CPO strengthens with increasing temperature or strain, but weakens and rotates with increasing second-phase content, These fabric changes suggest differing rates of grain growth, grain size reduction, and development of CPO owing to variations in the deformation conditions and, perhaps, mechanisms. To interpret natural mylonite structures or to extrapolate mechanical data to natural situations requires careful characterization of the microfabric, and, in particular, second-phase minerals. (c) 2007 Elsevier B.V, All rights reserved.

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.

Juxtaposition of melt impregnation and high-temperature shear zones in the upper mantle; field and petrological constraints from the Lanzo peridotite (Northern Italy)

Results of a field and microstructural study between the northern and the central bodies of the Lanzo plagioclase peridotite massif (NW Italy) indicate that the spatial distribution of deformation is asymmetric across kilometre-scale mantle shear zones. The southwestern part of the shear zone (footwall) shows a gradually increasing degree of deformation from porphyroclastic peridotites to mylonite, whereas the northeastern part (hanging wall) quickly grades into weakly deformed peridotites. Discordant gabbroic and basaltic dykes are asymmetrically distributed and far more abundant in the footwall of the shear zone. The porphyroclastic peridotite displays porphyroclastic zones and domains of igneous crystallization whereas mylonites are characterized by elongated porphyroclasts, embedded between fine-grained, polycrystalline bands of olivine, plagioclase, clinopyroxene, orthopyroxene, spinel, rare titanian pargasite, and domains of recrystallized olivine. Two types of melt impregnation textures have been found: (1) clinopyroxene porphyroclasts incongruently reacted with migrating melt to form orthopyroxene plagioclase; (2) olivine porphyroclasts are partially replaced by interstitial orthopyroxene. The meltrock reaction textures tend to disappear in the mylonites, indicating that deformation in the mylonite continued under subsolidus conditions. The pyroxene chemistry is correlated with grain size. High-Al pyroxene cores indicate high temperatures (11001030C), whereas low-Al neoblasts display lower final equilibration temperatures (860C). The spinel Cr-number [molar Cr/(Cr Al)] and TiO2 concentrations show extreme variability covering almost the entire range known from abyssal peridotites. The spinel compositions of porphyroclastic peridotites from the central body are more variable than spinel from mylonite, mylonite with ultra-mylonite bands, and porphyroclastic rocks of the northern body. The spinel compositions probably indicate disequilibrium and would favour rapid cooling, and a faster exhumation of the central peridotite body, relative to the northern one. Our results indicate that melt migration and high-temperature deformation are juxtaposed both in time and space. Meltrock reaction may have caused grain-size reduction, which in turn led to localization of deformation. It is likely that melt-lubricated, actively deforming peridotites acted as melt focusing zones, with permeabilities higher than the surrounding, less deformed peridotites. Later, under subsolidus conditions, pinning in polycrystalline bands in the mylonites inhibited substantial grain growth and led to permanent weak zones in the upper mantle peridotite, with a permeability that is lower than in the weakly deformed peridotites. Such an inversion in permeability might explain why actively deforming, fine-grained peridotite mylonite acted as a permeability barrier and why ascending mafic melts might terminate and crystallize as gabbros along actively deforming shear zones. Melt-lubricated mantle shear zones provide a mechanism for explaining the discontinuous distribution of gabbros in oceancontinent transition zones, oceanic core complexes and ultraslow-spreading ridges.

A computationally efficient shear deformable beam element for large deformation multibody applications

In this paper, a new two-dimensional shear deformable beam element based on the absolute nodal coordinate formulation is proposed. The nonlinear elastic forces of the beam element are obtained using a continuum mechanics approach without employing a local element coordinate system. In this study, linear polynomials are used to interpolate both the transverse and longitudinal components of the displacement. This is different from other absolute nodal-coordinate-based beam elements where cubic polynomials are used in the longitudinal direction. The accompanying defects of the phenomenon known as shear locking are avoided through the adoption of selective integration within the numerical integration method. The proposed element is verified using several numerical examples, and the results are compared to analytical solutions and the results for an existing shear deformable beam element. It is shown that by using the proposed element, accurate linear and nonlinear static deformations, as well as realistic dynamic behavior, can be achieved with a smaller computational effort than by using existing shear deformable two-dimensional beam elements.