Predição de propagação de fissuras através de modelos constitutivos locais e técnica de construção progressiva da trajetória de descontinuidade

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Calibración de Modelos Constitutivos de Hormigón

En este artículo se describen diversos modelos constitutivos del comportamiento dilatante del hormigón. Se describe una metodología para caracterizar el comportamiento mecánico del hormigón. En particular se calibra el modelo de plasticidad y daño para hormigón de Abaqus/Standard y Abaqus/Explicit a partir de ensayos uniaxiales y tria xiales en probetas cilíndricas. El modelo se enriquece con una subrutina de usuario para capturar mejor la ductilidad del material sometido a moderadas presiones de confinamiento. Además se define una variable de estado para controlar la eliminación de elementos. Finalmente se muestran algunos ejemplos de aplicación de impactos de proyectiles deformables y rígidos sobre muros de hormigón armado, incluyendo técnicas especiales como la conversión de elementos sólidos a partículas (SPH).

Reforço de pilares por confinamento com mantas de fibras de carbono: avaliação de modelos de comportamento

A presente dissertação incide sobre o estudo dos efeitos do confinamento com materiais compósitos de polímeros reforçados com fibras de carbono (CFRP) em pilares de estruturas de betão armado. A motivação para este estudo surge da necessidade de aprofundar conhecimentos acerca do comportamento dos pilares de betão reforçados por confinamento com CFRP, uma vez que a sua aplicação apresenta uma crescente importância, por exemplo, para aumento da resistência e da ductilidade de estruturas de betão armado. Fez-se, inicialmente, uma breve revisão das técnicas de reforço convencionais utilizadas em pilares de betão armado, com ênfase no reforço exterior com polímeros reforçados com fibras. A elevada resistência à tração, à corrosão e à fadiga, o baixo peso volúmico, a versatilidade e a diversidade dos sistemas comercializados com CFRP tornam este material muito competitivo para este tipo de aplicação. Na sequência desse estudo, realizou-se uma revisão bibliográfica acerca dos modelos de comportamento que permitem prever o desempenho de pilares de betão confinados com CFRP, sujeitos a esforços de compressão. Como forma de análise desses modelos, desenvolveu-se uma ferramenta numérica em ambiente Mathworks - Matlab R2015a, que permitiu a obtenção e posterior comparação dos diagramas de tensão-extensão descritos pelos modelos desenvolvidos por Manfredi e Realfonzo (2001), Ferreira (2007) e Wei e Wu (2011). Por fim, comparam-se os resultados experimentais de Paula (2003) e de Rocca (2007) com os dos modelos constitutivos referidos anteriormente, analisando-se também a influência de vários fatores na eficácia do confinamento, tais como o boleamento, o número de camadas de CFRP e a geometria da secção transversal. Foram ainda comparados e discutidos resultados relativos ao confinamento parcial de pilares. Os resultados obtidos indicam que os modelos analíticos representam relativamente bem o andamento das curvas do betão confinado para secções circulares, quadradas e retangulares, verificando-se as principais discrepâncias nestas duas últimas tipologias de secção transversal, dada a dificuldade associada à quantificação de parâmetros associados ao seu comportamento (por exemplo, boleamento de arestas). No entanto, verificou-se igualmente que com um adequado boleamento de arestas (e consequente aumento da relação entre o raio de canto e a largura da secção de betão), bem como com um aumento do número de camadas de material compósito, é possível aumentar a tensão resistente e a extensão axial na rotura do betão à compressão.

Reforço de pilares por confinamento com mantas de fibras de carbono: avaliação de modelos de comportamento

