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

книги из ГПНТБ / Образцов И.Ф. Методы расчета на прочность кессонных конструкций типа крыла

.pdf
Скачиваний:
30
Добавлен:
30.10.2023
Размер:
14.29 Mб
Скачать

(в направлении оси z) и двумя в своей плоскости. При этих степе­ нях свободы для рассматриваемой задачи получим три дифферен­ циальных уравнения равновесия:

aU" — b1 U— b2V — b^' -= 0;

b2U' + ^9" 4-62х" = 0;

(1.242)

—ct*=o,

 

где

 

— с -j— t.

(1.243)

При получении этих уравнений мы старались

придерживаться

обозначений, принятых В. 3. Власовым в цитированной здесь книге. Поэтому уравнения (1.242) отличаются от полученных им уравне­

ний только

значением коэффициента ct. Здесь

ct— обобщенная

жесткость оболочки и нервюры.

 

Введем в

рассмотрение новую функцию f(z)

и выразим через

эту функцию и ее производные искомые обобщенные перемещения U(z), 9 (z) и х (z) так, чтобы первое и третье уравнения (1.242) удовлетворялись тождественно при любом выборе функции f(z).

Тогда будем иметь

(1.244)

Имея в виду условия (1. 187) и (1. 188), а также наличие новой

обобщенной силы (поперечного бимомента)

Р, которая находится

по формуле

 

 

 

 

 

 

(j) тф3 dF,

 

(1.245)

 

 

 

 

и раскрывая эти выражения,

получим

 

 

 

B——aUr;

 

 

 

Н=Ь2и + Ь^ + Ь2Р-

 

(1.246)

Согласно (1.244)

 

 

 

 

ctb2 1

277

b2 7

b2

(1.247)

 

P=af”.

88

Систему уравнений

(1.242)

при

помощи (1.246)

приводим

к одному уравнению относительно

новой функции f(z) вида

где

 

 

 

 

(1.248)

 

 

 

 

 

 

2 (&?-&!)

(1.249)

 

 

 

Д=^-.

 

(1.250)

 

 

а

 

 

Представим общий

интеграл

дифференциального

уравнения

(1. 248) в виде

 

 

 

 

 

/= С1Ф1-|-С2Ф2-}- С3Ф3 + С4Ф44- С-г-|- С6.

(1.251)

Здесь

®j = ch аг-sin г;

 

 

 

 

Ф2=ch аг-cos 8г;

(1.252)

 

®3 = sh аг-cos £г;

 

 

 

®4=sh аг-sin г;

 

 

 

 

 

 

(1.253)

Производные от функций (1.252) даны в табл. 5.

через на-

Выражая постоянные интегрирования Су, С2 ... С6

чальные параметры Uo,

0О> х0 ,

Но,

Во,

Ро, относящиес;I

к сечению

г=0, решение задачи приводим к интегралам, записанным в очень удобной форме. Эти интегралы В. 3. Власов представил в виде табл. 14.

В пересечениях строк и столбцов таблицы вписаны функции, , зависящие от аргумента z и от обобщенных геометрических и упругих характеристик кессона.

Определение постоянных интегрирования

При помощи интегралов, записанных в табл. 14, решение задачи приводится к отысканию только трех постоянных интегрирования, остальные же три в каждой частной задаче будут определяться из граничных условий.

Для рассматриваемого кессона (см. фиг. 35) начало осей коор­

динат поместим на свободном конце. Тогда

 

при

z=0 Д0 = О; хо = О и HQ=H,

(1.254)

т. е. на свободном конце бимомент Во==0, крутящий момент Н0=Н, а Н — постоянный по всей длине оболочки. Принимаем, что крутя­ щий момент приложен к жесткой в своей плоскости передней нер­ вюре, поэтому поперечная деформация в плоскости этой нервюры Хо='О.

89'

 

 

 

 

 

 

 

 

 

 

Таблица li

 

*0

 

и0

00

 

Яо

 

Во

Ро

 

X

