[pt] Materiais não Newtonianos reagem de maneira diferente quando expostos a esforços cisalhantes ou extensionais. Uma equação constitutiva em que a tensão é função não só da taxa de deformação, mas também do tipo de escoamento é proposta e analisada teoricamente. Ela combina as informações obtidas em cisalhamento, extensão e movimento de corpo rígido em todas as regiões de um escoamento complexo. Esta equação é testada em códigos numéricos para duas geométricas axissimétricas diferentes: a contração abrupta 4:1 e o canal convergente R (2) z = C. Ambas as geometrias têm sido usadas para efetuar medições da viscosidade extensional. Para validação do modelo, foi construída uma seção de testes com uma contração abrupta 4:1. Foram feitos experimentos com fluidos Newtonianos e viscoelásticos utilizando a técnica de Velocimetria por Imagens de Partículas para determinação do campo de velocidade. / [en] Non-Newtonian materials exhibit different behavior if submitted to shear or extension. A constitutive equation in wich the stress is a function not only of the rate of deformation, but also of the type of the flow is proposed and analyzed theorecticaly. It combines information obtained in shear, extension and rigid body motion in all regions of a complex flow. This equation is tested with numerical codes for two different axissymmetrical geometries, namely, the 4:1 abrupt contraction and the R(2)z=C convergent channel. Both geometries have been employed to measure the extensional viscosity. In order to validate the model, a test section was buit with a 4:1 contraction. Newtonian and viscoelastic fluids were tested using the Particle Image Velocimetry technique to determine the velocity field.

[pt] Este estudo tem como objetivo avaliar a influência do desenvolvimento das propriedades mecânicas da pasta de cimento na integridade do poço de petróleo. Um modelo numérico foi desenvolvido para simular as diferentes fases de construção de poços, tais como perfuração, ação da fluência, reconstrução do revestimento e a hidratação do cimento. O comportamento da fluência do sal foi modelado através uma metodologia que combina o Duplo Mecanismo e a Lei da Potência. Duas fases foram consideradas para a pasta de cimento. Em primeiro lugar, a pasta de cimento foi idealizada como um fluido com um campo de tensão hidrostática. Durante as primeiras cinco horas após a cimentação, uma diminuição do campo de tensão hidrostática foi utilizada a fim de simular a retração autógena e a perda de água no processo. Quando a pasta de cimento foi considerada como elástico com um módulo de rigidez crescente ao longo do tempo quando solidificada. Na mudança de fase do fluido para sólido, um cuidado especial foi dado aos elementos que representam a pasta de cimento hidratada a fim de evitar a sobreposição de malhas. A nova lei de fluência do sal apresentou resultados precisos quando comparado aos valores obtidos em campo e a outros dados encontrados na literatura. No processo de hidratação, a redução das dimensões anulares foi observada. A taxa de deformação diminuiu com o aumento do módulo de rigidez da pasta de cimento. Mais importante ainda, esta deformação foi insignificante quando comparada com as dimensões do poço. Assim, a solidificação do cimento pode ser desconsiderada em simulações expeditas. / [en] The quality of the cement sheath plays a crucial role in the integrity of the oil well. This study aims at assessing the influence of cement hardening on deformations and stresses to which the sheath is submitted. A numerical model was developed to simulate the different stages of well construction such as drilling, creep action, casing reconstruction and - the focus herein - cement hardening. The salt creep behavior was modeled through a new methodology that combines the Double Mechanism with the Power Law. Two stages were considered for the cement. First, the cement was idealized as a fluid with a hydrostatic stress field. During the first five hours after cementation, a decrease in the hydrostatic stress field was employed in order to simulate the chemical shrinkage and loss of water in the hardening process. When hardened, the cement was considered to behave elastically with an increasing stiffness modulus over time. In the phase change from fluid to solid, special care was given to fitting the cement finite element mesh to the annulus, avoiding mesh overlapping. The new salt creep law showed accurate results when compared to field and other creep data. In the hardening process, the annular dimension reduction was expected and noticed. The deformation rate decreased with the increasing cement stiffness modulus. Most importantly, this deformation was negligible when compared with the dimensions of the well. Moreover, for practical and quick simulations, the hardening step can be avoided.

