MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade...

34
1 Silvestre T Pinho Mestrado em Engenharia Mecânica Mecânica não linear Dra. Lúcia Dinis Viscoelasticidade 08 de Novembro de 2005 Conclusões 08 de 08 de Novembro Novembro de 2005 de 2005 2 Introdução Elementos Maxwell Kelvin-Voigt + Complexos Hiper- elasticidade Objectivos Entender o que é viscosidade e visco- elasticidade Conceitos de fluência e relaxação de tensões Conhecer os elementos viscoelásticos básicos (Maxwell e Kelvin-Voigt) Reconhecer elementos básicos em modelos mais complexos Modelos viscoelásticos para grandes deformações (visco-hiperelasticidade)

Transcript of MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade...

Page 1: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

11

Silvestre T Pinho

Mestrado em Engenharia MecânicaMecânica não linear

Dra. Lúcia Dinis

Viscoelasticidade

08 de Novembro de 2005

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 22

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Objectivos

• Entender o que é viscosidade e visco-elasticidade

• Conceitos de fluência e relaxação de tensões

• Conhecer os elementos viscoelásticos básicos (Maxwell e Kelvin-Voigt)

• Reconhecer elementos básicos em modelos mais complexos

• Modelos viscoelásticos para grandes deformações (visco-hiperelasticidade)

Page 2: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 33

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Sumário da aula

• Introdução

• Elementos e modelos viscoelásticos

•Elemento de Maxwell

•Elemento de Kelvin-Voigt

•Modelos mais complexos

• Modelos materiais visco-hiperelásticos

• Recapitulação e conclusões

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 44

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Alguns conceitos e modelos

ViscosidadeViscosidade

••A viscosidade A viscosidade éé uma quantidade uma quantidade que descreve a resistência que que descreve a resistência que um material oferece ao um material oferece ao escoaescoa--mento.mento. O material O material éé forforççado a ado a escoarescoar por uma forpor uma forçça a FF..

••Quando a forQuando a forçça a éé retirada, o retirada, o material para de material para de escoar.escoar.

••Este fenEste fenóómeno meno éé representado representado pelo elemento viscoso, de acordo pelo elemento viscoso, de acordo com a lei de Newton da com a lei de Newton da viscosidade:viscosidade:

FF

FF

σ η γ=.

Page 3: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 55

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Alguns conceitos e modelos

ElasticidadeElasticidade

••Aos materiais que não Aos materiais que não escorrem, mas que mostram escorrem, mas que mostram propriedades tpropriedades tíípicas de spicas de sóólidos, lidos, éé usada uma mola para usada uma mola para representrepresentáá--loslos

••A mola A mola éé caracterizada pela lei caracterizada pela lei de de HookHook

••Quando a forQuando a forçça a FF éé retirada, a retirada, a mola retorna mola retorna àà configuraconfiguraçção ão inicialinicial

FF

FF

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 66

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Alguns conceitos e modelos

ViscoelasticidadeViscoelasticidade

••Um material Um material viscoelviscoeláásticostico apresenta apresenta simultaneamente caractersimultaneamente caracteríísticas viscosas e sticas viscosas e eleláásticassticas

••A A viscoelasticidadeviscoelasticidade depende da escala de depende da escala de tempo da experiência, quando comparada tempo da experiência, quando comparada com a escala de tempo do material (tempo com a escala de tempo do material (tempo de relaxade relaxaçção)ão)

Page 4: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 77

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Materiais viscoelásticos

••Não sNão sóó os materiais usados em os materiais usados em engenharia (engenharia (egeg polpolíímeros) são meros) são viscoelviscoeláásticossticos

(Mow 1991)

FromFromBioE515Lecture2.pptBioE515Lecture2.ppt

Osso: material poroOsso: material poro--viscovisco--eleláástico stico

anisotranisotróópico gradativopico gradativo

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 88

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Materiais viscoelásticos

ε

σ

1ε&12 εε && >

23 εε && >34 εε && >↑ε&••As tensões As tensões

dependem das dependem das deformadeformaçções e ões e das taxas de das taxas de defomadefomaççãoão

Page 5: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 99

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Materiais viscoelásticos – relaxação de tensões

ε

Tempo

T

Input

>

≤≤=Tt

TtT

t

0

0 0

ε

εε

Output

Tempo

Tensão

T

30 εε =

20 εε =

