[pt] O presente trabalho apresenta a Barragem de Terzaghi em um contexto da evolução tecnológica da construção de barragens de terra, colocando em foco sua construção e desempenho ao longo de 65 anos de operação bem como a verificação atual de suas condições físicas, analisando o impacto que a carga de sedimentos causa sobre a face montante da barragem, em termos de percolação, tensões e deformações. Tendo como base os métodos observacionais aplicados ao longo da operação da barragem e modelagem por elementos finitos, foram traçadas expectativas do fator de segurança para a seção principal da barragem ao longo do tempo. / [en] This work presents the Terzaghi Dam in a technological evolutionary context of the construction of earthfill dams, focusing on its own construction a performance after 65 years in operation as well as a scrutinized assessment of its physical conditions, implying the sediment load on its upstream slope, regarding seepage, stresses and strains. Taking into account observational methods applied throughout the operation of the dam and finite element simulations, the evolution of factor of safety over the dam lifetime and on was evaluated for the main section of the dam.

[pt] O design otimizado de estruturas considerando não-linearidades tem sido amplamente pesquisado nas décadas recentes. A análise de elementos finitos aplicada à otimização topológica é prejudicada pela deformação excessiva de elementos de baixa densidade sob alta compressão, o que impede o processo de encontrar uma solução ótima. Dois métodos, o esquema Interpolação de Energia e a técnica de Hiperelasticidade Aditiva, são implementados para superar essa dificuldade no problema de minimização da flexibilidade, e modelos de materiais hiperelásticos são usados para investigar suas influências na topologia otimizada. O Método das Assíntotas Móveis é usado para atualizar as variáves de projeto cujas sensibilidades foram calculadas pelo método adjunto. A equação de estado é resolvida através do método de Newton-Raphson com um incremento de carga ajustável para reduzir o custo computacional. Resultados de dois problemas de referência são comparado com aqueles já estabelecidos na literatura. O uso de diferentes modelos hiperelásticos apresentou pouca influência no design final da estrutura. O método de Interpolação de Energia foi capaz de convergir para cargas muito maiores que o método padrão, enquanto a Hiperelasticidade Aditiva apresentou dificuldades de convergência em estado plano de deformação. / [en] The optimized design of structures considering nonlinearities has been widely researched in the recent decades. The finite element analysis applied to topology optimization is jeopardized by the excessive deformation of low-density elements under high compression, which hinders the process of finding an optimal solution. Two methods, the Energy Interpolation scheme and the Additive Hyperelasticity technique, are implemented to overcome this difficulty in the minimum compliance problem, and hyperelastic material models are used to investigate their influence on the optimized topology. The Method of Moving Asymptotes is used to update the design variables whose sensitivities were calculated from the adjoint method. The state equation is solved through the Newton-Raphson method with an adjusting load step to reduce computational cost. Results for two benchmark problems are compared with those already established in the literature. The use of different hyperelastic models presented little influence on the final design of the structure. The Energy Interpolation method was able to converge for much higher loads than the default method, while the Additive Hyperelasticity presented convergence difficulties in plane strain.

