02 Aeroelasticidade Estática - Divergência Da Asa PDF

Aeroelasticidade Esttica Divergncia da Asa
A aeroelasticidade esttica o estudo de deflexo da estrutura flexvel da aeronave

sujeita s cargas aerodinmicas, quando se considera que as cargas e o movimento
no dependem de tempo.
Estas cargas induzem flexo e toro na asa, mudando o ngulo de incidncia e em
consequencia o escoamento aerodinmico, quais em torno mudam as cargas e as
deflexes at se chega ao equilbrio.
Existem dois fenmenos de aeroelasticidade esttica que so crticos:
divergncia ( divergence )
reverso de controle ( control reversal )
Divergncia ocorre quando os momentos causados por cargas aerodinmicas superam
os momentos restauradores da estrutura.
Aeroflio rgido bi-dimensional

Anlise Iterativa
Aeroflio simtrico ( sem arqueamento )
c a corda
k rigidez em toro da mola
ec distancia entre centro aerodinmico e eixo elstico
0 ngulo inicial do aeroflio
ngulo de toro devido s cargas aerodinmicas
a1 inclinao da curva de sustentao
V velocidade do ar (TAS)
O momento de arfagem :
M V 2ca10 ec qec 2 a10

2
q a presso dinmica:
1
V 2
2
A energia potencial ou energia de deformao da mola :
1
K 2
2
Equaes de Lagrange
A equao do Lagrange para um sistema com um grau de liberdade
W
d T T D U
Qx
dt x x x x
x
T a energia cintica
U a energia potencial ou de deformao

D a funo dissipadora
Qx a fora generalizada
W o trabalho da fora generalizada
2
W qec a10
Q
qec 2 a10

A equao do Lagrange neste caso:
U
Q
qec 2 a1 0
K qec a1 0
qR 0
K
2
ec 2 a1
R
K
Devido carga aerodinmica aplicada, o aeroflio torcionou um ngulo , ento o
momento de arfagem muda com o novo ngulo de incidncia.
Primeira iterao
O momento de arfagem com o novo ngulo de incidncia :
M qec 2 a1 0 qR0
A energia de deformao tem a mesma frmula.
K qec 2 a1 0 qR0
1 qR
qec a1
0 qR 1 qR 0
K
2
Continuando do mesmo jeito:

2
qR 1 qR qR 0
Usando a expanso em serie Taylor

2
1
No limite, o ngulo de toro vai ser
qR
0
1 qR
com x 1
nico passo Analise Direta

Considerando o ngulo ainda no determinado como incgnita, o momento de
arfagem :
M qec2 a1 0
A energia de deformao a mesma de antes.
O momento generalizado, baseado no trabalho incremental feito pelo momento

de arfagem por ngulo incremental :
2
W qec a1 0
Q
qec 2 a1 0

Aplicando as equaes de Lagrange para coordenada :
K qec 2 a1 0 K qec 2 a1 qec 2 a10

qec 2 a1
qR
0
0
2
K qec a1
1 qR
O ngulo de toro se tende a infinito quando q tende a 1/R e este define a
velocidade de divergncia.
qdiv
qdiv
K
1
2
R ec a1
2qdiv
1
2
Vdiv Vdiv
2
q / qdiv
0
1 q / qdiv
Asa engastada flexvel

Consideramos uma asa retangular com metade de envergadura s, corda c com
aeroflio simtrico e sem deformao inicial. O eixo elstico ( eixo de flexo )
posicionado a uma distancia ec atrs de centro aerodinmico ( em quarto da corda ).
A rigidez em toro GJ.
A inclinao da curva de sustentao aW e o ngulo de ataque a raiz fixo ( 0 ).
Se assume uma relao linear do ngulo de toro.
y
T
s
Considerando teoria das faixas, a sustentao incremental :
dL qcaW 0 T dy
s
A sustentao total :
y
s
L qcaW 0 T dy qcaW s0 T
s
2
0
s
A energia de deformao em toro :

2
s
d
1
1
GJ 2
T
U GJ
dy
GJ
dy
20
20
2s
s
dy
s
A coordenada generalizada T e o incremento do ngulo em funo do T :
y
T
s
O trabalho feito pelas cargas aerodinmicas determinado considerando o trabalho
do momento de arfagem em cada faixa pelo ngulo de toro da faixa. O trabalho total
determinado por integrao ao longo da semi-envergadura.
W dL ec qcaW 0 T dy ec
s
0
0
s
y
s 0 sT
y
2
qc aW 0 T dye T qec aW
T
s
3
s
2
0
s
Usando a equao do Lagrange com coordenada generalizada T :
U
QT
T
GJT
s s
qec 2 aW 0 T
s
3
2
s 0
s
GJ
2
2
qec
a
qec
a
W
W
T

