Finite-strain analysis in orthogneiss of the Gran Paradiso massif, Western Alps, Italy

A finite-strain study in the Gran Paradiso massif of the Italian Western Alps has been carried out to elucidate whether ductile strain shows a relationship to nappe contacts and to shed light on the nature of the subhorizontal foliation typical of the gneiss nappes in the Alps. The Rf/_ and Fry methods used on feldspar porphyroclasts from 143 augengneiss and 11 conglomerate samples of the Gran Paradiso unit (upper tectonic unit of the Gran Paradiso massif), as well as, 9 augengneiss (Erfaulet granite) and 3 quartzite conglomerate samples from the underlying Erfaulet unit (lower unit of the Gran Paradiso massif), and 1 sample from mica schist. Microstructures and thermobarometric data show that feldspar ductility at temperatures >~450°C occurred only during high-pressure metamorphism, when the rocks were underplated beneath the overriding Adriatic plate. Therefore, the finite-strain data can be related to high-pressure metamorphism in the Alpine subduction zone. The augen gneiss was heterogeneously deformed and axial ratios of the strain ellipse in XZ sections range from 2.1 to 69.8. The long axes of the finite-strain ellipsoids trend W/WNW and the short axes are subvertical associated with a subhorizontal foliation. The strain magnitudes do not increase towards the nappe contacts. Geochemical work shows that the accumulation of finite strain was not associated with any significant volume strain. Hence, the data indicate flattening strain type in the Gran Paradiso unit and constrictional strain type in the Erfaulet unit and prove deviations from simple shear. In addition, electron microprobe work was undertaken to determine if the analysed fabrics formed during high-P metamorphism. The chemistry of phengites in the studied samples suggests that deformation and final structural juxtaposition of the Gran Paradiso unit against the Erfaulet took place during high-pressure metamorphism. On the other hand, nappe stacking occurred early during subduction probably by brittle imbrication and that ductile strain was superimposed on and modified the nappe structure during high-pressure underplating in the Alpine subduction zone. The accumulation of ductile strain during underplating was not by simple shear and involved a component of vertical shortening, which caused the subhorizontal foliation in the Gran Paradiso massif. It is concluded that this foliation formed during thrusting of the nappes onto each other suggesting that nappe stacking was associated with vertical shortening. The primary evidence for this interpretation is an attenuated metamorphic section with high-pressure metamorphic rocks of the Gran Paradiso unit juxtaposed against the Erfaulet unit. Therefore, the exhumation during high-pressure metamorphism in the Alpine subduction zone involved a component of vertical shortening, which is responsible for the subhorizontal foliation within the nappes.

A staggered approach for the coupling of Cahn–Hilliard type diffusion and finite strain elasticity

We develop an algorithm and computational implementation for simulation of problems that combine Cahn–Hilliard type diffusion with finite strain elasticity. We have in mind applications such as the electro-chemo- mechanics of lithium ion (Li-ion) batteries. We concentrate on basic computational aspects. A staggered algorithm is pro- posed for the coupled multi-field model. For the diffusion problem, the fourth order differential equation is replaced by a system of second order equations to deal with the issue of the regularity required for the approximation spaces. Low order finite elements are used for discretization in space of the involved fields (displacement, concentration, nonlocal concentration). Three (both 2D and 3D) extensively worked numerical examples show the capabilities of our approach for the representation of (i) phase separation, (ii) the effect of concentration in deformation and stress, (iii) the effect of Electronic supplementary material The online version of this article (doi:10.1007/s00466-015-1235-1) contains supplementary material, which is available to authorized users. B P. Areias pmaa@uevora.pt 1 Department of Physics, University of Évora, Colégio Luís António Verney, Rua Romão Ramalho, 59, 7002-554 Évora, Portugal 2 ICIST, Lisbon, Portugal 3 School of Engineering, Universidad de Cuenca, Av. 12 de Abril s/n. 01-01-168, Cuenca, Ecuador 4 Institute of Structural Mechanics, Bauhaus-University Weimar, Marienstraße 15, 99423 Weimar, Germany strain in concentration, and (iv) lithiation. We analyze con- vergence with respect to spatial and time discretization and found that very good results are achievable using both a stag- gered scheme and approximated strain interpolation.

Ductile fabrics in the zone of active oblique convergence near the Alpine Fault, New Zealand: identifying the neotectonic overprint