A presente dissertação incide sobre o estudo dos efeitos do confinamento com materiais compósitos de polímeros reforçados com fibras de carbono (CFRP) em pilares de estruturas de betão armado. A motivação para este estudo surge da necessidade de aprofundar conhecimentos acerca do comportamento dos pilares de betão reforçados por confinamento com CFRP, uma vez que a sua aplicação apresenta uma crescente importância, por exemplo, para aumento da resistência e da ductilidade de estruturas de betão armado. Fez-se, inicialmente, uma breve revisão das técnicas de reforço convencionais utilizadas em pilares de betão armado, com ênfase no reforço exterior com polímeros reforçados com fibras. A elevada resistência à tração, à corrosão e à fadiga, o baixo peso volúmico, a versatilidade e a diversidade dos sistemas comercializados com CFRP tornam este material muito competitivo para este tipo de aplicação. Na sequência desse estudo, realizou-se uma revisão bibliográfica acerca dos modelos de comportamento que permitem prever o desempenho de pilares de betão confinados com CFRP, sujeitos a esforços de compressão. Como forma de análise desses modelos, desenvolveu-se uma ferramenta numérica em ambiente Mathworks - Matlab R2015a, que permitiu a obtenção e posterior comparação dos diagramas de tensão-extensão descritos pelos modelos desenvolvidos por Manfredi e Realfonzo (2001), Ferreira (2007) e Wei e Wu (2011). Por fim, comparam-se os resultados experimentais de Paula (2003) e de Rocca (2007) com os dos modelos constitutivos referidos anteriormente, analisando-se também a influência de vários fatores na eficácia do confinamento, tais como o boleamento, o número de camadas de CFRP e a geometria da secção transversal. Foram ainda comparados e discutidos resultados relativos ao confinamento parcial de pilares. Os resultados obtidos indicam que os modelos analíticos representam relativamente bem o andamento das curvas do betão confinado para secções circulares, quadradas e retangulares, verificando-se as principais discrepâncias nestas duas últimas tipologias de secção transversal, dada a dificuldade associada à quantificação de parâmetros associados ao seu comportamento (por exemplo, boleamento de arestas). No entanto, verificou-se igualmente que com um adequado boleamento de arestas (e consequente aumento da relação entre o raio de canto e a largura da secção de betão), bem como com um aumento do número de camadas de material compósito, é possível aumentar a tensão resistente e a extensão axial na rotura do betão à compressão.

Análise numérica de escorregamento em encostas

O perfeito entendimento de deslizamentos de terra tem sido objeto de grande interesse na engenharia geotécnica há muitos anos. Encostas naturais representam um desafio constante para os engenheiros, mais do que os taludes artificiais. Não é incomum que os escorregamentos causam perdas de vidas humanas e também danos materiais. As informações obtidas a partir de retro-análise de casos históricos, têm servido como uma ferramenta valiosa para a compreensão da estabilidade de taludes. Vale a pena notar que o primeiro método de cálculo de estabilidade de taludes, usando divisão de fatia, foi desenvolvido após um acidente na Noruega. Métodos de estabilidade por Equilíbrio Limite, desenvolvidos no século passado são atualmente utilizados na prática da engenharia. Apesar de suas deficiências, eles provaram ser confiáveis na determinação do Fator de Segurança. A geometria 2D e a abordagem completa da condição saturada são provavelmente as questões mais relevantes dos métodos atuais. A condição do solo não saturado está geralmente presente em taludes naturais e os métodos ignoram a influência da sucção do solo sobre a estabilidade do talude. O desenvolvimento de computadores no meio do século passado, tornaram os cálculos Equilíbrio Limite fáceis de serem executados. Por outro lado, eles não podem explicar a natureza do problema ou prever o comportamento do solo antes da ruptura. A década de 1960 marcou o início do desenvolvimento de programas de computador baseados em métodos numéricos que, entre outros, destaca-se, em geotecnia, a aplicação dos métodos de Elementos Finitos (FEM). Este método é ideal para muitos fins, uma vez que permite a simulação da sequência de construção e também a incorporação de diferentes modelos constitutivos. Este trabalho mostra a influência da sucção sobre a estabilidade de encostas naturais. Uma retroanálise de um escorregamento no Morro dos Cabritos, no Rio de Janeiro, foi feita através de um programa de análise de tensão (Plaxis 2D). Os resultados foram comparados com o FS obtido por programa de Equilíbrio Limite (Slide v5). O efeito de sucção do solo foi incorporada nas análises, através de sua influência sobre os parâmetros de resistência do solo e por diferentes possibilidades de infiltração de água na encosta. Este caso histórico foi estudado anteriormente por Gerscovich (1994), que na época, fez analise de fluxo e estabilidade utilizando um abordagem de Equilíbrio Limite 3D. Os resultados mostraram boa concordância entre a análise de Elementos Finitos e a simulação de Equilíbrio Limite, sem a presença de nível de água. Os resultados confirmaram que foi necessário haver pressões positivas para explicar a ruptura do talude.

