Neste trabalho é apresentado a modelagem mecânica e numérica de um escoamento elastoviscoplástico tixotrópico, em termos dos campos de velocidade, pressão, tensão e parâmetro de estrutura. A discretização numérica é feita pelo Método de Elementos Finitos Estabilizado, também conhecido como Galerkin Mínimos Quadrados (GMQ), através de elementos quadrangulares bilineares. O clássico problema da cavidade é utilizado nas simulações, a fim de comparar a formulação e o código utilizados com os resultados conhecidos na literatura. Além disso, apresenta-se o estudo de materiais elasto-viscoplástico tixotrópico em uma contração abrupta na escala 4:1, utilizando a formulação descrita e uma nova função viscosidade para fluidos viscoplásticos, denominada função Viscoplástica Harmônica (VPH). Resultados envolvendo a função VPH são introduzidos e discutidos pela primeira vez nesta tese e apresenta um ótimo ajuste de curva, quando comparada com outras funções disponíveis na literatura. Além da fácil implementação, essa função também apresenta um platô para as altas e baixas viscosidades que são fisicamente realistas, visto que não é possível uma viscosidade infinita ou nula. O menor tempo computacional é também uma característica perceptível nas simulações usando a nova função viscoplástica, isso é um atributo do seu equacionamento que não depende de um termo exponencial, como outros modelos. O estudo de qualidade de malha também é apresentado a fim de garantir a escolha do domínio discreto adequado. Apesar do uso de elementos de ordem inferior, o método GMQ mostrou-se estável na aproximação numérica de todos os problemas dispostos, garantindo até mesmo a análise sobre os efeitos da cinemática, da elasticidade e da tixotropia no escoamento dos fluidos dentro da contração abrupta. / In this work the mechanical and numerical modeling of a thixotropic elasto-viscoplastic flow in terms of the velocity, pressure, stress and structure parameter is presented. Numerical discretization is done by the Stabilized Finite Element Method, also known as Galerkin Least Squares (GLS), through bilinear quadrangular elements. The classical liddriven cavity problem is used in the simulations in order to compare the formulation and code used with the results in the literature. In addition, the study of elasto-viscoplastic thixotropic materials in a 4:1 abrupt contraction using the described formulation and a new viscosity function for viscoplastic fluids, called the Harmonic Viscoplastic Function (HVP), is presented. Results involving the HVP function are introduced and discussed for the first time in this thesis and present a better curve fit when compared to other functions available in the literature. Besides to easy implementation, this function also features a plateau for high and low viscosities that are physically realistic, since infinite or zero viscosity is not possible. The shortest computational time is also a perceptible feature in the simulations using the new viscoplastic function, this is an attribute of its equation that does not depend on an exponential term like other models. The mesh quality study is also presented in order to ensure the choice of the appropriate discrete domain. Despite the use of lower-order elements, the GLS method proved to be stable in the numerical approximation of all the problems, guaranteeing even the analysis of the effects of kinematics, elasticity and thixotropy on fluid flow within the abrupt contraction.