[pt] Na análise da predição de vida de estruturas a presença de defeitos por trincas é de fundamental importância. Estudos iniciais utilizando a Mecânica da Fratura Linear Elástica mostram que o Fator de Intensidade de Tensão controla o crescimento de trinca por fadiga. Contudo, sobrecargas podem induzir efeitos de memória no material ao reduzir, parar ou acelerar taxas de propagação, comportamento não previsto considerando-se um único parâmetro elástico. Devido às hipóteses consideradas, as soluções fornecidas pelos métodos analíticos não verificam condições importantes de equilíbrio e compatibilidade, especialmente quando são usados campos singulares idealizados de tensão e de deformação na região da ponta da trinca. Neste trabalho uma extensa revisão dos resultados analíticos recentemente publicados para materiais nos regimes elástico e elastoplástico são apresentadas: Williams e HRR com campo de distribuição singular das tensões e Creager-Paris com campo não-singular das tensões com o material elástico. Neste estudo, estas soluções são comparadas com resultados da análise numérica utilizando a discretização por elementos finitos tridimensionais elastoplásticos com o critério de escoamento do material de von Mises na representação do cegamento da ponta da trinca e importantes conclusões são apresentadas com referência aos limites de aplicação dos modelos teóricos bem como a extensão dos requerimentos para o modelo numérico. Além de carregamentos monotonicamente crescentes, o descarregamento foi também considerado na análise numérica, apesar desta condição não estar presente nas referências consideradas para as soluções analíticas. Para as simulações numéricas consideradas algoritmos de elastoplasticidade foram implementados em um framework com arquitetura baseada em plugins. / [en] In structure life prediction analysis the occurrence of crack defects are of paramount importance to be considered. Basic studies using Linear Elastic Fracture Mechanics approach shows that the Stress Intensity Factor (SIF) parameter controls Fatigue crack growth (FCG). However, overloads may induce material memory effects that delay, arrest or accelerate the FCG rate, a behavior that is not described by considering a single elastic parameter. To account for service variable amplitude loadings, the use of a prescribed stress-strain distribution has been proposed in recent research studies, as the driving force of FCG, using the critical damage approach. Due to the assumptions considered in the analytical derivations, important equilibrium and compatibility conditions are violated in the obtained solutions, since an idealized singular stress-strain field at the crack front region is used. In this work a comprehensive review of the recently published analytical results for solutions under the elastic and elastoplastic material behavior regimens are presented: Williams and HRR with singular stress distribution field and Creager-Paris with a non singular stress field, but with the material elastic approach. These solutions are compared, throughout the study, to results obtained from the numerical analysis using a 3D finite element discretization with elastoplastic von Mises yielding criteria employed for the modeling of the crack tip blunt and some comprehensive conclusions are derived regarding application limits of the theoretical models as well as the required extent of the numerical model representation. In addition to monotonic increasing applied loads, unloading conditions were also considered in the numerical analysis, although no reference to this condition is available in the literature, considering the analytical approach. For the numerical simulations presented elastoplastic algorithms were implemented in a plugin based framework.

[pt] Entre materiais renováveis de baixo impacto ambiental o bambu se destaca devido ao seu crescimento rápido, alta resistência à tração, baixo peso específico e flexibilidade. É um material compósito natural com fibras longas alinhadas unidirecionalmente, anisotrópico, que não está sendo usado amplamente nas construções civis por razão de falta de informações técnico-científicas. O bambu é um material inteligente com a distribuição das fibras longitudinais como functionally graded composite material (FGM). O objetivo deste trabalho é continuação dos estudos sobre o comportamento de bambu para aplicação nas obras de engenharia, que inicio na PUC-Rio em 1979. Neste estudo específico foi observado comportamento de bambu da espécie Phyllostachys pubescens (Mosô) (através dos ensaios experimentais) sujeito a carga de torção e flexão pura (realizados com utilização do sistema de reconhecimento de imagem) e simulações computacionais em programa de elementos finitos Abaqus. Os ensaios foram feitos para os corpos de prova inteiros (de 1 metro de comprimento) e os mesmos corpos de prova foram posteriormente fissurados no meio, alargando a fissura para cada ensaio até a ruptura. / [en] Among renewable materials with low environmental impact bamboo stands out because of rapid growth, high tensile strength, low specific weight and flexibility. It is a natural composite anisotropic material with unidirectionally aligned long fibers, which is not widely used in civil construction yet due to lack of technical and scientific information. Bamboo is an intelligent material with the distribution of longitudinal fibers as functionally graded composite material (FGM). The objective of this work is to continue the study the behavior of bamboo for application to engineering works, which started at PUC-Rio in 1979. In this specific study was observed the behavior of Moso bamboo (Phyllostachys pubescens) through the experimental observation of pure torsion and bending (realized with utilization of image recognition system) and computer simulations in finite element program Abaqus. The tests were done for the whole bamboo specimens (1 meter long) and the same specimens were subsequently cracked in the middle extending the crack for each test until failure.

