The shear strength of soils containing undulating shear zones—a numerical study

1988 ◽  
Vol 25 (3) ◽  
pp. 550-558 ◽  
Author(s):  
George T. Dounias ◽  
David M. Potts ◽  
Peter R. Vaughan

This paper investigates the behaviour of a clay layer containing an undulating shear surface, when sheared across the undulations. A relatively long soil block containing an undulating weak zone of finite thickness is assumed. A finite element study is undertaken, examining the effect of the thickness and the amplitude of the weak zone on the overall strength and deformation of the block. Also examined is the behaviour of the block when either only the weak zone or both the weak zone and the intact soil are strain softening. Key words: undulating shear surface, finite elements, strain softening.

2021 ◽  
Author(s):  
Pritam Ghosh ◽  
Kathakali Bhattacharyya

<p>We examine how the deformation profile and kinematic evolutionary paths of two major shear zones with prolonged deformation history and large translations differ with varying structural positions along its transport direction in an orogenic wedge. We conduct this analysis on multiple exposures of the internal thrusts from the Sikkim Himalayan fold thrust belt, the Pelling-Munsiari thrust (PT), the roof thrust of the Lesser Himalayan duplex (LHD), and the overlying Main Central thrust (MCT). These two thrusts are regionally folded due to growth of the LHD and are exposed at different structural positions. The hinterlandmost exposures of the MCT and PT zones lie in the trailing parts of the duplex, while the foreland-most exposures of the same studied shear zones lie in the leading part of the duplex, and thus have recorded a greater connectivity with the duplex. The thicknesses of the shear zones progressively decrease toward the leading edge indicating variation in deformation conditions. Thickness-displacement plot reveals strain-softening from all the five studied MCT and the PT mylonite zones. However, the strain-softening mechanisms varied along its transport direction with the hinterland exposures recording dominantly dislocation-creep, while dissolution-creep and reaction-softening are dominant in the forelandmost exposures. Based on overburden estimation, the loss of overburden on the MCT and the PT zones is more in the leading edge (~26km and ~15km, respectively) than in the trailing edge (~10km and ~17km, respectively), during progressive deformation. Based on recalibrated recrystallized quartz grain thermometer (Law, 2014), the estimated deformation temperatures in the trailing edge are higher (~450-650°C) than in the leading edge (350-550°C) of the shear zones. This variation in the deformation conditions is also reflected in the shallow-crustal deformation structures with higher fracture intensity and lower spacing in the leading edge exposures of the shear zones as compared to the trailing edge exposures.</p><p>The proportion of mylonitic domains and micaceous minerals within the exposed shear zones increase and grain-size of the constituent minerals decreases progressively along the transport direction. This is also consistent with progressive increase in mean R<sub>s</sub>-values toward leading edge exposures of the same shear zones. Additionally, the α-value (stretch ratio) gradually increases toward the foreland-most exposures along with increasing angular shear strain. Vorticity estimates from multiple incremental strain markers indicate that the MCT and PT zones generally record a decelerating strain path. Therefore, the results from this study are counterintuitive to the general observation of a direct relationship between higher Rs-value and higher pure-shear component. We explain this observation in the context of the larger kinematics of the orogen, where the leading edge exposures have passed through the duplex structure, recording the greatest connectivity and most complete deformation history, resulting in the weakest shear zone that is also reflected in the deformation profiles and strain attributes. This study demonstrates that the same shear zone records varying deformation profile, strain and kinematic evolutionary paths due to varying deformation conditions and varying connectivity to the underlying footwall structures during progressive deformation of an orogenic wedge.</p>


2021 ◽  
Author(s):  
Nikola Stanković ◽  
Vesna Cvetkov ◽  
Vladica Cvetković

<p>We report updated results of our ongoing research on constraining geodynamic conditions associated with the final closure of the Vardar branch of the Tethys Ocean by means of application of numerical simulations (previous interim results reported in EGU2020-5919).</p><p>The aim of our numerical study is to test the hypothesis that a single eastward subduction in the Jurassic is a valid explanation for the occurrence of three major, presently observed geological entities that are left behind after the closure of the Vardar Tethys. These include: ophiolite-like igneous rocks of the Sava-Vardar zone and presumably subduction related Timok Magmatic Complex, both Late Cretaceous in age as well as Jurassic ophiolites obducted onto the Adriatic margin. In our simulations we initiate an intraoceanic subduction in the Early/Middle Jurassic, which eventually transitions into an oceanic closure and subsequent continental collision processes.</p><p>In the scope of our study numerical simulations are performed by solving a set of partial differential equations: the continuity equation, the Navier-Stokes equations and the temperature equation. To this end we used I2VIS thermo-mechanical code which utilizes marker in cell approach with finite difference discretization of equations on a staggered grid [Gerya et al., 2000; Gerya&Yuen, 2003].</p><p>The 2D model consists of two continental plates separated by two oceanic slabs connected at a mid-oceanic ridge. Intraoceanic subduction is initiated along the ridge by assigning a weak zone beneath the ridge. Time-dependent boundary conditions for velocity are imposed on the simulation in order to model a transient spreading period. The change of sign in plate velocities is found to be useful for both obtaining obduction / ophiolite emplacement [Duretz et al., 2016] and causing back-arc extension. Changes in velocities are linear in time. Simulations follow a three-phase evolution of velocity boundary conditions consisting of two convergent phases separated by a single divergent phase where spreading regime is dominant. Effect of duration and magnitude of the second phase on model evolution is also explored.</p><p>Our so far obtained simulations were able to reproduce the westward obduction and certain extension processes along the active (European) margin, which match the existing geological relationships. However, the simulations involve an unreasonably short geodynamic event (cca 15-20 My) and we are working on solving this problem with new simulations. </p>


