CC BY
15
УДК 539.3
Хамитов Т.К. - кандидат технических наук, старший преподаватель
E-mail: tagirkx@yandex.ru
Фатыхова Р.Р. - студент
E-mail: regishka_f93@mail.ru
Казанский государственный архитектурно-строительный университет
Адрес организации: 420043, Россия, г. Казань, ул. Зеленая, д. 1
Об устойчивости упругопластической цилиндрической оболочки при продольном ударе
Аннотация
Рассматривается задача устойчивости упругопластической цилиндрической оболочки при продольном ударе. Известно, что для металлической оболочки при определенных соотношениях толщины к радиусу потеря устойчивости при продольном сжатии происходит за пределом упругости. Для двух видов нагружения оболочки получены критические напряжения с учетом конечности скорости распространения продольной волны. Задача решена в осесимметричной постановке при первом прохождении продольной волны сжатия вдоль оболочки. Используется статический критерий устойчивости. Показано, что при продольном ударе значения критических напряжений более чувствительны к характеру напряженного состояния за фронтом продольной волны, чем к учету конечности скорости распространения последней.
Ключевые слова: цилиндрическая оболочка, потеря устойчивости, продольная волна, критическое напряжение.
В настоящее время вопросам динамического поведения стержней, пластин и оболочек при ударных и быстро возрастающих нагрузках посвящено большое количество исследований. Обзор некоторых из них приводится в [1, 2], где в основном изучалось процесс возрастания прогибов во времени. К задачам, которые решаются в данной статье, непосредственное отношение имеют работы [3, 4], где рассматривалась устойчивость стержня и цилиндрической оболочки при продольном ударе в упругопластической постановке с учетом распространения продольной волны. В работе [5] изучалась устойчивость упругопластических стержней при различных граничных условиях. Здесь также определяются критические напряжения при первом прохождении продольной волны сжатия вдоль оболочки. Задача в упругой постановке рассматривалась в [6, 7]. Относительно методики решения и критерия потери устойчивости сохраняются те же предположения, что в вышеуказанных работах.
Постановка задачи
Рассмотрим цилиндрическую оболочку длиной l, радиусом R и толщиной h, на торце (незакрепленном в продольном направлении) которой прикладывается напряжение о0, превышающее предел пропорциональности материала. Диаграмма сжатия материала с линейным упрочнением показана (рис. 1), где E, Ek, Ec, - упругий, касательный и секущий модули соответственно. При этом рассмотрим два вида нагружения торца оболочки:
1) напряжение о0 прикладывается мгновенно (рис. 2);
2) напряжение линейно возрастает во времени от нуля до максимального значения о0 (рис. 3).
iT
£
а
ь-
f
0
V
т
t
о г, ь
Рис. 1. Диаграмма Рис. 2. «Мгновенное»
сжатия материала приложение напряжения
с линейным упрочнением на ударяемом торце
Рис. 3. Напряжение на ударяемом торце возрастает по линейному закону
Напряженное состояние оболочки по длине известно [8]. В первом случае от ударяемого торца распространяются упругая волна со скоростью а0 Е / р с напряжением за фронтом волны, равным пределу пропорциональности о., и пластическая волна со скоростью а1 =^Ек / р с напряжением о0 за фронтом волны (рис. 4а). Перед фронтом продольной волны невозмущенный участок стержня с нулевым уровнем напряжений. Во втором случае в начальный момент времени от ударяемого торца оболочки распространяется упругая волна с нулевым уровнем напряжений и со скоростью а0 (рис. 4б). Вслед за ней в промежутке времени 0<4<т1, начинают распространяться волны, несущие на своем фронте напряжение 0<о(0<о. (где т1 - момент времени, при котором напряжения на торце о(0=о.). В момент ^=т1 одновременно начинают распространяться две волны с разными скоростями а0 и а1. По длине оболочки образуется участок с постоянным уровнем напряжений о.. При напряжениях на торце оболочки, сколь угодно мало превосходящих напряжение о., распространяются пластические волны со скоростью а1. Будем изучать устойчивость оболочки вплоть до т2, когда напряжения достигают на торце величины о(0 = о0. Таким образом, в отличие от первого вида нагружения, по длине оболочка разбивается на четыре участка с различным уровнем напряжений (рис. 5б). Если устремить т1, т2^0, то приходим к картине напряжений, показанной на рис. 5а. В дальнейшем предполагается, что потеря устойчивости происходит при первом прохождении волны вдоль оболочки.
Рис. 4. Распространение продольных волн вдоль оболочки: а) при мгновенном приложении напряжения; б) при линейном возрастании напряжения
а о* 1 о Я.
и и 1 1 1 : 1 1, Ь и Л
ь Ъ
1 1
а) б)
Рис. 5. Напряженное состояние оболочки по длине: а) при мгновенном приложении напряжения; б) при линейном возрастании напряжения на ударяемом торце
Эпюры продольных напряжений по длине оболочки, соответствующие двум видам нагружения, показаны на рис. 5. На рис. 5а: Ь - полная длина сжатой части оболочки, 11 -длина участка, где напряжения равны о0, 12 - длина участка, где напряжения равны пределу пропорциональности материала о.. Между длинами участков имеют место соотношения:
Ь = а^, ¡1 = а/, Ь = ¡1 +/2, 11 = Ь^ = Ь ^, /2 = Ь
1 Е.
(1)
Выражения для длин по рис. 3-5б запишутся в виде:
/1 = а1(т2 -т1); /2 = (а0 -а1)(т2 -Т1); /3 = а0т1; Ь = /1 + /2 + /3 = а0т2. (2)
Учитывая, что г1=г2о6./о0 (рис. 3), выражения (2) можно представить в иной форме:
/1=VII1 - -
/2 = Ь
1 -ё1
1 - ^ |; /3 = Ь ^.
(3)
Решение
Уравнение возмущенного движения оболочки в осесимметричной постановке запишем в виде [9]:
52 Un ГЛ\д 2w I д\( В2 ,)dW 1 Щ ,d2w —-\ (D2 - D3)—- 1+—I \ — + s h I— 1+—2w + ph—- = 0.
дх2 Г 2 з! дх2 I дх Ц 2R ' )дх I 4R dt2
(4)
Здесь все обозначения общепринятые. При выводе уравнения (4) использована теория деформаций [10]. Коэффициенты В1, В2, 02, 03, характеризующие жесткость оболочки при пластических деформациях, зависят от параметра ¿0 - относительной координаты, определяющую толщину слоя разгрузки. В монографии [9] показано, что «...в достаточно большом числе случаев потери устойчивости оболочек, с которыми приходится иметь дело в инженерной практике, .наименьшие значения критических нагрузок будут реализовываться при г0 = 1, т.е. при исчезающей толщине слоя разгрузки».
1. Рассмотрим первый вид нагружения. При первом прохождении волны вдоль оболочки первая и вторая зоны оболочки будут изгибаться как пластическая среда, а третья зона - как упругая [3]. Запишем выражения для коэффициентов В1, В2, Ю2, 03:
D2 - D3 =
E„ + 3E,
36 E + 3Ek 36 Eh3
12(1 - m2)
h3 при 0 < x < l1
h при l1 < x < b,
при b < x < l.
(5)
где ц - коэффициент Пуассона.
n. í 2R -sh 1=
4 h2
3Eh — (Ec - Ek) + soh при 0 < x < l1
4 h2
3Eh, ~^(E -Ek) + ssh при l1 < x < b,
(6)
- Eh, 3
0
при b < x < l.
Пусть оболочка на торцах шарнирно оперта. Решение уравнения (4) ищем в виде:
кк
w(x,t) = ewtX fk sin(lkx), l k= —, (k = 1,2,3...). k =1 l Коэффициенты (5) и (6) также представим в виде рядов:
¥
(D2 - D3 ) = X dk C0S (Я kx)
BR+sh|=¿ bkcos(i kx),
2R 0 k=0
—R2=x ckcos(i kx),
4 R k=0
(7)
(8)
(9)
(10)
где
36l
3E
(Ec -E)l1 +(3Ek + E)b + —^((-b)
1 m
18n k
(Ec - E)sin (ll) +
3Ek + E -
3E
1-У
sin (l kb)
b0 = l
b = — k nk
— (Ec - E) + s 0 - Ss 111 +I 4R (E - Ek ) + Ss | b
4R(Ec - E) + S0 - ss 0sin (ll ) + ^—(E - Ek) + ss |sin (lkb) 2h
c0 = [E(l -h ) + El ], ck = ~ñk~k (Ec - E)sin(ll).
Подставляя (7)-(10) в (4), после некоторых преобразований получим бесконечный ряд, равный нулю:
к=1
к-1
/к °к+х -/а-+х ./а
г=к+1 г=1
2
= 0,
(11)
где акк=&к=Ь №с^2к)~Лк (2Ь()+Ь2к)+2с(ГС2к-2рИю при г=к, а1к=аы=Х2Хк2(ф_к-ф+к)-Х1Хк(Ь1_к+Ь1+к)+с1_к-с1+к при \фк.
Условие существования отличного от нулевого решения уравнения (11) приводит к
бесконечной системе однородных алгебраических уравнений для неизвестных /к:
¥ к-1
/к Ок + X 1,ак, +Х 1,ак, = 0, (к = 1,2,3...). (12)
г=к+1 г=1
Согласно статическому критерию устойчивости [3-7] (равенство нулю частоты колебаний) критический параметр определяется как наименьший корень уравнения:
ае1|\ай\\ = 0. (13)
2. Теперь рассмотрим второй вид нагружения. В этом случае оболочка по длине разбивается на четыре зоны с разными значениями напряжений в каждой зоне (рис. 5б). Коэффициенты Бь В2, Б2, Б3 имеют вид:
(А - А), В =
Е + 3ЕП3, 36 4 Еск, 3 с при 0 < х < /1
Е + ^ 36 4 Ек 3 при /1 < х < /1 + /2,
Е-12(1 - у2)' 4 Ек 3 при /1 + /2 < х < Ь,
Е-[ 12(1 - у2)' 4 Ек 3 при Ь < х < /.
В
— + а к = 2Я '
-1(Ес -Ек) + о- + (а0 ? а 1) (К - х)к при 0 < х < /1, 41 /
к1
—(Е - Ек) + а -41 1
(Ь - х) а, ——-к
при /1 < х < /1 + /2,
при /1 + /2 < х < Ь, при Ь < х < /.
Коэффициенты рядов (8)-(10) примут вид:
36/
(Ес - Е) +
3Ек + Е -
3Е
1-У
(Ь -13)+ Т
3Е
У2
18п к
(Ес - Е)81п(1 к/1) +
^ 3Е ^ 3Ек + Е - 3Е
1 - У2
81П (Я к (Ь -13 ))
Ь0 = — 0 2/
—(Ес - Е) + а0 - а 1 1/ +1 —(Е - Ек) + а^ |(Ь - /3) + а^
Ь = — к пк
41
(Ес -Е)81п(Як/1) +—(Е -Ек)п(Я к (Ь -/3)) +
Я ,Ь\ Ек
(1 - ос8 (Я Д)) + Оь ( (Я к (Ь - /3))-(Я кЬ ))
2к
с° = 1ё1 [Е/ + (Ес - Е)/1 ], ск = (Ес -Е)8Ш(Як/1).
(14)
(15)
Результаты расчета
Для примера рассмотрена оболочка из дюралюминия Д16Т с данными Е=7,5-106 кг/см2, Ек=3,72-106 кг/см2, ох=2000 кг/см2 [11]. В табл. 1-2 приведены критические
3
напряжения о0 зависимости от параметра ц=Ь/1 (п - отношение длины сжатого участка к полной длине оболочки) при некоторых значениях R/h и l/R. В табл. 1 даны результаты расчета для первого случая нагружения оболочки, в табл. 2 - для второго. Для обеспечения точности вычислений критические напряжения находились из определителей двухсотого порядка. В скобках даны значения о0/ост, где ост -критическое напряжение при статической потере устойчивости [9]:
=hE R c
1 + ^--
1
(16)
Таблица 1
l R 9h\ 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1
2 100 3215,67 (1,155) 3052 (1,059) 2984,25 (1,035) 2950,77 (1,024) 2932,35 (1,017) 2921,08 (1,013) 2913,38 (1,011) 2907,64 (1,009) 2903,17 (1,007)
4 2984,25 (1,035) 2932,36 (1,017) 2913,38 (1,011) 2903,17 (1,007) 2897,05 (1,005) 2893,52 (1,004) 2891,19 (1,003) 2889,41 (1,002) 2888,15 (1,002)
2 50 6898,06 (1,196) 6460,75 (1,121) 6164,7 (1,069) 6036 (1,047) 5978,35 (1,037) 5913,62 (1,026) 5890,85 (1,022) 5863,78 (1,017) 5845,43 (1,014)
4 6164,7 (1,069) 5978,37 (1,037) 5890,85 (1,022) 5845,43 (1,014) 5822,06 (1,01) 5809,06 (1,008) 5800 (1,006) 5792,7 (1,005) 5787,23 (1,004)
Таблица 2
l R X 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1
2 100 4595,17 (1,59) 4128,27 (1,43) 3847,44 (1,33) 3702,46 (1,28) 3599,95 (1,25) 3517,64 (1,22) 3457,1 (1,2) 3410,82 (1,18) 3371,26 (1,17)
4 3847,44 (1,33) 3599,95 (1,25) 3457,1 (1,2) 3371,26 (1,17) 3310,26 (1,148) 3265,25 (1,133) 3230,4 (1,121) 3202,37 (1,111) 3179,37 (1,103)
2 50 10410,3 (1,806) 9048,83 (1,569) 8338,96 (1,446) 7923,82 (1,374) 7631,64 (1,324) 7421,21 (1,287) 7260,08 (1,259) 7132 (1,237) 7027,86 (1,219)
4 8338,96 (1,446) 7631,64 (1,324) 7260,08 (1,259) 7027,86 (1,219) 6867,54 (1,191) 6749,28 (1,171) 6657,92 (1,154) 6584,91 (1,142) 6525,02 (1,132)
По результатам расчета можно сделать следующие выводы:
1. Из табл. 1 видно, что при «мгновенном» нагружении оболочки критические напряжения незначительно превышают статическое критическое значение и практически не зависят от длины сжатой части и общей длины. При этом продольная волна должна пройти расстояние в несколько статических полуволн потери устойчивости. Для короткой оболочки с длиной, сравнимой с длиной статической полуволны, критические напряжения сильно зависят от длины сжатой части и в несколько раз могут превышать статическое критическое значение. Этот вывод полностью согласуется с результатами работ [5, 6], где рассматривалась упругая оболочка. В этих же работах делается вывод, что граничные условия не влияют на значения критических напряжений для оболочек средней длины. Таким образом, для определения критической нагрузки при «мгновенном» нагружении оболочек средней длины с некоторой погрешностью можно использовать формулу (16).
2. Расчеты для второго случая нагружения оболочки (табл. 2) показали, что критические напряжения всегда получаются больше, чем для «мгновенного» нагружения при тех же исходных данных. С уменьшением длины сжатой части и общей длины оболочки значения о0 возрастают. Таким образом, при продольном ударе значения критических напряжений цилиндрической оболочки более чувствительны к характеру напряженного состояния за фронтом продольной волны, чем к учету конечности скорости распространения последней.
Список библиографических ссылок
1. Вольмир А.С. Нелинейная динамика пластинок и оболочек. - М.: Наука, 1972. -432 с.
2. Абросимов Н.А., Баженов В.Г. Нелинейные задачи динамики композитных конструкций. - Н. Новгород: Изд-во ННГУ, 2002. - 400 с.
3. Аракелян А.Е., Мовсисян Л.А. Об упругопластической устойчивости стержня и цилиндрической оболочки при продольных интенсивных нагружениях // Известия АН Армянской ССР. Сер. физ.-мат. наук, 1987, т. 40, № 5. - С. 3-13.
4. Мовсисян Л.А. К устойчивости упругопластических стержней при ударных нагрузках // Известия АН Армянской ССР. Сер. физ.-мат. наук, 1986, т. 39, № 2. -С. 15-23.
5. Хамитов Т.К., Шигабутдинов Ф.Г. Об устойчивости упругопластических стержней при продольном ударе // Аналитическая механика, устойчивость и управление движением. Труды X Международной Четаевской конференции. Т.2. Секция 2. Устойчивость. - Казань, 12-16 июня 2012. - С. 560-569.
6. Хамитов Т.К., Шигабутдинов Ф.Г. К вопросу о потере устойчивости упругих цилиндрических оболочек при продольном сжатии усилием ударного типа // Вестник Иркутского государственного технического университета. - Иркутск, 2010, № 3 (43). -С. 70-76.
7. Шигабутдинов Ф.Г., Хамитов Т.К. Определение критических усилий потери устойчивости упругих цилиндрических оболочек при продольном сжатии силами ударного типа // Вестник КГТУ им. А.Н. Туполева, 2011, № 2. - С. 85-92.
8. Рахматуллин Х.А., Демьянов Ю.А. Прочность при интенсивных кратковременных нагрузках. - М.: Физматгиз, 1961. - 369 с.
9. Королев В.И. Упругопластические деформации оболочек. - М.: Машиностроение, 1971. - 304 с.
10. Ильюшин А.А. Пластичность. - М.: Гостехиздат, 1948. - 376 с.
11. Вольмир А.С. Устойчивость упругих систем. - М.: Физматлит, 1963. - 880 с.
Khamitov T.K. - candidate of technical sciences, senior lecturer E-mail: tagirkx@yandex.ru Fatykhova R.R. - student E-mail: regishka_f93@mail.ru
Kazan State University of Architecture and Engineering
The organization address: 420043, Russia, Kazan, Zelenaya st., 1
On stability of elastic-plastic cylindrical shell under longitudinal impact Resume
We consider the problem of stability of elastic-plastic cylindrical shell under longitudinal impact. It is known that for a metal shell with certain ratios of thickness to the radius of the buckling under longitudinal compression occurs beyond the elastic limit. Elastic-plastic material properties are characterized by a diagram with linear hardening, we use the theory of deformation without regard to unloading. For the two types of loading of the shell critical load obtained taking into account the finite speed of propagation of longitudinal waves. Under a critical load refers to the load at which the frequency of the bending oscillations is equal to zero. The problem is solved in the axisymmetric formulation of the first passage of the longitudinal compression waves along the shell. It is shown that if the longitudinal impact values of critical loads are more sensitive to the nature of the stress state behind the front of the longitudinal wave than to the consideration of the finite speed of propagation of the latter.
Keywords: cylindrical shell, buckling, longitudinal wave, critical voltage.
Reference list
1. Vol'mir A.S. Nonlinear dynamics of plates and shells. - M.: Nauka, 1972. - 432 p.
2. Abrosimov N.A., Bazhenov V.G. Nonlinear problems of dynamics of composite structures. - N. Novgorod: publishing house of NNSU, 2002. - 400 p.
3. Arakelyan A.E., Movsisyan L.A. On elastic-plastic stability of rod and cylindrical shell under longitudinal intensive loads // Proceedings of the Academy of Sciences of the Armenian SSR. Series of physico-mathematical Sciences. 1987, vol. 40, № 5. - P. 3-13.
4. Movsisyan L.A. On the stability of elastic-plastic bars for impact loads // Proceedings of the Academy of Sciences of the Armenian SSR. Series of physico-mathematical sciences. 1986, vol. 39, № 2. - P. 15-23.
5. Khamitov T.K., Shigabutdinov F.G. On stability of elastic-plastic rods under longitudinal impact // Analytical mechanics, stability and control. 10th International Chetaev Conference. Vol. 2. Section 2. Stability. - Kazan, June 12-16, 2012. - P. 560-569.
6. Khamitov T.K., Shigabutdinov F.G. On the buckling of elastic cylindrical shells under longitudinal compression force impact type // Vestnik Irkutskogo gosudarstvennogo tekhnicheskogo universiteta. Irkutsk, 2010, № 3 (43). - P. 70-76.
7. Shigabutdinov F.G., Khamitov T.K. Because the determination of the critical efforts of the buckling of elastic cylindrical shells under longitudinal compression forces impact type // Vestnik KSTU n.a. A.N. Tupolev, 2011, № 2. - P. 85-92.
8. Rakhmatullin K.A., Demyanov Y.A. Strength under the intensive short-time loads. - M.: Fizmatgiz, 1961. - 369 p.
9. Korolev V.I. Elastic-plastic deformation of shells. - M.: Mashinostroenie, 1971. - 304 p.
10. Ilyushin A.A. Plasticity. - M.: Gostehizdat, 1948. - 376 p.
11. Vol'mir A.S. Stability of elastic systems. - M.: Fizmatlit, 1963. - 880 p.