10 εε =

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 1010

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Materiais viscoelásticos – fluência

Tensão

Tempo

T

≤≤=Tt

TtT

t

0

0 0

σ

σσ

Input

Tempo

ε

T

10 σσ =

30 σσ =

20 σσ =

Output

Page 6: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 1111

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Diferentes elementos

ElasticidadeElasticidade

Em elasticidade, a Em elasticidade, a mola mola éé a unidade a unidade principal, e assumeprincipal, e assume--se se que responde que responde instantaneamente e instantaneamente e de forma reversde forma reversíível a vel a cargas e deformacargas e deformaçções ões aplicadasaplicadas

MolaMola

εσ E=

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 1212

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

••A unidade viscosa fundamental A unidade viscosa fundamental ééo amortecedor, onde a tensão o amortecedor, onde a tensão ééassumida ser uma funassumida ser uma funçção da taxa ão da taxa de deformade deformaçção e não da ão e não da deformadeformaççãoão

••O amortecedor O amortecedor éé normalmente normalmente assumido ser linear, e a constante assumido ser linear, e a constante de proporcionalidade de proporcionalidade éé chamada chamada viscosidadeviscosidade

ViscosidadeViscosidade

AmortecedorAmortecedor

( )εσ &f=

εησ &=

Diferentes elementos

DashpotDashpot

Page 7: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 1313

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Um ponto chave em modelar o comportamento viscoelástico de diferentes materiais está em escolher a forma adequada das componenteselástica e viscosa (eg linear) bem comocombinar os elementos no melhorsistema possível por forma a que o comportamento em função do tempo seja previsto adequadamente

Diferentes elementos

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 1414

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elementos em sElementos em séérie:rie:

εεε

σσσ

=+

==

21

21

εσ ,

22 ,εσ11,εσ

1 2 0== εσ

1 2

Diferentes elementos

Page 8: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 1515

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elementos em paralelo:Elementos em paralelo:

1

2

1

2

11,εσ

22 ,εσ

εσ ,

εεε

σσσ

==

=+

21

21

Diferentes elementos

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 1616

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elemento de Maxwell

1 2

22 εσ E=11 εησ &=

12

21

21 εηεεεε

σσσ&=⇒

=+

==E

•• O elemento de O elemento de MaxwellMaxwell éé o modelo o modelo bbáásico de um flusico de um fluíído do ((flufluíídodo, porque a , porque a resposta a longo prazo resposta a longo prazo correspondente a uma correspondente a uma tensão constante tensão constante aplicada aplicada éé idêntica idêntica àà de de um um flufluíídodo――ieie continua a continua a deformardeformar--se)se)

Page 9: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 1717

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Tempo

Def

orm

ação 0ε

( )tH0εε =

Estado inicial:

Respostainstantânea:

Estado final:

Elemento de Maxwell - relaxação de tensões

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 1818

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

1 2

22 εσ E=11 εησ &=

12

21

21 εηεεεε

σσσ&=⇒

=+

==E ( )tHEE 011 εεεη =+&

Tempo

0ε( )tH0εε =

Elemento de Maxwell - relaxação de tensões

Page 10: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 1919

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elemento de Maxwell - relaxação de tensões

( )tHEE 011 εεεη =+&Tempo

0ε( )tH0εε =

( ) ηεσ EteEt

−= 0

Deformação viscosa:

Tensões:

( ) ( )ηεε Etet

−−= 101

Tempo

Ten

são

0εE

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 2020

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elemento de Maxwell – fluência

Tempo

Ten

são

Estado inicial

Respostainstantânea

Estado final

0σ( )tH0σσ =

Page 11: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 2121

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elemento de Maxwell – fluência

Tempo

Ten

são 0σ

( )tH0σσ =

( )( ) 20

10

εσ

εησ

EtH

tH

=

= &

1 2

22 εσ E=11 εησ &=

12

21

21 εηεεεε

σσσ&=⇒

=+

==E

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 2222

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elemento de Maxwell – fluência

Tempo

0σ( )tH0σσ =

( ) ( )η

σσε

ttH

Et 00 +=

Def. elástica:

( ) ( )tHE

t 02

σε =

Def. viscosa:

Def. total:

( )η

σε

tt 0

1 =

Tempo

E

def

orm

ação

Page 12: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 2323

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elemento de Kelvin-Voigt

