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

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

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

Для расчета сечения крыла на изгиб вычерчивается профиль расчетного сечения крыла, на котором размещаются пронумерованные стрингеры и лонжероны (рис.12). В носике и хвостике профиля следует располагать стрингеры с большим шагом, чем между лонжеронами. Расчет сечения крыла на изгиб проводится методом редукционных коэффициентов и последовательных приближений.

y

 

 

y 1

 

xi

Fi

 

 

 

3

i

1

2

x1

 

yi O

O

 

x n

Рис. 12

6.1. Порядок расчета первого приближения

Определяются в первом приближении приведенные площади поперечного сечения продольных ребер (стрингеров, поясов лонжеронов) с присоединенной обшивкой

Fiпр. = Fi Σ × j1i , Fi Σ = Fi + Fiобшприс. . , (9)

где F - действительная площадь сечения i-го ребра; F обш.

-

присоединенная

i

i прис.

 

 

 

площадь обшивки ( F обш

= 0,8 × d × ℓ- для растянутой панели,

F обш

= 30 × d2 -

i прис.

 

 

i прис.

 

для сжатой панели); j1i - редукционный коэффициент первого приближения. Если материал полок лонжеронов и стрингеров разный, то следует сде-

лать приведение к одному материалу через редукционный коэффициент по модулю упругости

ji E = E i ,

Eo

где Ei - модуль материала i-го элемента; Eo - модуль материала, к которому приводится конструкция (как правило, это материал пояса самого нагруженного лонжерона). Тогда

j1i Σ = j1i × ji E

21

Вслучае разных материалов поясов лонжеронов и стрингеров в формулу

(9)вместо j1i подставляется j1i Σ .

Определяем координаты xi и yi центров тяжести сечений продольных элементов профиля относительно произвольно выбранных осей x и y (рис. 12)

и вычисляем статические моменты элементов Fiпр. × xi и Fiпр. × yi .

Определяем координаты центра тяжести сечения первого приближения по формулам

 

 

n

 

 

 

Fiпр. × xi

 

x1

=

i =1

,

n

c

 

 

 

 

Fiпр.

 

 

 

i =1

 

 

 

n

 

 

 

Fiпр. × yi

 

y1

=

i =1

.

n

c

 

 

 

 

Fiпр.

 

 

 

i =1

 

Через найденный центр тяжести проводим оси x1и y1 (ось x1 удобно выбрать параллельной хорде сечения) и определяем координаты центров тяжести всех элементов сечения относительно новых осей.

Вычисляем моменты инерции (осевые и центробежный) приведенного сечения относительно осей x1 и y1 :

n

 

n

I пр1 .

n

× yi1 .

I пр1 . = Fiпр. × yi12 ,

I пр1 . = Fiпр. × xi12 ,

1 = Fiпр. × xi1

x

y

i =1

x y

i =1

 

i =1

 

 

 

Определяем угол поворота главных центральных осей сечения:

 

2 × I

пр1 .

1

 

tg 2a =

 

x y

 

.

I пр1 . - I пр1 .

 

y

 

x

Если угол α будет больше 5о, то оси x1и y1 следует повернуть на этот угол (положительное значение угла соответствует вращению осей по часовой стрелке) и далее вести расчет относительно главных центральных осей. В целях упрощения расчета угол α рекомендуется вычислять только при расчетах последнего приближения. Обычно, если ось x1выбрана параллельно хорде сечения , угол α оказывается незначительным и им можно пренебречь.

Определяем напряжения в элементах сечения в первом приближении

si

=

M

изг.

× yi

× ji × jE .

1

 

1

1

 

 

I

пр1 .

 

 

 

 

 

x

 

 

Полученные напряжения s1i сравниваем с σкр.л−н и σкр.стр. для сжатой панели и с σв.л−н и σв.стр. - для растянутой панели.

22

6.2. Определение критических напряжений стрингеров

