Neste trabalho, propõe-se uma formulação não linear baseada no método dos elementos de contorno, para representação de domínios poro-elasto-plásticos reforçados. Esta formulação é apresentada para os casos saturado e não saturado. Para o problema poroso enrijecido um acoplamento com o método dos elementos finitos é empregado, e a técnica de mínimos quadrados permite a regularização dos deslocamentos e do vetor de forças de superfície ao longo das interfaces de acoplamento. São empregadas expressões analíticas para o tratamento das integrais de contorno e de domínio presentes na formulação do método dos elementos de contorno. A formulação de Biot é empregada para a descrição de meios porosos saturados e uma formulação energética baseada nos trabalhos de Coussy é adaptada para a extensão ao caso não saturado. Neste caso, a pressão capilar e energia das interfaces são levadas em consideração. O nível de saturação é descrito pelo modelo de Van Genuchten e o comportamento do esqueleto é descrito ou pelo modelo de Drucker-Prager ou pelo modelo de Cam-Clay modificado. O problema não linear obtido por uma descrição temporal associada a discretização espacial é resolvido pelo método de Newton-Raphson. No caso saturado, o operador tangente consistente é definido e utilizado para obtenção da solução do sistema. Exemplos numéricos são apresentados para validar a formulação proposta. / In this work a nonlinear formulation of the boundary element method (BEM) is proposed to deal with saturated and unsaturated poro-elasto-plastic 2D reinforced domains. To model reinforced porous domains a BEM/FEM (Finite Element Method) modified coupling technique is employed. The coupling is made by using the least square method to regularize the displacement and traction distributions along the interfaces. Analytical expressions have been derived for all boundary and domain integrals required for the formulation. The Biot formulation is used for the description of the saturated porous environments and an energetic consistent formulation based on work of Coussy is adopted for its extension to the framework of unsaturated porous media. In this case, the capillar pressure and the interface energy are taken into account. The Van Genuchten model is used for the determination of saturation level in non-saturated poro-elasto-plastic problems. The Drucker-Prager modified model if used for the saturated poro-elasto-plastic problems and the modified Cam-Clay model for the representation of non-saturated poro-elasto-plastic problems. For the saturated case, the consistent tangent operator is derived and employed inside a Newton procedure to solve non-linear problems. Numerical solutions are presented to validate the proposed models.

