The thesis provides the required knowledge base for establishing Laser Induced Stress Wave Thermometry (LISWT) as a viable alternative to current infrared technologies for temperature measurement up to 1000°C with ±1°C resolution. The need for a non-contact, high resolution thermal measurement methodology applicable to Rapid Thermal Processing (RTP) motivated the work. A stress wave propagation model was developed and a complex, temperature-dependent elasto-viscoplastic constitutive law was identified. A stagger-grid finite difference scheme was followed to approximate the solution field subject to temperature and plate thickness variations. Extensive numerical experiments were conducted to identify the proper time and spatial steps. A Gabor wavelet transform scheme was also employed for the extraction of wafer thermal and geometric information from exploring wave attenuation and dispersion. Researched results concluded that wave group velocity is a nonlinear function of temperature. Nonlinearity became more prominent at high temperatures and low frequencies. As such, for LISWT to achieve better thermal resolution at high temperatures, low frequency components of the induced stress wave should be exploited. The results also showed that the influence of temperature on attenuation is relatively small. It is not recommended to use attenuation for resolving temperature variation as small as several degrees Celsius. In addition to temperature, geometry also was found to have an impact on wave dispersion and attenuation. The results showed that the influence of thickness on wave velocity is significant, thus suggesting that for LISWT to achieve high temperature resolution, wafer thickness must be accurately calibrated in order to eliminate all possible errors introduced by thickness variation. The study established the basic framework for LISWT to be applicable to silicon wafer RTP at elevated temperatures. The model and methods developed for the course of the research can be easily adapted to account for other nondestructive evaluation applications involving the use of surface, plate or bulk waves for material characterization and thermal profiling.

The thesis provides the required knowledge base for establishing Laser Induced Stress Wave Thermometry (LISWT) as a viable alternative to current infrared technologies for temperature measurement up to 1000°C with ±1°C resolution. The need for a non-contact, high resolution thermal measurement methodology applicable to Rapid Thermal Processing (RTP) motivated the work. A stress wave propagation model was developed and a complex, temperature-dependent elasto-viscoplastic constitutive law was identified. A stagger-grid finite difference scheme was followed to approximate the solution field subject to temperature and plate thickness variations. Extensive numerical experiments were conducted to identify the proper time and spatial steps. A Gabor wavelet transform scheme was also employed for the extraction of wafer thermal and geometric information from exploring wave attenuation and dispersion. Researched results concluded that wave group velocity is a nonlinear function of temperature. Nonlinearity became more prominent at high temperatures and low frequencies. As such, for LISWT to achieve better thermal resolution at high temperatures, low frequency components of the induced stress wave should be exploited. The results also showed that the influence of temperature on attenuation is relatively small. It is not recommended to use attenuation for resolving temperature variation as small as several degrees Celsius. In addition to temperature, geometry also was found to have an impact on wave dispersion and attenuation. The results showed that the influence of thickness on wave velocity is significant, thus suggesting that for LISWT to achieve high temperature resolution, wafer thickness must be accurately calibrated in order to eliminate all possible errors introduced by thickness variation. The study established the basic framework for LISWT to be applicable to silicon wafer RTP at elevated temperatures. The model and methods developed for the course of the research can be easily adapted to account for other nondestructive evaluation applications involving the use of surface, plate or bulk waves for material characterization and thermal profiling.

A verified elasto-viscoplastic finite element model is used to develop a better understanding of the performance of embankments with geosynthetic reinforcement constructed over rate-sensitive soil. The interaction between reinforcement and prefabricated vertical drains (PVDs) and their effects on time-dependent behaviour of embankments are examined. For rate-sensitive soils, the generation of creep-induced pore pressures following the end of construction is evident along the potential slip surface. As a result, the minimum factor of safety with respect to embankment stability occurs after the end of construction. The combined use of reinforcement and PVDs are shown to provide an effective means of minimizing creep-induced excess pore pressure, increasing overall stability, and decreasing deformation of the embankments. The combined effects of the viscoelastic properties of geosynthetic reinforcement (polyester, polypropylene and polyethylene) and the rate-sensitive nature of foundation soils on the performance of embankments are examined. The effect of various factors, including reinforcement type (i.e., stiffness and viscosity), soil viscosity, construction rate and allowable long-term reinforcement strain, on the time-dependent behaviour of embankments are considered. The long-term performance of reinforced embankments is investigated for different maximum allowable long-term reinforcement strains. From a series of finite element analyses, the ideal allowable reinforcement strains to minimize embankment deformation while providing optimum long-term service height of the embankment, considering the effect of soil and reinforcement viscosity, are proposed for soils similar to those examined in this study. The currently proposed design methods for embankments with creep-susceptible reinforcement over rate-sensitive soils appear to be overly conservative. This study proposes a refined approach for establishing the allowable long-term reinforcement strains that are expected to provide adequate performance while reducing the level of conservativeness of reinforced embankment design. Finally, a previously developed elasto-viscoplastic constitutive model is modified to incorporate the effect of soil structure using a state-dependent fluidity parameter and damage law. The model was evaluated against data from a well-documented case study of a reinforced test embankment constructed on a sensitive Champlain clay deposit in Saint Alban, Quebec. The benefit of basal reinforcement and the effect of reinforcement viscosity are then discussed for these types of soil deposits. / Thesis (Ph.D, Civil Engineering) -- Queen's University, 2011-01-21 22:26:40.133