[pt] A performance de protocolos geralmente empregados na caracterização reológica de pastas de cimento é investigada. Para este fim duas pastas de cimento foram utilizadas, Pasta A e Pasta B, com a mesma densidade e razão água-cimento, diferindo quanto aos aditivos empregados. O procedimento API (American Petroleum Institute), tal qual, sua modificação é estudada. Curvas de escoamento são determinadas com o auxílio de um reômetro rotacional através de três diferentes métodos (ramp up/ramp down, ramp down e viscosidade mínima). O procedimento API gerou curvas de escoamento que superestimaram (Pasta A) ou subestimaram (Pasta B) a viscosidade dentro da faixa de cisalhamento estudada. Os desvios tornam-se particulamente maiores em baixas taxas, e para Pasta B, por causa de efeitos tixotrópicos, o regime permanente não é atingido. A modificação do protocolo API investigada resultou em curvas de escoamento similares ao método de viscosidade mínima. Um novo modelo constitutivo fenomenológico irreversível é proposto para descrição reológica de pastas de cimento (voltadas para poços de petróleo) sob escoamento de cisalhamento simples. Para esse objetivo uma série de experimentos foram realizados, como taxa de deformação constante, construção e destruição da microestrutura. O modelo proposto utiliza a própria fluidez, recíproca da viscosidade, para descrever a evolução da estrutura da pasta. O efeito macroscópico do fenômeno de hidratação na reologia é introduzido por meio de uma função exponencial para tensão limite de escoamento. Todos os dez parâmetros necessários no modelo são experimentalmente obtíveis. Os resultados mostraram boa concordância entre os dados experimentais e a modelagem, embora pequenas diferenças tenham sido observadas para taxas de deformação abaixo de 7,1 s(−1). / [en] The performance of reological protocols usually employed for the rheological characterization of oil well cement slurries is investigated. To this end two cement slurries were employed, namely Paste A and Paste B. The API procedure is addressed and a modification to it is investigated. Flow curves are determined with the aid of a rotational rheometer using three different methods (ramp up/down, ramp down, and minimum-viscosity). The API procedure yielded flow curves that significantly overestimate (Paste A) or underestimate (Paste B) the viscosity throughout the whole range of shear rate. The deviations become particularly dramatic in the low end of this range and for Paste B, clearly because the steady state is not attained due to a thixotropic behavior. The investigated modification to the API procedure resulted in flow curves similar to the ones based on the minimum-viscosity method. A novel irrevesible phenomenological constitutive model is proposed to describe oil well cement slurry rheology under simple shear flow. To this end several rheological experiments were conducted, namely constant shear rate, construction and microstructural destruction experiment. The proposed model employ the fluidity itself, reciprocal of viscosity, to describe the paste s structure evolution. The macroscopical rheological effect of hydration phenomena is introduced in the model by means of a yield stress exponential function. All ten parameters needed in the model are experimentally obtainable. The present results have shown a good agreement between experimental data and modeling, albeit minor differences were observed for low shear rates below 7.1 s(−1).