f*

 

 

0

 

72

 

 

1

 

 

Ф2 +----- Фд

4

“й(аФ1+ Фз)

 

-Й(«Ф1+?Ф3)

 

-------- Фл

2аМ4[(“3-3^Ф> +

 

2Г2а?

 

 

 

 

 

 

 

 

 

 

 

2afi«

 

 

 

 

 

 

 

 

 

 

 

 

+ (За2 - 3) ф3]

и

S2

 

г2

 

-^(2а -2а Ф2 + г2ф4)

 

 

 

 

^(«Ф1Фз)

Ф2 — — Ф4

0

2а аз2[(“3-3а 2)Ф1-

— Ф4

 

 

 

 

 

 

 

 

2 2а 4

 

 

 

 

 

 

2а л 4

 

 

 

 

 

 

 

 

 

 

-(За2р -Р2)ф3]

 

 

0

— Д272®4

 

Д172(«Ф1Ф3)

1

д172 (аф1 — Зфз) + 71-г

-^W-2«3®2+^4)

- («*1 + Ф3)

 

 

 

 

 

 

 

 

 

 

zip

 

н

0

 

0

0

 

 

1

 

0

0

 

в

— Д2Фд

 

1(а®1— ?Ф3)

0

 

А172 (a®i — Ф3)

 

«’=-Д«4

_Х(аФ1 + Фз)

р

Д2(®Ф1 + Ф3)

— Д2Ф4

0

 

— Д272$4

 

S2

г2

Фд

 

 

М(аФ1 “ Фз)

Фо+---

 

 

 

 

 

 

 

 

 

 

2~2сф

4

 

 

 

Здесь 71 — ,2

. 2 ;

72 — .2

b2

as2

as^

 

 

 

 

 

. 2 ; Д1 —

9 R >

Д2 = 9 „ ;

 

 

 

 

 

 

*1~ *2

о1

— b2

2з|3

2сф

 

 

 

 

 

« = ^ф22(Л + /?2 + бДЛ);