Analisam-se problemas de interação solo-estrutura através de uma formulação tridimensional obtida da combinação entre o Método dos Elementos de Contorno (MEC) e o Método dos Elementos Finitos (MEF). Os elementos estruturais que interagem com o solo são modelados pelo Método dos Elementos Finitos. E o solo, considerado como um meio semi-infinito, homogêneo, elástico linear e isotrópico, é modelado pelo Método dos Elementos de Contorno, empregando a solução fundamental de Mindlin. A solução fundamental de Mindlin é particularmente adequada para o tipo de problemas em análise, ou seja, problemas envolvendo sólidos tridimensionais semi-infinitos, já que é necessário discretizar apenas a superfície carregada do solo e/ou a linha de carga, e não todo o sólido tridimensional. A discretização da estaca em vários elementos finitos de viga tridimensional, permitindo a consideração de estacas de qualquer tamanho e submetidas a qualquer tipo de carregamento é uma das principais contribuições do trabalho. Outra contribuição diz respeito à consideração da ação horizontal no sistema placa-estaca-solo, diretamente aplicada na placa, e não redistribuída no topo das estacas. Por isso, os elementos estruturais laminares (radiers, sapatas e blocos de fundação) são modelados por elementos finitos de casca plana, possibilitando a consideração dos efeitos de flexão e de membrana. Essa abordagem permite a análise tanto de um grupo de estacas com bloco de capeamento rígido quanto de um radier estaqueado de qualquer rigidez. O acoplamento entre o MEC e o MEF é feito através de uma formulação mista em que a matriz dos coeficientes de influência do solo obtida pelo MEC é adicionada à matriz de rigidez dos elementos estruturais obtida pelo MEF, resultando em uma matriz de rigidez equivalente. Exemplos numéricos de interação estaca-solo, placa-solo e casca plana-estaca-solo são resolvidos para verificar, validar e demonstrar a eficiência das formulações desenvolvidas e implementadas. / Soil-structure interaction problems are analyzed by a tridimensional formulation obtained combining the Boundary Element Method (BEM) and the Finite Element Method (FEM). Structural elements that interact with the soil are modeled by the Finite Element Method. And the soil, considered as a semi-infinite, homogeneous, linear elastic and isotropic medium, is modeled by the Boundary Element Method, using Mindlin\'s fundamental solution. Mindlin\'s fundamental solution is particularly suitable for the type of problems under analysis, viz. problems involving semi-infinite three-dimensional solids, since it is necessary to discretize only the loaded surface of the soil and/or the line-load, and not all the three-dimensional solid. The discretization of the pile in several three-dimensional beam finite elements, allowing the consideration of piles of any size and subjected to any type of loading, is one of the main contributions of this work. Another contribution is about the consideration of the horizontal load in the plate-pile-soil system, directly applied to the plate, and not redistributed on the top of the piles. Therefore, the laminar structural elements (rafts, footings, and foundation blocks) are modeled by flat shell finite elements, making possible the consideration of the effects of flexion and membrane. This approach allows the analysis of both a capped pile group and a piled raft of any stiffness. The coupling between the BEM and the FEM is performed through a mixed formulation in which the matrix of the soil\'s influence coefficients obtained by the BEM is added to the stiffness matrix of the structural elements obtained by the FEM, resulting in an equivalent stiffness matrix. Numerical examples of pile-soil, plate-soil and flat shell-pile-soil interaction are solved to verify, validate and demonstrate the efficiency of the developed and implemented formulations.

Cette thèse s'articule autour de deux thématiques : les méthodes numériques pour la propagation d'ondes acoustiques sous écoulement et les méthodes de réduction de modèles. Dans la première thématique, nous développons une méthode de couplage d'éléments finis et d'éléments de frontière pour résoudre l'équation d'Helmholtz convectée, lorsque l'écoulement est uniforme à l'extérieur d'un domaine borné. En particulier, nous proposons une formulation bien posée à toutes les fréquences de la source. Dans la deuxième thématique, nous proposons une solution au problème classique d'accumulation d'arrondis machine qui survient en calculant l'estimateur d'erreur a posteriori dans la méthode des bases réduites. Par ailleurs, nous proposons une méthode non intrusive pour calculer une approximation sous forme séparée des systèmes linéaires résultant de l'approximation en dimension finie de problèmes aux limites dépendant d'un ou plusieurs paramètres / This thesis has two topics : numerical methods for acoustic wave propagation in a flow and reduced order models. In the first topic, we develop a coupled finite element and boundary element method to solve the convected Helmholtz equation, when the flow is uniform outside a bounded domain. In particular, we propose a formulation that is well-posed at all the frequencies of the source. In the second topic, we propose a solution to the classical problem of round-off error accumulation that occurs when computing the a posteriori error bound in the reduced basis method. Furthermore, we propose a non intrusive method for the approximation, in a separated representation form, of linear systems resulting from the finite-dimensional approximation of boundary-value problems depending on one or several parameters

