Effect of the Flow Rule on the Bearing Capacity of Strip Foundations on Sand by the Upper-Bound Limit Analysis and Slip Lines

The influence of the flow rule on the bearing capacity of strip foundations placed on sand was investigated using a new kinematic approach of upper-bound limit analysis. The method of stress characteristics was first used to find the mechanism of the failure and to compute the stress field by using the Mohr-Coulomb yield criterion. Once the failure mechanism had been established, the kinematics of the plastic deformation was established, based on the requirements of the upper-bound limit theorem. Both associated and nonassociated plastic flows were considered, and the bearing capacity was obtained by equating the rate of external plastic work to the rate of the internal energy dissipation for both smooth and rough base foundations. The results obtained from the analysis were compared with those available from the literature. (C) 2014 American Society of Civil Engineers.

Lower-Bound Axisymmetric Formulation for Geomechanics Problems Using Nonlinear Optimization

A lower-bound limit analysis formulation, by using two-dimensional finite elements, the three-dimensional Mohr-Coulomb yield criterion, and nonlinear optimization, has been given to deal with an axisymmetric geomechanics stability problem. The optimization was performed using an interior point method based on the logarithmic barrier function. The yield surface was smoothened (1) by removing the tip singularity at the apex of the pyramid in the meridian plane and (2) by eliminating the stress discontinuities at the corners of the yield hexagon in the pi-plane. The circumferential stress (sigma(theta)) need not be assumed. With the proposed methodology, for a circular footing, the bearing-capacity factors N-c, N-q, and N-gamma for different values of phi have been computed. For phi = 0, the variation of N-c with changes in the factor m, which accounts for a linear increase of cohesion with depth, has been evaluated. Failure patterns for a few cases have also been drawn. The results from the formulation provide a good match with the solutions available from the literature. (C) 2014 American Society of Civil Engineers.

Intrinsic correlation between dilatation and pressure sensitivity of plastic flow in metallic glasses

Taking shear-induced dilatation into consideration in shear transformation zone (STZ) operations, we derive a new yield criterion that reflects the pressure sensitivity in plastic flow in metallic glasses (MGs), which agrees well with experiments. Furthermore, an intrinsic theoretical correlation between the pressure sensitivity coefficient and the dilatation factor is revealed. It is found that the pressure sensitivity of plastic flow of MGs originates in the dilatation of microscale STZs.

Avaliação de ligações soldadas de perfis tubulares em estruturas de aço através do método dos elementos finitos.

A necessidade de se obter soluções para os diversos problemas estruturais na engenharia, associado ao amplo crescimento e difusão da construção metálica fez surgir, por volta dos anos 60, os perfis tubulares. As ligações soldadas entre perfis tubulares são as ligações de maior difusão nas construções com este tipo de perfil, sendo necessário portanto, estudos mais detalhados que tornem a utilização de ligações soldadas mais seguras e otimizadas. O Eurocode 3 possui um item específico para as ligações entre perfis tubulares, em especial as ligações soldadas devido a suas particularidades. Infelizmente, a norma brasileira NBR8800 não aborda as ligações envolvendo perfis tubulares. Para a execução da análise numérica no presente trabalho, faz-se necessário uma revisão bibliográfica. Com estes resultados, foram feitas modelagens de ligações tubulares com o método dos elementos finitos de forma a otimizar os modelos a serem utilizados em uma análise paramétrica futura. Desta forma, foram desenvolvidos dois modelos numéricos, um considerando ligação tipo T entre perfis tubulares quadrados e outro para uma ligação tipo K entre perfis circulares. Estes modelos foram caracterizados através de elementos de casca com seis graus de liberdade por nó considerando-se adicionalmente o efeito de membrana. A análise não-linear realizada considerou a não-linearidade do material através do critério de plastificação de Von Mises através de uma lei constitutiva tensão versus deformação bi-linear e a não-linearidade geométrica foi atribuída através da formulação de Lagrange atualizado. Dentre as principais conclusões obtidas no presente trabalho, pode-se citar que os resultados para as ligações tipo T, o Eurocode 3, fornece resultados que precisam ser observados com cautela. Todavia, para as ligações do tipo K, os resultados numéricos mostraram-se sempre inferiores aos valores através do Eurocode 3, representando um dimensionamento a favor da segurança.

Optimal displacement mechanisms beneath shallow foundations on linear-elastic perfectly plastic soil