Критическое напряжение стрингера σкр.стр. вычисляется из условия об-

щей и местной форм потери устойчивости. Для вычисления σкр.стр. общей фор-

мы потери устойчивости используем выражение

 

 

 

s

 

= s

 

×

 

1 + ν

,

(10)

 

кр.стр.

 

 

 

 

 

ув

 

 

 

в

1 + n + n2

 

 

гден=

. Здесь σкр. э - критическое напряжение, вычисленное по формуле

 

 

укр. э

 

 

 

 

 

 

 

 

 

 

Эйлера:

 

 

 

 

 

 

 

 

 

 

 

 

 

sкр. э

=

p2

× E

, l =

n × a

,

(11)

 

 

 

l2

 

 

ix

 

 

 

 

 

 

 

 

 

 

где μ - коэффициент, зависящий от условий опирания концов стрингера; a - шаг нервюр; λ - гибкость стрингера с присоединенной обшивкой; ix - радиус инер-

ции относительно центральной оси сечения.

В формуле (11) под ix следует понимать imin , но в целях упрощения по-

ложение главной инерциальной оси считаем совпадающим с осью x. В свою очередь

i x

=

I

x

,

F пр.

 

 

 

 

 

 

стр.

 

где I x - момент инерции стрингера с присоединенной обшивкой относительно оси x (рис.13); Fстрпр.. - площадь сечения стрингера с присоединенной обшивкой. Ширина присоединенной обшивки берется рав-

ной 30 δ (рис.13).

F пр.

= F

+ F прис.

стр.

стр.

обш.

При этом

I

x

= I

x

+ F

прис. × a2

+ I

x

+ F

× a

2

,

 

 

обш. 1

 

стр.

 

2

 

 

 

 

 

1

 

 

 

2

 

 

 

23

siл− на

где I x1 - момент инерции присоединенной обшивки относительно собственной центральной оси x1 (обычно значения I x1 -малы); I x2 - момент инерции стрин-

гера относительно собственной центральной оси x2 .

Для вычисленияσкр.стр. местной формы потери устойчивости рассмотрим

потерю устойчивости свободной полки стрингера как пластины, шарнирно опертой по трем сторонам (рис.14). На рис. 14 обозначено: а – шаг нервюр; b1 – высота свободной полки стрингера (рис.13). Для рассматриваемой пластинки σкр. вычисляется по асимптотической формуле (10), в которой

sкр. э

 

δ

 

2

= 0,9 × kσ × E

 

с

 

,

 

 

b1

 

 

где kσ – коэффициент, зависящий от условий нагружения и опирания пластины, dс – толщина свободной полки стрингера.

Для рассматриваемого случая

b

 

2

kσ = 0,425 +

1

.

a

 

 

 

Для сравнения с действительными напряжениями, полученными в результате редуцирования, выбирается меньшее напряжение, найденное из расчетов общей и местной потери устойчивости.

В процессе редуцирования необходимо обратить внимание на следующее: если напряжения в сжатой полке лонжерона окажутся больше или равными разрушающим в любом из приближений, то конструкция крыла не способна выдержать расчетную нагрузку и ее надо усилить. Дальнейшие приближения в этом случае делать не следует. Если в каком-либо сжатом стрингере с номером

"k" (с присоединенной обшивкой) напряжение sik окажется меньше σкр.стр. k , то редукционный коэффициент для него и в последующем приближении следует оставить прежним jik+1 = jik ; если в каком-либо сжатом стрингере (с присоединенной обшивкой) с номером "m" напряжение s im окажется больше σ кр.стр. m то в последующем

приближении редукционный коэффициент следует вычислять по формуле

ji +1 = sкр.стр. m × yлi − на m × ymi ;

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

высит σкр.стр. j , то конструкция явно перетяже-

лена и требует облегчения.

24

В растянутой зоне уточнение редукционных коэффициентов в процессе последовательных приближений ведется так же, но сравнение расчетных напряжений ведется не с σкр.стр. , а с σв .