Τα δομικά στοιχεία των μηχανολογικών κατασκευών υπόκεινται σε επαναλαμβανόμενες κυκλικές καταπονήσεις, από τις οποίες δημιουργούνται και διαδίδονται ρωγμές. Οι καταπονήσεις αυτές, οι οποίες προκαλούν κόπωση στις κατασκευές, μπορεί να είναι είτε καθαρά μηχανικές είτε θερμικές ή να προκύπτουν σα συνδυασμός θερμικής και μηχανικής φόρτισης. Τυπικές περιπτώσεις θερμικών και θερμομηχανικών φορτίσεων εμφανίζονται σε κατασκευές, όπως σωλήνες κυκλωμάτων ψύξης, πιεστικά δοχεία, συνιστώσες ηλεκτρικών κυκλωμάτων, θάλαμοι μηχανών εσωτερικής καύσης και πτερύγια στροβιλοκινητήρων. Η κυκλική μεταβολή του θερμικού φορτίου στις προαναφερθείσες περιπτώσεις, συνιστά συνθήκες θερμικής κόπωσης. Επίσης, λόγω της σχετικά υψηλής συχνότητας του φορτίου η θερμοκρασία παρουσιάζει έντονη μεταβολή στο χώρο και στο χρόνο. Ο προσδιορισμός της διάρκειας ζωής ενός δομικού στοιχείου κατά τη φάση του σχεδιασμού μπορεί να γίνει με τη βοήθεια πειραματικών διαδικασιών. Τα πειράματα όμως κόπωσης είναι δαπανηρά και χρονοβόρα και προφανώς απαιτούνται περισσότερες από μια πειραματικές δοκιμές. Οπότε, είναι εύλογο να υπάρχουν υπολογιστικά εργαλεία που να δίνουν τη δυνατότητα στο μηχανικό να εκτιμήσει την διάρκεια ζωής ή τη σοβαρότητα της βλάβης ενός εξαρτήματος. Τα περισσότερα υπολογιστικά μοντέλα αναφέρονται σε καθαρά μηχανικές καταπονήσεις. Έτσι υπάρχει πρόσφορο έδαφος για την ανάπτυξη υπολογιστικών εργαλείων για την ανάλυση προβλημάτων θερμικής και θερμομηχανικής κόπωσης. Τέτοιου είδους εργαλεία θα πρέπει να λαμβάνουν υπόψη το κλείσιμο των ρωγμών, που συμβαίνει λόγω των θερμικών παραμορφώσεων, διότι είναι δυνατόν να επηρεάζεται τοπικά το θερμοκρασιακό πεδίο. Επομένως, χρειάζεται επαναληπτική διαδικασία για τον προσδιορισμό του θερμικού και τασικού πεδίου που αλληλεπιδρούν. Είναι προφανές ότι η ανάλυση της θερμικής κόπωσης εξελίσσεται σε συνθέτη διαδικασία, που θα πρέπει να συμπεριλαμβάνει τον υπολογισμό της κατανομής της θερμοκρασίας, την τοπική επίδραση του άκρου της ρωγμής στο τασικό πεδίο καθώς και την επαφή των επιφανειών της ρωγμής. Η μέθοδος των συνοριακών στοιχείων είναι ικανή να αντιμετωπίζει τέτοιου είδους τοπικές επιδράσεις. Η παρούσα διατριβή επικεντρώνεται στην ανάπτυξη υπολογιστικού εργαλείου βασισμένου στα συνοριακά στοιχεία, για την πρόβλεψη της διάδοσης ρωγμών και την εκτίμηση της διάρκειας ζωής, εξαρτημάτων υπό θερμική και θερμομηχανική κόπωση. Έμφαση δίνεται σε περιπτώσεις που το θερμικό φορτίο προκαλεί κλείσιμο της ρωγμής και σε περιπτώσεις διεπιφανειακών ρωγμών, όπου το θερμοκρασιακό πεδίο επηρεάζεται από την θερμική αντίσταση ανάμεσα στις επιφάνειες της ρωγμής. Στο πρώτο κεφάλαιο γίνεται βιβλιογραφική ανασκόπηση σε εργασίες που εστιάζουν σε φαινόμενα κόπωσης και διάδοσης ρωγμών, καθώς και στην ανάπτυξη υπολογιστικών μοντέλων για την πρόβλεψη της διάδοσης ρωγμών. Επιπλέον, προσδιορίζεται λεπτομερώς το αντικείμενο της παρούσας διατριβής και εξηγείται η συνεισφορά της και τα καινοτόμα σημεία της. Στο δεύτερο κεφάλαιο περιγράφεται η ιδιόμορφη συμπεριφορά του άκρου της ρωγμής, δίνονται οι διατυπώσεις των μεγεθών θραύσης που χρησιμοποιούνται στην ανάλυση της κόπωσης και αναφέρονται τρόποι με τους οποίους μελετάται η διάδοση ρωγμών. Στο τρίτο κεφάλαιο περιγράφονται λεπτομερώς οι ολοκληρωτικές συνοριακές διατυπώσεις για την επίλυση προβλημάτων θερμοελαστικότητας. Στο τέταρτο κεφάλαιο περιγράφονται οι υπολογιστικές διαδικασίες που ακολουθούνται στην παρούσα εργασία για τον προσδιορισμό του πεδίου θερμοκρασιών και μετατοπίσεων, καθώς και ο τρόπος που προσομοιώνεται η διάδοση ρωγμής. Στο πέμπτο κεφάλαιο παρατίθενται τα αποτελέσματα που προέκυψαν από τις αναλύσεις για διάφορες περιπτώσεις, ενώ στο έκτο κεφάλαιο εξάγονται συμπεράσματα και διατυπώνονται προτάσεις για μελλοντική έρευνα. / The prediction of fatigue life is essential for the integrity and reliability of a structure when designing engineering components that undergo cyclic loading. In most cases, the mechanical cyclic loads are taken into account in order to evaluate the life and damage tolerance of structures with existing cracks. However, there exists a category of structures that experience severe thermal cycling that acts alongside the mechanical loads. Such structures include cooling system pipes, pressure vessels, pistons and combustion chambers of internal combustion engines, gas turbine blades and components of electrical circuits. Interfacial crack growth is of paramount importance when designing components that are protected by thermal barrier coatings in order to increase their endurance and efficiency. These types of structures are exposed to very intense thermo-mechanical cycling, which gradually causes delamination and eventually leads to spallation of the coating Numerical simulations, via the finite element method, are a common trend, when analysing the endurance of coated components. However, important aspects such as the heat exchange between the contacting faces and friction are not taken into account in fracture assessments of these components. The boundary element method is very attractive for crack-growth analyses because only the boundary is meshed, rather than the whole domain of the problem. In the present thesis, the boundary integral equations of uncoupled, time-dependent thermo-elasticity are employed to account for the time-varying nature of the thermal load. Our study discusses the influence of crack closure on quasi-static, sub-critical crack extension in the presence of thermo-mechanical cyclic loading. Appropriate thermal and mechanical boundary conditions are imposed on the numerical model to account for the contact state. The validity of the code to compute the temperature distribution under thermal cycling is checked through analytical solutions. Afterwards, a pure mode-I and mixed mode fracture problems in homogeneous material are analysed and the results are compared to other boundary element solutions. The singularity resulting from tractions and heat flux around the crack tip is effectively captured by singular quarter-point elements, while the fracture magnitudes can be computed using appropriate traction formulas. In these problems, the fatigue life is evaluated in terms of load cycle when the crack closure is considered. The number of cycles required for an existing crack to grow a certain length can be empirically predicted using the Paris’ law. The crack extension angle is evaluated by means of the maximum circumferential stress. The results are discussed, clearly indicating the impact of crack closure on fatigue life evaluation. The main conclusion is that crack closure should be incorporated into the analysis whenever the contact effect is inevitable. Otherwise, the fatigue life may be underestimated, leading to a conservative design. Finally, the sub-domain boundary element procedure is applied to interfacial cracks where the crack closure is more pronounced. Specifically, a case of a thermal barrier coating system is investigated. The thermal resistance between the contacting crack faces is incorporated into the procedure and it is assumed to be dependent on the contact pressure. If crack closure due to thermal distortion takes place, then the displacement and traction field may affect the heat flux between the crack faces, and the thermal and mechanical parts of the problem will need to be solved repeatedly until thermo-mechanical convergence is achieved. The results suggest that there are significant effects on the behaviour of stably growing cracks and the evaluation of failure capacity, emanating from crack closure, the amount of thermal resistance and the phase angle between the mechanical and thermal loads.