An energy method for a linear-elastic perfectly plastic method utilising the von Mises yield criterion with associated flow developed in 2013 by McMahon and co-workers is used to compare the ellipsoidal cavity-expansion mechanism, from the same work, and the displacement fields of other research by Levin, in 1995, and Osman and Bolton, in 2005, which utilise the Hill and Prandtl mechanisms respectively. The energy method was also used with a mechanism produced by performing a linear-elastic finite-element analysis in Abaqus. At small values of settlement and soil rigidity the elastic mechanism provides the lowest upper-bound solution, and matches well with finite-element analysis results published in the literature. At typical footing working loads and settlements the cavity-expansion mechanism produces a more optimal solution than the displacement fields within the Hill and Prandtl mechanisms, and also matches well with the published finite-element analysis results in this range. Beyond these loads, at greater footing settlements, or soil rigidity, the Prandtl mechanism is shown to be the most appropriate.

Generalized Viscoplastic Modeling of Debris Flow

Various concepts have been proposed or used in the development of rheological models for debris flow. The earliest model developed by Bagnold was based on the concept of the “dispersive” pressure generated by grain collisions. Bagnold’s concept appears to be theoretically sound, but his empirical model has been found to be inconsistent with most theoretical models developed from non-Newtonian fluid mechanics. Although the generality of Bagnold’s model is still at issue, debris-flow modelers in Japan have generally accepted Takahashi’s formulas derived from Bagnold’s model. Some efforts have recently been made by theoreticians in non-Newtonian fluid mechanics to modify or improve Bagnold’s concept or model. A viable rheological model should consist both of a rate-independent part and a rate-dependent part. A generalized viscoplastic fluid (GVF) model that has both parts as well as two major rheological properties (i.e., the normal stress effect and soil yield criterion) is shown to be sufficiently accurate, yet practical, for general use in debris-flow modeling. In fact, Bagnold’s model is found to be only a particular case of the GVF model. Analytical solutions for (steady) uniform debris flows in wide channels are obtained from the GVF model based on Bagnold’s simplified assumption of constant grain concentration.

小湾坝基卸荷岩体参数研究

The Xiao-wan Power Station is the second highest arch dam in the world under construction. The height of the dam is about 292m. Large-scale excavation in the dam foundation of Xiao-wan Power Station has brought intensive unloading phenomenon. We collected a large number of firsthand data on unloaded rock mass in dam foundation, which supplies a natural testing ground for researching unloaded rock mass after excavation. Detailed study was carried out on the parameters of unloaded rock mass in the dam foundation of Xiao-wan Power Station. The study is not only importance for the Xiao-wan Power Station, but also has important instruction significance to similar projects in the Southwest of China. In order to study the mechanical parameters of unloaded rock mass, large field and laboratory tests were carried out. The test results showed the size effect of the sample is obvious. The change of deformation modulus of rock is not obvious. However, the Poisson's ratio of rock is increased under unloaded condition, its value is even more than 0.5. The theoretical forecasted results is accordance to the field tests including sound wave data and deformation monitoring data, which shows the forecasted results were reasonable. The soften yield criterion was adopted in the thesis to study the characteristics of the brittle rock mass in order to simulate their brittle failures. Based on the study results above, the transform of the structural plane network model to the numerical one was carried out, which made it feasible to consider the influence of large amount joints on the mechanical characteristic of rock mass in the numerical analyses. Using a factor, the degree of the damage or strengthen of rock mass can be determined rapidly, which proposed a rapid and feasible method for the determination of the parameters of rock mass.

Modeling Coulombic failure of sea ice with leads

[1] Sea ice failure under low-confinement compression is modeled with a linear Coulombic criterion that can describe either fractural failure or frictional granular yield along slip lines. To study the effect of anisotropy we consider a simplified anisotropic sea ice model where the sea ice thickness depends on orientation. Accommodation of arbitrary deformation requires failure along at least two intersecting slip lines, which are determined by finding two maxima of the yield criterion. Due to the anisotropy these slip lines generally differ from the standard, Coulombic slip lines that are symmetrically positioned around the compression direction, and therefore different tractions along these slip lines give rise to a nonsymmetric stress tensor. We assume that the skewsymmetric part of this tensor is counterbalanced by an additional elastic stress in the sea ice field that suppresses floe spin. We consider the case of two leads initially formed in an isotropic ice cover under compression, and address the question of whether these leads will remain active or new slip lines will form under a rotation of the principal compression direction. Decoupled and coupled models of leads are considered and it is shown that for this particular case they both predict lead reactivation in almost the same way. The coupled model must, however, be used in determining the stress as the decoupled model does not resolve the stress asymmetry properly when failure occurs in one lead and at a new slip line.

