Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O presente trabalho apresenta um estudo referente à aplicação da abordagem Bayesiana como técnica de solução do problema inverso de identificação de danos estruturais, onde a integridade da estrutura é continuamente descrita por um parâmetro estrutural denominado parâmetro de coesão. A estrutura escolhida para análise é uma viga simplesmente apoiada do tipo Euler-Bernoulli. A identificação de danos é baseada em alterações na resposta impulsiva da estrutura, provocadas pela presença dos mesmos. O problema direto é resolvido através do Método de Elementos Finitos (MEF), que, por sua vez, é parametrizado pelo parâmetro de coesão da estrutura. O problema de identificação de danos é formulado como um problema inverso, cuja solução, do ponto de vista Bayesiano, é uma distribuição de probabilidade a posteriori para cada parâmetro de coesão da estrutura, obtida utilizando-se a metodologia de amostragem de Monte Carlo com Cadeia de Markov. As incertezas inerentes aos dados medidos serão contempladas na função de verossimilhança. Três estratégias de solução são apresentadas. Na Estratégia 1, os parâmetros de coesão da estrutura são amostrados de funções densidade de probabilidade a posteriori que possuem o mesmo desvio padrão. Na Estratégia 2, após uma análise prévia do processo de identificação de danos, determina-se regiões da viga potencialmente danificadas e os parâmetros de coesão associados à essas regiões são amostrados a partir de funções de densidade de probabilidade a posteriori que possuem desvios diferenciados. Na Estratégia 3, após uma análise prévia do processo de identificação de danos, apenas os parâmetros associados às regiões identificadas como potencialmente danificadas são atualizados. Um conjunto de resultados numéricos é apresentado levando-se em consideração diferentes níveis de ruído para as três estratégias de solução apresentadas. / This work presents a study on the application of Bayesian approach as a technique for solving the inverse problem of structural damage identification, where the integrity of the structure is continuously described by a structural cohesion parameter. The structure chosen for analysis is a simply supported Euler - Bernoulli beam. The damage identification is based on changes in the impulse response of the structure caused by the presence thereof. The direct problem is solved by the finite element method (FEM), which, in turn, is parameterized by the cohesion parameter of the structure. The problem of identifying damages is formulated as an inverse problem, whose solution, from the Bayesian framework, is a posteriori probability distribution of the cohesion parameter, obtained using the sampling methodology of Monte Carlo with Markov Chain. The uncertainties inherent to the measured data will be included in the likelihood function. Three solution strategies are presented. In the Strategy 1, the cohesion parameters of the structure are sampled from probability density functions a posteriori that have the same standard deviation. In the Strategy 2, after a previous analysis of the damage identification process, are determined potentially damaged regions and the cohesion parameters associated with these regions are sampled from probability density functions a posteriori that have different deviations. In the Strategy 3, after a preliminary analysis of the damage identification process, only the parameters associated with regions identifed as potentially damaged are updated. A set of numerical results are presented taking into account different noise levels for the three considered strategies.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / O presente trabalho apresenta um estudo referente à aplicação da abordagem Bayesiana como técnica de solução do problema inverso de identificação de danos estruturais, onde a integridade da estrutura é continuamente descrita por um parâmetro estrutural denominado parâmetro de coesão. A estrutura escolhida para análise é uma viga simplesmente apoiada do tipo Euler-Bernoulli. A identificação de danos é baseada em alterações na resposta impulsiva da estrutura, provocadas pela presença dos mesmos. O problema direto é resolvido através do Método de Elementos Finitos (MEF), que, por sua vez, é parametrizado pelo parâmetro de coesão da estrutura. O problema de identificação de danos é formulado como um problema inverso, cuja solução, do ponto de vista Bayesiano, é uma distribuição de probabilidade a posteriori para cada parâmetro de coesão da estrutura, obtida utilizando-se a metodologia de amostragem de Monte Carlo com Cadeia de Markov. As incertezas inerentes aos dados medidos serão contempladas na função de verossimilhança. Três estratégias de solução são apresentadas. Na Estratégia 1, os parâmetros de coesão da estrutura são amostrados de funções densidade de probabilidade a posteriori que possuem o mesmo desvio padrão. Na Estratégia 2, após uma análise prévia do processo de identificação de danos, determina-se regiões da viga potencialmente danificadas e os parâmetros de coesão associados à essas regiões são amostrados a partir de funções de densidade de probabilidade a posteriori que possuem desvios diferenciados. Na Estratégia 3, após uma análise prévia do processo de identificação de danos, apenas os parâmetros associados às regiões identificadas como potencialmente danificadas são atualizados. Um conjunto de resultados numéricos é apresentado levando-se em consideração diferentes níveis de ruído para as três estratégias de solução apresentadas. / This work presents a study on the application of Bayesian approach as a technique for solving the inverse problem of structural damage identification, where the integrity of the structure is continuously described by a structural cohesion parameter. The structure chosen for analysis is a simply supported Euler - Bernoulli beam. The damage identification is based on changes in the impulse response of the structure caused by the presence thereof. The direct problem is solved by the finite element method (FEM), which, in turn, is parameterized by the cohesion parameter of the structure. The problem of identifying damages is formulated as an inverse problem, whose solution, from the Bayesian framework, is a posteriori probability distribution of the cohesion parameter, obtained using the sampling methodology of Monte Carlo with Markov Chain. The uncertainties inherent to the measured data will be included in the likelihood function. Three solution strategies are presented. In the Strategy 1, the cohesion parameters of the structure are sampled from probability density functions a posteriori that have the same standard deviation. In the Strategy 2, after a previous analysis of the damage identification process, are determined potentially damaged regions and the cohesion parameters associated with these regions are sampled from probability density functions a posteriori that have different deviations. In the Strategy 3, after a preliminary analysis of the damage identification process, only the parameters associated with regions identifed as potentially damaged are updated. A set of numerical results are presented taking into account different noise levels for the three considered strategies.