2020 ◽  
Vol 6 (22) ◽  
pp. eaba1529
Author(s):  
Christopher J. Tulley ◽  
Åke Fagereng ◽  
Kohtaro Ujiie

The rheology of the metamorphosed oceanic crust may be a critical control on megathrust strength and deformation style. However, little is known about the strength and deformation style of metamorphosed basalt. Exhumed megathrust shear zones exposed on Kyushu, SW Japan, contain hydrous metabasalts deformed at temperatures between ~300° and ~500°C, spanning the inferred temperature-controlled seismic-aseismic transition. Field and microstructural observations of these shear zones, combined with quartz grain-size piezometry, indicate that metabasalts creep at shear stresses <100 MPa at ~370°C and at shear stresses <30 MPa at ~500°C. These values are much lower than those suggested by viscous flow laws for basalt. The implication is that relatively weak, hydrous, metamorphosed oceanic crust can creep at low viscosities over a wide shear zone and have a critical influence on plate interface strength and deformation style around the seismic-aseismic transition.


1999 ◽  
Vol 10 (2) ◽  
pp. 151-159
Author(s):  
K. B. Murashkin ◽  
V. M. Paskonov ◽  
S. V. Fortova

1977 ◽  
Vol 14 (11) ◽  
pp. 2510-2523 ◽  
Author(s):  
P. R. Cobbold

This paper offers a generalized mechanical explanation for the origin and development of bandlike deformation structures such as shear zones, mylonite zones, kink bands, 'pressure-solution' seams, extension gashes, and similar folds.Methods of continuum mechanics are used to examine permissible variations in strain rate, stress, and rheological properties across a region containing ideal banded perturbations. For bands to develop, the rheological properties must vary across the banding. The physical basis for this variation is a corresponding variation in microstructure or chemical composition, influenced in turn by finite deformation, stress, and temperature. Many rocks are likely to soften or harden during progressive deformation and these changes may be enhanced by thermal or other agents. Deformation softening (including strain softening and rotation softening) is a cause of instability and has two effects: first, the deformation tends to accelerate under constant stress; second, the deformation tends to become locally perturbed. Deformation hardening has compensatory effects.Banded perturbations do not appear spontaneously in a deforming rock, but evolve towards an ideal banded form by processes of nucleation and propagation. Evidence for these processes comes from theoretical analysis, experimental data, and observation of bandlike structures that have formed as a result of natural processes of deformation.


2006 ◽  
Vol 129 (6) ◽  
pp. 717-726 ◽  
Author(s):  
Kamil Kahveci

This numerical study looks at laminar natural convection in an enclosure divided by a partition with a finite thickness and conductivity. The enclosure is assumed to be heated using a uniform heat flux on a vertical wall, and cooled to a constant temperature on the opposite wall. The governing equations in the vorticity-stream function formulation are solved by employing a polynomial-based differential quadrature method. The results show that the presence of a vertical partition has a considerable effect on the circulation intensity, and therefore, the heat transfer characteristics across the enclosure. The average Nusselt number decreases with an increase of the distance between the hot wall and the partition. With a decrease in the thermal resistance of the partition, the average Nusselt number shows an increasing trend and a peak point is detected. If the thermal resistance of the partition further declines, the average Nusselt number begins to decrease asymptotically to a constant value. The partition thickness has little effect on the average Nusselt number.


Author(s):  
Rafiq Manna ◽  
Patrick H. Oosthuizen

Abstract Simultaneous natural convective heat transfer from the top, bottom and side surfaces of two-sided inclined square plates having various thicknesses has been numerically investigated. The aim of this work is to determine whether the plate thickness has a significant influence on the heat transfer rates from the plate surfaces when the plate is inclined to the horizontal and to determine how the heat transfer rate varies with this angle of inclination. The upper, lower and side surfaces of the plate have been assumed to be isothermal and at the same temperature which is higher than that of the surrounding fluid. The range of conditions considered is such that laminar, transitional, and turbulent flow occur over the plate. The numerical solution has been obtained using the commercial CFD solver ANSYS FLUENT©. In this study, results have only been obtained for the case where the plate is exposed to air. Inclination angles of between 0 and 40 degrees from the horizontal and plate dimensionless thicknesses (thickness-to-side length ratios) of between 0 and 0.3 have been considered. Variations of the mean Nusselt number with Rayleigh number for the top surface, bottom surface, side surface and that averaged over all heated surfaces of the plate for various inclination angles and for various plate dimensionless thicknesses have been obtained.


1987 ◽  
Vol 24 (3) ◽  
pp. 430-440 ◽  
Author(s):  
D. H. Chan ◽  
N. R. Morgenstern

A finite element analysis was performed of the deformations observed during the excavation of the Edmonton Convention Centre. Local geology in the Edmonton area consists of layers of shale with weak bentonite seams overlain by glacial deposits. The presence of the bentonite seams, which possess strain-softening characteristics, controls the stability of the excavation, which is located in a valley wall. To simulate the influence of the bentonite seams a strain-softening finite element model is used to estimate the amount of deformation in the foundation of the excavation. Field measurements indicated that localized progressive straining had occurred during the excavation process, and substantial heave of the foundation floor was observed. The finite element results show progressive deformation of the excavation and propagation of shear zones. Good agreement between the finite element results and the field observations is obtained. Key words: progressive failure, strain softening, finite element analysis, shear band, excavation stability.


2019 ◽  
Vol 194 ◽  
pp. 322-343 ◽  
Author(s):  
Hussein M. Elsanadedy ◽  
Yousef A. Al-Salloum ◽  
Tarek H. Almusallam ◽  
Abdulhafiz O. Alshenawy ◽  
Husain Abbas

Sign in / Sign up

Export Citation Format

Share Document