Добавил:
timofeev.9@mail.ru Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

1503_Конструкция и расчет самолета на прочность крыло

.pdf
Скачиваний:
137
Добавлен:
05.01.2020
Размер:
293.91 Кб
Скачать

сужения, длины центроплана и т.д.). Можно воспользоваться данными приведенными в таблице 2.

Таблица 2

Распределение циркуляции посечениямдля трапецевидныхкрыльев

 

 

 

~

пл . при 2 ≤

л ≤ 5

 

 

z

 

 

Г

 

 

 

 

 

 

 

 

 

η=2

 

η=3

 

η=4

η=5

 

 

 

 

 

0

1,2604

 

1,321

 

1,3553

1,3792

0,1

1,254

 

1,3124

 

1,3454

1,3669

0,2

1,2242

 

1,2858

 

1,3152

1,3413

0,3

1,1989

 

1,2395

 

1,2625

1,2794

0,4

1,1463

 

1,1713

 

1,1857

1,198

0,5

1,0763

 

1,0811

 

1,0845

1,0884

0,6

0,9911

 

0,9727

 

0,9631

0,9577

0,7

0,8946

 

0,8622

 

0,8287

0,8137

0,8

0,7865

 

0,7241

 

0,6875

0,6624

0,9

0,6345

 

0,5664

 

0,5236

0,4669

0,95

0,4933

 

0,447

 

0,3979

0,3817

1,0

0,0

 

0,0

 

0,0

0,0

Расчетная погонная аэродинамическая нагрузка (направление qаэр. приближенно можно считать перпендикулярным плоскости хорд крыла) для плоского крыла при

0о ≤ χ0,25 ≤ 10о

 

 

qаэр. =

Pразр.

~

 

 

× Гпл.

(1)

L

 

 

 

Для крыльев со стреловидностью 10о ≤ χ0,25 ≤ 35о

где

~ч = [~пл.

Г

Г

qаэр. =

Pразр.

~

 

 

 

 

 

× Гч

, (2)

 

 

 

L

 

 

 

 

 

 

 

 

- 2 × (1 -

 

)× (1 - cosч)]×

1

(3)

z

cosч

 

 

 

 

 

 

 

При учете стреловидности не принимается во внимание крутка крыла. Для крыльев со стреловидностью χ › 35 о формула (3) дает ошибку в значениях циркуляции до 20 %.

Методика расчета для неплоских крыльев любой формы изложена в рабо-

те [1].

По эпюре распределенных нагрузок qаэр, вычисленных формулам (1) или (2), строятся последовательно эпюры Qаэр. известные дифференциальные зависимости, находим

для 12 сечений по и Mаэр.. Используя

Qаэр. =

z

аэр.

q

dz

L /

2

 

M аэр. =

 

z

аэр.

 

Q

dz

L / 2

Интегрирование проводится численно, используя метод трапеций (рис.2). По результатам вычислений строятся эпюры изгибающих моментов и перерезывающих сил.

11

qаэр

qi

qi+1

z

O

zi 0,5

 

0,0

1,0

Рис.2

4.2Определение массовых и инерционных сил

4.2.1Определение распределенных сил от собственного веса конструкции крыла. Распределение массовых сил по размаху крыла с незначительной погрешностью можно считать пропорциональным аэродинамической нагрузке

qкр. = Gкр. × qаэр. ,

G

или пропорционально хордам

 

G

× n

э

× f

 

qкр. =

кр.

 

y

 

× b

 

S

 

 

 

 

 

 

 

Погонная массовая нагрузка приложена по линии центров тяжести сечений, расположенной, обычно, на 40-50% хорды от носка. По аналогии с аэродинамическими силами определяются Qкр. и Mкр.. По результатам вычислений строят эпюры.

4.2.2. Определение распределенных массовых сил от веса баков с топливом. Распределенная погонная массовая нагрузка от баков с топливом

q топл. = c × b × B × g × nэy × f ,

где γ – удельный вес топлива; B – расстояние между лонжеронами, являющимися стенками бака (рис.3).

Относительная толщина профиля в сечении

c = cкорн. - (cкорн. - cконц. )× z (4)

Центр тяжести условно можно считать расположенным на середине между передним и задним лонжеронами. В целях упрощения расчетов пренебре-

12

гаем кривизной крыла, то есть форму баков принимаем в виде усеченных пирамид (рис. 4).

Далее находятся Qтопл. и Mтопл. и строятся их эпюры. При вычислении Qтопл. и Mтопл. следует вводить дополнительные сечения границ расположения

Sбака1

Sбака2

B – текущие хорды баков