Estudo paramétrico para um muro de gravidade submetido a carregamentos sísmicos.

A presente dissertação apresenta um estudo do comportamento de um muro de gravidade flexível submetido a carregamentos sísmicos. A influência do carregamento sísmico, e mais precisamente da variação da aceleração horizontal de pico é avaliada num estudo paramétrico, através da análise dos resultados obtidos para os deslocamentos e acelerações horizontais que ocorreram no paramento do muro de gravidade. Optou-se pela análise ao longo do tempo, introduzindo o carregamento sísmico ao modelo numérico do qual fazia parte o muro de gravidade através de 5 históricos temporais de acelerogramas horizontais normalizados em 0,05g, 0,10g, 0,15g, 0,20g e 0,25g. O evento sísmico de origem dos acelerogramas normalizados foi o terremoto ocorrido no Chile em 27 de fevereiro de 2010 e medido na estação em San Jose de Maipo, a 332,7km do epicentro do terremoto. Utiliza-se um software computacional aplicado à Geotecnia, o Plaxis, para a calibração de um modelo numérico em situação estática e posterior desenvolvimento das análises paramétricas em condições sísmicas. As análises realizadas para avaliação do comportamento do muro de gravidade em condição sísmica foram a verificação da influência da variação do acelerograma horizontal normalizada introduzido na base do modelo nos valores das acelerações e deslocamentos horizontais atuantes no muro e a verificação da influência do carregamento sísmico em comparação com a situação estática, comparando-se os valores dos deslocamentos horizontais obtidos na análise estática e nas análises dinâmicas. Os fatores que influenciam potencialmente nos resultados numéricos em condições sísmicas foram apresentados, destacando-se os aspectos relacionados à modelagem numérica em condição de sismo. Os modelos constitutivos oferecidos, a possibilidade de incorporação do carregamento sísmico na base do modelo e a possibilidade de consideração do amortecimento visco-elástico de Rayleigh nos materiais são as principais qualidades do Plaxis, utilizado na modelagem numérica. Avaliou-se positivamente a potencialidade do Plaxis, visto que o programa mostrou-se uma ferramenta capaz de simular o comportamento de muros de gravidade sujeitos a carregamentos sísmicos.

Simulação numérica de dutos enterrados, submetidos à perda de apoio e elevação localizada.

Os recentes desastres ocorridos no país, como o rompimento da adutora em Campo Grande e os desastres relacionados às enchentes urbanas, mostram a necessidade de desenvolvimento de pesquisas científicas que auxiliem na compreensão e no dimensionamento das estruturas projetadas para atender a demanda da população. Os métodos analíticos e experimentais mais utilizados possuem algumas limitações de ordem teórica ou prática. Por outro lado, os métodos numéricos, capazes de simular etapas construtivas e envolver materiais com diferentes modelos constitutivos numa mesma análise, buscam atender às necessidades práticas dos projetos de geotecnia e, ao mesmo tempo, complementam os modelos analíticos e experimentais. Nesse trabalho foram realizadas comparações entre resultados obtidos em ensaios experimentais e resultados extraídos do modelo computacional, buscando aumentar a compreensão sobre a interação solo-estrutura em relação à distribuição de tensões mobilizadas e aos deslocamentos e deformações provocados. A simulação numérica foi feita com a utilização do PLAXIS/3D, software de análise geotécnica baseado no método dos elementos finitos. Os ensaios foram confeccionados na Escola de Engenharia de São Carlos/USP por Costa (2005) e envolveram dutos enterrados submetidos à perda de apoio ou elevação localizada. O estudo experimental foi realizado através de modelos físicos compostos por um maciço de areia pura, contendo um tubo repousando sobre um alçapão no centro do vão. Os modelos físicos foram equipados com instrumental capaz de medir as deflexões e as deformações específicas ao longo do duto, além das tensões totais no maciço de solo circundante e na base do equipamento.

