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

книги из ГПНТБ / Давыдов, В. В. Динамические расчеты прочности судовых конструкций учебник

.pdf
Скачиваний:
20
Добавлен:
22.10.2023
Размер:
12.6 Mб
Скачать

Подставив (3.41), в котором неизвестными являются только пара­ метры а.;, в уравнение (3.40), получим

00

2 а! [(1 + «О EIf)v (х) -f Rmfj(x)— ma)2fi (х)] еш = р (х) еш .

Заменив четвертую производную f)v (х) самой функцией по урав­ нению (3.6) Е Ifjv(x) = mtijfj ),

ОО

 

 

получим еш 2 А/ [(1 +

xt) тК2— та2-f- Ria] f,■(х) = р (х) еш

 

или

tfmfj(x) = p(x) ,

(а)

если ввести обозначение

 

 

 

2

(3-42)

Для того чтобы избавиться от ряда и найти параметры а-п умножим ряд (а) на любую форму колебаний /) (х) dx и проинтегрируем его по длине балки. В силу первого из условий ортогональности (3.20) все интегралы, кроме одного, в котором перемножаются одинаковые формы колебаний (одного тона), исчезнут * и мы получим

„ ч 2 ,

со2 . . / 2гш

J mf2(х) dx — Pfj (а).

(б)

Поясним получение в уравнении (б) правой части. Приложенную в сечении х= а сосредоточенную силу Реш можно рассматривать как распределенную на небольшом участке вблизи х= а нагрузку боль­ шой интенсивности (см. рис. 26). Поэтому после умножения (а) на /у (х) интегрирование по длине балки сведется к интегрированию на участке Ах

I

а+Д*

J Р (*) // (х) dx=

| pfj (х) dx = pfj (а) Ах = Pf, (а).

О

а

Обозначив в соответствии с (3.27) и (3.22) через Nj обобщенную

жесткость

 

Nj = X21 т [fj (х)]2 dx

(3.43)

о

 

* Полезно напомнить, что ряд (а) содержит все функции / х,

f 2, f 3, . . . ,

