Universidade Federal de Santa Catarina
Centro Tecnolgico
Departamento de Engenharia Mecnica
EMC5352 Introduo ao Projeto Aeronutico
Agro Aeronave
Professor: Juliano T. Micheletto
Alunos: Chaiane Wiggers de Souza
Fbio do Monte Sena (capito)
Fernando Jnior Groff
Guilherme Tozzi Marques
Marco Aurlio Stimamiglio
Timmermann
Florianpolis, dezembro de 2010.
Introduo
O presente trabalho constitui-se como finalizao da disciplina de Introduo ao
Projeto Aeronutico. No incio do semestre, o desafio foi proposto como projetar uma
aeronave dentre as opes oferecidas pelo professor, ou outras sugeridas pelo grupo. No
nosso caso, escolhemos uma aeronave agrcola, cujo uso seria destinado a aplicao de
inseticidas e pesticidas.
O projeto era desenvolvido medida que os checks eram seguidos. No final, o
contedo foi sintetizado, resultando na aeronave projetada.
Desenvolvimento
Primeiramente, a equipe escolheu o desafio: construir uma aeronave agrcola. Tal escolha
deu-se ao perceber a falta de aeronaves agrcolas no Brasil, ou seja, uma possibilidade de
mercado. O capito foi escolhido por demonstrar pr-atividade, Fbio Sena. Em seguida,
entendendo o nosso pblico-alvo, escolhemos trs diretrizes de projeto:
Custo: o consumidor deseja uma boa relao custo-benefcio. Uma aeronave muito
cara no atende a esse requisito.
Segurana: o avio voar em baixos altitudes, perto de fios eltricos e construes,
por exemplo. O risco de acidentes alto. A proteo ao piloto fundamental.
Manuteno: o custo para manter a aeronave em operao deve ser baixo para que o
agricultor possa maximizar o lucro.
Para melhor basear nosso projeto, fizemos uma pesquisa de algumas aeronaves com uso
agrcola, com destaque Ipanema (Neiva), um avio totalmente brasileiro movido a lcool. As
aeronaves escolhidas foram:
Neiva Ipanema
American Champion Citabria - 7GCAA Citabria "A" Package
Air Tractor AT502-B
Aero Boero AB-115
Piper Cherokee PA-28-140
Socata TB-9
Tabela 1: Anlise histrica
Fabrica
nte
NEIVA
Nome
Ipanema (lcool)
Enverga
dura
(m)
Altura
(m)
Compri
mento
(m)
Aerofli
o
Peso
vazio
(kg)
11.69
AMERICAN
CHAMPION
CITABRIA
7GCAA
Citabria "A"
Package
10.2
AIR TRACTOR
AERO
BOERO
PIPER
CHEROKEE
SOCATA
AT-502B
AB-115
PA-28-140
TB9
15,84
10.78
9.2
9.89
2.22
2.4
2.1
2.25
3.02
7.43
6.9
7.23
7.16
7.72
NACA 4412
503.5
NACA 65415
1949
557
544
647
Peso de
carga
(kg)
Peso
til (kg)
Mxim
o peso
de
decolag
em (kg)
Peso de
Pouso
(kg)
Capacid
ade de
combus
tvel (L)
Velocid
ade
mxima
(km/h)
Velocid
ade de
cruzeiro
Velocid
ade de
stall
(km/h)
1800
Range
(alcanc
e) (km)
Tipo de
trem de
pouso
Distnci
a de
decolag
em (m)
Autono
mia
Capacid
ade de
carga
Motoriz
ao
Hlice
748
244.5
2451
748
4272
770
430
413
975
1060
3629
264 utilizveis
644
57.5
260.7
220.48
230
309
204 a 222
202.8
248
144.84
201
196
88
82.1
67-84
87
108
610
804.7
Usually landed: 85.
UP: 132 (8000ps
3629kg). DOWN:
109 (mesmas
condies)
998
pneus 8.50 x 10-10
lonas
rodas: 5.00 x 5
354
347 (4264kg)
867
4,5 a 5h
15.3
Lycoming IO-540K1J5, 6 cilindros
(Motor a lcool)
1 Lycoming
O-320-A2B
P&W PT6A-15AG
1 AVCO
Lycoming 0
235 C2A,
2400 Horas
TBO
1 Lycoming O320-E2A
Sensenich
M74DM
320 hp, 2.700 rpm
150 hp (111.9
kW)
680 @ 2200 RPM
115 hp a
2800 rpm
150 hp (113
kW)
1.9m de
dimetro
15.14
Hartzell trip
rea de
asa
(m)
Consum
o
Links de
refern
cia
975
15.3
Hartzell HC-B3TN3D/T10282N+4
29,01
15.3
1 Lycoming O320-D2A 4cylinder
horizontallyopposed direct
drive
119 kW (160 hp)
11.9
25L/h
http://www.aeronei
va.com.br/site/conte
nt/produtos/produt
os_ipanema_dados.a
sp
NEIVA
http://en.wiki
pedia.org/wiki
/American_Ch
ampion_Citabr
ia
http://www.airtract
or.com/502a502b
AMERICAN CHAMPION
CITABRIA
AIR
TRACTOR
http://pt.wi
kipedia.org/
wiki/Aero_B
oero_AB115
AERO
BOERO
http://en.wikip
edia.org/wiki/Pi
per_Cherokee#
cite_note-POH16
PIPER
CHEROKEE
http://www.soc
ata.org/
SOCA
TA
Agro
Project
Alongamento AR
6,15
6,80
8,65
7,60
5,59
8,22
6,67
Relao peso
potncia (kg/hp)
Relao de massas
5,63
4,99
6,28
6,70
6,50
6,63
6,41
47,2%
67,3%
45,6%
72,3%
56,8%
61,0%
29,3%
Carga alar (kg/m)
90,33
48,89
147,26
50,33
64,40
89,08
136,67
Para melhor visualizao dos resultados, foram elaborados os seguintes grficos:
Alongamento
9
8
7
6
5
4
3
2
1
0
Alongamento AR
Relao peso potncia (kg/hp)
8,00
6,00
4,00
2,00
0,00
Relao peso potncia (kg/hp)
Relao de massas
80,0%
70,0%
60,0%
50,0%
40,0%
30,0%
20,0%
10,0%
0,0%
Relao de massas
Carga alar (kg/m)
150,00
100,00
50,00
0,00
Carga alar (kg/m)
Fotos das aeronaves pesquisadas:
Figura 1: Ipanema
Figura 2: Desenhos tcnicos do Ipanema
Figura 3: 7GCAA Citabria A
Figura 4: AT-502B
Figura 5: Desenhos tcnicos do AT-502B
Figura 6: PA-28-140
Figura 7: TB9
Com base na anlise histrica e levando em considerao que nosso projeto no seria
to otimizado quanto as aeronaves pesquisadas, fizemos uma estimativa pessimista em
relao ao peso da aeronave:
Tabela 2: Estimativa de massas da aeronave
Tanque de combustvel
Tanque de produtos qumicos
Motor
Fuselagem
Piloto
Peso inicial estimado
150kg
1200kg
100kg
500kg
100kg
2050kg
Para decidir a configurao da aeronave, comparamos as vantagens e desvantagens de
alguns itens:
Tabela 3: Comparao das configuraes possveis
Item
Superfcies de
sustentao
Formas de controle
Subitem
Asa retangular
Vantagens
Barato;
Fcil construo.
Asa baixa
Melhor aplicao do
produto;
Maior
manobrabilidade de
rolamento;
Menor efeito solo;
Exige menor pista de
decolagem.
Menor peso
Configurao
Desvantagens
Baixa eficincia
aerodinmica devida
ao arrasto induzido.
Menor estabilidade
lateral.
Propulso
Sistema de pouso
Estrutura
convencional
Motor a lcool
Motor localizado
na parte frontal
estrutural.
Mais barato;
Ecologicamente
correto.
CG localizado mais na
frente.
Hlice localizada
na frente do
motor
Escoamento sem
perturbaes na
hlice.
Trem
convencional
Mais barato;
Menor necessidade
de manuteno.
Mais fcil construo;
Mais barato.
Resistncia a
corroso;
Relativamente leve;
Facilidade de
manipulao
Resistente;
Baixo peso;
Adequa-se ao formato
desejado.
Alta resistncia.
Resistente;
Duro.
Trem de pouso
fixo
Materiais:
Alumnio
Fibra de vidro
Compsitos
Ao
Exige maior tanque de
combustvel.
Esteira de vrtices na
asa;
Aumento do arrasto
total (aumenta o
arrasto de atrito sobre
a aeronave).
Mais susceptvel a
acidentes no nariz;
Pior visibilidade.
Perda de eficincia
aerodinmica.
Caro.
Denso.
Com base na tabela e nas diretrizes de projeto, a concepo escolhida foi:
Propulso:
Motor a lcool, localizado na parte frontal do avio, atrs da hlice
Sistema de pouso: Trem convencional
Estrutura:
Materiais:
Alumnio;
Fibra de vidro;
Ao (trem de pouso);
Compsitos.
Forma:
Semelhante aos esboos a mo (em anexo).
Nossa anlise aerodinmica partiu da frmula do nmero de Reynolds, que um
nmero adimensional que define vrias propriedades do escoamento de fluidos.
Para a nossa asa, retangular, que ter corda de 1,5m e ter velocidade de cruzeiro de
aproximadamente 200 km/h e considerando a densidade de 1,225 kg/m e viscosidade do ar
17,89
Pa.s (dados retirados do livro ROSA, Edison da. Introduo ao Projeto Aeronutico:
uma contribuio contribuio SAE Aerodesign; colaborao Juliano Toporoski. Florianpolis:
UFSC/GRANTE, 2006):
= 1,225 kg/m
c = 1,5 m
= 17,89 *
Pa.s
Tabela 4: Nmeros de Reynolds
Cl
0,7
0,8
0,9
1,0
1,1
1,2
1,3
1,4
1,5
1,6
1,7
1,8
1,9
2,0
2,1
Velocidade de estol Velocidade
Re para a velocidade atual
(m/s)
atual (m/s)
30,0
3,08E+06
55,9
32,5
3,34E+06
52,3
35,0
3,59E+06
49,3
37,5
3,85E+06
46,8
40,0
4,11E+06
44,6
42,5
4,37E+06
42,7
45,0
4,62E+06
41,0
47,5
4,88E+06
39,5
50,0
5,14E+06
38,2
52,5
5,39E+06
37,0
55,0
5,65E+06
35,9
57,5
5,91E+06
34,9
60,0
6,16E+06
33,9
62,5
6,42E+06
33,1
65,0
6,68E+06
32,3
Para a velocidade estimada de cruzeiro, v = 200/3,6 m/s,
Re = 5,70*106.
Para essa mesma velocidade, tambm calculamos o nmero de Mach usando
vsom=340m/s:
M=0,1634.
Tendo em vista as premissas de projeto definidas no projeto informacional,
escolhemos seis itens para a seleo de perfil. Em seguida, eles foram valorados de acordo
com essas premissas, em ordem de importncia.
Tabela 5: Critrios para seleo de perfil e respectivo peso
Critrio
Arrasto
Fabricao
Coeficiente de momento
Relao Cl/Cd
Resistncia Estrutural
Sustentao
Valor (1 a 5)
2
5
5
1
3
3
Alguns perfis foram escolhidos para passar por essa avaliao, baseados na anlise
histrica e em perfis semelhantes queles. Os candidatos foram:
NACA 0010
NACA 4412
NACA 65-210
SD 7062 (14%)
SD 8040 (10%)
Utilizando o software XFLR5 para anlise aerodinmica, os dados dos perfis foram
avaliados, como mostra a tabela a seguir:
Tabela 6: Comparao dos perfis
Peso
NACA
4412
Nota
Raio
SD
7062
Nota
NACA
0010
Nota
NACA
65-210
Nota
SD
8040
Nota
-
Espessura
(%)
Camber (%)
11,94
14
9,95
9,99
10,01
3.97
0.04
1.1
2.65
Cl mx
1,485
1,936
1,737
1,685
1,899
Cl/Cd mx
177,2
170,0
111,6
109,6
142,0
0,044
0,039
0,028
0,035
0,035
-0,030
-0,039
0,028
-0,003
-0,014
Cd (para Cl
mx)
Cm (para Cl
mx)
Fabricao
Resistncia
Estrutural
Total
64
86
73
37
O perfil selecionado foi o SD 7062, visto que ele apresentou maior pontuao segundo os
critrios estabelecidos.
Para o clculo do coeficiente de sustentao da asa, primeiramente foi feita a
simulao da curva Cl x do perfil, fixando Re=5,7.106 e Mach=0,163. Ento procedemos com
o clculo de CL segundo os mtodos conhecidos.
Segundo da Rosa,
Em que a0 obtido por:
Figura 8: Grfico Cl x alfa do perfil
Ento:
Resolvendo CL por Nicolai,
Aps definidas as equaes, prosseguimos com a simulao da asa no software XFLR5,
com Re=5,6*106, M=0,16, =1,225 kg/m3 e =1,5*10-5m2/s pelo mtodo do 3D pannels.
A tabela abaixo apresenta valores de CL em funo de Cl (Nicolai) e de (linearizao),
sendo os valores de Cl obtidos em simulao no XFLR5. Por comparao, os valores obtidos
com a simulao da asa encontram-se na ltima coluna.
Tabela 7: Valores de CL
Alpha
()
-10,0
Cl
CD
CDp
0,0058
CL por
linearizao
-0,1672
CL por
Nicolai
-0,5152
CL na
simulao
-0,4484
-0,6924
0,0101
-9,5
-0,6355
0,0097
0,0053
-0,1531
-0,4728
-0,4097
-9,0
-0,5782
0,0092
0,0048
-0,1391
-0,4302
-0,3708
-8,5
-0,5204
0,0089
0,0043
-0,1250
-0,3872
-0,3319
-8,0
-0,4627
0,0085
0,0038
-0,1109
-0,3443
-0,2930
-7,5
-0,4047
0,0081
0,0034
-0,0968
-0,3011
-0,2539
-7,0
-0,3462
0,0078
0,0030
-0,0827
-0,2576
-0,2149
-6,5
-0,2875
0,0075
0,0027
-0,0686
-0,2139
-0,1757
-6,0
-0,2286
0,0073
0,0025
-0,0546
-0,1701
-0,1366
-5,5
-0,1697
0,0072
0,0022
-0,0405
-0,1263
-0,0974
-5,0
-0,1108
0,0069
0,0020
-0,0264
-0,0824
-0,0581
-4,5
-0,0517
0,0068
0,0018
-0,0123
-0,0385
-0,0189
-4,0
0,0075
0,0067
0,0016
0,0018
0,0056
0,0204
-3,5
0,0665
0,0066
0,0015
0,0159
0,0495
0,0597
-3,0
0,1256
0,0065
0,0013
0,0299
0,0935
0,0989
-2,5
0,1848
0,0064
0,0012
0,0440
0,1375
0,1382
-2,0
0,2435
0,0063
0,0011
0,0581
0,1812
0,1775
-1,5
0,3021
0,0061
0,0010
0,0722
0,2248
0,2167
-1,0
0,3608
0,0060
0,0010
0,0863
0,2684
0,2559
-0,5
0,4194
0,0058
0,0009
0,1004
0,3120
0,2951
0,0
0,4778
0,0057
0,0010
0,1144
0,3555
0,3342
0,5
0,5362
0,0057
0,0010
0,1285
0,3989
0,3733
1,0
0,5944
0,0056
0,0011
0,1426
0,4423
0,4123
1,5
0,6527
0,0057
0,0012
0,1567
0,4856
0,4513
2,0
0,7109
0,0057
0,0013
0,1708
0,5289
0,4902
2,5
0,7688
0,0058
0,0014
0,1849
0,5720
0,5291
3,0
0,8265
0,0060
0,0016
0,1989
0,6149
0,5678
3,5
0,8837
0,0061
0,0017
0,2130
0,6575
0,6065
4,0
0,9406
0,0063
0,0019
0,2271
0,6998
0,6450
4,5
0,9972
0,0064
0,0021
0,2412
0,7419
0,6835
5,0
1,0532
0,0066
0,0024
0,2553
0,7836
0,7219
5,5
1,1077
0,0067
0,0026
0,2694
0,8242
0,7601
6,0
1,1607
0,0069
0,0029
0,2834
0,8636
0,7982
6,5
1,2152
0,0073
0,0032
0,2975
0,9041
0,8362
7,0
1,2695
0,0077
0,0035
0,3116
0,9445
0,8741
7,5
1,3231
0,0081
0,0038
0,3257
0,9844
0,9118
8,0
1,3757
0,0085
0,0042
0,3398
1,0236
0,9494
8,5
1,4268
0,0090
0,0046
0,3539
1,0616
0,9868
9,0
1,4757
0,0096
0,0051
0,3679
1,0980
1,0241
9,5
1,5230
0,0102
0,0056
0,3820
1,1332
1,0612
10,0
1,5696
0,0107
0,0061
0,3961
1,1678
1,0982
10,5
1,6124
0,0114
0,0067
0,4102
1,1997
1,1349
11,0
1,6524
0,0121
0,0074
0,4243
1,2294
1,1715
11,5
1,6779
0,0129
0,0081
0,4384
1,2484
1,2079
12,0
1,7014
0,0138
0,0090
0,4524
1,2659
1,2441
12,5
1,7200
0,0151
0,0102
0,4665
1,2797
1,2802
13,0
1,7369
0,0166
0,0117
0,4806
1,2923
1,3160
13,5
1,7499
0,0185
0,0136
0,4947
1,3020
1,3516
14,0
1,7587
0,0210
0,0161
0,5088
1,3085
1,3870
14,5
1,7596
0,0243
0,0195
0,5228
1,3092
1,4221
15,0
1,7538
0,0286
0,0238
0,5369
1,3049
1,4571
15,5
1,7345
0,0346
0,0299
0,5510
1,2905
1,4918
16,0
1,7171
0,0413
0,0367
0,5651
1,2776
1,5263
16,5
1,7177
0,0468
0,0423
0,5792
1,2780
1,5605
17,0
1,7147
0,0530
0,0486
0,5933
1,2758
1,5945
17,5
1,7089
0,0599
0,0557
0,6073
1,2715
1,6283
18,0
1,7003
0,0676
0,0635
0,6214
1,2651
1,6618
18,5
1,6869
0,0764
0,0724
0,6355
1,2551
1,6950
19,0
1,6708
0,0860
0,0821
0,6496
1,2431
1,7280
19,5
1,6530
0,0961
0,0924
0,6637
1,2299
1,7607
20,0
1,6337
0,1068
0,1033
0,6778
1,2155
1,7932
20,5
1,6141
0,1178
0,1144
0,6918
1,2009
1,8253
21,0
1,5949
0,1287
0,1255
0,7059
1,1867
1,8572
21,5
1,5767
0,1396
0,1364
0,7200
1,1731
1,8889
22,0
1,5613
0,1499
0,1469
0,7341
1,1617
1,9202
22,5
1,5465
0,1601
0,1572
0,7482
1,1506
1,9512
23,0
1,5316
0,1702
0,1675
0,7623
1,1396
1,9820
23,5
1,5199
0,1797
0,1771
0,7763
1,1309
2,0125
24,0
1,5076
0,1894
0,1869
0,7904
1,1217
2,0426
24,5
1,4965
0,1988
0,1963
0,8045
1,1134
2,0725
25,0
1,4877
0,2077
0,2053
0,8186
1,1069
2,1021
25,5
1,4766
0,2170
0,2148
0,8327
1,0986
2,1313
26,0
1,4689
0,2257
0,2236
0,8468
1,0929
2,1603
26,5
1,4622
0,2342
0,2322
0,8608
1,0879
2,1889
27,0
1,4535
0,2431
0,2412
0,8749
1,0814
2,2172
27,5
1,4477
0,2514
0,2495
0,8890
1,0771
2,2453
28,0
1,4426
0,2595
0,2577
0,9031
1,0733
2,2729
28,5
1,4361
0,2680
0,2663
0,9172
1,0685
2,3003
29,0
1,4321
0,2759
0,2743
0,9313
1,0655
2,3274
29,5
1,4279
0,2838
0,2823
0,9453
1,0624
2,3541
30,0
1,4237
0,2918
0,2903
0,9594
1,0593
2,3805
3,0000
2,5000
2,0000
1,5000
Cl do perfil
1,0000
CL por
linearizao
0,5000
CL por
Nicolai
0,0000
-20,0
-10,0
0,0
10,0
20,0
30,0
-0,5000
-1,0000
Figura 9: Grfico comparativo dos CL
40,0
CL na
simulao
Em seguida, foi plotada a curva de distribuio de CL na asa para =0.
Figura 10: Distribuio de CL na asa
Percebe-se que o coeficiente de sustentao nas pontas de asa menor, diferena
esta causada pelas perdas aerodinmicas. Percebeu-se tambm que a soma das duas reas (da
curva de CL e do downwash) uma constante para cada trecho da asa.
A figura a seguir ilustra o escoamento na asa, onde percebe-se pela corrente de ar a
ocorrncia do efeito dos vrtices de ponta de asa.
Figura 11: Distribuio de Cp na asa e vrtices de ponta de asa
Para o clculo do coeficiente de Oswald, interpolamos na tabela de Glauert, para = 1
(asa retangular) encontra-se = 0,052, e como ew = 1/(1+ ):
ew = 0,95087
O coeficiente de arrasto induzido em funo de CL calculado da forma:
A tabela a seguir mostra os valores de CDi em funo de CL e de :
Tabela 8: CDi em funo de CL e de alfa
Alfa ()
-10,0
-9,5
-9,0
-8,5
-8,0
-7,5
-7,0
-6,5
-6,0
-5,5
-5,0
-4,5
-4,0
-3,5
-3,0
-2,5
-2,0
-1,5
-1,0
-0,5
0,0
0,5
1,0
1,5
2,0
2,5
3,0
3,5
4,0
4,5
5,0
5,5
6,0
6,5
7,0
7,5
8,0
8,5
CL
-0,4484
-0,4097
-0,3708
-0,3319
-0,2930
-0,2539
-0,2149
-0,1757
-0,1366
-0,0974
-0,0581
-0,0189
0,0204
0,0597
0,0989
0,1382
0,1775
0,2167
0,2559
0,2951
0,3342
0,3733
0,4123
0,4513
0,4902
0,5291
0,5678
0,6065
0,6450
0,6835
0,7219
0,7601
0,7982
0,8362
0,8741
0,9118
0,9494
0,9868
CDi
0,0100964
0,0084271
0,0069055
0,0055328
0,0043100
0,0032382
0,0023182
0,0015507
0,0009365
0,0004760
0,0001696
0,0000179
0,0000209
0,0001788
0,0004916
0,0009593
0,0015815
0,0023581
0,0032886
0,0043725
0,0056091
0,0069978
0,0085377
0,0102278
0,0120671
0,0140545
0,0161887
0,0184683
0,0208920
0,0234582
0,0261652
0,0290113
0,0319947
0,0351135
0,0383656
0,0417491
0,0452617
0,0489011
9,0
9,5
10,0
10,5
11,0
11,5
12,0
12,5
13,0
13,5
14,0
14,5
15,0
15,5
16,0
16,5
17,0
17,5
18,0
18,5
19,0
19,5
20,0
20,5
21,0
21,5
22,0
22,5
23,0
23,5
24,0
24,5
25,0
25,5
26,0
26,5
27,0
27,5
28,0
28,5
29,0
29,5
30,0
1,0241
1,0612
1,0982
1,1349
1,1715
1,2079
1,2441
1,2802
1,3160
1,3516
1,3870
1,4221
1,4571
1,4918
1,5263
1,5605
1,5945
1,6283
1,6618
1,6950
1,7280
1,7607
1,7932
1,8253
1,8572
1,8889
1,9202
1,9512
1,9820
2,0125
2,0426
2,0725
2,1021
2,1313
2,1603
2,1889
2,2172
2,2453
2,2729
2,3003
2,3274
2,3541
2,3805
0,0526651
0,0565513
0,0605571
0,0646801
0,0689175
0,0732668
0,0777251
0,0822898
0,0869579
0,0917265
0,0965926
0,1015531
0,1066053
0,1117457
0,1169714
0,1222791
0,1276655
0,1331277
0,1386621
0,1442657
0,1499349
0,1556666
0,1614574
0,1673038
0,1732026
0,1791505
0,1851442
0,1911800
0,1972548
0,2033651
0,2095077
0,2156795
0,2218767
0,2280964
0,2343352
0,2405898
0,2468571
0,2531340
0,2594173
0,2657036
0,2719903
0,2782742
0,2845521
Para o clculo do rendimento na empenagem, necessrio conhecer a geometria da
aeronave.
Vamos ento basear a geometria da aeronave com o intuito de obter t 0,9. O estado
utilizado para calcular o t da aeronave quando em cruzeiro.
Primeiramente, calcula-se o CL da aeronave atravs da frmula:
Quando em cruzeiro, o avio mantm uma altitude constante, ento a sustentao
igual ao peso da aeronave. A velocidade de cruzeiro de 200 km/h = 55,5 m/s. A rea da
superfcie da asa 15 m2.
Possuindo o valor de CL , utilizando a frmula:
Encontra-se = 0,067. Sendo a tg (0,067) = 0,00117.
Com isso, calcula-se Ht da aeronave com a seguinte frmula:
Sendo H =lref . tg e lref = 5,375m. Verifica-se que o valor de t = 0,9 esta compreendido
entre os valores de 0,5 Z 0,9. Sendo Z = 0,9 , calcula-se o Ht da aeronave, que Ht = 0,98m.
O prximo passo foi o clculo das resistncias da agro aeronave.
Utilizando a tabela da pgina 132 do livro, estimam-se valores de CD:
Parte
Fuselagem
Nacele
Flaps
Bequilha
Trem de pouso
Descrio
Forma aerodinmica,
sem salincias
Acima da asa em
avies pequenos
60% da envergadura,
defletido 30
Roda da bequilha e
sua estrutura
Rodas e estrutura
expostas
CD
0,06
A
SC
0,1
SC
0,025
S1
0,25
Spr
0,4
Spr
Figura 12: CL x CD da aeronave (verde) e da asa (vermelho)
Em seguida, consideramos o efeito solo. Este efeito de especial importncia na nossa
aeronave, uma vez que a mesma voa em baixas altitudes para aplicao do veneno. O grfico
abaixo mostra a influncia do efeito solo a uma altitude de 5m.
Figura 13: Grfico CL x CD da aeronave sem levar em conta efeito solo (verde) e levando em conta o
efeito solo (azul)
Estipulamos os seguintes dados para as rodas:
Raio
0,10 m
0,075 m
Rodas dianteiras
Roda traseira
Dimetro do pneu
0,05 m
0,04 m
Calculamos o arrasto na aeronave utilizando a seguinte tabela:
Tabela 9: Clculo do arrasto para diferentes fases do vo
Corrida (=0)
Cd0 para =0
0,006
Decolagem (=20)
Arrasto de perfil (Cd0)
Cd0 para =20
(aproximadamente de
estol)
0,012
Arrasto induzido (CDi)
0,0051
Cruzeiro (=2)
Cd0 para =2
0,006
0,2703
Arrasto parasita (CDP)
0,0111
0,026918
Arrasto de rolagem (Q)
Q=0
Q=0
16892,8
0
Arrasto Total
8606,2 N
1871,2 N
362,9 N
83,69 m/s
100,43 m/s
117,17 m/s
Por semelhana com o tipo de aeronave, o motor selecionado foi um motor genrico
de uma aeronave de transporte ttico areo 350HP, cuja curva foi fornecida pelo professor:
10000,00
TD [N]
8000,00
TD
[N]
6000,00
4000,00
2000,00
0,00
0
50
100
150
Figura 14: Curva de potncia disponvel do motor escolhido
200
Tabela 10: Potncia disponvel e requerida
v [m/s]
10
20
30
40
50
60
70
80
90
100
110
120
130
140
150
160
170
180
CL
CD
22,345
5,58625
2,482778
1,396563
0,8938
0,620694
0,45602
0,349141
0,275864
0,22345
0,184669
0,155174
0,132219
0,114005
0,099311
0,087285
0,077318
0,068966
TR [N]
0,038906
0,022189
4504,158694
2249,249594
0,012359
0,010689
0,009951
0,009457
0,009131
0,009437
0,008838
831,7877035
586,1839221
450,6862428
367,6899719
314,8261856
291,7391438
257,8883437
0,008713
230,6402471
TD [N]
7835,50
7652,00
7449,50
7228,00
6987,50
6728,00
6449,50
6152,00
5835,50
5500,00
5145,50
4772,00
4379,50
3968,00
3537,50
3088,00
2619,50
2132,00
PD
[N.m/s]
78355
153040
223485
289120
349375
403680
451465
492160
525195
550000
566005
572640
569335
555520
530625
494080
445315
383760
PR [N.m/s]
225207,9347
134954,9756
66543,01628
52756,55299
45068,62428
40445,89691
37779,14227
37926,0887
36104,36812
36902,43953
700000
600000
500000
400000
Series1
300000
200000
100000
0
0
O volume de cauda
definido da seguinte forma:
50
100
150
200
um parmetro, fundamental para a estabilidade longitudinal,
A partir de um valor referncia obtido da Tabela 10.1 pg. 187 para a aeronave Piper
Cherokee, que possui S = 16,26m e b = 10,16m, valores semelhantes aos estimados para a
nossa aeronave. Analisando a tabela obtivemos um valor para o volume de cauda
Estimou-se ainda um valor para = 5m.
Com isso determinamos o valor de
Analogamente ao valor de
= 0,61.
obtido a seguir, onde o valor de
foi obtido da tabela
10.1 e
Com os valores das reas dos estabilizadores, pode-se determinar o dimensionamento
final dos mesmos, conforme se pode observar na tabela:
Profundor
2,745
3
0,915
-
S (m)
b (m)
c (m)
c (m)
Leme
1,11
1,5
1
0,48
Para fazer o clculo de hn, utiliza-se a frmula:
Sendo a = 4.437, aII = 1.8, h0 = 0.375, = 0.067, VH = 0.61 e t = 0.9.
Com isso, calcula-se hn = 0,582795.
Para fazer o clculo de MS, utiliza-se:
Sendo hn = 0,0582795 e h = 0.51, chega-se a:
MS = 0,072795
A curva CM x foi plotada no software XFLR5, tendo como base o perfil selecionado e o
comprimento da corda e a envergadura da asa.
Figura 15: Curva CM x alfa da asa
Sendo a massa da aeronave 2050 kg, seu peso 20110 N. Com este peso, calcula-se n1
=2,179.
Para a velocidade de cruzeiro, q = 3788,7 Pa. Sendo S = 15 m, encontra-se a
sustentao para o maior CL da aeronave, que L = 40217 N. Pensando nisso, no diagrama CM
x obtivemos o valor de ngulo de ataque para essa carga, 14,5.
Em seguida, plotamos a curva da distribuio de momento fletor ao longo da asa para
este ngulo, com a finalidade de calcular os esforos a que ela est submetida. Escolhemos o
formato da viga, uma viga caixo, para melhor suportar a toro.
Figura 16: Distribuio de Momento Fletor na asa para velocidade de cruzeiro a =14,5
Tendo em vista a segurana, arredondamos o momento fletor mximo da longarina
para 5.104 N.m. A partir deste valor, e levando em considerao o camber do perfil
selecionado, iniciamos o dimensionamento da longarina.
Figura 17: Esquematizao do perfil
O dimensionamento da longarina se dar levando em conta a tenso das partes
superior e inferior da viga, visto que as partes laterais no sero to afetadas pelo momento
fletor.
Portanto, a altura mxima que a longarina pode ter de 0,21m. Tomando uma altura
menor para que a aba caiba no perfil, escolhemos uma altura de 0,10m e uma espessura para
a alma e para a aba de 0,02m. O prximo passo ser encontrar a largura da viga. Relacionando
tenso e o momento fletor, temos:
Para a viga caixo:
O material escolhido foi o Alumnio 2014-T6, cujas propriedades esto listadas a seguir:
Tabela 11: Propriedades do Alumnio 2014-T6
Densidade (kg/m3)
2800
Mdulo de Elasticidade (GPa)
Tenso de Escoamento (MPa)
Tenso de Cisalhamento (MPa)
72,4
662
290
Utilizamos um fator de segurana 1,2, de forma que aumentasse a segurana mas no
fosse muito acrscimo de peso.
Ao compararmos esta frmula com as anteriores, obtemos:
b=0,0583 m
O peso da longarina foi calculado pela densidade do material:
m=132 kg
Figura 18: Esquematizao da longarina
Concluso
Ao final desenvolvimento do projeto, percebeu-se que o avio ainda no encontra-se
pronto para comercializao. Para atingir esse nvel, um estudo mais aprofundado deveria ser
feito, complementado com mais simulaes, testes e otimizaes.
Esse curso, foi, entretanto, muito instrutivo na aprendizagem de aeronutica. Apesar
de algumas dificuldades, como obteno de dados de aeronaves existentes, utilizao dos
programas de simulao e falta de literatura aprofundada no assunto em portugus, muito foi
aprendido, no s em sala de aula, mas principalmente na elaborao do projeto de um avio.
Referncias
ROSA, Edison da; TOPOROSKI, Juliano. Introduo ao projeto aeronutico: uma
contribuio competio SAE Aerodesign. Florianpolis, SC: UFSC, Centro Tecnolgico, 2006.
283p.
http://www.airtractor.com/ - Acesso em setembro de 2010.
http://www.globoaviacaoagricola.com.br/ - Acesso em setembro de 2010.
http://www.aeroneiva.com.br/ - Acesso em setembro de 2010.
http://www.engbrasil.eng.br/ - Acesso em setembro de 2010.
http://www.ganum.unifei.edu.br/ - Acesso em outubro de 2010.
http://www.inf.unioeste.br/~rogerio/02f-Aplicacao-Aerodinamica2.pdf - Acesso em
novembro de 2010.
http://www.matweb.com acesso em 08 de dezembro de 2010.