download completo - SET - USP
download completo - SET - USP
download completo - SET - USP
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
146<br />
Larissa Driemeier & Sergio Persival Baroncini Proença<br />
Para a fase de correção do modelo C de dano, adaptou-se uma formulação mista (Comi<br />
& Driemeier 13 ), de modo que apenas continuidade C 0 é requerida para aproximação do<br />
campo de dano, mesmo quando termos de quarta ordem estão presentes no funcional.<br />
A formulação variacional mista pode ser obtida introduzindo o gradiente de dano Θ<br />
como uma variável independente e adicionando ao funcional a restrição Θ =∇D ,<br />
ponderada por um multiplicador de Lagrange μ.<br />
O funcional discretizado assume a forma:<br />
min<br />
max<br />
ΔD, Θ μ<br />
+<br />
*<br />
{ } com ΔD<br />
≥ 0<br />
L c<br />
* ⎡1<br />
2 T<br />
2 ⎤<br />
L c ( ΔD,<br />
Θ;<br />
μ) ≡<br />
∫ ( 1−<br />
N DD<br />
) ε Eε<br />
+ 2k<br />
N DD<br />
+ b3<br />
( N DD<br />
) dΩ<br />
(40)<br />
⎢⎣ 2<br />
⎥⎦<br />
+<br />
2<br />
T<br />
( N Θ) dΩ + c ( ∇N<br />
Θ) ( ∇N<br />
Θ)<br />
∫c1<br />
D ∫ 2 D<br />
Ω<br />
Ω<br />
T<br />
∫<br />
2( N μ μ ) ( ∇N<br />
DD<br />
− N DΘ)<br />
Ω<br />
Ω<br />
dΩ<br />
D<br />
dΩ<br />
O problema de complementaridade da forma discreta da fase de correção é dado por:<br />
( ) ( ( ) )( Δ ) ( Δ )<br />
f ≡ϕ ε − a ε + b D + D + d ε, D + D −k − gμ<br />
≤0 (41)<br />
n<br />
ΔD ≥ 0 ,<br />
T<br />
f ΔD=<br />
0<br />
(42)<br />
c1Θ+ c2Θ<br />
− hμ =0<br />
(43)<br />
T T<br />
g D −h<br />
Θ =0<br />
(44)<br />
n<br />
onde:<br />
g ≡<br />
c<br />
2<br />
∫<br />
Ω<br />
≡<br />
( ∇N<br />
)<br />
∫<br />
Ω<br />
c<br />
2<br />
D<br />
T<br />
N<br />
μ<br />
dΩ;<br />
≡<br />
T<br />
( ∇N<br />
Θ<br />
) ∇N<br />
ΘdΩ;<br />
h ≡<br />
∫<br />
c<br />
1<br />
∫<br />
Ω<br />
c N<br />
1<br />
Ω<br />
N<br />
T<br />
Θ<br />
T<br />
Θ<br />
N<br />
N<br />
Θ<br />
μ<br />
dΩ;<br />
dΩ<br />
(45)<br />
Considera-se n D , n Θ e n μ o número de variáveis generalizadas D , Θ e m,<br />
e<br />
respectivamente, considerando todo a malha de elementos finitos e n D , n<br />
e<br />
Θ e n e μ o<br />
número das mesmas variáveis por elemento.<br />
Na modelagem o multiplicador de lagrange μ é descontínuo através dos elementos. Tal<br />
opção permite desacoplar a equação 44 ao nível de elemento, podendo-se condensar<br />
Cadernos de Engenharia de Estruturas, São Carlos, v. 10, n. 46, p. 127-155, 2008