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

книги / Статические и динамические проблемы теории упругости

..pdf
Скачиваний:
6
Добавлен:
19.11.2023
Размер:
27.66 Mб
Скачать

моментов показаны на рис. 9, а. Представление этих моментов в векторной

форме дано на рис. 9, б. Так как т[ и М\ на одной стороне, а т и М\ на дру­ гой представляют один и тот же изгибающий момент, то равнодействующая

/72/ и М\ должна быть равна и противоположна по направлению равнодей­ ствующей /п/ и М/, как это показано на рис. 9, б. Далее, обозначая через

Ус угол между соседними коленами в точке i и раскладывая моменты на со­ ставляющие по направлениям ///_i и z/_i, получаем

Mi +

Mricos у, — /72/ sin ус =

0;

z

M ri sin уI + /72/ COS yt =

(22)

/72/ +

0.

Еще два уравнения получаются из условий совместности деформаций. Рассмотрим колено, расположенное слева от i-й опоры. Положительные на­ правления всех сил и моментов, действующих на это колено, показаны

на рис. 10. Обозначим через <р{ и ф/ углы поворота поперечного сечения в 2-й опоре, соответственно относительно осей Z/_i и г/7-—\. Для знака углов

моментов /72/_1 и /72/ и

сил

Q/—1 и P z-i, получаем

 

ф/ =

772/OCj

/72/—iOC

P i—l^i

Qi—1^>>

(23)

где aly a2, tx и /2 берутся согласно формулам (3),

(4), (5) и (6). Аналогично,

используя соотношения (7), (8), (16) и (17), находим

 

ф/ =

Af/Pi — TH/_iP2

Т /—is -f- 5/_it22,

(24)

где piT p2, s и v2определяются выражениями (9), (10), (14) и (19). Рассматривая теперь колено, расположенное справа от 2-й опоры и обо­

значая через ф/ и ф/ углы поворота поперечного сечения, расположенного в опоре 2 относительно осей zi и ус этого колена, получаем тем же способом,

что и раньше, следующие зависимости:

Ф; = mj»! — ml+ia2 +

PJi + Qik\

ф£ = MriPi — Mi-L1P2 +

7\-S— SiVL-

Здесь !Tt- = Ti- 1 + 5/_ir представляет собою крутящий момент в t-й опо­ ре, а коэффициент берется из выражения (18). Малые углы поворотов, подсчитанные согласно формулам (23), (24) и (25), могут быть представлены в виде векторов, чьи направления совпадают с направлениями соответствую­ щих углов поворота1 и чья величина пропорциональна соответствующим

углам поворота (рис. 11). Теперь, так как % и ф/, с одной стороны, а ф' и ф/, с другой, представляют повороты одного и того же поперечного сече­ ния опорной шейки в i-й опоре, то результирующий поворот ф* и ф, дол­ жен быть соответственно равен результирующему повороту ср^ и ф/, как это указано на рис. 11. Отсюда получаем

9, = T,cosY ,4 ^ sin Y ,;

i

(Jo)

1|>£= — ф£Sin V,- + l|Jicos у(.

Э т и выражения вместе с зависимостями (22) образуют систему из четы­ рех уравнений для i-й опоры. Аналогичные уравнения можно записать

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

уравне­

ний, равное числу неизвестных величин

т\, т,

 

М\, где i = I, 2, 3, ..., /г,

а п — число опор.

 

 

 

 

(25) ввести в систему уравнений (26),

Если зависимости (23), (24) и

то эта система принимает следующий вид:

 

 

 

 

 

tTLiCL1

ГП{—1СС2

Р i—1^1

Q i—1/2 3=3

 

 

~

 

+ Pih +

Qih) cos у,- +

Mli+iP2 —

+ 7 »

sin y£;

 

Af'px —

 

 

P2 +

Ti—,s +

S£_ ,O2 = (Atfp, — Mi+ip2 —

 

— 5,-t»! +

7 »

cos уi — (/n/aj — m‘i+[a2+

P,./1+

sin yt.

(27)

Теперь выражения типа (22) дают

 

 

 

 

 

 

 

rrii =

-

Mj +

cos V/

 

m£=

 

+ Mri cos y£

 

 

sin V/

 

 

sin Y£

 

 

 

 

 

 

 

 

 

(28)

mt_i =

м Li + Щ-\ cos v,_i

 

mi+ 1 =

 

Mi+ 1 + ^£+1 cos Vj-t-i

 

 

 

 

 

sinV,_,

 

 

 

sinVt+i

 

 

 

 

 

 

 

 

 

 

Подставляя эти величины в (27), получаем

 

 

 

 

 

 

sm T,—! •Mf-i — a, ctg

 

 

— 2ax ctg у,М,- —

 

 

 

 

sin 7,

a xCOS2 7/

Pi sin Y,) Mri

 

 

 

 

 

siny,

+

 

 

£+

 

 

 

 

I

/"ft

*

 

<=08 7 . 0 0 5

7 , + ,

\

 

 

 

+

( p2 sin Yi - CC2 — sinTf+| +

 

JMl+, -

 

 

“ 2

 

 

 

 

 

 

 

+

(ЛЛ + <Ыг) cos Y( -

 

Между положительным направлением оси поворота и направлением вращения су­ ществует такая же связь, как и между вращением и перемещением правого винта.

Sfr sin vt- +

T(s sin Vil

 

M 'i-iP2 -I- К + Pi) M.\ + (ax — Рл) cos y (M rc +

sin Vi cos v i+ l

sin у,-

 

+ p2 cos Vi + a2 sin v i+ i

M'l+l =

Mi-)-i -f- a2 sin Vi+ i

= тi_iS + Tcs cos Vi — S i_,o , — S.Ox cos yt-— (/V i +

ФЛ) sin V.-

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

Особые случаи уравнений ( 2 9 ) . Когда у2 = Уз = V. то уравнения (29) принимают следующую упрощенную форму:

М‘- 1 — a 2 ctg yM 'i-i 2ах ctg уМ.\ —

+ «1 cos v ctg v + Pi sin v) M'i + (P2 sin у

 

 

— a2cosvctgv)M i+i — a 2ctgvMi+i = Л - 1/1 4-

 

+

Qi-\t2+ (Pit,. + Qit2) cos v — Sfl! sin v + 7\s sin v;

(30)

 

Mi—iP2 + (ax + PJ) Mi + (at — Pi) cos уM'i 4-

 

 

+ (Pa + аг) COS yMli+\ + a2Mi+i = — Ti_iS +

 

+

Tts cos v — Si-ittj — Sfl! cos v — (Ptt1 + QitJ sin y.

 

Если коленчатый вал нагружается только кручением, то уравнения (30) принимают вид, приводимый ниже. Подставляя в исходное уравнение Л _ i = = Pt = Qi- 1 = = Sc- 1 = S£ = 0 и полагая Tt-\ = Tf = T, получаем

MLi — a2 ctg yM i-i — 2ax ctg yAl' —

— (lE Y + ai C0SVctgY + Pi sin v) M] + (p3 siny

— a2 cos v ctg v) M'+i — a2 ctg yM'i+1 = Ts sin у;

(31)

— P2^ _ , + (ax Pi) Mi 4" (ai — Pi) cos yMi 4- 4- («2 4- P2) cos vMt+i 4- a^M'+i = Ts (cos v — 1).

ПРИМЕНЕНИЕ УРАВНЕНИЙ МНОГОКОЛЕНЧАТОГО ВАЛА

КСЛУЧАЮ ТРЕХКОЛЕНЧАТОГО ВАЛА

Вслучае трехколенчатого вала имеем

Уг =

Ys = Y = 120°;

sin у =

V^3;

cosy = ----- ctgy = — у = - .

Так

как вал предполагается

свободно

опертым, то изгибающие моменты

в концевых поперечных сечениях равны нулю, т. е.

м{ = M 'I = м{ = м\ = о.

Общая формула для трехколенчатого вала. При допущениях, введенных в предыдущем параграфе, уравнения (30) при двух промежуточных опорах

примут вид

 

 

 

 

 

 

^

(<*1 + Pi) М2 ----

Y (а1— Pi) ^ 2

----

(«2 + Рг)Мз + а2Мз =

= — 7 xs — 4

T2s — V

2 +

4 " s *wi — Щ - (*Vi + СУ2):

otiM -

( 4 «1 +

4 Pi) М2

+ ( 4

р2 -

4 а*) м^ + 4

«2м5=

= 4

Т2S-|--4 — |^1^1 4" Ql^2----

 

о” (Р2^14" QA) I

4“и1^2>

— р2Мо (осх 4 "

Pi) М 3-------

2~ (a i

Р1)

= ^ 2S

2~ ^ з5’

 

— S2U 4" 4 5з^1

4 — (^з^1 4" Фз^г) —

 

агм1+ 4

агМг 4" а1 Мз—

ai Н—4- Pi) М3 =

= 4

Т3s + 4 4

[^Vi +

Ог^г----

 

2 "

")■ Qs4)j ----

4- S3vv

Коленчатый вал подвергается кручению моментами Т, приложенными на концах. Уравнения (32), записанные для двух промежуточных опор, теперь приводятся к виду

К

4 - Pi) Ml -

4

(«1 -

 

Pi) M S -

4

^

+

P2) M i

+

a 2MS =

-------5 - 7 s ;

 

«1M i -

( 4

« 1 4 - 4

Pi) M S 4- ( 4

Рг —

r

« * ) M * +

4

« . M S =

4

7 s ;

 

 

 

 

 

 

 

/

 

1

4

 

 

Px) MS

 

 

з

 

 

 

 

f

(33)

-

psMS 4 - ( « 1 4- Pi) M i

-

 

(a , -

=

-------4

T s;

 

 

 

 

 

-

cc2M i 4- 4

 

a ’-M S +

a i M 3 -

 

( 4

“ i +

4

Pi) м з

= 4

Ts-

 

 

 

 

Числовой пример. Рассмотрим трехколенчатый вал дизельного

двигателя мощностью

900 л. с.,

все размеры которого даны в разобранном

 

ранее случае. Оттуда из зависимостей

(3),

(6),

(9), (10), (14),

(18)

и

(19)

найдем, что

 

 

 

 

 

 

 

 

 

 

 

 

а х =

0,851

2е +

^

:

а г =

0,478

2й +

 

f

 

s =

0,232

2e+ ( .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Si

 

 

 

Вг = 0,864

 

-i- f

; р2 = 0,381 -

2е 4- f

 

и, = 0,408

(2e +

/)r

.

(34)

 

 

 

Di

 

 

V

 

- ■;

 

 

 

 

 

 

 

 

 

 

 

 

D1

 

 

 

 

 

Si

 

 

 

 

 

=

0,640

(2e +

f)r

;

<

=

0,559

(2e4

 

)r

;

f2 =

1,013 (2e +

^ f .

 

 

 

 

 

 

 

£>1

 

 

 

 

 

 

 

w

 

 

 

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Si

 

 

Применяя эти зависимости к случаю, который обсуждался выше и описывался уравнениями (33), получаем

1,715М?, +

0.0065МГ, — 0.429Л4' +

0,478М$ =

— 0,3487;

0,85ЬМ/, — 1,712Л4о 4 0,1663Л4^ +

0.239М5 =

0.174Г;

— О.ЗвШ^-}- l,715Afi + 0,0065Mj = — 0,3487;

0,478М/, +

0,239М5 + 0 ,8 5 Ш з— 1,712М£ =

0.174Г

Отсюда имеем Л42 = —0,2107, МГ> = —0,2607, М13 — —0,2607, Мг3 = —0,2057.

Теперь легко подсчитать углы закручивания для каждого колена. Уравнение (12) дает величину закру­ чивания, когда защемление отсутству­ ет. Имеем б z= Tw. Уменьшение этой величины из-за влияния изгибающего момента М подсчитывается по форму­

ле

(15а).

=

Для первого колена найдем 6 '=

— 0,210 Ts. Для второго имеем

бо= — (0,260 + 0,260) Ts. Для третье­

го — 60 = — 0,205 Ts. Уменьшение

приведенной длины из-за влияния за­ щемления будет таким:

для первого колена

б '/б = 0,2\0s/w= 0,0305, т. е.

около 3%;

для второго колена

б0/б = 0,520s/w = 0,0775, т. е. около 7,8%;

для третьего колена

6Q/ 6 = 0,250s/cw = 0,0301, т. е. около 3% .

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

чительное. Для

среднего колена

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

половину

подсчитанного

для

разобранного на

стр.

120 числового примера при допущении,

что

колено

защемлено по

п

п

п

середине опорных шеек. Подобные резуль-

*1

У

таты могут быть получены для

любого

 

1f

 

трехколенчатого вала при у2 = у3 =

120°.

 

 

Теперь могут быть

рассчитаны изгиб

94сп

 

и кручение того же коленчатого

вала си­

 

лами Р, Q и S, приложенными к шатунным

 

 

 

шейкам.

Силы

Р и S можно

получить из

 

 

 

диаграмм давления

газа

и

инерционных

 

 

 

сил. Исследуя первый случай, изображен­

 

 

 

ный на

рис. 12,

где

максимальная

сила

 

 

 

действует на первое колено,

получаем

зна­

 

 

 

чения рассматриваемых сил.

представлен­

 

 

 

ных в табл. 1.

 

 

 

 

 

 

Масштаб

28,8'10*кг-сп

Номер

колена

1

2

3

Угловаяко­ ординатапо­ коле­ворота на,дагр

 

 

Номер

Р , кг

S . кг

Q. кг

колена

20

—29 500

13 300

3500

1

140

—4260

—2600

3500

2

260

— 1950

5000

3500

3

Угловаяко­ ординатапо­ воротаколе­ на,град

Р , кг

 

90

1370

210

— 31 750

330

— 3670

S, кг Q. кг

6010

3490

13 900

3490

2710

3490

Подставляя в уравнения (32) числовые значения (34) и значения сил, приведенные

в табл. 1, получаем

 

 

 

 

 

 

 

 

 

 

 

1,715М12 +

0,0065/^2 — 0,429М3' +

0,478A^ =

— 460

103 кг-см;

 

О.вбШ з — 1.712MJ +

0,1663Мз +

0.239MJ = — 237 103 кг-см;

 

 

 

 

 

 

 

— 0.38Ш1, + 1,715Л^ +

0,0065Л^ =

 

 

 

 

 

 

 

 

=

— 111

103 кг-см;

 

 

 

 

 

 

— 0,тм12+ 0,239^2 +

0,85 ш ' - 1,712МГ3=

 

 

 

 

 

mi Отсюда

=

7,45 • 103 кг-см.

 

 

 

 

 

имеем

 

103 кг-см;

 

 

 

 

 

 

 

 

М12 = — 297

 

 

 

 

 

 

 

 

Мг2 =

— 8,24 • 103

кг-см;

 

 

 

 

 

 

 

 

М 1Ъ= — 62,8 . Ю3

кг-см;

 

 

|

/ 1

 

 

 

 

М3 — 46,1

103 кг-см.

[

Л

 

 

После

подстановки этих значений в выраже­

 

 

ние (28) получим

103 кг - см;

 

 

 

 

 

 

т12 =

— 161

>СГ

 

Т

 

 

 

 

 

 

— 331

103 кг - см;

Tribal

 

^

К

 

 

т2 =

 

4*2

 

 

m3 =

— 89,3 • 103 кг •см;

z

>1

 

 

 

 

 

 

тг3 =

— 99,1

103 кг •смт

Масштаб

 

 

(

 

 

На рис.

13 представлено изменение из-

 

L.-

J

 

гибающего момента по длине первого колена,

 

 

 

28,8‘JO^KtCM

около которого, как предполагается, находит­

 

Рис.

15.

 

 

ся максимальное значение

изгибающего мо­

 

 

 

мента. На рис. 13, а даны моменты в плоскос­

ти колена, а на рис. 13, 6 — моменты в пло­ скости, перпендикулярной плоскости колена. Напряжения в любом элементе колена можно теперь подсчитать обычным способом.

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

Из уравнения (32) получим

М2 = 1 1 7 103 кг - см;

Мг2 = — 78,9 • 103 кг - см;

М2 =

— 251

103 кг •см;

М2 = — 27,4 • 103 кг •см.

Подставив эти величины в зависимости (28), получим

т12=

158 • 103

кг - см;

т2= 181 103 кг - см;

тз =

— И З-

103 кг см;

тг3 = — 274 • 103 кг - см.

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

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

The

Herbert

«pendulum» hardness tester.

The Engineer, 1923, vol. 136,

July, 6, p.

21.

Перепечатка: The pendulum

hardness tester. T i m o s h e n ­

k o

S. P.

The collected papers. New York — London — Toronto, McGraw-Hill

 

 

 

Publishing Company, Ltd, 1953, p. 363—365.

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

Пусть через е (рис. 1) будет обозначено расстояние от центра тяжести прибора Cgjxo геометрического центра шара О, 0 — переменный угол ко­ лебаний, Q — вес маятника, / — момент инерции прибора относительно оси, проходящей через О и перпендикулярной плоскости колебаний. При допущении о малости ей 0 уравнение колебаний будет иметь вид

/0 + eQQ = 0.

(1)

Используя обозначение для радиуса инерции

 

r = VlgiQ,

(2)

находим из уравнения (1) период одного качания в виде

 

Г = я - £ = - .

(3)

У eg

 

т. е. период качания обратно пропорционален У е. Можно ожидать,

что

тот же результат будет и в случае идеально упругой плоской поверхности, так как распределение давления на контактной поверхности в точке А сим­ метрично относительно вертикальной оси АО. В случае появления отпечатка в материале будут иметь место условия для колебаний, изображенные на рис. 2. Распределение контактного давления в этом случае уже не симметрично отно­ сительно вертикальной оси АО и плечо пары сил Q зависит не только от экс­ центриситета е, но также и от величины отпечатка в материале. Допустим,

1 Pendulum hardness tester. The Engineer, 1923, vol. 135, April, 13, p. 390—391.

что при колебании плечо пары сил изменяется и равно (еу+ б) 0, где б — по­ стоянная в ел и ч и н а , зависящая от величины отпечатка. Тогда уравнение ко­

лебаний будет

(4)

/0 + + б) Q0 = 0.

(5)

V(e + b)g

Период уменьшается с увеличением б, т. е. с увеличением отпечатка.

Можно также видеть, что колебания могут быть вызваны и при отри­ цательном знаке величины е, когда центр тяжести Cg находится выше гео­ метрического центра шара О, для этого необходимо только, чтобы е < б.

 

Т а б л и ц а 1

 

Т а б л и ц а 2

 

ю г

 

 

 

107

 

 

вычисленные

полученные

 

вы численные

полученные

 

по формуле

эксперимен­

 

по ф орм уле

эксперимен­

 

(5)

тально

 

 

(5)

тально

0 ,0 5

133

133

 

— 0,1

И З

И З

0,1

100

100

 

0

7 7

7 6

0 ,1 5

8 3

81

 

0,1

6 3

6 3

0 ,2

7 2

7 0

 

0 ,3

4 8

4 8

0 ,3

6 0

6 0

 

 

 

 

Применим теперь формулу (5) к результатам, полученным в экспери­

ментах со

стеклом1.

 

полученные для

е =

0,3 мм и е =

0,05 мм,

Сравнивая результаты,

по формуле (5) получаем б =

0,0137 мм, г =

10,6 см. Для тех же

значений

б и г по формуле (5) можно вычислить период качания Т для каждого зна­ чения е. Некоторые результаты, полученные подобным образом, даны в табл. 1. Можно видеть, что формула (5) находится в хорошем соответствии с результатами эксперимента. Применяя эту же формулу для случая быс­ трорежущей инструментальной стали2, получаем из эксперимента = 0,3 и е = — 0,1) б = 0,188 мм.

1См. рис. 3, стр. 390 работы, указанной в сноске на стр. 127.

2То же.

Результаты расчетов по формуле (5) при б = 0,188 мм и г = 10,6 см даются в табл. 2. Снова получено хорошее соответствие формулы (5) с эк­ спериментом.

Когда центр тяжести прибора расположен выше центра шара, формула для расчета периода Т примет вид

7 = л

г

(6)

V (6

 

e ) g

Легко видеть, что, когда е близко к б, прибор становится очень чувствитель­ ным. Небольшие изменения величины б, т. е. небольшие вариации размера отпечатка, будут оказывать очень большое влияние на период колебания Т.

МАТЕМАТИЧЕСКИЙ МЕТОД ОПРЕДЕЛЕНИЯ МОДУЛЯ УПРУГОСТИ

Mathematical determination

of the modulus of

elasticity.

Mechanical

Engi­

neering, 1923, vol. 45, N 4,

April, p. 259—260.

Перепечатка: Determination

of the modulus of elasticity.

T i m o s h e n k o

S. P. The

collected

papers.

New York — London — Toronto, McGraw-Hill Publishing Company, Ltd, 1953,

p.366—370.

Всвязи с публикацией Давида Гуэлбаума1 по проблеме математического метода определения модуля упругости, а также в связи с замечаниями к ней

В.Брайанса2 представляет интерес обсуждение способа решения этой задачи.

При изгибе стержня, например, с прямоугольным поперечным сечением

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

е = 6/pd0 = z/p.

(1)

Соответствующее напряжение равно (рис. 1, б)

 

р = еЕ = гЕ/р,

(2)

и, интегрируя по площади поперечного сечения, получаем хорошо известное уравнение упругой линии

М = EI/р.

(3)

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

жение получается из формулы (2) для z =

+ ft/2,

откуда ртах = hE/2p

или рlh = Е/ртах.

 

 

 

 

Тогда для железа Е = 2,1

106 кг/см2, а ртах =

2100 кг/см2*, Так что

//

2,1

Ю6

САА

 

Р/1~

2,1

10:* 2 ~

500'

 

1 G u e l b a u m D. Mathematical determination

of the modulus of elasticity. Mecha­

nical Engineering, 1922, vol. 44, N 12, p. 858—859.

 

 

- B r y a n s W. R. Mathematical determination of the modulus of elasticity. Mechani­ cal Engineering, 1923, vol. 45, N 2, p. 136.