Neste trabalho é apresentado a modelagem mecânica e numérica de um escoamento elastoviscoplástico tixotrópico, em termos dos campos de velocidade, pressão, tensão e parâmetro de estrutura. A discretização numérica é feita pelo Método de Elementos Finitos Estabilizado, também conhecido como Galerkin Mínimos Quadrados (GMQ), através de elementos quadrangulares bilineares. O clássico problema da cavidade é utilizado nas simulações, a fim de comparar a formulação e o código utilizados com os resultados conhecidos na literatura. Além disso, apresenta-se o estudo de materiais elasto-viscoplástico tixotrópico em uma contração abrupta na escala 4:1, utilizando a formulação descrita e uma nova função viscosidade para fluidos viscoplásticos, denominada função Viscoplástica Harmônica (VPH). Resultados envolvendo a função VPH são introduzidos e discutidos pela primeira vez nesta tese e apresenta um ótimo ajuste de curva, quando comparada com outras funções disponíveis na literatura. Além da fácil implementação, essa função também apresenta um platô para as altas e baixas viscosidades que são fisicamente realistas, visto que não é possível uma viscosidade infinita ou nula. O menor tempo computacional é também uma característica perceptível nas simulações usando a nova função viscoplástica, isso é um atributo do seu equacionamento que não depende de um termo exponencial, como outros modelos. O estudo de qualidade de malha também é apresentado a fim de garantir a escolha do domínio discreto adequado. Apesar do uso de elementos de ordem inferior, o método GMQ mostrou-se estável na aproximação numérica de todos os problemas dispostos, garantindo até mesmo a análise sobre os efeitos da cinemática, da elasticidade e da tixotropia no escoamento dos fluidos dentro da contração abrupta. / In this work the mechanical and numerical modeling of a thixotropic elasto-viscoplastic flow in terms of the velocity, pressure, stress and structure parameter is presented. Numerical discretization is done by the Stabilized Finite Element Method, also known as Galerkin Least Squares (GLS), through bilinear quadrangular elements. The classical liddriven cavity problem is used in the simulations in order to compare the formulation and code used with the results in the literature. In addition, the study of elasto-viscoplastic thixotropic materials in a 4:1 abrupt contraction using the described formulation and a new viscosity function for viscoplastic fluids, called the Harmonic Viscoplastic Function (HVP), is presented. Results involving the HVP function are introduced and discussed for the first time in this thesis and present a better curve fit when compared to other functions available in the literature. Besides to easy implementation, this function also features a plateau for high and low viscosities that are physically realistic, since infinite or zero viscosity is not possible. The shortest computational time is also a perceptible feature in the simulations using the new viscoplastic function, this is an attribute of its equation that does not depend on an exponential term like other models. The mesh quality study is also presented in order to ensure the choice of the appropriate discrete domain. Despite the use of lower-order elements, the GLS method proved to be stable in the numerical approximation of all the problems, guaranteeing even the analysis of the effects of kinematics, elasticity and thixotropy on fluid flow within the abrupt contraction.

A maioria dos líquidos encontrados na natureza são não newtonianos e o estudo do seu comportamento reológico tem uma importância significante em diferentes áreas da engenharia. Entre eles, existe uma classe de fluidos que exibem pequena deformação aparente quando sujeitos a um nível de tensões inferior a uma tensão limite de escoamento, referido como comportamento viscoplástico. Nesta classe de materiais, alguns apresentam também comportamento elástico quando submetidos a baixas taxas de cisalhamento. A presente Tese tem como objetivo o estudo numérico de escoamentos bidimensionais em regime permanente de fluidos elasto-viscoplásticos através de uma expansão-contração planar. O modelo mecânico é definido pelas equações de conservação de massa e de balanço de momentum acopladas ao modelo elasto-viscoplástico proposto nesta Tese. Esta modelagem é aproximada por um método de elementos finitos multi-campos estabilizado baseado na metodologia de Galerkin mínimos-quadrados que possui como variáveis primais os campos de tensão extra polimérica, velocidade e pressão. As condições de compatibilidade entre os sub-espaços de elementos finitos para tensão extra-velocidade e velocidade-pressão são violadas, permitindo assim a utilização de interpolações de igual ordem. O método estabilizado 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. Nesta Tese é adotada uma metodologia alternativa para a definição das zonas rígidas do escoamento como sendo a posição onde a taxa de cisalhamento é igual a um valor dado pela relação de parâmetros reológicos do fluido, especificamente a tensão limite de escoamento e a viscosidade newtoniana para baixas taxas de cisalhamento. Nas simulações numéricas realizadas, o tempo de relaxação adimensional, o número de salto, o coeficiente power-law, a vazão adimensional e a massa específica adimensional são variados de forma a avaliar de que modo influenciam na dinâmica de escoamentos elastoviscoplásticos. Os resultados obtidos estão qualitativamente de acordo com a literatura, atestando a estabilidade da formulação empregada. / Non-Newtonian fluids are the majority of liquids found in nature and the study of their rheological behavior has a significant importance on different areas of engineering. Among them, there is a class of materials that exhibits little apparent deformation when subjected to a stress level behind an yield stress, referenced as viscoplastic material. In this class of materials, some fluids also exhibit elastic behavior at low shear rates. The present work aimed to a numerical study of two-dimensional steady state laminar flows of elasto-viscoplastic fluids through a planar expansion-contraction cavity. The mechanical model was defined by the mass conservation and momentum balance equations coupled to the elasto-viscoplastic model porposed in this work. This modeling has been approximated by a stabilized multi-field finite element method based on the Galerkin least-squares methodology, having as primal variables the elastic extra-stress component, velocity and pressure fields. In this way, the compatibility conditions between the extra-stress-velocity and pressure-velocity (Babuška- Brezzi condition) finite element subspaces are violated, allowing to use equal-order finite element interpolations. The stabilized method has been implemented in the finite element code for non-Newtonian fluids under development at the Laboratory of Applied and Computational Fluid Mechanics (LAMAC) of the Federal University of Rio Grande do Sul. An alternative methodology is adopted to define the yield surface as the position where the strain rate is equal to a value given by the relation of the rheological parameters of the fluid, namely the yield stress and the viscosity at low shear rates. In the performed numerical simulations, the non-dimensional relaxation time, the jump number, the power-law coefficient, the non-dimensional flow rate, and the non-dimensional density are varied in order to evaluate their influence on the elasto-viscoplastic fluid dynamics. All results found are in qualitatively accordance with the affine literature, and attesting the good stability features of the formulation.

