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

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

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

удовлетворить, если принять следующие выражения:

 

Ху = -----^-Еу3хе~уг{а2 — х2){Ь2 — у2)-,

(11)

Хг =

-^-Еу2хе~уг(а2 — х2)у — 2цху;

(12)

¥ г =

± -Е ч Н е-^ ф *-у*)х.

(13)

При увеличении г выражения (9) — (13) стремятся к зависимостям (4), ко­ торые были найдены ранее. Значение у должно быть выбрано таким образом, чтобы сделать величину потенциальной энергии кручения V минимальной.

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

v=ЖIГГ[х>+*=+к- +

z‘)

О—а —b '

'

ипредположить, что |e~yzdz = — , то получим

оY

 

V =

Ех2а3Ь31— Зу +

(1 +

v) Г-|р а2Ь2уь +

(а2+ Ь2) у3+

 

 

 

 

 

+

_ 1 2

____L U

 

(14)

 

 

 

 

 

Ml+v)* а3 J}

 

 

Уравнение, из которого можно определить величину у, будет иметь вид

 

 

(1

+ v )

2

а*Ь»-у* + ± ( а 2 + Ь2)у2

= 3.

(15)

Для

случая очень

узкого

прямоугольника получим

 

 

 

 

 

 

 

a V = T q h r .

 