Análise tensão deformação de uma barragem de concreto em solo residual preponderantemente anisotrópico.

Estudos de fundações de obras de arte, como barragens de concreto-gravidade e barragens de terra, devem contemplar todos os esforços atuantes no maciço de fundação, principalmente as tensões e as deformações esperadas durante todo o processo construtivo e no período pós-construtivo. Quando essas estruturas são apoiadas sobre rochas de boa sanidade, a escolha do barramento geralmente favorece a implantação de barragens de concreto. Entretanto, quando os maciços de fundação são formados por solos, a opção técnica geralmente mais bem aceita é quanto à utilização de barragens de terra. Em ambos os casos, as análises de estabilidade e de deformação são desenvolvidas por métodos analíticos bem consolidados na prática. Nas condições mais adversas de fundação, seja em rochas ou em solos com marcante anisotropia e estruturas reliquiares herdadas da rocha-mãe, a utilização de modelos constitutivos anisotrópicos em análises por elementos finitos propicia simulações mais realistas dessas feições estruturais, contribuindo para o seu melhor conhecimento. O presente trabalho teve por objetivo demonstrar a utilização de um modelo constitutivo anisotrópico no estudo da fundação da Barragem San Juan, localizada na República Dominicana, que foi concebida como uma estrutura tipo concreto-gravidade apoiada sobre solos residuais jovens altamente anisotrópicos. Nessa obra, apesar de sua pequena altura, a presença marcante de descontinuidades ensejou um estudo mais detalhado do comportamento tensão-deformação da fundação, levando em conta o levantamento detalhado da atitude das descontinuidades presentes no maciço e os resultados de ensaios de resistência em planos paralelos e normais às descontinuidades. Para a estimativa de deformações, os respectivos módulos de Young (Es) foram estimados com base em correlações com a resistência à penetração SPT desenvolvidas nesta dissertação, a partir de um estudo estatístico baseado em várias publicações disponíveis na literatura. As análises numéricas por elementos finitos foram desenvolvidas através do programa Plaxis 2D, utilizando-se como modelo constitutivo aquele denominado The Jointed Rock Model, que é particularmente recomendado para análises de estabilidade e deformação de materiais anisotrópicos. Os resultados das análises numéricas foram comparados com as análises de equilíbrio limite elaboradas para o projeto executivo da referida barragem, pelo programa Slope, utilizando o método rigoroso de Morgenstern e Price, que se mostrou conservador. Os resultados das análises numéricas mostraram sua inequívoca versatilidade para a escolha de opções de reforço da fundação, através de dentes que objetivavam o aumento das condições de estabilidade da barragem.

Development of numerical methodologies for parameter identification and shape optimization in metal forming simulations