[pt] Este trabalho apresenta o desenvolvimento/montagem de um sistema computacional para análise de fluxo em meios porosos, meios fraturados, porosos fraturados e em combinações destes meios, considerando regime permanente ou transiente, sob condições saturadas e não saturadas. O sistema consiste de quatro programas, três programas de funções específicas interligadas por rotinas de programação feitas na linguagem Cmaismais e o quarto é um visualizador de resultados. O FracGen 3D (Telles, 2006) gera fraturas ou famílias de fraturas de forma determinística ou probabilística. O programa ICEM CFD v.14 divide o domínio de interesse em sub-dominios, através da geração de malha de elementos finitos. O programa FTPF-3D (Telles, 2006) utiliza o método de elementos finitos para discretizar as equações governantes no espaço e em diferenças finitas no tempo, e para resolver a não linearidade, utiliza o método iterativo de Picard ou o método iterativo BFGS e finalmente O Pos3D é o responsável pela visualização dos resultados. Neste trabalho foram desenvolvidos cinco exemplos, dois deles para a validação deste procedimento, e três aplicados a um talude típico do Rio de Janeiro, os quais incluem fraturas verticais e juntas de alívio. Estes casos estudados verificam a influência das fraturas nos meios porosos em termos de carga de pressão, totais e campo de velocidades, para a verificação do comportamento hidráulico dos maciços e de eventuais instabilidades. / [en] This work presents the development/assembly of a computational system for flow analysis in porous media, fractured and fractured porous media and in combination of both media, considering steady or transient states under saturated and unsaturated conditions. The system comprehends four computational programs, three of them of specific functions interconnected by Cplusplus programing routines and the last program is an output viewer. FracGen 3D program (Telles, 2006) generates fractures or fracture families in a determinist or probabilistic way. ICEM CFD v.14 program divides the interest domain in sub-domains by means of the element finite mesh generation. FTPF-3D program (Telles, 2006) uses the element finite method to discretize the governing equations in the space domain and the difference finite method for the time domain and for solving the nonlinearity is used the iterative Picard or BFGS method, so that, finally, Pos3D viewer program is answerable by visualization of the results. In the present dissertation five examples were developed, two of them for the validation of this procedure and the three others applied to a typical slope in Rio de Janeiro, which include vertical fractures and relief joints on their slopes. All those studied cases evaluate the influence of the fractures on porous media in terms of pressure and total heads and velocity fields for verifying of the hydraulic behavior of solid masses and eventual instabilities.

[pt] Emulsões podem ser usadas como agentes de controle de mobilidade em métodos de recuperação avançada de óleo de forma a conseguir um varrido mais eficiente do reservatório. A aplicação de tal técnica requer um entendimento completo de como as emulsões escoam em um material poroso. O comportamento macroscópico do escoamento deve ser determinado baseado no escoamento na escala de poros, em que as emulsões não podem ser tratadas como líquidos não Newtonianos monofásicos, uma vez que os diâmetros das gotas são da mesma ordem de magnitude das gargantas de poros. Experimentos realizados por Cobos et al.(1) sobre o escoamento de emulsões através de micro-capilares com garganta, que serviu como modelo para a geometria de uma garganta de poro conectando dois poros adjacentes, mostrou o efeito da fase dispersa na queda de pressão para diferentes condições de escoamento e características das emulsõs. De forma a ampliar a faixa de número de capilaridade e tamanhos de gota explorados nos experimentos, determinar o efeito da tensão interfacial, razão de viscosidades, geometria do capilar e determinar características do escoamento com mais de uma gota, é estudado o escoamento de gotas suspensas em uma fase contínua escoando através de um micro-capilar. A evolução da interface das gotas é determinada pelo método de level-set, que é acoplado a solução das equações de Navier Stokes, baseada no método de elementos finitos. Foi investigado o efeito da distância inicial entre duas gotas, assim como o efeito da tensão interfacial na velocidade de cada uma das gotas. Os resultados obtidos fornecem uma descrição mais detalhada do escoamento de emulsões em gargantas de poros. / [en] Emulsions can be used as mobility control agent in enhanced oil recovery methods in order to achieve a more efficient sweep of the reservoir. The application of such technique requires full understanding of how emulsions flow in a porous material. The macroscopic flow behavior can be determined based on the pore scale flow, at which emulsions cannot be treated as a single phase non Newtonian liquid, since the drop diameters are in the same order of magnitude of the pore throats. Experiments by Cobos et al.(1) of flow of emulsions through constricted micro capillaries, wich served as a model to the geometry of a pore throat connecting two adjacent pore bodies, have shown the effect of the dispersed phase on the pressure drop for different flow conditions and emulsions characteristics. In order to widen the range of capillary number and drop size explored in the experiments and to determine the effect of interfacial tension, viscosity ratio, capillary geometry and determine the characteristics of the flow with more the one drop, we study the flow of a drops suspended in a continuous phase flowing through a constricted micro capillary. The evolution of the drop interface is determined by the level set method which is incorporated to a fully coupled implicit Navier Stokes solver based finite element method. We investigate the effect of an initial distance between the two drops, as well as the effects of the interfacial tension on the drops speed. The results provide a more detailed description of the flow of emulsions through pore throats.