3
2
s
3q ec 2 s 2 aW
T
6GJ 2q ec 2 s 2 aW
A presso dinmica na divergncia :
3GJ
qW 2 2
ec s aW
A soluo exata da equao diferencial
GJ
2, 47GJ
2 2 2 2
ec s aW 2
ec s aW
2
qW _ exato
Variao da sustentao ao longo da asa engastada
3q ec 2 s 2 aW
dL
y
y
qcaW 0 qcaW 1
0
2 2
dy
s
6GJ 2q ec s aW s
Usando a presso dinmica de divergncia:
3 q / qW y
dL
3q ec 2 s 2aW
y
q caW 1
0 qcaW 1
0
2 2
dy
6GJ 2q ec s aW s
2 1 q / qW s
A sustentao total :
3 q / qW
dL
dy q csaW 1
0
dy
0
4 1 q / qW
s
LTotal
Formulao com matrizes

A matriz de rigidez da viga: [K]. A equao de equilbrio elstico:
[M elast ] [ K ]
O momento de arfagem aerodinmico:
M aero q aW ec 2 y
A aW ec 2 y matriz diagonal
[M elast ] [M aero ] K q A 0
Chegamos a um problema de autovalores:
1
A I 0
q
Exemplo numrico
Um modelo de asa tem rigidez em toro GJ = 23 Nm2
Dados adicionais:
s = 0,5 m ; semi-envergadura
c = 0,15 m; corda
aW = 6 per radiano
= 1,225 Kg/m3 ; densidade do ar ao nvel do mar
e = 0,25; eixo elstico localizado o,25c atrs de a.c. ou localizado na metade da corda
Usando a formula analtica:
3GJ
3 23
2
qW 2 2
8177,8
N/m
ec s aW 0, 25 0,152 0,52 6
VDW
qWexato
2 8177,8
115,5 m/s
1, 225
2qW
2, 47GJ
2, 47 23
2
2 2
6733,0
N/m
ec s aW 0, 25 0,152 0,52 6
VDWexato
2qWexato
2 6733,0
104,8 m/s
1, 225
A soluo com elementos finitos

Velocidade de divergencia
106,00
104,00
102,00
m/s
100,00
98,00
96,00
94,00
92,00
90,00
88,00
1
No. de elementos
10
O efeito do trim sobre a aeroelasticidade esttica

O exemplo apresentado mostra que o aumento da velocidade leva ao aumento do
ngulo de toro e em consequencia o aumento da sustentao. Na prtica a mudana
da velocidade requer o trim da aeronave via o ajuste de profundor para manter o
equilbrio de foras aerodinmicas e de inrcia.
Efeito do trim sobre a velocidade de divergncia de um modelo simples
A aeronave voa horizontalmente em vo uniforme. A fuselagem rgida e pode se

deslocar em sentido vertical e girar. A asa tem flexibilidade em toro. O efeito de
downwash da asa sobre a empenagem horizontal ignorado. As sees da asa e da
empenagem so simtricas.
Cargas
T empuxo ( thrust )
D arrasto ( drag )
W peso ( weight )
LW sustentao de uma asa ( Wing Lift )
LT sustentao da cauda ( Tail Lift )
Coordenadas Generalizadas ( graus de liberdade )

z deslocamento vertical
0 ngulo de arfagem da aeronave
T ngulo de toro na ponta da asa
Equaes de equilbrio da aeronave rgida:
F
F
0T D 0
0 2 LW W LT 0
2q csaW 0 T
2
W LT 0
(1)
lW LT lT 0
(2)
M 0
T
2q csaW 0
2
A energia elstica de deformao:

2
d
1
GJ 2
T
U 2 GJ
T
dy GJ dy
20
s
s
dy
0
s
U
2GJ
O trabalho virtual duas vezes o anterior ( da asa engastada ), e a terceira equao :
GJ

T 2qec 2 saW 0 T
s
2 3
(3)
Podemos eliminar LT usando as equaes (1) e (2), obtendo a seguinte equao:
2q csaW 0 lW lT WlT
2
Arrumando as equaes (4) e (3) em forma matricial:
2q csaW
q ec 2 saW
WlT
0 l l
GJ W T
2
2
q ec saW 2
T 0
s
3
q csaW
(4)
As solues so:
WlT / lW lT
WlT / lW lT
T
4GJ / (ecs) q cs aW / 3 4GJ