Por parte da indústria de estampagem tem-se verificado um interesse crescente em simulações numéricas de processos de conformação de chapa, incluindo também métodos de engenharia inversa. Este facto ocorre principalmente porque as técnicas de tentativa-erro, muito usadas no passado, não são mais competitivas a nível económico. O uso de códigos de simulação é, atualmente, uma prática corrente em ambiente industrial, pois os resultados tipicamente obtidos através de códigos com base no Método dos Elementos Finitos (MEF) são bem aceites pelas comunidades industriais e científicas Na tentativa de obter campos de tensão e de deformação precisos, uma análise eficiente com o MEF necessita de dados de entrada corretos, como geometrias, malhas, leis de comportamento não-lineares, carregamentos, leis de atrito, etc.. Com o objetivo de ultrapassar estas dificuldades podem ser considerados os problemas inversos. No trabalho apresentado, os seguintes problemas inversos, em Mecânica computacional, são apresentados e analisados: (i) problemas de identificação de parâmetros, que se referem à determinação de parâmetros de entrada que serão posteriormente usados em modelos constitutivos nas simulações numéricas e (ii) problemas de definição geométrica inicial de chapas e ferramentas, nos quais o objetivo é determinar a forma inicial de uma chapa ou de uma ferramenta tendo em vista a obtenção de uma determinada geometria após um processo de conformação. São introduzidas e implementadas novas estratégias de otimização, as quais conduzem a parâmetros de modelos constitutivos mais precisos. O objetivo destas estratégias é tirar vantagem das potencialidades de cada algoritmo e melhorar a eficiência geral dos métodos clássicos de otimização, os quais são baseados em processos de apenas um estágio. Algoritmos determinísticos, algoritmos inspirados em processos evolucionários ou mesmo a combinação destes dois são usados nas estratégias propostas. Estratégias de cascata, paralelas e híbridas são apresentadas em detalhe, sendo que as estratégias híbridas consistem na combinação de estratégias em cascata e paralelas. São apresentados e analisados dois métodos distintos para a avaliação da função objetivo em processos de identificação de parâmetros. Os métodos considerados são uma análise com um ponto único ou uma análise com elementos finitos. A avaliação com base num único ponto caracteriza uma quantidade infinitesimal de material sujeito a uma determinada história de deformação. Por outro lado, na análise através de elementos finitos, o modelo constitutivo é implementado e considerado para cada ponto de integração. Problemas inversos são apresentados e descritos, como por exemplo, a definição geométrica de chapas e ferramentas. Considerando o caso da otimização da forma inicial de uma chapa metálica a definição da forma inicial de uma chapa para a conformação de um elemento de cárter é considerado como problema em estudo. Ainda neste âmbito, um estudo sobre a influência da definição geométrica inicial da chapa no processo de otimização é efetuado. Este estudo é realizado considerando a formulação de NURBS na definição da face superior da chapa metálica, face cuja geometria será alterada durante o processo de conformação plástica. No caso dos processos de otimização de ferramentas, um processo de forjamento a dois estágios é apresentado. Com o objetivo de obter um cilindro perfeito após o forjamento, dois métodos distintos são considerados. No primeiro, a forma inicial do cilindro é otimizada e no outro a forma da ferramenta do primeiro estágio de conformação é otimizada. Para parametrizar a superfície livre do cilindro são utilizados diferentes métodos. Para a definição da ferramenta são também utilizados diferentes parametrizações. As estratégias de otimização propostas neste trabalho resolvem eficientemente problemas de otimização para a indústria de conformação metálica.

Caracterização experimental e modelação numérica de espumas metálicas

As espumas de alumínio são materiais ultraleves, o que as torna atractivas para um largo espectro de aplicações comerciais na área da defesa, na indústria automóvel e aeroespacial, entre outras. Actualmente, há um grande interesse na utilização de espumas de alumínio em componentes estruturais ultraleves, incorporados em sistemas de absorção de energia para protecção contra o impacto. O recurso à simulação numérica para resolver problemas de engenharia em várias áreas é cada vez mais comum. A modelação numérica dos materiais assume vital importância quando o problema envolve a análise de processos tecnológicos como, por exemplo, a conformação plástica de materiais, ou a análise de estruturas. Deste modo, torna-se imprescindível garantir que a modelação dos materiais é de tal forma rigorosa que permite simular o melhor possível o seu comportamento real nas condições concretas da análise a realizar. A forma mais comum de garantir o rigor dos modelos utilizados é a validação dos modelos numéricos tendo por base resultados experimentais. Neste trabalho, fez-se a caracterização do comportamento mecânico das espumas de alumínio com nome comercial ALPORAS!, obtidas pelo processo de fabrico denominado expansão directa do metal fundido por adição de um agente expansor. Esta caracterização consistiu num conjunto de ensaios experimentais quer no regime quasi-estático, quer no regime dinâmico. No regime quasi-estático realizaram-se ensaios de compressão uniaxial e de compressão multiaxial. Para a caracterização no regime dinâmico foram realizados ensaios em barras de Hopkinson de polimetil-metacrilato (PMMA). Com base nos resultados experimentais obtidos determinaram-se os parâmetros dos dois modelos constitutivos para espumas metálicas implementados no programa comercial Abaqus™/Explicit. Estes modelos, e os respectivos parâmetros determinados, foram validados reproduzindo numericamente alguns ensaios experimentais quasi-estáticos e dinâmicos. Assim, verificou-se a adequabilidade dos modelos em diversas condições quer em termos de esforços quer em termos de regime de taxa de deformação Por último, desenvolveu-se uma estrutura inovadora para absorção de energia durante um impacto, constituída por componentes perfilados em liga de alumínio e por componentes em espumas de alumínio. Esta estrutura foi testada exclusivamente com recurso à simulação numérica, utilizando os modelos constitutivos validados anteriormente.