[pt] Problemas de engenharia geotécnica onde a resistência ao cisalhamento é fator determinante para o projeto são geralmente avaliados com base no critério de ruptura de Mohr-Coulomb; no entanto, modelos constitutivos avançados são necessários para abordar grande parte de problemas geotécnicos onde estimativas de deformação são necessárias. Lade e Kim propuseram um modelo constitutivo elastoplástico avançado para investigar o comportamento tensão vs. deformação de materiais com atrito, sob condição drenada ou não drenada. Nesta dissertação, este modelo foi implementado no programa computacional de diferenças finitas FLAC3D com o objetivo de estudar o comportamento elasto-plástico de aterros. As equações governantes do modelo constitutivo de Lade-Kim são apresentadas e um algoritmo explícito de integração de tensões foi implementado na presente pesquisa. Os resultados numéricos obtidos com o modelo de Lade e Kim foram comparados com aqueles publicados na literatura, envolvendo ensaios triaxiais de laboratório realizados em areia e análises numéricas 2D do comportamento de barragem. Finalmente, após validação da implementação do modelo, análises do comportamento 2D e 3D de aterros foram executadas, verificando a influência destes tipos de análise nos valores calculados, bem como resultados numéricos foram também comparados com aqueles determinados pelo tradicional modelo elasto-plástico de Mohr-Coulomb. / [en] Geotechnical engineering problems, where the design is mainly influenced by the shear strength, can be successfully analyzed considering a traditional constitutive model; however, advanced constitutive models must be used when the estimative of deformation is an important factor to be considered. Lade and Kim proposed the Lade-Kim model in order to model a granular material stress-strain constitutive relation under different combinations of effective stress and drained or undrained conditions. In this dissertation, this model is implemented and validated in the FLAC3D finite volume program, in order to study the elasto-plastic behavior of embankments. The governing equations of the constitutive model of Lade-Kim are presented and, using an appropriate stress integration algorithm, it is implemented in FLAC3D. This procedure is validated by using experimental results of triaxial tests performed on the sand of the Sacramento River and performing two-dimensional deformation analyzes of embankments and comparing the results with the literature. Finally, the Lade-Kim model is used in two-dimensional and three-dimensional deformation analyses of embankments in order to analyze the influence of the three-dimensional condition. On the other hand, the same analyzes were carried out using the Mohr- Coulomb constitutive model in order to compare the results and analyze the efficiency of the Lade-Kim constitutive model. Thus, it is concluded that the Lade-Kim constitutive model was successfully implemented and that it provides accurate results regarding modeling geomechanical behavior under different combinations of effective stresses.

