Применение уравнений нейтрального равновесия к задаче о цилиндре, сжатом боковым давлением Текст научной статьи по специальности «Физика»

УДК 519.6

ПРИМЕНЕНИЕ УРАВНЕНИЙ НЕЙТРАЛЬНОГО РАВНОВЕСИЯ К ЗАДАЧЕ О ЦИЛИНДРЕ, СЖАТОМ БОКОВЫМ ДАВЛЕНИЕМ

Г.И. ВОЛОКИТИН, КВ. КОЗЛОВ, В.Б. ФЕДОСЕЕВ

(Донской государственный технический университет)

Рассматривается задача потери устойчивости полого кругового цилиндра. Используется теория малых деформаций, наложенных на конечную. Докритическое начальное напряженное состояние неоднородное. Закон состояния определен соотношением Блейтц-Ко.

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

Введение. Согласно феноменологическому подходу, принятому в механике сплошных сред, базовые соотношения для нелинейно-упругих тел включают конкретное представление удельной потенциальной энергии деформации. Предлагаемое выражение должно допускать описание основных нелинейных эффектов и приводить к доступному математическому рассмотрению простых деформаций. На свойства резиноподобных материалов сильно влияют химический состав и процессы изготовления. Для некоторых сортов резины подходит модель несжимаемого упругого материала, конкретизированная потенциалами Муни, Трелоара, Г.М. Бартенева и Т.Н. Хазановича, К.Ф. Черных и И.М. Шубиной [1]. Однако часто необходим учет сжимаемости материала. Как практический, так и теоретический интерес в этом случае представляет потенциал Блейтца и Ко. Целью данной работы является исследование устойчивости равновесия толстостенного цилиндра из этого материала при гидростатическом сжатии на основе трехмерных уравнений нелинейной теории упругости.

Определяющие соотношения и начальная деформация нелинейно упругого цилиндра.

Закон состояния задается выражением плотности потенциальной энергии деформации, предложенным в работах Блейтца и Ко [2, 3]:

¡V =

