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

книги из ГПНТБ / Теория ветвления и нелинейные задачи на собственные значения

..pdf
Скачиваний:
8
Добавлен:
24.10.2023
Размер:
9.31 Mб
Скачать

190 К. С Е Н СЕ Н И Г

Обозначая через еп то значение г, при котором может возникнуть п-я мода выпучивания, имеем

 

 

 

 

(5.6)

 

1 +<а + « т З г -

Если п -*■ оо, то е„ -н>------ '

> так

чт0 все Ф°РМЫ выпучива­

ния могут возникнуть прежде,

чем деформированный радиус умень­

шится до величины

^ R —

' где° ~ К0ЭФФиРиент Пуассо­

на. В частности, все формы

выпучивания могут возникнуть рань­

ше, чем деформированный

радиус станет

равным R/2.

Величина еп, вычисленная

согласно

формуле (5.6), согласуется

с точностью до членов низшего порядка по knh со значением, которое

дает теория тонких пластин. Из (5.6) мы видим, что еп есть функция

от 4klh2 и коэффициента Пуассона а. На рис. 3

изображен график

величины — е„ в зависимости от 4k„h2 для о = 0,3.

Прямая

линия —

соответствующий график, полученный из теории

тонких

пластин.

Эта прямая касается кривой в начале координат.

 

 

6. Решение задачи 2

Ради простоты мы ищем решения, для которых недеформированная ось цилиндра переходит в плоскую кривую. Тогда можно без потери общности предположить, что ось переходит в кривую, выходящую из начала координат на плоскости Ѳ = 0. Предположим, что и1 — четная функция от Ѳ, и2 — нечетная функция от Ѳ и и3 — четная функция от Ѳ. Предположим также, что и1, и2 и и3 могут быть гладким образом продолжены ниже плоскости z = 0 так, чтобы функция ы1 была четной по г, и2 — четной по г, а и? — нечетной функцией от г.

XI. Н Е К О Т О Р Ы Е ЗА Д А Ч И О В Ы П У Ч И В А Н И И ТЕО РИ И УПРУГОСТИ J9J

Поэтому запишем следующие разложения функций и1, и2, и3.

 

и1 = аг-\-'2і

2

amn(r)cosmQcosknz,

 

7 7 1=0 7 1 = 0

 

foo

оо

Ьтп (г) sin тѲ cos knz,

 

«2= 2

2

(6. 1)

771=0 7 1 = 0

оооо

w3 = f e + 2 2 cmn{r)cosmQsmknz,

то=0 n= 0

kn ■—Я П

Подставляя (6.1) в (4.5) и (4.7), убеждаемся в том, что переменные разделяются. Из (4.5) получаем

 

 

 

 

 

г @тп

ri

O-mn^

 

r2

®тп

a-\-b АукпИГ!

 

 

 

 

( Alt +

~2 Аз )

tn b m n А3 г Ьтп'

(•^5 Н

а-\-Ь

) ^n^mn =