[pt] Atualmente, o escoamento de fluidos não Newtonianos está presente em diversos processos industriais e também no dia- a-dia de todos. Na presente dissertação, inicialmente, é apresentada uma nova equação constitutiva capaz de modelar fluidos que possuam tanto elasticidade como viscoplasticidade, sendo este novo modelo definido como elasto-viscoplástico. Para a melhor compreensão do modelo, anteriormente, é realizada, uma breve revisão dos fluidos não Newtonianos, com maior ênfase para o modelo viscoplástico SMD e o modelo viscoelástico de Oldroyd-B, que possuem grande influência no novo modelo constitutivo. Em seguida é apresentada uma nova metodologia de adimensionalização levando-se em conta apenas os parâmetros reológicos neste processo, garantindo a independência das equações com relação a vazão. Por fim, é realizada a simulação deste novo modelo em um escoamento com uma expansão e contração abrupta. Nesta simulação, é analisado o comportamento do fluido ao passar por esta cavidade com relação às principais variáveis da equação constitutiva. A simulação foi realizada pelo método dos elementos finitos e os resultados mostram a influência do número de Deborah reológico, da velocidade, do expoente power-law e da razão entre os tempos de retardo e de relaxamento, na perda de carga e no padrão de escoamento. / [en] Currently, the flow of non-Newtonian fluids is present in many industrial processes and also in day-to-day of all. In this dissertation, initially it is exposed a new constitutive equation capable of modeling fluids that present as much elasticity as viscoplasticity, and this new model is called as Elasto-viscoplastic. For better understanding of this model, previously, a brief review of non- Newtonian fluids is done, with a focus on the SMD viscoplastic model and on the Oldroyd-B viscoelastic model, both which have a great influence in this new constitution model.Afterward, a new methodology of non- dimensionalization is presented, taking into account only the rheological parameters in this process and ensuring the independence of equations with respect to flow. Finally, the simulation of this new model is performed in a flow within a abrupt expansioncontraction geometry. In this simulation, it is studied the behavior of the fluid through this cavity with respect to key variables of the constitutive equation. The simulation was carried out by the finite element method and, by the analysis of the results, it is shown the influence of the rheological Deborah, flow speed, the power-law exponent and the ratio of the time of retardation and relaxation on the head loss and the flow pattern.