Este trabalho propõe uma técnica de modelagem multiescala concorrente do concreto considerando duas escalas distintas: a mesoescala, onde o concreto é modelado como um material heterogêneo, e a macroescala, na qual o concreto é tratado como um material homogêneo. A heterogeneidade da estrutura mesoscópica do concreto é idealizada considerando três fases distintas, compostas pelos agregados graúdos e argamassa (matriz), estes considerados materiais homogêneos, e zona de transição interfacial (ZTI), tratada como a parte mais fraca entre as três fases. O agregado graúdo é gerado a partir de uma curva granulométrica e posicionado na matriz de forma aleatória. Seu comportamento mecânico é descrito por um modelo constitutivo elástico-linear, devido a sua maior resistência quando comparado com as outras duas fases do concreto. Elementos finitos contínuos com alta relação de aspecto em conjunto com um modelo constitutivo de dano são usados para representar o comportamento não linear do concreto, decorrente da iniciação de fissuras na ZTI e posterior propagação para a matriz, dando lugar à formação de macrofissuras. Os elementos finitos de interface com alta relação de aspecto são inseridos entre todos os elementos regulares da matriz e entre os da matriz e agregados, representando a ZTI, tornando-se potenciais caminhos de propagação de fissuras. No estado limite, quando a espessura do elemento de interface tende a zero (h ?0) e, consequentemente, a relação de aspecto tende a infinito, estes elementos apresentam a mesma cinemática da aproximação contínua de descontinuidades fortes (ACDF), sendo apropriados para representar a formação de descontinuidades associados a fissuras, similar aos modelos coesivos. Um modelo de dano à tração é proposto para representar o comportamento mecânico não linear das interfaces, associado à formação de fissuras, ou até mesmo ao eventual fechamento destas. A fim de contornar os problemas causados pela malha de elementos finitos de transição entre as malhas da macro e da mesoescala, que, em geral, apresentam diferenças expressivas 5 de refinamento, utiliza-se uma técnica recente de acoplamento de malhas não conformes. Esta técnica é baseada na definição de elementos finitos de acoplamento (EFAs), os quais são capazes de estabelecer a continuidade de deslocamento entre malhas geradas de forma completamente independentes, sem aumentar a quantidade total de graus de liberdade do problema, podendo ser utilizados tanto para acoplar malhas não sobrepostas quanto sobrepostas. Para tornar possível a análise em multiescala em casos nos quais a região de localização de deformações não pode ser definida a priori, propõe-se uma técnica multiescala adaptativa. Nesta abordagem, usa-se a distribuição de tensões da escala macroscópica como um indicador para alterar a modelagem das regiões críticas, substituindo-se a macroescala pela mesoescala durante a análise. Consequentemente, a malha macroscópica é automaticamente substituída por uma malha mesoscópica, onde o comportamento não linear está na iminência de ocorrer. Testes numéricos são desenvolvidos para mostrar a capacidade do modelo proposto de representar o processo de iniciação e propagação de fissuras na região tracionada do concreto. Os resultados numéricos são comparados com os resultados experimentais ou com aqueles obtidos através da simulação direta em mesoescala (SDM). / A concurrent multiscale analysis of concrete is presented, in which two distinct scales are considered: the mesoscale, where the concrete is modeled as a heterogeneous material and the macroscale that treats the concrete as a homogeneous material. The mesostructure heterogeneities are idealized as three phase materials composed of the coarse aggregates, mortar matrix and the interfacial transition zone (ITZ). The coarse aggregates are generated from a grading curve and placed into the mortar matrix randomly. Their behavior is described using an elastic-linear constitutive model due to their significant higher strength when compared with the other two phases of the concrete. Special continuum finite elements with a high aspect ratio and a damage constitutive model are used to describe the nonlinear behavior associated to the propagation of cracks, which initiates in the ITZ and then propagates to the mortar matrix given place to a macro-crack formation. These interface elements with a high aspect ratio are inserted in between all regular finite elements of the mortar matrix and in between the mortar matrix and aggregate elements, representing the ITZ. In the limit case, when the thickness of interface elements tends to zero (h ?0) and consequently the aspect ratio tends to infinite, these elements present the same kinematics as the continuous strong discontinuity approach (CSDA), so that they are suitable to represent the formation of discontinuities associated to cracks, similar to cohesive models. A tensile damage model is proposed to model the nonlinear mechanical behavior of the interfaces, associated to the crack formation and also to the possible crack closure. To avoid transition meshes between the macro and the mesoscale meshes, a new technique for coupling non-matching meshes is used. This technique is based on the definition of coupling finite elements (CFEs), which can ensure the continuity of displacement between independent meshes, without increasing the total number of degrees of freedom of the problem. This technique can be used to couple non-overlapping and overlapping meshes.To make possible the concurrent multiscale analysis, where the strain localization region cannot be defined a priori, an adaptive multiscale model is proposed. In this approach the macroscale stress distribution is used as an indicator to properly change from the macroscale to the mesoscale modeling in the critical regions during the analysis. Consequently, the macroscopic mesh is automatically replaced by a mesoscopic mesh where the nonlinear behavior is imminent. A variety of tests are performed to show the ability of the proposed methodology in predicting the behavior of initiation and propagation of cracks in the tensile region of the concrete. The numerical results are compared with the experimental ones or with those obtained by the direct simulation in mesoscale (DSM).