[pt] Nos últimos 40 anos, tem-se observado um crescente interesse por compósitos formados pela adição de fibras a matrizes de concreto, os chamados concretos reforçados com fibras. Esse interesse é justificado, sobretudo, pelo significativo ganho em tenacidade que as fibras proporcionam, atenuando as características frágeis do concreto. De fato, em virtude do mecanismo de reforço promovido pelas fibras, o concreto com fibras é capaz de absorver muito mais energia de deformação até a ruptura, apresentando, no regime pós-fissuração, um comportamento muito mais suave que o concreto simples. Esse comportamento é acompanhado por um processo de fissuração mais uniforme, no qual observam-se fissuras mais finas e menos espaçadas. Além disso, registram-se aumentos nos valores de resistência do material e nos níveis de deformação que ele atinge até seu completo esgotamento. Tendo em vista os benefícios que as fibras aportam ao desempenho do concreto, seu emprego seria recomendável a estruturas em que a ductilidade é um dos parâmetros principais de projeto, ou naquelas feitas com concretos de alta resistência, uma vez que estes apresentam um comportamento ainda mais frágil que os concretos de resistência normal. Além disso, a utilização de fibras no combate aos esforços de cisalhamento mostra-se extremamente vantajosa e promissora. Neste trabalho, apresenta-se um modelo constitutivo para concreto reforçado com fibras baseado na formulação hipoelástica de ELWI E MURRAY (1979), originalmente proposta para concreto simples. As especificidades do comportamento do concreto com fibras frente às mais diversas solicitações, tais como, tração, compressão e cisalhamento, são incorporadas ao modelo através de relações tensão-deformação adequadas a esse material. Tais relações provêm de estudos analíticos e experimentais sobre o assunto, publicados na literatura técnica especializada. O modelo assim obtido é implementado no programa FEPARCS (ELWI E MURRAY, 1980), capaz de realizar análises númericas não-lineares através do método dos elementos finitos. Por fim, utiliza-se esse program para simular a resposta de uma estrutura de concreto com fibras, cujo ensaio experimental aparece minuciosamente descrito em (CRAIG, 1987). Os resultados numéricos obtidos são comparados com os experimentais correspondentes, em termos da curva carga versus deslocamento, desenvolvimento e distribuição de fissuras, progressão do escoamento da armadura longitudinal (convencional) e modo de ruptura. Avalia-se, então, a eficiência da implementação realizada na descrição do comportamento de estruturas de concreto com fibras. / [en] Along the past forty years, an increasing interest on composite materials formed by the addition of discrete fibers to a concrete matrix is being observed. These composites are known as fiber reinforced concretes. The interest on the use of fibers as reinforcement is justified by their significative contribution to concrete thoughness, as they reduce the brittle characteristics of that material. In fact, due to fiber reinforcement mechanism, fiber reinforced concrete can absorb much more strain energy until failure, in comparison to ordinary concrete. The cracking process seems to be more uniform, as the distance between cracks are reduced. Besides that, the material strength and the deformation levels at cracking and rupture are greater, on the case of fiber reinforced concrete. The benefits that fibers bring to concrete behavior indicate that they could be used as complementary reinforcent for concrete structures when ductility is a major design concern, or when high strength concrete is employed, as this class of material tends to be much more brittle then normal strength concrete. Fibers are also effective as shear reinforcement, and they could even replace stirrups in this function. In this work, a constituive model for fiber reinforced concrete is presented. This model is based on the formulation originally proposed by ELWI AND MURRAY (1979) for the case of ordinary concrete. The behavior characteristcs of fiber reinforced concrete are incorporated as adaquated uniaxial stress-strain relations in tension and compression. The behavior under shear stress is also considered. The model is then implemented in the finite element program FEPARCS (ELWI AND MURRAY, 1980). A numerical analysis on the response of a fiber reinforced concrete structure is conducted. Results reported in technical literature (CRAIG, 1987) are compared to those obtained by the finite element analysis. The efficiency of the model is then verified. / [es] En los últimos 40 anos, se ha observado un creciente interés por compuestos formados por la adición de fibras a matrizes de concreto, los llamados concretos reforzados con fibras. Ese interés se debe a la significativa ganancia en tenacidad que las fibras proporcionan, atenuando las características frágiles del concreto. De hecho, en virtud del mecanismo de refuerzo promovido por las fibras, el concreto con fibras es capaz de absorver mucha más energía de deformación hasta la ruptura, presentando, en el régimen posfisuración, un comportamiento mucho más suave que el concreto simple. Este comportamiento se ve acompañado por un proceso de fisuración más uniforme, en el cual se observan fisuras más finas y menos espaciadas. Además, se registran aumentos en los valores de resistencia del material y en los niveles de deformación que alcanza hasta su completa destrucción. Teniedo en cuenta los beneficios que las fibras aportan al desempeño del concreto, sería recomendable su empleo en extructuras donde la ductilidad es uno de los parámetros principales de proyecto, o en aquellas hechas con concreto de alta resistencia, ya que éstos presentan un comportamiento más frágil que los concretos de resistencia normal. En este trabajo, se presenta un modelo constitutivo para concreto reforzado con fibras que tiene como base la formulación hipoelástica de ELWI Y MURRAY (1979), originalmente propuesta para concreto simple. Las especificidades del comportamiento del concreto con fibras frente a las más diversas solicitudes, tales como, tracción, compresión y cisallamiento, se incorporan al modelo a través de relaciones tensión-deformación adecuadas a ese material. Tales relaciones provienen de estudios analíticos y experimentales sobre el asunto, publicados en la literatura técnica especializada. La implementación del modelo obtenido fue realizada a través del programa FEPARCS (ELWI Y MURRAY, 1980), capaz de realizar análisis númerico no lineal a través del método de los elementos finitos. Por fin, se utiliza ese programa para simular la respuesta de una extructura de concreto con fibras, cuyo ensayo experimental aparece minuciosamente descrito en (CRAIG, 1987). Los resultados numéricos obtenidos se comparan con los experimentales correspondientes, considerando la curva carga versus deslocamiento, desarrollo y distribuición de fisuras, progresión del deslizamiento de la armadura longitudinal (convencional) y modo de ruptura. Se evalúa entonces, la eficiencia de la implementación en la descrición del comportamiento de extructuras de concreto con fibras.