The mid-crustal Alpine Schist in central Southern Alps, New Zealand has been exhumed during the past similar to3 m.y. on the hanging wall of the oblique-slip Alpine Fault. These rocks underwent ductile deformation during their passage through the similar to 150-km-wide Pacific-Australia plate boundary zone. Likely to be Cretaceous in age, peak metamorphism predates the largely Pliocene and younger oblique convergence that continues to uplift the Southern Alps today. Late Cenozoic ductile deformation constructively reinforced a pre-existing fabric that was well oriented to accommodate a dextral-transpressive overprint. Quartz microstructures below a recently exhumed brittle-ductile transition zone reflect a late Cenozoic increment of ductile strain that was distributed across deeper levels of the Pacific Plate. Deformation was transpressive, including a dextral-normal shear component that bends and rotates a delaminated panel of Pacific Plate crust onto the oblique footwall ramp of the Alpine Fault. Progressive ductile shear in mylonites at the base of the Pacific Plate overprints earlier fabrics in a dextral-reverse sense, a deformation that accompanies translation of the schists up the Alpine Fault. Ductile shear along that structure affects not only the 12-km-thick section of Alpine mylonites, but is distributed across several kilometres of overlying nonmylonitic rocks. (C) 2001 Elsevier Science Ltd. All rights reserved.

A mathematical framework for finite strain elastoplastic consolidation. Part 1: Balance laws, variational formulation, and linearization

A mathematical formulation for finite strain elasto plastic consolidation of fully saturated soil media is presented. Strong and weak forms of the boundary-value problem are derived using both the material and spatial descriptions. The algorithmic treatment of finite strain elastoplasticity for the solid phase is based on multiplicative decomposition and is coupled with the algorithm for fluid flow via the Kirchhoff pore water pressure. Balance laws are written for the soil-water mixture following the motion of the soil matrix alone. It is shown that the motion of the fluid phase only affects the Jacobian of the solid phase motion, and therefore can be characterized completely by the motion of the soil matrix. Furthermore, it is shown from energy balance consideration that the effective, or intergranular, stress is the appropriate measure of stress for describing the constitutive response of the soil skeleton since it absorbs all the strain energy generated in the saturated soil-water mixture. Finally, it is shown that the mathematical model is amenable to consistent linearization, and that explicit expressions for the consistent tangent operators can be derived for use in numerical solutions such as those based on the finite element method.

Elastoplastic consolidation at finite strain. Part 2: finite element implementation and numerical examples

A mathematical model for finite strain elastoplastic consolidation of fully saturated soil media is implemented into a finite element program. The algorithmic treatment of finite strain elastoplasticity for the solid phase is based on multiplicative decomposition and is coupled with the algorithm for fluid flow via the Kirchhoff pore water pressure. A two-field mixed finite element formulation is employed in which the nodal solid displacements and the nodal pore water pressures are coupled via the linear momentum and mass balance equations. The constitutive model for the solid phase is represented by modified Cam—Clay theory formulated in the Kirchhoff principal stress space, and return mapping is carried out in the strain space defined by the invariants of the elastic logarithmic principal stretches. The constitutive model for fluid flow is represented by a generalized Darcy's law formulated with respect to the current configuration. The finite element model is fully amenable to exact linearization. Numerical examples with and without finite deformation effects are presented to demonstrate the impact of geometric nonlinearity on the predicted responses. The paper concludes with an assessment of the performance of the finite element consolidation model with respect to accuracy and numerical stability.

A finite-strain solid–shell using local Löwdin frames and least-squares strains

A finite-strain solid–shell element is proposed. It is based on least-squares in-plane assumed strains, assumed natural transverse shear and normal strains. The singular value decomposition (SVD) is used to define local (integration-point) orthogonal frames-of-reference solely from the Jacobian matrix. The complete finite-strain formulation is derived and tested. Assumed strains obtained from least-squares fitting are an alternative to the enhanced-assumed-strain (EAS) formulations and, in contrast with these, the result is an element satisfying the Patch test. There are no additional degrees-of-freedom, as it is the case with the enhanced-assumed-strain case, even by means of static condensation. Least-squares fitting produces invariant finite strain elements which are shear-locking free and amenable to be incorporated in large-scale codes. With that goal, we use automatically generated code produced by AceGen and Mathematica. All benchmarks show excellent results, similar to the best available shell and hybrid solid elements with significantly lower computational cost.

Phase-field analysis of finite-strain plates and shells including element subdivision

With the theme of fracture of finite-strain plates and shells based on a phase-field model of crack regularization, we introduce a new staggered algorithm for elastic and elasto-plastic materials. To account for correct fracture behavior in bending, two independent phase-fields are used, corresponding to the lower and upper faces of the shell. This is shown to provide a realistic behavior in bending-dominated problems, here illustrated in classical beam and plate problems. Finite strain behavior for both elastic and elasto-plastic constitutive laws is made compatible with the phase-field model by use of a consistent updated-Lagrangian algorithm. To guarantee sufficient resolution in the definition of the crack paths, a local remeshing algorithm based on the phase- field values at the lower and upper shell faces is introduced. In this local remeshing algorithm, two stages are used: edge-based element subdivision and node repositioning. Five representative numerical examples are shown, consisting of a bi-clamped beam, two versions of a square plate, the Keesecker pressurized cylinder problem, the Hexcan problem and the Muscat-Fenech and Atkins plate. All problems were successfully solved and the proposed solution was found to be robust and efficient.