Análise não linear física e geométrica de pórticos espaciais

Apresenta-se uma formulação do tipo incrementaliterativa destinada a análise não linear de pórticos espaciais. Considera-se os efeitos não lineares introduzidos pelas mudanças de configuração geométrica da estrutura e também pela combinação destes efeitos com aqueles inerentes ao comportamento plástico exibido pelo material. As relações cinemáticas empregadas permitem a consideração de deslocamentos arbitrariamente grandes, acompanhadas de pequenas deformações . A modelagem do comportamento plástico do material é efetuada através do conceito de rótula plástica, estabelecido a partir de um critério de plastificação generalizado. Adota-se uma matriz de rigidez geométrica de barra baseada em momentos semitangenciais. Para elementos com extremos plastificados, é deduzida uma matriz de rigidez elasto-plástica. Emprega-se um método numérico do tipo incremental-iterativo, que utiliza como condição básica de controle da análise a constância do trabalho realizado pelos incrementos de cargas, em cada passo incremental (Método de Controle por Trabalho).A formulação permite uma descricão completa do desempenho mecânico da estrutura, inclusive em estágio de deformação pós-crítico em que ocorre regressão do carregamento com aumento de deslocamentos, ou vice-versa. A formulação foi implementada em um programa computacional elaborado em linguagem FORTRAN. Vários exemplos numéricos são apresentados para mostrar a eficiência das procedimentos propostos.

Análise elastoplástica de grelhas planas de concreto armado

O objetivo deste trabalho é apresentar um programa computacional para analisar problemas elastoplásticos com pequenas deformações em grelhas planas de concreto armado, utilizando o método dos deslocamentos. O material é definido através de um critério de escoamento generalizado que leva em consideração o efeito dos esforços de corte, momento fletor e momento torsor nas seções críticas e apresenta endurecimento isótropo. Uma solução incremental empregando o método da matriz de rigidez tangente foi adotado. Algumas comparações foram feitas com resultados experimentais, numéricos e teóricos.

Avaliação de um modelo matemático para a determinação do coeficiente de atrito no processo de conformação por dobramento

Pós-graduação em Engenharia Mecânica - FEG

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

Avaliação do processamento por atrito linear em chapas da liga de titânio Ti-6Al-4V.

Esta tese tem por objetivo a aplicação do processamento por atrito linear na liga de titânio Ti-6Al-4V. Derivado da solda por atrito linear, é um processo recente desenvolvido na década de 90 para união de alumínio. Sua aplicação em outros tipos de materiais como aços e ligas de alto desempenho, em especial o titânio, tem interessado a industria. A metodologia utilizada nesta tese para avaliar o processamento por atrito linear, consistiu na execução de ensaios mecânicos de tração em condições mistas em chapas da liga de titânio Ti-6Al-4V. A máquina utilizada para o processamento das chapas foi um centro de usinagem CNC convencional, adaptado com dispositivos especiais. Além dos ensaios de tração em condições mistas, foram executadas medições de microdurezas nas regiões atingidas pelo processo, avaliação das microestruturas resultantes e medições de tensão residual para uma caracterização mais ampla do processo. As microestruturas na região processada são caracterizadas por uma estrutura totalmente transformada. As temperaturas de pico na região processada excederam a temperatura -transus durante o processamento e a transformação da fase + ocorreu durante a fase de resfriamento. A transformação da fase para resultou na formação de agulhas de fase nos contornos e pelo interior dos grãos da fase . Pequenas regiões com estrutura equiaxial de grãos ( globular) foram observados na zona de processamento. A abordagem dos resultados quantitativos foi feita de forma estatística, visando identificar os parâmetros de maior interação com os resultados observados. Foi identificado nesta tese que a rotação da ferramenta apresentou a maior influência nos resultados de tensão residual, microdureza e tensão de escoamento. Uma importante contribuição à modelagem da tensão de escoamento para materiais anisotrópicos é proposta, baseado em um critério de escoamento ortotrópico. Equações complementares baseadas nos testes mistos de tração e cisalhamento são propostas para modificar o modelo ortotrópico. O intuito deste modelo é indicar em que condições o material tem seu regime de escoamento atingido, podendo servir de base para simulações práticas de peças em condições similares.

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.