~2 ^3 (

+

 

— frnm j — (И3 -f- 7Ц)

Ь,тп

а + Ь A a fln b m n '

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

т

 

 

 

 

 

 

— (A + Y ^ s) у-гömTl— ( Y

+ ^4 ) 7з~

 

 

 

 

 

 

 

 

 

-

( Л5 +

 

а

 

m k n

'■0,

(6.2)

 

 

 

 

 

 

 

- ^ х Л )

 

'

Ь

(

н

1

/

trfi

 

\

 

 

 

п

 

 

 

 

-^Э ( сяш + — Cmn----^

Cmn j — (Л6 -f- A7) k„^п^тп

 

 

 

 

 

 

 

— ^ ^5 ”1"

a - \ - b

 

) ^

11

( ^mn“I- ~ a m n ~bt n b m n j = 0.

 

Из нетривиальных краевых условий (4.7) выводим следующие

соотношения:

 

 

 

 

 

 

 

 

 

 

 

 

{А3-\~ A^) йтп-\-А^^тЬтпА"— \атп^

AbknCmn — 0

при г R,

 

 

 

 

 

Ьтп ^ а

т п = о

при

r = R,

 

 

(6.3)

 

 

 

 

 

bcmyi

akncimn = 0

при

 

г = R.

 

 

 

Требование, чтобы существовало ненулевое решение однородных уравнений (6.2) с однородными граничными условиями (6.3), опре­ деляет значения а, Ь, е, 6, при которых может возникнуть тп-я форма выпучивания. Хотя очевидно, что данная процедура может быть про­ изведена при любой функции плотности энергии деформации, мы снова приведем результат лишь для функции плотности энергии деформации частного вида (5.5). В этом случае мы обнаруживаем, что ненулевые

192

К. СЕ НСЕ НИГ

решения существуют только при т = 1 и п Ф 0. Условие существова­ ния ненулевого решения имеет вид

«11

«12

«13

«21

«22

«23 = 0,

Ö31

«32

«33

где

а і і = - ^ Г В Д ) - ^ / і ( г £ ) 1 ,

 

 

П

*-

 

' П

 

-1

« 1 2

 

/ З.Л-[-20

1

.

А

\ г/

 

 

 

 

 

 

 

 

 

 

 

 

Л+ 2В

 

 

 

 

+

(

2[х

«13=

 

А

-

 

 

 

 

 

g~ 11{гп)і

 

 

 

 

«гі —

 

-рг 11(rn) +

( т іг ^ т )

022 —

 

ЗЛ + 2В 1

 

 

 

 

 

L

 

 

Гг,

«23—

0

Гп

ГЛ (Гп)

 

~

(Гл)1 >

 

L

 

т„

 

J

«31=

--- р — I 1 ігп),

 

 

 

 

 

1п

 

 

 

 

_ А, + 2[д, а г- / ч

 

 

 

«32—

 

~

^ ІіѴп),

 

 

 

«33—

2Л ,, .

 

 

 

 

g

I l ( r n)-

 

 

 

 

При этом

 

 

 

 

 

 

 

 

 

 

 

A = 1 + - ^ L - ,

 

 

 

 

 

 

'

a - \- b

В

2fib a+ b

Гп — k nr ,

(6.4)

I

 

ЗЛ+2В

1

-k )h (r n ),

Гп >

(6.5)

Л

-I

2В Ii {fn)i

( 6-6)

* в - гХ = — г„.

Здесь / 0 и / і — модифицированные функции Бесселя первого рода.

Снова над величинами а, Ь, е, 6, фигурирующими в (6.5), (6.6), нужно

о

было бы написать кружочек. Обозначим через 8П величину б, опреде-

XI. Н Е К О Т О Р Ы Е З А Д А ЧИ О В Ы П У Ч И В А Н И И Т Е О Р И И УПР УГ ОСТИ I9 3

ляемую уравнением (6.4), которое можно решить относительно —6П

численно. На рис. 4 дан график —ö„ в зависимости от (— Для

ряда значений а. Прямолинейный график служит графиком —6„,

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

7. Заключение

Приведенные здесь результаты частично основываются на статье [7]. Другие статьи, посвященные точной теории подобных задач, перечис­ лены в литературе.

ЛИТЕРАТУРА

[1]Biot М. A., Exact theory of buckling of a thick slab, Appl. Sei. Res., 12, Section A (1963), 183—198.

[2]John F., Plane strain problems for a perfectly elastic material of harmonic

type, Comm. Pure Appl. Math., 13 (1960), 230—296.

[3]Lubkin S., Determination of buckling criteria by minimization of total ener­ gy, N.Y.U., CIMS, Res. Rep. IMM-NYU, 1957, pp. 241.

[4]Lur’e A. I., Bifurcation of equilibrium of a perfectly elastic body, Jour. Appl. Math. Mech., 30 (1966), 855—869.

[5]Read H. E., On the stability of equilibrium of thick and thick-walled isotro­ pic elastic solids, Ph. D. Thesis, University of Delaware, 1964.

[6]Sensenig С. B., The buckling of a thick circular plate using a non-linear theory, N.Y.U., CIMS, Res. Rep. IMM-NYU., 1959, pp. 262.

[7]Sensenig С. B., Instability of thick elastic solids, Comm. Pure Appl. Math.,

17

(1964),

451—491.

stability of a circular tube and under end thrust, Quart.

[8] Wilks E. W., On the

J.

Mech. and

Appl.

Math., 8 (1955), 88—100.

[9*]Лурье A. H.,

Теория упругости, «Наука», M., 1970.

ІЮ*]Савин Г. Н., Койфман Ю- И., Общая нелинейная теория упругости, При­ кладная механика, 6, № 12 (1970).

[ 11*]Лурье А. И., Теория упругости для полулинейного материала, ПММ, 32, № 6 (1968).

13 0 1 2 8 5

X I I

ФОРМЫ РАВНОВЕСИЯ НЕЛИНЕЙНО УПРУГИХ СТЕРЖНЕЙ

Стюарт Антман

1. Введение

Теория Эйлера — Бернулли плоского изгиба нерастяжимой эластики характеризуется определяющим соотношением

М — const-(k К)-

(1-1)

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

Рис. 1.

что решение соответствующих уравнений можно найти в квадратурах, которые выражаются при помощи эллиптических функций типа Якоби и Вейерштрасса (см. лекции I и IV).

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

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

Исследовать такие задачи было бы удобно, если бы их решения можно было получить при помощи квадратур. Этого, однако, не удается сделать в том случае, когда в дополнение к (1.1) выполняется

определяющее соотношение

 

. N = const -е,

(1.2)

,где N ■— результирующее осевое напряжение, а г — либо растяжение (изменение длины на единицу длины), либо материальная деформация кривой, задающей стержень.

XII. ФОР МЫ Р АВН О ВЕ С ИЯ Н Е Л И Н Е Й Н О У ПР У Г И Х С Т Е Р Ж Н Е Й

195

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

(i)Существует ли простая форма определяющих уравнений, кото­ рая приводит к задачам, интегрируемым в квадратурах?

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

Как показано в [1] и [2], ответ на эти вопросы утвердительный. В самом деле, одни и те же уравнения дают ответ на оба вопроса (і), (іі). Если т — длина дуги недеформированного стержня, s (т) — длина дуги деформированного стержня, Ф (т) — угол наклона каса­ тельной к недеформированному стержню и <р (т) — угол наклона касательной к деформированному стержню, то эти определяющие уравнения имеют вид

М =

Ы = Ж '

W =

б, т),

(1.3)

причем

 

 

 

 

б = s' — 1, |і = s'k К = (1 +

б) k — К = —(Ф — Ф)'.

(1.4)

Здесь штрих означает дифференцирование по т.

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

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

для энергии деформации W положительно определен.

Можно показать [1], что (1.1) и (1.2) обеспечивают хорошую аппрок­ симацию уравнений вида (1.3), когда растяжение мало.2

2. Формулировка основных уравнений

Уравнения равновесия можно получить из рис. 2, где изображено равновесие сил, действующих на выделенный элемент стержня. Здесь через Q обозначена результирующая срезающая сила, а через q — нормальная нагрузка.

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

13*

196

С. АНТМАН

с использованием недеформированной длины дуги т, то мы получим

+

+

=

(2.1)

(р + М)ЛГ + (1+б)М ' =

0.

(2.2)

К этим двум уравнениям добавим определяющие уравнения

М = Wß, N = W6

(2.3)

Рис. 2.

и геометрические соотношения

 

 

Ф' = ц +

К,

(2-4)

х' = (1 + б) cos ф, у'

= (1 + б) sin ф.

(2.5)

Правила выбора знака для этих величин указаны на рис. 2. Урав­ нения (2.1)—(2.5) составляют замкнутую систему, к которой надо присоединить подходящие граничные условия. Функции К (т) и q (т) считаются заданными.

3. Решение уравнений

Сейчас мы изучим тот случай, когда К и q постоянны, а W не зави­ сит явно от т. Мы ищем общее представление решения системы (2.1)— (2.5). Умножая (2.1) на М7(1 + б) и интегрируя полученное выраже­ ние, находим

М ' 2 + (1 + б)2 N2 — 2 qM (1 + б)2 = а (1 + б)2,

(3.1)

где а — произвольная постоянная интегрирования. Это уравнение справедливо для любого определяющего соотношения. Подстанов­ ка (2.3) в (2.1) дает

W 2= а (1 + б)2 + (Н - б)2 — (1 + б)2 Щ ,

(3.2)

а подстановка (2.3) в (2.2) дает

(ц + К ) ^ + ( 1 + б ) Щ = 0.

(3.3)

Так как величина К постоянна, (3.3) можно переписать в виде

[(м-+ К) W J + [(1 + б) W d' - W >' - №6ö' = о,

(3.4)

XI I. Ф ОР МЫ Р А В НО В ЕС И Я Н Е Л И Н Е Й Н О У П Р У Г И Х С Т Е Р Ж Н Е Й

197

 

 

 

 

 

 

 

а так как W не зависит явно от т, то

 

 

 

 

 

 

i r = IP > ' +

lF66'.

 

 

(3.5)

Таким образом, (3.4) можно проинтегрировать,

что дает

 

 

 

 

([l + K)Wil+ (l+ ö )W 6- W = b.

 

(3.6)

Здесь b — произвольная постоянная

интегрирования. Уравнение

(3.6) есть просто алгебраическое соотношение между р и б.

 

 

Теперь упростим (3.2). Сначала запишем (3.3) в виде

 

 

(р +

К) [IV P ' + ^ ß ö '] + (1 + б) [WVp' + Wbbб'] = 0.

 

(3.7)

Предполагая,

что

 

 

 

 

 

 

 

 

(р + /С) ^ 6 + ( 1 +б) І^бб^О,

 

(3.8)

из (3.7) получим

 

 

 

 

 

 

 

_А'.... (ц+К) WW+U + б)

,

 

 

Поэтому

 

 

(ц+ Ю^ив + О + в)^««

>*•

 

 

 

 

( 1 + б ) [ ^ г б0_ ^ б] ,

 

 

 

 

 

 

(3.10)

 

 

 

+ У ) ^ 6+ + б)\Гб6 И- •

 

Следовательно,

(3.2)

принимает вид

 

 

-

(3^.11) j

 

 

+

Г (ц + /С) ^мб + (1 + б ) Г 6б1 2

р '2 =

2qWtx — Wl) L

V

s a

Знаменатель правой части (3.11) не обращается в нуль ввиду поло­ жительной определенности гессиана W. Если из уравнения (3.6) выразить б как функцию от р, что возможно, если выполняется (3.8),

то правую часть (3.11)

можно

считать функцией

одного только р,

и, следовательно, (3.11)

дает уравнения траектории р в фазовой пло­

скости. Кроме того, если мы

обозначим

правую

часть (3.11) через

/ (р), то (3.11) имеет интеграл

 

 

 

т = ± j

/ -1/2 dp +

const,

(3.12)

который дает неявное уравнение для р как функции от т. Если эта функция известна, то можно найти б из (3.6), а затем получить ф (т), X (т), у (т) дальнейшими интегрированиями.

В случае когда условие (3.8) нарушается, можно предположить, что числитель в (3.9) не обращается в нуль. Тогда мы можем проделать аналогичные операции, которые приводят к следующему уравнению для б в фазовой плоскости:

б'2 = (а + 2qWil— ITб) (и+*)У и|1 + (1 + а)Гмв

(3.13)

 

] •

Правая часть этого уравнения теперь рассматривается как функ­ ция от б.

198

С. АНТМАН

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

4. Пример

Положим W = Е Іу 2/2 + ЕАЬ2!2, где Е І и ЕА — заданные постоян­ ные. Тогда определяющие уравнения линейны относительно дефор­

маций р и б: М = E ly, N = ЕА8. Так как s' =

1 + б >- 0 (условие

непрерывности), (3.6) можно переписать в виде

 

 

1 + б = -

(ц + К ?И А \Х/\

 

(4.1)

где В — произвольная постоянная. Тогда (3.11) принимает вид

 

(£ /р ')2 = — (р + К)2НА] X

 

 

 

X + 2£/<7р — (ЕА)2 [(В — (р + К)2 НА)42 — I]2}.

(4.2)

Замена переменных р + К =

(АВ/І)112 sin ф,

w = tg (ф/2)

пре­

образует (4.2) в уравнение, которое можно проинтегрировать при помощи эллиптических интегралов Вейерштрасса, но этот факт не помо­ гает, по-видимому, ни вычислениям, ни физической интерпретации. Поэтому дальше мы займемся изучением одной частной краевой задачи с применением методики, допускающей непосредственное обобщение. Разберем задачу, соответствующую рис. 1. Пусть — угол, на кото­ рый опирается недеформированная дуга, причем —а/К ^ т ^ а/К- Концы стержня шарнирно закреплены, так что в них отсутствует момент; следовательно, имеем граничные условия

М (±а/К ) = 0.

Отсюда следует, что

 

р (±а//С) = 0.

(4.3)

Кроме того, мы имеем геометрическое условие неподвижности концов. Это дает условия

X ( а ! К) X (—a l l С) = 2sin а ! К,

(4.4)

у ( а ! К) - у ( - а ! К) = 0.

С применением (2.5) условия (4.4) можно записать в виде некоторых интегральных выражений от <р.

Вводя новые безразмерные переменные

I = (р + К)а/К, t = Кх!а, ß2 = КЧ /а2А

(4.5)

р = a3q!K)3EI,

XI I.

Ф О РМЫ Р А В Н О В Е С И Я Н Е Л И Н Е Й Н О

У П Р У Г И Х

С Т Е Р Ж Н Е Й

199

можно написать формулы

(4.1),

(4.2),

(4.3)

в виде

 

 

 

1 +

8 =

(у2 — ß2£2)'/2,

 

(4.6)

( - § ) 2 = (1+ б)2 [а0 + 2Р (S-

а) -

82/ß4] =

 

 

 

=

(у2 - П 2) [«о+ 2р (£ -

а) - [(у2 -

ß2£2) -

1 ]2/ß4] =

 

 

 

 

=

A(Q,

 

(4.7)

 

 

£ (±1) =

а.

 

 

(4-8)

Здесь у2 и а0 — произвольные постоянные. Из (4.6) и (4.8) мы видим, что растяжение при t = ±1 одинаково; обозначая его через 80, из (4.6)

находим

(4.9)

У2 = (1 + 60)2 + ß2a 2.

Первоначальная система дифференциальных уравнений, соответ­ ствующая (2.1)—(2.5), допускает постоянные решения. Однако при помощи этих уравнений можно легко показать, что постоянное реше­ ние £ = а, удовлетворяющее (4.8), тривиально, так как оно может удовлетворять остальным граничным условиям (4.4) только при

^5 = 0 и р = 0.

(Для некоторых других типов граничных условий тем не менее существуют нетривиальные постоянные решения.)

Теперь изучим функцию h (£), определенную в (4.7), чтобы выяс­ нить ограничения, налагаемые граничными условиями (4.8) на при­

роду решений нашей системы. Из (4.7) имеем

 

h (а) = (1 + 60)2 (а0 - W ) -

(4.10)

Тогда из уравнения (4.7) следует, что для нетривиальных решений h (а) > 0. Кроме того, формула (4.7) показывает, что h (£) имеет нули при £ = rfcy/ß. Далее, так как £ должна иметь одно и то же значение а на обоих концах стержня, из теоремы Ролля вытекает, что сущест­ вует значение t = Т в открытом интервале (—1, 1), для которого d£/dt = 0, но ö ? X J d t 0. (Если бы здесь d%!dt2 = 0, то из теоремы единственности для уравнения второго порядка, соответствующего (2.1), следовало бы, что решение тривиально.) Если условие 1 + 8 > 0 не нарушается, то эти условия означают, что h (£) имеет простой нуль в £ (Т), причем | £ (Т) | < y/ß. Если £t — наибольший нуль функции h (£), меньший а, в котором h (£) меняет знак, а £2 — наименьший нуль h (£), больший а, в котором h (£) меняет знак, то по крайней мере один из нулей £ь £2 находится в открытом интервале (—y/ß, y/ß) оси £ и h (£) ведет себя, как показано на рис. 3.

Когда h (£) имеет такой вид, качественная природа решения £ (^)

хорошо известна (например, из анализа в фазовой плоскости). В част­ ію

ности, функция £ (t) периодическая с полупериодом ы = \ /i(£)_1/2d£,

V

U

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