В результате мы получаем новые уточненные редукционные коэффициенты последующего приближения jij+1 . Далее рассчитываем следующее приб-

лижение в том же порядке и снова уточняем редукционные коэффициенты. Расчет продолжается до тех пор, пока редукционные коэффициенты двух последующих приближений практически совпадут (в пределах 5%).

7. Расчет сечения крыла на сдвиг.

Расчет сечения крыла на сдвиг ведется без учета влияния кручения (поперечная сила QΣ считается приложенной в центре жесткости сечения, полагая, что на сдвиг работают стенки лонжеронов и обшивка).

7.1. Порядок расчета

Для расчета многоконтурного сечения на сдвиг делаются продольные разрезы в панелях таким образом, чтобы контур стал открытым. Для сечения крыла разрезы удобно делать в плоскости хорд в носке крыла и в стенках лонжеронов (рис. 15). В местах разрезов прикладываются неизвестные замыкающие погонные касательные усилия qI , qII , qIII .

Погонные касательные усилия qi в обшивке панелей сечения крыла определяются как сумма погонных касательных усилий qoi в незамкнутом контуре и замыкающих усилий qI , qII , qIII . Усилия qoi определяются формулой

 

 

Q

× S

отс.

 

qoi

= -

Σ

 

1, i 1

,

(12)

I

 

 

 

 

xред.

 

гдеQΣ -расчетная перерезываю-

 

k =i−1

 

 

 

щая сила; S1отс. =

y

k

× F пр.

-

,i−1

k =1

k ред.

 

 

 

 

 

статический момент площади части сечения, ограниченного 1-м и (i-1) – м ребрами (принятый порядок нумерации ребер очевиден из

рис. 14); I xред. - главный момент инерции всего сечения, причем положение центра тяжести берется

25

из последнего приближения расчета на изгиб.

В формуле (12) направление поперечной силы QΣ считается положительным при его совпадении с положительным направлением оси y, т.е. вверх. Положительные направления потоков касательных усилий совпадают с направлением обхода начала координат по часовой стрелке.

Для определения замыкающих потоков погонных касательных усилийqI , qII , qIII составляем канонические уравнения

[A]×{q}+ {Ao } = {0}

Коэффициенты канонических уравнений (элементы матрицы [A]и вектора {Ao }) определяются выражениями:

a11

=

1

×

li

, a22

=

1

×

li

, a33

=

1

×

li

,

Go

di ред.

Go

di ред.

Go

di ред.

 

 

ABCDEFA

 

 

BCDEFB

 

 

CDEC

 

(здесь суммирование ведется по панелям, где qI , qII , qIII не равны нулю соответственно),

a12

= a21

=

1

×

 

li

 

 

,

a23

= a32

=

1

×

li

,

Go

 

di ред.

Go

di ред.

 

 

 

BCDEF

 

 

 

 

 

 

 

CDE

 

 

 

 

 

a13

= a31

=

1

×

 

li

,

 

 

 

 

 

 

 

 

 

G

 

d

i ред.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CDE

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

o

 

 

 

 

 

 

 

 

(здесь суммирование ведется по панелям, где qI , qII ; qII , qIII ; qI , qIII - соответственно не равны нулю),

a10

=

1

×

li × qoi

, a20

=

1

×

li × qoi

, a30

=

1

×

li × qoi

,

 

 

 

 

 

 

di ред.

 

 

 

 

Go ABCDEFA

di ред.

 

Go BCDEF

 

 

Go CDE di ред.

(здесь суммирование ведется по панелям, где qI ,q0i ; qII , q0i ; qIII , q0i - соответственно не равны нулю). Здесь li -длина i-той панели; Go - приведенный модуль

сдвига (для обшивки из дюраля G = 105

кг × см− 2 ) ; d

i ред.

-

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

o

 

 

 

 

щина

обшивки di ред. = di × Y ;

 