[pt] O rebaixamento e controle das águas subterrâneas são atividades implementadas e monitoradas de forma contínua ao longo da vida dos projetos de mineração subterrânea ou a céu aberto. A implementação apropriada e eficiente destas atividades depende de estudos hidrogeológicos de grande porte, que permitem avaliar os sistemas de controle mais adequados. A procura da eficiência técnico-econômica destes processos demanda análises numéricas de fluxo tridimensionais de toda a região em estudo, caracterizada por profundas e complexas estratificações de materiais permeáveis abaixo do lençol freático, como normalmente abrangem projetos de mineração, onde as soluções analíticas não podem mais ser aplicadas. O presente trabalho de pesquisa contribui na melhor compreensão das formulações numéricas que representam o comportamento do fluxo subterrâneo, através de dois estudos de caso, o primeiro em uma mina subterrânea e segundo em uma mina superficial. No caso da mina subterrânea foram incorporadas feições cársticas através de elementos discretos 1D dentro de um modelo tridimensional de elementos finitos com o intuito de representar caminhos preferenciais de fluxo. Foram discutidas as vantagens de incorporar tais feições de forma explicita, quantificando o fluxo que passam por estas, que alimentam à mina através de conexões diretas com um rio adjacente. Estes elementos discretos permitem uma representação mais realista do meio hidrogeológico e ao mesmo tempo, uma avaliação mais aprimorada dos efeitos no comportamento do fluxo subterrâneo devido à impermeabilização superficial do rio, como a solução mais coerente para este problema de infiltração. Também foi elaborado um modelo hidrogeológico conceitual para representar o comportamento hidrogeológico de uma mina a céu aberto, desenvolvendo uma sistemática de uso adequado das condições de contorno e de restrição, a calibração deste modelo e a verificação de diferentes cenários de fluxo, como resultado da incorporação das diferentes técnicas de controle das águas avaliadas em regime transiente, que abrange poços de bombeamento, paredes cut-off assim como ponteiras filtrantes e drenos horizontais. As diferentes técnicas modeladas mostraram resultados satisfatórios, sendo que arranjos de várias técnicas, configuradas de forma localizada, resultam ser mais recomendáveis e eficientes para tratar problemas particulares. / [en] Dewatering and groundwater control are activities continuously implemented and monitored throughout the duration of underground mining or open pit projects. The proper and efficient implementation of these activities depends on large hydrogeological studies, for assessing the most appropriate control systems. The demand for technical and economic efficiency of these processes requires three-dimensional flow numerical analysis of the entire study area, characterized by deep and complex stratifications of permeable materials below the water table, as usually cover mining projects, where analytical solutions cannot be applied. The present research contributes to a better understanding of the numerical formulations that represent the behavior of groundwater flow via two case studies—the first in an underground mine and the second in a surface mine. In the case of the underground mine, 1D discrete elements were incorporated within a finite-three dimensional model in order to represent preferential flow paths. The advantages of incorporating such features explicitly to quantify the flow passing through them, that feed the mine through direct connections with an adjacent river, were discussed. The use of discrete elements allows for a more realistic representation of the hydrogeological environment and, at the same time, a more refined assessment of the effects on the behavior of groundwater flow due to surface sealing of the river, as the most coherent solution to this infiltration problem. Furthermore, a conceptual hidrogeologic model representing the hydrogeological behavior of an open pit mine was created, developing a system of the appropriate use of boundary and constraint conditions, calibrating the model and verifying different flow scenarios, as a result of the incorporation of different water control techniques assessed in transient regime, such as covering pumping wells, cut-off walls as well as wellpoints and horizontal drains. The different techniques modeled showed satisfactory results, with various arrangements of techniques configured in a localized form resulting in the most desirable and efficient treatments for particular problems.