O presente estudo trata do problema da corrosão de armaduras em estruturas de concreto armado submetidas à penetração de cloretos. É amplamente reconhecido nos meios técnico e científico que, dentre os diversos processos causadores de patologias nas estruturas de concreto armado, destacam-se aqueles que desencadeiam a corrosão das armaduras. O processo de corrosão é iniciado quando a concentração de cloretos na interface armadura/concreto atinge um valor limite, despassivando a camada química protetora ao redor da armadura. A utilização de uma abordagem probabilística é capaz de tratar o problema de forma mais consistente, uma vez que as variáveis envolvidas no fenômeno possuem um alto grau de aleatoriedade associado. A formulação utilizada para a análise não linear mecânica emprega o método dos elementos finitos, utilizando o modelo de dano de Mazars para descrever o comportamento do concreto e o modelo elastoplástico para descrever o comportamento do aço. A lei de Fick é empregada para descrever o mecanismo de difusão dos íons cloretos no interior dos elementos de concreto armado e leis empíricas são utilizadas para representar a penalização das áreas de aço das armaduras bem como a redução de sua tensão de escoamento. A formulação utilizada para a análise probabilística emprega o FORM (First Order Reliability Method), o método de superfície de resposta e a simulação de Monte Carlo. Enfoque especial é dado na análise de estruturas hiperestáticas, com o objetivo de determinar a probabilidade de falha do sistema, sendo esta uma das contribuições desse estudo. Os resultados obtidos indicam de que forma os efeitos corrosivos considerados contribuem para a redução da durabilidade estrutural. / This work deals with the problem of reinforcement corrosion of concrete structures subjected to chloride penetration. It is widely recognized in the technical and scientific communities that, among the different processes which cause pathologies in reinforced concrete structures, reinforcement corrosion is one of the most relevant. The corrosion process starts when chloride concentration at the reinforcement/concrete interface reaches a threshold value, causing depassivation of the protecting chemical layer surrounding the armor. A probabilistic approach is capable of dealing with the problem in a more consistent manner, since the variables involved in the phenomenon have a high degree of randomness. The formulation used in the mechanical non-linear analysis utilizes the finite element method, employing Mazars\' damage model to describe the concrete behavior and the elastoplastic model to describe the steel behavior. Fick\'s law is used to describe the diffusion mechanism of chloride ions inside the concrete elements and empirical laws are used to represent the steel area reduction as well as the reduction of steel yield stress. The formulation used for the probabilistic analysis employs the FORM (First Order Reliability Method), the response surface method and the Monte Carlo simulation. A contribution of this study concerns the analysis of hyper-static structures, in order to evaluate the failure probability of the system, giving special attention to this type of problem. The obtained results show how the considered corrosive effects influences the structural durability reduction.

