[pt] A resposta dinâmica de sistemas estruturais não lineares tem sido um item de grande interesse nas pesquisas em engenharia civil. Problemas onde há interação base flexível-estrutura são de grande importância na análise estrutural, já que a maioria das estruturas civis é apoiada sobre sistemas flexíveis (solo ou sistemas de apoio com dissipação de energia). Nesta área, o estudo de sistemas submetidos a ações sísmicas é um tópico relevante, já que estas solicitações têm um grande conteúdo de frequências, o que pode influenciar consideravelmente as respostas da estrutura. Neste contexto, o conhecimento da resposta em frequência de estruturas não lineares sob uma excitação de base é uma ferramenta útil para avaliar os potenciais efeitos de ações sísmicas sobre estes sistemas. Na presente tese é desenvolvida uma metodologia de análise não linear dinâmica de sistemas estruturais reticulados sob excitações de base, considerando não linearidade geométrica e apoios flexíveis, representados por molas unidimensionais, com comportamento elasto-plástico. Através de uma análise paramétrica é avaliada a variabilidade das respostas de sistemas esbeltos submetidos a ações sísmicas reais, sismos artificiais, assim como ações sísmicas sucessivas. O problema no espaço é resolvido pelo método dos elementos finitos. Para a análise em frequência, é apresentada uma metodologia baseada no método do balanço harmônico e no método de Galerkin, juntamente com técnicas de continuação para a obtenção das curvas de ressonância não lineares. O problema no tempo é abordado através da integração das equações de movimento pelos métodos de Runge-Kutta e Newmark, associado ao método de Newton-Raphson. / [en] The dynamic response of nonlinear structures has been a topic of interest in civil engineering research. Problems in which base-structure interaction is present have a great importance in structural analysis, since most structures rests on flexibel systems (soil or supports with dissipation). In this research area, the study of structures under the action of seismic loads represent a relevant topic, since this kind of excitations may excite several vibration modes and thus influence strongly the dynamic response. In this context, the prediction of the nonlinear structural behavior in frequency domain of structures under base excitation is a useful resource to assess the potential effects of sismic loads on these systems. In this thesis, a methodology for nonlinear dynamic analysis of plane frame structures under base excitation is presented considering geometric nonlinearity and elastic supports represented by elasto-plastic unidimensional springs. Trough a parametric analysis, the variability of the dynamic responses of slender structural systems under the actions of real earthquakes, synthetics earthquakes, as well as the action of multiple earthquakes is assessed. The structural systems here analyzed are discretized in space using a nonlinear finite element formulation. For the response in frequency domain, a scheme based on the Balance Harmonic Method and the Galerkin method, in conjunction with continuation methods, is formulated to obtain the nonlinear resonance curves. The nonlinear dynamic response in the time domain is calculated by direct integration of the equations of motion. For this, the Runge-Kutta method and the Newmark method in association with the iterative Newton-Raphson scheme are employed.

