Anisotropic viscous models of large-deformation Mohr-Coulomb failure

We have developed a way to represent Mohr-Coulomb failure within a mantle-convection fluid dynamics code. We use a viscous model of deformation with an orthotropic viscoplasticity (a different viscosity is used for pure shear to that used for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation, neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane, which is formally equivalent in any symmetric formulation. We also allow for strain-weakening of the material. The material can remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-integration-point finite element code and show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion. The formulation itself is general and applies to 2- and 3-dimensional problems.

An anisotropic viscous representation of Mohr-Coulomb failure for use is modelling coupled mantle-continent dynamics

In mantle convection models it has become common to make use of a modified (pressure sensitive, Boussinesq) von Mises yield criterion to limit the maximum stress the lithosphere can support. This approach allows the viscous, cool thermal boundary layer to deform in a relatively plate-like mode even in a fully Eulerian representation. In large-scale models with embedded continental crust where the mobile boundary layer represents the oceanic lithosphere, the von Mises yield criterion for the oceans ensures that the continents experience a realistic broad-scale stress regime. In detailed models of crustal deformation it is, however, more appropriate to choose a Mohr-Coulomb yield criterion based upon the idea that frictional slip occurs on whichever one of many randomly oriented planes happens to be favorably oriented with respect to the stress field. As coupled crust/mantle models become more sophisticated it is important to be able to use whichever failure model is appropriate to a given part of the system. We have therefore developed a way to represent Mohr-Coulomb failure within a code which is suited to mantle convection problems coupled to large-scale crustal deformation. Our approach uses an orthotropic viscous rheology (a different viscosity for pure shear to that for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation - neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane which is formally equivalent in any symmetric formulation. It is also important to allow strain-weakening of the material. The material should remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-Integration-point finite element code and will show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion, and comparisons between mantle convection models using the von Mises versus the Mohr-Coulomb yield criteria. The formulation itself is general and applies to 2D and 3D problems, although it is somewhat more complicated to identify the slip plane in 3D.

Comparison of (120)Sn((6)He, (6)He)(120)Sn and (120)Sn(alpha,alpha)(120)Sn elastic scattering and signatures of the (6)He neutron halo in the optical potential

Cross sections of (120)Sn(alpha,alpha)(120)Sn elastic scattering have been extracted from the alpha-particle-beam contamination of a recent (120)Sn((6)He,(6)He)(120)Sn experiment. Both reactions are analyzed using systematic double-folding potentials in the real part and smoothly varying Woods-Saxon potentials in the imaginary part. The potential extracted from the (120)Sn((6)He,(6)He)(120)Sn data may be used as the basis for the construction of a simple global (6)He optical potential. The comparison of the (6)He and alpha data shows that the halo nature of the (6)He nucleus leads to a clear signature in the reflexion coefficients eta(L) : The relevant angular momenta L with eta(L) >> 0 and eta(L) << 1 are shifted to larger L with a broader distribution. This signature is not present in the alpha-scattering data and can thus be used as a new criterion for the definition of a halo nucleus.

Reaction cross sections for (8)B, (7)Be, and (6)Li+(58)Ni near the Coulomb barrier: Proton-halo effects

Elastic scattering of (8)B, (7)Be, and (6)Li on a (58)Ni target has been measured at energies near the Coulomb barrier. Optical-model fits were made to the experimental angular distributions, and total reaction cross sections were deduced. A comparison with other systems provides striking evidence for proton-halo effects on (8)B reactions. As opposed to the situation for the neutron-halo nucleus (6)He, for which particle transfer dominates, the ""extra"" cross section observed for (8)B appears to result entirely from projectile breakup.

Control of the persistent currents in two interacting quantum rings through the Coulomb interaction and interring tunneling

The persistent current in two vertically coupled quantum rings containing few electrons is studied. We find that the Coulomb interaction between the rings in the absence of tunneling affects the persistent current in each ring and the ground-state configurations. Quantum tunneling between the rings alters significantly the ground state and the persistent current in the system.

A numerical study of flexural buckling of foliated rock slopes

The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter-layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat-type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic-perfectly plastic. Condition of slip at the interfaces are determined by a Mohr-Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14:87-104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub-vertical slopes. Copyright (C) 2001 John Wiley & Sons, Ltd.

Simulation of the micro-physics of rocks using LSMearth

The particle-based Lattice Solid Model (LSM) was developed to provide a basis to study the physics of rocks and the nonlinear dynamics of earthquakes (MORA and PLACE, 1994; PLACE and MORA, 1999). A new modular and flexible LSM approach has been developed that allows different microphysics to be easily included in or removed from the model. The approach provides a virtual laboratory where numerical experiments can easily be set up and all measurable quantities visualised. The proposed approach provides a means to simulate complex phenomena such as fracturing or localisation processes, and enables the effect of different micro-physics on macroscopic behaviour to be studied. The initial 2-D model is extended to allow three-dimensional simulations to be performed and particles of different sizes to be specified. Numerical bi-axial compression experiments under different confining pressure are used to calibrate the model. By tuning the different microscopic parameters (such as coefficient of friction, microscopic strength and distribution of grain sizes), the macroscopic strength of the material and can be adjusted to be in agreement with laboratory experiments, and the orientation of fractures is consistent with the theoretical value predicted based on Mohr-Coulomb diagram. Simulations indicate that 3-D numerical models have different macroscopic properties than in 2-D and, hence, the model must be recalibrated for 3-D simulations. These numerical experiments illustrate that the new approach is capable of simulating typical rock fracture behaviour. The new model provides a basis to investigate nucleation, rupture and slip pulse propagation in complex fault zones without the previous model limitations of a regular low-level surface geometry and being restricted to two-dimensions.

Determination of stresses around a cylindrical single pile caused by horizontal movements of soils with the finitive element method

This study was carried out with the aim of modeling in 2D, in plain strain, the movement of a soft cohesive soil around a pile, in order to enable the determination of stresses resulting along the pile, per unit length. The problem in study fits into the large deformations problem and can be due to landslide, be close of depth excavations, to be near of zones where big loads are applied in the soil, etc. In this study is used an constitutive Elasto-Plastic model with the failure criterion of Mohr-Coulomb to model the soil behavior. The analysis is developed considering the soil in undrained conditions. To the modeling is used the finite element program PLAXIS, which use the Updated Lagrangian - Finite Element Method (UL-FEM). In this work, special attention is given to the soil-pile interaction, where is presented with some detail the formulation of the interface elements and some studies for a better understand of his behavior. It is developed a 2-D model that simulates the effect of depth allowing the study of his influence in the stress distribution around the pile. The results obtained give an important base about how behaves the movement of the soil around a pile, about how work the finite element program PLAXIS and how is the stress distribution around the pile. The analysis demonstrate that the soil-structure interaction modeled with the UL-FEM and interface elements is more appropriate to small deformations problems.

Dimensionamento de muros de suporte considerando diferentes metodologias de cálculo

A presente dissertação consiste em verificar a tendência do coeficiente de segurança quando se varia certos parâmetros (ângulo de atrito interno do terreno, inclinação do terrapleno no tardoz do muro e o ângulo que a massa de solo faz quando se comporta como parte integrante do muro numa situação limite) considerados no cálculo. Para atingir os objectivos anteriormente referidos, dividiu-se o trabalho em duas fases, a primeira fase teve como objectivo verificar qual a tendência do coeficiente de segurança quando sujeito à variação de dois parâmetros, o ângulo de atrito interno do terreno e a inclinação do terrapleno no tardoz do muro que varia de 5° até ao valr do ângulo de atrito interno do terreno. A segunda fase consiste em analisar qual a tendência do coeficiente de segurança quando sujeito à variação de três parâmetros, o ângulo de atrito interno do terreno, de 20° a 45°, a inclinação do terrapleno no tardoz do muro, de 10° até ao valor do ângulo de atrito interno do terreno, e o ângulo que a massa de solo faz quando se comporta como parte integrante do muro numa situação limite. Para ambas as situações efectuaram-se os respectivos cálculos quer pela Teoria de Rankine e quer pela Teoria de Mohr – Coulomb, havendo casos em que foi necessário conjugar estas duas Teorias.

Fractional describing function of systems with Coulomb friction

This paper studies the describing function (DF) of systems constituted by a mass subjected to nonlinear friction. The friction force is decomposed into two components, namely, the viscous and the Coulomb friction. The system dynamics is analyzed in the DF perspective revealing a fractional-order behavior. The reliability of the DF method is evaluated through the signal harmonic contents.

Cargas de colapso de sistemas de alvenaria

Dissertação apresentada na Faculdade de Ciências e Tecnologia da Universidade Nova de Lisboa para obtenção do grau de Mestre em Engenharia Civil na especialidade de Estruturas

Análise limite aplicada ao estudo tridimensional de paredes de alvenaria

Dissertação para obtenção do Grau de Mestre em Engenharia Civi

Aplicação de retroanálise em obras subterrâneas

Dissertação para obtenção do Grau de Mestre em Engenharia Civil - Perfil Geotecnia

Avaliação da capacidade resistente última de estacas em compressão usando análises em tensão/deformação

Dissertação para obtenção do Grau de Mestre em Engenharia Civil - Ramo de Estruturas e Geotecnia

Comportamento e modelação da ligação GFRP/betão em elementos de betão armado expostos a ambientes agressivos

Dissertação para a obtenção do Grau de Doutor em Engenharia Civil, na Especialidade de Estruturas pela Universidade Nova de Lisboa