![](/user_photo/_userpic.png)
книги из ГПНТБ / Теория ветвления и нелинейные задачи на собственные значения
..pdf190 К. С Е Н СЕ Н И Г
Обозначая через еп то значение г, при котором может возникнуть п-я мода выпучивания, имеем
|
|
|
|
(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В |
2[х |
«13= |
|
А |
„ - |
|
|
|
|
|
g~ 11{гп)і |
|
|
|
|
||
«гі — |
|
-рг 11(rn) + |
( т іг ^ т ) |
||||
022 — |
|
ЗЛ + 2В 1 |
|
|
|
|
|
|
2ß |
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 |
2Ö |
-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д (Н - б)2 W» — (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. Пусть 2а — угол, на кото рый опирается недеформированная дуга, причем —а/К ^ т ^ а/К- Концы стержня шарнирно закреплены, так что в них отсутствует момент; следовательно, имеем граничные условия
М (±а/К ) = 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