#11(728) 2020 ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
3
Машиностроение и машиноведение
УДК 534.014.4 doi: 10.18698/0536-1044-2020-11-3-12
Динамическая устойчивость трубопровода с протекающей по нему жидкостью
В.П. Радин, В.П. Чирков, А.В. Щугорев, В.Н. Щугорев, О.В. Новикова
НИУ«МЭИ»
Dynamic Stability of Pipelines with Fluid Flow
V.P. Radin, V.P. Chirkov, A.V. Shchugorev, V.N. Shchugorev, O.V. Novikova
National Research University — Moscow Power Engineering Institute
Рассмотрена классическая задача об устойчивости участка трубопровода с протекающей в нем жидкостью. Уравнение возмущенного движения получено методом разложения по формам собственных колебаний с дальнейшим применением процедуры метода Бубнова — Галеркина. Граница области устойчивости на плоскости параметров потока жидкости найдена с использованием критерия Рауса — Гурвица применительно к неконсервативным задачам устойчивости. При фиксированных значениях относительной массы построены траектории характеристических показателей как параметрически зависимые функции от скорости протекающей жидкости с определением частоты колебаний трубопровода при потере устойчивости по типу флаттера. Изучены формы флаттера в разных точках границы области устойчивости. Формы флаттера изображены пучком изогнутых осей трубопровода в дискретные моменты времени на протяжении одного периода.
Ключевые слова: устойчивость трубопровода, метод разложения по формам, граница и частота флаттера, формы колебаний трубопровода
The classical problem of stability of a pipeline section with fluid flow is considered in this paper. The equation of perturbed motion is solved by a method of expansion by forms of natural oscillations with further application of the Bubnov — Galerkin method. The boundary of the stability domain on the plane of fluid flow parameters is determined using the Raus — Hurwitz criterion for non-conservative stability problems. For fixed values of the relative mass, the trajectories of characteristic indicators are constructed as functions parametrically dependent on the velocity of the fluid flow. The frequency of pipeline oscillations in the event of loss of stability is determined by the flutter type. Flutter modes at various points of the boundary of the stability domain are examined. Flutter modes are represented by a beam of curved axes of the pipeline at discrete points of time throughout one period.
Keywords: pipeline stability, method of decomposition by form, boundary and frequency of flutter, pipeline oscillation forms
4
ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
#11(728) 2020
Среди классических неконсервативных задач теории упругой устойчивости наибольший интерес представляет задача об устойчивости участка трубопровода с протекающей по нему жидкостью [1-10]. Это связано с такими особенностями данной задачи, как двухпараметри-ческая неконсервативная нагрузка от протекающей по трубопроводу жидкости (масса и скорость), немонотонность границы флаттера и необъяснимое поведение характеристических показателей при изменении параметров неконсервативной нагрузки.
В последние годы интерес к этой задаче вырос в связи с расчетом углеродных нанотрубок [11, 12]. Отметим, что кроме традиционных для механики методов, изложенных в работах [1316], в зарубежной литературе достаточно часто для решения задач неконсервативной устойчивости применяют метод дифференциальных квадратур [8-12, 17, 18].
Суть этого метода состоит в аппроксимации производных искомой функции взвешенной суммой значений этой функции в точках дискретизации области, в которой ищут решение. Весовые коэффициенты рассчитывают с помощью интерполяционных многочленов Лан-гранжа [19, 20]. Такой метод позволяет получить решение с достаточной точностью на сетке с меньшим количеством узлов, чем у методов конечных разностей и конечных элементов.
Однако для дополнительных исследований указанной задачи применим классический метод разложения ее решения по формам собственных колебаний с дальнейшим использованием процедуры метода Бубнова — Галеркина.
Цель работы — построение траекторий характеристических показателей и наглядные изображения форм флаттера в различных точках границы области устойчивости.
Рассмотрим прямолинейный участок трубопровода, наполненный невязкой жидкостью, которая движется с невозмущенной скоростью и (рис. 1). Расчетную схему трубопровода представим в виде защемленного на одном конце стержня с изгибной жесткостью Е1, погонной массой т0 и длиной I. Массу жидкости,
УА
£7, от0 _
и
I
т
Рис. 1. Участок трубопровода с протекающей по нему жидкостью
приходящуюся на единицу длины трубопровода, обозначим через т, а перемещения сечений стержня при отклонении от прямолинейной формы — через м(х, :), где х — продольная координата; t — время.
Будем считать движение жидкости одномерным, а условия на выходе из трубопровода такими, что осевые усилия от потока отсутствуют. При этом можно трактовать жидкость как подвижную нагрузку с интенсивностью сил инерции, равной -тй2м/й:2, где й2/Л2 — полная производная по времени. С учетом того, что м = м[х(:), t], а и = йх/й:, получаем следующие выражения:
• для первой производной
ймГх(:),: 1 дм дм дх дм Тдм
-!=-=^ =-+--=-+ и-;
й: д: дх д: д: дх
• для второй производной
й2м д2м
+ 2и-
м
+ и2
м
(1)
й:2 д:2 дхд: дх2
Первый член в правой части выражения (1) учитывает обычную нормальную даламберовую силу. Второй член соответствует кориолисово-му ускорению из-за поворота потока частиц жидкости. Третий член связан с центробежными силами, определяемыми кривизной потока при колебаниях трубопровода. Тогда уравнение малых колебаний около невозмущенной формы равновесия можно записать в виде
4м , чд 2м „ ТТ д 2м
Е1--ь(т0 + т)--ъ 2ти--ъ
дх4 ^ ' д:2 дхд:
м
+ ти2-= 0.
дх2
(2)
Уравнение (2) необходимо дополнить граничными условиями
дм
м = — = 0 при х = 0; дх
д 2м д 3м
-= 0, -= 0 при х = I.
дх2 дх3
Введем безразмерные параметры:
(3)
с х 1 Е1
\ =—; Т = Ю0:; Ю0 = —-;
I I2 \т0+т
т12
т
а = Ш-; Р= .
V Е1 т0+т
где Е, — безразмерная продольная координата; х — безразмерное время; ю0 — характерная
#11(728) 2020 ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
5
частота; а — параметр, характеризующий скорость протекания жидкости и ее расход; Р — погонная масса протекающей по трубопроводу жидкости, отнесенная к их суммарной погонной массе.
С учетом безразмерных параметров уравнение (2) и граничные условия (3) запишутся следующим образом:
Э„ ггт Э. ЭЭ
- + 2а^Р
Э£Эт Эw
- + a2
w = -
Э£2 Эт2 = 0 при £ = 0;
+ —= 0; (4)
w
w
= 0, чтт = 0 при £ = 1.
Э£2 Э£3
Представим решение уравнения (4) в виде ряда
w(£, т) = X qk (т) Фк (£),
к=1
(5)
где п — число членов ряда; qk (х) — обобщенные координаты; ф к © — формы собственных колебаний консольного стержня.
Выбор числа членов ряда п определяется точностью представления решения х). Как показали многочисленные вычисления для неконсервативных задач подобного рода, достаточная точность достигается уже при п = 8 [15]. Если ввести в рассмотрение вектор обобщенных координат q(х) и вектор собственных форм ф(£), то разложение (5) может быть переписано в векторной форме
х) = q т(х)ф©;
"q1(T)" " Ф1(£)"
q(T) = q2(T) ; ф(£) = Ф2(£)
qn(т)_ Фп (£)_
Формы собственных колебаний выразим через функции Крылова
й(х£) = 1 [еЬ(х^) + ес8(х^)];
^2 (х£) = 1 [зЬ(х^) + 81и(х^)];
Яз(х£) = 2 [сЬ(х^)" сс8(х^)];
^4 (х^) = 2 [*Мх^) - З1п(х^)]. Здесь х — частотный параметр.
Для защемленного на одном конце стержня имеем
Фк (£) = S 4 (х к £) - S3(x к £),
S1(x к)
где частотные параметры хк определяются из уравнения
1 + cos xch х = 0.
Подставим ряд (5) в уравнение (4) и применим процедуру метода Бубнова — Галеркина. Для вектора обобщенных координат q(T) получим обыкновенное дифференциальное уравнение в матричном виде
Aq + 2«VpBq + [C + a 2D] q = 0, (6)
где A, B, C и D — матрицы размера nxn,
A = j ф(£) ф T(£)d£;
dy(£) "" d4 ф(£)
0
1
B = j ф(£)
0 1
C = j ф(£)
0 1
D = j ф(£)
d£4
d£2
d£;
d£;
d£.
Представляя вектор обобщенных координат q(f) как
q(í) = qо ехр(Л£),
для характеристических показателей X имеем проблему собственных значений в виде матричного полинома
Р2 X2 + РД + Ро = 0.
Здесь Р2, Р1, Р0 — матрицы-коэффициенты характеристического полинома,
Р2 = А; Р1 = 2аТРв; Ро = С + а 2Б.
В качестве варьируемых параметров системы примем относительную массу жидкости Р и безразмерную скорость ее течения а. Совокупность значений параметров Р и а, содержащихся в элементах этих матриц и соответствующих характеристическим показателям с отрицательными действительными частями, определяет область устойчивости тривиального положения равновесия.
Для трубопроводов — защемленного на одном конце и свободного на другом — потеря
6
ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
#11(728) 2020
С(а, Р) =
устойчивости происходит только колебательным образом по типу флаттера. Границу флаттера удобно построить с использованием одного из критериев отрицательности действительных частей характеристических показателей — критерия Рауса — Гурвица. Для этого составим характеристический полином Гурвица Н(а, Р) из элементов матрицы
" 0 Е "
-А -1 (С + а2 В) -2а^РА -1Б
где 0 и Е — нулевая и единичная матрицы размера п х п.
Матрица в(а, Р) получена в результате преобразования уравнения (6) к нормальной форме Коши. Равенство нулю главного минора матрицы Гурвица порядка 2п -1, т. е. А2п-1(а, Р) = 0, определяет границу области флаттера на плоскости параметров нагруже-ния а, Р [14, 15].
Отметим, что в неконсервативных системах, где возможен квазистатический тип потери устойчивости — дивергенция, соответствующая граница ее области может быть найдена из условия равенства нулю определителя матрицы Гурвица или (что то же самое) равенства нулю свободного члена характеристического полинома.
На рис. 2 сплошная линия есть граница области устойчивости прямолинейной формы равновесия трубопровода, полученная при удержании восьми членов ряда в разложении решения возмущенного движения в ряд по формам собственных колебаний. Пунктирная линия определяет частоту флаттера ю $, отне-
(X, Юу/С01
20
15
10
><
Флаттер
а* в
А, -------Т
1 1 | Устойчивость
0,25
0,50
0,75
Рис. 2. Граница области устойчивости и частота флаттера
сенную к первой собственной частоте трубопровода ю1.
Такая кривая, впервые построенная М.П. Пайдуссисом [1] в 1966 г., встречается во многих работах, посвященных устойчивости трубопроводов. Граница флаттера при Р^0 исходит из окрестности точки V при критическом значении скорости потока а* = V20,05 = = 4,48, а точнее несколько ниже этой точки. Вычисления показали, что при Р^ 0 параметр а* ^ 4,20.
Этот факт иллюстрирует известный парадокс Циглера, показывающий снижение параметров неконсервативных нагрузок при исче-зающе малом демпфировании [13-15]. Действительно, если масса протекающей по трубопроводу жидкости весьма мала, то уравнение (6) совпадает по структуре с уравнением возмущенного движения для известной задачи Бека — сжатия консольного стержня следящей силой Р. Решение этой задачи дает Р12/Е1 = 20,05.
Обсуждаемая точка V также соответствует классическому поведению характеристических показателей. Здесь имеет место явное взаимодействие первой и второй собственных форм и частот. С увеличением параметра а первая собственная частота ю1 как мнимая часть первого характеристического показателя Х1 увеличивается, а вторая ю2 уменьшается (рис. 3). Действительные части характеристических показателей Х1 и X 2 равны нулю.
Такая ситуация сохраняется вплоть до достижения параметром а критического значения а*, при котором происходит слияние частот ю1 =ю2 =ю $ и превращение их в частоту флаттера ю$, а один из характеристических показателей получает положительную действительную часть, что свидетельствует о возрастании амплитуд колебаний трубопровода.
Аналогично ведут себя характеристические показатели и в других классических задачах неконсервативной устойчивости (задача Реута, флаттер пластин и оболочек в потоке газа и др.). Безусловно, наличие в системе демпфирования вносит коррективы в зависимости характеристических показателей от параметров неконсервативной нагрузки [13-15].
Однако наиболее замысловато X к ведут себя в рассматриваемой системе при изменении параметров протекающей по трубопроводу жидкости. Проведем исследование траекторий характеристических показателей на комплексной
#11(728) 2020 ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
7
Яе X, 1т X
Рис. 3. Зависимость действительных и мнимых частей характеристических показателей от скорости потока а при Р^ 0
плоскости ИеX, 1тX при фиксированной относительной массе Р и увеличении параметра а до значения а*, при которых происходит пересечение одним (первым по очереди) из характеристических показателей мнимой оси и X переходит в правую полуплоскость, а его действительная часть становится положительной.
Этот факт соответствует динамической потери устойчивости, т. е. флаттеру. Мнимая часть X определяет частоту флаттера ю . Значения Р и а* образуют границу области устойчивости — границу флаттера. С применением метода дифференциальных квадратур такие траектории, построенные в ряде зарубежных работ, носят название диаграмм Ар-ганда [5, 8-11].
На рис. 4 построены траектории характеристических показателей X1, X 2 и фрагмент траектории X 3 при Р = 0,2. Траектории всех характеристических показателей начинаются на мнимой оси X со значений юк, соответствующих частотам собственных колебаний консольной трубы с неподвижной жидкостью как корни уравнения |С-ю2А = 0. При параметре а = 5,6 и 1тX2 =юf = 13,75 характеристический показатель X2 в точке А (см. рис. 2) переходит в правую полуплоскость и ИеX2 > 0, что соответствует флаттеру в системе.
Проведем аналогичное исследование для случая Р = 0,5 (рис. 5). Здесь неустойчивость в системе обусловлена пересечением мнимой оси в точке В характеристическим показателем X 3 при а* = 9,32 (см. рис. 2). Частота флаттера составляет 1т X 3 =ю f = 26,6.
Для пояснения немонотонности границы флаттера, обусловленной выходом и временным возвращением характеристического показателя в правую полуплоскость, проведем вычисление траекторий X к для параметра относительной массы Р = 0,88 (см. рис. 2 и 6).
Динамическая неустойчивость системы обусловлена поведением характеристического показателя X1, который при а* = 14,2 в точке С переходит в правую полуплоскость, но с ростом а до значения а* = 15,7 он возвращается в левую полуплоскость с отрицательной действительной частью в точке Л плоскости ИеX, 1тX. Дальнейшее увеличение параметра а приводит к тому, что при а* = 16,9 в точке Е характеристический показатель X1 опять переходит в правую полуплоскость уже безвозвратно.
1т X
20
10
0
-15
Ху/
( (^2
1 А
Юу
-10
-5
0 ЯеХ
Рис. 4. Траектории характеристических показателей при Р = 0,2
1т X
20
10
0
-15
4^3 в
Х-у
---"Ю!
-10
-5
0 ЯеХ
Рис. 5. Траектории характеристических показателей при Р = 0,5 и а = уаг
8
ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
#11(728) 2020
1т X
Е/
СО3
Ч / / V су ( м
со2
\ /С с? Ух ®1
-50 -25 0 Ле А,
Рис. 6. Траектории характеристических показателей при Р = 0,88
Таким образом, граница области устойчивости в зависимости от параметра, характеризующего относительную массу протекающей жидкости, определяется выходом различных характеристических показателей в правую полуплоскость. Так, если параметр Ре [0; 0,384], то мнимую ось слева направо пересекает характеристический показатель X 2, определяя при этом значения параметра а*, т. е. границу области устойчивости, а также частоту флаттера ю f = 1т X 2.
Для параметра Ре [0,384; 0,530] границу и частоту флаттера определяет показатель X 3. В диапазоне параметра Ре [0,530; 0,614] вновь инициатива по смене устойчивости на неустойчивость принадлежит показателю X 2. И, наконец, при Ре [0,614; 1,000] кривую а*(Р) определяет характеристический показатель Х1.
Аналогичные результаты с применением метода дифференциальных квадратур для определения характерных участков на плоскости параметров Р, а получены в работе [17] и сделаны выводы о преобладающих модах в колебательном движении трубопровода при потере устойчивости. Влияние скорости потока жидкости на изменение модовых форм рассмотрено в статье [10].
Проведем исследование формы флаттера на основе метода разложения решения уравнения возмущенного движения по формам собственных колебаний трубопровода.
Зависимость частоты флаттера, отнесенная к первой собственной частоте пустого трубопровода, показана на рис. 2 пунктирной линией. Из рисунка следует, что характеры кривой ю f (Р) и
границы области устойчивости а*(Р) вполне согласованы. Поведение характеристических показателей, а точнее их расположение в момент динамической потери устойчивости трубопровода, определяют форму флаттера.
Наглядную картину этих форм флаттера дают фигуры, построенные на рис. 7, а-в. Здесь для трех характерных точек на границе флаттера и для р = 30 моментов времени, на которые поделен период колебаний трубопровода Т = 2 л/ю!, построены формы изгиба оси трубопровода для моментов времени тк = кДТ, к = 0,1, 2, ...,30, где ДТ = Т/р.
Причем выбран фрагмент движения, соответствующий установившимся колебаниям трубопровода, возникающим после некоторых весьма малых начальных условий для обобщенных координат qk (т). Деформации трубопровода при безразмерной продольной координате ^ = 0...1 в каждый момент времени определяются выражением (6).
На рис. 7, а на границе области устойчивости для а* = 5,6 и Р = 0,2, где при малых начальных возмущениях реализуется флаттер с частотой ю! = 13,72, изображены формы изгиба трубопровода в течение одного периода. Здесь основной вклад в форму флаттера вносят первая и вторая собственные формы, что наблюдается визуально по рисунку, а также подтверждено соотношениями между обобщенными координатами.
Потеря устойчивости связана с выходом в правую полуплоскость характеристического показателя X 2 (см. рис. 4). Увеличение относительной массы жидкости Р повышает критическое значение а*, частоту флаттера ю f, а также инициирует смену характеристического показателя с X2 на X3 (см. рис. 5), ответственного за потерю устойчивости, действительная часть которого становится положительной.
На рис. 7, б для точки на границе флаттера с параметрами а* = 9,32 и Р = 0,5 с траекториями X к, изображенными на рис. 5, построена форма изгиба трубопровода при динамической потере устойчивости. Здесь уже явно проявляется вклад третьей собственной формы в форму флаттера с частотой ю f = 26,55. Наиболее любопытными и неожиданными являются случаи значений Р, близких к единице. Эта область характеризуется сложными траекториями характеристических показателей.
На рис. 7, в показана форма флаттера при а* = 13,53 и Р = 0,8. Здесь наступление флатте-
#11(728)2020
ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
9
Рис. 7. Формы флаттера на границе области устойчивости при различных значениях
безразмерных параметров: а — а = 5,60, Р = 0,2; б — а = 9,32, Р = 0,5; в — а = 13,53, Р = 0,8
ра связано с увеличением 1т X и выходом X в правую полуплоскость. Причем этот выход происходит при значениях, близких к ю3, при существенном снижении X 4. Видимо эти причины и приводят к тому, что формы изгиба трубопровода (форма флаттера) содержат формы собственных колебаний с первой по четвертую.
Выводы
1. Полученные данные существенно дополняют результаты в области расчета устойчивости упругих систем, находящихся в условиях некон-
сервативного нагружения. В частности, на основе метода разложения решения уравнения возмущенного движения по формам собственных колебаний построены траектории характеристических показателей в зависимости от параметров нагружения, определены частоты динамической потери устойчивости. Получены наглядные рисунки, характеризующие формы флаттера на границе области устойчивости.
2. Результаты работы могут быть использованы при проектировании и расчете на устойчивость технических объектов, содержащих трубопроводы с протекающей по ним жидкостью.
10
ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
#11(728) 2020
Литература
[1] Gregory R.W., Paidoussis M.P. Unstable Oscillation of Tubular Cantilever Fluid. Proceedings
of the Royal Society of London. Series A, Mathematical and Physical Series, 1966, vol. 293, iss. 1435, pp. 528-542.
[2] Феодосьев В.И. Избранные задачи и вопросы по сопротивлению материалов. Москва,
Наука, 1973. 400 с.
[3] Elishakoff I., Vittori P. A paradox of non-monotonicity in stability of pipes conveying fluid.
Theoretical and Applied Mechanics, 2005, vol. 32, pp. 235-282.
[4] Каган-Розенцвейг Л.М. О механизме потери устойчивости равновесия консольной
трубы с протекающей жидкостью. Вестник гражданских инженеров, 2012, № 1(30), с. 102-107.
[5] Marzani A., Mazzotti M., Viola E., Vittori P., Elishakoff I. FEM Formulation for Dynamic In-
stability of Fluid-Conveying Pipe on No uniform Elastic Foundation. Mechanics Based Design of Structures and Machines: An International Journal, 2012, vol. 40(1), pp. 83-95, doi: 10.1080/15397734.2011.618443
[6] Bahaadini, R., Hosseini, M. Flow-induced and mechanical stability of cantilever carbon
nanotubes subjected to an axial compressive load. Applied Mathematical Modelling, 2018, vol. 59, pp. 597-613, doi: https://doi.org/10.1016/j.apm.2018.02.015
[7] Paidoussis M.P. The canonical problem of the fluid-conveying pipe and radiation of the
knowledge gained to other dynamics problems across. Journal of Sound and Vibration, 2008, vol. 310, pp. 462-492, doi: https://doi.org/10.1016/j.jsv.2007.03.065
[8] Bellman R.E., Kashef B.G., Casti J. Differential quadrature: A technique for the rapid solution
of nonlinear partial differential equations. Journal of Computational Physics, 1972, vol. 10, iss. 1, pp. 40-52, doi: 10.1016/0021-9991(72)90089-7
[9] Shu C. Differential Quadrature and Its Application in Engineering. Springer-Verlag, 2000. 340 p.
[10] Wang L., Dai H.L., Ni Q. Nonconservative pipes conveying fluid: evolution of mode shapes with increasing flow velocity. Journal of Vibration and Control, 2015, vol. 21(16), pp. 33593367, doi: 10.1177/1077546314522490
[11] Bahaadini R., Mohammad H. Flow-induced and mechanical stability of cantilever carbon nanotubes subjected to an axial compressive load. Applied Mathematical Modelling, 2018, vol. 59, pp. 597-613, doi: 10.1016/j.apm.2018.02.015
[12] Bahaadini R., Mohammad R. D., Mohammad H., Zahra K. Stability analysis of composite thin-walled pipes conveying fluid. Ocean Engineering, 2018, vol. 160, pp. 311-323, doi: 10.1016/j.oceaneng.2018.04.061
[13] Болотин В.В. Неконсервативные задачи теории упругой устойчивости. Москва, Физ-матгиз, 1961. 339 с.
[14] Радин В.П., Чирков В.П., Щугорев А.В., Щугорев В.Н. Методы определения критических значений неконсервативных нагрузок в задачах устойчивости механических систем. Известия высших учебных заведений. Машиностроение, 2019, № 10, с. 3-13, doi: 10.18698/0536-1044-2019-10-3-13
[15] Радин В.П., Самогин Ю.Н., Чирков В.П., Щугорев А.В. Решение неконсервативных задач теории устойчивости. Москва, Физматлит, 2017. 240 с.
[16] Seyranian A.R., Elishakoff I. Modern problem of structural stability. New York, SpringerVerlag Wien, 2002. 394 p.
[17] Tornabene F., Marzani A., Viola E., Elishakoff I. Critical Flow Speeds of Pipes Conveying Fluid Using the Generalized Differential Quadrature Method. Advances in Theoretical and Applied Mechanics, 2010, vol. 3, no. 3, pp. 121-138.
[18] Marzani A., Tornabene F., Viola E. Nonconservative stability problems via Generalized Differential Quadrature method. Journal of Sound and Vibration, 2008, vol. 315, pp. 176-196, doi:10.1016/j.jsv.2008.01.056
[19] Барулина М.А. Применение обобщенного метода дифференциальных квадратур к решению двумерных задач механики. Известия Саратовского университета. Новая сер. Сер. Математика. Механика. Информатика, 2018, т. 18, вып. 2, с. 206-216, doi: 10.18500/1816-9791-2018-18-2-206-216
[20] Вержбицкий В.М. Численные методы (математический анализ и обыкновенные дифференциальные уравнения). Москва, Директ-Медиа, 2013. 400 с.
#11(728) 2020 ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
11
References
[1] Gregory R.W., Paidoussis M.P. Unstable Oscillation of Tubular Cantilever Fluid. Proceedings
of the Royal Society of London. Series A, Mathematical and Physical Series, 1966, vol. 293, iss. 1435, pp. 528-542.
[2] Feodos'yev V.I. Izbrannyye zadachi i voprosy po soprotivleniyu materialov [Selected tasks and
questions on the resistance of materials]. Moscow, Nauka publ., 1973. 400 p.
[3] Elishakoff I., Vittori P. A paradox of non-monotonicity in stability of pipes conveying fluid.
Theoretical and Applied Mechanics, 2005, vol. 32, pp. 235-282.
[4] Kagan-Rozentsveyg L.M. On the stability loss mechanism of tubular cantilever conveying flu-
id. Vestnik grazhdanskikh inzhenerov, 2012, no. 1(30), pp. 102-107 (in Russ.).
[5] Marzani A., Mazzotti M., Viola E., Vittori P., Elishakoff I. FEM Formulation for Dynamic In-
stability of Fluid-Conveying Pipe on No uniform Elastic Foundation. Mechanics Based Design of Structures and Machines: An International Journal, 2012, vol. 40(1), pp. 83-95, doi: 10.1080/15397734.2011.618443
[6] Bahaadini, R., Hosseini, M. Flow-induced and mechanical stability of cantilever carbon
nanotubes subjected to an axial compressive load. Applied Mathematical Modelling, 2018, vol. 59, pp. 597-613, doi: https://doi.org/10.1016/jj.apm.2018.02.015
[7] Paidoussis M.P. The canonical problem of the fluid-conveying pipe and radiation of the
knowledge gained to other dynamics problems across. Journal of Sound and Vibration, 2008, vol. 310, pp. 462-492, doi: https://doi.org/10.1016/jjsv.2007.03.065
[8] Bellman R.E., Kashef B.G., Casti J. Differential quadrature: A technique for the rapid solution
of nonlinear partial differential equations. Journal of Computational Physics, 1972, vol. 10, iss. 1, pp. 40-52, doi: 10.1016/0021-9991(72)90089-7
[9] Shu C. Differential Quadrature and Its Application in Engineering. Springer-Verlag, 2000.
340 p.
[10] Wang L., Dai H.L., Ni Q. Nonconservative pipes conveying fluid: evolution of mode shapes with increasing flow velocity. Journal of Vibration and Control, 2015, vol. 21(16), pp. 33593367, doi: 10.1177/1077546314522490
[11] Bahaadini R., Mohammad H. Flow-induced and mechanical stability of cantilever carbon nanotubes subjected to an axial compressive load. Applied Mathematical Modelling, 2018, vol. 59, pp. 597-613, doi: 10.1016/j.apm.2018.02.015
[12] Bahaadini R., Mohammad R. D., Mohammad H., Zahra K. Stability analysis of composite thin-walled pipes conveying fluid. Ocean Engineering, 2018, vol. 160, pp. 311-323, doi: 10.1016/j.oceaneng.2018.04.061
[13] Bolotin V.V. Nekonservativnyye zadachi teorii uprugoy ustoychivosti [Non-conservative problems of the theory of elastic stability]. Moscow, Fizmatgiz publ., 1961. 339 p.
[14] Radin V.P., Chirkov V.P., Shchugorev A.V., Shchugorev V.N. Methods for Determinning Critical Values of Nonconservative Loads in Problems of Stability of Mechanical Systems. Proceedings of Higher Educational Institutions. Machine Building, 2019, no. 10, pp. 3-13 (in Russ.), doi: 10.18698/0536-1044-2019-10-3-13
[15] Radin V.P., Samogin Yu.N., Chirkov V.P., Shchugorev A.V. Resheniye nekonservativnykh zadach teorii ustoychivosti [Solution of non-conservative problems of stability theory]. Moscow, Fizmatlit publ., 2017. 240 p.
[16] Seyranian A.R., Elishakoff I. Modern problem of structural stability. New York, SpringerVerlag Wien, 2002. 394 p.
[17] Tornabene F., Marzani A., Viola E., Elishakoff I. Critical Flow Speeds of Pipes Conveying Fluid Using the Generalized Differential Quadrature Method. Advances in Theoretical and Applied Mechanics, 2010, vol. 3, no. 3, pp. 121-138.
[18] Marzani A., Tornabene F., Viola E. Nonconservative stability problems via Generalized Differential Quadrature method. Journal of Sound and Vibration, 2008, vol. 315, pp. 176-196, doi:10.1016/j.jsv.2008.01.056
[19] Barulina M.A. Application of Generalized Differential Quadrature Method to Two-dimensional Problems of Mechanics. Izvestiya Saratovskogo universiteta. Novaya ser. Ser. Matematika. Mekhanika. Informatika, 2018, vol. 18, iss. 2, pp. 206-216, doi: 10.18500/1816-9791-2018-18-2-206-216
12
ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
#11(728) 2020
[20] Verzhbitskiy V.M. Chislennyye metody (matematicheskiy analiz i obyknovennyye differentsi-al'nyye uravneniya) [Numerical methods (mathematical analysis and ordinary differential equations)]. Moscow, Direkt-Media publ., 2013. 400 p.
Информация об авторах
РАДИН Владимир Павлович — кандидат технических наук, профессор кафедры «Робототехника, мехатроника, динамика и прочность машин». НИУ «МЭИ» (111250, Москва, Российская Федерация, Красноказарменная ул., д. 14, e-mail: radinvp@gmail.com).
ЧИРКОВ Виктор Петрович — доктор технических наук, профессор кафедры «Робототехника, мехатроника, динамика и прочность машин». НИУ «МЭИ» (111250, Москва, Российская Федерация, Красноказарменная ул., д. 14, e-mail: victor.p.chirkov@gmail.com).
ЩУГОРЕВ Алексей Владимирович — кандидат технических наук, ассистент кафедры «Робототехника, мехатро-ника, динамика и прочность машин». НИУ МЭИ (111250, Москва, Российская Федерация, Красноказарменная ул., д. 14, e-mail: alexey.salut@gmail.com).
ЩУГОРЕВ Владимир Николаевич — кандидат технических наук, доцент кафедры «Робототехника, мехатроника, динамика и прочность машин». НИУ «МЭИ» (111250, Москва, Российская Федерация, Красноказарменная ул., д. 14, e-mail: shchugorev.vn@gmail.com).
НОВИКОВА Ольга Валерьевна — кандидат технических наук, доцент кафедры «Робототехника, мехатроника, динамика и прочность машин». НИУ «МЭИ» (111250, Москва, Российская Федерация, Красноказарменная ул., д. 14, e-mail: NovicovaOV@mpei.ru).
Статья поступила в редакцию 28.06.2020 Information about the authors
RADIN Vladimir Pavlovich — Candidate of Science (Eng.), Professor, Department of Robotics, Mechatronics, Dynamics and Strength of Machines. National Research University — Moscow Power Engineering Institute (111250, Moscow, Russian Federation, Krasnokazarmennaya St., Bldg. 14, e-mail: radinvp@gmail.com).
CHIRKOV Victor Petrovich — Doctor of Science (Eng.), Professor, Department of Robotics, Mechatronics, Dynamics and Strength of Machines. National Research University — Moscow Power Engineering Institute (111250, Moscow, Russian Federation, Krasnokazarmennaya St., Bldg. 14, e-mail: victor.p.chirkov@gmail.com).
SHCHUGOREV Aleksei Vladimirovich — Candidate of Sciences (Eng.), Teaching Assistant, Department of Robotics, Mechatronics, Dynamics and Strength of Machines. National Research University — Moscow Power Engineering Institute (111250, Moscow, Russian Federation, Krasnokazarmen-naya St., Bldg. 14, e-mail: alexey.salut@gmail.com).
SHCHUGOREV Vladimir Nikolaevich — Candidate of Sciences (Eng.), Associate Professor, Department of Robotics, Mechatronics, Dynamics and Strength of Machines. National Research University — Moscow Power Engineering Institute (111250, Moscow, Russian Federation, Krasnokazarmennaya St., Bldg. 14, e-mail: shchugorev.vn@gmail.com).
NOVIKOVA Olga Valerievna — Candidate of Sciences (Eng.), Associate Professor, Department of Robotics, Mecha-tronics, Dynamics and Strength of Machines. National Research University — Moscow Power Engineering Institute (111250, Moscow, Russian Federation, Krasnokazarmennaya St., Bldg. 14, e-mail: NovicovaOV@mpei.ru).
Просьба ссылаться на эту статью следующим образом:
Радин В.П., Чирков В.П., Щугорев А.В., Щугорев В.Н., Новикова О.В. Динамическая устойчивость трубопровода с протекающей по нему жидкостью. Известия высших учебных заведений. Машиностроение, 2020, № 11, с. 3-12, doi: 10.18698/0536-1044-2020-11-3-12
Please cite this article in English as: Radin V.P., Chirkov V.P., Shchugorev A.V., Shchugorev V.N., Novikova O.V. Dynamic Stability of Pipelines with Fluid Flow. BMSTU Journal of Mechanical Engineering, 2020, no. 11, pp. 3-12, doi: 10.18698/0536-1044-2020-11-3-12