[pt] A Modelagem de Bacias Sedimentares e Sistemas Petrolíferos (BPSM) integra técnicas sofisticadas das disciplinas de geologia, engenharia e desenvolvimento de software, com o objetivo de mitigar os riscos econômicos da atividade de exploração de hidrocarbonetos. Entretanto, ainda que o processo de formação de bacias sedimentares apresente alta complexidade, envolvendo diversos processos geológicos, a técnica empregada majoritariamente em análises numéricas para representar o comportamento tensãodeformação das rochas sedimentares consiste em uma lei empírica desenvolvida pelo geofísico Lawrence F. Athy na década de 1930, fundamentada em uma série de simplificações sobre o problema mecânico. Neste contexto, o presente trabalho avalia a capacidade de modelos constitutivos baseados na mecânica do contínuo: elástico, elastoplástico com superfície de escoamento aberta e elastoplástico com superfície de escoamento fechada de representar o comportamento mecânico de rochas sedimentares associado aos processos de deposição e compactação observados na BPSM, usando os dados do modelo empírico de Athy como referência. Após a investigação inicial, os modelos constitutivos são comparados à solução tradicional (Lei de Athy) na representação de outros processos geológicos (erosão, compressão tectônica e extensão tectônica) usando cenários simplificados, com a finalidade de projetar o impacto desse tipo de representação em uma análise convencional de BPSM. / [en] Basin and Petroleum System Modeling integrates sofisticated techniques from the fields of geology, engineering and software development, aiming to mitigate the economic risks presented in the exploration of hydrocarbons. Even though the formation of sedimentary basins presents a high level of complexity, involving several geological processes, the main technique applied to represent the tension-deformation behavior of sedimentary rocks in numerical analyses is an empirical law developed by geophysicist Lawrence F. Athy in the 1930s, based on a series of simplifications about the mechanical problem. In that context, this work evaluates the capability of constitutive models based on Continuum Mechanics: elastic, elastoplastic with an uncapped yield surface and elastoplastic with a capped yield surface to represent the mechanical behavior of sedimentary rocks associated with the deposition and compaction processes observed in BPSM, using the data from Athy s empirical model as reference. After the initial investigation, the constitutive models are compared to the traditional solution (Athy s Law) in the representation of other geological processes (erosion, tectonic compression and tectonic extension) using simplified scenarios to predict the impact of such models in conventional BPSM analyses.

[pt] Essa pesquisa usa o Método do Ponto Material (MPM), para avaliar diferentes aspectos de deslizamentos de encosta. Esse tema é de suma importância, visto que os deslizamentos de terra são o desastre natural que mais causa perdas humanas no Brasil. Esse método numérico foi verificado a fim de que fosse encontrado o fator de segurança e a superfície de ruptura em um talude infinito com solo representado pelo modelo constitutivo de Drucker-Prager. Além disso, foi validado para calcular a força de impacto em um anteparo, sendo o volume impactante um objeto qualquer elástico ou um solo com modelo de Drucker-Prager. Os cálculos foram executados com o código desenvolvido pela PUC-Rio e produziram excelentes resultados. Foi observada uma grande dependência dos resultados com a malha, e assim como no Método dos Elementos Finitos, o refinamento da malha gera convergência para um resultado. O método foi considerado satisfatório para cálculo de uma parede de retenção em locais de risco. / [en] This research uses a numerical method, the Material Point Method, to evaluate different aspects of slope landslides. This theme is of paramount importance since landslides are the natural disaster that shows the highest number of deaths in Brazil. The method was verified so that it could find the safety factor and rupture surface in an infinite slope with soil represented by the Drucker-Prager constitutive model. In addition, the method was validated to calculate impact force in a bulkhead; the impacting volume could be any elastic object or soil with the Drucker-Prager model. The calculations were performed using the code developed by PUC-Rio and showed excellent results. It was observed high dependence on the mesh discretization, as well as the Finite Element Method, the refinement of the mesh generates convergence for a result. The method was considered satisfactory for calculating a retention wall at risk locations.