q
1
ecs 4qW
6GJ 2q ec s a Wl l l Wl l l 1 q / q
2q cs a 1 q / 4q
q cs a 12GJ q ec s a
2 2
2 2
A presso dinmica na divergncia :
qA 4qW
Efeito do trim vertical sobre a variao da sustentao

ao longo da asa
Substituindo as expresses achadas para 0 e T , a derivada da sustentao :
dL WlT / lT lW 2 q / qW 3 y / s 2
dy
4s 1 q / 4qW
Efeito do trim sobre a sustentao da asa e da

empenagem horizontal
Para manter o trim, o equilbrio de foras e momentos deve ser mantido.
No caso de seo simtrica da asa, a sustentao da asa e da cauda permanece
constante durante mudanas de velocidade.
No caso de uma seo arqueada, vai aparecer um momento de arfagem e as foras
vo mudar com a velocidade, mas ainda vo ficar em equilbrio.
Efeito de enflechamento sobre a aeroelasticidade

esttica
A maioria das aeronaves tem asa com enflechamento positivo ( swept-back ). O motivo
principal aerodinmico ( reduo do arrasto ). Um efeito semelhante pode ser obtido
com enflechamento negativo ( swept-forward ). O enflechamento negativo tem um
efeito prejudicial sobre a divergncia.
Efeito de enflechamento da asa sobre o ngulo efetivo de incidncia
A flexo tem uma influencia maior que a toro sobre o ngulo de incidncia, mas o
modo de divergncia ainda ocorre em toro.
a.
b.
c.
Para caso sem enflechamento (AC), o ngulo de incidncia no muda por causa da
flexo
Para caso com enflechamento positivo ( AD), o ngulo de incidncia efetivo
diminui visto que o ponto D se levanta mais que o A, em flexo.
Para caso com enflechamento negativo (AB), o ngulo de incidncia efetivo
aumenta visto que o ponto A se levanta mais que o B, em flexo.
Em consequencia, para asa com enflechamento negativo, a velocidade de divergncia

diminui em comparao asa sem enflechamento, enquanto que a asa com
enflechamento positivo tem um aumento da velocidade de divergncia.
ngulo de incidncia efetivo devido ao

batimento e arfagem
Consideramos uma asa no afilada com corda c, semi-envergadura s, ngulo de
enflechamento ( positivo atrs ) e ngulo de incidncia 0 na raiz, numa corrente do
ar com velocidade V.
Para ter um modelo mais realista temos que considerar a flexo junto com a toro.
Duas molas de rotao: K que controla o batimento e K que contm a toro.
O batimento e a toro ocorrem ao longo e em torno do eixo elstico y.

De novo vamos usar a teoria das faixas para modelo aerodinmico.
c cos
pitch
cos
c
flap
y ' c sin y ' sin

c
A arfagem diminui o ngulo

efetivo de incidncia.
Batimento para baixo aumenta o

ngulo de incidncia de uma asa
com enflechamento positivo.
No vo com velocidade constante e nivelado a asa deflete para cima ( negativo ) e

em consequencia enflechamento positivo diminui o ngulo de incidncia efetivo.
Em realidade, as deflexes devidas flexo influenciam o ngulo de incidncia efetivo
mais que a toro.
Efeito do ngulo de enflechamento sobre a

velocidade de divergncia
Considerando o batimento e a toro junto:
dL qaW cdy 0 cos sin

O trabalho virtual:
c sin
y
W qaW cdy 0 cos sin

4
cos
0
c cos
qaW cdy 0 cos sin
4
cs cos
s2
cs sin
qaW cdy 0 cos sin

4
2 cos
4
s
A energia de deformao das duas molas :
1
1
2
U K K 2
2
2
Aplicando as equaes de Lagrange para as duas coordenadas generalizadas: e ,
i.e.
Obtivemos duas equaes:
cs 2
c 2 s sin
K qaW 0 cos sin
2
cos
c 2 s cos
K qaW
0 cos sin
4
Arranjando em formato matricial :
s 2 tan cs sin 2
s 2 cs sin cos
K q aW c
q aW c
2
4
2
4

