04.07.2014 Views

download completo - SET - USP

download completo - SET - USP

download completo - SET - USP

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Estudo da localização de deformações associada ao emprego de modelo constitutivo de dano<br />

143<br />

L c<br />

+<br />

⎡1<br />

2 T<br />

T<br />

T<br />

( ΔD) ≡<br />

∫ ( 1 − D) ε Eε<br />

− 2( b1P<br />

ε − k) D − ( b1<br />

( b2<br />

− 1)<br />

P ε − b3<br />

)<br />

∫<br />

3<br />

Ω<br />

2<br />

b b<br />

1<br />

2<br />

P<br />

Ω<br />

T<br />

⎢⎣ 2<br />

3<br />

εD<br />

dΩ +<br />

∫<br />

Ω<br />

c<br />

1<br />

T<br />

2<br />

( ∇D) ∇DdΩ +<br />

∫c2<br />

( ∇ D)<br />

Ω<br />

T<br />

∇<br />

2<br />

DdΩ<br />

D<br />

2<br />

⎤<br />

dΩ<br />

⎥⎦<br />

(29)<br />

para o modelo C apresentado no item anterior, são também solução da fase de correção.<br />

P é o vetor de projeção isótropa ( P T ≡ { 1 1 1 0 0 0 } ) e K é o cone<br />

convexo de funções não-negativas e contínuas em suas primeiras derivadas. Para<br />

soluções fisicamente admissíveis, a solução deve ser restrita a D = D + ΔD<br />

≤1 .<br />

Uma formulação mista alternativa pode ser obtida para o modelo C pela introdução do<br />

gradiente do dano Θ como uma variável independente, adicionando ao funcional a<br />

restrição Θ = ∇D<br />

, ponderada por um mutiplicador de lagrange μ.<br />

Esta formulação requer continuidade C 0 para o campo de dano e é mais conveniente<br />

para a formulação em elementos finitos.<br />

* *<br />

*<br />

0<br />

min max{ L c } ΔD<br />

∈K<br />

; ; K ≡<br />

μ<br />

{ ΔD<br />

∈C<br />

ΔD<br />

≥ 0 em Ω}<br />

ΔD<br />

*<br />

L c<br />

+<br />

⎡1<br />

2 T<br />

T<br />

T<br />

( ΔD,<br />

Θ;<br />

μ) ≡<br />

∫ ( 1−<br />

D) ε Eε<br />

− 2( b1P<br />

ε − k) D − ( b1<br />

( b2<br />

−1)<br />

P ε −b3<br />

)<br />

∫<br />

Ω<br />

2<br />

b1b<br />

2P<br />

3<br />

T<br />

εD<br />

Ω<br />

3<br />

⎢⎣ 2<br />

dΩ +<br />

∫<br />

Ω<br />

c Θ<br />

1<br />

2<br />

dΩ +<br />

∫<br />

Ω<br />

c<br />

2<br />

T<br />

( ∇Θ) ∇ΘdΩ + 2μ( ∇D<br />

− Θ)<br />

∫<br />

Ω<br />

n<br />

dΩ<br />

Uma importante característica da formulação é que o incremento de dano nos problemas<br />

de mínimo 28-29 e de ponto-sela 30 possui a restrição somente de não ser negativo.<br />

Além disso, pode-se provar que as condições de contorno (equações 24 e 25) são<br />

satisfeitas implicitamente para os funcionais acima, conforme demonstram Comi &<br />

Perego 37 para o caso de modelo elastoplástico com gradiente de segunda ordem, e<br />

Comi 11 e Comi & Driemeier 13 para modelo elástico com dano com gradiente segundo e<br />

gradientes segundo e quarto, respectivamente.<br />

D<br />

2<br />

⎤<br />

dΩ<br />

⎥⎦<br />

(30)<br />

3.3 Implementação em elementos finitos com variáveis generalizadas<br />

3.3.1 Fase de previsão<br />

O funcional da fase de previsão discretizado tem a seguinte forma:<br />

min max<br />

δu,<br />

δε<br />

δσ<br />

{ L p } com<br />

δu<br />

= 0 on Γu<br />

Cadernos de Engenharia de Estruturas, São Carlos, v. 10, n. 46, p. 127-155, 2008

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!