Orientador: Gabriela Rezende Fernandes / Resumo: Neste trabalho é apresentada uma formulação do MEC (Método dos Elementos de Contorno) considerando-se não-linearidade física para analisar microestruturas de materiais heterogêneos no contexto da análise em multi-escala. A microestrutura, também denominada como EVR (Elemento de Volume Representativo), é modelada como uma chapa em sub-regiões onde vazios ou inclusões podem ser considerados dentro da matriz, sendo diferentes propriedades elásticas e modelos constitutivos definidos para cada sub-região. A equação integral para o deslocamento é obtida a partir do Teorema de Betti, onde para considerar o fenômeno dissipativo, um campo de esforços iniciais é considerado. A equação algébrica da chapa é obtida após a discretização do contorno externo e interface em elementos e do domínio das subregiões em células. Na análise multi-escala cada ponto da estrutura (macrocontínuo) é representado por um EVR, onde o comportamento do material não é definido por um modelo constitutivo, mas através da solução do problema de equilíbrio do EVR quando sujeito à deformação referente ao ponto do macrocontínuo. O problema de equilíbrio do EVR é definido em termos da flutuação dos deslocamentos, sendo o mesmo satisfeito quando seu campo de forças se encontra em equilíbrio. Após a solução do EVR, os deslocamentos no contorno e as forças dissipativas são atualizados e as forças de superfície sobre o contorno recalculadas para se obter a tensão homogeneizada. O custo computacional obtido com a presente... (Resumo completo, clicar acesso eletrônico abaixo) / Mestre

