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

Фролов ЭM.Динамика и прочность машин.Теория механизмов и машин

.pdf
Скачиваний:
28
Добавлен:
23.02.2023
Размер:
26.85 Mб
Скачать

310

Глава 5.2. МЕХАНИКА СЛОИСТЫХ КОМПОЗИТОВ

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

Коэффициенты D определяют изгибные жесткости стенки. В частности, D\\ и D22 соот­ ветствуют изгибу в плоскостях XZ и yz, а D^^ - кручению. Коэффициент D\2 отражает связь между изгибными деформациями в плоскостях XZ и yZy обусловленную эффектом Пуассона, а Z)i6 и D2e - связь между изгибом и кручением. Из третьего равенства (5.2.5) следует, что рас­ сматриваемые коэффициенты зависят от коор­ динаты е, т.е. от положения базовой плоскости, к которой приведены моменты. Для того чтобы получить истинную изгибную жесткость, следует рассмотреть действие только одного момента и задать координату е с помощью соответствующе­ го уравнения (5.2.9). В результате изгибные (т/и=11, 22) и крутильная (тт—66) жесткости будут иметь вид

(1)1 /г(0)

тттт

Инаконец, коэффихщенты К (5.2.7) опре­ деляют жесткость стенки при межслоевом сдвиге

вплоскостях XZ и yz- Для однородной (однослойной) стенки эти жесткости имеют слЬ-

дующий вид: К-^ = G^h; К2 = Gy^h, Для трех­ слойной стенки, у которой податливость при сдвиге определяется слоем изотропного легкого

заполнителя, К^ = К2 = G^h jh^ , где GQ И

Ло - модуль сдвига и толщина слоя заполнителя. Соотношения (5.2.5) и (5.2.6) в совокупно­ сти с уравнениями равновесия или движения, а

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

L В результате решения исходной системы уравнений находят перемещения и деформации

влюбой точке конструкции.

2.По формулам

^1

- ^ 1

1^1

+^12^2

Л')