L'objectif de ce travail est la mise en oeuvre de la mécanique de la rupture et d'un modèle d'endommagement afin de vérifier l'intégrité du réservoir cryotechnique d'ARIANE 5 en présence de défauts. Ces derniers sont généralement assimilés à des fissures surfaciques semi-elliptiques. L'approche globale basée sur l'analyse FAD fournit une prédiction trop conservative lorsque cette méthode utilise le facteur d'intensité K et la ténacité KIc. Ce conservatisme peut toutefois être partiellement levé à condition d'introduire la notion d'intégrale J et d'utiliser une ténacité dérivée de la courbe JR - ∆a afin d'autoriser une légère avancée de la fissure dans le cas d'un matériau ductile. Cette courbe JR - ∆a s'obtient à partir d'essais sur des éprouvettes fissurées de différentes géométries (CT, SENT et SCT). L'approche locale de la rupture constitue cependant la seule alternative possible afin de statuer sur l'acceptabilité d'un défaut suffisamment important pour générer une plasticité étendue. Cette approche nécessite alors de décrire précisément les champs de contrainte et déformation en pointe de fissure. Par conséquent il est nécessaire de tenir compte du caractère anisotrope de la plasticité, lié au procédé de laminage, lors de la modélisation du comportement élasto-plastique del'alliage d'aluminium 2219 T87. Dans cette étude un critère phénoménologique est utilisé pour modéliser la plasticité dont l'anisotropie est décrite par différents coefficients. L'optimisation de ces paramètres repose sur la comparaison entre les résultats issus de simulations numériques et ceux obtenus par des essais sur des éprouvettes lisses et entaillées. L'approche local de la rupture repose également sur la compréhension des micromécanismes de rupture. L'alliage d'aluminium 2219 T87 étudié dans le cadre de cette thèse présenteune rupture ductile dont la phase de germination s'effectue sur les inclusions Al2Cu et constitue le mécanisme prépondérant de la rupture. L'endommagement est modélisé par le modèle de Gurson modifié par Needleman et Tvergaard (GTN). L'anisotopie de rupture est simulée par deux lois de germination différentes selon le sens de laminage et le sens transverse. / The purpose of this work is to use fracture mechanics and damage model to assess the structural integrity of the ARIANE 5 tank containing defects. These flaws are usually assumed to be semi-elliptical surface cracks.This situation can be analysed using the FAD (Failure Assessment Diagram) approach. This method is based on the stress intensity factor K and a material fracture toughness (Kmat) which can be the plane strain fracture toughness (KIC) for a conservative prediction. A more realistic prediction is achieved by employing the integral J and a material parameter derived from crack resistance curve JR - ∆a . This curve is obtained by using cracked test samples with different geometries (CT, SENT and SCT). The local approach to fractureis an alternative method to assess failure in particular in the case of large scale plasticity. In this case, the previous method is not able to give accurate predictions. Analyses based on local approach concepts requirean accurate evaluation of strain and stress fields near the crack tip. It is therefore important to take plastic anisotropy into account to model the plastic behavior of the rolled sheets. In this work this is done using amacroscopic phenomenological model. The parameters of this model are determined by minimizing of the difference between the simulation results and the test measurements for smooth and notched tensile test samples. The local approach to fracture is also based on the understanding of the failure micromechanismsof materials. Ductile fracture of metallic materials can be described as a three stage process. The first is void initiation at inclusions, the second is void growth and the third is void coalescence. Nucleation of defectsat Al2Cu inclusions is the main damage mechanism. In this work, the proposed model for ductile fracture is based on the Gurson model modified by Needleman and Tvergaard (GTN). This model is extended to account for fracture anisotropic by considering two nucleation gaussian laws respectively for rolling andtransverse directions.