[pt] Durante a gelificação por parafina, a rede de cristais de parafina modifica o comportamento do óleo cru. Ele muda de um material newtoniano de baixa viscosidade para um material com dependência temporal e de alta viscosidade com tensão limite do escoamento. Com isto, é totalmente desafiador descobrir a pressão minima de início do fluxo de petróleo gelificado com uma microestrutura tão complexa. Através da minha dissertação de mestrado, investigamos dois materiais viscoplásticos: um gel de cabelo com uma tixotropia desprezível e uma suspensão aquosa tixotrópica 2 por cento de laponita para imitar o início de fluxo de óleos gelificados. Para ambos os materiais, o gradiente de pressão axial mínimo necessário para o início do fluxo foi medido e os valores medidos estavam de acordo com a previsão do balanço de força convencional. Por outro lado, os casos da indústria exibiram que o balanço de força mencionado acima leva a uma superestimação da pressão mínima inicial. Em alguns estudos, uma explicação elicitada é o comportamento tixotrópico do petróleo gelificado, mas nossos resultados acima mencionados serviram para refutá-lo. Durante a primeira parte da minha tese de doutorado, buscamos verificar por laboratório porque o balanço de força não se aplica ao petróleo gelificado, em seguida, buscamos uma explicação fisicamente adequada para essa discrepância e também uma maneira confiável de prever a pressão mínima de início do escoamento. Nesta linha, mostramos o efeito do encolhimento do petróleo gelificado na discordância entre a tensão limite estática de escoamento e o gradiente de pressão mínimo necessário para iniciar o fluxo, através da reometria e do fluxo de fluido em um tubo. Em seguida, introduzimos uma equação de balanço de força modificada com o efeito de encolhimento incluído para obter o melhor estimação da pressão minima de início de escoamento. Outro elemento essencial sobre o início de escoamento de petróleo gelificado é descobrir uma estratégia confiável para modelar matematicamente a reologia do material. Na maioria dos modelos que visam prever o comportamento reológico de petróleo gelificado, as mudanças da microestrutura durante o fluxo são consideradas tixotrópicas; dependência temporal reversível. Porém, observamos em nossos experimentos com histórias de fluxo e térmicas bem controladas que o caráter irreversível da dependência temporal é bastante evidente. Assim, na segunda parte da tese, propomos um modelo baseado no desenvolvimentos anteriores de Souza Mendes e colaboradores que considera a dependência temporal irreversível observada experimentalmente para petróleo gelificado. A capacidade preditiva do modelo proposto é então avaliada através de comparações com dados experimentais. / [en] Throughout the wax gelation, the network of parafinn crystals modifies the behavior of waxy crude oil. It changes from a low viscosity Newtonian to a high viscosity time-dependent material with yield strength. Now, it is totally challenging to find out the restart pressure for gelled crude oil ow with such a complex microstructure. Through my Master s dissertation, we investigated two viscoplastic materials, namely a hair gel with a negligible thixotropy and a quite thixotropic 2 percent aqueous suspension of Laponite to mimic the startup flow of waxy crude oils. For both materials, the minimum axial pressure gradient required for the onset of flow was measured, and the measured values were in good agreement with the prediction of a conventional force balance. On the other hand, industry cases have exhibited that the just mentioned force balance leads to an overestimation of the minimum startup pressure gradient. In some studies, an elicited explanation is the thixotropic behavior of the gelled crude, but our results above-mentioned served to falsify it. Over the first part of my PhD thesis, we aimed to verify in the laboratory that why the force balance does not hold for gelled crude oil and then we sought a physically proper explanation for this discrepancy and also a reliable way to predict the minimum startup pressure gradient. Along these lines, we show the role of gelled crude oil s shrinkage in the discordance between static yield strength and required minimum pressure gradient to onset the flow, through rheometry and fluid flow in a tube. Then, we introduce a modified force balance equation with the role of shrinkage included to best estimate the minimum restart pressure gradient. Another essential element through the restart flow of gelled waxy crude is to find out a reliable strategy to mathematically model the material s rheology. In most models that aim at predicting the rheological behavior of gelled waxy crude oil, the microstructure changes during ow are assumed to be thixotropic (reversible time dependent). But, we observed in our experiments with well-controlled flow and thermal histories that the irreversible character of time dependence is quite evident. Thus, in the second part of thesis we propose a model based on previous developments by Souza Mendes and co-workers that accounts for the irreversible time dependence observed experimentally in a waxy crude oil. The predictive capability of the proposed model is then assessed via comparisons with experimental data.