Передний лонжерон

Задний лонжерон

l бака

z = 0,8

O z

0.0

1.0

Рис. 3

топлива, если они не совпадают с координатами таблицы 1.

4.2.3. Построение эпюр от сосредоточенных сил. Сосредоточенные инерционные силы от агрегатов и грузов, расположенных в крыле и присоединенных к крылу, приложены в их центрах тяжести и принимаются направленными параллельно аэродинамическим силам. Расчетная сосредоточенная нагрузка

Pсоср. = Gсоср. × nэy × f

H1'

H2'

l бака

H1

H2

Рис.4

Результаты приводятся в виде эпюр Qсоср. и Mсоср.. Строятся суммарные эпюры QΣ и MxΣ от всех сил, приложенных к крылу, с учетом их знаков:

QΣ = Qаэр. + Qкр. + Qтопл. + Qсоср.

M xΣ = M аэр. + M кр. + M топл. + M соср.

13

4.3 Вычисление моментов, действующих относително условной оси

4.3.1. Определение M

аэр. от аэродинамических сил. Аэродинамические силы

 

z усл.

действуют по линии центров давления, положение которой считается извест-

ным. Вычертив крыло в плане, отметим положение ∆Qаэрi на линии центров давления и по чертежу определим hаэрi (рис.5).

Линия четвертей хорд

 

DQi

 

 

аэр.

Линия центров давления

 

DQiкр.

 

Расчетное сечение

 

 

 

 

 

zусл.

 

 

 

 

 

 

 

 

 

 

 

 

 

hаэр.

 

hкр.

 

 

 

 

 

 

 

 

 

 

 

 

i

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Линия центров тяжести крыла

O z

0.0

1.0

Рис. 5

Далее вычисляем M Zаэрycл.i иM Zаэрycл. по формулам

DM

аэр

= DQаэр. × hаэр

 

zусл. i

i

i

n

M Zycаэрл. = M Zycаэрл.i i =1

и строим эпюруM Zаэрycл. .

4.3.2. Определение M Zycл. от распределенных массовых сил крыла

( M Zycкр. л. и M Zycтоплл.. ). Массовые силы, распределенные по размаху крыла, действу-

ют по линии центров тяжести его конструкции (см. рис. 5).

DM

кр.

= DQкр. × hкр.

,

 

Zycл.

i

i

 

где DQiкр. - расчетная сосредоточенная сила от веса части крыла между двумя соседними сечениями; hiкр. - плечо от точки приложения силы DQiкр. до оси zусл. .

14

Аналогично вычисляются значения DM Zтоплycл.. . По расчетам строятся эпюры

M

кр.

 

и M

топл. .

 

 

 

 

 

 

Zycл.

 

Zycл.

 

 

 

 

 

 

 

 

4.3.3. Определение M Zycл. от сосредоточенных сил.

 

 

 

 

 

 

 

M соср.

= Pсоср. × hсоср. ,

 

 

 

 

 

 

 

Zycл.i

i

i

 

где

Pсоср.

= P × n

э

× f , расчетный

вес

каждого агрегата или груза; hсоср.

-

 

 

 

i

i

y

 

 

i

 

расстояние от центра тяжести каждого агрегата или груза до оси.

После вычисления M

соср.

определяется суммарный момент M ZΣycл. от

 

Zycл.i

 

всех сил, действующих на крыло, и строится эпюра M ZΣycл. (имеется ввиду ал-

гебраическая сумма).

Σ

аэр.

кр.

топл.

соср.

M Zycл.

= M Zycл.

+ M Zycлc

+ M Zycлc

+ M Zycл.

4.4. Определение расчетных значений M изг. и M кр. для заданного сечения крыла

Для определения M изг. и M кр. следует:

- найти приближенное положение центра жесткости (рис. 6)

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

Hi2 ×li

 

 

 

x

 

 

=

 

i=1

 

,

 

 

 

m

 

 

ц.ж.

 

 

 

 

 

 

 

 

 

 

 

Hi2

 

 

 

 

 

 

 

 

 

i=1

 

 

где H i - высота i-го лонжерона;

l i - расстояние от выбранного полюса А до

стенки i-го лонжерона; m – количество лонжеронов;

- вычислить момент относительно оси Z, проходящей через приближен-

ное положение центра жесткости и параллельной оси Zусл.

 

M ZΣ = QΣ × a - M ZΣycл. ;

 

a

 

 

 

 

 

 

Xц.ж.

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H1

 

 

 

 

 

 

А

 

 

 

 

 

 

 

H2

 

Zусл.

 

 

 

 

 