Les céramiques sont couramment utilisées depuis les années 1960 comme matériaux constitutifs de blindages. En effet, grâce à leurs très bonnes propriétés physiques et mécaniques, elles permettent, pour un même niveau de protection, un gain de masse important par rapport aux blindages métalliques. Cependant, la microstructure d’une céramique peut avoir une forte influence sur sa résistance à l’impact. Le but de cette thèse est, à partir d’essais de caractérisation et en se basant sur l’utilisation de quatre nuances de carbure de silicium présentant des microstructures différentes, d’éclaircir les liens entre microstructure et performance à l’impact. Les campagnes expérimentales de compression dynamique et d’écaillage sont basées sur une utilisation innovante du moyen GEPI installé au CEA Gramat. Pour l’étude du comportement en compression dynamique des céramiques, il a permis d’utiliser la technique d’analyse lagrangienne et ainsi de remonter à l’évolution de la résistance des matériaux au cours du chargement. Pour les essais d’écaillage, il a permis, entre autres, une étude de la sensibilité à la vitesse de déformation de la résistance en traction dynamique. La caractérisation de la fragmentation dynamique est quant à elle basée sur des essais d’impact sur la tranche. Un essai innovant d’impact sur céramique préalablement fragmentée a également été dimensionné et réalisé. Ces différents essais expérimentaux ont permis de mettre en évidence et de comprendre l’influence de la microstructure du matériau sur son comportement face aux différents types de sollicitations. L’ensemble des résultats expérimentaux a été comparé à des simulations numériques permettant de valider les lois de comportement utilisées. Le modèle de fragmentation des matériaux fragiles DFH (Denoual-Forquin-Hild) a ainsi montré de très bonnes capacités à simuler le comportement des céramiques sous chargement de traction dynamique (écaillage et fragmentation) / Since the sixties, ceramics are commonly used as armour materials. Indeed, thanks to their interesting physical and mechanical properties, they allow a significant weight benefit in comparison to monolithic steel plate armours. However, the microstructure of the ceramic may have a strong influence on its penetration resistance. Based on characterisation tests and on the use of four silicon carbide grades, this work aims to highlight the links between the microstructure and the ballistic efficiency. Experimental compressive and spalling tests are based on the use of the GEPI device. For studying the compressive dynamic behaviour, it allows using the lagrangian analysis method and characterising the yield strength of the material. For studying the tensile dynamic behaviour, it allows assessing the strain-rate sensitivity of the spall strength. An analysis of the fragmentation process is performed based on Edge-On Impact tests. Moreover, an innovating impact test on fragmented ceramics has been designed and performed. The different experimental results allow a better understanding of the influence of the ceramic microstructure on its behaviour under the different loadings. All the experimental data have been compared to numerical results allowing validating the constitutive models. The DFH (Denoual-Forquin-Hild) damage model of brittle materials showed very good capacities to simulate the tensile dynamic behaviour of ceramics (spalling and fragmentation)

