CC BY
42
УДК 539.3:534.1 В01: 10.18287/2541-7533-2016-15-3-209-222
ВЛИЯНИЕ АСИММЕТРИЧНЫХ НАЧАЛЬНЫХ НЕСОВЕРШЕНСТВ ФОРМЫ НА СВОБОДНЫЕ КОЛЕБАНИЯ ТОНКИХ ОБОЛОЧЕК
© 2016
С. В. Серёгин кандидат технических наук, старший научный сотрудник управления научно-исследовательской деятельностью, Комсомольский-на-Амуре
государственный технический университет, Seregin-komsHome@yandex.ru
В рамках теории пологих оболочек изучается влияние начальных отклонений от идеальной круговой цилиндрической формы на изгибные колебания тонких изотропных оболочек. В работе впервые показана возможность возникновения дополнительной зоны расщепления изгибного частотного спектра, обусловленная наличием у оболочки начальных несовершенств формы. Установлено, что расщепление частотного спектра имеет место не только в случаях, когда число волн окружных динамических деформаций равно числу волн несовершенств формы оболочки, как это принято считать в настоящее время, но и в случаях, когда число формообразующих волн в два раза меньше волн несовершенств формы. В первом случае обе расщеплённые частоты ниже соответствующих частот колебаний идеальной оболочки. Установлено, что при некоторых геометрических параметрах оболочки частоты радиальных колебаний могут быть соизмеримыми с изгибными. Решения, учитывающие взаимодействие сопряжённых изгибных и радиальных форм колебаний, дополнены новыми результатами и выводами. Сделан вывод, что при действии периодических нагрузок, неизбежно возникающих в условиях эксплуатации оболочечных конструкций, нерезонансные зоны, определённые согласно традиционному подходу, по существу могут оказаться резонансными. Резонансно опасные режимы колебаний оболочечных конструкций с неправильностями такого рода могут возникать на частотах, намного меньших или больших, чем это предсказывает традиционная теория. Такое поведение несовершенных оболочек может привести к высокому напряжённо-деформируемому состоянию конструкции и сложным видам динамической неустойчивости в условиях эксплуатации. Полученные в этой работе результаты и выводы требуют уточнения уже решённых ранее задач динамики оболочек с несовершенствами формы.
Оболочка, начальные неправильности, несовершенства, формы, изгибные, радиальные колебания, расщепление, частотный спектр, параметр волнообразования, резонанс.
Цитирование: Серёгин С.В. Влияние асимметричных начальных несовершенств формы на свободные колебания тонких оболочек // Вестник Самарского университета. Аэрокосмическая техника, технологии и машиностроение. 2016. Т. 15, № 3. С. 209-222. Б01: 10.18287/2541-7533-2016-15-3-209-222
Введение
Тонкие круговые цилиндрические оболочки широко используются в качестве основных несущих элементов, применяемых в авиастроении, кораблестроении, ракетостроении, космической технике. Такие конструкции имеют неизбежные начальные отклонения от идеальной круговой формы, обусловленные изготовлением их на производстве или конструктивными особенностями при проектировании судов, ракет и космических аппаратов.
В условиях эксплуатации оболочки подвергаются интенсивным динамическим нагрузкам, в результате которых могут происходить сложные, нестационарные процессы и специфические явления, приводящие к высоким напряжённо-деформируемым состояниям (НДС) и динамической неустойчивости конструкций [1].
Большая часть работ в научной литературе посвящена изучению устойчивости геометрически несовершенных оболочек [2-5], определению их НДС [6-8] и несущей способности [9-12]. Меньшая их часть посвящена изучению динамических процессов в тонких оболочках [13-15]. Краткий обзор работ по расчёту оболочек с несовершенствами формы представлен в [16-19].
В настоящее время принято считать, что осесимметричная часть динамического прогиба оказывает существенное влияние только на нелинейное поведение оболочек, а при колебаниях с малыми амплитудами, при решении таких задач, в математической модели достаточно учитывать лишь сопряжённые изгибные формы [18-20].
Однако в работе [21] на примере более простой (предельной) задачи колебаний бесконечно длинной круговой цилиндрической оболочки - кольца, находящегося в условиях плоской деформации, показано, что уже в линейной постановке учёт радиальных колебаний приводит как к качественно, так и количественно лучшим результатам, чем взаимодействие лишь сопряжённых изгибных форм. Результаты этой работы качественно согласуются с численными исследованиями [22; 23]. Предложенная математическая модель получила свое распространение в задачах динамики оболочек конечной длины, несущих присоединённую массу, ослабленных отверстием, в которых позволила устранить ряд качественных противоречий и установить новые эффекты и особенности взаимодействия сопряжённых изгибных форм с радиальными [24-26]. Результаты этих работ хорошо согласуются с численными исследованиями [27-30], в которых также обозначены противоречия и погрешности аналитических решений.
Настоящая работа является продолжением исследования [21; 22] в том плане, что полученные результаты для предельной задачи - колебаний изолированного, геометрически несовершенного кольца, обобщаются на случай колебаний оболочки уже конечной длины, имеющей начальные отклонения от идеальной круговой формы.
Математическая модель
Пусть шарнирно опёртая по торцам круговая цилиндрическая оболочка массой М0, радиусом Я, длиной ' и толщиной И совершает малые изгибные колебания вблизи зоны основного резонанса.
Линеаризованные уравнения движения. Анализ основывается на известных уравнениях теории пологих оболочек, которые для изотропной несовершенной оболочки имеют вид [32]:
Б ^4 д2wn д2Ф 1 д2Ф д2я
—V V =-0--1----р-
И дх2 ду2 Я дх2 дг2 '
^^ д2я0 д2я 1 д2я
—V 4Ф =--0----
Е дх2 ду2 Я дх2
Здесь Б = ЕИ3/(12(1 - /2 )) цилиндрическая жёсткость, где Е - модуль Юнга; /л -коэффициент Пуассона; V4 - бигармонический оператор Лапласа; х, у, г) - динамический прогиб; Ф(х, у,г) функция напряжений; р- массовая плотность; г - время.
Граничные условия и условие периодичности решения. Как и в подавляющем большинстве случаев, задача решается при удовлетворении тангенциальным граничным условиям Ы1 = Т = 0 «в среднем» [19; 32]:
И с д 2Ф Иг д 2Ф
0 аИ * =0; - п 0 1дхду у =0 при х=а х=' (2)
Условие непрерывности окружного перемещения имеет вид:
2nR ~ 2nR
f-^ = f 0 ^ 0
Сд2Ф д2Ф^
E
Kdx2 М dy2 J R dy dy
w dw0 dw +--- 0
dy = 0. (3)
Конечномерная модель. Для идеальной оболочки каждая точка окружности может стать узлом колебаний. Экспериментально установлено, что несовершенства формы исключают эту неопределённость и приводят к фиксации узлов в окружном направлении. Положение узлов определяется формой начальных несовершенств [17; 19]. Поэтому принято считать, что динамический прогиб несовершенной оболочки в линейной постановке может быть аппроксимирован выражением, учитывающим только взаимодействие сопряжённых sinPy sinax, cosPy sinax (a = mn/1, P = n¡R, n - число окружных волн динамических деформаций, m - число продольных полуволн, n, m = 1,2,3..., N) изгибных форм, которые являются формами собственных колебаний идеальной оболочки и соответствуют одному и тому же волновому параметру n . Этим формам, как известно, соответствует одна и та же собственная частота [19; 32]. Такой подход позволяет установить качественные особенности динамического деформирования идеальных оболочек, обнаруженных в экспериментах. Однако, для реальных (имеющих несовершенства формы) оболочек традиционная математическая модель приводит к результатам, не всегда согласующимся с экспериментальными данными [18; 19].
В настоящей работе математическая модель, предложенная для более простой (предельной) задачи колебаний бесконечно длинной круговой цилиндрической оболочки - кольца, находящегося в условиях плоской деформации [21], обобщается на случай колебаний геометрически несовершенной оболочки конечной длины.
Уточнение конечномерной модели предполагает, что динамическая асимметрия уже в линейной постановке приводит не только к связанности сопряжённых изгибных форм, но и к взаимодействию низкочастотных изгибных колебаний оболочки с высокочастотными радиальными колебаниями. Динамический прогиб оболочки в новом подходе имеет вид [21, 24-26]:
1
N
w(x, y, t) = 2 (fr sin Py + fn,m cos Py + fm )sin ax, (4)
n,m
в котором ,т ^) и /2",т ^) - обобщённые координаты; дополнительно введённая к
[19] обобщённая координата /3т (1) отвечает радиальным колебаниям.
В [17] авторы используют аналогичный подход к построению конечномерной модели, но для оболочки, имеющей отклонения в окружном направлении (несовершенства распределены не по всей длине оболочки). Отличия в том, что в (4) координата
/т ^) , описывающая радиальные формы колебаний, зависит от количества продольных полуволн т, а в [17] - от количества продольных полуволн т и от количества волн в окружном направлении " , что, по мнению автора, противоречит физическим представлениям о динамическом поведении реальной оболочки. Для рассматриваемой в этой статье задачи - оболочки, имеющей отклонения, соответствующие характеру ее волнообразования, авторы работы [17] динамический прогиб аппроксимировали более сложным выражением. В связи с этим некоторые результаты в данной статье будут несколько схожи с теми, что получены в [17], однако они будут уточнены и дополнены новыми выводами.
Начальные несовершенства. Известно, что наибольшее влияние на статические [1-12] и динамические [13-16; 17-19; 32] характеристики оболочек оказывают те начальные несовершенства, которые имеют такую же форму, как и прогибы, возника-
ющие при колебаниях или потере устойчивости оболочек. Поэтому с теоретической и практической сторон этот случай наиболее интересен.
Пусть оболочка имеет начальные отклонения от идеальной круговой формы, изменяющиеся по закону:
я0 (х У) = /0 ^п (ДУ + ф0) ^п ах = (/ю ^п РУ + /20 с°8 РУ) ^п ах,
(5)
где /10 = /0соъу0; /20 = /0^т^0; /0 - амплитуда начальных отклонений; (р0 - начальный угол отсчёта (рис. 1).
Рис. 1. Схематическое изображение рассматриваемых несовершенств оболочки
Собственные частоты и формы колебаний геометрически несовершенных оболочек
Подстановка функции динамического прогиба (4) и функции начальных отклонений от идеальной круговой формы (5) в первое уравнение линеаризованной системы (1) приводит к неоднородному дифференциальному уравнению, из которого определяются функция напряжений в срединной поверхности:
Ф = Е[Ф0 б1П(ох)вт(вУ)+ Ф1 вт(ах)СОБ^У) + Ф2 §т(/У) + Ф3 сов(в) +
+ Ф4 СОБ(2ОХ)вт(ву) + Ф5 СОБ(2ОХ)сов(вУ) + Ф6 Б1П(ОХ) + Ф7 СОБ(2ОХ)+ (6) + Ф8 ып(2ву) + Ф9 соБ(2вУ) + Ф01 х72].
Коэффициенты в (6) определяются по формулам:
Ф 0 =
Иа
Иа
! 2 2 Иа
Я (а2 +в / Я(а2 +в ) 2в
+ ,
22 Иа
Ф 3 ^20^3, Ф 4 =
2Р
„2„.2 о2
+
И 2а2в2 2(4а2 +в2 )
а10а3,
И2а£р£ И И2в2 ( \
Ф 5 = -Та ^^ Ф 6 = а3, Ф 7 = 7Т"Г (а10а1 + а20а2 А
о1/1„,2 , о2т а Я 16а
(4а2 +в2 )2
24а +
И2а2
а2 Я
16а2
22 Иа
. пи,/ ч _ Л/ СС /- \
Ф8 = (а10а2 + а20а1 А Ф 9 = (а10а1 - а20а2 ),
16в
16в
Ф 01 определяется из условия непрерывности окружного перемещения у(х, у, г) при удовлетворении краевым условиям (2): Ф 01 = в2 (/10 /1 + /20 /2 )/4.
Подстановка (4), (5), (6) и ортогонализация уравнения движения (1) в первом приближении приводит к системе трёх связанных динамических уравнений относительно безразмерных координат аг — /1/И (г — 1,2,3):
^ъ спа1 С12аша20а2 — 0,
а 2 + С21а10а20а1 + С22а2 + С23а20а3 — 0, а3 + С31а10а1 + С32а20а2 + С33а3 — 0,
(7)
где точками обозначено дифференцирование по безразмерному времени т — Л (Л-собственная частота), а коэффициенты равны:
С11 = 1 +
8®
-[(3 + в4 ) а20 + в4 а^0
3а
8®
С22 = 1 +
8®
-[(3 + в4 ) а^0 + в4 <
2
а —
К У
в—ПК.
п1
С33 — 2 ®
2 ав
Р +-
8
2 +
(1 + 462 )2
(а10 + а20 )
8а
12
С13 — С23 —
3п®
1 +
в'
(1 + 462 У
4а
12
С31 — С32 —
3ж®'
1 +
в4
(1 + 62 )2
Квадраты безразмерных частот, соответствующие изгибным и радиальным формам колебаний круговой цилиндрической оболочки, определяются по формулам:
а(1 + 62 )2
К Л — V— )
п Е 12(1 -,
в4
-; р2 — 1+
ав4
12 (1 -Ц2 )"
5(1 - Ц ) (1 + в2 )2
В традиционном решении система модальных уравнений (7) имеет только первые два уравнения и а3 — 0. Безразмерные расщеплённые собственные частоты изгиб-ных колебаний Оп — ®пг/ ®п (г —1,2;®пг - частоты колебаний оболочки с динамической асимметрией) равны [19, 32]:
О П, — 1 + а6 </!
О2 — 1 + ^0, "п 2 „
(в' + 3).
(8)
Таким образом, начальные несовершенства вида (5), соответствующие характеру волнообразования оболочки, расщепляют изгибный частотный спектр. При этом частоты колебаний Оп1 и Оп2 всегда больше единицы. Безразмерные частоты оболочки, имеющей идеальную геометрическую поверхность, равны О п1 — О п 2 — 1.
В новом решении система динамических уравнений определяет три частоты О пг (г — 1,2,3 ). Первым двум частотам Оп1 и Оп2 соответствуют преимущественно изгиб-ные формы колебаний, а третьей частоте О п3 - преимущественно радиальные формы. Большая из расщеплённых собственных частот Оп1 , как и в традиционном решении
[19], определяется точно, её квадрат равен О^ — 1 + ав4а120 /8. Меньшая из расщеплён
1
1
ных собственных частот, в отличие от традиционного решения, ниже частоты колебаний идеальной оболочки.
Значения частот, полученные при решении задачи, учитывающей взаимодействие сопряжённых изгибных форм (традиционное решение) при различных параметрах а, в и а10, представлены на рис. 2, а; значения частот, полученные при решении задачи, учитывающей взаимодействие сопряжённых изгибных и радиальных форм, - на рис. 2, б. На представленных графиках две верхние линии соответствуют частотам Оп2, две нижние - О п1.
1.1
1.05
а
0.95
Е = = 0.5
= 0.8 V г!--
"Л1 е в = 02 У = 0.3 /
а.
1.05
0.95
0.9
0.85
0.2 0.4 0.6 0.8 1 0.8
■^г- > \ V £■ = 02 ? в = 0,3 V г = 0.5
^ / 1е=оя
"«1
0 0.2 0.4 0.6 0.8
а б
Рис. 2. Зависимость расщеплённых частот Оп1 и Оп2 от амплитуды несовершенств а10
Видим, что в и в новом, и в традиционном решениях эффект расщепления изгиб-ного частотного спектра зависит от геометрических и волновых параметров оболочки. Вторая частота нового решения равна первой частоте традиционного решения. Одна из расщеплённых собственных частот выше соответствующей частоты оболочки с идеальной круговой формой, другая - ниже, как и отмечается в работе [17]. В традиционном решении обе расщеплённые частоты больше частоты колебаний идеальной оболочки. С увеличением амплитуды несовершенств в обоих решениях расщепление спектра Оп2 - Оп1 увеличивается.
Следует отметить, что частота, равная Оп1 = ^ 1 + ав4а^0 / 8 , для оболочек, размеры которых соответствуют реальным конструкциям, применяемым в промышленности (а < 0,5; в < 0,8 и а10 = 1), практически равна единице, то есть частоте колебаний идеальной оболочки.
1.4
1.3
О
«з
1.2
1.1
е = 02 а10=2
Е =0, 5^ = 1
я = 0 2 = 1
0 0.2 0.4 0.6 0.8
в
Рис. 3. Зависимость частоты преимущественно радиальных колебаний Оп2 от параметра волнообразования в
Частоты преимущественно радиальных колебаний представлены на рис. 3. Видим, что частоты радиальных колебаний с ростом е, в и а10 увеличиваются. При определённых значениях е, в и а10 они могут быть соизмеримы с частотами преимущественно изгибных колебаний. Отметим, что о близости частот, соответствующих изги-бным и радиальным формам колебаний, говорится и в задачах, близких к рассматриваемой [24-26].
Численное решение
Рассмотрим оболочку со следующими геометрическими и физическими характеристиками: I / Я = 1; Я/И = 200; Я = 1; Е = 2 -1011 МПа; р = 7800 кг/м3; ¡ = 0,3. Минимальной частоте спектра оболочки с такими параметрами соответствует птЬ = 8 . Пусть оболочка имеет начальные неправильности формы /0 = И, изменяющиеся по закону (5), и число волн несовершенств - п0 = 8. Схематическое изображение такой оболочки представлено на рис. 1, результаты расчёта - на рис. 4. Сплошной линией обозначена большая из расщеплённых безразмерных собственных частот Оп2, штриховой -
меньшая Оп1, пунктирной - безразмерная частота колебаний идеальной оболочки Оп .
1,15
1,1
1,05 ^ 1
0,95 0,9 0,85
? ™ ^
^ / \ 25%
\
Л'....... \
\ \ \ Д = 5%
----
Рис. 4. Зависимость расщеплённых частот преимущественно изгибных колебаний от количества формообразующих волн п
Из рис. 4 видим, что расщепление изгибного частотного спектра возникает не только при п = п0, как это принято считать в настоящее время [17-19; 32-33], но и при
п = 2п0 . В случае, когда п = п0, обе расщеплённые частоты меньше частоты колебаний геометрически совершенной оболочки, что противоречит традиционному решению [19], где частоты больше, и согласуются с новым аналитическим решением и решением, представленным в [17]. Снижение частоты (рис. 4) Оп2 -Оп1 составляет 5%. В аналитическом решении расстройка частотного спектра (в рассматриваемом примере е = 0,102 , в = 0,393 ) Оп2 -Оп1 равна 9%. В случае, когда п = 2п0 , одна из расщеплённых частот Оп2 больше частоты колебаний соответствующей идеальной оболочки, другая Оп1 - меньше. Расстройка спектра может быть довольно существенной и в рассматриваемом примере превосходит известный случай ( п = п0) в пять раз.
Следует отметить, что аналогичный эффект был обнаружен автором на примере колебаний несовершенного кольца, находящегося в условиях плоской деформации [22].
Формы колебаний и соответствующие им расщеплённые собственные частоты для случая п = п0 и п = 2п0 представлены на рис. 5.
Если оболочка имеет идеальную геометрическую поверхность, то обеим сопряжённым изгибным формам для каждых п соответствуют равные частоты: Л4 = 301,5 Гц, Л = 139,25 Гц.
Рис. 5. Формы колебаний и соответствующие им расщепленные частоты: а - при п = п0 = 8: Л81 = 132 Гц, Л82 = 139,2 Гц; б - при п0 = п /2 = 4 : Л41 = 265,27Гц, Л42 = 348,7 Гц
Видим, что при п = п0 = 8 большая из расщеплённых собственных частот несовершенной оболочки Л82 практически равна частоте колебаний оболочки без ^0(х, У) - Л. Меньшая из расщеплённых собственных частот несовершенной оболочки Л81 ниже соответствующей частоты колебаний оболочки без ^0(х, у) - Л8 более чем на 5%. Расстройка частотного спектра Л82 -Л81 « 5% . В зоне п0 = п /2 = 4 большая из расщеплённых собственных частот колебаний оболочки с ^0(х, у) - Л42 превышает соответствующую частоту идеальной оболочки Л4 на 15%. Меньшая частота Л41 - ниже 12% по сравнению с Л4. Расстройка частотного спектра Л42 — Л41 ~ 27% .
Стоит отметить, что при анализе динамического поведения оболочек с несовершенствами формы наблюдались и более сложные (неоднозначные) формы колебаний [27; 28; 34].
Полученные в этой статье результаты позволяют уточнить раздел механики в области колебаний оболочек и могут потребовать уточнения уже решённых ранее задач динамики оболочек с начальными несовершенствами формы.
В работе обнаружена возможность возникновения зоны расщепления изгибного частотного спектра, обусловленная наличием у оболочки начальных отклонений от идеальной круговой формы. Асимметричные несовершенства формы приводят к расщеплению изгибного частотного спектра. При этом расстройка частотного спектра наблюдается не только при п = п0, как это принято считать в настоящее время, но и при
п = 2п0. В случае, когда п = п0, обе расщеплённые частоты меньше, а не больше частоты колебаний геометрически совершенной оболочки. В случае, когда п = 2п0, одна из расщеплённых частот Оп2 больше частоты колебаний соответствующей идеальной оболочки, другая Оп1 - меньше. Расстройка спектра может быть довольно существенной. Эффект расщепления спектра как при п = п0, так и при п = 2п0 , зависит от геометрических и волновых параметров оболочки. Частоты преимущественно радиальных ко-
а
б
Основные выводы
лебаний при определённых значениях параметров волнообразования и амплитуд несовершенств формы могут быть соизмеримы с частотами преимущественно изгибных колебаний.
Библиографический список
1. Кузнецов А.П., Кузнецов С.П. Нелинейные колебания, катастрофы, бифуркации, хаос // Известия высших учебных заведений. Прикладная нелинейная динамика. 2013. Т. 21, № 3. С. 1-10.
2. Анисимов М.К., Поникаров С.И., Черенков А.В. Численный анализ устойчивости цилиндрических оболочек, нагруженных осевой сжимающей силой, с учётом геометрических несовершенств // Машиностроение и инженерное образование. 2006. № 4. С. 49-56.
3. Бегичев М.М. Численный анализ устойчивости стержневых систем и оболочек при упругих и пластических деформациях с учётом начальных несовершенств. Автореферат дис. ... кандидата технических наук. Москва, 2013. 24 с.
4. Лопаницын Е.А., Матвеев Е.А. Устойчивость цилиндрических оболочек с начальными несовершенствами под действием внешнего давления // Известия Российской академии наук. Механика твёрдого тела. 2011. № 2. С. 16-25.
5. Олевский В.И., Мильцын А.М., Плетин В.В. Многофакторные исследования устойчивости тонкостенных оболочек с несовершенствами // Вестник Нижегородского университета им. Н.И. Лобачевского. 2011. № 4(4). С. 1663-1665.
6. Ефремова Г.И. Определение напряжённо-деформированного состояния подкреплённых цилиндрических оболочек с начальной погибью произвольной формы // Труды Центрального научно-исследовательского института имени академика А Н. Крылова. 2009. № 42. С. 19-26.
7. Попов О.Н., Моисеенко М.О., Трепутнева Т.А. Влияние симметричной общей начальной погиби на напряжённо-деформированное состояние и устойчивость пологих цилиндрических оболочек // Строительная механика и расчёт сооружений. 2010. № 4. С. 34-39.
8. Товстик П.Е., Черняев С.П. Нелинейное деформирование тонких оболочек с учётом несовершенств формы срединной поверхности // Вестник Санкт-Петербургского университета. Серия 1. Математика. Механика. Астрономия. 2004. № 3. С. 88-95.
9. Варяничко М.А., Карасев А.Г., Лихачева О.В., Красовский В.Л. Влияние начальных несовершенств геометрии на критическое давление замкнутых упругих гладких пологих конических оболочек // Вюник Придшпровсько! державно! академп будiвництва та архггектури. 2010. № 6 (147). С. 20-31.
10. Лопашев П.А. Метод автоматизации учёта несовершенств геометрической формы цилиндрических оболочек при расчёте несущей способности по методу конечных элементов // Труды Центрального научно-исследовательского института им. академика А Н. Крылова. 2009. № 46. С. 173-180.
11. Лукин Я.А., Ефименко Е.С. Оптимизация несущих систем кузовов цельнометаллических вагонов типа замкнутой оболочки с учётом начальных несовершенств // Вестник Ростовского государственного университета путей сообщения. 2008. № 1 (29). С. 9-13.
12. Мамай В.И. Несущая способность тонкостенных оболочек с локальными несовершенствами // Известия Российской академии наук. Механика твёрдого тела. 2011. № 2. С. 26-31
13. Коротков А.В., Куликов Ю.А. Вынужденные параметрические колебания криволинейных композитных труб с начальными технологическими несовершенствами // Механика композиционных материалов и конструкций. 2014. Т. 20, № 1. С. 16-33.
14. Мильцын А.М., Олевский В.И., Плетин В.В. О формах закритического волнообразования неоднородно нагруженных цилиндрических оболочек с технологическими несовершенствами // Восточно-Европейский журнал передовых технологий. 2011. Т. 5, № 7 (53). С. 44-48.
15. Снисаренко С.И., Хвацков Б.Е. Динамическое деформирование цилиндрических композитных оболочек с наличием структурных и технологических несовершенств // Научный вестник Новосибирского государственного технического университета. 2009. № 3. С. 115-120.
16. Моисеенко М.О., Песцов Д.Н., Попов О.Н., Трепутнева Т.А. Обзор работ по расчету пластин и оболочек с малыми начальными несовершенствами формы срединной поверхности: депонированная рукопись № 501-В2011; 22.11.2011. 28 с.
17. Тарануха Н.А., Лейзерович Г.С. Динамика «неправильных» оболочек. Владивосток: Дальнаука, 2005. 422 с.
18. Amabili M. Nonlinear vibrations and stability of shells and plates. New York: Cambridge University Press, 2008. 392 p.
19. Кубенко В.Д., Ковальчук П.С., Краснопольская Т.С. Нелинейное взаимодействие форм изгибных колебаний цилиндрических оболочек. Киев: Наукова думка, 1984. 220 с.
20. Варадан Т.К., Пратхап Дж., Рамани Х.В. Нелинейные свободные изгибные колебания тонкостенных круговых цилиндрических оболочек // Аэрокосмическая техника. 1990. № 5. С. 21-24.
21. Серёгин С.В., Лейзерович Г.С. Свободные колебания бесконечно длинной круговой цилиндрической оболочки с начальными неправильностями и малой присоединённой массой // Учёные записки Комсомольского-на-Амуре государственного технического университета. 2014. Т. 1, № 4 (20). С. 36-43.
22. Лейзерович Г.С., Приходько Н.Б., Серёгин С.В. О влиянии малой присоединённой массы на расщепление частотного спектра кругового кольца с начальными неправильностями // Строительная механика и расчёт сооружений. 2013. № 6 (251). С. 49-51.
23. Лейзерович Г.С., Приходько Н.Б., Серёгин С.В. О влиянии малой присоединённой массы на колебания разнотолщинного кругового кольца // Строительство и реконструкция. 2013. № 4(48). С. 38-42.
24. Серёгин С.В. Свободные изгибно-радиальные колебания тонкой круговой цилиндрической оболочки, несущей присоединённую массу // Вестник МГСУ. 2014. № 11. С. 74-81.
25. Серёгин С.В. Численное и аналитическое исследование свободных колебаний круговых цилиндрических оболочек, несущих присоединённую массу, линейно распределённую вдоль образующей // Вычислительная механика сплошных сред. 2014. Т. 7, № 4. С. 378-384.
26. Серегин С.В. Свободные колебания тонкой круговой цилиндрической оболочки, ослабленной отверстием // Известия высших учебных заведений. Авиационная техника. 2015. № 3. С. 9-13.
27. Серёгин С.В. Влияние площади контакта и величины линейно распределённой и сосредоточенной массы с круговой цилиндрической оболочкой на частоты и формы свободных колебаний // Вестник МГСУ. 2014. № 7. С. 64-74.
28. Серёгин С.В. Об эффекте расщепления изгибного частотного спектра тонких круговых цилиндрических оболочек, несущих присоединённую массу // Строительная механика и расчет сооружений. 2014. № 6 (257). С. 59-61.
29. Серёгин С.В., Лейзерович Г.С. Влияние присоединённой массы на динамические характеристики тонкой оболочки // Проблемы машиностроения и автоматизации. 2015. № 4. С. 83-89.
30. Серегин С.В. Исследование динамических характеристик оболочек с отверстиями и присоединённой массой // Вестник МГСУ. 2014. № 4. С. 52-58.
31. Серёгин С.В. Влияние присоединённого тела на частоты и формы свободных колебаний цилиндрических оболочек // Строительная механика и расчёт сооружений. 2014. № 3(254). С. 35-38.
32. Вольмир А.С. Нелинейная динамика пластинок и оболочек. М.: Наука, 1972.
432 с.
33. Гришин Н.Ф. Свободные колебания цилиндрических оболочек, имеющих начальную погибь // В кн.: Динамика судовых оболочек (Материалы по обмену опытом). Вып. 208. Л.: Судостроение, 1974. С. 31-40.
34. Серёгин С.В. Динамика тонких цилиндрических оболочек с присоединённой массой. Комсомольск-на-Амуре: Комсомольский-на-Амуре государственный технический университет, 2016. 175 с.
INFLUENCE OF ASYMMETRICAL INITIAL FORM IMPERFECTIONS
ON FRE OSCILLATIONS OF THIN SHELLS
© 2016
S. V. Seregin Candidate of Science (Engineering), Senior Research Associate, Research Management Department, Komsomolsk-na-Amure State Technical University, Komsomolsk-on-Amur, Russian Federation, Seregin-komsHome@yandex.ru
In the framework of the theory of shallow shells we studied the effect of the initial deviations from the ideal circular cylindrical shape on flexural vibrations of thin isotropic shells. For the first time the possibility of an additional zone of splitting the flexural frequency spectrum due to the presence of initial shell form imperfections is shown in the paper. It is found that the splitting of the frequency spectrum occurs not only in cases when the number of waves of circumferential dynamic strain is equal to the number of waves of shell form imperfections as generally agreed at present, but also in cases when the number of shape-generating waves is half that of wave imperfections. In the first case, both the split frequencies are lower than the corresponding frequency of the ideal shell. It was established that for some geometrical parameters of the shell the frequencies of radial vibrations may be commensurate with the flexural ones. Solutions that take into account the interaction of bending and radial mode shapes are complemented with new results and conclusions. It was concluded that under the influence of periodic loads that inevitably arise in the operating conditions of shell structures, non-resonant zones defned according to the traditional approach can, in fact, be resonant. Resonance hazardous vibration conditions of shell structures with irregularities of this kind can occur at frequencies much lower or higher than the traditional theory predicted. This behavior of imperfect shells may lead to high stress-strain state of the structure and complex types of dynamic instability under operating conditions. The results and conclusions obtained in this work call for the revision of problems of the dynamics of shells with form imperfections that were solved earlier.
Shell, initial irregularities, imperfections, shape, flexural and radial vibrations, splitting, frequency spectrum, wave generation parameter, resonance.
Citation: Seregin S.V. Influence of asymmetrical initial form imperfections on fre oscillations of thin shells. Vestnik of Samara University. Aerospace and Mechanical Engineering. 2016. V. 15, no. 3. P. 209-222. DOI: 10.18287/2541-7533-2016-15-3-209-222
References
1. Kuznetsov A.P., Kuznetsov S.P. Nonlinear oscillations, catastrophe, bifurcation, chaos. Izvestiya VUZ. Applied Nonlinear Dynamics. 2013. V. 21, no. 3. P. 1-10. (In Russ.)
2. Anisimov M.K., Ponikarov S.I., Cherenkov A.V. Numerical analysis of the stability of cylindrical shells loaded with axial compressive force, taking into account the geometric imperfections. Mechanical Engineering and Engineering Education. 2006. No. 4. P. 49-56. (In Russ.)
3. Begichev M.M. Chislennyy analiz ustoychivosti sterzhnevykh sistem i obolochek pri uprugikh i plasticheskikh deformatsiyakh s uchetom nachal'nykh nesovershenstv. Avtoreferat dis. kand. tekhnicheskikh nauk [Numerical analysis of stability of frame structures and shells under elastic and plastic deformations with account for initial imperfections. Extended abstract of Cand. Sci. (Eng.) Diss.]. Moscow, 2013. 24 p.
4. Lopanitsyn E.A., Matveev E.A. Stability of cylindrical shells with initial imperfections under the action of external pressure. Mechanics of Solids. 2011. V. 46, Iss. 2. P. 170-178. DOI: 10.3103/s0025654411020038
5. Olevs'kyy V.I., Mil'tsyn A.M., Pletin V.V. Multifactorial study of the stability of thin-walled shells with imperfections. Vestnik of Lobachevsky University of Nizhni Novgorod. 2011. No. 4 (4). P. 1663-1665. (In Russ.)
6. Efremova G.I. Defining stressed & strained state of stiffened cylindrical shells with initial bending of arbitrary shape. Proceedings of the Krylov State Research Centre. 2009. No. 42. P. 19-26. (In Russ.)
7. Popov O.N., Moiseenko M.O., Treputneva T.A. Influence of symmetric general initial flexure on mode of deformation and stability of flat cylindrical shells. Structural Mechanics and Analysis of Constructions. 2010. No. 4. P. 34-39. (In Russ.)
8. Tovstik P.E., Chernjaev S.P. Nonlinear deformation of thin shells with the middle surface imperfections. Vestnik Sankt-Peterburgskogo Universiteta. Ser. 1. Matematika, Mekhanika, Astronomiya. 2004. No. 3. P. 88-95. (In Russ.)
9. Varyanichko M.A., Karasev A.G., Lykhachova O.V., Krasovsky V.L. The influence of the initial imperfections of geometry on the critical pressure of the closed elastic plain flatsloped conical shells. Bulletin of Prydniprovs 'ka State Academy of Civil Engineering and Architecture. 2010. No. 6 (147). P. 20-31. (In Russ.)
10. Lopashev P.A. Method for automatic account of geometric shape defects for cylinder shells in FEM analysis of bearing capacity. Proceedings of the Krylov State Research Centre. 2009. No. 46. P. 173-180. (In Russ.)
11. Lukin Ya.A., Efimenko E.S. Optimization of bearing systems of all-metal bodies of cars such as the closed shell taking into account the initial imperfections. Vestnik Rostovskogo Gosudarstvennogo Universiteta Putey Soobshcheniya. 2008. No. 1 (29). P. 9-13. (In Russ.)
12. Mamai V.I. Load-carrying capacity of thin-walled shells with local imperfections. Mechanics of Solids. 2011. V. 46, Iss. 2. P. 179-173. DOI: 10.3103/S002565441102004X
13. Korotkov A.V., Kulikov Yu.A. Forced Parametric Vibration of Curved Composite Pipes with Initial Technological Imperfections. Mekhanika kompozitsionnykh materialov i konstruktsii. 2014. V. 20, no. 1. P. 16-33. (In Russ.)
14. Miltsin A.M., Olevsky V.I., Pletin V.V. On the forms of buckling waves of non uniformly loaded cylindrical shells with technological imperfections. Eastern-European Journal of Enterprise Technologies. 2011. V. 5, no. 7 (53). P. 44-48. (In Russ.)
15. Snisarenko S.I., Hvaskov B.E. Dynamic deformation of cylindrical composite membranes with structural and technological defects. Scientific Bulletin of NSTU. 2009. No. 3. P. 115-120. (In Russ.)
16. Moiseenko M.O., Pestsov D.N., Popov O.N., Treputneva T.A. Obzor rabot po raschetu plastin i obolochek s malymi nachal'nymi nesovershenstvami formy sredinnoy poverkhnosti: deponirovannaya rukopis' [A review of studies on the calculation of plates and shells with small initial imperfections of the middle surface]. Deposited manuscript № 501-B2011; 22.11.2011. 28 p.
17. Taranukha N.A., Leizerovich G.S. Dinamika ««nepravil'nykh» obolochek [Dynamics of irregular shells]. Vladivostok: Dal'nauka Publ., 2005. 422 p.
18. Amabili M. Nonlinear vibrations and stability of shells and plates. New York: Cambridge University Press, 2008. 392 p.
19. Kubenko V.D., Koval'chuk P.S., Krasnopol'skaya T.S. Nelineynoe vzaimodeystvie form izgibnykh kolebaniy tsilindricheskikh obolochek [Nonlinear interaction of forms of flexural vibrations of cylindrical shells]. Kiev: Naukova Dumka Publ., 1984. 220 p.
20. Varadan T. K., Pratkhap Dzh., Ramani Kh.V. Nonlinear free bending vibrations of thin-walled circular cylindrical shells. Aerokosmicheskaya tekhnika. 1990. No. 5. P. 21-24. (In Russ.)
21. Seregin S.V., Leizerovich G.S. Free vibration of infinite length circular cylindrical shell with initial imperfections and a small added mass. Scholarly Notes of Komsomolsk-na-Amure State Technical University. 2014. V. 1, no. 4 (20). P. 36-43. (In Russ.)
22. Lejzerovich G.S., Prihod'ko N.B., Seregin S.V. Influence of low additional mass on splitting of the frequency spectrum of a circular ring with initial irregularities. Structural Mechanics and Analysis of Constructions. 2013. No. 6. P. 49-51. (In Russ.)
23. Lejzerovich G., Prihod'ko N., Serjogin S. On the influence of small additional masson vibration of a ring with different thicknesses. Building and Reconstruction. 2013. No. 4. P. 38-42. (In Russ.)
24. Seregin S.V. Free flexural radial vibrations of a thin circular cylindrical shell bearing added mass. VestnikMGSU. 2014. No. 11. P. 74-81. (In Russ.)
25. Seregin S.V. Numerical and analytical investigation of free vibrations of circular cylindrical shells with added mass linearly distributed along generatrix. Computational Continuum Mechanics. 2014. V. 7, no. 4. P. 378-384. (In Russ.)
26. Seregin S.V. Free vibrations of a thin circular cylindrical shell weakened by a hole. Russian Aeronautics. 2015. V. 58, Iss. 3. P. 258-262. DOI: 10.3103/S1068799815030022
27. Seregin S.V. Influence of the contact area and value of the linearly distributed and concentrated mass with a circular cylindrical shell on the frequency and modes of natural oscillations. Vestnik MGSU. 2014. No. 7. P. 64-74. (In Russ.)
28. Seregin S.V. Regarding the splitting of the frequency spectrum of the bending of thin circular cylindrical shells carrying the associated mass. Structural Mechanics and Analysis of Constructions. 2014. No. 6 (257). P. 59-61. (In Russ.)
29. Seregin S.V., Leizerovich G.S. Influence of the added mass on the dynamic characteristics of thin shells. Engineering and Automation Problems. 2015. No. 4. P. 83-89. (In Russ.)
30. Seregin S.V. Investigation of dynamic characteristics of shells with holes and added mass. Vestnik MGSU. 2014. No. 4. P. 52-58. (In Russ.)
31. Seregin S.V. The influence attached body on the frequency and forms free vibrations cylindrical shells. Structural Mechanics and Analysis of Constructions. 2014. No. 3 (254). P. 35-38. (In Russ.)
32. Vol'mir A.S. Nelineynaya dinamika plastinok i obolochek [Nonlinear dynamics of plates and shells]. Moscow: Nauka Publ., 1972. 432 p.
33. Grishin N.F. Svobodnye kolebaniya tsilindricheskikh obolochek, imeyushchikh nachal'nuyu pogib'. V kn.: Dinamika sudovykh obolochek (Materialy po obmenu opytom). Iss. 208. Leningrad: Sudostroenie Publ., 1974. P. 31-40. (In Russ.)
34. Seregin S.V. Dinamika tonkikh tsilindricheskikh obolochek s prisoedinennoy massoy [Dynamics of thin cylindrical shells with added mass]. Komsomolsk-on-Amur: Komsomolsk-na-Amure State Technical University Publ., 2016. 175 p.