•• O elemento de O elemento de KelvinKelvin--VoigtVoigt éé o modelo o modelo bbáásico de um ssico de um sóólido lido ((ssóólidolido, porque a , porque a resposta a longo prazo resposta a longo prazo correspondente a uma correspondente a uma tensão constante tensão constante aplicada aplicada éé idêntica idêntica àà de de um um ssóólidolido――ieie deixa de deixa de deformardeformar--se)se)

σεεηεεε

σσσ=+⇒

==

=+E&

21

21

2

22 εσ E=

11 εησ &=

1

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 2424

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elemento de Kelvin-Voigt – relaxação de tensões

Tempo

Def

orm

ação

Estado inicial:

Respostainstantânea:

Estado final:

0ε( )tH0εε =

Page 13: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 2525

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elemento de Kelvin-Voigt – relaxação de tensões

Tempo

Def

orm

ação

0ε( )tH0εε =

σεεηεεε

σσσ=+⇒

==

=+E&

21

21

2

22 εσ E=

11 εησ &=

1

( ) ( ) σεδηε =+ tHEt 00

0εE

Ten

são

Tempo

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 2626

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elemento de Kelvin-Voigt – fluência

Tempo

Ten

são 0σ

( )tH0σσ =

Estado inicial:

Respostainstantânea:

Estado final:

Page 14: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 2727

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elemento de Kelvin-Voigt – fluência

Tempo

Ten

são 0σ

( )tH0σσ =

σεεηεεε

σσσ=+⇒

==

=+E&

21

21

2

22 εσ E=

11 εησ &=

1

( )tHE 0σεεη =+&

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 2828

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Elemento de Kelvin-Voigt – fluência

Tempo

0σ( )tH0σσ =( )tHE 0σεεη =+&

( ) ( )ησε Et

eE

t−−= 10

Soluçãohomogénea: ( ) ηε Et

cet−=

Solução particular:( )

Et 0σ

ε =

Solução total com c.i.

( )( ):00 ==tεTempo

Def.

E

Page 15: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 2929

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos mais complexos

•• NemNem o o modelomodelo de Maxwell de Maxwell nemnem o o modelomodelo de Kelvinde Kelvin--Voigt Voigt produzproduzumauma respostaresposta viscoelviscoeláásticastica quequecapture capture qualitativamentequalitativamente o o comportamentocomportamento de de muitosmuitosmateriaismateriais reaisreais; ;

•• PortantoPortanto, , modelosmodelos maismaiscomplexoscomplexos devemdevem ser ser usadosusados

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 3030

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos mais complexos: 3 elementos

MaxwellMaxwell em paralelo com uma molaem paralelo com uma mola

εεεηεεε

εεηεεσ1111

21

21221EE

EEE=+⇒

+=

+=+=&

&

1ε 2ε

η

2E

1E

Page 16: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 3131

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos mais complexos: 3 elementos

1ε 2ε

η

1E2E

KelvinKelvin--VoigtVoigt em sem séérie com uma rie com uma mola (mola (three three

parameter solidparameter solid))

( ) εεηεεεε

εεηεσ21121

21

22111EEE

EE=++⇒

+=

=+=&

&

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 3232

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Kelvin-Voigt em série com uma mola –relaxação de tensões

Tempo

ε0ε

( )tH0εε =

η

1E2E

Estado inicial:

Respostainstantânea:

Estado final:

Page 17: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 3333

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Kelvin-Voigt em série com uma mola –relaxação de tensões

Tempo

ε0ε

( )tH0εε =1ε 2ε

η

1E2E

( ) εεηεεεε

εεηεσ21121

21

22111EEE

EE=++⇒

+=

=+=&

&

( ) ( )tHEEEE 0221121 εεεηε ==++ &

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 3434

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Kelvin-Voigt em série com uma mola –relaxação de tensões

Tempo

ε0ε

( )tH0εε =( ) ( )tHEEEE 0221121 εεεηε ==++ &

Elemento elástico

( ) ( ) ( )[ ]ηεεε tEE

eEE

EtHt 211

21

0202

+−−+

−=

Tensão total

( ) ( ) ( )[ ]ηεεεσ tEE

eEE

EEtEt 211

21

0

2

20222

+−−+

−==

Tempo

Tensão

02εE21

021

EE

EE

+

ε

( ) ( )[ ]ηεε tEE

eEE

Et 211

21

021

+−−+

=

Elemento de Kelvin