2
2
2
q aW sc sin cos
q aW sc cos
4
4
s 2 cs sin cos
q aW s
2
4
0
2
2
q aW sc cos
4
A divergncia ocorre quando o determinante igual a zero.

s 2 tan cs sin 2
q aW sc 2 cos 2
K q aW c
K
2
4
4
q aW
sin
cos
s
cs sin cos
3
sc

0
4
4
2
A velocidade de divergncia :
Vdiv
2 K K
aW K sc 2 cos2 / 4 K cs 2 tan / 2 c 2 sin 2 / 4
Para asa sem enflechamento ( = 0 ) a velocidade de divergncia :
Vdiv
8K
aW c 2 s
A velocidade da divergncia normalizada com respeito velocidade de asa

sem enflechamento:
Efeito do enflechamento sobre a distribuio da

sustentao
No caso analisado, o efeito do enflechamento sobre a alterao da sustentao para
cada faixa o mesmo.
No grfico abaixo, a sustentao por unidade de envergadura, normalizada com
respeito ao caso = 0 apresentada em funo do ngulo de enflechamento.

02 Aeroelasticidade Estática - Divergência Da Asa PDF

Enviado por

Direitos autorais:

Formatos disponíveis

02 Aeroelasticidade Estática - Divergência Da Asa PDF

Enviado por

Dados do documento

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

02 Aeroelasticidade Estática - Divergência Da Asa PDF

Enviado por

Direitos autorais:

Formatos disponíveis

Aeroelasticidade Esttica Divergncia da Asa

A aeroelasticidade esttica o estudo de deflexo da estrutura flexvel da aeronave

Aeroflio rgido bi-dimensional

Aeroflio simtrico ( sem arqueamento )

M V 2ca10 ec qec 2 a10

A energia potencial ou energia de deformao da mola :

U a energia potencial ou de deformao

A equao do Lagrange neste caso:

momento de arfagem muda com o novo ngulo de incidncia.

Continuando do mesmo jeito:

Usando a expanso em serie Taylor

No limite, o ngulo de toro vai ser

nico passo Analise Direta

O momento generalizado, baseado no trabalho incremental feito pelo momento

Aplicando as equaes de Lagrange para coordenada :

K qec 2 a1 0 K qec 2 a1 qec 2 a10

Asa engastada flexvel

Considerando teoria das faixas, a sustentao incremental :

A energia de deformao em toro :

A coordenada generalizada T e o incremento do ngulo em funo do T :

Usando a equao do Lagrange com coordenada generalizada T :

A presso dinmica na divergncia :

Variao da sustentao ao longo da asa engastada

Formulao com matrizes

Chegamos a um problema de autovalores:

Usando a formula analtica:

A soluo com elementos finitos

O efeito do trim sobre a aeroelasticidade esttica

Efeito do trim sobre a velocidade de divergncia de um modelo simples

A aeronave voa horizontalmente em vo uniforme. A fuselagem rgida e pode se

Coordenadas Generalizadas ( graus de liberdade )

Equaes de equilbrio da aeronave rgida:

A energia elstica de deformao:

O trabalho virtual duas vezes o anterior ( da asa engastada ), e a terceira equao :

Podemos eliminar LT usando as equaes (1) e (2), obtendo a seguinte equao:

Arrumando as equaes (4) e (3) em forma matricial:

4GJ / (ecs) q cs aW / 3 4GJ

A presso dinmica na divergncia :

Efeito do trim vertical sobre a variao da sustentao

Efeito do trim sobre a sustentao da asa e da

Efeito de enflechamento sobre a aeroelasticidade

Efeito de enflechamento da asa sobre o ngulo efetivo de incidncia

Em consequencia, para asa com enflechamento negativo, a velocidade de divergncia

ngulo de incidncia efetivo devido ao

O batimento e a toro ocorrem ao longo e em torno do eixo elstico y.

y ' c sin y ' sin

A arfagem diminui o ngulo

Batimento para baixo aumenta o

No vo com velocidade constante e nivelado a asa deflete para cima ( negativo ) e

Efeito do ngulo de enflechamento sobre a

dL qaW cdy 0 cos sin

A energia de deformao das duas molas :

Obtivemos duas equaes:

Arranjando em formato matricial :

aW K sc 2 cos2 / 4 K cs 2 tan / 2 c 2 sin 2 / 4

Para asa sem enflechamento ( = 0 ) a velocidade de divergncia :

A velocidade da divergncia normalizada com respeito velocidade de asa

Efeito do enflechamento sobre a distribuio da

Você também pode gostar