F(Or (О .. (0J0

 

^12

- "^12

^12

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

Более полная информация о методах расче­ та композитных элементов конструкций пред­ ставлена в работе [6].

5.2.2. ДИССИПАТИВНЫЕ СВОЙСТВА СЛОИСТЫХ КОМПОЗИТОВ

Воспользуемся результатами, представлен­ ными в п. 5.1.7 для монослоя, и рассмотрим диссипативные свойства слоистых материалов.

При идеальной взаимосвязи слоев потери энергии в многослойном композите при цикли­ ческом нагружении равны сумме потерь в моно­ слоях. Величину потерь за цикл нагружения в /-М слое определяют с помощью матрищя

^монослоя через амплитудные значения

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

(5.1.100):

 

/=1

vV)Y'Y'

.W-^^\o

 

l

^к f

(/^'^

где h^ '

относительная толщина /-го слоя

[й^^-н^Чн],

или

i»'=i{<,j'[»|,j,

(О 2

(/)

2

(/•)

4 " .

4'»,2 .

(/)

2

('•)

S') т. .»?)

Щ"1 +^ху

С учетом равенств (5.2.1), где принимают r=Z/, находят деформации в координатах, связанных с направлением армирования /-го слоя (см. рис. 5.1.1).

3. С помощью закона Гука

где

[*]=[s]t[e<"p'>p"p'>p]

(5.2.10) матрицы жесткости и податливос-

• n - l'

ти пакета слоев

СВОЙСТВА КОНСТРУКЦИОННЫХ композиционных МАТЕРИАЛОВ

311

Для частных случаев структуры пакета сло­ ев многослойного материала формула (5.2.10) может быть упрощена. Так, при одноосном на­ гружении перекрестно армированного композита

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

/ / ; /

/2

,M^llUll^22-^12

^22

2^16^16^22 --^26^1^2в^п)

+

^

;

-^

 

^22

 

4^66^61^22

• ^26^12 J

/

/ /

/

Л2

^22(^^11^22

-^12

Здесь штрихами отмечены компоненты матриц жесткости и диссипативных характерис­ тик монослоя. На рис. 5.2.2 представлены зави­ симости от угла армирования коэффихщента

рассеяния ^х ~^1 ^ модуля упругости

*

Е^ - Е-^ при одноосном нагружении симмет­ рично армированного композита.

ных материалах. Величины Ei и Е2 являются модулями упругости при нагружении соответ­ ственно в продольном и поперечном направле­ ниях. Величины Qj и G2 определяют пределы прочности при растяжении в этих направлениях, а а2 и а2 - пределы прочности при сжатии.

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

Х(Т]1Х(20''С}

W

^7

0,5

\ Л / 4

 

 

 

 

 

 

5 А ^

 

 

20

100

200

ЗОО^С

Рис. 5.2.3. Зависимость относительных механических характеристик однонаправленного стеклопластика от температуры:

7-^1(7); 2-GÎ2( 7); J-a 1^(7);

^ai(7);J^a2(7);(^i2(7)

К{1)1Х(20Ч)

у1

О

15

30

45

60

75

(р"

Рис. 5.2.2. Зависимость коэффицисЕгга диссипации i|/^ и модуля Ех от угла ориентации волокон ф при £*!=103,4 ГПа; ^=7,2 ГПа; Gi2=3,8 ГПа; vi2=0,29; vi/j =0,45 %;

V|/2=5,5 %; \|/g=6,7 %

5.2.3. СВОЙСТВА КОНСТРУКЦИОННЫХ КОМПОЗИЦИОННЫХ МАТЕРИАЛОВ

В табл. 5.2.1 представлены типовые меха­ нические и физические характеристики наиболее распространенных композиционных материалов. Индекс 1, как и ранее, соответствует направле­ нию волокон в однонаправленном материале или направлению основного армирования в тка­

0,5

\ х ^

\

 

 

 

 

8

^

 

 

20

100

200

ЗОО'С

N

Рис. 5.2.4. Зависимость относительных механических характеристик однонаправленного органопластика от температуры:

l-Ei(T); 2-Б2(Т); 3-Gi2(T); 4-^1(1);

5-aï(7); 6-a2{T); l-G^iT)- »^niT)

 

Вид

Марка

Материал

арми­

напол­

 

рования

нителя

Стекло­

Тканые

1СУ8/

пластики

 

З-ВМ-78

 

 

Т-25(ВМ)

 

 

-78

 

 

Т-10

 

 

НП-750

 

Одно­

ВМПС

 

направ­

 

 

ленные

РВМН

 

 

РБН

 

Хаоти­

БС

 

чески

 

 

армиро­

БС

 

ванные

БС

 

 

Угле­

Одно­

ВМН-4

пластики

направ­

УКН-

 

ленные

 

 

5000

 

 

ЭЛУР

 

 

ЛУ-П

Органо­

Одно­

СВМ-У

пластики

направ­

жсвм

 

ленные

 

 

Армое

 

Тканые

Арт.

 

 

56313

 

 

Арт/

 

 

56334

 

5.2.1. Типовые свойства композиционных материалов

 

 

 

Марка

El

El

^1

^2

^1

 

^2

т

а10б,

 

Ср.

Р,

связую­

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

щего

 

ГПа

 

 

МПа

 

 

 

к-1 1

Вт/(мК)

кДж/(кг-К)

кг/м^

ЭХД-У

22

28

400

400

280

 

350

32

6,2

0,32

1,5

1850

ЭП-5122

50

18

800

90

500

 

200

30

12,0

0,40

1,1

1850

ЭП-5122

26

20

560

310

260

 

220

28

8,0

0,33

1,1

1900

ЭХД-У

26

22

240

260

220

 

240

45

6,5

0,29

1,1

1800

эхд-

55

3,0

1800

12

600

 

100

32

20

0,4

1,2

2100

МК

60

3,5

1500

12

600

1

100

35

15

1 0,4

1,4

2050

ЭДТ-10

ЭДТ-10

40

3,0

1300

10

400

1

100

30

20

0,4

1,4

2100

Р-2

18

18

80

80

130

 

130

45

20

0,35

1,2

1950

ФН

18

18

500

500

140

 

140

50

25

0,36

1,3

1950

Р-2М

18

18

100

100

130

,

130

50

18

0,40

1,3

2000

УП-632

150

3,5

900

10

600

 

100

25

-0,4

0,50

0,72

1500

эхд-

120

3,0

1300

10

800

 

70

30

0,2

0,53

0,77

1540

МК

130

3,0

800

И

90

 

80

30

1 -0,1

0,47

0,76

1480

ПАИС

 

ЭНУП

130

2,5

800

11

600

 

80

25

0,2

0,48

0,75

1490

эхд-

85

3,0

2500

10

210

 

100

32

1 0,8

0,15

0,80

1300

мк

80

3,0

2200

^\%

170

 

100

35

1,0

0,15

0,80

1320

эхд-

 

мк

95

3,5

2700

12

346

 

105

39,2

-0,5

0,17

0,95

1300

эхд-

 

мд

 

 

 

 

 

 

 

 

 

 

 

 

эдт-

34

34

600

500

150

 

160

44

0,77

0,13

1,1

1250

10П

30

45

550

700

140

 

180

42

0,75

0,13

1,1

1250

ЭДТ-

 

10П

 

 

 

 

 

 

 

 

 

 

 

 

5.2.2. Механические характеристики конструкционных композитов при температуре 20° С

 

 

Материал

р-10-3,

El

El

Gn

Vl2

^r

a j

^2

42

ai-106

а2-10б

 

 

 

 

кг/м^

 

ГПа

 

 

 

ГПа

 

МПа

 

 

К-1

 

 

 

 

 

 

 

 

 

Однонаправленный

2,05

60-70

11

0,5

0,28

1,75-2,1

0,4

27

72

36

8

100

стеклопластик (v=0,75)

 

 

 

 

 

 

 

 

 

 

 

 

Однонаправленный

 

 

 

 

 

1,4-2,2

 

12-28

97-138

20-44

 

 

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

1,30-1,38

78-95

4,1-5,5

2,1

0,29

0,28-0,31

-2,5

100-200

ных волокон (v=0,65)

 

 

 

 

 

 

 

 

 

 

 

 

Однонаправленный

1,5

140-

5-6

5,5

0,32

1,1-1,6

0,5-1

20

75

30-50

-1,5

40

углепластик (v=0,5)

 

160

 

 

 

 

 

 

 

 

 

 

Однонаправленный

1,9

250

25

9,8

0,22

1,2

1,16

20

70

60

2,4

100

боропластик

 

 

 

 

 

 

 

 

 

 

 

 

Однонаправленный

2,6

218

141

42

0,23

1,52

1,76

141

211

162

0,65

2,78

бороалюминий (v=0,4)

 

 

 

 

 

 

 

 

 

 

 

 

Однонаправленный

 

108

 

2

 

 

 

 

 

 

-2

 

углеорганопластик (50 %

1,5

4-6

0,4

1,2

0,7

25

100

56

40

углеродных и 50 %

 

 

 

 

 

 

 

 

 

 

 

 

органоволокон) (v=0,6)

 

 

 

 

 

 

 

 

 

 

 

 

Стеклотекстолит на основе ткани

1,8

19,2

24,2

7

0,14

0,38

0,22

50

270

85

9,5

9,5

сатинового переплетения (v=0,65)

 

 

 

 

 

 

 

 

 

 

 

 

Углепластик на основе

1,39

12,4

11

5,8

0,16

0,15

0,165

95

120

20

13,3

16,6

углеткани (v=0,57)

 

 

 

 

 

 

 

 

 

 

 

 

1Свазиизотропный хаотически

2

16,3

16,3

6,5

0,27

0,13

0,25

130

250

127

-

-

армировашпой стекловолокнис-

тый пресс-материал (v=0,62)

Способ

изготов­

ления

 

О

 

со

Намотка

о

О

 

H

»о

H

»s

о

д

д

Вык­

ладка

д

Прессо­ о

вание со

д

д

д

о

Намотка дд

Прессо­

вание

»

»

314

СПИСОК ЛИТЕРАТУРЫ

Х(Т)/Х(20''С)

метров окружающей среды наибольшее влияние

на прочность композитов с полимерной матри­

 

цей оказывает температура. Соответствующие

 

зависимости для типовых композиционных ма­

\ / ^

териалов (см. табл. 5.2.2) представлены на рис.

5,2.3 - 5.2.6.

0,5

 

 

 

 

2^

N ^ ^

 

 

СПИСОК ЛИТЕРАТУРЫ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/4

/

NJ

 

1.

Бакулин В. Н., Ряссоха А, А. Метод ко­

 

 

 

 

 

 

 

 

нечных элементов и голографическая интерфе­

 

 

 

 

 

 

 

 

рометрия в механике композитов. М.: Машино­

 

 

 

 

 

 

 

 

строение, 1987. 311 с.

 

 

20

 

100

 

200

 

500"С

 

2. Булаве Ф. Я., Радинып И. Г. Деформа-

 

 

 

тивные свойства однонаправленно армированно­

 

 

 

 

 

 

 

 

го пластика при трансверсальном нагружении //

Рис. 5.2.5. Зависимость относительных механических

Вопросы динамики и прочности. Рига: Зинатне,

1980.

С. 73-81.

 

 

характеристик стеклотекстолита от температуры :

 

 

7-^1(7); 2~oliT);

3-G^(T);

^а^СТ);

5-GI(T);

 

3. Булаве Ф. Я., Радиньш И. Г. Микроме­

ханика ползучести однонаправленно армирован­

 

 

 

6-1^2^7)

 

 

 

 

ных пластиков при продольном сдвиге //

Меха­

 

 

 

 

 

 

 

 

ника армированных пластиков. Рига: Рижский

 

 

 

 

 

 

 

 

политех, ин-т, 1981. С. 19-26.

 

Х(Т)/Х(20''С)

 

 

 

 

 

 

 

4. Ван Фо Фы Г. А. Теория армированных

 

 

 

 

 

 

 

 

материалов. Киев: Наукова Думка, 1971. 232 с.

Ю

 

 

 

 

 

 

 

 

5. Ванин Г. А. Микромеханика композици­

 

 

 

 

 

 

 

онных материалов. Киев: Наукова Думка, 1985.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

302 с.

 

 

 

 

 

 

 

 

 

 

 

 

6. Ваеильев В. В. Механика конструкций

 

 

 

 

 

 

 

 

из. композиционных материалов. М.: Маши­

 

 

 

 

 

 

 

 

ностроение, 1988. 270 с.

 

 

0,5

 

 

 

 

 

 

 

 

7. Викарио А., Толавд Р. Критерии прочно­

 

 

 

 

 

 

 

сти и анализ разрушения конструкций из компо-

^

 

 

^

 

 

 

—d

JA4

 

 

 

ЗИ1Д10ННЫХ материалов // Композиционные ма­

 

 

 

 

 

 

 

 

териалы / Под ред. Л. Браутмана и Л. Крока.

 

 

 

 

 

 

 

 

Г.7. М.: Машиностроение, 1978. С. 62-107.

 

 

 

 

 

 

 

 

 

8. Гольденблат И. И., Копнов В. А. Крите-

 

 

 

 

 

 

 

 

рии

прочности и пластичности конструкцион-

 

 

500

1000

 

1500

2000'C ных

материалов. М.:

Машиностроение,

1968.

 

 

 

 

 

 

 

 

191с.

9. Зиновьев П. А., Ермаков Ю. Н. Анизот­

 

 

 

 

 

 

 

 

 

Рис. 5.2.6. Зависимость относительных механических

ропия диссипативных свойств волокнистых ком­

характеристик термостойкого материала

позитов // Механика композитных материалов.

 

на основе углеткани от температуры:

 

1985,№ 5. С. 816-825.

 

 

l-Ei(T); 2-Е2(Т); З-а^СГ); 4-G^(T);

S-UjiT);

 

10.

Каламкаров А. Л., Кудрявцев

Б. А.,

 

 

 

6-G~2(T)

 

 

 

 

Партон В. 3. Асимптотический метод осред­

 

 

 

 

 

 

 

нения в механике композитов регулярной струк­

 

 

 

 

 

 

 

 

туры // Итоги науки и техники. ВИНИТИ. Сер.

В

табл.

5.2.1

приведены

харакгеристики

Механика деформируемого твердого тела. 1987.

конкретных

композиционных

материагюв. В

N^ 19. С. 78-147.

 

 

 

И. Кильчинский А. А. Об одной модели для

связи с тем, что номенклатура этих

материалов

 

определения термоупругих характеристик

мате­

непрерывно расширяется, в табл. 5.2.2 представ­

риалов, армированных волокнами // Прикладная

лены средние харакгеристики

основных

классов

механика. 1965. Т. 1. №

12. С. 65.

 

отечественных композитов. Индекс 12 соот­

 

 

12.

Кристенсен Р. Введение в механику

ветствует характеристикам в плоскости слоя; Vj2

композитов. М.: Мир, 1982. 334 с.

 

- коэффициент Пуассона при нагружении в на­

 

13.

Максимов Р. Д., Плуме Э. 3., Понома-

правлении армирования; (г^2 ^

'^12

"

^ОДУЛЬ

ров в. М. Прочностные свойства однонаправ-

сдвига и предел прочности при сдвиге. Из пара-

ленно армированных гибридных композитов //

с п и с о к ЛИТЕРАТУРЫ

315

Механика композитных материалов. 1984, N 1.

С.35-41.

14.Москвитвн В. В. Циклические нагружения элементов конструкций. М.: Наука, 1981. 344 с.

15.Победря Б. Е. Механика композицион­ ных материалов. М.: Изд-во МГУ, 1984. 336 с.

16.Рабинович А. Л., Верховский И. А. Об упругих постоянных ориентированных стекло­ пластиков // Инж. журнал. 1964. N1 1. С. 90-100.

17.Радиньш И. Г., Гурвич М. Р. Вязкоупругие свойства слоистых армированных пласти­ ков при длительном плоском нагружении // Механика армированных пластиков. Рига: Риж­ ский политех, ин-т, 1981. С. 76-89.

18.Роулацдс Р. Течение и потеря несущей способности композитов в условиях двухосного напряженного состояния: Сопоставление расчета

иэкспериментальных данных // Сб. переводов. Механика, 16. М.: Мир, 1978. С. 140-179.

19.Скудра А. М., Булаве Ф. Я. Структур нал теория армированных пластиков. Рига: Зинатне, 1978. 192 с.

20.Скудра А. М., Булаве Ф. Я. Прочность армированных пластиков. М.: Химия, 1982. 216с.

21. Скудра А. М., Булаве Ф. Я., Гур­ вич М. Р., Круклиш»ш А. А. Элементы строи­ тельной механики стержневых систем из компо­ зитных материалов. Рига: Зинатне, 1989. 250 с.

22.Хашин 3., Розен В. Упругие модули ма­ териалов, армированных волокнами // Приклад­ ная механика, 1964, № 2. С. 71-82.

23.Хилл Р. Теория механических свойств волокнистых композиционных материалов. Ч. 1. Упругое поведение // Механика. Вып. 2. М.: Мир, 1966. С. 131-143.

24.Хорошун Л. П. О методике определения упругих модулей армир<5ванных тел // Механика полимеров. 1968. № 1. С. 78.

25.Цай С , Хан X, Анализ разрушения композитов // Механика. Вьш. 16. М.: Мир, 1978. С. 104-139.

26.Шаффер Б. Соотношения между на­ пряжениями и деформациями для армированных 1шастиков при действии внешних сил парал­ лельно и нормально их внутренним волокнам //

Ракетная техника и космонавтика. 1964. № 2. С. 163-169.

27.Щермергор Т. Д. Теория упругости микронеоднородных сред. М*.: Наука, 1977. 399с.

28.Вах 1. Déformation behaviour and failure of glassfibrereinforced resin material // Plastics and Polymers. 1970. Vol. 38. N 133. P. 27-30.

29.Chou T. K., Fukuda H. Stiffness and strength of hybrid composites // Composite Materials Japan. - US conf. Tokyo. 1981. P. 87.

30.Foye R. L. An Evaluation on Varions Engineering Estimâtes of Transverse Properties of

Unidirectional Composites // SAMPE Journal. 1966. vol. 10. P. 16-31.

31.Gruber M. В., Overbeeke I. L., Chou T.

A reasonable sandwich beam concept for composite compression test // J. Composite Materials. 1982. Vol. 16. P. 162-171.

32.НШ R. llieory of Mechanical Properties of Fibre Strengthened Materials. 111. Self-consistent Models // J. Mech. Phys. SoUds. 1965. 13. P. 189198.

33.Kalnin I. L. Evaluation of unidirectional glass-graphite fiber epoxy resin composites // Com­ posite Materials: Testing and Design (Second confér-

ence), ASTM STP 497, American Society for Testing and Materials, 1972. P. 551-563.

34.Knappe W., Schneider W. Bruchkriterien fur unudirektionelen Glasfaser / Kunststoffe unter ebenen Kurzzeit und Langzeit Beansprunchung // Kunststoffe. 1972. Bd. 62. H. 12. S. 864.

35.Pinselli R. La fibre aramide et ses appUcation dans les composites hybrides / / Matériaux et techniques. 1984. Jan. Fev. 72 Année, NS 1-2. P. 43-48.

36.Short D., Summerscales I. A new theoretical approach to the strength of fibre composite hybrid materials // 13th Reinforced Plastic Congress, London, 1982. P. 225-229.

37.Whithney J. M. Elastic ModuU of Unidirectional Composites with Anisotropic Filaments // J. Composite Materials. 1967. Vol. 1. P. 188-193.

38.Wu E. M., Ruhmann D. C. Stress rupture of glass-epoxy composites: environmetal and stress effects // Composite Rehability. ASTM STP 580. 1975. P. 263-287.

316

Р а з д е л 6 КОЛЕБАНИЯ МЕХАНИЧЕСКИХ СИСТЕМ

Глава 6.1

КОЛЕБАНИЯ ЛИНЕЙНЫХ СИСТЕМ С КОНЕЧНЫМ ЧИСЛОМ СТЕПЕНЕЙ СВОБОДЫ

6.1.1. ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ

.ЛИНЕЙНЫХ СИСТЕМ С КОНЕЧНЫМ ЧИСЛОМ СТЕПЕНЕЙ СВОБОДЫ

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

дТ

дТ

дФ

дП ^^,.

^^

дд,

 

+ _ +

= Q.(t)

(/ = 1,2,...,4

dq^

dq^

dq^

 

 

 

 

(6.1.1)

Линейными называют

системы,

 

колебания кото­

рых описываются линейными

дифференциаль­

ными уравнениями, что имеет место, в частно­ сти, если кинетическая энергия системы Г, диссипативная функция Рэлея Ф и потенциальная энергия П могут быть представлены в виде квад­ ратичных форм обобщенных координат qi и ско­ ростей q^\

-п п

Потенциальная энергия системы приводит­ ся к форме (6.1.2), если обобщенные координа­ ты qi отсчитываются от какого-либо из положе­ ний равновесия системы, в которьгх ôlJ/dq^ = О

(/—I, 2,...,«). Система соверщает колебательные движения около положения равновесия, если оно является устойчивым. Достаточным услови­ ем устойчивости является минимум потенциаль­ ной энергии и, как следствие, положительная определенность квадратичной формы (6.1.2) потенциальной энергии, которая может быть установлена критерием Сильвестра:

^11

^12

:

^\т

^21

^22

 

^2т >0 (т = 1,2,...,/î).

^ml

^т2

С

тт

Пример 1. Определить, при каких значени­ ях жесткости пружины положение равновесия в поле тяжести системы, представленной на рис. 6.1.1, будет устойчивым.

t-J = 1

содержащих постоянные коэффициенты Л/,-, Ь^р Су, называемые соответственно обобщенными инерционными коэффициентами, коэффициен­ тами вязкого сопротивления и квазиупругими коэффициентами.

Часто для получения линейных дифференциальньгх уравнений дополнительно вводят по­ нятие малых колебании, при которых обобщен­ ные координаты qi и скорости q^ рассматривают как величины первого порядка малости, а в функциях Т, Ф и П ограничиваются малыми второго порядка. Таким образом, некоторые системы при малых колеба}{иях ведут себя, как, линейные.

Рис. 6.1.1. к определению устойчивости положения равновесия

 

 

 

 

 

ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ЛИНЕЙНЫХ СИСТЕМ

 

 

 

 

 

 

 

 

317

 

Потенциальная

энергия

системы

 

выражает­

тивная функция, вычисляемая как половина

мощ­

ся формулой

 

 

 

 

 

 

 

 

 

 

 

 

ности,

развиваемой

 

силами

 

сопротивления

П = 3mgl(l

- со8ф) - 2mgl{\

-

cosvj/)

 

+

(вязкого трения),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

O-HjJCj

/ 2 + 1 ^ 2 ( ^ : 2 - Х | ) ^ / 2

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ cl ( ф - х | / )

/ 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая,

 

что

 

при

 

 

а « 1

Потенциальная энергия системы равна энергии

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c o s a

« 1 -

а

/ 2 ,

запишем

/7(н/,

ф)

в

форме

упругих деформаций пружин

 

 

 

 

 

 

 

 

(6Л.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

77 = CjXj

/

2 + С2 (^2

-

Xj )

/

2

+

 

 

 

Я =

с^iM/^ +2с12Ч/ф+С22Ф^1/2,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где cii=cP-2mgl;

ci2'=C2\=-cP;

C22—cP+3mgL

 

 

+ С3 (Х3 - Х2 )

/

2 + С4Х3

/

2.

 

 

быть

Согласно

критерию

Сильвестра

должно

 

Учитывая соотношения

между

mj^^ Цу^ и Су^

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и сравнивая

полученные

выражения

Г, Ф и 77 с

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>0;

и

42 >0,

 

 

 

 

(6.1.2) и (6.1.4),

находим

 

 

 

 

 

 

 

 

 

 

 

 

41

 

 

 

 

 

Ч^

 

О

 

 

О '

 

 

 

 

 

 

 

 

'21

^22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О

4w / 3

 

О

 

 

 

 

 

откуда следует, что положение равновесия v|/=0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(р=0 является устойчивым, если

c>6mg/L

 

 

 

V^3y

 

 

О

 

О

 

 

w^

 

 

 

 

 

 

Дифференциальные

уравнения

колебаний

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

линейной

системы около положения

равновесия

 

 

2|я

-\х

 

О ^

 

 

Г2с

 

- с

 

О ^

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

(6ЛЛ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6Л.2) имеют вид

 

 

 

 

 

 

 

 

 

 

В

 

- ц

2|л

 

,

С -

 

с

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О

- ц

 

 

 

 

 

О

 

 

 

 

 

7=1

 

 

 

 

 

 

 

 

 

 

 

 

(6Л.З)

Записывая работу 5 ^ возмущающей силы

F{î) на

или в матричной

форме

 

 

 

 

 

 

 

возможньЕх

перемещениях

ÔXj,

8x2,

6x3

грузов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

системы

0^=7^(00^:1,

определяем

компоненты

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aq + Bq + Cq

- Q ( 0 ,

 

 

 

 

(6.1.4)

вектора

 

столбца

 

внешних

 

возмущений

где q=(<7/) - вектор (матрица-столбец)

 

обобщен­

Q(0=(f(0,O,O)T.

 

 

 

 

 

 

 

 

 

 

 

 

ных координат;

А=(ау);

 

В=(Ьу);

C=(c^y)

 

- со­

 

2. Задачу, сформулированную в примере 2,

ответственно

симметричные

матрицы

 

обобщен­

легко решить, используя уравнения

поступатель­

ных инерционных коэффициентов, коэффици­

ного

движения

твердого

тела

W х

= У^^А;/ •

ентов вязкого сопротивления и квазиупругих

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к

 

коэффициентов;

Q(t)

-

вектор обобщенных

сил,

Проекции сил Fj^p

приложенных к у-му

грузу,

обусловленных

возмущающимися

воздействия­

легко

определить с

помощью

рис. 6.1.2,

б.

Для

ми.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример

2.

Составить

дифференциальные

трех грузов

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

WjJCj

= F{t) -

HjXj

-

CjXj

-f 1^2 (^2

-

^1 ) +

уравнения колебаний около положения равнове­

сия системы, состоящей из трех грузов

 

(рис.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6.1.2,

а),

если

их

массы

т\=т,

 

/W2=4m/3,

 

 

"T C'y (X-^

— Xi ^5

 

 

 

 

 

 

 

 

 

 

гпу^т,

сила сопротивления А;-го демпфера

про­

/«2^2

=-M2(-^2

- - i i ) - C 2 ( X 2

- X i )

+

 

 

порциональна

скоростям

движения

 

"поршня"

 

 

относительно "цилиндра" Fj^^^^^=\xj^Vfç^'^^, при­

 

 

Ч-Цз(Хз

- Х 2 ) + С з ( Х з

- X 2 ) ;

 

 

 

 

чем )Л1=Ц2~ЦЗ~Й4~Ц' ^ сила упругости

/:-й

 

 

 

 

 

 

пружины

пропорциональна

ее

 

удлинению

тзХз

= - Цз(хз - Х 2 ) - С з ( х з

- Х 2 ) -

 

 

Fl^y^P=CfçA/^, где С1=С2=Сз=С4=с; Лу^ - деформа­

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ция пружины.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- ^4^3 ~ ^4^3*

 

 

 

 

 

 

 

 

 

 

 

1. Решим

задачу

с

помощью

 

уравнений

 

 

 

 

 

 

 

 

 

 

 

 

Лагранжа

П

рода.

 

Выбираем

 

в

качестве

Сравнивая полученные уравнения с (6.1.3),

обобщенных

координат

 

перемещения

Xj,

Х2,

(6.1.4), приходим к тому же результату,

что

и в

Х3

соответствующих

грузов

из

 

положения

п. 1.

Форма

(6.1.4)

записи

дифференциальных

равновесия.

 

Кинетическая

 

 

 

энергия

 

 

2

 

 

 

2

 

 

-2

 

 

 

 

 

 

уравнений движения, называемая прямой [9, 67],

Т - m|Xj

/ 2 + ^«2^2 /

-^ "^ Щ^ъ

I

^'

 

Диссипа-

не является

единственной.

 

 

 

 

 

 

 

 

 

318 Глава 6.1. КОЛЕБАНИЯ ЛИНЕЙНЫХ СИСТЕМ С КОНЕЧНЫМ ЧИСЛОМ СТЕПЕНЕЙ СВОБОДЫ

^WV. N

>^^ччччччччччч\ччччччччччччч\\ччччччч\ч\чччччччч;ч

а)

/?;=^^/

 

Rif.=Z^ Х^

— ^ -^

 

 

Fz^c(xM

Fj-C(XfX2)

Р^=~СХз

— ^ -^—Î I

^Ч\ЧЧ\ЧЧЧЧЧЧ\ЧЧЧЧЧЧ\\\Ч\\ЧЧЧЧЧЧЧ\\ЧЧЧЧ\\\ЧЧЧ\Ч^Ч lit

ю

Рис. 6.1,2. Система с тремя степенями свободы

КОЛЕБАНИЯ ЛИЫЕЙНОЙ СИСТЕМЫ С ОДНОЙ СТЕПЕНЬЮ СВОБОДЫ

319

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

п

<ii =Ya^ij^j

0' = 1Д-..,/^),

(6.1.5)

7=1

где ô^y - приращение /-Й обобщенной координа­ ты Qi от действия единичной у-й обобщенной силы Gy~l- Симметрическую матрицу Ъ={Ъф называют матрицей коэффициентов влияния Ъу.

В обобщенных силах, кроме возмущающих

В качестве обобщенных координат примем прогиб балки в месте крепления груза у и угол поворота фуза ф. Приложив единичные силу и момент по правилу Верещагина [86], вычислим коэффициенты влияния оц—4Р/9^/,

012—021=2/2/9£7, b^T^ll^EJ. Тогда ма1рицы, входящие в уравнение (6.1.4), будут иметь вид

(уЛ

q =

в

10,

 

О,

 

 

 

О

\

1

11

Q ( 0 -

 

ь =

 

 

 

9£7

2 / ^ 3/

Qj{t) и сил сопротивления среды (внещнее тре­ ние) -ЬцС:, учитывают даламберовы силы инер­

ции -djjÇj:

 

Qj=Qj{t)-bjjqj-ajjqj

( у = 1 , 2 , . . . , 4

Подстановка Qj в (6.1.5) приводит к уравнениям колебаний в обратной форме, которые в матрич­ ной записи имеют вид

q = 5 ( 0 ( 0 - Bq - Aq).

(6-1.6)

Матрицы \—{аф ' и B=(^^y) коэффициентов инерции и сопротивления в обратной форме записи являются диагональными.

Для перехода от уравнений (6.1.6) в обрат­ ной форме к прямой форме (6.1.4) необходимо (6.1.6) домножить слева на матрицу 5"^=С квазиупругих коэффициентов.

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

6.1.3).

/77. /

С = 5"

3 / 8 / 3

- I / 4 /

Л I М\

\/21

 

6.1.2.КОЛЕБАНИЯ ЛИНЕЙНОЙ СИСТЕМЫ

СОДНОЙ СТЕПЕНЬЮ СВОБОДЫ

Всоответствии с (6.1.3) дифференциальное уравнение колебаний имеет вид

aq + bq-\-cq=:Q(t).

(6.1.7)

Свободные (собственные) колебания

происходят

при отсутствии

возмущающих

воздействий

[G(/)~0]> если в начальный момент

(t=0) сис­

тема выведена из состояния покоя.

 

Свободные колебания системы без трения

(консервативная

система). Дифференциальное

уравнение колебаний согласно (6.1.7) имеет вид

aq + cq = О,

общее

рещение

которого

q(t) = As,m(cDt+ а) (рис. 6.1.4, а),

где круговая

(циклическая) частота со=(с/т)^'^,

а амплитуда

А и начальная фаза а

колебаний

зависят

от на­

чальных

условий q(0)=q(),

q(0)

-

q^,

причем

л =\qo

+(qo /

(X)) j

; a

- arctg(^oco

/ ^Q).

Периодом

колебании

Т=2%/а) называют

продолжительность одного полного цикла коле­ баний (рис. 6.1.4). Используется также частота в герцах v=l/7'==(o/27C, равная численно числу колебаний в секунду.

Свободные колебания системы при наличии вязкого трения (диссипативная система). Диффе­

ренциальное уравнение колебаний

согласно

(6.1.7) имеет вид aq + bq -^ cq = О или

 

2

(6.1.8)

q + 2zq + С0 q = О,

где величину г=Ь/2а называют коэффициентом демпфирования.

 

При '"MajiOM" демпфировании (s<oc))

реше­

шш. ние (6.1.8)может бьггь записано в виде

 

 

 

q(t) = Ае"''^ sm((D^t + а),

(6.1.9)

 

где (Oi=(ci)2-s2)^'^, а /4 и а определяются

из

на­

Рис. 6.1.3. Система с двумя степенями свободы

чальных условий