[pt] Esta pesquisa tem como objetivo o desenvolvimento de uma metodologia com base no método de elementos finitos para o estudo do problema de contato com atrito que surge na interação solo-duto no comportamento de dutos enterrados. O tratamento das restrições de contato e a incorporação da lei do atrito (tipo Coulomb) é feita através do Método da Penalidade, tomando como referência o trabalho de Laursen e Simo (2002), onde as restrições de contato são impostas de maneira aproximada com o uso de parâmetros de penalidade para as forças normal e tangencial no contato. As relações cinemáticas são descritas em termos de uma função diferenciável da distancia entre os corpos (gap). Para a discretização são empregados elementos hexaedricos de 8 nós com uma formulação híbrida - Enhanced Assumed Strain (EAS) e os efeitos de comportamento não-linear dos materiais envolvidos são considerados na presença de grandes deslocamentos e grandes deformações. A teoria da plasticidade é utilizada para modelar a natureza não-linear das relações constitutivas do duto. Aplicações considerando o problema de contato com atrito que surge na interação solo-duto são apresentadas. / [en] This research aims the development of a methodology based on the finite element method for 3D frictional contact problems such as the frictional contact problem that arises from the soil-pipe interaction of buried pipelines. The treatment of the contact restrictions and the incorporation of the friction law of the Coulomb type are carried out through a penalty method, where the contact restrictions are imposed in an approximated manner using penalty parameters for both normal and tangential forces. The cinematic relations are established in terms of the diferential funtion of the gap between bodies. Hexahedral eight-node elements are employed based on the Enhanced Assumed Strain (EAS) concept and the effects of the non-linear behavior of the materials are considered in presence of large displacements and large deformations. The theory of plasticity is used to model the non-linear nature of the constitutive relations of the pipe. Applications are presented considering the frictional contact problem that arises on the interaction surfaces of a buried structure such as an oil pipeline.

[pt] Reservatórios de petróleo e gás estruturalmente compartimentados por falhas geológicas selantes são encontrados em diversas regiões do mundo. Durante a fase de explotação, a integridade do selo destas falhas pode ser comprometida pelas deformações decorrentes dos processos de depleção e/ou injeção de fluidos. Estas deformações, em conjunto com as propriedades físicas e geométricas das rochas e falhas presentes, podem alterar significativamente o estado de tensões do maciço rochoso fazendo com que uma falha reative e se torne hidraulicamente condutora. A esse fenômeno estão associados riscos de exsudação, perda de integridade de poços e outros potencias problemas geomecânicos. Na literatura, diversas modelagens numéricas têm sido utilizadas a fim de caracterizar e prever os fenômenos de reativação e/ou abertura de falhas geológicas. A maior parte de estas abordagens faz uso de modelos bidimensionais considerando seções críticas na hipótese de estado plano de deformação. Essas simplificações são adotadas a fim de evitar a complexidade geométrica e o alto custo computacional de uma modelagem tridimensional. No entanto, a configuração tridimensional dos planos de falha pode induzir a reativação em direção a zonas mais críticas do que aquelas contidas numa única seção. Neste trabalho apresenta-se uma metodologia para análise de reativação de falhas geológicas e discute-se a importância do uso dos modelos 3D na previsão do comportamento geomecânico de reservatórios compartimentados por falhas geológicas. São apresentados 3 modelos diferentes. O primeiro exemplo traz um modelo bidimensional apresentado na literatura, faz-se uma comparação dos resultados com representação por meio do elemento de interface, por meio do continuo equivalente e por meio de um elemento solido com fraturas embutidas. O segundo exemplo faz-se um comparativo entre a utilização de elementos quadrilaterais e triangulais para a representação da falha em modelos 3D. Para o terceiro modelo foram realizadas simulações numéricas considerando modelos 2D e 3D em um simulador in-house baseado no método dos elementos finitos. Para a representação do meio continuo foram utilizados elementos quadrilaterais para o caso 2D, e elementos hexaédricos e tetraédricos para o caso 3D. Para a representação das falhas geológicas foram utilizados elementos de interface de espessura nula segundo o critério de ruptura de Mohr-Coulomb. Da comparação dos resultados, constata-se que as análises 2D e 3D forneceram previsões de reativação similares. No entanto, as previsões de pressões de abertura foram distintas em ambos os modelos devido às diferentes trajetórias de migração de fluido. Particularmente em modelos com geometria irregular confirma-se a importância do emprego de modelo 3D. / [en] Oil and gas reservoirs that are structurally compartmented by sealing geological faults are common in several areas around the world. During production, the deformations from the processes of fluid depletion and/or injection can compromise the integrity of the seal of the faults. This deformation, together with the physical and geometrical properties from the rocks and faults can significantly change the stress state. Therefore, it might cause fault reactivation, turning it in a hydraulic conduit. Related to this phenomenon, are the exudation, loss of wellbore integrity and other potential geomechanical problems. There are several numerical modelling techniques available in literature to characterize and predict the reactivation and/or opening of geological faults. In most of these modelling approaches, bi-dimensional models are used for critical sections through the assumption of plane strain conditions. The reason for using 2D models is to avoid the geometrical complexity and the high computational costs associated to three-dimensional modeling. On the other hand, the fault planes in the three-dimensional approach can show fault reactivation in a more critical direction e than the one represented by the bi-dimensional model. In this work, a methodology is presented in order to assess geological fault reactivation. In addition, the importance of using 3D models in the prediction of the geomechanical behavior of reservoirs compartmented by geological faults is discussed. Three different models are presented. The first example is based on a two dimensional model from the literature. A comparison between approaches using interface elements, equivalent continuum elements and solid element with fractures is carried out in the first example. The second example brings a comparison between the quadrilateral and triangular elements to represent faults in a 3D model. In addition, an analysis was carried out considering 2D and 3D models using an in house software based on the finite element method. To simulate the continuum medium, quadrilateral elements are used in the 2D case and in the 3D case hexahedral and tetrahedral elements are employed. In addition, to represent the geological faults, interface elements with zero thickness are used in association with the Mohr-Coulomb failure criterion. In the case study, predictions of fault reactivation were similar in the 2D and 3D models. However, fault opening pressures were different in both models, due to the 3D fluid migration path. It also confirmed the importance of using 3D models when simulating irregular geometries.