[pt] O presente trabalho tem como objetivo a implementação e avaliação de um modelo constitutivo para solos reforçados com fibra (compósito). A principal característica do modelo constitutivo implementado é que cada material (matriz de solo e fibra) segue sua própria lei constitutiva e ao mesmo tempo interagem entre si. Utilizando um algoritmo explícito, são implementados os modelos Cam Clay Modificado e Lade-Kim para a matriz de solo, cuja verificação é feita mediante o programa PLAXIS 2D e curvas tensão-deformação obtidas da literatura, respectivamente. Posteriormente, é adicionado o comportamento da fibra no desenvolvimento das tensões no compósito e verificado mediante a comparação das curvas tensão-deformação com as apresentadas por Diambra et al. (2013). As linguagens de programação utilizadas neste trabalho foram duas, a primeira é a utilizada no programa MATLAB, onde os códigos dos modelos são verificados e validados em relação à um conjunto de ensaios triaxiais de areia reforçada com fibra. Posteriormente foi usada a linguagem de programação FORTRAN para incluir o modelo constitutivo para solo reforçado com fibras no programa de elementos finitos ABAQUS, através da sub-rotina UMAT. Porém, para a implementação na sub-rotina UMAT os códigos dos modelos implementados no MATLAB sofrem algumas modificações com a finalidade de que o ABAQUS consiga compilar e representar adequadamento o comportamento do modelo constitutivo, mediante a correta utilização de vetores e propriedades desta. Finalmente, são modelados ensaios triaxiais drenados para verificar que a implementação mediante a sub-rotina UMAT é satisfatória. / [en] The present work aims to implement and evaluate a constitutive model for fiber-reinforced soils (composite). The main characteristic of the constitutive model implemented is that each material (soil and fiber matrix) follows its own constitutive law and at the same time interact with each other. Using an explicit algorithm, the Cam Clay Modified and Lade-Kim models are implemented for the soil matrix, verified by the PLAXIS 2D software and stress-strain curves obtained from the literature, respectively. Later, it is included the behavior of the fiber in the development of the stresses in the composite and verified by the comparison of the stress-strain curves with those presented by Diambra et al. (2013). The programming languages used in this work were two, the first one is the one used in the MATLAB program, where the codes of the models are verified and validated in relation to a set of triaxial tests of fiber-reinforced sands. Later the programming language was converted into FORTRAN to include the constitutive model for fiber reinforced soil in the ABAQUS finite element software, through the UMAT subroutine. However, for the implementation in the UMAT subroutine the codes of the models implemented in MATLAB undergo some modifications in order that ABAQUS can compile and represent adequately the behavior of the constitutive model through the correct use of vectors and its properties. Finally, drained triaxial tests are modeled to verify that the implementation through the UMAT subroutine is satisfactory.