Ψ - редукционный коэффи-

 

циент обшивки.

 

 

Модуль сдвига обшивки

 

панели крыла не равен модулю

 

сдвига материала обшивки, а

 

зависит еще от ее кривизны,

 

толщины, шага нервюр и

 

стрингеров (размеров подкреп-

 

ляющей

 

клетки), подкрепля-

ющих профилей, характера нагружения пластины. Значения модуля сдвига бо-

26

лее или менее точно определяются опытным путем для данной конструкции. В расчете приходится большей частью пользоваться средними величинами G, полученными из испытаний аналогичных конструкций. Так как

Y = G ,

Go

то при вычислении мы будем пользоваться значениями редукционных коэффициентов, приведенными на рис. 15. Значения коэффициента Ψ для обшивки из другого материала следует умножить на ϕE .

Определяем погонные сдвиги в панелях как соответствующую сумму погонных касательных усилий в открытом контуре и замыкающих интенсивностей:

qΣсдв. = qoi + qI + qII + qIII

По результатам расчета строим схему потоков погонных касательных усилий по контуру сечения.

8.Расчет сечения крыла на кручение

8.1.Определение положения центра жесткости сечения крыла

Положение центра жесткости определяется по формуле

 

2

n

 

xцуточн.ж. . =

× (qI × FI + qII × FII + qIII × FIII + qoi

× wi ),

QΣ

 

i =1

 

где FI , FII , FIII - площади контуров ABCDEFA, BCDEFB, CDEC соответственно

i

Уточненное

QΣ положение ц.ж.

i-1 i

ri

 

 

 

 

 

ц. ж.

A

 

x

уточ.

 

x

 

 

 

ц.ж.

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 16

(рис. 15), которые подсчитываются по чертежу сечения крыла; qI , qII , qIII - потоки погонных касательных усилий, полученные в результате расчета сечения на сдвиг от силы QΣ ; qoi - потоки погонных касательных усилий в открытом контуре сечения крыла от сдвига; ωi - секториальная площадь, соответствующая i – той панели (рис. 16). Приближенно значения ωi

можно вычислить как площадь треугольника

27

wi » 1 × ℓi × ri ,

2

где ri - длина перпендикуляра, опущенного из произвольно выбранного полюса А (в качестве такого полюса можно взять координату предварительно определенного центра жесткости) на хорду дуги элемента контура сечения i .

 

 

8.2. Определение потока касательных усилий от кручения

 

При

расчете на кручение замыкающие потоки

касательных усилий

qкр.

, qкр. , q

кр. определяются из системы уравнений

 

I

II

III

 

 

 

[ A]×{q} = 2 ×x{F} ,

 

 

 

2 ×{F}т ×{q}= M кр.уточн.

(13)

 

 

 

Здесь компоненты вектора {F} - FI , FII , FIII - площади

контуров ABCDEFA,

BCDEFB, CDEC (см. рис. 15); ξ - относительный угол закручивания сечения;

M кр.уточн. - крутящий момент относительно уточненного положения центра же-

сткости сечения, определяемый из выражения

 

 

 

M кр.уточн. = M кр. ± QΣ × d

 

Здесь d

расстояние между приближенным и уточненным положениями цен-

тров жесткости сечения крыла.

 

 

Значения коэффициентов aij канонических уравнений (13) те же, что при

расчете на сдвиг. После определения потоков замыкающих касательных усилий qIкр. , qIIкр. , qIIIкр. при кручении, суммарные погонные сдвиговые усилия находим подобно расчету на сдвиг, положив qoi = 0 , т.е.

qiкрΣ . = qIкр. + qIIкр. + qIIIкр.

По результатам расчета строим суммарную эпюру потоков погонных ка-

сательных усилий от сдвига и кручения qiсдвΣ + кр по контуру расчетного сечения крыла. При построении суммарной эпюры положительные значения потоков откладываем внутрь контура сечения.

28