Passive safety is a well-known term. This term can be further categorized into different topics of the car passive safety, restraint systems, safety assistants (ABS, ESP, ASR, etc.). One of these topics is passive safety of the fuel system. Safety and tightness of the fuel system must be guaranteed even under non-standard conditions, for example a collision against a fixed obstacle. This issue is not often mentioned in the field of car safety. It is considered a standard. Passive safety of the fuel system is often ensured using various interesting technical solutions and devices, usually patented ones. The development of these solutions is supported by numerical simulations in different stages of development process. The doctoral thesis deals with impact loading of the plastic components of the fuel system, in particular Fuel Supply Module (FSM), which is mounted inside the fuel tank. The flange is the most important part of the fuel supply module from the car safety point of view. The flange closes FSM on the external side of the fuel tank. The thesis focuses on the finite element analysis of the complete or partial FSM, and the flange itself during impact loading. The main objective of this thesis are numerical material models, taking into account important aspects of the mechanical behavior of polymer materials during impact loading. There are a lot of ad hoc invented or standardized experiments described in this thesis. These experiments are used for estimation of the material parameters or comparison of numerical analysis vs real conditions, or tests. The solver LS-DYNA was mainly used for numerical simulations. The final results of this thesis brings new quantified knowledge about behavior of the Typical Semi-Crystal Polymer (TSCP), not only for impact loading. The practical part of this thesis defines new methodology for the numerical simulation approach of impact loading for FSM. This methodology is directly usable for new product development. A lot of numerical material models were developed and tested. The best results were achieved using numerical material model *MAT_24 with combination of *MAT_ADD_EROSION card. The limits and parameters for this numerical material model was estimated empirically during conducting experiments. The numerical material model SAMP-1 was partly solved in this doctoral thesis, but more detail study will be given in future works.

Als Großdemonstrator des Sonderforschungsbereichs 532 soll auf dem Gelände der Fakultät für Bauingenieurwesen der RWTH Aachen ein Ausstellungspavillon mit einer Dachkonstruktion aus Textilbeton errichtet werden. Der Beitrag zeigt die automatisierte Bemessung dieser komplexen Tragstruktur mit Hilfe eines hierfür entwickelten numerischen Bemessungsstools. Die Auswertung erfolgt für alle Lastfallkombinationen nach DIN 1055-100 und unter Berücksichtigung der Reduktion der Tragfähigkeit der textilen Bewehrung infolge einer Umlenkung im Riss sowie einer Erhöhung der aufnehmbaren Textilspannungen infolge Biegebeanspruchung. Zur Ermittlung der Tragfähigkeit der Textilbetonschale wurden am Institut für Massivbau Dehnkörper- und Biegeversuche durchgeführt, die in Bewehrungsgrad und Dicke dem realen Bauwerk am Schalenrand entsprachen. / Within the collaborate research center 532 at RWTH Aachen University the construction of an exposition hall with a double-curved roof structure consisting of textile reinforced concrete is currently planned. A numerical tool has been developed to calculate the necessary number of reinforcement layers. Further, the tool evaluates the characteristic stresses of the load case combinations in the ultimate limit state taking into account a reduction of the load-bearing capacity of the textile reinforcement due to alignment of the rovings in the crack bridge and a better activation of the inner filaments for bending induced tension. The resistance of the material has been determined in experimental investigations of tensile and bending specimens of the same thickness and reinforcement ratio as the planned shell structure.