[j, . . . . а умножаем мы его на одну определенную |у.

 

125

из (б), получим

 

 

 

 

 

 

 

 

P f j (а)

1

аГ

 

 

P fi (а)

 

 

1

1

 

----------1 !

 

 

 

 

 

 

 

 

 

 

! 1

1

^ | з

+

^ 1

+

1 .

Гй

7

 

N /

 

 

 

 

 

 

 

 

 

 

 

Подставив а;-

в (3.41)

и заменив в ней

 

 

 

 

 

 

 

еш cos соt +

i sin соt,

 

найдем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2rw + X

w =

 

Р fj (а) fj (x)

 

 

 

 

*7

 

 

 

 

 

 

 

 

Л -

i *

 

 

 

 

 

 

 

 

1

 

ч

 

 

з

 

3

 

+

1

« |

1

I

 

 

Г 2

 

—■, М(

 

 

1 '

 

+

>> is? 1е

+

S

IS 1---------- 1

(cos соt + i sin coO

+ X

или

°°

 

 

 

I^ 1---- — ^ cos соt +

2rm

x ) sin wt + i [

 

 

 

 

T

 

W =

P fi (а) fj (х)

'

 

 

 

 

 

J m A

 

 

Nj

1 1 - i *

 

 

 

 

H i

 

 

 

 

 

 

 

 

 

 

{

 

 

*7

 

 

*7

 

 

 

 

 

 

 

 

 

 

Отбросив в

числителях членов

ряда

мнимые слагаемые ( [ . . . ] ,

в соответствии

с вещественными

значениями возмущающей

силы,

Р cos с о и

объединив косинус и синус в одну функцию (см. примеча­

ния к с. 20),

получим

 

 

 

 

 

 

 

 

w(x,

t) =

Pfj (a) f; (х) cos (соt — р;)

(3.44)

^

 

 

со2 \

2 ,

/ 2гсо

 

 

 

 

 

 

 

 

 

 

3

 

 

>2

 

 

*7

 

 

 

 

 

 

 

 

 

 

 

 

 

5 i=1

 

 

 

 

 

 

где

 

 

 

2лсо

 

 

 

 

 

 

 

 

 

sin Ру

 

*7

 

 

 

 

 

 

 

 

 

 

/ ( ‘- f N\

7f -

(3.45)

 

 

 

 

 

1 - ^ 1

 

 

 

 

 

 

 

COS Р ; =

___ ±

 

 

 

 

 

 

 

/e-t

2

+ I —

 

-(- К

 

 

 

 

 

 

 

 

 

 

 

 

 

I

 

 

,

 

При отсутствии сопротивления формула (3.44) переходит в (3.28),

а р;- = 0.

126

Если бы колебания вызывались возмущающим моментом М соsat, приложенным в сечении х = а, то вместо (3.44) мы получили бы

Wj (х,

(3.46)

при тех же значениях фазовых узлов (3.45).

Решение в замкнутой форме с учетом сопротивления. Интеграл

уравнения (3.40) при правой части, равной нулю,

можно искать и

в замкнутой комплексной

форме в

виде некоторой

функции формы

F (х), умноженной на еш

= cos at

+ i sin соt:

 

 

w{x, t) = F(x)eiat

(3.47)

или ввести две функции формы

w (х, t) = Fx(х) cos at + F2 (x) sin at.

Взяв решение в последнем виде, подставив его в (3.40) и приравняв нулю коэффициенты при гармонических функциях cos at и sin at, получим два дифференциальных уравнения четвертого порядка, в ко­ торые войдут обе функции Fx (х) и F2 (х) и их четвертые производные. Путем некоторых преобразований можно одну из функций исключить, приведя задачу к интегрированию уравнения восьмого порядка, во­ семь постоянных которого определяют из граничных условий и усло­ вий сопряжения. Процесс разыскания восьми постоянных довольно громоздок.

Подстановка решения (3.47) в уравнение (3.40) после сокращения на неравный нулю множитель eiat приводит к уравнению вида

(1 + XI) EIFlv (х) — та2 [1 — — ) F (х) — 0.

\

пиа /

Процесс выделения вещественной части после интегрирования этого уравнения получается также достаточно сложным. Поэтому мы ограничимся детальным рассмотрением лишь случая, когда суще­ ствует один вид сопротивления, а именно внутреннее (гистерезисное) сопротивление, внешнее же отсутствует. В этом случае уравнение (3.40) оказывается несколько проще

(1 -f- xt) EI

d*w

т d2w

0.

(3.48)

 

дх4

~dt*

 

 

Подставив (3.47) в (3.48) и сокращая на неравный нулю множитель ё~ш*, получим уравнение относительно неизвестной функции формы

F™(x)

тш‘

F(x) = 0.

(3.49)

(1 +

х/) El

 

 

 

127

Общий интеграл этого дифференциального уравнения

F(x) — A sin v-y- + Bcos v-y- + Csh v ^ - - f D chv-y-

(3.50)

отличается от (3.31) с. 117 лишь иным значением безразмерного ар гумента

mcoi

(1 -|-х < )£ / -j/ j _j_ xt

где v — определяется формулой (3.32).

Постоянные А, В, С, D находятся из граничных условий, а при действии на балку сосредоточенных усилий в пролете еще и из усло­ вий сопряжения.

Комплексный аргумент v можно преобразовать следующим обра­ зом:

v = v ( 14->«) = V [ 1 ------ X,i----- — и а 4 - . . .

\4 32

Учитывая относительную малость коэффициентов неупругого со­ противления,* можно ограничиться удержанием лишь двух членов, т. е. принять

v = v 11---- — ш ].

(3.51)

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

гумент V.

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

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

Например:

. -

. /

. XV \

.

.XV

.

. XV

sin V =

s in

V— l

= s in V COS I ---------

COS V s in

t -------

 

\

4 )

 

4

 

4

sinv = sinvch ----------

1 COS V sn

----- .

4

 

4

Даже для днищевых перекрытий ф < 0,6 и согласно формуле (1.6) х =

* 0, 1.

128

Таким образом получаются и остальные формулы перехода

 

sin v = sinvch

XV

i cos v sh

XV

 

 

4

 

 

4

 

cosv = cosvch

xv

i sin v sh

XV

 

 

4

 

 

4

(3.52)

 

XV

 

 

XV

sh v = sh v cos

-i

ch v sin

 

 

4

 

 

4

 

ch v = ch v cos

XV

-i

sh v sin

XV

 

 

4

 

 

4

 

Найденную функцию F (x)

надо

подставить в формулу

(3.47)

w (х, t) = F (x) еш = F (x) (cos cat + i sin соt)

(3.53)

иосвободиться от мнимых членов (отбросить их).

Пр и м е р . Найти вынужденные колебания свободно опертой по концам призматической балки под действием сосредоточенной силы Р cos at, приложен­

ной к середине ее пролета (см рис. 25 при а = Ь) со следующими данными:

Длина балки

/,

м ..................................

 

 

 

 

10

 

 

Погонная масса т, т/м .......................

 

 

 

 

0,0086

[8 ,75-10 5* кгс-с2/см2]

Момент инерции сечения /,

м4

. .

.

. 180-10

8[180 см4]

Модуль упругости Е, кН/м2

...............

 

 

1,96-108

[2-10е кгс/см2]

Амплитуда силы Р, кН .......................

 

 

 

 

0,225

[23

кгс]

Круговая частота <о, с

 

 

 

 

18

 

 

Коэффициент

неупругого

сопротивле­

 

 

ния .........................................................

 

 

 

 

 

 

0 , 2

 

 

Решение в форме ряда.

Формы главных

колебаний рассматриваемой балки

(см. случай 5 табл. 3) являются синусоидами

 

 

 

 

. . .

.

х

=

.

/лх

 

 

fj (х)

= sin p.j

sin

J - j - ,

 

 

в связи с чем в ряде (3.44)

 

 

 

 

 

 

 

 

 

w =

 

P f l ( a ) h

(х) cos (at — Р/)

(в)

 

 

 

 

1—

+ X2

 

2

j

 

 

 

 

 

 

 

 

 

 

 

i=i, з, 5, . . .

останутся лишь слагаемые с нечетными /, все четные слагаемые равны нулю,

так

как а = —

и f ; (а)

= sin —

=

0 при

четных /'.

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

равны

 

 

Собственные частоты колебаний балки* по формуле (3.8)

 

к■. = $ _ л / ~ Е /

Лгс2

j

/ ~ Ё 7

_

/2я 2

f

1 , 9

6 - 1 0

^ - 1

8 80 -

1 Q — .

7

/ 2 V

т

[2

у

т

 

Ю2

V

 

0,0086

 

1 '

 

 

 

/

Г у'2 я 2

 

2• 10е -1180

= 20/2

 

 

 

 

 

 

 

L10002

 

8,75по-1C-5

 

 

 

 

*Формы и

частоты

первых трех тонов этой

балки

приведены

на рис. 23.

5

В. В. Давыдов,

Н. В. Маттес

 

 

 

 

 

 

 

 

129

Входящие в расчетную формулу обобщенные жесткости по формуле (3.23)

равны

I

 

 

I

 

 

 

 

 

 

 

 

 

Ni = EI (

[f"j(x)]2dx = E f

I

( 'L ^ s i n i y - j

dx =

 

n*EIj*__

97,4 • 1,96 ■108 • 180• 10“ 8/ 4

, 2

/ k4 H m/ ;

 

2/3

2 - 103

 

 

1 7

 

 

 

 

 

 

 

,,

97,4-2-10e-180-/4

 

 

.

 

 

Ni =

-------------------------------

= 1 7 ,5 i4 КГС/СМ

 

 

 

2 -1000s

 

 

 

 

 

Вычисление

применительно к приведенной

выше

расчетной формуле (в)

и формулам (3.45) для фазового угла

 

 

 

 

 

sin р,- =

;

cos Р/ =

 

ч

 

 

 

 

 

Vi1 —-»!

выполнены в табл. 1 0 .

Заимствуя амплитуды колебаний и фазовый угол из граф 8 и 11 таблицы, получим решение задачи

 

 

 

w (х , /)мм =

 

 

ггу

 

 

 

 

 

 

 

47,41 sin -у- cos (о)/ — 46°30') —

 

 

— 0,16 sin

cos (о / — ll°30') +

0 ,0 2 sin ^ p co s(c o / — 11°30') + . . .

Амплитуды

колебаний середины балки (под силой)

могут

быть

найдены

следующим

,

 

 

/

имеем

 

 

 

 

 

образом:

при х = - у

 

 

 

 

 

w ( — , Л

= 47,41 cos (со/ — pt) +

0,16 cos (со/— рз) + 0,02 cos (со/ — р5) + • • ■;

' 2

'ММ

 

 

 

 

 

 

 

 

 

 

 

w (— ,

/'j

= [ 4

7 , 4 cos1

Рх + 0

, 1

cos6

рз + 0 , 0 2cos р5 + . . .] cos со/+

'

2

'мм

 

 

 

 

. . .] sin со/ = 32,8 cos со/ + 34,4 sin со/.

+ [47,41 sin рх +

0,16 sin рз + 0,02 sin рб +

Амплитуды при х

вычислены в графе 12, а выражения

в квадратных

скобках

в графах 13

и 14 табл. 10.

Результирующая

амплитуда,

очевидно,

равна

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а макс =

32,82 + 34,42 = 47,5

мм.

 

 

 

Ввиду близости к резонансу с первым тоном (со =

18,

a A,j =

20 с-1 ) сопро­

тивление влияет довольно сильно. Расчет без сопротивления дал бы амплитуду

69 мм.

Вдали от резонанса даже довольно большое внутреннее сопротивление (х = 0,2) почти не сказывается. Расчет той же балки под действием возмущаю

щей силы той же амнлитуды (0,225 кН), но с частотой со — 30 с- 1 дал бы ампли туды колебаний под силой 10,2 мм с учетом сопротивления и 10.5 мм без учета сопротивления.

Решение в замкнутой форме. По формуле (3.37) при х = — при дейст-

,

/

получим при отсутствии

вии силы в середине пролета, т. е. при а = Ь =

 

130

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

 

 

 

 

О

 

 

 

 

 

 

 

 

 

а

 

 

 

 

 

 

 

 

 

а-

 

 

 

 

Р/3

 

 

 

а

 

 

 

 

 

 

 

ч

 

 

 

 

 

 

 

 

 

VO

 

 

 

 

' 2EIv3

sin V

а

 

 

 

 

Е~>

 

 

sh2

-

 

 

 

 

 

 

 

 

sh v

I COS COt

=

 

 

 

 

 

 

 

 

 

Р13

 

/,

v

 

v

\

cos сot.

 

----------

 

tg ------- th —

/

 

4£/v3

 

\

2

 

2

 

 

Аргумент v равен по формуле

 

(3.32)

 

 

 

 

 

 

 

 

 

 

 

 

,

/

ягсо2

 

 

 

«d

 

' -

'

V

т т =

 

 

К

 

 

 

0S

 

 

 

 

 

 

 

 

 

о.

= ш

V

 

 

0,0086- !8 2

 

_

«

 

 

 

Е

 

 

1,96-108-180-10—8 _

О .

 

 

о

 

 

 

 

= 2,98.

 

 

 

’в*

 

 

 

 

 

 

 

 

Подставив значения тригоно­

ж

метрических

и

гиперболических

ч

cd

функций,

 

получим:

 

 

 

ю

 

 

 

 

 

w --- ■

0,225-103

(12,350 —

«d

 

 

 

 

 

4-1,96-180-2,983

'

 

 

о

 

 

 

 

 

 

 

 

 

VO

 

 

 

 

 

 

 

 

 

4

— 0,903) cos cot 0,069 cos соt,

о

ж

т. e. амплитуда

колебаний сере­

м

3

дины балки равна 69 мм.

 

 

X

 

 

X

Для

 

учета

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

4)

 

ЕС

нужно в то же

выражение (3.37)

к

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

>>

5

v по формуле (3.51)

 

 

 

2

 

 

 

а

 

 

 

 

/

 

 

 

 

ж

 

 

 

 

 

 

 

 

 

н

 

 

 

 

 

 

 

 

 

а>

 

v =

v

1 ----- -- v.i

 

 

sr

 

 

 

о

 

 

 

 

\

4

 

 

 

 

Cl

w =

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

ww'fluis [si]

ки S03 [31]

ww _ uis [д]

1г/

шс,

Й= ' Й50Э

=$ UIS

.А! N

(о) 'id

г\Усст / /А

г1

./г‘л

'jv

' егг

 

1гU!S “ (0>^

rf О О ^ o' o'

—• СО OJ

т}4 »—■ о h- о о

О О О О

СО СО СО СО

о о о о

СО —1— —

h- О © О

00 00 00 оо

СО© © ©

юоо со «о

© © ©

о"© о"o'

со© © О ^ —CN <МО О О

©© © ©

©© © ©

—« © © ©

сосо

© © — ©

— — © ©

оо © © ©

о о о ©

(МCN) —*Ю

оо со сч о

© ~ © ©

© о о ©

СЧ00 О 00

—<ю©

5*

131

 

У ч и т ы в а я , что

 

 

 

 

 

 

 

 

 

 

з

= 1 Н------хг,

 

 

 

1

------- хг

 

4

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

получим

 

 

 

 

 

 

 

Р1з j 1 +

щ

v

 

xv .

, I V

XV .

sin | -------------г

sh —

------- г

4EIv3

 

2

 

8

2

8

егшг

 

V

 

XV

ch | —-----—

i

 

 

cos

 

 

 

 

 

 

 

2

8

 

Далее следует воспользоваться формулами (3.52) и избавиться от мнимости в знаменателях. После довольно длительных преобразований получим

PI3 j 1 + - j - хг ]

8EIV3

, . . XV sh v — г sin —

, ,

V

,

„XV

sh2

------h cos2 —

 

2

 

8

sin v

i sh

cos*2

• +

sh2 —

\

 

 

I (cos at +

i sin at).

Выполнив перемножения и оставляя лишь действительную часть, оконча­ тельно найдем

 

 

 

 

 

 

 

 

3

, x v

 

 

 

 

 

 

 

 

sin v -4------х sh —

 

 

 

2

/

 

8Д/v3

 

 

4

4

 

 

 

 

COS'1

- +

sh2 XV

 

 

 

 

 

.

 

, 3

.X V

 

 

 

 

 

 

 

sh v H------x sin —

 

 

 

 

 

 

 

________4_______ 4_

I cos at +

 

 

 

 

 

 

, „

V .

„XV

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sh2 ------ b cos2 —

 

 

 

 

 

 

 

 

 

 

2

8

 

 

 

 

 

^

.

 

,

,

xv

3

,

 

. xv \

 

 

------- x sin v +

 

s h —

— x sh v — sin —

 

 

4

 

 

 

 

4

4

 

 

4

sin

at

v .

, „ xv

 

, ,

v

,

„XV

 

 

 

cos2

------ h sh2 —

 

sh2 ------f- cos2 —

 

 

 

2

 

 

8

 

 

2

 

8

 

 

Подставив численные значения входящих в последнюю формулу величин

* Значения тригонометрических и гиперболических функций для аргумен-

v =

2,98;

V

=

1,49;

vv

 

 

VV

 

 

— = 0,149 и — = 0,0745

 

 

2

 

 

4

 

 

8

 

 

заимствуем из «Пятизначных

математических

таблиц» Б. И.

С е г а л а и

К. А. С е м е н д я е в а ,

 

М.,

Физматгиз,

1959.

 

 

 

132

и выполнив вычисления во избежание потери точности с завышенным числом значащих цифр, получим

w А = 0,0402 cos со/ + 0,0320 sin <в/ м.

Амплитуда колебаний середины балки равна

^0,04022-f 0,03202 =0,0514 м или 51 мм.

§ 12

ВЛИЯНИЕ УПРУГОГО ОСНОВАНИЯ и ОСЕВЫХ с и л НА СВОБОДНЫЕ ПОПЕРЕЧНЫЕ КОЛЕБАНИЯ БАЛКИ

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

EI dxi

d2w

■kw=*q.

dx2

Этим выражением можно воспользоваться для составления диффе­ ренциального уравнения свободных колебаний балки, необходимо лишь учесть, что в случае колебаний прогиб w (х , t) является функ­ цией двух переменных х и /, а под интенсивностью загрузки ее пони­ мают загрузку силами инерции

d2w

v dt2

Тогда уравнение упругих колебаний будет иметь вид

Е1 d*w

d2w . -kw... +, m д2т

= 0.

(3.54)

дхл

дх2

3/2

 

 

Решение этого дифференциального уравнения ищем в виде беско­

нечного ряда

ОО

 

 

 

 

(х) cos (Xjt -f a.j).

(3.55)

w (х, 0 = 2 Л

 

/=i

 

 

 

Подстановка его в (3.54) приводит * к дифференциальному урав­

нению относительно функций формы колебаний:

 

 

f7 w -

j t f"i{х) ~

n^ r r ^ f>{x)= °-

(3-56)

Общий интеграл последнего уравнения имеет вид:

 

fl (х) = aj sin г)/ -j- +

bj cos г]/ - у + cysh £y - у +

d,-ch - у ,

(3.57)

где a;-, bjt Cj, dj произвольные

постоянные интегрирования,

а gу

* Порядок рассуждений полностью аналогичен приведенному выше при подстановке (3.5) в (3.2).

133

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

(3.58)

В справедливости проведенного решения легко убедиться, подста­

вив интеграл (3.57) в уравнение (3.56) и заменив в нем

т)у- их зна­

чениями (3.58).

г)у, а также

Для получения постоянных а/, bj, Cj, djt параметров

связанной с ними частоты Яу следует воспользоваться граничными ус­ ловиями на концах балки. Четыре граничных условия дают систему четырех однородных алгебраических уравнений относительно dj, bjt

с,- и dj.

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

Равенство нулю определителя представляет собой транцендентное уравнение относительно неизвестных £/ и т]у, или, что то же, уравне­ ние относительно частоты Яу (от которой £у и ту зависят), т. е. может быть названо уравнением частоты.

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

П р и м е р . Найти собственные частоты поперечных колебаний балки, лежащей на упругом основании жесткости k и свободно опертой по концам на жесткие ножевые опоры, при условии, что балка растягивается силой Т.

Составим, используя выражение (3.57), граничные условия, выбрав начало координат на левой опоре

х = 0 /у = 0

bj + dj — 0;

 

х =

0

/■ = 0

+

= 0 ;

 

х =

l

fj — 0

я/ sin т|/ +

bj cos r\j + Cj sh £y + dj ch £y =

0 ;

x = l

f" = 0

— £ZyT]y sin T)y—&/Пу cos ’I/+ c 3 sh £/+ d $

ch = °-

Умножив первое из этих уравнений на т] 2 и сложив со вторым, получим

(S/+ -П/) rf/ = °- ■

Поскольку £у и т]у — величины положительные, d/ = 0, следовательно, и bj = 0. Из двух оставшихся граничных условий

составим определитель

(l2 4 - T|/) sh sin г), = 0 .

134

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