O estudo do comportamento reológico de fluidos não Newtonianos tem grande importância em diversas áreas da engenharia. O aumento na demanda destes fluidos por exemplo, no uso doméstico, pessoal e em processos químicos acarreta em dificuldades que vão desde o processo de sua mistura ao seu manuseio. Dentre os fluidos não Newtonianos estão os viscoelásticos, os quais exibem deformação aparente quando os níveis de tensões são inferiores ao limite de escoamento do material e, dentro desta classe, alguns ainda podem apresentam comportamento elástico quando submetidos à baixa taxa de cisalhamento. Juntamente com os efeitos elasto-viscoplásticos, os materiais podem apresentar comportamento tixotrópico, onde devido as tensões sua reestruturação não é instantânea. Na presente Tese, fez-se um estudo numérico a fim de analisar o problema especifico de escoamentos de fluidos elasto-viscoplásticos com comportamento tixotrópico em uma expansão planar abrupta na razão de aspecto 4 : 1. O modelo mecânico aplicado consiste de uma equação viscoelástica para o campo de tensões, uma evolutiva para o parâmetro de estrutura do material, bem como as equações de conservação de massa e momento. O modelo mecânico aplicado mostrou-se capaz de prever o comportamento tixotrópico. A aproximação numérica do modelo aplicado foi feita através do método estabilizado de elementos finitos, especificamente o método Galerkin Mínimos-Quadrados (GLS), o qual foi implementado no código de elementos finitos para fluidos não Newtonianos em desenvolvimento no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) da Universidade Federal do Rio Grande do Sul. Os fenômenos reológicos presentes no problema foram analisados a partir da influência da cinemática do escoamento, da elasticidade e da tixotropia, no nível de estruturação do material, na posição e tamanho das regiões não escoadas e na deformação elástica do material. Os resultados mostraram-se satisfatórios, pois condizem com os apresentados na literatura, apontando a boa predição do modelo mecânico aplicado bem como a robustez de sua implementação computacional. / The rheological behavior of non-Newtonian uids study is of great then importance in many areas of engineering. The increase in demand of these uids - for example, domestic use, personal and processes chemical - causes di culties ranging from the process of mixing it to handling. Among the non-Newtonian uids are viscoelastic, which exhibit apparent deformation when stress levels are lower than the yield limit of the material and, within this class, some still have elastic behavior when subjected to low shear rates. Together with the elastic-viscoplastic e ects, materials may exhibit thixotropic behavior, ie, due the restructuring strain is not instantaneous. In this thesis was made a numerical study to simulate the speci c problem of ow of elastic-viscoplastic uids with thixotropic behavior in abrupt planar expansion { common geometry in industrial systems associated with elastic-viscoplastic uids, whose ratio the aspect is 4:1. The mechanical behavior of most of these structured materials, are highly non-Newtonian, with this, there is need to make them more pseudoplastic, causing undesirable behaviors such as thixotropic which is a phenomenon somewhat characterized and modeled in the literature. The mechanical model applied is able to predict thixotropic behavior and is composed of a viscoelastic equation for the stress eld and an evolutionary to the material structure parameter in addition to the mass and momentum conservation equations. This mechanical model is approximated by a stabilized nite element model, called the Galerkin method of least squares (GLS). In order to study the rheological phenomena present in the problem is analyzed the in uence of the ow kinematics, elasticity and thixotropy in the level of structure of the material, in position and size of unyielded regions and the elastic deformation of the material. The results were satisfactory, since the the study of intensity the U the results agree with those reported in the literature pointing to good prediction of the mechanical model applied well as the hardiness of their computational implementation.The results of the elasticity showed quite spectacular behavior of unyielded regions of the material since for high values of relaxation time, the unyielded region of channel larger has the form of " ngers"because of the high exibility of the material along the line of symmetry of the channel. The results of the analysis of the elastic deformation show, that model correctly dosing the elasticity in the unyielded regions. The thixotropic e ects reported a slower structuring of the material with increasing characteristic time, in response to strain change caused by the expansion, results indicate that for high values of relaxation time, higher the distance between the unyielded regions of smaller channel and larger.

O estudo do comportamento reológico de fluidos não Newtonianos tem grande importância em diversas áreas da engenharia. O aumento na demanda destes fluidos por exemplo, no uso doméstico, pessoal e em processos químicos acarreta em dificuldades que vão desde o processo de sua mistura ao seu manuseio. Dentre os fluidos não Newtonianos estão os viscoelásticos, os quais exibem deformação aparente quando os níveis de tensões são inferiores ao limite de escoamento do material e, dentro desta classe, alguns ainda podem apresentam comportamento elástico quando submetidos à baixa taxa de cisalhamento. Juntamente com os efeitos elasto-viscoplásticos, os materiais podem apresentar comportamento tixotrópico, onde devido as tensões sua reestruturação não é instantânea. Na presente Tese, fez-se um estudo numérico a fim de analisar o problema especifico de escoamentos de fluidos elasto-viscoplásticos com comportamento tixotrópico em uma expansão planar abrupta na razão de aspecto 4 : 1. O modelo mecânico aplicado consiste de uma equação viscoelástica para o campo de tensões, uma evolutiva para o parâmetro de estrutura do material, bem como as equações de conservação de massa e momento. O modelo mecânico aplicado mostrou-se capaz de prever o comportamento tixotrópico. A aproximação numérica do modelo aplicado foi feita através do método estabilizado de elementos finitos, especificamente o método Galerkin Mínimos-Quadrados (GLS), o qual foi implementado no código de elementos finitos para fluidos não Newtonianos em desenvolvimento no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) da Universidade Federal do Rio Grande do Sul. Os fenômenos reológicos presentes no problema foram analisados a partir da influência da cinemática do escoamento, da elasticidade e da tixotropia, no nível de estruturação do material, na posição e tamanho das regiões não escoadas e na deformação elástica do material. Os resultados mostraram-se satisfatórios, pois condizem com os apresentados na literatura, apontando a boa predição do modelo mecânico aplicado bem como a robustez de sua implementação computacional. / The rheological behavior of non-Newtonian uids study is of great then importance in many areas of engineering. The increase in demand of these uids - for example, domestic use, personal and processes chemical - causes di culties ranging from the process of mixing it to handling. Among the non-Newtonian uids are viscoelastic, which exhibit apparent deformation when stress levels are lower than the yield limit of the material and, within this class, some still have elastic behavior when subjected to low shear rates. Together with the elastic-viscoplastic e ects, materials may exhibit thixotropic behavior, ie, due the restructuring strain is not instantaneous. In this thesis was made a numerical study to simulate the speci c problem of ow of elastic-viscoplastic uids with thixotropic behavior in abrupt planar expansion { common geometry in industrial systems associated with elastic-viscoplastic uids, whose ratio the aspect is 4:1. The mechanical behavior of most of these structured materials, are highly non-Newtonian, with this, there is need to make them more pseudoplastic, causing undesirable behaviors such as thixotropic which is a phenomenon somewhat characterized and modeled in the literature. The mechanical model applied is able to predict thixotropic behavior and is composed of a viscoelastic equation for the stress eld and an evolutionary to the material structure parameter in addition to the mass and momentum conservation equations. This mechanical model is approximated by a stabilized nite element model, called the Galerkin method of least squares (GLS). In order to study the rheological phenomena present in the problem is analyzed the in uence of the ow kinematics, elasticity and thixotropy in the level of structure of the material, in position and size of unyielded regions and the elastic deformation of the material. The results were satisfactory, since the the study of intensity the U the results agree with those reported in the literature pointing to good prediction of the mechanical model applied well as the hardiness of their computational implementation.The results of the elasticity showed quite spectacular behavior of unyielded regions of the material since for high values of relaxation time, the unyielded region of channel larger has the form of " ngers"because of the high exibility of the material along the line of symmetry of the channel. The results of the analysis of the elastic deformation show, that model correctly dosing the elasticity in the unyielded regions. The thixotropic e ects reported a slower structuring of the material with increasing characteristic time, in response to strain change caused by the expansion, results indicate that for high values of relaxation time, higher the distance between the unyielded regions of smaller channel and larger.

O estudo do comportamento reológico de fluidos não Newtonianos tem grande importância em diversas áreas da engenharia. O aumento na demanda destes fluidos por exemplo, no uso doméstico, pessoal e em processos químicos acarreta em dificuldades que vão desde o processo de sua mistura ao seu manuseio. Dentre os fluidos não Newtonianos estão os viscoelásticos, os quais exibem deformação aparente quando os níveis de tensões são inferiores ao limite de escoamento do material e, dentro desta classe, alguns ainda podem apresentam comportamento elástico quando submetidos à baixa taxa de cisalhamento. Juntamente com os efeitos elasto-viscoplásticos, os materiais podem apresentar comportamento tixotrópico, onde devido as tensões sua reestruturação não é instantânea. Na presente Tese, fez-se um estudo numérico a fim de analisar o problema especifico de escoamentos de fluidos elasto-viscoplásticos com comportamento tixotrópico em uma expansão planar abrupta na razão de aspecto 4 : 1. O modelo mecânico aplicado consiste de uma equação viscoelástica para o campo de tensões, uma evolutiva para o parâmetro de estrutura do material, bem como as equações de conservação de massa e momento. O modelo mecânico aplicado mostrou-se capaz de prever o comportamento tixotrópico. A aproximação numérica do modelo aplicado foi feita através do método estabilizado de elementos finitos, especificamente o método Galerkin Mínimos-Quadrados (GLS), o qual foi implementado no código de elementos finitos para fluidos não Newtonianos em desenvolvimento no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) da Universidade Federal do Rio Grande do Sul. Os fenômenos reológicos presentes no problema foram analisados a partir da influência da cinemática do escoamento, da elasticidade e da tixotropia, no nível de estruturação do material, na posição e tamanho das regiões não escoadas e na deformação elástica do material. Os resultados mostraram-se satisfatórios, pois condizem com os apresentados na literatura, apontando a boa predição do modelo mecânico aplicado bem como a robustez de sua implementação computacional. / The rheological behavior of non-Newtonian uids study is of great then importance in many areas of engineering. The increase in demand of these uids - for example, domestic use, personal and processes chemical - causes di culties ranging from the process of mixing it to handling. Among the non-Newtonian uids are viscoelastic, which exhibit apparent deformation when stress levels are lower than the yield limit of the material and, within this class, some still have elastic behavior when subjected to low shear rates. Together with the elastic-viscoplastic e ects, materials may exhibit thixotropic behavior, ie, due the restructuring strain is not instantaneous. In this thesis was made a numerical study to simulate the speci c problem of ow of elastic-viscoplastic uids with thixotropic behavior in abrupt planar expansion { common geometry in industrial systems associated with elastic-viscoplastic uids, whose ratio the aspect is 4:1. The mechanical behavior of most of these structured materials, are highly non-Newtonian, with this, there is need to make them more pseudoplastic, causing undesirable behaviors such as thixotropic which is a phenomenon somewhat characterized and modeled in the literature. The mechanical model applied is able to predict thixotropic behavior and is composed of a viscoelastic equation for the stress eld and an evolutionary to the material structure parameter in addition to the mass and momentum conservation equations. This mechanical model is approximated by a stabilized nite element model, called the Galerkin method of least squares (GLS). In order to study the rheological phenomena present in the problem is analyzed the in uence of the ow kinematics, elasticity and thixotropy in the level of structure of the material, in position and size of unyielded regions and the elastic deformation of the material. The results were satisfactory, since the the study of intensity the U the results agree with those reported in the literature pointing to good prediction of the mechanical model applied well as the hardiness of their computational implementation.The results of the elasticity showed quite spectacular behavior of unyielded regions of the material since for high values of relaxation time, the unyielded region of channel larger has the form of " ngers"because of the high exibility of the material along the line of symmetry of the channel. The results of the analysis of the elastic deformation show, that model correctly dosing the elasticity in the unyielded regions. The thixotropic e ects reported a slower structuring of the material with increasing characteristic time, in response to strain change caused by the expansion, results indicate that for high values of relaxation time, higher the distance between the unyielded regions of smaller channel and larger.

Kyoto University (京都大学) / 0048 / 新制・課程博士 / 博士(工学) / 甲第16375号 / 工博第3456号 / 新制||工||1523(附属図書館) / 29006 / 京都大学大学院工学研究科社会基盤工学専攻 / (主査)教授 岡 二三生, 教授 勝見 武, 准教授 木元 小百合 / 学位規則第4条第1項該当

Cette thèse vise à proposer des outils innovants pour le calcul de structures aéronautiques évoluant à haute température. En effet, les régimes de fonctionnement des moteurs actuels conduisent à des évolutions élasto-viscoplastiques généralisées dans les pièces métalliques et l’utilisation de modèles simplifiés (élastiques) n’est plus totalement satisfaisante en terme de précision, même en phase de préconception. De même, la géométrie complexe permettant le refroidissement continu des pièces (micro-perforations) doit être prise en compte de manière exacte. Les techniques de calcul standard pour ce genre de problème conduiraient à des simulations lentes et peu flexibles (la moindre modification entraînant une remise en œuvre complète de la chaîne de calcul). Plus précisément, cette thèse étend les méthodes de type global/local non-intrusives au cas de la viscoplasticité généralisée en utilisant deux échelles de temps et d'espace, chacune adaptée aux phénomènes locaux et globaux à capturer. La méthode est ensuite étendue au calcul de nombreux cycles complexes de chargement, par des techniques de saut de cycles. Le schéma de couplage en temps permet alors une adaptation locale du pas de temps par sous-domaine. Des techniques d’accélération de convergence sont proposées, à l’échelle d’un incrément puis à celle de la succession de cycles (sauts de cycles). Ces développements permettent d’obtenir rapidement et précisément une estimation du cycle limite qui alimente un modèle de durée de vie. Le couplage non-intrusif est réalisé dans un script de programmation pilotant un code commercial (dans notre cas le langage Python et Abaqus/Standard). La méthode a été appliquée sur des plateformes de calculs industrielles, en réutilisant directement des maillages et les mises en données issues de modèles intervenant plus tôt dans la chaîne de calcul. Un cas métier, issu d’un bureau d’études de Safran Aircraft Engines, a pu être traité. / This thesis consists in developing innovating tools destined to the simulation of aeronautical structures evolving at high temperature. Indeed, working rates of current engines lead to an elasto-viscoplastic evolution generalized in metallic parts and the use of simplified models (linear elastic) are no longer totally satisfying in term of accuracy, even in initial design process. Likewise, the complex geometry allowing the continuous cool down process of parts (micro-perforations) has to be exactly taken into account. The standard computation techniques dedicated to this kind of models would lead to slow simulations with a lack of flexibility (the slightest modifications leading to restart the whole design process of the computation chain).More precisely, this thesis extends the noninvasive global/local methods to the framework of viscoplasticity generalized to the whole structure, using two scales in time and space, each one adapted to global and local phenomena to capture. The method is then extended to the computation of high number of complex load cycles, by skipped cycles techniques. The time coupling scheme lets then a local adaptation of time steps per subdomain. Convergence acceleration techniques are also set up, first for one time step and then through several load cycles (skipped cycles). These developments conduct to obtain quickly an evaluation of the limit cycle providing data to a lifetime expectancy model.The noninvasive coupling is realized in a programming language script managing the commercial software (respectively in our case Python and Abaqus/Standard). The method has been applied on industrial computational platforms, by reusing directly meshes and data from previous engineering tasks appearing earlier in the computational chain. A genuine test case from a Safran Aircraft Engines design office, was performed successfully.

[pt] Sendo a exploração de petróleo cada vez mais profunda, atingindo vários quilômetros abaixo dos leitos marítimos, uma simulação numérica adequada para prever o comportamento de rochas antes da perfuração está se tornando cada vez mais importante para a indústria de petróleo, devido aos altos custos operacionais para garantir a estabilidade do poço. Um dos problemas típicos que os engenheiros de petróleo enfrentam é a instabilidade de poços em rochas evaporíticas, que geralmente são sujeitas a deformações excessivas e possuem características de comportamento mecânico dependentes do tempo. Grandes reservas de petróleo que foram encontradas ao redor do mundo em rochas evaporíticas encorajaram aos engenheiros a investigar mais profundamente essa resposta mecânica dependente do tempo (creep). Por outro lado, a concentração de tensões em torno dos furos pode torná-lo um projeto de engenharia inviável. A desestabilização devido a grandes deformações pode levar a problemas como colapso e fechamento do poço. Assim, para estudar a estabilidade de poços e projetar um sistema de produção de petróleo adequado, a influência do creep deve ser considerada adequadamente no comportamento esperado do material. Nesta pesquisa, as equações que definem o modelo elásto-viscoplástico de Sterpi e Gioda (2007) são implementadas no programa computacional FLAC 3D - ITASCA, código baseado no método das diferenças finitas. O procedimento de implementação consiste em gerar uma DLL (Dynamic Link Library) escrita em linguagem de programação C (mais mais) e integrar as tensões ao longo de um incremento de deformações usando um algoritmo de integração explícito. O modelo elásto-viscoplástico foi validado com resultados experimentais envolvendo rochas salinas disponíveis na literatura. Finalmente, é realizada uma simulação numérica com o programa FLAC 3D, considerando o modelo de Sterpi e Gioda, com o objetivo de estudar o fechamento de um poço no pré-sal brasileiro. / [en] Since oil exploration is reaching several kilometers below seafloor, an adequate numerical simulation for prediction of rock behavior prior to drilling is becoming increasingly important for the oil industry due to the high operating costs to ensure well stability. One of the typical problems that oil engineers face is the instability of wells in evaporitic rocks, which are often subject to excessive deformation and have time-dependent mechanical behavior characteristics. Huge oil reserves that were found around the world in evaporitic reservoirs encouraged engineers to further investigate this mechanical time-dependent response (creep). On the other hand, the concentration of stresses around the well can make it an unviable engineering project. Destabilization due to large deformations may be involved, leading to problems such as collapse and well closure. Thus, in order to study the stability of the boreholes and design an appropriate oil production system, creep must be properly considered in the expected behavior of the material. In this research, the equations that define the elasto-viscoplastic model of Sterpi and Gioda (2007) will be implemented in the computer program FLAC 3D - Itasca, code based on the finite difference method. The implementation procedure consists of generating a Dynamic Link Library (DLL) written in C (plus plus) programming language and integrating the stress along a strain increment by using an explicit integration algorithm. The elasto-viscoplastic model has been validated with experimental results from triaxial creep tests in salt rocks available in the literature, indicating a successful implementation. Finally, a numerical simulation in FLAC 3D considering the model of Sterpi and Gioda is performed in order to study a pre-salt wellbore closure.