Modelización constitutiva y computacional del daño y la fractura de materiales compuestos

En el trabajo se definen modelos constitutivos que permiten reproducir el proceso de fallo de estructuras de materiales compuestos en distintas escalas bajo cargas estáticas. Se define un modelo constitutivo para determinar la respuesta de estructuras de materiales compuestos mediante la teoría de laminados. El modelo es validado mediante un programa de ensayos experimentales con probetas con un agujero central geométricamente similares. Se muestra la capacidad del modelo de detectar el efecto tamaño. Se define un modelo constitutivo para materiales transversalmente isótropos bajo estados tridimensionales de tensión. El modelo se valida analizando numéricamente el proceso de agrietamiento de la matriz. Finalmente se desarrolla un modelo analítico para determinar el agrietamiento de la matriz y la delaminación entre las capas.

Simulation of interlaminar and intralaminar damage in polymer-based composites for aeronautical applications under impact loading

La aplicación de materiales compuestos de matriz polimérica reforzados mediante fibras largas (FRP, Fiber Reinforced Plastic), está en gradual crecimiento debido a las buenas propiedades específicas y a la flexibilidad en el diseño. Uno de los mayores consumidores es la industria aeroespacial, dado que la aplicación de estos materiales tiene claros beneficios económicos y medioambientales. Cuando los materiales compuestos se aplican en componentes estructurales, se inicia un programa de diseño donde se combinan ensayos reales y técnicas de análisis. El desarrollo de herramientas de análisis fiables que permiten comprender el comportamiento mecánico de la estructura, así como reemplazar muchos, pero no todos, los ensayos reales, es de claro interés. Susceptibilidad al daño debido a cargas de impacto fuera del plano es uno de los aspectos de más importancia que se tienen en cuenta durante el proceso de diseño de estructuras de material compuesto. La falta de conocimiento de los efectos del impacto en estas estructuras es un factor que limita el uso de estos materiales. Por lo tanto, el desarrollo de modelos de ensayo virtual mecánico para analizar la resistencia a impacto de una estructura es de gran interés, pero aún más, la predicción de la resistencia residual después del impacto. En este sentido, el presente trabajo abarca un amplio rango de análisis de eventos de impacto a baja velocidad en placas laminadas de material compuesto, monolíticas, planas, rectangulares, y con secuencias de apilamiento convencionales. Teniendo en cuenta que el principal objetivo del presente trabajo es la predicción de la resistencia residual a compresión, diferentes tareas se llevan a cabo para favorecer el adecuado análisis del problema. Los temas que se desarrollan son: la descripción analítica del impacto, el diseño y la realización de un plan de ensayos experimentales, la formulación e implementación de modelos constitutivos para la descripción del comportamiento del material, y el desarrollo de ensayos virtuales basados en modelos de elementos finitos en los que se usan los modelos constitutivos implementados.

Um modelo elasto-viscoplástico para análise de peças de concreto estrutural, submetidas a estados planos de tensão, através do método dos elementos finitos

O objetivo deste trabalho é desenvolver um modelo computacional, baseado no método dos elementos finitos, para o estudo de peças de concreto armado e protendido submetidas a estados planos de tensão. O estudo abrange situações de carga de curta e longa duração, onde consideram-se fluência e retração do concreto e relaxação do aço. São utilizados modelos constitutivos elasto-viscoplásticos para descrever o comportamento dos materiais. Implementou-se um modelo de camadas superpostas para melhor representar o comportamento do concreto, onde o material é composto de diversas camadas que sofrem a mesma deformação. Cada camada possui diferentes características materiais e a tensão total é obtida pela soma das diferentes contribuições de cada camada. Para a fissuração da concreto, utilizou-se um modelo de fissuras distribuídas, que leva em conta a contribuição do concreto entre fissuras. Tanto a amadura passiva como a de pratensão são introduzidas no modelo como uma linha de material mais rígido dentro do elemento de concreto. Os deslocamentos ao longo da armadura são referenciados aos deslocamentos nodais do elemento de concreto. Deste modo, obtém-se uma matriz de rigidez para a armadura com as mesmas dimensões que a matriz de rigidez do elemento de concreto, A matriz de rigidez do elemento concreto-aço é a soma das duas matrizes. Considera-se aderência perfeita entre o concreto e o aço. Os resultados obtidos com esse programa computacionai são comparados com valores experimentais disponíveis.