A maioria dos líquidos encontrados na natureza são não newtonianos e o estudo do seu comportamento reológico tem uma importância significante em diferentes áreas da engenharia. Entre eles, existe uma classe de fluidos que exibem pequena deformação aparente quando sujeitos a um nível de tensões inferior a uma tensão limite de escoamento, referido como comportamento viscoplástico. Nesta classe de materiais, alguns apresentam também comportamento elástico quando submetidos a baixas taxas de cisalhamento. A presente Tese tem como objetivo o estudo numérico de escoamentos bidimensionais em regime permanente de fluidos elasto-viscoplásticos através de uma expansão-contração planar. O modelo mecânico é definido pelas equações de conservação de massa e de balanço de momentum acopladas ao modelo elasto-viscoplástico proposto nesta Tese. Esta modelagem é aproximada por um método de elementos finitos multi-campos estabilizado baseado na metodologia de Galerkin mínimos-quadrados que possui como variáveis primais os campos de tensão extra polimérica, velocidade e pressão. As condições de compatibilidade entre os sub-espaços de elementos finitos para tensão extra-velocidade e velocidade-pressão são violadas, permitindo assim a utilização de interpolações de igual ordem. O método estabilizado 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. Nesta Tese é adotada uma metodologia alternativa para a definição das zonas rígidas do escoamento como sendo a posição onde a taxa de cisalhamento é igual a um valor dado pela relação de parâmetros reológicos do fluido, especificamente a tensão limite de escoamento e a viscosidade newtoniana para baixas taxas de cisalhamento. Nas simulações numéricas realizadas, o tempo de relaxação adimensional, o número de salto, o coeficiente power-law, a vazão adimensional e a massa específica adimensional são variados de forma a avaliar de que modo influenciam na dinâmica de escoamentos elastoviscoplásticos. Os resultados obtidos estão qualitativamente de acordo com a literatura, atestando a estabilidade da formulação empregada. / Non-Newtonian fluids are the majority of liquids found in nature and the study of their rheological behavior has a significant importance on different areas of engineering. Among them, there is a class of materials that exhibits little apparent deformation when subjected to a stress level behind an yield stress, referenced as viscoplastic material. In this class of materials, some fluids also exhibit elastic behavior at low shear rates. The present work aimed to a numerical study of two-dimensional steady state laminar flows of elasto-viscoplastic fluids through a planar expansion-contraction cavity. The mechanical model was defined by the mass conservation and momentum balance equations coupled to the elasto-viscoplastic model porposed in this work. This modeling has been approximated by a stabilized multi-field finite element method based on the Galerkin least-squares methodology, having as primal variables the elastic extra-stress component, velocity and pressure fields. In this way, the compatibility conditions between the extra-stress-velocity and pressure-velocity (Babuška- Brezzi condition) finite element subspaces are violated, allowing to use equal-order finite element interpolations. The stabilized method has been implemented in the finite element code for non-Newtonian fluids under development at the Laboratory of Applied and Computational Fluid Mechanics (LAMAC) of the Federal University of Rio Grande do Sul. An alternative methodology is adopted to define the yield surface as the position where the strain rate is equal to a value given by the relation of the rheological parameters of the fluid, namely the yield stress and the viscosity at low shear rates. In the performed numerical simulations, the non-dimensional relaxation time, the jump number, the power-law coefficient, the non-dimensional flow rate, and the non-dimensional density are varied in order to evaluate their influence on the elasto-viscoplastic fluid dynamics. All results found are in qualitatively accordance with the affine literature, and attesting the good stability features of the formulation.