A finite-strain solid–shell using local Löwdin frames and least-squares strains

A finite-strain solid–shell element is proposed. It is based on least-squares in-plane assumed strains, assumed natural transverse shear and normal strains. The singular value decomposition (SVD) is used to define local (integration-point) orthogonal frames-of- reference solely from the Jacobian matrix. The complete finite-strain formulation is derived and tested. Assumed strains obtained from least-squares fitting are an alternative to the enhanced-assumed-strain (EAS) formulations and, in contrast with these, the result is an element satisfying the Patch test. There are no additional degrees-of-freedom, as it is the case with the enhanced- assumed-strain case, even by means of static condensation. Least-squares fitting produces invariant finite strain elements which are shear-locking free and amenable to be incorporated in large-scale codes. With that goal, we use automatically generated code produced by AceGen and Mathematica. All benchmarks show excellent results, similar to the best available shell and hybrid solid elements with significantly lower computational cost.

Semi-implicit finite strain constitutive integration and mixed strain/stress control based on intermediate configurations

A new semi-implicit stress integration algorithm for finite strain plasticity (compatible with hyperelas- ticity) is introduced. Its most distinctive feature is the use of different parameterizations of equilibrium and reference configurations. Rotation terms (nonlinear trigonometric functions) are integrated explicitly and correspond to a change in the reference configuration. In contrast, relative Green–Lagrange strains (which are quadratic in terms of displacements) represent the equilibrium configuration implicitly. In addition, the adequacy of several objective stress rates in the semi-implicit context is studied. We para- metrize both reference and equilibrium configurations, in contrast with the so-called objective stress integration algorithms which use coinciding configurations. A single constitutive framework provides quantities needed by common discretization schemes. This is computationally convenient and robust, as all elements only need to provide pre-established quantities irrespectively of the constitutive model. In this work, mixed strain/stress control is used, as well as our smoothing algorithm for the complemen- tarity condition. Exceptional time-step robustness is achieved in elasto-plastic problems: often fewer than one-tenth of the typical number of time increments can be used with a quantifiable effect in accuracy. The proposed algorithm is general: all hyperelastic models and all classical elasto-plastic models can be employed. Plane-stress, Shell and 3D examples are used to illustrate the new algorithm. Both isotropic and anisotropic behavior is presented in elasto-plastic and hyperelastic examples.

Finite-strain laminates: Bending-enhanced hexahedron and delamination

With a new finite strain anisotropic framework, we introduce a unified approach for constitutive model- ing and delamination of composites. We describe a finite-strain semi-implicit integration algorithm and the application to assumed-strain hexahedra. In a laminate composite, the laminae are modeled by an anisotropic Kirchhoff/Saint-Venant material and the interfaces are modeled by the exponential cohesive law with intrinsic characteristic length and the criterion by Benzeggagh and Kenane for the equivalent fracture toughness. For the element formulation, a weighted least-squares algorithm is used to calculate the mixed strain. Löwdin frames are used to model orthotropic materials without the added task of per- forming a polar decomposition or empirical frames. To assess the validity of our proposals and inspect step and mesh size dependence, a least-squares based hexahedral element is implemented and tested in depth in both deformation and delamination examples.

Least-squares finite strain hexahedral element/constitutive coupling based on parametrized configurations and the Löwdin frame

Two novelties are introduced: (i) a finite-strain semi-implicit integration algorithm compatible with current element technologies and (ii) the application to assumed-strain hexahedra. The Löwdin algo- rithm is adopted to obtain evolving frames applicable to finite strain anisotropy and a weighted least- squares algorithm is used to determine the mixed strain. Löwdin frames are very convenient to model anisotropic materials. Weighted least-squares circumvent the use of internal degrees-of-freedom. Het- erogeneity of element technologies introduce apparently incompatible constitutive requirements. Assumed-strain and enhanced strain elements can be either formulated in terms of the deformation gradient or the Green–Lagrange strain, many of the high-performance shell formulations are corotational and constitutive constraints (such as incompressibility, plane stress and zero normal stress in shells) also depend on specific element formulations. We propose a unified integration algorithm compatible with possibly all element technologies. To assess its validity, a least-squares based hexahedral element is implemented and tested in depth. Basic linear problems as well as 5 finite-strain examples are inspected for correctness and competitive accuracy.

Folding processes and solitary waves in structural geology

Evolution of localized folding patterns in layered elastic and visco-elastic materials is reviewed in the context of compressed geological systems. The thin strut or plate embedded in a visco-elastic medium is used as an archetypal example to describe localized buckles which, in contrast to those from earlier formulations, appear in the absence of triggering imperfections. Structural and material effects are surveyed and important nonlinear characteristics are identified. A brief review of possible methods of analysis is conducted.