УДК 539.3
Решение задачи бифуркации цилиндрической оболочки при комбинированном
нагружении
д.т.н. проф. Охлопков Н.Л., Нигоматулин Ф.В. Тверской государственный технический университет 8 (920) 696-49-95, Fedor.nigomatulin@mail.ru 8 (4822) 52-63-63, kafsm@yandex.ru
Аннотация. Рассматривается задача бифуркации тонкостенной круговой цилиндрической оболочки с учетом сложного характера деформирования в момент потери устойчивости при пропорциональном докритическом нагружении осевой сжимающей силой, крутящим моментом и внутренним давлением в девиаторном пространстве деформаций А.А. Ильюшина Э(3). Связь напряжений и деформаций принята в соответствии с определяющими соотношениями гипотезы компланарности. Для определяющих функций пластичности использованы несколько вариантов аппроксимаций, предложенных В.Г. Зубчаниновым[1].
Ключевые слова: пластичность, устойчивость, сложное нагружение, оболочка.
Решение задачи строится на основе теории устойчивости неупругих систем В.Г. Зубча-нинова [1]. Используется условие несжимаемости материала и условие однородности напряженного состояния в оболочке до момента потери устойчивости. Задача решается в геометрически линейной постановке.
Уравнения связи напряжений и деформаций в момент потери устойчивости оболочки принимаем в соответствии с определяющими соотношениями гипотезы компланарности, которые в скоростях принимают вид [2]:
= + (о' - Nt)S —, (г, j = 1,2,3), (1)
а
где: S = ds/ dS = Pt; t = cos j1; 3y = eiJ; SJj - компоненты тензора-девиатора напряжений;
3JJ - компоненты тензора-девиатора деформаций,
а - модуль вектора напряжений; ds / dS, N - определяющие функции пластичности; J1 - угол сближения ( cos j1 = s • p1 );
S - длина дуги траектории деформации. Символ с точкой наверху означает дифференцирование по обобщенному параметру времени d / dt = d / dS • dS / dt.
Зависимость s = Ф(3) = Ф(£) полагаем универсальной для простого нагружения. Дифференциальные уравнения равновесия элемента цилиндрической оболочки, потерявшей устойчивость, и уравнения совместности деформаций имеют вид [1]:
+ = о ^22 , дТп = Q
ЭХ, ЭХ2 ' ЭХ2 дХ} 9
Щ + 2 ЁМ^ + ^ + T22i22 + 2Тихи +±Т22= 0, (2)
дХ2 дХхдХ2 дХ\ ххкхх llKl1 xlKxi R 22 V ;
d2sn д\2 д%2 _ 1 д2 W дХ22 дХ2 дХЖ2 ~ R' дХ2
где: е - скорости деформирования срединной поверхности, % - скорости изменения кривизны и кручения срединной поверхности, W - функция прогибов оболочки. Первые два уравнения системы (3) будут удовлетворены, если положить:
ТП=ЕИЦ
11 8Х22
Т22 = Eh
д2ф
дХ*
Ти = -Eh
а2Ф дХхдХ2
(3)
где: ф - функция скоростей усилий, E - модуль Юнга, Ь - толщина оболочки. Решение основных уравнений задачи представляем в виде рядов Фурье
00 Л
к=1 ^
)
mpR
¥ 1
Ф = Х Bk si^(1rnXl - nX2 )? 1 m = r r l
(4)
к=l
где: ш, п - целые числа, определяющие число полуволн в направлениях Х1? X2 соответственно (XI - в направлении образующей, Х2 - в окружном направлении), Ь - длина рабочей зоны оболочки, Я - радиус срединной поверхности. Оболочку принимаем «длинной», шарнирно подкрепленной по торцам. В уравнениях (4) сохраняем по одному члену ряда [1].
В результате, окончательно получаем систему алгебраических уравнений задачи о собственных числах [1]:
s
-K
EN* 2sv12
+1-
.3N*0 2 giv
K*S* — 1 2
(l+r2 )2 -1
1 - &
. gl.
k2 +
9
8 liv
K* S*
"Г" _ -» -r* — *
^—-3 n* K*S*,
21 m 4 i
V — (1 + r2)2 -3(1 -Si
2 P
(5)
*
2
* )S* ,
l
где: i — 3R/h - гибкость оболочки.
СУ су* 2 СУ * ТУ * * 2 *
S* — S11r + S22 + 2ст12г, K* — s11 + s22r -2s12r,
P2* N2* . А * Д 3. ЛГ* N2*2 D2 1 D* P*2 n Eh3 Ф — -4--1, Ф— -Ф, gi — —— ---H g2 — -— -(P3—V), D —-.
P N* 2 D 2 3 N* 2 D 2 P1 9
Интегралы принимают вид:
Z — -. p h
2Z (6)
i i 2GP* — f Pz*n-1 dz*,2GN* — f Nz*m-1dz*.
m J n J
(7)
-1
-1
Полагаем, что в зоне пластической догрузки 31 = 0°(г = 1), в зоне упругой разгрузки = 180°(т = -1) и зона разгрузки примыкает к границе = -1 , где = • Для определяющих функций пластичности N и йа/dS принимаем аппроксимации предложенные В.Г. Зубчаниновым [1].
р р
(8)
Р Р
N — 2G(1 -ш) 0 <J £ p, P — 2G(1 -1) 0 <J <-, P — 2G -<J1 <p,
где: G - модуль сдвига, ш - параметр пластичности А.А. Ильющина, X - параметр разупрочнения.
Разбивая в (7) интегрирование по зонам и полагая, что зона разгрузки примыкает к поверхности Zp = — 1, получим:
2
К = -[1 - (-l)m-w(l - zpm )], m
Pi = -[1 - (-1)m -1(1 - z*pm)], (m = 1,2,3). m
(9)
Для вычисления координаты границы раздела зон z*, имеем уравнения [1]:
f — 2CS«iK;1; f — P*Z* -P*. (10)
При использовании модифицированной теории устойчивости A.A. Ильюшина [1] интегралы P* рассчитываются по формуле (9),
а для Л/т имеем выражение [1]^
1
К = -[1 - (-1)m ](1 -w).
m
(11)
При расчете по модифицированной теории устойчивости В.Г. Зубчанинова имеем систему алгебраических уравнений [1]:
а „ Ер
Eg2
- K - -
2sv1
1 m (1+г2)
cv =
2 m
(12)
21 m
v = (1 + r2)2
Интегралы P*m рассчитываются по формуле (9).
Решение бифуркационной задачи позволяет для заданной комбинации полуволн m, п вычислить критическую гибкость оболочки i в зависимости от значения модуля вектора напряжений о в момент потери устойчивости.
Рассматриваются траектории пропорционального докритического деформирования оболочки осевой сжимающей силой, крутящим моментом и внутренним давлением в девиа-торном пространстве деформаций A.A. Ильюшина Э(3) (рисунок 1).
Рисунок 1. Траектории Рисунок 2. Аппроксимация диаграммы
деформирования деформирования стали 9Х2
В качестве материала оболочек принимаем сталь 9Х2 [3]. Диаграмма деформирования материала при осевом сжатии и ее аппроксимация в зоне упруго пластических деформаций представлены на рисунке 2.
В работах [4, 5] показано, что для процессов пропорционального нагружения, реализуемых в плоскости Э1 - Э3 (при сжатии и закручивании оболочки) расчеты по теории устойчивости А.А. Ильюшина, с учетом разгрузки материала, [1] позволяют, для ряда
конструкционых сталей, получить физически достоверные результаты. Предположительно, для данных процессов можно использовать и более простые аппроксимации определяющих функций пластичности, что сделано в настоящей работе.
Расчеты выполнены для нескольких траекторий пропорционального нагружения оболочки в плоскости Э1 - Э3. На рисунках 3, 4, 5 представлены кривые наименьших гибкостей, построенные как огибающие кривых устойчивости, построенные при различных комбинациях полуволн ш,п.
Рисунок 3. Огибающие кривых устойчивости полученные по модифицированной
теории устойчивости А.А.Ильюшина
Рисунок 4. Огибающие кривых устойчивости полученные по модифицированной
теории устойчивости В.Г.Зубчанинова
Рисунок 5. Огибающие кривых устойчивости полученные по теории устойчивости
А.А.Ильюшина
На рисунках 6, 7, 8, 9 кривая 1 соответствует модифицированной теории устойчивости А.А. Ильюшина, кривая 2 - модифицированной теории устойчивости В.Г. Зубчанинова, кривая 3 - теории устойчивости А.А. Ильюшина.
Анализ полученных результатов показал, что из рассмотренных вариантов решения наименьшие значения критических напряжений для оболочек реализуется в случае использования модифицированной теории устойчивости В.Г. Зубчанинова, что вполне ожидаемо. Расчеты по теории устойчивости А.А. Ильюшина и ее модифицированного варианта практически совпадают. Наименьшие значения критических напряжений реализуются в случае чистого кручения оболочки.
Рисунок 6. Огибающие кривых устойчивости при ai=ö и a2=0
Рисунок 7. Огибающие кривых устойчивости при a1=3ö и a2=0
Рисунок 8. Огибающие кривых устойчивости при a1=6ö и а2=0
Рисунок 9. Огибающие кривых устойчивости при a1=9ö и а2=0
Выполненные тестовые расчеты в дальнейшем будут распространены на спектр конструкционных сталей [3] и сопоставлены с результатами расчетов, полученных при использовании более сложных аппроксимаций определяющих функций пластичности [4, 5].
Литература
1. Зубчанинов В.Г. Устойчивость и пластичность. Т. 1. Устойчивость / В.Г. Зубчанинов. -М.: Физматлит, 2007. - 448 с.
2. Зубчанинов В.Г. Математическая теория пластичности: Монография / В.Г. Зубчанинов. -Тверь: ТГТУ, 2002. - 300 с.
3. Зубчанинов В.Г. Экспериментальная пластичность: Монография. Книга 1. Процессы сложного деформирования / В.Г. Зубчанинов, Н.Л. Охлопков, В.В. Гараников. - Тверь: ТГТУ, 2003. - 172 с.
4. Охлопков Н.Л. О влиянии сложного характера деформирования в момент потери устойчивости на критические параметры напряжений круговой цилиндрической оболочки / С.А. Соколов, Н.Л. Охлопков // Вестник Тверского государственного технического университета: Научный журнал. - Тверь: ТГТУ, 2008. - Вып. 13. - С. 229-234.
5. Охлопков Н.Л. О предельных поверхностях критических напряжений и деформаций материала в решении задачи устойчивости круговой цилиндрической оболочки при простых процессах / Н.Л. Охлопков, С.В. Черемных // Фундаментальные и прикладные проблемы техники и технологии. - Орел: ОрГТУ, 2012. - Вып. 5 (295). - С. 30-36.