A maioria dos líquidos encontrados na natureza são não newtonianos e o estudo do seu comportamento reológico tem uma importância significante em diferentes áreas da engenharia. Entre eles, existe uma classe de fluidos que exibem pequena deformação aparente quando sujeitos a um nível de tensões inferior a uma tensão limite de escoamento, referido como comportamento viscoplástico. Nesta classe de materiais, alguns apresentam também comportamento elástico quando submetidos a baixas taxas de cisalhamento. A presente Tese tem como objetivo o estudo numérico de escoamentos bidimensionais em regime permanente de fluidos elasto-viscoplásticos através de uma expansão-contração planar. O modelo mecânico é definido pelas equações de conservação de massa e de balanço de momentum acopladas ao modelo elasto-viscoplástico proposto nesta Tese. Esta modelagem é aproximada por um método de elementos finitos multi-campos estabilizado baseado na metodologia de Galerkin mínimos-quadrados que possui como variáveis primais os campos de tensão extra polimérica, velocidade e pressão. As condições de compatibilidade entre os sub-espaços de elementos finitos para tensão extra-velocidade e velocidade-pressão são violadas, permitindo assim a utilização de interpolações de igual ordem. O método estabilizado 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. Nesta Tese é adotada uma metodologia alternativa para a definição das zonas rígidas do escoamento como sendo a posição onde a taxa de cisalhamento é igual a um valor dado pela relação de parâmetros reológicos do fluido, especificamente a tensão limite de escoamento e a viscosidade newtoniana para baixas taxas de cisalhamento. Nas simulações numéricas realizadas, o tempo de relaxação adimensional, o número de salto, o coeficiente power-law, a vazão adimensional e a massa específica adimensional são variados de forma a avaliar de que modo influenciam na dinâmica de escoamentos elastoviscoplásticos. Os resultados obtidos estão qualitativamente de acordo com a literatura, atestando a estabilidade da formulação empregada. / Non-Newtonian fluids are the majority of liquids found in nature and the study of their rheological behavior has a significant importance on different areas of engineering. Among them, there is a class of materials that exhibits little apparent deformation when subjected to a stress level behind an yield stress, referenced as viscoplastic material. In this class of materials, some fluids also exhibit elastic behavior at low shear rates. The present work aimed to a numerical study of two-dimensional steady state laminar flows of elasto-viscoplastic fluids through a planar expansion-contraction cavity. The mechanical model was defined by the mass conservation and momentum balance equations coupled to the elasto-viscoplastic model porposed in this work. This modeling has been approximated by a stabilized multi-field finite element method based on the Galerkin least-squares methodology, having as primal variables the elastic extra-stress component, velocity and pressure fields. In this way, the compatibility conditions between the extra-stress-velocity and pressure-velocity (Babuška- Brezzi condition) finite element subspaces are violated, allowing to use equal-order finite element interpolations. The stabilized method has been implemented in the finite element code for non-Newtonian fluids under development at the Laboratory of Applied and Computational Fluid Mechanics (LAMAC) of the Federal University of Rio Grande do Sul. An alternative methodology is adopted to define the yield surface as the position where the strain rate is equal to a value given by the relation of the rheological parameters of the fluid, namely the yield stress and the viscosity at low shear rates. In the performed numerical simulations, the non-dimensional relaxation time, the jump number, the power-law coefficient, the non-dimensional flow rate, and the non-dimensional density are varied in order to evaluate their influence on the elasto-viscoplastic fluid dynamics. All results found are in qualitatively accordance with the affine literature, and attesting the good stability features of the formulation.