Este trabalho apresenta uma estratégia numérica desenvolvida usando o método dos elementos finitos para simular o processo de falha de compósitos cimentícios reforçados com fibras de aço. O material é descrito como um compósito composto por três fases: matriz cimentícia (pasta, argamassa ou concreto), fibras descontínuas discretas, e interface fibra-matriz. Um novo esquema de acoplamento para malhas de elementos finitos não-conformes foi desenvolvido para acoplar as malhas geradas independentes, da matriz cimentícia e de uma nuvem de fibras de aço, baseado na utilização de novos elementos finitos desenvolvidos, denominados elementos finitos de acoplamento. Utilizando este esquema de acoplamento, um procedimento não-rígido é proposto para a modelagem do complexo comportamento não linear da interface fibra-matriz, utilizando um modelo constitutivo de dano apropriado para descrever a relação entre a tensão de cisalhamento (tensão de aderência) e deslizamento relativo entre a matriz e cada fibra de aço individualmente. Este esquema também foi adotado para considerar a presença de barras de aço para as análises de estruturas de concreto armado. As fibras de aço são modeladas usando elementos finitos lineares com dois nós (elementos de treliça) com modelo material elastoplástico. As fibras são posicionadas usando uma distribuição randômica uniforme isotrópica, considerando o efeito parede. Uma abordagem contínua e outra descontínua são investigadas para a modelagem do comportamento frágil da matriz cimentícia. Para a primeira, é utilizado um modelo de dano isotrópico com duas variáveis de dano para descrever o comportamento de dano à tração e à compressão. A segunda emprega uma técnica de fragmentação de malha que utiliza elementos finitos degenerados, posicionados entre todos os elementos finitos que formam a matriz cimentícia. Para esta técnica é proposto um modelo constitutivo à tração, compatível com a abordagem descontínua forte contínua, para prever a propagação de fissura. Para acelerar o cálculo e aumentar a robustez dos modelos de dano contínuos para simular o processamento de falhas, um esquema de integração implícito-explícito é utilizado. Exemplos numéricos são apresentados ao longo do desenvolvimento desta tese. Inicialmente, exemplos numéricos com um único reforço são apresentados para validar a técnica desenvolvida e para investigar à influência das propriedades geométricas 7 das fibras e sua posição em relação à superfície de falha. Posteriormente, exemplos mais complexos são considerados envolvendo uma nuvem de fibras. Nestes casos, atenção especial é dada à influência da distribuição das fibras no comportamento do compósito relacionado ao processo de fissuração. Comparações com resultados experimentais demonstram que a aplicação da ferramenta numérica para modelar o comportamento de compósitos cimentícios reforçados com fibras de aço é muito promissora e pode ser utilizada como uma importante ferramenta para melhor entender os efeitos dos diferentes aspectos envolvidos no processo de falha deste material. / This work presents a numerical strategy developed using the Finite Element Method (FEM) to simulate the failure process of Steel Fiber Reinforced Cementitious Composites (SFRCCs). The material is described as a composite made up by three phases: a cementitious matrix (paste, mortar or concrete), discrete discontinuous fibers, and a fiber-matrix interface. A novel coupling scheme for non-matching finite element meshes has been developed to couple the independent generated meshes of the bulk cementitious matrix and a cloud of discrete discontinuous fibers based on the use of special finite elements developed, termed Coupling Finite Elements (CFEs). Using this approach, a nonrigid coupling procedure is proposed for modeling the complex nonlinear behavior of the fiber-matrix interface by adopting an appropriate constitutive damage model to describe the relation between the shear stress (adherence stress) and the relative sliding between the matrix and each fiber individually. This scheme has also been adopted to account for the presence of regular reinforcing bars in the analysis of reinforced concrete structural elements. The steel fibers are modeled using two-node finite elements (truss elements) with a one-dimensional elastoplastic constitutive model. They are positioned using an isotropic uniform random distribution, considering the wall effect of the mold. Continuous and discontinuous approaches are developed to model the brittle behavior of the bulk cementitious matrix. For the former, an isotropic damage model including two independent scalar damage variables for describing the composite behavior under tension and compression is considered. The discontinuous approach is based on a mesh fragmentation technique that employs degenerated solid finite elements in between all regular (bulk) elements. In this case, a tensile damage constitutive model, compatible with the Continuum Strong Discontinuity Approach (CSDA), is proposed to predict crack propagation. To increase the computability and robustness of the continuum damage models used to simulate the failure processes in both of the strategies, an implicit-explicit integration scheme is used. Numerical analyses are performed throughout the presentation of the work. Initially, numerical examples with a single reinforcement are presented to validate the technique and to investigate the influence of the fibers geometrical properties and its position relative to the crack surface. Then, more complex examples involving a cloud of steel fibers are considered. In these cases, special attention is given to the analysis of the influence of the fiber distribution on the composite behavior relative to the cracking process. Comparisons with experimental results demonstrate that the application of the numerical tool for modeling the behavior of SFRCCs is very promising and may constitute an important tool for better understanding the effects of the different aspects involved in the failure process of this material.

La méthode d'intégration produit a été proposée pour résoudre des équations linéaires de Fredholm de deuxième espèce singulières dont la solution exacte est régulière, au moins continue. Dans ce travail on adapte cette méthode à des équations dont la solution est juste intégrable. On étudie également son extension au cas non linéaire posé dans l'espace des fonctions intégrables. Ensuite, on propose une autre manière de mettre en oeuvre la méthode d'intégration produit : on commence par linéariser l'équation par une méthode de type Newton puis on discrétise les itérations de Newton par la méthode d'intégration produit / The product integration method has been proposed for solving singular linear Fredholm equations of the second kind whose exact solution is smooth, at least continuous. In this work, we adapt this method to the case where the solution is only integrable. We also study the nonlinear case in the space of integrable functions. Then, we propose a new version of the method in the nonlinear framework : we first linearize the eqaution by a Newton type method and then discretize the Newton iterations by the product integration method

