İdari Hizmetler

O grupo de Modelagem Computacional do Laborat´orio de Mecˆanica dos Pavimen- tos da Universidade Federal do Cear´a (LMP/UFC) decidiu desenvolver um novo sistema

baseado no MEF para an´alises de pavimentos. Este sistema, denominado CAP3D, foi desen- volvido com os conceitos de orientac¸˜ao a objetos descritos neste cap´ıtulo e tem como base o c´odigo aberto FEMOOP (Martha e Parente-Jr, 2002). Os requisitos b´asicos deste novo sistema ´e a capacidade de lidar com problemas bi e tridimensionais de elementos finitos com materias de comportamento n˜ao-linear e dependentes do tempo.

A seguir encontra-se uma descric¸˜ao em forma resumida das classes e func¸˜oes im- plementadas no sistema CAP3D. Para um maior detalhamento o leitor ´e remetido `a literatura fonte desta revis˜ao constante nos trabalhos de Martha e Parente-Jr (2002) e Holanda et al. (2006b).

A organizac¸˜ao geral das classes do CAP3D pode ser vista na Figura 2.17, onde todas as relac¸ ˜oes s˜ao do tipo derivac¸˜ao (has a). As principais classes do programa s˜ao: Con-

trol, Node, Element, Shape, Analysis Model, Material, Integration Point, Constitutive Model

e Load. A classe Control representa o n´ıvel global discutido na sec¸˜ao anterior. Esta classe

Figura 2.17: Estrutura de classes do CAP3D (Holanda et al., 2006b).

fornece uma interface comum para os algoritmos de an´alise do problema. Atualmente ´e com- posta de 4 subclasses: LinearStatic, EquilibriumPath, QuasiStatic e LinearNewmark como ilustrado no esquema da Figura 2.18. A primeira ´e respons´avel por an´alises el´asticas lineares sem a considerac¸˜ao das forc¸as inerciais, a segunda ´e relativa a algoritmos iterativos usados no equil´ıbrio est´atico de estruturas de comportamento n˜ao-linear. As classes QuasiStatic e

LinearNewmark s˜ao respons´aveis pelas an´alises dependentes do tempo gerando estados de

equil´ıbrio baseados em intervalos de tempos anteriores. A subclasse LinearNewmark im- plementa o m´etodo de integrac¸˜ao direta de Newmark para an´alises dinˆamicas (considerac¸˜ao de forc¸as inerciais e amortecimento) onde os deslocamentos podem ser calculados de forma total ou incremental. Esta classe ´e proposta e foi implementada no ˆambito dos objetivos desta dissertac¸˜ao e as formulac¸˜oes necess´arias a sua implementac¸˜ao encontram-se descrititas nas sec¸˜oes seguintes. A subclasse QuasiStatic implementa an´alises lineares considerando forc¸as e deslocamentos dependentes do tempo, mas sem a considerac¸˜ao de forc¸as inerciais. Tipicamente, este tipo de an´alise tem sido usada na determinac¸˜ao de tens˜oes e deformac¸ ˜oes em pavimentos cuja camada de superf´ıcie tem comportamento viscoel´astico linear, mas as an´alises quasi-est´aticas n˜ao s˜ao restritas somente a este tipo de material.

Figura 2.18: Classe Control do CAP3D (Holanda et al., 2006b).

mento) de dados relativos aos n´os do modelo (coordenadas, condic¸˜oes de suporte, etc.). Ela tamb´em lida com algumas vari´aveis durante a execuc¸˜ao do programa, tais como os graus de liberdade (g.d.l.) nodais e os deslocamentos correntes.

Element ´e a subclasse que define o comportamento geral de um elemento finito.

As principais tarefas da classe s˜ao: (i) a indicac¸˜ao do n´umero e direc¸˜ao dos g.d.l. nodais

ativos; (ii) c´alculo dos vetores do elemento (g e f); (iii) c´alculo das matrizes (K, M e C) e

(iv) determinac¸˜ao das respostas do modelo (tens˜oes e deformac¸˜oes). Uma importante carac-

ter´ıstica desta classe ´e a capacidade de lidar com modelos multi-dimensionais em elementos finitos de uma maneira gen´erica.

A subclasse Shape em linhas gerais trabalha com a geometria e os aspectos de interpolac¸˜ao do elemento (dimens˜ao, topologia, n´umero de n´os, conectividade e ordem de interpolac¸˜ao). Atrav´es de suas classes derivadas para elementos uni, bi e tridimensionais fornece as func¸ ˜oes de forma e suas derivadas em relac¸˜ao a coordenadas param´etricas e carte- sianas e avalia a matriz Jacobiana do elemento.

Para uma modelagem eficiente dos materiais envolvidos nas an´alises, o CAP3D disp˜oe de 2 classes: Material e Constitutive Model. A classe Material ´e uma classe abstrata que lida com os dados dos diferentes modelos constitutivos implementados no programa, que atualmente incluem o el´astico linear, viscoel´astco linear e o modelo resiliente. O obje- tivo b´asico desta classe ´e armazenar as propriedades dos materias lidas a partir do arquivo de entrada de dados. A classe abstrata Constitutive Model tem como func¸˜ao b´asica o c´alculo do vetor de tens˜oes (σ) a partir de um dado vetor de deformac¸˜oes (ε) e a avaliac¸˜ao da matriz con- stitutiva (C), a ser usada na Equac¸˜ao 2.12 a partir do corrente estado de tens˜ao/deformac¸˜ao.

Deve ser notado que nos modelos (viscoel´asticos e elasto-pl´asticos) em que as tens˜oes n˜ao somente dependem do estado de tens˜ao corrente e sim de todo o hist´orico de carregamento, uma vez que cada ponto de integrac¸˜ao tem um hist´orico de tens˜oes, o sistema automaticamente cria um objeto do tipo Constitutive Model para cada ponto de integrac¸˜ao da malha de elementos finitos. Obviamente, o tipo de cada um destes objetos depende do material associado ao elemento correspondente. Devido a esta dependˆencia do hist´orico de

tens˜oes, algumas classes derivadas s˜ao necess´arias para o armazenamento das vari´aveis in- ternas (deformac¸ ˜oes em passos de tempo anteriores) e atualizac¸˜ao destas vari´aveis ap´os a convergˆencia do algoritmo global ser atingida.

Finalmente, a classe Load lida com as condic¸˜oes de contorno do problema e as forc¸as de corpo. Assim como classes anteriores, a classe abstrata Load permite o trata- mento generalizado para diferentes condic¸˜oes de carregamento dispon´ıveis. Al´em do mais, esta classe tamb´em avalia os vetores de carregamentos nodais e a aplicac¸˜ao destes em seus respectivos g.d.l. Assim como implementado na avaliac¸˜ao da Equac¸˜ao 2.12 na classe Con-

stitutive Model, a classe Load avalia a Equac¸˜ao 2.8 de uma maneira geral, independente da

dimens˜ao do elemento, forma, ordem de interpolac¸˜ao e equac¸˜ao diferencial do problema. Esta classe tamb´em inclui carregamentos vari´aveis no tempo, onde cada carregamento dis- tribu´ıdo e nodal tˆem uma func¸˜ao de tempo associada a ele. Desta forma, o termo q da Equac¸˜ao 2.8 pode ser definido por:

q= ¯qh(t) (2.140)

ondeq representa a variac¸˜ao espacial do carregamento distribu´ıdo dentro do elemento e¯ h(t)

´e uma func¸˜ao de tempo (constante, harmˆonica, ...) associada ao carregamento.

Para um maior detalhamento das definic¸˜oes das classes no sitema CAP3D o autor encoraja a leitura de Holanda et al. (2006b).