l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис.6

- для стреловидного крыла сделать поправку на стреловидность (рис.7) по формулам

15

 

M кр.

= M zΣ

= M zΣ ×cosc + M xΣ × sinc

 

 

 

1

 

 

 

 

 

 

M изг.

= M xΣ

= M zΣ × sinc - M xΣ ×cosc

 

 

 

1

 

 

 

 

 

 

Линия четвертей хорд

 

 

 

 

 

 

 

x

x1

Положение расчетного

 

 

 

Σ

 

сечения после поправки

 

 

 

 

 

 

 

 

 

M x

 

на стреловидность

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

χ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

90o

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M Σz χ

 

 

 

 

 

 

 

z1

 

 

 

 

 

 

Расчетное сечение

 

 

 

 

Ось жесткости

 

 

 

z

 

 

 

 

 

O

 

 

 

 

1.0

 

 

 

 

 

 

 

 

0.5

 

 

0.0

 

 

 

Рис. 7

5. Выбор конструктивно-силовой схемы крыла, подбор параметров расчетного сечения

5.1. Выбор конструктивносиловой схемы крыла

Тип конструктивно-силовой схемы крыла выбирается с использованием рекомендаций, изложенных в лекциях и работах [1,2,3,7].

5.2. Выбор профиля расчетного сечения крыла

Относительная толщина профиля расчетного сечения определяется по формуле (4). Из работы [9] выбирается симметричный (для простоты) профиль, соответствующий по толщине c рассматриваемому типу самолета и составляется таблица 3. Подобранный профиль вычерчивается на миллиметровой бумаге в масштабе (1:10, 1:25). В случае отсутствия в справочнике профиля необходимой толщины можно взять из справочника наиболее близкий по толщине профиль и все данные пересчитать по формуле

Таблица 3.

16

x%

yв%

yн%

(см)

(см)

(см)

(см)

 

yв

yн

h=yв

+yн

0

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

100

 

 

 

 

 

 

y = yсправ.

×

 

 

 

с

 

,

 

 

 

 

 

 

 

ссправ.

где y – расчетное значение ординаты; yсправ. - табличное значение ординаты; cсправ. - таб-

личное значение относительной толщины профиля крыла.

Для стреловидного крыла следует сделать поправку на стреловидность по формулам

b= b ×cos c,

c ¢ = c cos c

5.3.Подбор параметров сечения (ориентировочный расчет)

5.3.1.Определение нормальных усилий, действующих на панели крыла.

Для последующих расчетов будем считать положительными направления Mизг. , и QΣ в расчетном сечении (рис. 8). Пояса лонжеронов и стрингеры с присоединенной обшивкой воспринимают изгибающий момент Mизг. . Усилия,

нагружающие панели, можно определить из выражения

N = ± M изг. ,

H ср.

где Hср. = F ; F – площадь поперечного сечения крыла, ограниченная крайними

B

лонжеронами; B - расстояние между крайними лонжеронами; (рис. 9).

Для растянутой панели усилие N принять со знаком плюс, для сжатой - со знаком минус.

17

На основе статистических данных в расчете следует принять усилия, воспринимаемые полками лонжеронов - N лна = a × N , Nстр. = b × N , Nобш. = g × N .

Значения коэффициентов a, b, g даны в таблице 4 и зависят от типа кры-

ла. Таблица 4.

Типкрыла

Значения коэффициента

 

α

β

γ

Кессонное

0,70

0,20

0,1

моноблочное

0,30

0,45

0,25

5.3.2. Определение толщины обшивки. Толщину обшивки d для растянутой зоны определяют по 4-ой теории прочности:

 

 

 

1

 

 

M изг.

× g 2

 

M кр.

 

2

 

d

раст.

=

 

×

 

 

 

 

+ 3 ×

 

 

 

,

 

 

 

 

 

 

 

sв

 

 

 

 

 

 

2 × Hср.

 

 

 

 

 

 

 

 

Hср. × B

 

× B

 

где sв - напряжение предела прочности материала обшивки; g - коэффициент, значение которого приведено в таблице 4. Для сжатой зоны толщину обшивки следует принять равной dсж. =1,3 ×d раст. .

5.3.3.Определение шага стрингеров и нервюр. Шаг стрингеров и нервюр а выбирают с таким расчетом, чтобы поверхность крыла не имела недопустимой волнистости.

Для расчета прогибов обшивки считаем ее свободно опертой на стринге-

y

p

 

F

c

 

 

 

 

 

 

B

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hср.

 

 

 

F

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

O

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 9

 

 

Рис. 10

 

 

 

 

 

 

 

 