Page 18: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 3535

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Kelvin-Voigt em série com uma mola –fluência

Tempo

Ten

são

0σ( )tH0σσ =

η

1E2E

Estado inicial:

Respostainstantânea:

Estado final:

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 3636

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Kelvin-Voigt em série com uma mola –fluência

( ) εεηεεεε

εεηεσ21121

21

22111EEE

EE=++⇒

+=

=+=&

&

( )21

221110

εεε

εεηεσσ

+=

=+== EEtH &

1ε 2ε

η

2E

Tempo

Ten

são

0σ( )tH0σσ =

Page 19: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 3737

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Kelvin-Voigt em série com uma mola –fluência

Elemento de Kelvin

( ) ( ) ( )ησεεηεσ tE

eE

tEtH 111

011110

−−=⇒+= &

Elemento elástico

( ) ( ) ( )2

02220

E

tHtEtH

σεεσ =⇒=

( )21

221110

εεε

εεηεσσ

+=

=+== EEtH &

Defomação total

( )

+

−=

21

0

11 1

EE

et

tE η

σε

Tempo

ε

2

0

E

σ

+

21

0

11

EEσ

Tempo

Ten

são

0σ( )tH0σσ =

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 3838

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Kelvin-Voigt em série com uma mola –funções de relaxação e de fluência

A A funfunççãoão de de relaxarelaxaççãoão de de tensõestensões

e a e a funfunççãoão de de fluênciafluência sãosãodeterminadasdeterminadas tomandotomando um valor um valor unitunitááriorio dada varivariáávelvel de de entradaentrada(input), (input), sejaseja estaesta tensãotensão ououdeformadeformaççãoão

Page 20: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 3939

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Kelvin-Voigt em série com uma mola –funções de relaxação e de fluência

( ) ( )[ ]21

0

2

202 1 21e

EE

EEt

tEE ⇒−+

−= +−εεσ η

FunFunçção de relaxaão de relaxaççãoão

( )211

21

21

21 ,EEE

eEE

EE

EEtG

t

+=

++

=− η

ττ

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 4040

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Kelvin-Voigt em série com uma mola –funções de relaxação e de fluência

( )21

0

11 1

EE

et

tE

σεη

+

−=

FunFunçção de fluênciaão de fluência

( )112

,11

EE

e

EtJ

t ητ

τ

=−

+=−

Page 21: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 4141

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 4242

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Page 22: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 4343

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Comparação - Fluência

Tempo

Ten

são

0σ( )tH0σσ =

Input

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 4444

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Comparação - Fluência

TempoKelvinKelvin--VoigtVoigt

ε

E

KelvinKelvin--VoigtVoigt em em sséérie com uma rie com uma mola mola Tempo

ε

2

0

E

σ

MaxwellMaxwell

Tempo

E

0σ ε

Output

Page 23: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 4545

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Comparação – Relaxação de tensões

Tempo

ε0ε

( )tH0εε =

Input

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 4646

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Comparação – Relaxação de tensões

MaxwellMaxwell

Tempo

Tensão0εE

KelvinKelvin--VoigtVoigt 0εE

Tensão

Tempo

KelvinKelvin--VoigtVoigt em em sséérie com uma rie com uma mola mola Tempo

Tensão

02εE

Output

Page 24: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 4747

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Modelos constitutivos com variModelos constitutivos com variááveis internasveis internas

• Muitos materiais usados em engenharia são inelásticos

• Os modelos hiperelásticos estudados na última aula não são adequados para estudar esses materiais, em que pode haver dissipação

• No restante desta aula, vamos estudar materiais inelásticos com base no conceito de variáveis internas

• Aplicações: Modelação de dano e materiais viscoelásticos

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 4848

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Conceito de variConceito de variááveis internasveis internas

• O estado termodinâmico corrente (com T=cte) de um

material (hiper)elástico pode ser determinado com base apenas no gradiente de deformação F

• O gradiente de deformação pode ser medido e por • O gradiente de deformação F pode ser medido e por isso é chamado variável externa

•isso é chamado variável externa

• O estado termodinâmico corrente de materiais que involvam dissipação pode ser determinado a partir de um número de variáveis internas

Page 25: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 4949

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Conceito de variConceito de variááveis internasveis internas

• Essas variáveis descrevem aspectos da estrutura interna de materiais, associados com efeitos irreversíveis (dissipação)

•• As tensões e deformações dependem das variáveis internas internas