Structure-property quantification of an A36 steel alloy was the focus of this study in order to calibrate and validate a plasticity-damage model. The microstructural parameters included grain size, particle size, particle number density, particle nearest neighbor distances, and percent of ferrite and pearlite. The mechanical property data focused on stress-strain behavior under different applied strain rates (0.001/s, 0.1/s, and 1000/s), different temperatures (293 K and 573 K), and different stress states (compression, tension, and torsion). Notch tension tests were also conducted to validate the plasticity-damage model. Also, failure of an A36 I-beam was examined in cyclic loads, and the crack growth rates were quantified in terms of fatigue striation data. Dynamic strain aging was observed in the stress-strain behavior giving rise to an important point that there exists a critical temperature for such behavior.

[pt] A representação matemática de comportamentos complexos em materiais exige formulações constitutivas sofisticada, como é o caso de modelos com múltiplas superfícies de plastificação. Assim, um modelo elastoplástico complexo demanda um procedimento robusto de integração das equações de evolução plástica. O desenvolvimento de esquemas de integração para modelos de plasticidade é um tópico de pesquisa importante, já que estes estão diretamente ligados à acurácia e eficiência de simulações numéricas de materiais como metais, concretos, solos e rochas. O desempenho da solução de elementos finitos é diretamente afetado pelas características de convergência do procedimento de atualização de estados. Dessa forma, este trabalho explora a implementação de modelos constitutivos complexos, focando em modelos genéricos com múltiplas superfícies de plastificação. Este estudo formula e avalia algoritmos de atualização de estado que formam uma estrutura robusta para a simulação de materiais regidos por múltiplas superfícies de plastificação. Algoritmos de integração implícita são desenvolvidos com ênfase na obtenção de robustez, abrangência e flexibilidade para lidar eficazmente com aplicações complexas de plasticidade. Os algoritmos de atualização de estado, baseados no método de Euler implícito e nos métodos de Newton-Raphson e Newton-Krylov, são formulados utilizando estratégias de busca unidimensional para melhorar suas características de convergência. Além disso, é implementado um esquema de subincrementação para proporcionar mais robustez ao procedimento de atualização de estado. A flexibilidade dos algoritmos é explorada, considerando várias condições de tensão, como os estados plano de tensões e plano de deformações, num esquema de integração único e versátil. Neste cenário, a robustez e o desempenho dos algoritmos são avaliados através de aplicações clássicas de elementos finitos. Além disso, o cenário desenvolvido no contexto de modelos com múltiplas superfícies de plastificação é aplicado para formular um modelo elastoplástico com dano acoplado, que é avaliado através de ensaios experimentais em estruturas de concreto. Os resultados obtidos evidenciam a eficácia dos algoritmos de atualização de estado propostos na integração de equações de modelos com múltiplas superfícies de plastificação e a sua capacidade para lidar com problemas desafiadores de elementos finitos. / [en] The mathematical representation of complex material behavior requires a sophisticated constitutive formulation, as it is the case of multisurface plasticity. Hence, a complex elastoplastic model demands a robust integration procedure for the plastic evolution equations. Developing integration schemes for plasticity models is an important research topic because these schemes are directly related to the accuracy and efficiency of numerical simulations for materials such as metals, concrete, soils and rocks. The performance of the finite element solution is directly influenced by the convergence characteristics of the state-update procedure. Therefore, this work explores the implementation of complex constitutive models, focusing on generic multisurface plasticity models. This study formulates and evaluates state-update algorithms that form a robust framework for simulating materials governed by multisurface plasticity. Implicit integration algorithms are developed with an emphasis on achieving robustness, comprehensiveness and flexibility to handle cumbersome plasticity applications effectively. The state-update algorithms, based on the backward Euler method and the Newton-Raphson and Newton-Krylov methods, are formulated using line search strategies to improve their convergence characteristics. Additionally, a substepping scheme is implemented to provide further robustness to the state-update procedure. The flexibility of the algorithms is explored, considering various stress conditions such as plane stress and plane strain states, within a single, versatile integration scheme. In this scenario, the robustness and performance of the algorithms are assessed through classical finite element applications. Furthermore, the developed multisurface plasticity background is applied to formulate a coupled elastoplastic-damage model, which is evaluated using experimental tests in concrete structures. The achieved results highlight the effectiveness of the proposed state-update algorithms in integrating multisurface plasticity equations and their ability to handle challenging finite element problems.