9.Проверка обшивки и стенок лонжеронов на прочность

иустойчивость

Врезультате проверочного расчета должно быть дано заключение о прочности подобранного сечения крыла. Для этого обшивка и стенки лонжеронов проверяются на прочность и устойчивость.

Максимальные нормальные напряжения, действующие на соответствую-

щую панель обшивки (или стенки лонжерона) с учетом ϕE

smax = s

×

E

обш

,

 

обш

i max

E

 

 

 

 

стр

а значения редукционного коэффициента обшивки находятся по выражению

ϕобш = bприс. ,

b

где bприс . = 30 × d - для сжатой зоны; bприс. = 0,8 × ℓ- для растянутой зоны; - ширина рассматриваемой панели обшивки (шаг стрингеров). Тогда средние нормальные напряжения в панелях обшивки

sср. обш. = jобш × smax . обш

Касательные напряжения, действующие в обшивке (или стенке лонжерона) от сдвига и кручения, вычисляются как

tiΣ

=

qiсдвΣ

+ кр

 

 

.

 

 

 

 

di

Критические касательные напряжения tкр вычисляются аналогично sкр по

формуле

 

 

 

 

 

 

 

 

 

σср.обш.

 

1 + n*

 

 

 

 

 

tкр

 

 

 

 

 

 

= tв ×1 + n* + n* 2 ,

 

τ

 

 

где tв

= 0,7 × sв , n*

=

τв

,

 

a

 

 

 

 

tкр.э

 

 

 

 

tкр.э

d

 

2

 

 

 

 

= 0,9 × E × kτ ×

 

i .

 

 

 

 

 

 

b

 

 

 

 

Значения коэффициента kτ берется из ра-

 

 

 

боты [8] для пластины шарнирно опертой по

 

Рис. 17

 

контуру.

 

 

 

 

 

 

 

Для оценки устойчивости элементов кры-

 

 

 

 

 

 

29

 

 

 

 

 

ла вычисляется коэффициент

 

 

sср.обш.

 

t

iΣ

2

k

 

=

 

 

+

 

.

 

 

 

 

 

кр

 

sкр

 

 

 

 

 

 

 

 

tкр

Значения коэффициента kкр позволяют судить о работе обшивки (стенки

лонжерона) на устойчивость при одновременном действии сжатия и сдвига (рис. 17). Растянутая панель находится в этом случае в облегченных условиях и в данном случае не рассматривается. Если kкр < 1, то потери устойчивости не

произойдет, при kкр > 1 пластинка потеряет устойчивость. Потеря устойчивости

не является критерием потери несущей способности конструкции.

При проверке обшивки на прочность вычисляются значения коэффициента kпр с использованием четвертой теории прочности:

kпр

=

s2max + 3 × ti2Σ

обш

sв

 

 

£ 1 ,

(14)

где s2max + 3 × ti2Σ = sэквIV . .

обш

Для стенок лонжеронов σmax = 0 (чистый сдвиг) коэффициент kпр вычис-

ляется по формуле (14) до потери устойчивости и, если стенка потеряла устойчивость, то по выражению

kпр =

2 × tΣ

£ 1.

sв

 

 

Значения коэффициента kпр < 1 позволяют сделать вывод о том, что условие прочности соблюдается.

Ли т е р а т у р а

1.Бадягин А.А., Егер С.М. и др. Проектирование самолетов. – М.: Машиностроение, 1972. 516с.

2.Зайцев В.Н. Рудаков В.Л. Конструкция и прочность самолетов. – Киев: Вища школа, 1976. 400с.

3.Шульженко М.Н. Конструкция самолетов. – М.: Машиностроение, 1971.

4.Кравец А.С. Характеристики авиационных профилей. – М.: Оборонгиз, 1939.

5.Макаревский А.И., Корчемкин Н.Н., Француз Т.А., Чижов В.М. Прочность самолета. – М.: Машиностроение, 1975. 280с.

30