Submitted by LUIS HENRIQUE DE REZENDE CROZARIOL null (luiscrozariol@gmail.com) on 2017-09-26T22:06:47Z No. of bitstreams: 1 Dissertação Luis 26-09.pdf: 3642072 bytes, checksum: 6bd5d174cfe596f9dc2def41f1270185 (MD5) / Approved for entry into archive by Monique Sasaki (sayumi_sasaki@hotmail.com) on 2017-09-28T13:33:18Z (GMT) No. of bitstreams: 1 crozariol_lhr_me_ilha.pdf: 3642072 bytes, checksum: 6bd5d174cfe596f9dc2def41f1270185 (MD5) / Made available in DSpace on 2017-09-28T13:33:18Z (GMT). No. of bitstreams: 1 crozariol_lhr_me_ilha.pdf: 3642072 bytes, checksum: 6bd5d174cfe596f9dc2def41f1270185 (MD5) Previous issue date: 2017-07-31 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Neste trabalho é apresentada uma formulação do MEC (Método dos Elementos de Contorno) considerando-se não-linearidade física para analisar microestruturas de materiais heterogêneos no contexto da análise em multi-escala. A microestrutura, também denominada como EVR (Elemento de Volume Representativo), é modelada como uma chapa em sub-regiões onde vazios ou inclusões podem ser considerados dentro da matriz, sendo diferentes propriedades elásticas e modelos constitutivos definidos para cada sub-região. A equação integral para o deslocamento é obtida a partir do Teorema de Betti, onde para considerar o fenômeno dissipativo, um campo de esforços iniciais é considerado. A equação algébrica da chapa é obtida após a discretização do contorno externo e interface em elementos e do domínio das subregiões em células. Na análise multi-escala cada ponto da estrutura (macrocontínuo) é representado por um EVR, onde o comportamento do material não é definido por um modelo constitutivo, mas através da solução do problema de equilíbrio do EVR quando sujeito à deformação referente ao ponto do macrocontínuo. O problema de equilíbrio do EVR é definido em termos da flutuação dos deslocamentos, sendo o mesmo satisfeito quando seu campo de forças se encontra em equilíbrio. Após a solução do EVR, os deslocamentos no contorno e as forças dissipativas são atualizados e as forças de superfície sobre o contorno recalculadas para se obter a tensão homogeneizada. O custo computacional obtido com a presente formulação é menor que aquele referente ao modelo desenvolvido pelo Método dos Elementos Finitos, sendo a resposta homogeneizada do EVR comparada ao modelo de elementos finitos a fim de validar a formulação apresentada nesse trabalho. / A BEM formulation, considering dissipative phenomena, to analyze microstructures of heterogeneous materials in the context of multi-scale analysis is presented. The microstructure, also denoted as RVE (Representative Volume Element), is modelled as a zoned plate where voids or inclusions can be considered inside a matrix, being different elastic properties and constitutive models defined for each sub-region. The integral representation for displacement is obtained from Betti’s Theorem, where to consider the dissipative phenomena, an initial forces field is considered. The plate algebraic equation is obtained after discretizing the external boundary and interfaces into elements and the sub-regions domain into cells. In the multi-scale analysis, each macro-continuum point is represented by a RVE, being the material behaviour not governed by a phenomenological constitutive model, but defined after the solution of the RVE equilibrium problem due to the macro strain. The RVE equilibrium problem is defined in terms of displacement fluctuations, being satisfied when the forces field is in equilibrium. After the RVE solution, the boundary displacements and dissipative forces are updated and the boundary tractions recalculated to obtain the homogenized stress. The computational cost obtained with the proposed formulation is smaller than the formulation developed by the Finite Element Method. Besides, the homogenized response is compared to the finite element model to show its accuracy.