(15')

Для

того

чтобы определить

угол

закручивания 0,

приравняем

значение

потенциальной энергии, выражаемой формулой (14), величине работы крутя­ щего момента М. Тогда с учетом формулы (8') имеем

е . ф

 

или, положив V = 0,3, получим

 

 

е =

т(/ — 0,425а).

(16)

Влияние местной нерегулярности в поперечном сечении с координатой г = = 0 на величину угла закручивания 0 такое же, как и уменьшение длины I на величину 0,425а. Если же потребуется принять во внимание границы х = ± а , то выражения для Х 2, У2 должны быть заменены следующими:

Хг = -\ - Еу2те~уг (а2— х2)у — 2рту(\— е~к {а~~х)); Уг = -^ - Еу2хе~уг (Ь2у2) х + yak (b2у2) е~к(а“ х).

Величины а2у2 и уменьшение длины I на величину б/, соответствующее уменьшению угла закручивания 0, даны в таблице.

Решение обсуждаемой задачи может быть получено иначе, если зада­ ваться выражениями для перемещений. Эти выражения должны быть вы­ браны таким образом, чтобы удовлетворить условию, что поперечное сече­ ние с координатой г = 0 остается плоским. Они будут содержать одну или две постоянные величины, которые представляют собой координаты

alb

asY2

6/

системы.

Эти постоянные величины

можно

 

 

 

определить с помощью

вариационных

урав­

оо

3,846

0,425 а

нений равновесия.

прямоугольного по­

10

3,604

0,428 а

 

Для

случая

узкого

5

3,047

0,390 а

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

 

 

 

может рассматриваться

как тонкая пласти­

чае можно

предположить,

на,

заделанная

на крае

z = 0. В этом слу-

что перемещение v в направлении у выражает-

ся в такой форме:

 

Мх

 

 

ч

 

 

 

 

v =

z ----- —(1 — е

 

(17)

 

 

 

 

a

v

 

 

 

 

 

 

 

 

 

 

где М — крутящий момент; а — постоянная величина, которая должна быть определена; С — жесткость призмы при кручении.

Видно, что при этом условия (v)z=о = 0, (dv/dz)z=о = 0 удовлетво­ ряются. Далее видно, что с увеличением z относительный угол закручива­

ла

ния d2v/dxdz = — (1 е-*2) стремится к величине М/С. Пренебрегая малой

величиной e~al, получаем угол закручивания в виде

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

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

производит крутящий момент М. Таким образом, получаем

2aD аа£

+ 2(1 _ v)(; ----- =

с

 

где D — изгибная жесткость пластины. В случае узкого прямоугольного поперечного сечения имеем 2aD/С = 1/2 (1 — v), а из приведенного выше соотношения получим следующую формулу:

_1_

а

0,488а.

(19)

а

У 6(1

— v)

 

 

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

Вычисление, проделанное для случая, когда выражение для перемеще­ ния принималось в форме

v =

Мх

~ а

~ е ') + РЛАГ®*]

(20)

С

и содержало две константы а и р ,

дало результат, отличающийся от при­

веденного в формуле (19) только в последнем десятичной знаке.

В дальнейшем, для того, чтобы установить влияние жесткости мест­ ной нерегулярности на характер устойчивости плоского листа, изгибаемого в своей плоскости1 (рис. 2), воспользуемся выражением (17). Известно, что, если увеличивать изгибающий момент, может наступить такое состояние, когда плоская форма изгиба становится неустойчивой и появляется боковое выпучивание, как это показано на рис. 2. Такое выпучивание сопровожда­ ется кручением.

Для того чтобы установить критическую величину изгибающего момен­ та Мкр, создающего возможность для такого вида неустойчивости, возьмем дифференциальное уравнение равновесия. При допущении о малости пере­ мещений можно сделать вывод, что крутящий момент в поперечном сечении тп будет равен Mdy/dz и в соответствии с выражением (17) получим диффе­ ренциальное уравнение кручения в следующей форме:

« T - C (- S — l l r - S - ) .

Изгибающий момент в плоскости наименьшей изгибной жесткости Л40, и соответствующее уравнение равновесия будет

 

Л40 = — В

d2y

 

 

 

1ПГ

 

 

где В — изгибная жесткость.

 

 

 

 

Из уравнений (21) и (22) получим

 

 

 

d40

о

dW

М2а2

0 =

0.

dz*

а*

dz2

ВС

Интегрируя это уравнение и замечая, что 0 и d2Q/dz2 равны нулю при и х = I, находим

0 = A sin Рг,

где

P = / - - f « ’ + y - f + ^

- Т -

С учетом формулы (24) получим

(2 1)

равен

(22)

(23)

х= 0

<24>

м „ - - Ц З . j / ( и - ^ т )

или, учитывая формулу (19) для 1/а, имеем

Мкр= nVl BC (l + U 8 - ^ - ) .

(25)

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

1 M i c h e l l A. G. М. On the elastic stability of long beams under transverse forces. The London, Edinburgh and Dublin Philosophical Magazine and Journal of Sciences, ser. 5, 1899, vol. 48, N 292, p. 298— 309; Т и м о ш е н к о С. П. К вопросу о явлениях резонанса в валах. Изв. С.-Петербургского политехнического института, 1905, т. 3, вып. 1— 2, стр. 55— 106, Отд. оттиск: С.-Петербург, тип. А. Э. Коллинса, 1905, 52 стр.

МЕМБРАННАЯ АНАЛОГИЯ ПРИ ИЗГИБЕ

A membrane analogy to flexure. Proceedings of the London Mathematical Society, 1922 [1921], vol. 20, ser. 2, pt 5, p. 398—400; pt 6, p. 401—407. Перепечат­ ка: 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. 321— 328.

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

Предполагаем, что центральная линия балки длиной / располагается горизонтально и один конец ее заделан* а также, что силы приложены к по­ перечному сечению, совпадающему с другим концом. Они распределяются таким образом, чтобы быть статически эквивалентными вертикальной на­ грузке Р, действующей вертикально вниз вдоль линии, проходящей через центр тяжести поперечного сечения. Расположим начало координат на за­ деланном конце и направим ось z вдоль центральной линии, а ось х — вер­ тикально вниз. Далее предположим, что оси хи у являются главными осями инерции поперечного сечения. В этом случае в соответствии с решением Сен-Венана имеем

Xx = Ye = Xy = 0;

 

Zz = - P ( l - z ) - f .

(,)

Составляющие напряжений X и Y будут функциями только от х и у. Урав­ нение равновесия

дХг

,

дУ2 ,

dZz

= 0

дх

г

ду

дг

 

будет удовлетворено, если положить

X =

дер

Рх2

( 2)

ду

2/ + /(!/). у , — -£■

где ср обозначает функцию напряжения, а / является произвольной функцией только от у.

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

Условие того, что цилиндрическая поверхность стержня свободна от напряжений, имеет вид Xz cos (xv) + Уг cos (yv) = 0, оно может быть за­ писано следующим образом:

Подставляя выражения (2) в уравнения совместности деформаций V2KZ =

= О, V2X 2 = ----- j + v)/

, получаем

 

 

 

а2Ф

а2ф

v

Ру

- Г ( у ) + С.

(4)

дх2

ду2

1 + v

/

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

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

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

- • и - а д + ( * л * *

и работа нормального давления

6 Э Д ф [т р 7 ‘Т ' — f'(y) + c dxdU'

откуда следует, что функция <р может быть найдена из условия, что инте­ грал

S =

1

J+

ф( т Т 7 -7~ — Г (У) + с)} dxdiJ

(5)

(-|г)2+ (“tr)

минимален.

 

 

примем

 

Применяя метод Релея — Ритца,

 

 

ф = Яо'Фо +

Д1Ф1 + Я2Ф2 +

(6)

где ф0, фь ... — функции, обращающиеся в нуль на границе. Коэффициен­ ты а0, а1э ...— могут быть определены из условия минимума, который имеет вид

6S

= 0.

(7)

дап

 

 

Точность приведенного решения будет зависеть от числа членов ряда, пред­ ставленного выражением (6).

Если граница поперечного сечения задается уравнением F (х, у) = О, а функция F отлична от нуля внутри поперечного сечения, то решение (6)

можно взять в следующей форме:

т—со п=оо

(8)

<P= F(X, у) 2

2 атпхтуп

т= 0

п=0

 

Теперь покажем, как находить функцию <р, если граница

сечения балки

имеет ту или иную определенную форму.

Эллипс. Уравнение ограничивающей эллипс кривой имеет вид х2/о2 + + у21Ь2 = 1. Правая часть уравнения (3) будет равна нулю, если положить

Дифференциальное уравнение (4) будет таким:

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

выражения (8) дает точное решение уравнения (4а), если принять

 

_

Ра2

(1 +

v) а2 +

vb2

 

а°1

/

2 (1 +

v) (За2 + Ь2) 9

При а = b получим

решение для

круга.

 

Прямоугольник.

В случае прямоугольника границы задаются уравне­

ниями х = ± а, у =

±Ь. Члены, стоящие в правой части уравнения (3),

будут равны нулю, если положить / =

Ра2121. Соответствующее уравнение

для мембраны будет иметь вид

 

 

 

 

д2<р

д2ф __

v

Ру

 

дх2

ду2

1 + v

/

В этом случае ф также является четной функцией относительно х и нечет­ ной функцией относительно у. Эти условия и условия на границе будут удовлетворены, если взять выражение (6) в такой форме:

о )

т= 0 п=1

Из уравнения (7) можно получить

iyn+л—1

_

v

Р

Я 2 ш + 1 .л —

1 + v

/

,

1 а2 (2т + I)2 + п2

 

 

 

п(2т+ 1)

где а = Ыа. Касательные напряжения, выражаемые формулой (2), будут равны

У _

/4т

,

U

/ 2 ___ у2\

у ___

Ят

 

дф

Р

_

( Ю )

Лг

ду

+

2/

Х ^

Уг

дх 9

 

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

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

 

д2ф

__

v

 

Ру

Откуда получим

ду2

 

1 + v

 

/ •

 

 

 

 

 

ф =

1 + V -& -W -

Ь*у).

Тогда

 

 

6/

 

 

 

 

 

 

 

Х г

 

 

 

 

(П)

В центре поперечного сечения имеем

 

 

(Хг)х=у=о

 

2/

|l

3 ( i + v ) а2)

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

 

 

 

32ф _

v

Ру

*

 

 

 

Отсюда

 

 

дх1

Г+Т ~ Г

 

 

 

__

р

 

 

 

 

Р

 

 

Х 2

2- х

2),