[pt] O escoamento de cápsulas suspensas em uma fase líquida através de canais e capilares micrométricos representa um problema complexo que ocorre em diferentes aplicações, de glóbulos vermelhos em hemodinâmica até escoamento em meios porosos. Em aplicações de meios porosos, a compreensão da dinâmica na microescala é fundamental para avaliar o comportamento macroscópico do escoamento. Canais e capilares com constrição podem ser usados para modelar uma garganta conectando dois poros adjacentes. O escoamento de uma cápsula suspensa através de tais modelos foi analisado para avaliar as características do escoamento considerando os efeitos inerciais (isto é, número de Reynolds finito), incluindo a máxima diferença de pressão necessária para empurrar uma cápsula através da constrição em função do raio da cápsula, a tensão inicial e o material da membrana, geometria do canal e do capilar, assim como as condições de escoamento. De fato, neste estudo, a resposta da pressão é fundamental para avaliar o efeito de bloqueio da cápsula. As fases líquidas internas e externas foram descritas pelas equações de Navier-Stokes, enquanto que a dinâmica da membrana da cápsula foi modelada por uma estrutura flexível 1-D tipo mola. O problema de interação fluido-estrutura foi resolvido usando o método de elementos finitos acoplado ao método de fronteira imersa. Os resultados mostraram a redução da mobilidade da fase contínua devido à presença da cápsula através da constrição. Tais resultados podem ser usados para projetar microcápsulas para bloquear caminhos preferenciais de fluxo da água no processo de deslocamento de óleo em meios porosos. / [en] The flow of capsules suspended in a liquid phase through small channels and capillaries poses a complex problem presented in different applications, from red blood cells on hemodynamics to flow in porous media. In applications of porous media, the understanding of microscale dynamics is fundamental to assess the macroscopic flow behavior. Constricted channels and capillaries can be used to model a pore throat connecting two adjacent pore bodies. The flow of a suspended capsule through such models was analyzed to evaluate the flow characteristics considering inertial effects (i.e. finite Reynolds numbers), including the maximum pressure difference required to push a capsule through the constriction as a function of capsule radius, initial membrane tension, membrane material, channel and capillary geometries, as well as flow conditions. In fact, in this study, the pressure response is fundamental in order to assess the capsule blocking mechanism. Inner and outer liquid phases were described by the Navier-Stokes equations and capsule membrane dynamics was modeled by a 1-D spring-like flexible structure. The fluid-structure interaction problem was solved using the finite element method coupled with the immersed boundary method. Results showed the mobility reduction of the continuous phase due to the presence of a capsule as it flows through the constriction. Such results can be used to design microcapsules to block preferential water flow paths in oil displacement process in porous media.