Analisam-se problemas de interação solo-estrutura através de uma formulação tridimensional obtida da combinação entre o Método dos Elementos de Contorno (MEC) e o Método dos Elementos Finitos (MEF). Os elementos estruturais que interagem com o solo são modelados pelo Método dos Elementos Finitos. E o solo, considerado como um meio semi-infinito, homogêneo, elástico linear e isotrópico, é modelado pelo Método dos Elementos de Contorno, empregando a solução fundamental de Mindlin. A solução fundamental de Mindlin é particularmente adequada para o tipo de problemas em análise, ou seja, problemas envolvendo sólidos tridimensionais semi-infinitos, já que é necessário discretizar apenas a superfície carregada do solo e/ou a linha de carga, e não todo o sólido tridimensional. A discretização da estaca em vários elementos finitos de viga tridimensional, permitindo a consideração de estacas de qualquer tamanho e submetidas a qualquer tipo de carregamento é uma das principais contribuições do trabalho. Outra contribuição diz respeito à consideração da ação horizontal no sistema placa-estaca-solo, diretamente aplicada na placa, e não redistribuída no topo das estacas. Por isso, os elementos estruturais laminares (radiers, sapatas e blocos de fundação) são modelados por elementos finitos de casca plana, possibilitando a consideração dos efeitos de flexão e de membrana. Essa abordagem permite a análise tanto de um grupo de estacas com bloco de capeamento rígido quanto de um radier estaqueado de qualquer rigidez. O acoplamento entre o MEC e o MEF é feito através de uma formulação mista em que a matriz dos coeficientes de influência do solo obtida pelo MEC é adicionada à matriz de rigidez dos elementos estruturais obtida pelo MEF, resultando em uma matriz de rigidez equivalente. Exemplos numéricos de interação estaca-solo, placa-solo e casca plana-estaca-solo são resolvidos para verificar, validar e demonstrar a eficiência das formulações desenvolvidas e implementadas. / Soil-structure interaction problems are analyzed by a tridimensional formulation obtained combining the Boundary Element Method (BEM) and the Finite Element Method (FEM). Structural elements that interact with the soil are modeled by the Finite Element Method. And the soil, considered as a semi-infinite, homogeneous, linear elastic and isotropic medium, is modeled by the Boundary Element Method, using Mindlin\'s fundamental solution. Mindlin\'s fundamental solution is particularly suitable for the type of problems under analysis, viz. problems involving semi-infinite three-dimensional solids, since it is necessary to discretize only the loaded surface of the soil and/or the line-load, and not all the three-dimensional solid. The discretization of the pile in several three-dimensional beam finite elements, allowing the consideration of piles of any size and subjected to any type of loading, is one of the main contributions of this work. Another contribution is about the consideration of the horizontal load in the plate-pile-soil system, directly applied to the plate, and not redistributed on the top of the piles. Therefore, the laminar structural elements (rafts, footings, and foundation blocks) are modeled by flat shell finite elements, making possible the consideration of the effects of flexion and membrane. This approach allows the analysis of both a capped pile group and a piled raft of any stiffness. The coupling between the BEM and the FEM is performed through a mixed formulation in which the matrix of the soil\'s influence coefficients obtained by the BEM is added to the stiffness matrix of the structural elements obtained by the FEM, resulting in an equivalent stiffness matrix. Numerical examples of pile-soil, plate-soil and flat shell-pile-soil interaction are solved to verify, validate and demonstrate the efficiency of the developed and implemented formulations.