bl = Y G (d2F2 +

 

Ь2 = 1 О (dfa-d2^.

 

В силу сказанного выше решение рассматриваемой задачи при­ водится к трем интегралам (выписываем их из табл. 14):

£/ — UJ^uu (г) НКин & + Р(fiup (z);

 

* = — ЦЛ* и (z) — HKzh (z) 4- РйК. р (z);

- (1.255)

^ = ийК.и (z) 4-604- НК. я (г) — РаК.р (г).

 

Для сечения z=l коэффициенты будут иметь вид

 

Кии (0 = ch «Zcos pZ —— sh aZsin Z;

 

 

KUH

(2a —2ap ch aZcos pZ4-r2sh aZsin pZ);

 

 

2сф

 

 

 

 

Л'Ур(/) = -^-5Ь aZsin pz;

 

 

 

 

Kz

(ach »Zsin pZ4- p sh aZcos pZ);

(1.256)

 

Kz h(.1)—~ (“ ch al sin pZ4- p sh aZ cos BZ);

 

 

2ap

 

 

 

 

Kz p

4 l(a3—-Зар2) ch al sin pZ 4-

 

 

4- (3a2p — p3) sh al cos pZ].

 

J

 

 

 

 

 

Постоянные интегрирования, входящие в

(1.255), определяем

из граничных условий при z=l:

 

 

 

 

 

U = Q-,

 

 

 

 

z = 0;

 

 

 

(1-257)

 

0 = 0.

 

 

 

 

Указанный край оболочки закреплен от продольных и попереч­

ных перемещений,

поэтому граничные

условия (1.257)

являются

чисто геометрическими.

 

 

совместно

 

Имея в виду условия (1.257)

и

решая

равенства

(1.255), получим для постоянных

интегрирования выражения:

Р<г=

И [Кин (0 Кг v (/) + Кг н (/) Кии (01

(1.258)

 

 

 

 

КирУЖ^ЦУ + К^УЖии (/)

и0 _ НКин (0 — PqKup (/) .

(1.259)

Кии(!)

 

'Кцр (0 Ки е (0

l-A^ep(Z)

Kuu (!)

 

Kuh (!) Ks и (I)

(1.260)

+ К*нМ ■

Kuu (!)

91

Формула для определения нормальных напряжений

Нормальные напряжения определяем из равенства (1. 185). Дифференцируя первое из выражений (1.255) и имея в виду (1.258) и (1.259), получим формулу для определения нормальных напряжений в любой точке кессона с учетом деформируемости кон­ тура поперечного сечения:

a(z, 5)=Яъ£[фз(£--а)+ф1(^-+ ) Ъ^) + иоЕ X

х[ф3('а-Л-Ф]( + -^)1?3(5)+Р0^(,^+^)?3(5). (1-261)

Геометрические и упругие характеристики этой формулы пред­ ставлены равенствами (1.249), (1.250), (1.252), (1.253) и табл. 14.

Фиг. 37. Графики (распределения бимоменггных нор­ мальных напряжений по размаху кессона при раз­ личных жесткостях нервюр.

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

Пусть кессон, показанный на фиг. 37, закручивается моментом Н=50 000 кгсм. Материал оболочки — дуралюмин.

92

1. Вначале будем считать,

что

оболочка не имеет

поперечных

подкреплений (нервюр).

В

этом

случае

обобщенная жесткость

ct = c и определяется по формуле

(1.241). График

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

нормальных бимоментных

напряжений вдоль поясов

кессона,, вы­

численных ПО1 этой формуле,

показан на

фиг. 37

сплошной ли­

нией 1.

 

 

 

.

 

2. Примем теперь контур

поперечного

сечения

оболочки неде-

формируемым, т. е. положим, что Л=/2=о°. Нормальные напряже­ ния определим по формуле (1.209). График распределения напря­ жений вдоль поясов показан на фигуре сплошной линией 4.

3. Пусть оболочка подкреплена нервюрами в виде пластин тол­ щиной 8* =0,05 см, шаг между нервюрами Ь = 20 см. Обобщенную жесткость ct находят из выражения (1.243). Напряжения, вычис­

ленные по формуле (1.261), представлены на фиг. 37 штрихпунктирной линией 2.

4. Примем толщину нервюр S*=0,l см, шаг между

ними &=

= 20 см. График распределения нормальных напряжений,

получен­

ных для этого случая по формуле (1. 261), нанесен на фиг.

37 пунк­

тирной линией 3.

 

При толщине стенок нервюр 8*=0,05 см жесткость поперечных

сечений оболочки незначительная, поэтому бимоментные напряже­ ния, показанные кривой 2, имеют небольшую величину. С увеличе­ нием толщины стенок нервюр до 0,1 см увеличивается жесткость кессона и бимоментные напряжения сильно возрастают. При этом график бимоментных напряжений близко подходит к кривой 4, по­ лученной для оболочки с жестким контуром поперечного сечения.

В заключение отметим, что пластинчатые нервюры при сравни­ тельно небольших толщинах придают оболочке большую попереч­ ную жесткость. Жесткость кессона зависит также от шага нервюр.

Глава VI

ИЗГИБ И КРУЧЕНИЕ КЕССОНОВ С УЧЕТОМ УПРУГОСТИ ЗАДЕЛКИ

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

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

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

93

Фиг. 38, Расчетная модель кессонов, соединенных с оболочкой центроплана.

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

Ввиду симметричности расположения оболочек относительноплоскости симметрии самолета расчетную модель можно предста­ вить состоящей из прямоугольной оболочки — центроплана длиной. li— и прикрепленного к ней кессона.

§ 15. Изгиб и кручение кессона с упругой заделкой сосредоточенными на свободном конце поперечной силой, изгибающим и крутящим моментами

Пусть на кессон (фиг. 38, б) действуют поперечная сила Q и изги­ бающий и крутящий моменты М и Н, приложенные к передней: жесткой в своей плоскости нервюре. Так как кессон связан с обо­ лочкой центроплана, то перемещения, возникающие от действия внешних сил в кессоне, будут передаваться на оболочку центро­ плана.

Представим искомые обобщенные продольное и поперечное пе­ ремещения точки Mi(z, s) элементарной полоски кессона в форме

и(г,

s) = ^1(2)®1(s) + t72(z)<?2(s)4-Z73(z)<P3(s);

(1.262)

 

v(z, s) =

(z)^j(s)+ n2(z)<p2(s),

(1.263);

а искомые продольное и поперечное

перемещения точки

A42(z, s)

элементарной полоски оболочки центроплана в виде

 

u(z, ^) = t/j (г) cpj (s)-f-t/2(z) <p2(s)-f-(73 (z) tp3(s);

(1.264>

 

ц(г, s)=V1(z)^(s)-b/2(z)^(s).

(1.265)

Начало координат для кессона

поместим на свободном, конце-

в плоскости

передней нервюры, а для оболочки центроплана — в.

плоскости симметрии самолета.

 

 

 

Выбор аппроксимирующих

функций

Функции

(s), <p2(s),

срз (s), <J»i

(s) и ф2(«)

описаны

в гл. I и

представлены на фиг. 3.

Оболочка центроплана может иметь иные размеры и геометри­ ческие характеристики, чем кессон, поэтому ее аппроксимирующие функции будут отличаться по величине от функций, принятых для кессона.

95-

.Эти функции выберем так: ?i(s) ==?(«);

?2 (S)

(1.266)

?3(s) = ^(s)^(s);

Ш=У(«);

ip2(s)=A(s),

где ?! (s) —функция, соответствующая перемещению сечения по закону плоскости (см. фиг. 39, а);

?2(s) —функция, относящаяся к депланации сечения при изги­ бе поперечной силой (см. фиг. 39, б).

Фиг. 39. Эпюры аппроксимирующих функций для оболочки центроплана.

Коэффициент ортогональности с определяется по формуле

с —

(1.267)

6-Zr

96

Здесь

 

 

 

 

/72 = 82t/2;

 

 

 

J

x

1 \ 6

2

+ дА

 

 

 

 

/

(1.268)

iZb

d2 — высота и

ширина оболочки центроплана;

 

®з($) — функция,

соответствующая обобщенной

депланации,

 

 

возникающей при кручении (фиг. 39, в);

 

% (5) —функция,

соответствующая поступательному переме-'

_

_

щению элементарной полоски в вертикальном направ-

лении (фиг.

39, г);

 

 

.

Ф2($) —функция,

отражающая

поворот поперечных сечений

_

_

оболочки под действием крутящего момента (фиг. 39, д);

h(s')—длина перпендикуляра,

опущенного из

начала коор­

динат на соответствующую пластину оболочки центро­ плана.

О вычислении коэффициенте!в дифференциальных уравнений

Коэффициенты дифференциальных уравнений для кессона и оболочки центроплана вычисляются по формулам (1.8). Они пред­ ставлены выражениями (1.26), (1.28).

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

Дифференциальные уравнения для консольной

ицентропланной оболочек

Всоответствии с искомыми обобщенными продольными и попе­ речными перемещениями, представленными в форме разложений (1.262) — (1.265), элементарные поперечные полоски рассматривае­

мых оболочек в отношении перемещений из плоскости поперечного сечения будут обладать тремя степенями свободы, соответствующи­ ми трем обобщенным координатам cpi(s), <p2(s) и <p3(s), а в отно­

шении контурных перемещений

точек полоски в плоскости попе­

речного

сечения — двумя степенями свободы, соответствующими

обобщенным координатам ipi(s)

и ф2(«).

 

При этих степенях свободы

(от=3, га=2), заданных

внешних

нагрузках

и вычисленных коэффициентах (1.2'6), (1.28)

система

дифференциальных уравнений (1.8) примет вид:

7

428

97

Соседние файлы в папке книги из ГПНТБ