• Portanto, o conceito de variável interna postula que o estado termodinâmico corrente num ponto de um material dissipativo é definido pelo gradiente de deformação F (variável externa; e também eventualmente a temperatura) e variáveis internas ξ

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 5050

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Conceito de variConceito de variááveis internasveis internas

• Ou seja, o estado termodinâmico pode ser entendido como um estado fictício de equilíbrio termodinâmico

•como um estado fictício de equilíbrio termodinâmico

• O estado termodinâmico corrente é descrito pelos valores actuais (e não pela história) do gradiante de deformação (e temperatura), bem como de um determinado número de variáveis internas

• A natureza das variáveis internas pode ser física • A natureza das variáveis internas pode ser física (mecânica, térmica, química ou eléctrica), descrevendo processos microscópicos sob um ponto de vista macroscópico

Page 26: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 5151

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

EquaEquaçções constitutivasões constitutivas

• Energia livre de Helmholtz para caso isotérmico (energia de deformação): ( )mξξF ,...,, 1Ψ=Ψ

• Desigualdade de Clausius-Plank (Clausius-Duhem): ( ) 0,...,, 1intint ≥Ψ−= mwD ξξF&

( ) ∑= ∂

Ψ∂+

∂Ψ∂

=Ψm

m

1

1 ::,...,,α

αα

ξξ

FF

ξξF &&&

( )( )( ) ( )

:,...,,

:,...,,

,...,,

1

11

1intint

Ψ∂−

Ψ∂−=

=Ψ−=

∑=

mmm

mwD

αα

α

ξξ

ξξFF

F

ξξFP

ξξF

&&

&

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 5252

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

EquaEquaçções constitutivasões constitutivas

• Equações constitutivas ( )

F

ξξFP

Ψ∂= m,...,, 1

•F∂

• Dissipação:

0:1

int ≥Ξ=∑=

m

αα ξ&

• Equações constitutivas internas • Equações constitutivas internas ( )

αα

ξ

ξξF

Ψ∂−=Ξ m,...,, 1

Page 27: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 5353

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

EquaEquaçções constitutivasões constitutivas

Equações de evolução e equilíbrio termodinâmico •

Equações de evolução e equilíbrio termodinâmico • É necessário complementar as equações

constitutivas + equações constitutivas internas com uma relação cinemática que descreva como as variáveis αξ variam no tempo―equações de evolução:

( )( )mAt ξξFξ ,...,,)( 1αα =&

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 5454

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Materiais viscoelásticos

• Postulado:

( ) ∑=

∞∞ ΓΥ+Ψ+Ψ=ΓΓΨm

m J1

isovol1 ),()()(,...,,α

αα CCC

com com

0),(,0)(,0)1(1

isovol =Υ=Ψ=Ψ ∑=

∞∞m

αα ΙII

)();( isovol C∞∞ ΨΨ J � Resposta volumétrica e

isocórica quando ∞→t

∑=

ΓΥm

1

),(α

αα C � Potencial dissipativo =1α

αΥ � Energia livre configuracional

Page 28: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 5555

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Materiais viscoelásticos• Cada variável interna αΓ caracteriza o

comportamento viscoso do material • Γ• αΓ � medida de deformações inelásticas

• Γ• αΓ � C • Para um material hiperelástico (em regime

isotérmico), o estado de equilíbrio termodinâmico é definido por C .

• Para um material viscoelástico (em regime

isotérmico), o estado (fictício) de equilíbrio termodinâmico é definido por ( )mΓΓΓ ,...,,, 21C

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 5656

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Materiais viscoelásticos

Tensões Segundas tensões de Piola-Kirchoff

( )C

CS

∂ΓΓΨ∂

= m,...,,2 1

Dissipação

02

1:

),(2

1

int ≥ΓΓ∂

ΓΥ∂−= ∑

α α

αα &m

DC

Page 29: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 5757

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Materiais viscoelásticos

Decomposição das segundas tensões de Piola-Kirchoff:

( )=

ΓΓΨ∂= m1,...,,

2C

CS

∑=

∞∞

=∂

ΓΥ+Ψ+Ψ∂

m

J1

isovol ),()()(

αα

C

CC

=

∞∞

ΓΥ∂+

∞∂

Ψ∂+

∞∂

Ψ∂=

mJ

1

iso

iso

vol

vol ),(2)(2

)(2

α

α

αα

44344214342143421Q

C

C

S

C

C

S

C

• αQ � tensões de desiquilíbrio

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 5858

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Materiais viscoelásticos

• É possível mostrar que as tensões de desiquilíbrio

αQ são conjugadas energéticas de αΓ , e que as equações constitutivas internas se podem escrever:

α

ααα Γ∂

ΓΥ∂−=

),(2

CQ

Page 30: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 5959

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Exemplo

ε

1q

0>∞E

01 >E

0>mE0>mη

01 >η

mq

mγ0>=

α

αα

ητ

E

σ

qm

1

=

+==

∞ ∑

εσ

σσα

α

E

==

= ∞∞

αγη

εσ

mq ,...,1, == αγη ααα &

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 6060

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Exemplo

αααα

α γηεσσσ &==+= ∞∞=

∞ ∑ qEqm

,,1

( )ααα γε −= EqTemos também

Derivando a expressão anterior: ( )ααα γε &&& −= Eq

O que leva a:•

=+ ετ α

α

αα E

qq&

E a dissipação é: ( ) 02

11

int ≥== ∑∑==

αα

ααα

α γηγ &&mm

qD

Page 31: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 6161

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Exemplo

( ) ( ) ( )∑=

∞ +=m

m

1

1 ,,...,,α

αα γευεψγγεψ

A energia de deformação pode ser definida como:

Com: ( ) ( ) ( )22

2

1,e

2

1ααα γεγευεεψ −== ∞∞ EE

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 6262

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Exemplo

( ) ( ) ( )

( )

σ

σ

γεεε

γευε

εψε

γγεψ

α

αα

αα

=

+=

=−+=

=∂

−∂+=

q

EE

m

d

d,...,, 1

( ) ( ) ( )22

2

1,e

2

1ααα γεγευεεψ −== ∞∞ EE

A derivada de ψ em ordem à deformação ε permite obter a tensão σ:

Page 32: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 6363

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Modelos visco-hiperelásticos

Exemplo

Derivando ψ relativamente às variáveis internas γα, obtem-se:

( ) ( ) αααα

αα γεγ

γευqE =−=

∂−

,

E a dissipação pode ser expressa como:

αα α

α γγυ

&∑= ∂

∂−=

m

D1

int

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 6464

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Recapitulação

• Vimos o que é viscosidade e visco-elasticidade

• Conceitos de fluência e relaxação de tensões

• Vimos quais os elementos viscoelásticosbásicos (Maxwell e Kelvin-Voigt)

• Vimos modelos complexos como combinação de elementos básicos

• Modelos viscoelásticos para grandes deformações (visco-hiperelasticidade)

Page 33: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 6565

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Conclusões

• Vários materiais reais têm comportamento viscoelástico

• De um modo geral, o modelo mais adequado para cada material―o qual pode ser bastante complexo―pode ser obtido por combinação de elementos simples

• Para grandes deformações, é adequado usar uma formulação baseada em hiperelasticidade

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 6666

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Referências

• G A Holzapfel. Nonlinear Solid Mechanics, A continuumapproach for enginneers. John Wiley & Sons Ldt, England, 2000.

• L E Malvern. Introduction to the Mechanics of a ContinuousMedium. Prentice-Hall, Inc, USA, 1969.

• I Doghri. Mechanics of Deformable Solids, Linear and nonlinear, analytical and computational aspects. Springer-Verlag, Germany, 2000.

• V A Lubarda. Elastoplasticity theory. CRC Press LLC, 2002

• G T Mase, G E Mase. Continuum Mechanics for Engineers, Second Edition, CRC Press LLC, 1999

• J Bischoff, Viscoelasticity―Constitutive response of tissue at large strain rates. Esci 274 Mechanics of Biomaterials, The University of Auckland, 2002

Page 34: MestradoemEngenhariaMecânica Mecânicanãolinearstpinho/teaching/feup/y0506/...Viscoelasticidade •Um material viscoelástico apresenta simultaneamente caracteríísticas viscosas

Conclusões

08 de 08 de NovembroNovembro de 2005de 2005 6767

Introdução

Elementos

Maxwell

Kelvin-Voigt

+ Complexos

Hiper-elasticidade

Apontamentos

• Esta apresentação (em formato pdf) encontra-se em:

http://www.fe.up.pt/~ldinis

e em

http://www.fe.up.pt/~stpinho

• Perguntas?