(хр

И

(1-М

(1)

Здесь /р /2, /3 - главные инварианты меры деформации Коши-Грина, коэффициенты

а, р, |а, - модули упругости. Константа а выражается через постоянные Ламе и коэффициент Пуассона:

X V

а = — = -

2|а, 1 - 2\

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

Далее будет использован закон состояния для напряжений в форме Фингера [1]:

Т =

81У

~дГ

(

\

81

77 дЖ 772 Р--------------^

81-,

где

'1 2 /

Т - тензор напряжений Коши; Е - единичный тензор; ^ - мера деформации Фингера. Выражение тензора Т, найденное с учетом потенциала (1), принимает вид:

Т = ^[^1Е + с°^-с1^2] . (2)

Обобщенные модули задаются равенствами:

т-о. , и(!-Р)

1 _ _ 7-

^ I ^

ТО. _ 22_

3 I

V -'з у

К‘-РК

с =

Рассматривается осесимметричная деформация кругового полого цилиндра.

Материальные координаты - цилиндрические координаты точки недеформированного цилиндра -г, ф, г. Координаты точки после деформации, возникающей за счет поверхностного нагружения равномерно распределенным давлением по наружной боковой поверхностям цилиндра, обозначены соответственно Я, Ф, Z. Место точки в актуальной конфигурации определяется выражениями:

Я = Я(г)ег +сЫъ, Ф = ф. (3)

Здесь ег, е 73 - базисные векторы цилиндрических координат; г, ф, г -радиус, азимутальный угол и высота соответственно. Введем обозначения:

о = д', Ь = — . (4)

г

Тогда градиент деформации, определяемой соотношениями (3), принимает вид

0 С д I д

УД =

Є,----------Ь—Є-----------------1-7,--------

Гг г Г (р &

Я = аегег + 6ефеф + ёцц .

V 7 Ч' ^ у

о

где V - набла-оператор в отсчетной конфигурации. Отсюда находим меру Фингера и главные инварианты:

^ = а2егег +62ефеф +ё\ц, 11=11{р) = а2 +Ъ2 +ё2, 12= а2 (Ь2 + ё2} + Ъ2ё2, 1ъ=а2Ъ2ё2 . (5)

Тензор напряжений является соосным мере Фингера и с помощью физических компонент [1] выражается в виде:

Т = С в в +0 в в +С 7,7, ,

г г г ффф г 3 3'

где, учитывая формулы (2) и (3), для физических компонент можно получить следующие выражения:

ст =

=

Із

3/2 (/31+“-У1/з1-“-/2+а2(у1/3+/1)-а4),

^ ~12 +Ь2Ыз +/:)-64) ,

Ф т 3/2

13

ц(1-р)

а, =

3/2 (1э1+а --72+^2(у173+71)-^4)- (6)

±3

В соотношениях (6) инварианты вычислены по формулам (5), а у1 = у^~р ■

Если отсутствуют массовые силы, условие равновесия в объеме имеет вид

д 1 д . д . . ч

е-----1—е--------1-7, — 1а е е +а е е +с 7,7, = 0.

г дК Я ф<ЗФ дгг г г ф ф ф '

Учитывая зависимости (6) и деривационные формулы [1], можно это условие свести к одному уравнению

с/ст а

- + ^К-С<р) = °- (7)

с/г гЪ

Поскольку в этом уравнении напряжения аг, сф определены формулами (6), то ясно, что (7) -

это нелинейное дифференциальное уравнение второго порядка. Приведем вариант уравнения, отвечающий упрощенному варианту потенциала (1), - материалу Ноулса-Стернберга

(а = 1/2, Р = О, V = 1/4):

„ъ\

= 0. (8)

йа а-Ъ а а~ а"

сіг Зг

Ь Ъ2 Ьъ

где

Условия равновесия на цилиндрических поверхностях имеют вид:

ГО, г = г

аг= в, (9)

{-р, г = гн

гв - внутренний радиус трубы; гн - наружный радиус; р - интенсивность давления. В ча-

1

стности, для материала Ноулса-Стернберга о. =ц 1-

а Ьс1

Краевое условие на торце удовлетворяется интегрально, в «смысле Сен-Венана»:

гнЪ{гн)

<2 = ^ о2аЬгс1г,

(Ю)

где <2 - продольная сила, приложенная к торцам цилиндра для равновесия. Окончательно задача Ламе для цилиндра определяется краевой задачей (7), (9). Причем параметр осевого смещения точек с/ задается в соответствии с условием (10). В упрощенном случае материала Ноулса-Стернберга задача Ламе для нелинейно-упругого цилиндра рассмотрена в [4].

Для решения нелинейной краевой задачи (7), (9) предложен численный метод, использующий конечные разности [5]. Выбираются равноотстоящие узлы, производные аппроксимировались центральными разностями. Узловые значения аргумента определяются формулами:

Г]=гъ+{]~ 2)к 7 = 2,3,...,Ы + 2, ^ = ^-.

Используем также еще две «законтурные» точки гх =гв-И, гм+ъ =гн+И . Обозначим х. значения

искомой функции Я(г) на выбранной сетке значений аргумента. Тогда, в частности, из (7) имеем систему Ж + 1 нелинейных алгебраических уравнений:

х,=^-(х,+1

3 г.

,3

Ъ ¥

где а =

'7+1

V-!

2/2

X,

Ъ =

Еще два условия определим из краевых условий:

Х{=ХЪ- 2/2

( У2/3

х2

кг2;

-2/2

1-

Рн

И

ч-1/З, у2/3

Таким образом, имеем систему N + 3 алгебраических уравнений с N + 3 неизвестными. Решение системы получено методом простых итераций [5]. В качестве нулевого приближения итерационного метода решения краевой задачи на сетке принимались узловые значения соответствующего решения задачи Ламе для цилиндра в линейной теории упругости. Такой выбор обусловлен тем, что при малой относительной толщине оболочки и небольших давлениях уравнения (7), (8) и (9) переходят в соотношения линейной теории, для которых известны точные решения [6].

Следует отметить, что в упрощенном варианте материала Ноулса-Стернберга допускается решение дифференциального уравнения (8) в аналитическом виде. Его можно переписать, учи-с/ а-Ъ й

тывая тождество — =---------—, следующим образом

йг

л

йа

Ж

1

з Л

= 0.

„ а _ с1а , & с1Ь а-Ь

Применим подстановку t = — .C помощью тождеств — = t + b— и — =-----

Ь с1Ь ай Ф г

ренциальное уравнение можно представить в виде:

с/г Л

г (^ - 1)г^2 + t + А)

последнее диффе-

Отсюда находим

г = С,

з

№

\t~l\2 (t2 +t + 4)8 где Cj - произвольная постоянная.

4л/Ї5 л/Ї5

(И)

_ с1г \с1Я

Далее, учитывая равенство — =-------, получаем другое дифференциальное уравнение

г I Я

dR

= -3-

Отсюда имеем

п ,At2 + t + 4)4 { 3 2t + l ,

R = Ln-------------;----exp —= arctan .— = C2C,(t).

2лУІ5 л/15

(12)

Выражения (11) и (12) дают общее решение исходного дифференциального уравнения (8) в параметрической форме. Для решения краевой задачи необходимо удовлетворить краевым условиям (9) и найти границы изменения параметра (в и (н(1в\->гв, \->гн). Эти значения пара-

метра t определяются из системы двух нелинейных уравнений, вытекающих из краевых условий (9):

^ У (О

з

У

^н)’

і+Еу

с(0

с(0'

Далее, изменяя параметр t от tB до tH, с помощью (11) и (12) можно определить зависимость R(r). В [4] указаны полученные ранее другими авторами параметрические формулы общего решения дифференциального уравнения (8). Однако соотношения (11), (12) имеют отличия, связанные с выбором параметра. Также в [4] не прояснен вопрос решения краевой задачи. Точное решение сходной нелинейной краевой задачи Ламе для сферы из материала Ноулса-Стернберга приведено в [7].

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

w = u(r, ф, г)ег + v(r, ф, г)е<р + w(r, ф, z)i3.

Применяя деривационные формулы, для градиента вектора добавочной деформации имеем равенство

Vw = ■

1 du а дг

-ее, +-

1

R

1 du

dv ^ 5ф

1 dw

1 dv

1

л-----------ui, л-------------eem +

ф 9 d dz а дг ф R

du

5ф

Л

--v

е е +

ф г

-h ег+-

1 dv

¥>+’

1 ÖW

-еггъ+-

1 ÖW

-Уз ■

(13)

d дг d дг а дг R дф Смежные формы равновесия цилиндра исследуем с помощью оператора А.И. Лурье © [1,6], который для материала Блейтца и Ко определен соотношением (10), полученным в работе [7]. Учитывая это соотношение и равенство (13), для тензора 0 имеем следующее представление:

0 =

л ^ л

А, — + А,

1 дг 2

(

и ■

ду

5ф

(

„ ди

С,-----------\~

1 дг 2

V у

ґ дул и

V

л дм А3 —

3 &

Л

(

егег +

/

д—+в2

1 дг 2

(

ду

5ф

■ вД

3 &

Л

е е +

ф ф

^— + ^2 1 дг 2

5ф

У

„ дм

+ С-о —

3

Л

(

/

ду

Ц — + Д

1 дг 2

ди

5ф

У

--V

/

е е +

Г ф

у У

5ф

--V

е е +

ф г

у У

дм ди Л

+ Сг2~ ЄЛ

дг дг)

(14)

Я,

дм

ди

1 - 1 2 - 3 г

бг

Н\ — \ие, + дг

Р,

дм

-д

Л

Уз +

Сдм

йт-+£2 —

5ф

дг

¥>•

1 5ф ^ дг г

Коэффициенты при неизвестных и, V, м и их частных производных известны после решения начальной задачи Ламе функции переменной г. Приведем выражение одного из этих коэффициентов:

И(1-Р),

Д =■

а1.

3/2

-((1 + 2а)у1/з1_“) + (2а- 1)/3(1+а) +Уга213 +ЗЬ2сі2).

Уравнения нейтрального равновесия в объеме [1] (У© = 0) приводят к системе трех линейных дифференциальных уравнений второго порядка с переменными коэффициентами:

а©

гг

дг Я

а©

©_ -©.

дг Я

\

д@г2 а

-----—+ —

дг Я

© +©

Гф фГ

а© ^

фГ

5ф

а© ^ фф 5ф

д д® г

ё дг а а©

= 0,

г<р

сі дг

= 0,

(15)

а©,„Л а а©

©

аф

сі дг

- = 0.

В уравнениях системы (15) компоненты тензора © выражены в соответствии с представлением (14). Нейтральное равновесие на боковой поверхности определяется условием [1,6]

Н® = -р{умЫ-Ы^мт), (16)

где V - набла-оператор в метрике деформированного состояния; N - вектор внешней нормали к деформированной поверхности, нагруженной гидростатическим давлением интенсивности р. На торцах принимаем граничные условия вида при г = 0 и г = Ь( ¿-длина цилиндра):

и, = 0, ©гг =©гг = 0 . (17)

Соотношения (17) означают недопустимость добавочного перемещения в осевом направлении и отсутствие трения на торцах [8]. Эти равенства не противоречат характеру нагружения цилиндра

наружным давлением и сжатия с торцов абсолютно жесткими гладкими плитами.

Разыскиваем решение краевой задачи (15), (16)

и =Хшп(г)со^п^со^Кг,

(18)

™ = 2т(.г)С0*пу$тХтг,

где п и X - параметры волнообразования, в соответствии с (17) определяется равенством

пт

Хт = —^~, а числа п, т - целые.

Подстановка (18) в уравнения (15) и краевые условия (16) приводит к разделению переменных. В результате получаем однородную систему обыкновенных дифференциальных уравнений:

АХ 1 +

Х'пт +П\ А + | Л + —Н\ \2'пт

А+^-в2-п2г2)-х2т^н2

А'2+-{А2-В2-Р2)

+х,„

А+—{А~въ)

г =о,

пт ’

Д7'

1 пт

А'-А+-(А+^)

71-«I А+-А \Х'т-

+п2ВЛ + -Х2ОЖт->

і' (Ц

Д +~(^2 +Р2+В2)

К

(19)

-пк \ -в^+-а \гт = о,

аг" +\о;+—а\г' -х \о,+-с,\х' -х \о!}+—о2+-с2\х -

1 пт | 1 ^ II Ш7 т I 2 (Ц ) ^ ж I 2 ^ (Л ^

-Х_и| -Р2+^-С, |7_-а

-І +ХІ-С, 12 =0

ті 1 т т 3 ят

/? а

Из краевых условий (16) после разделения переменных получим при г = гв и г = гн:

4^1 + ^4 +-^(Хпт +П^„т) + ^т^А + ^ = 0>

+7ит) = 0,

Д7' - Д -

1 пт 2

(20)

^1-Яи^2-^ = 0.

В условиях (20) ц = 0, если г = гв, и ц = р , если г = гн.

Однородная система (19) и краевые условия (20) выражают нейтральное равновесие цилиндра. При бифуркационных значениях параметров р и с1 возможны нетривиальные решения

краевой задачи (19), (20). Они соответствуют возмущенным равновесным формам цилиндра. Критические значения давления находим как наименьшие бифуркационные значения р при надлежащем выборе п и Хт. Бифуркационные значения - собственные числа задачи (19), (20). Отметим, задача на собственные значения является нелинейной, поскольку искомый параметр входит нелинейно через функцию Я (г). Для решения такой задачи на собственные значения использован численный метод [8].

Выводы. Численный алгоритм определения бифуркационных нагрузок, включающий решение задачи Ламе о начальной деформации цилиндра, вычисление коэффициентов системы (19) и по-

лучение собственных значений, реализован на ЭВМ. Расчеты проводились как для упрощенного варианта материала - материала Ноулса-Стернберга, так и для полного закона Блейтца и Ко. В качестве примера применения полученных уравнений нейтрального равновесия рассматривались некоторые случаи потери устойчивости, изучаемые в теории оболочек [9]. Для очень длинных цилиндрических оболочек под внешним боковым давлением критическая нагрузка разыскивалась в виде соотношения Грасгофа-Бресса: /; = 2ц(1 + у)р,г', где в - относительная толщина оболочки, р, - безразмерный параметр давления. Обнаружено, что критическому (минимальному) давлению отвечают плоские формы выпучивания: п = 2, т = О. Представления о влиянии нелинейности на величину верхнего критического давления в зависимости от относительной толщины отражены в табл. 1. Данные второй строки согласуются с результатами, приведенными в [8]. Численный анализ также показал, что для не слишком толстостенных цилиндров влияние на критическое значение параметра давления точности решения начальной задачи (задачи Ламе) и кон-

Г

станты р незначительное. Изучался случай оболочек средней длины, где принималось -^- = 0,38 .

Для таких оболочек критическая нагрузка разыскивалась в соответствии с формулой Саусвелла-Папковича: р = 2\1(\ + \)р„г2’5. Длина оболочки неизменная (¿/ = 1), цилиндр нагружен внешним

боковым давлением. Критические значения параметра давления и соответствующие значения параметров волнообразования содержатся в табл. 2, параметр т = 1.

Таблица 1

8 = 2(г„-г„)/(гн+гв) 0,05 0,10 0,20

р, (v = 0,272) 0,26 0,27 0,29

р, (v = 0,250) 0,34 0,36 0,38

Таблица 2

8 = 2(r„-г„)/(гн+гв) 0,05 0,10 0,20

п 4 3 3

р, (v = 0,272) 0,38 0,39 0,44

р, (v = 0,250) 0,49 0,52 0,56

Библиографический список

1. Лурье А.И. Нелинейная теория упругости / А.И. Лурье. - М.: Наука, 1980.

2. Blatz P.J., Ко W.L. Application of finite elasticity theory to deformation of rubbery materials //Trans. Soc. Rheol. 1962, v.6. - P. 223-251.

3. Ко W.L. Application of finite elasticity theory to the behavior of rubber-like materials. Thesis in partial fulfillment of requirements for the degree of doctor of philosophy. California Institute of Technology, 1963.

4. Fu Y.B., Ogden R.W. Nonlinear elasticity: theory and applications. Cambridge University Press, 2001.

5. Калиткин H.H. Численные методы / H.H. Калиткин. - М.: Наука, 1978.

6. Лурье А.И. Теория упругости / А.И. Лурье. - М.: Наука, 1970.

7. Волокитин Г.И. Условия бифуркации равновесия сферы / Г.И. Волокитин, Д.В. Моисеев. // Современные проблемы механики сплошной среды: тр. XV междунар. конф. - Ростов н/Д: Изд-во ЮФУ, 2011. - Т. 1. - С. 69-73.

8. Волокитин Г.И. Устойчивость нелинейно-упругого цилиндра при боковом давлении и осевом сжатии / Г.И. Волокитин // ПММ. - 1982. - Т. 42. - Вып. 2. - С. 289-295.

9. Григолюк Э.И. Устойчивость оболочек / Э.И. Григолюк, В.В. Кабанов. - М.: Наука, 1978.

Материал поступил в редакцию 26.03.2011.

References

1. Lur'e A.I. Nelinejnaya teoriya uprugosti / A.I. Lur'e. - M.: Nauka, 1980. - In Russian.

2. Blatz P.J., Ko W.L. Application of finite elasticity theory to deformation of rubbery materials // Trans. Soc. Rheol. 1962, v.6. - P. 223-251.

3. Ko W.L. Application of finite elasticity theory to the behavior of rubber-like materials. Thesis in partial fulfillment of requirements for the degree of doctor of philosophy. California Institute of Technology, 1963.

4. Fu Y.B., Ogden R.W. Nonlinear elasticity: theory and applications. Cambridge University Press, 2001.

5. Kalitkin N.N. Chislenny'e metody' / N.N. Kalitkin. - M.: Nauka, 1978. - In Russian.

6. Lur'e A.I. Teoriya uprugosti / A.I. Lur'e. - M.: Nauka, 1970. - In Russian.

7. Volokitin G.I. Usloviya bifurkacii ravnovesiya sfery' / G.I. Volokitin, D.V. Moiseev // Sovre-menny'e problemy' mexaniki sploshnoj sredy': tr. XV mezhdunar. konf. - Rostov n/D: Izd-vo YuFU, 2011. - T. 1. - S. 69-73. - In Russian.

8. Volokitin G.I. Ustojchivost' nelinejno-uprugogo cilindra pri bokovom davlenii i osevom szhatii / G.I. Volokitin // PMM. - 1982. - T. 42. - Vy' p. 2. - S. 289-295. - In Russian.

9. Grigolyuk E' .1. Ustojchivost' obolochek / E' .1. Grigolyuk, V.V. Kabanov. - M.: Nauka, 1978. -In Russian.

APPLICATION OF NEUTRAL EQUILIBRIUM EQUATIONS

TO PROBLEM ON CYLINDER WITH LATERAL PRESSURE COMPRESSION

G.I. VOLOKITIN, K.V. KOZLOV, V.B. FEDOSEYEV

(Don State Technical University)

The buckling problem of the hollow circular cylinder is considered. The small-strain theory superimposed on the finite is used. The subcritical initial stress state is inhomogeneous. The law of state is defined by Blatz-Ko. Keywords: finite deformation, stresses, circular cylinder, tensor, buckling, stability.