Escoamentos sem inércia de fluido elasto-viscoplástico, dentro de uma cavidade, são numericamente analisados. As soluções visam compreender a influência dos efeitos viscosos e elásticos na topologia de superfícies de escoamento. Assumindo-se que o colapso da microestrutura do material é instantâneo, o modelo mecânico é constituído pelas equações governantes de massa e momentum para fluidos incompressíveis, associado a uma equação hiperbólica para o tensor tensão extra, baseado na equação do modelo Oldroyd-B (Nassar et al, 2011). A principal característica do modelo é considerar a viscosidade e o tempo de relaxação como função da taxa de deformação, permitindo a pseudoplasticidade de viscosidade e restringindo os efeitos elásticos para as regiões não deformadas do material. As simulações numéricas são realizadas através do método de Galerkin mínimos quadrados a três campos: tensor tensão extra, pressão e velocidade. Os resultados mostram que as superfícies de escoamento do material são fortemente influenciadas pela ação combinada entre os efeitos elásticos e viscosos, estando em conformidade com a recente visualização experimental dos fluxos elasto-viscoplásticos. / Elasto viscoplastic uid ows without inertia, within a cavity, are numerically analyzed. The solutions aim to understand the in uence of viscous and elastic e ects on the topology of yield surfaces. Assuming that the collapse of the material microstructure is instantaneous, the mechanical model consists of the governing equations of mass and momentum for incompressible uids, associated with a hyperbolic equation for the extra stress tensor, based on the equation of the Oldroyd-B model (Nassar et al, 2011). The main feature of model is to consider the viscosity and the relaxation time as a function of shear rate, allowing the shear-thinning of viscosity and restricting elastic e ects for regions not deformed material. Numerical simulations are performed by the method of Galerkin Least Squares to three elds: extra stress tensor, pressure and velocity. The results show that the yield surfaces of material are strongly in uenced by the combined action between the elastic and viscous e ects, complying with the recent experimental visualization of elasto-viscoplastic ows.

Escoamentos sem inércia de fluido elasto-viscoplástico, dentro de uma cavidade, são numericamente analisados. As soluções visam compreender a influência dos efeitos viscosos e elásticos na topologia de superfícies de escoamento. Assumindo-se que o colapso da microestrutura do material é instantâneo, o modelo mecânico é constituído pelas equações governantes de massa e momentum para fluidos incompressíveis, associado a uma equação hiperbólica para o tensor tensão extra, baseado na equação do modelo Oldroyd-B (Nassar et al, 2011). A principal característica do modelo é considerar a viscosidade e o tempo de relaxação como função da taxa de deformação, permitindo a pseudoplasticidade de viscosidade e restringindo os efeitos elásticos para as regiões não deformadas do material. As simulações numéricas são realizadas através do método de Galerkin mínimos quadrados a três campos: tensor tensão extra, pressão e velocidade. Os resultados mostram que as superfícies de escoamento do material são fortemente influenciadas pela ação combinada entre os efeitos elásticos e viscosos, estando em conformidade com a recente visualização experimental dos fluxos elasto-viscoplásticos. / Elasto viscoplastic uid ows without inertia, within a cavity, are numerically analyzed. The solutions aim to understand the in uence of viscous and elastic e ects on the topology of yield surfaces. Assuming that the collapse of the material microstructure is instantaneous, the mechanical model consists of the governing equations of mass and momentum for incompressible uids, associated with a hyperbolic equation for the extra stress tensor, based on the equation of the Oldroyd-B model (Nassar et al, 2011). The main feature of model is to consider the viscosity and the relaxation time as a function of shear rate, allowing the shear-thinning of viscosity and restricting elastic e ects for regions not deformed material. Numerical simulations are performed by the method of Galerkin Least Squares to three elds: extra stress tensor, pressure and velocity. The results show that the yield surfaces of material are strongly in uenced by the combined action between the elastic and viscous e ects, complying with the recent experimental visualization of elasto-viscoplastic ows.

Escoamentos sem inércia de fluido elasto-viscoplástico, dentro de uma cavidade, são numericamente analisados. As soluções visam compreender a influência dos efeitos viscosos e elásticos na topologia de superfícies de escoamento. Assumindo-se que o colapso da microestrutura do material é instantâneo, o modelo mecânico é constituído pelas equações governantes de massa e momentum para fluidos incompressíveis, associado a uma equação hiperbólica para o tensor tensão extra, baseado na equação do modelo Oldroyd-B (Nassar et al, 2011). A principal característica do modelo é considerar a viscosidade e o tempo de relaxação como função da taxa de deformação, permitindo a pseudoplasticidade de viscosidade e restringindo os efeitos elásticos para as regiões não deformadas do material. As simulações numéricas são realizadas através do método de Galerkin mínimos quadrados a três campos: tensor tensão extra, pressão e velocidade. Os resultados mostram que as superfícies de escoamento do material são fortemente influenciadas pela ação combinada entre os efeitos elásticos e viscosos, estando em conformidade com a recente visualização experimental dos fluxos elasto-viscoplásticos. / Elasto viscoplastic uid ows without inertia, within a cavity, are numerically analyzed. The solutions aim to understand the in uence of viscous and elastic e ects on the topology of yield surfaces. Assuming that the collapse of the material microstructure is instantaneous, the mechanical model consists of the governing equations of mass and momentum for incompressible uids, associated with a hyperbolic equation for the extra stress tensor, based on the equation of the Oldroyd-B model (Nassar et al, 2011). The main feature of model is to consider the viscosity and the relaxation time as a function of shear rate, allowing the shear-thinning of viscosity and restricting elastic e ects for regions not deformed material. Numerical simulations are performed by the method of Galerkin Least Squares to three elds: extra stress tensor, pressure and velocity. The results show that the yield surfaces of material are strongly in uenced by the combined action between the elastic and viscous e ects, complying with the recent experimental visualization of elasto-viscoplastic ows.