O objetivo do presente trabalho é o desenvolvimento de formulações e de códigos computacionais que possibilitem a análise bidimensional estática de meios contínuos anisotrópicos viscoplásticos reforçados ou não por fibras. Especificamente, as análises numéricas envolvem aplicações dos métodos dos elementos de contorno (MEC) e dos elementos finitos (MEF), comparando-se os resultados obtidos com respostas analíticas e experimentais, disponíveis na literatura, buscando-se assim, subsídios teóricos que permitam o entendimento de problemas mais gerais envolvendo meios anisotrópicos. Para tanto são empregados elementos finitos triangulares com aproximações cúbica e quadrática para os deslocamentos na modelagem dos domínios. Na consideração do reforço com fibras, elementos finitos de barras simples são empregados. A formulação desenvolvida proporciona também a consideração de distribuição randômica das fibras imersas no meio sem qualquer aumento dos graus de liberdade do problema analisado, diferindo-se assim, das formulações conhecidas até o momento. Com o MEC, a análise de plasticidade e viscoplasticidade em meios com anisotropia geral é feita de maneira original no trabalho, destacando-se a consideração de lei de fluxo plástico não-associativa e o tratamento de viscosidade apenas com integrais de contorno, sem a utilização de aproximações de domínio. Uma quantidade significativa de exemplos é apresentada, possibilitando a verificação da eficiência das formulações e dos códigos desenvolvidos / The objective of the present work is the development of formulations and computational codes that enable the static bidimensional analysis of the viscoplastic anisotropic medias reinforced, or not, by fibers. Specifically, the numerical analysis involve applications of the boundary elements method (BEM) and finite elements (FEM), comparing the results obtained with analytical and experimental solutions available in the literature, allowing the understanding of general problems in anisotropic media. Two-dimensional finite elements with cubic and quadrate approximations for the displacements are used to model domains. Reinforcements are modeled by truss finite elements. The developed formulation provides the consideration of random distribution of the fibers, without any additional degree of freedom of the problem. With the BEM, the plasticity and viscoplasticity analysis in general anisotropic medias is originally developed in the present work, emphasizing the consideration of non-associative plastic flow and the treatment of viscosity just with boundary integrals, without domain approximation. Various examples are shown in order to verify the efficiency of the proposed formulation and developed computational codes

Neste trabalho, propõe-se uma formulação não linear baseada no método dos elementos de contorno, para representação de domínios poro-elasto-plásticos reforçados. Esta formulação é apresentada para os casos saturado e não saturado. Para o problema poroso enrijecido um acoplamento com o método dos elementos finitos é empregado, e a técnica de mínimos quadrados permite a regularização dos deslocamentos e do vetor de forças de superfície ao longo das interfaces de acoplamento. São empregadas expressões analíticas para o tratamento das integrais de contorno e de domínio presentes na formulação do método dos elementos de contorno. A formulação de Biot é empregada para a descrição de meios porosos saturados e uma formulação energética baseada nos trabalhos de Coussy é adaptada para a extensão ao caso não saturado. Neste caso, a pressão capilar e energia das interfaces são levadas em consideração. O nível de saturação é descrito pelo modelo de Van Genuchten e o comportamento do esqueleto é descrito ou pelo modelo de Drucker-Prager ou pelo modelo de Cam-Clay modificado. O problema não linear obtido por uma descrição temporal associada a discretização espacial é resolvido pelo método de Newton-Raphson. No caso saturado, o operador tangente consistente é definido e utilizado para obtenção da solução do sistema. Exemplos numéricos são apresentados para validar a formulação proposta. / In this work a nonlinear formulation of the boundary element method (BEM) is proposed to deal with saturated and unsaturated poro-elasto-plastic 2D reinforced domains. To model reinforced porous domains a BEM/FEM (Finite Element Method) modified coupling technique is employed. The coupling is made by using the least square method to regularize the displacement and traction distributions along the interfaces. Analytical expressions have been derived for all boundary and domain integrals required for the formulation. The Biot formulation is used for the description of the saturated porous environments and an energetic consistent formulation based on work of Coussy is adopted for its extension to the framework of unsaturated porous media. In this case, the capillar pressure and the interface energy are taken into account. The Van Genuchten model is used for the determination of saturation level in non-saturated poro-elasto-plastic problems. The Drucker-Prager modified model if used for the saturated poro-elasto-plastic problems and the modified Cam-Clay model for the representation of non-saturated poro-elasto-plastic problems. For the saturated case, the consistent tangent operator is derived and employed inside a Newton procedure to solve non-linear problems. Numerical solutions are presented to validate the proposed models.