ры и нервюры (рис. 10). Наибольшее значение прогиба достигается в центре рассматриваемой пластины:

18

Nстр.
Nстр.

 

w = d ×

 

p × ℓ

4

,

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

где p = P

-удельная нагрузка на крыло;

D =

 

 

 

E × d3

 

 

 

 

 

-цилиндрическая же-

 

 

 

 

S

 

12

× (1 - n2 )

 

сткость обшивки. Значения коэффициентов d в зависимости от a

приведены в

 

 

 

 

 

 

 

 

 

работе [8]. Обычно это отношение равно 3.

Расстояние между стрингерами и нервюрами следует выбирать так, что-

бы w £ 0,02 ¸ 0,002 .

Число стрингеров в сжатой панели

nсж. = Bo ,

где Bo - длина дуги обшивки сжатой панели.

Количество стрингеров в растянутой панели следует уменьшить на 20%. Как отмечалось выше, расстояние между нервюрами a = 3 × ℓ.

5.3.4. Определение площади сечения стрингеров. Площадь сечения стрингера в сжатой зоне в первом приближении

Fстр.сж. = sкр.стр. × nсж. ,

где σкр.стр. - критическое напряжение стрингеров в сжатой зоне ( в первом при-

ближении σкр.стр. = 0,5 ÷ 0,6 σв ).

Площадь сечения стрингеров в растянутой зоне

Fстр. раст. = sв × 0,9n раст. ,

где σв - предел прочности материала стрингера при растяжении.

5.3.5. Определение площади сечения лонжеронов. Площадь полок лонжеронов в сжатой зоне

F

=

N лн

,

 

л. сж.

 

sкр. л-н

 

 

где σкр. л-н - критическое напряжение при потере устойчивости полки лонжеро-

на. σкр. л-н ≈ 0,8σв (берется предел прочности материала лонжерона).

Площадь каждой полки двухлонжеронного крыла находится из условий

 

 

 

 

Fлн1

=

H12

 

, F

 

= F

+ F

, (5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fл- н 2 H 22

 

л- н сж.

 

лн1

 

лн 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а для трехлонжеронного крыла

 

 

 

 

 

 

 

 

 

 

 

 

Fлн1

=

H12

,

Fлн 2

=

 

H 22

, F

 

= F

 

+ F

н 2

+F

лн 3

(6)

 

 

 

 

 

 

 

 

 

Fл- н 2 H 22

 

Fл- н 3

 

H32

л- н сж.

лн1

л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

19

 

 

 

 

 

 

 

Площадь лонжеронов в растянутой зоне

F

=

N лн

,

 

лн раст.

 

sв × k

 

 

где k – коэффициент, учитывающий ослабление поясов лонжеронов крепежными отверстиями; при заклепочном соединении k = 0,9 ÷ 0,95.

Площадь каждой полки находится аналогично площади в сжатой зоне из условий (5) или (6).

5.3.6. Определение толщины стенок лонжеронов. Предполагаем, что вся перерезывающая сила воспринимается стенками лонжеронов

n

QΣ = Qi ,

i =1

где Qi - сила, воспринимаемая стенкой i-го лонжерона. Для трехлонжеронного крыла (n=3)

Q1 = H12 , Q2 = H 22 ,

Q2 H 22 Q3 H32

где H1, H 2 , H3 - высоты стенок лонжеронов в расчетном сечении крыла.

Толщина стенки

dст. i

=

 

Qi

.

(7)

Hi

× tкр. i

 

 

 

 

Здесь τкр. i - критическое напряжение потери ус-

тойчивости стенки лонжерона крыла от сдвига (рис. 11). Для вычислений τкр. i следует принять все четыре стороны стенки свободно опертыми:

 

 

 

2

 

 

 

dст.i

 

 

tкр. i = kτi

× 0,9 × E ×

 

 

,

(8)

Hi

 

 

 

 

 

 

 

H

i

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

где kτi = 5,60 + 3,78 ×

 

 

 

 

при

a > Hi , при a < H i следует заменить в (8) Hi на a,

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а в формуле для k

- (Hi

 

 

) на (a

)

. Формула (8) справедлива для τ

кр

 

≤ τ

пц.

 

τi

 

 

 

a

 

 

Hi

 

 

 

 

 

i

 

Подставляя значения τкр.i

из (8) в (7), находим толщину стенки i-го лон-

жерона

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dст.i

= 3

 

 

Qi × Hi

.

 

 

 

 

 

 

 

 

 

 

 

 

 

0,9 × kτi × E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6. Расчет сечения крыла на изгиб

20