Yz =

V

ху.

( 12)

 

1 + V

2/

 

 

 

1 + V

/

 

 

По сравнению с обычным элементарным решением касательные напряжения

понижаются в отношении 1/(1 + v).

dx/Xz = dy/Yz ли­

Можно указать, что дифференциальное уравнение

ний касательных напряжений в соответствии с формулами (12) дает

у = С(а2 — х2Г.

(13)

Выражение (12) будет точным решением задачи изгиба, если граница опи­ сывается уравнением (13) или, что то же самое, уравнением (y/b)l/v = 1 —

— х2/а2. В действительности члены, стоящие в правой части уравнения (3), будут равны нулю, если положить

Ра2

f(y) = 21

Соответствующее решение уравнения (4) будет иметь вид

V

Р

у(хг — а2) + а2Ь

/ и \!+1/v

1 —v

2/

( ^

а выражения (2) дают решение (12).

Для случая, когда а и b величины одного порядка, было использовано полное решение (9). Учитывая следующие соотношения:

V _ L - _ L TT2

9

V

( - 1)"

 

2 J Л 2 “ 6 п

2 J

Па

12. ’

( - 1)”

 

 

яз

2 кл-\

 

 

 

 

(2т+ 1 ) [(2т +

I)2 + fea]

32

 

4 ( i H ‘ '

получаем

 

 

 

 

 

 

 

 

 

 

(Xz)x=o,t/=o

=

ab

____ ^ _ а2 {J _

.

_ i_ у

( - ^ J

j

 

8

1 +

v

| 3

*

л12

2ш^

n2ch nnla

(X ;) I = 0,JI=J

 

 

1

 

0& 1 з

 

л2

S

1

|

g

 

v

 

n2 ch ял /а

Эти формулы совпадают с хорошо известным решением Сен-Венанаг. Используя метод Релея-Ритца, можно получить решение задачи в дру­

гой форме, более удобной для числового расчета. Ограничивая общее вы­

Координаты точек

hIf*

х = 0.

у = о

х =

0, у = ь

OJC

 

 

 

Приближенн ое

 

Точное

Приближен­

Точное

 

решение

ное решение

решение

решение

0,5

0 983

0,981

1,033

1,040

1

0,940

0,936

1,126

1,143

2

0,856

0,856

1,396

1,426

4

0,805

0,826

1,988

1,934

ражение (8) только двумя членами и полагая

<р =

2 — а2) (у2b2) х

X (Ау By3), из уравнений (7) получаем

 

 

 

 

 

А

v

р

 

 

 

 

11

+

а 2

 

 

 

T + V W

/ 1 , 3

1 И 1 I

8 \ , ~ т

9 *

 

 

 

 

^ 7

^

5

а 2 у ( 11

 

a 2

J

21

35а2

В = — 1 + v

Ш*

I 1

 

5

Н —

+ —

V

21

9

 

 

 

\ 7

^

35а2

 

 

 

а 2

) [ 11

^

а 2

)

Соответствующие составляющие касательного

напряжения (10) равны

(Хг)х=о.у=о = - ^ г +

Аа2Ь2,

(Хг)х=0.у=ь =

 

 

 

2a2b2 + Bb2). (15)

Для того чтобы определить точность приближенного решения, были вычислены напряжения в центре прямоугольника (х = 0, у = 0) и в точке с координатами = 0, у = Ь). При этом v бралось равным 1/4. Величина выражений в квадратных скобках (14) и соответствующие величины реше­ ния (15) приведены в таблице.

Видно, что если а и &— величины одного порядка, то приближенное решение (15) достаточно точное. В случае необходимости можно всегда по­ высить точность решения, увеличивая число членов в общем выражении (8).

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

1 S a i n t - V e n a n t B . Memoire sur la flexion des prismes sur les glissements transversaux et longitudinaux qui l’accompagnent lorsqu’elle nes’ opere pas uniformement ou en, arc de cercle, et sur la forme courbe affectee alors par leurs sections transversales primitivement planes. Journal de mathematiques pures et appliquees, Liouville, ser. 2, 1856, t. 1, p. 89— 189 (см. стр. 163). [Перевод на русский язык: С е н - В е н а н Б. Мемуар о кручении призм. Мемуар об изгибе призм. М., Физматгиз, 1961, 518 стр. (см. стр. 463].

2 В таком случае точное решение (14) неудобно для численного счета.

напряжений в следующем виде (удовлетворяя условиям на границе):

Ф T q r v ~1Г ^ — (Ь~У) к), (16)

где k — константа, которую надо определить из уравнения (7). Положив

e~bk =

0, получим уравнение

 

 

 

0,8а2(akf =

0,8а (akf + (0,4 — 2а2) (akf + 6аak — 7 = 0.

(17)

Если а

очень велико,

то имеем ak =

-^ -j/lO . Если а = 4, из уравнения (17)

получим ak = 1,298, а формула (10) даст

 

 

 

(Х2)х=о1У=ь =

2,038.

 

Этот результат только на 25% меньше точного решения, приведенного в таб­ лице. Можно сделать вывод, что в случае узкого прямоугольника прибли­ женное решение (16) является достаточно точным.

Методом, примененным для прямоугольника, можно получить прибли­ женное решение и в некоторых других случаях. Например, если уравнения границ таковы (рис. 1), что у = ±Ь , х2 + у2 + г2 = 0, то положим

f(y) = - % г ( г г — У*)-

Приближенное выражение для функции ср можно взять в таком виде:

Ф = (У2 — Ь2) (х2+ у2— г2) (Ay + By3),

где А и В могут быть вычислены из уравнения (7). Таким же путем задача может быть решена в других случаях, когда поперечное сечение ограничено вертикальными линиями у = ±Ь и двумя окружностями, симметрично рас­ положенными относительно оси у.

Треугольник. Если треугольник симметрично расположен относитель­ но оси у (рис. 2), то уравнение границ будет у + а = 0, х = =h tg а (2а — у). Если положить

f{y) = - J - t g a ( 2a — yf,

то члены, стоящие справа в уравнении (3), будут равны нулю и надо будет решать следующее уравнение равновесия мембраны, закрепленной по краям:

д2ср .

д2<р

v

Ру

tg2a -^ -(2 а — у) + с.

(18)

~дх*~ + ~di?

Т+ 7 “Т

 

 

Используем общее решение (8), подбирая константу с таким образом,

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

 

Если tg2a = v/1

— v, то решение задачи очень простое. Получим его

наложением на напряжения

 

 

 

х ; =

р_

Х г

+ -i-(2 a — у )'

(19)

 

2/

 

 

 

 

напряжений от кручения, вычисленных из хорошо известной функции на­

пряжения

 

 

ухт

 

 

 

 

 

 

Ф =

(у + а)\х*-----\ -(2 a -y f\

 

 

 

 

Соответствующие напряжения

равны

 

 

у"

5<Р —

ИТ

(х2 +

 

2ау — у2), Y' =

д(р

g -2 x (y + a). (20)

 

 

2а

 

 

 

дх

2а

Надо только подобрать значение fix. Крутящий момент, соответствующий напряжениям (19), будет

X'zt/dxdy =

Ра.

Напряжению (20) соответствует крутящий момент, равный

2 Ц q>dxdy =

 

27

рта*

 

 

 

 

5 /~ з

 

Условием, что крутящий момент равен нулю, будет

2

D , 27

■рта4=

0,

—F-

Ра Н------

 

5

 

5 / 3

 

 

откуда получим

 

 

 

 

 

 

 

рт =

2

 

/ 3

Р

 

 

 

27

а3

 

 

 

Подставляя значения для рт в формулы (20) и складывая выражения (19) и (20), получаем

Хг = Х'2 +

Х’г = Щ - Р \ - х 2 + а(2а-у)\;

Y, = у ; + Y\ = Щ - Р х { у + а).

Напряжения Хг в точках,

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

функцией вида

 

(**)*=<> = Щ - Р (2 а ~ у ) .

Наибольшее значение этого напряжения будет равно

(*z)х=0,у=’—а

2 /~ 3

р

9а*