Modelo da ligação entre concreto e armadura na análise de estruturas de concreto pelo método dos elementos finitos

Este trabalho apresenta um modelo computacional para análise tridimensional de peças de concreto armado, utilizando o método dos elementos finitos. São utilizados modelos constitutivos elasto-viscoplásticos para representar o comportamento dos materiais, quando submetidos a cargas de curta ou longa duração. Se a estrutura está submetida a cargas de curta duração, o modelo desenvolvido fornece, quando alcançado o estado estável, a solução do problema elastoplástico. Em contrapartida, se a estrutura está submetida a cargas de longa duração, o modelo pode representar o comportamento viscoelástico com envelhecimento do concreto. Foi utilizado um modelo de camadas superpostas para assegurar uma representação adequada do comportamento real do concreto ao longo do tempo. Deste modo, admite-se que o concreto é constituído por um número conveniente de camadas, que apresentam a mesma deformação. Para a fissuração do concreto, utilizou-se um modelo de fissuras distribuídas, que leva em consideração a contribuição do concreto entre fissuras O objetivo deste estudo é a modelagem da armadura e do reforço nas peças de concreto armado, através de um modelo de aderência entre o concreto e o material adjacente. A armadura é introduzida no modelo como uma linha de material mais rígido dentro do elemento de concreto. E, a consideração da degradação da aderência é realizada através de um modelo onde os efeitos da aderência imperfeita são incluídos pela introdução de graus de liberdade para os deslocamentos relativos entre os materiais. Os resultados obtidos com este programa computacional são comparados com valores experimentais disponíveis na literatura, visando validar o modelo matemático e a metodologia numérica.

Modelagem de situações de punção em lajes de concreto armado através do método dos elementos finitos

O objetivo deste trabalho é o desenvolvimento de um programa computacional, elaborado através do método dos elementos finitos, para análise tridimensional de situações de punção em lajes de concreto armado com e sem armadura de cisalhamento. São utilizados modelos constitutivos elasto-viscoplásticos para representar o comportamento do concreto. Para a simulação do concreto submetido a um carregamento de longa duração, empregou-se o modelo de camadas superpostas. Este modelo admite que o concreto é constituído por um número conveniente de camadas com a mesma deformação, porém características materiais diferentes. A tensão total é obtida pela soma das contribuições de cada camada. A armadura é introduzida no modelo como uma linha de material mais rígido no interior do elemento de concreto. Considera-se aderência perfeita entre o concreto e o aço. Assim, os deslocamentos ao longo das barras de aço são determinados a partir dos deslocamentos nodais dos elementos de concreto. Desta forma, a matriz de rigidez da armadura tem a mesma dimensão da matriz de rigidez do elemento de concreto, as quais são somadas para a obtenção da matriz de rigidez total Para representar o comportamento após à fissuração do concreto, é utilizado o modelo de fissuras distribuídas, onde apenas as propriedades dos materiais são modificadas e a malha de elementos finitos não necessita de atualização. Sugerem-se duas formulações para a consideração da contribuição do concreto entre fissuras e da redução na capacidade de transferência de corte devido ao aumento da abertura da fissura. Este último é incluído através de uma simplificação, que consiste em atribuir ao módulo de elasticidade transversal, correspondente ao plano fissurado, um valor reduzido. Para comprovar a validade do programa computacional desenvolvido, comparam-se os resultados numéricos com valores experimentais de ensaios realizados por Coelho (1999), Trautwein (2001) e Fusco (1988). Obtiveram-se boas aproximações para lajes com e sem armadura de cisalhamento. Desta forma, o programa computacional desenvolvido pode ser empregado para generalizar resultados experimentais em lajes cogumelo com variação na resistência do concreto, na bitola das barras de armadura, no espaçamento entre as barras e no número de camadas de armadura de cisalhamento.