Submitted by Renata Lopes (renatasil82@gmail.com) on 2017-02-23T19:26:50Z No. of bitstreams: 1 francianeconceicaopeters.pdf: 4730497 bytes, checksum: 201c60342a8bf9edc9b308fa50fafa54 (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2017-02-24T12:08:43Z (GMT) No. of bitstreams: 1 francianeconceicaopeters.pdf: 4730497 bytes, checksum: 201c60342a8bf9edc9b308fa50fafa54 (MD5) / Made available in DSpace on 2017-02-24T12:08:43Z (GMT). No. of bitstreams: 1 francianeconceicaopeters.pdf: 4730497 bytes, checksum: 201c60342a8bf9edc9b308fa50fafa54 (MD5) Previous issue date: 2010-01-27 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A construção de imagens associadas à distribuição de condutividades no interior de um meio condutor a partir de injeção de corrente elétrica e medidas de potencial no contorno externo do corpo é uma técnica conhecida como tomografia por impedância elétrica. É um problema inverso que tem sido estudado visando aplicações biomédicas, monitoramento de processos industriais e investigação geofísica. Em alguns casos, é possível levar em consideração informações conhecidas sobre o domínio do corpo no processo de construção da imagem, recaindo no problema da detecção de inclusões que é o problema efetivamente tratado neste trabalho. Este problema pode ser resolvido por meio da minimização de uma função da diferença entre potenciais medidos no contorno e calculados para uma dada distribuição de condutividades. O presente trabalho desenvolve uma estratégia para a solução deste problema baseada na parametrização da geometria do contorno das inclusões cujas formas e dimensões se pretende determinar. O problema de minimização é resolvido por meio do Método de Levenberg-Marquardt e o problema direto via Método dos Elementos de Contorno. Para avaliar o desempenho da estratégia proposta são apresentados resultados numéricos envolvendo contornos definidos por splines, problemas com a presença de ruído nas medidas, avaliação de protocolos de injeção de corrente e medição de potencial elétrico e ainda uma aplicação voltada ao monitoramento cardíaco. / The images reconstruction of the conductivity distribution inside a conductive body based on electrical current injection and potential measurements on the outer boundary of this body is a technique known as electrical impedance tomography. This is an inverse problem that has been studied in biomedical applications, industrial process monitoring and geophysics investigation. In some cases, it is possible to take into account in the reconstruction process, informations about the body, leading to the problem of identifying inclusions, that is the problem actually treated in this work. This inverse problem can be solved by the minimization of a function, defined as the difference between the measured potentials and the computed ones for a given conductivity distribution. The present work describes a strategy to solve this problem based on the parametrization of the inclusions boundary, whose shape and size is intended to be determined. The minimization problem is solved via Levenberg-Marquardt Method and the forward one is solved via Boundary Elements Method. In order to evaluate the performance of the proposed strategy, numerical experiments with inclusions of boundaries defined by splines, problems with noisy data, current injection and potential measurement protocols and an application of the strategy to the cardiac function monitoring are presented.