Este trabalho apresenta uma estratégia numérica desenvolvida usando o método dos elementos finitos para simular o processo de falha de compósitos cimentícios reforçados com fibras de aço. O material é descrito como um compósito composto por três fases: matriz cimentícia (pasta, argamassa ou concreto), fibras descontínuas discretas, e interface fibra-matriz. Um novo esquema de acoplamento para malhas de elementos finitos não-conformes foi desenvolvido para acoplar as malhas geradas independentes, da matriz cimentícia e de uma nuvem de fibras de aço, baseado na utilização de novos elementos finitos desenvolvidos, denominados elementos finitos de acoplamento. Utilizando este esquema de acoplamento, um procedimento não-rígido é proposto para a modelagem do complexo comportamento não linear da interface fibra-matriz, utilizando um modelo constitutivo de dano apropriado para descrever a relação entre a tensão de cisalhamento (tensão de aderência) e deslizamento relativo entre a matriz e cada fibra de aço individualmente. Este esquema também foi adotado para considerar a presença de barras de aço para as análises de estruturas de concreto armado. As fibras de aço são modeladas usando elementos finitos lineares com dois nós (elementos de treliça) com modelo material elastoplástico. As fibras são posicionadas usando uma distribuição randômica uniforme isotrópica, considerando o efeito parede. Uma abordagem contínua e outra descontínua são investigadas para a modelagem do comportamento frágil da matriz cimentícia. Para a primeira, é utilizado um modelo de dano isotrópico com duas variáveis de dano para descrever o comportamento de dano à tração e à compressão. A segunda emprega uma técnica de fragmentação de malha que utiliza elementos finitos degenerados, posicionados entre todos os elementos finitos que formam a matriz cimentícia. Para esta técnica é proposto um modelo constitutivo à tração, compatível com a abordagem descontínua forte contínua, para prever a propagação de fissura. Para acelerar o cálculo e aumentar a robustez dos modelos de dano contínuos para simular o processamento de falhas, um esquema de integração implícito-explícito é utilizado. Exemplos numéricos são apresentados ao longo do desenvolvimento desta tese. Inicialmente, exemplos numéricos com um único reforço são apresentados para validar a técnica desenvolvida e para investigar à influência das propriedades geométricas 7 das fibras e sua posição em relação à superfície de falha. Posteriormente, exemplos mais complexos são considerados envolvendo uma nuvem de fibras. Nestes casos, atenção especial é dada à influência da distribuição das fibras no comportamento do compósito relacionado ao processo de fissuração. Comparações com resultados experimentais demonstram que a aplicação da ferramenta numérica para modelar o comportamento de compósitos cimentícios reforçados com fibras de aço é muito promissora e pode ser utilizada como uma importante ferramenta para melhor entender os efeitos dos diferentes aspectos envolvidos no processo de falha deste material. / This work presents a numerical strategy developed using the Finite Element Method (FEM) to simulate the failure process of Steel Fiber Reinforced Cementitious Composites (SFRCCs). The material is described as a composite made up by three phases: a cementitious matrix (paste, mortar or concrete), discrete discontinuous fibers, and a fiber-matrix interface. A novel coupling scheme for non-matching finite element meshes has been developed to couple the independent generated meshes of the bulk cementitious matrix and a cloud of discrete discontinuous fibers based on the use of special finite elements developed, termed Coupling Finite Elements (CFEs). Using this approach, a nonrigid coupling procedure is proposed for modeling the complex nonlinear behavior of the fiber-matrix interface by adopting an appropriate constitutive damage model to describe the relation between the shear stress (adherence stress) and the relative sliding between the matrix and each fiber individually. This scheme has also been adopted to account for the presence of regular reinforcing bars in the analysis of reinforced concrete structural elements. The steel fibers are modeled using two-node finite elements (truss elements) with a one-dimensional elastoplastic constitutive model. They are positioned using an isotropic uniform random distribution, considering the wall effect of the mold. Continuous and discontinuous approaches are developed to model the brittle behavior of the bulk cementitious matrix. For the former, an isotropic damage model including two independent scalar damage variables for describing the composite behavior under tension and compression is considered. The discontinuous approach is based on a mesh fragmentation technique that employs degenerated solid finite elements in between all regular (bulk) elements. In this case, a tensile damage constitutive model, compatible with the Continuum Strong Discontinuity Approach (CSDA), is proposed to predict crack propagation. To increase the computability and robustness of the continuum damage models used to simulate the failure processes in both of the strategies, an implicit-explicit integration scheme is used. Numerical analyses are performed throughout the presentation of the work. Initially, numerical examples with a single reinforcement are presented to validate the technique and to investigate the influence of the fibers geometrical properties and its position relative to the crack surface. Then, more complex examples involving a cloud of steel fibers are considered. In these cases, special attention is given to the analysis of the influence of the fiber distribution on the composite behavior relative to the cracking process. Comparisons with experimental results demonstrate that the application of the numerical tool for modeling the behavior of SFRCCs is very promising and may constitute an important tool for better understanding the effects of the different aspects involved in the failure process of this material.