Le présent travail est consacré à l'étude du comportement thermo-hydro-mécanique linéaire et non linéaire des roches poreuses de type argilites par approche de changement d'échelle. A partir des observations microstructurales de ces matériaux, un modèle conceptuel a été proposé. Dans ce modèle, le volume élémentaire représentatif du milieu hétérogène est composé d'une phase matricielle argileuse contenant des inclusions sphériques de minéraux de quartz et de calcite et des inclusions ellipsoïdales aplaties représentant l'espace poreux. Dans un premier temps, le procédé de la modélisation a été exploité par la détermination des propriétés effectives isotropes et isotrope transverses des argilites : la conductivité thermique et les propriétés thermo-hydro-mécanique croisées. En outre, de nombreuses études numériques ont mis en évidence l'influence de la morphologie de l'espace poreux, de la minéralogie et des schémas d'estimation sur les résultats prédictifs. Dans un deuxième temps, nous avons modélisé le comportement mécanique non linéaire (élasto-plastique, élasto-viscoplastique) des roches argileuses. La comparaison entre les simulations numériques et les résultats expérimentaux disponibles (essai de compression triaxiale, essai de fluage) a confirmé la validation du modèle développé / The present work deals with the linear and non-linear thermo-hydro-mechanical behaviour of porous rocks such as the argillite by the multiscale modelling approach. Based on microstructure observations, a conceptual model was proposed. In this model, the representative elementary volume of a heterogeneous medium is composed of an argillaceous matrix containing spherical inclusions of minerals quartz and calcite and ellipsoidal inclusions representing the pore space. In a first step, the process of modelling has been exploited by determining the isotropic and transversely isotropic effective properties of the argillite: thermal conductivity and thermo-hydro-mechanical properties. Furthermore, many numerical studies have highlighted the influence of the morphology of the pore space, of the mineralogy and of the estimate schemes to the predictive results. In a second step, we modelled the non linear mechanical behaviour (elasto-plastic, elasto-viscoplastic) of argillaceous rocks. The comparison between numerical simulations and available experimental results (triaxial compression test, creep test) confirmed the validation of the model developed

A major instrumented geosynthetic reinforced approach embankment was constructed to 5.5 m elevation above ground, with prefabricated vertical drains, over a soft compressible clay deposit at Leneghan, Newcastle, Australia in May 1995. The field monitoring of settlements for over six years shows that the embankment manifests significant creep. The instrumentation, field performance and the finite element analyses for predicting the long-term performance of this embankment are described in this thesis. The maximum settlement of 1.1 m was observed one year after the completion of construction. However, the embankment continued to settle at a rate of 0.4 mm/day for the next 5 years. The horizontal displacements of 0.09-0.14 m at various locations and the maximum reinforcement strains of 0.67% were recorded. A numerical model was developed to perform a fully coupled large deformation elasto-viscoplastic finite element analysis for this performance prediction based on creep model proposed by Kutter and Sathialingam (1992). The foundation soil was modelled with creep material behaviour using six noded linear strain triangular elements. A well-documented case history ??? Sackville embankment, New Brunswick, Canada was analysed using this model as a benchmark problem and the model was found to predict all the behaviour characteristics reasonably well. The results obtained from finite element analysis using this model are shown to be in reasonable agreement with the observed performance of Leneghans embankment in terms of settlements, horizontal displacements, excess pore pressures and geosynthetic strains. But, the prediction of settlements was less than satisfactory beyond April 1999. Finite element analyses were performed to study the sensitivity of this embankment behaviour on the variation of hydraulic conductivity values and geosynthetic reinforcement properties. This sensitivity study indicated that the kv variation, the kh/kv ratio and the nominal values of geosynthetic properties adopted in the benchmark analysis are reasonable enough for the long-term behaviour prediction.