pontifícia universidade católica de minas gerais análise de novos ...

230
(158)
onde x1,...,x8 e y1,...,y8 são as coordenadas dos nós do elemento e xg, yg são as
coordenadas globais do ponto, n1,...,n8 são as funções de forma, α1,... α8 são os
valores das isoclínicas nos nós, αg é o valor de α para o presente ponto plotado,
δ1, δ8 são os valores das isocromáticas nos nós e δg é o valor das isocromáticas
correspondente ao ponto atual que está sendo plotado. A qualidade do contorno
depende do intervalo de mudança de fase. Quanto mais fino o intervalo, melhor é o
contorno resultante. Para cada um dos elementos plotados, a extensão das faces na
direção x (xext) e direção y (yext) é calculado a partir das coordenadas globais dos
nós (Figura 31(a)). Com estes, a mudança de fase em intervalos Δr e Δs (Figura 31
(b)) são calculados de tal forma que eles estão igual ao tamanho do pixel ou mais
finas do que o tamanho do pixel.
Figura 31: (a) Representação de um elemento típico para calcular incrementos
na direção x e y e (b) representação da digitalização incrementos no sistema
de coordenadas naturais.
Coordenadas globais Coordenadas locais
Fonte: Adaptado por Ashokan e Ramesh (2009).
Ashokan e Ramesh (2009) escolheram o problema de um disco circular sob
compressão diametral para explicar o método. Um modelo 2D de um disco circular
de 60 mm de diâmetro e 6 mm de espessura foi modelado, e a análise foi feita
usando ABAQUS V6.4. O modelo foi discretizado utilizando elementos quadrilaterais
com oito nós (Figura 32(a)) com 2064 elementos e 6339 nós. As condições de