Научная статья на тему 'ДИНАМИКА ВАЛОПРОВОДА ГРЕБНОГО ВИНТА ПРИ РАЗГОННЫХ РЕЖИМАХ'

ДИНАМИКА ВАЛОПРОВОДА ГРЕБНОГО ВИНТА ПРИ РАЗГОННЫХ РЕЖИМАХ Текст научной статьи по специальности «Механика и машиностроение»

CC BY
133
18
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ВАЛОПРОВОД / РЕЖИМ РАЗГОНА / МАТЕМАТИЧЕСКАЯ МОДЕЛЬ / ДИНАМИЧЕСКИЕ УСИЛИЯ / КРУТИЛЬНЫЕ КОЛЕБАНИЯ / СОБСТВЕННАЯ ЧАСТОТА / ДИССИПАТИВНЫЕ СИЛЫ / МЕТОД ФУРЬЕ

Аннотация научной статьи по механике и машиностроению, автор научной работы — Царенко Сергей Николаевич, Рак Александр Николаевич, Безлобенко Борис Николаевич

В работе приведены результаты исследования динамических усилий, возникающих в валопроводе при вводе главного двигателя в номинальный режим работы. Расчетная схема валопровода представляет собой упругий однородный стержень с распределенным моментом инерции. Левый торец стержня поворачивается по заданному закону движения, на правом торце расположен винт в виде инертного диска, на который действует момент сопротивления, пропорциональный квадрату угловой скорости. Математическая модель напряженно-деформированного состояния вала представлена нелинейным дифференциальным уравнением второго порядка в частных производных. Для практической реализации математической моделипредложен ряд допущений, которые позволяют ее упростить. Решение упрощенной математической модели получено методом Фурье для собственных функций, ортогональных с весом. Закон движения левого торца вала, имитирующий разгон главного двигателя, принят в виде экспоненциальной функции. В качестве примера для расчета использованы энергомеханические характеристики большого морозильного рыболовного траулера «Механик Ковтун». Результаты исследования показали, что динамические усилия могут многократно превышать номинальное значение, если разгон главного двигателя осуществляется в течение нескольких периодов главной формы собственных колебаний валопровода. При плавном разгоне главного двигателя с выходом на номинальную частоту вращения за время, превышающее двадцать периодов главной формы собственных колебаний, динамический эффект становится несущественным. Предложенная математическая модель также может быть использована для исследования напряженно-деформированного состояния валопровода при различных динамических воздействиях, таких как выход винта из воды, попадание в неоднородный набегающий поток, удар винта о лед, резонансные колебания при пульсирующей работе двигателя и др.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по механике и машиностроению , автор научной работы — Царенко Сергей Николаевич, Рак Александр Николаевич, Безлобенко Борис Николаевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

PROPELLER SHAFT LINE DYNAMICS AT ACCELERATION MODE

The results of the study of the dynamic forces that arise in the shafting when the main engine is brought into the nominal operating mode are presented in the work. The design scheme of the shafting is an elastic homogeneous rod with a distributed moment of inertia. The left end of the shaft rotates according to a given law of motion, on the right end there is a screw in the form of an inert disk, on which a moment of resistance acts proportional to the square of the angular velocity. The mathematical model of the stress-strain state of the shaft is represented by a nonlinear differential equation of the second order in partial derivatives. For the practical calculation of the mathematical model, assumptions that make it possible to simplify the mathematical model are introduced. The solution of a simplified mathematical model is obtained by the Fourier method for eigenfunctions orthogonal to weight. The law of motion of the left end of the shaft, simulating the acceleration of the main engine, is taken as an exponential function. For an example of the calculation, the energy-mechanical characteristics of a large freezing fishing trawler MECHANIC KOVTUN are taken. The results of the study have showed that the dynamic forces can be many times higher than the nominal value if the main engine is accelerated during several periods of the main form of natural vibrations of the shafting. The dynamic effect is small at a slow acceleration of the main engine, reaching the rated speed for more than twenty periods of the main mode of oscillation. The proposed mathematical model can also be used to study the stress-strain state of the shafting under various dynamic loads, such as: the propeller leaves the water, it gets into an inhomogeneous incoming flow of water, the propeller hits ice, resonant vibrations during pulsating engine operation, etc.

Текст научной работы на тему «ДИНАМИКА ВАЛОПРОВОДА ГРЕБНОГО ВИНТА ПРИ РАЗГОННЫХ РЕЖИМАХ»

*кВЕСТНИК

ш-Г-............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

VjyiOPCKOro И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

СУДОСТРОЕНИЕ И СУДОРЕМОНТ

DOI: 10.21821/2309-5180-2021-13-4-548-558

PROPELLER SHAFT LINE DYNAMICS AT ACCELERATION MODE

S. N. Tzarenko1, А. N. Rak2, B. N. Bezlobenko3

1 — Kamchatka State Technical University, Petropavlovsk-Kamchatskу, Russian Federation

2 — Donetsk National Technical University, Donetsk, Ukraine

3 — Petropavlovsk-Kamchatsky Shipyard, Petropavlovsk-Kamchatskу, Russian Federation

The results of the study of the dynamic forces that arise in the shafting when the main engine is brought into the nominal operating mode are presented in the work. The design scheme of the shafting is an elastic homogeneous rod with a distributed moment of inertia. The left end of the shaft rotates according to a given law of motion, on the right end there is a screw in the form of an inert disk, on which a moment of resistance acts proportional to the square of the angular velocity. The mathematical model of the stress-strain state of the shaft is represented by a nonlinear differential equation of the second order in partial derivatives. For the practical calculation of the mathematical model, assumptions that make it possible to simplify the mathematical model are introduced. The solution of a simplified mathematical model is obtained by the Fourier methodfor eigenfunctions orthogonal to weight. The law of motion of the left end of the shaft, simulating the acceleration of the main engine, is taken as an exponential function. For an example of the calculation, the energy-mechanical characteristics of a large freezing fishing trawler MECHANIC KOVTUN are taken. The results of the study have showed that the dynamic forces can be many times higher than the nominal value if the main engine is accelerated during several periods of the main form of natural vibrations of the shafting. The dynamic effect is small at a slow acceleration of the main engine, reaching the rated speed for more than twenty periods of the main mode of oscillation. The proposed mathematical model can also be used to study the stress-strain state of the shafting under various dynamic loads, such as: the propeller leaves the water, it gets into an inhomogeneous incoming flow of water, the propeller hits ice, resonant vibrations during pulsating engine operation, etc.

Keywords: shaft line, acceleration mode, mathematical model, dynamic forces, torsional vibrations, natural frequency, dissipative forces, Fourier method.

For citation:

Tzarenko, Sergey N., Alexander N. Rak, and Boris N. Bezlobenko. "Propeller shaft line dynamics at acceleration

mode." Vestnik Gosudarstvennogo universiteta morskogo i rechnogo flota imeni admirala S. O. Makarova

13.4 (2021): 548-558. DOI: 10.21821/2309-5180-2021-13-4-548-558.

УДК 539.3/.6:629.1

ДИНАМИКА ВАЛОПРОВОДА ГРЕБНОГО ВИНТА ПРИ РАЗГОННЫХ РЕЖИМАХ

^ С. Н. Царенко1, А. Н. Рак2, Б. Н. Безлобенко3

г

1 — Камчатский государственный технический университет, г. Петропавловск-Камчатский, Российская Федерация

2 — Донецкий национальный технический университет, г. Донецк, Украина

3 — Петропавловск-Камчатский судоремонтный завод, г. Петропавловск-Камчатский, Российская Федерация

В работе приведены результаты исследования динамических усилий, возникающих в валопроводе при вводе главного двигателя в номинальный режим работы. Расчетная схема валопровода представляет собой упругий однородный стержень с распределенным моментом инерции. Левый торец стержня поворачивается по заданному закону движения, на правом торце расположен винт в виде инертного диска, на который действует момент сопротивления, пропорциональный квадрату угловой скорости. Математическая модель напряженно-деформированного состояния вала представлена нелинейным дифференциальным уравнением второго порядка в частных производных. Для практической реализации математической модели

ВЕСТНИК«!

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ^^

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

предложен ряд допущений, которые позволяют ее упростить. Решение упрощенной математической модели получено методом Фурье для собственных функций, ортогональных с весом. Закон движения левого торца вала, имитирующий разгон главного двигателя, принят в виде экспоненциальной функции. В качестве примера для расчета использованы энергомеханические характеристики большого морозильного рыболовного траулера «Механик Ковтун». Результаты исследования показали, что динамические усилия могут многократно превышать номинальное значение, если разгон главного двигателя осуществляется в течение нескольких периодов главной формы собственных колебаний валопровода. При плавном разгоне главного двигателя с выходом на номинальную частоту вращения за время, превышающее двадцать периодов главной формы собственных колебаний, динамический эффект становится несущественным. Предложенная математическая модель также может быть использована для исследования напряженно-деформированного состояния валопровода при различных динамических воздействиях, таких как выход винта из воды, попадание в неоднородный набегающий поток, удар винта о лед, резонансные колебания при пульсирующей работе двигателя и др.

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

Для цитирования:

Царенко С. Н. Динамика валопровода гребного винта при разгонных режимах / С. Н. Царенко, А. Н. Рак, Б. Н. Безлобенко // Вестник Государственного университета морского и речного флота имени адмирала С. О. Макарова. — 2021. — Т. 13. — № 4. — С. 548-558. DOI: 10.21821/2309-5180-2021-13-4548-558.

Введение (Introduction)

Валопровод является одним из наиболее нагруженных узлов, входящих в состав судовых энергетических установок. Основной функцией валопровода является передача крутящего момента от главного двигателя (ГД) движителю. Помимо кручения он подвергается также продольным и поперечным нагрузкам, которые имеют как статическую, так и динамическую составляющие [1]. Расчет на статические нагрузки не вызывает трудностей и выполняется классическими методами сопротивления материалов. Учету динамических нагрузок, в силу многообразия условий их возникновения, посвящены исследования как признанных, так и современных ученых, краткий исторический обзор развития теории исследования крутильных колебаний валопроводов представлен в работе [2].

Как правило, систему подвижных узлов при исследовании ее динамики рассматривают изолированно, соединенной жестко [1], [3]-[6] или упруго [7], [8] с корпусом судна, при этом может дополнительно учитываться гироскопический эффект от движения судна в целом [4], [7]. Традиционно сформировалось несколько подходов к решению задачи крутильных колебаний, но общие требования к выполнению расчетов регламентированы Регистром судоходства1. Наиболее распространенным является метод дискретизации, суть которого состоит в том, что все конструктивные элементы валопровода, включая подвижные узлы ГД, соединительные муфты, гребной винт и другие механизмы, заменяются системой инертных дисков, соединенных между собой упругими связями [1]-[3], [5], [6]. Данная модель отработана для типовых схем и хорошо согласуется с результатами экспериментальных исследований [9]. Помимо упругих связей в нее могут быть включены демпфирующие элементы [3], [6]. Однако этот подход имеет свои недостатки, а именно точность определения динамических усилий зависит от уровня дискретизации модели, а кроме того, учитывая, что для ее расчета в основном используют численные методы, возникает проблема с верификацией результатов, особенно при появлении нелинейных членов в уравнениях.

Альтернативой метода дискретизации является рассмотрение валопровода как системы с распределенными и сосредоточенными массами. В этом случае динамическая задача сводится к решению дифференциального уравнения второго - четвертого порядка в частных производных, для решения которого используют точные и приближенные методы расчета [10]. Данный подход не получил широкого распространения при исследовании динамики крутильных колебаний, но успешно используется для решения задач изгибных колебаний валопровода [8], [11]. В работе [12]

2 о 2

1 Правила классификации и постройки морских судов. Ч. VII. Механические установки. НД № 2-020101-124. СПб.: Российский морской регистр судоходства, 2020. 106 с.

^ВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^уЮРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

модальный анализ роторной системы выполнен с использованием программного комплекса объектного моделирования на основе метода конечных элементов. Преимуществом такого подхода является возможность получения полного спектра частот собственных крутильных, продольных и изгибных колебаний на базе одной модели. Однако сложно судить о качестве такого расчета без четкого обоснования выбора типа элемента, плотности сетки при разбивке на элементы, назначения граничных условий и задания шага времени для динамических задач.

В существующих динамических расчетах собственные колебания, как правило, не учитывают, так как из-за множественного действия диссипативных сил они являются быстро затухающими [1], при этом исключение составляет исследование резонансных явлений [3], [6]. Тем не менее в теории строительных конструкций известно, что кратковременно динамический эффект, в котором значительную роль играют собственные колебания, в случаях удара или внезапного приложения нагрузки может приводить к увеличению внутренних усилий в два раза и более [13], [14]. В работе [15] при статистическом анализе повреждений валопроводов указывается факт значительных нагрузок в режиме запуска дизельного двигателя, а также при зигзагообразных маневрах. В трудах [16], [17] представлены результаты теоретических и экспериментальных исследований динамики гребного винта и валопровода для различных типов силового воздействия, включая импульсное. Тем не менее в специальной или отраслевой литературе в полной мере данный вопрос не раскрыт, а множественные повреждения и интенсивный износ узлов валопровода лишь подтверждают актуальность дальнейших исследований в этом направлении.

Целью работы является определение условий возникновения динамических усилий в вало-проводе в режиме разгона ГД.

г

Методы и материалы (Methods and Materials)

Систему, состоящую из гребного и промежуточных валов, будем рассматривать как однородный упругий стержень длиной l с крутильной жесткостью GJp и равномерно распределенным по длине моментом инерции масс yJ (рис. 1). С правой стороны вала располагается гребной винт с моментом инерции 1М, а также действует момент гидродинамического сопротивления, пропорциональный квадрату угловой скорости:

M = R.

, /)

где Rc — коэффициент, учитывающий физические свойства среды и геометрию винта, который ориентировочно определяется через номинальную угловую скорость и мощность на винте [18]; у (/в, — угол поворота винта; ^ — время.

Рис. 1. Расчетная схема валопровода

2

Левый торец стержня поворачивается по заданному закону движения: 0 = f ) (см. рис. 1), моделирующего разгон ГД. Дифференциальное уравнение движения сечений вала можно представить в виде [19]:

д2у д2у д3у

дх

д1

' д1дх2

и дУ В

'и 2 а7 " В

ду

а7

5(х-/) = т(х, ,

(1)

где у — угол поворота произвольного сечения вала; и ц 2 — коэффициенты, характеризующие внутреннее и внешнее затухание соответственно (для рассматриваемой задачи ограничимся случаем: |и12 = 0; т(х, 0 = 0); 5(х) — дельта-функция Дирака; т (х, — интенсивность внешнего скручивающего момента; х — продольная координата.

Движение произвольного сечения вала будем рассматривать как сложное, состоящее из переносного движения левого торца вала и относительного движения деформации вала. Тогда углы поворота для произвольного сечения будут определяться зависимостью

у = 0 + ф,

(2)

где ф — угол закручивания.

Подставляя выражение (2) в уравнение (1), а также учитывая, что угол 0 не зависит от продольной координаты, получим уравнение относительно углов закручивания:

п Г2Ф л

2 Л Г0Гф + *

Г2 0

Х 5(Х " 1 ) = РР й0 + Л

П г

2

Гф

а7

2

(3)

и* -').

Зависимость (3) является нелинейным неоднородным уравнением с переменным коэффициентом, аналитическое решение которого для общего случая неизвестно. С целью упрощения этого уравнения введем следующие допущения:

1) будем полагать, что ф << 0, это позволит исключить из уравнения нелинейный член;

2) учитывая, что первое слагаемое в скобках по своей сути соответствует локально-действующей диссипативной силе, и при выходе на номинальную скорость 0ном будет происходить гашение колебаний (ф ^ 0), вместо мгновенной будем использовать среднюю скорость в течение рассматри-

• 0(2 )-0(^ ) ваемого промежутка времени & = ¿2 - ^ , 0ср = ———;

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

В уравнении (3) для удобства перейдем к безразмерным величинам и с учетом принятых допущений получим

аг2 ах2 рах ^ ; с

а0

ах

\ 2

Ыь-1) +

а20

ах

2

(4)

Здесь = х /1 — безразмерная координата; т = ^ /1 — безразмерное время, с = .^оТ^ТуТР — скорость распространения крутильных волн в стержне [14]; ц = с ср ; Яс = С .

С-1?1

В справочнике [1] приведена ступенчатая диаграмма ввода дизельного двигателя в режим полного хода, в соответствии с которой наиболее подходящей для описания режима разгона при выводе на номинальную угловую скорость 0ном на каждой ступени будет экспоненциальная функция вида

е(т) = бном (1 - е-).

(5)

Здесь а — параметр, который зависит от времени выхода на номинальную скорость. Тогда закон движения левого торца стержня можно представить в виде

еМ = вном -

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

с

т — -

1 — е"

а

2 О 2

■Е

(6)

ES

со

В дальнейшем точкой будем обозначать производную по т, штрихом — по Для определения динамических усилий необходимо решить граничную задачу (4):

ф(0, т) = 0;

JФ'(1, т) = ^ ф(1, т), (7)

с однородными начальными условиями:

ф(0 ) = о; < ( 0 ) = 0. (8)

Согласно методу Фурье, решение уравнения (4) найдем в виде ряда произведения функций:

ф( . (*). (9)

n=1

Собственные функции Ф; возьмем такими, которые удовлетворяют уравнению вида [19]:

Ф;Ф , = о, (10)

где X — собственные значения.

^ n

С учетом первого граничного условия (7) решение уравнения (10) можно представить в виде

= ^ . (11)

к n

Из второго граничного условия (7) получим уравнение для нахождения собственных значений:

tan X; =-1-. (12)

Здесь £ = IM / [yJp/) — относительная инерционная нагрузка.

Собственные функции (11) будут ортогональны с весом р(^) = 1 + £8(1 — [14]. Подставив ряд (9) в выражение (4), умножив его обе части на рФn и взяв интеграл по длине, получим уравнение для нахождения функций

n+Í?T+-У ( (1)-Ф: (0)-№. (1))Й , (13)

w _____ . , . . ^ ___ ^ , _, ___ ^ ____

^ л 2 л 2 \ n \ f n \ f J n n V / / /-ч 2

n ^) a n x nv ;ax2

где АП = (2Xn - sin2Xn + 4£Xn sin2 Xn) / (4Xn) — квадрат нормы собственных функций [14]; цn =

= цф n (1) / А П.

С учетом зависимостей (5) и (6) общее решение уравнения (13) запишется в виде

w,, (т) = Ane^ sin р„ т+ Bne^ cos р„ т + Cne^ + Dne™ + E,, (14)

2 где Pn =^Х2 - Д2 ; An, Bn — произвольные постоянные, определяемые из начальных условий (8); о Cn, D , En — коэффициенты, определяемые из равенства значений при соответствующих функциях S левой и правой части уравнения (13) (в силу громоздкости выражений значения коэффициентов A ,B , C ,D,E в статье не приводятся).

И Н -С В n ny ny n n Г /

W¡UJ Принимая во внимание ранее введенные обозначения, спектр частот собственных колебаний

определяется по формуле

Ю, =V„c /1. (15)

С учетом зависимостей (2), (9), (11) и (14) выражения для углов поворота, угловых скоростей и крутящих моментов произвольных сечений вала можно представить в виде

1 - е-°х^ х--

v(4> х) = £ ^П^ (( X Sin х + Bne~>-х cos p„ х + Cne ~2"х + Dne^ + ^) + <U -

n=1 K n C

ч с ^ sin X n %

l n=1 Xn

x(eAn* cosPnx(Pn^n - AnBn) - e~An* sinPnт(ДЛ + PnB) - 2aQe"2a* - a^e) + ёном (l - e); GJ

M(¡j,x) = --iXcosS x

l n=1

X ((T sin Pn т + Bne"fn' cos Pn т + Cne-2a< + Dne^ + E„).

; (16)

(17)

(18)

Динамический коэффициент примем как отношение крутящего момента в сечении вала к номинальному моменту сопротивления:

GJ.

Кд ($,х) = Rgf-Zcos*яnТ sinPnх + Bne"ДnТ cosрих + Cne"2ах + Dne^ + ^).

(19)

Результаты (Results)

В качестве примера рассмотрим механические характеристики БМРТ большого морозильного рыболовного траулера «Механик Ковтун»: диаметр гребного винта D = 4 м; момент инерции винта с присоединенной водой IM = 3226 кг-м2; номинальная частота вращения пном = 156 об/мин ( 0ном = 16,3 рад/с), номинальная мощность на винте Л"ном = 2650 кВт, длина валопровода (промежуточного и гребного валов) l = 11,35 м; приведенная жесткость валопровода GJp = 1,33 • 108 Н-м2; распределенный момент инерции валопровода у J = 13 кг • м (относительная инерционная нагрузка £ = IM / (y/pl j « 22), скорость распространения волн крутильных колебаний в стальном стержне c = 3200 м/с; по номинальной мощности и угловой скорости определяем коэффициент Rc:

R = Nhom / еНом = 587 кг-м2.

В результате решения уравнения (12) значения первых трех собственных чисел составят:

Х1 = 0,212; = 3,16; Х3 = 6,29.

Рассмотрим случай, когда винт выходит на 95 % номинальной скорости вращения за время, сопоставимое с периодом главной формы колебаний:

Гр = 2к! / (V) = 0,105 с.

Для этого случая параметры модели будут принимать значения: a = 0,1; Эср = 11,2 рад/с; ц = 0,158; Rc = 3,98 ; Д1 = 7,11 • 10"3 ; Р1 « 0,212 , для последующих форм колебаний Дп убывает,

а вп мало отличается от Xп. с

На графике (рис. 2, а) штриховой линией показан угол поворота левого торца вала в соот- 1

ветствии с законом движения (6), сплошной линией — угол поворота правого торца вала соглас- I

но зависимости (16). Аналогичные зависимости для угловых скоростей (5) и (17) представлены I

к

на графике (рис. 2, б). На рис. 3 приведена эпюра крутящих моментов для момента времени t = 0,053 с. Как видно из графика, изменение крутящих моментов по длине вала является незначительным. Зависимость для динамического коэффициента (19) на левом торце вала (£, = 0) представлена на графике (рис. 4, а). Для рассматриваемых параметров модели Kд ~ 11. Из приведенных на рис. 2-4 графиков следует, что для принятого режима разгона динамический эффект, являющийся весьма существенным, приведет к поломке оборудования, однако практически реализовать такой режим разгона крайне сложно.

2 О 2

n=1

а)

б)

Рис. 2. Углы поворота (а) и угловые скорости торцевых сечений валопровода (б)

Рис. 3. Эпюра крутящих моментов

а)

г

со

б)

см о

Рис. 4. Динамические коэффициенты при разгоне вала за время Т = 0,105 с (а) и Т = 2 с (б)

С увеличением времени разгона динамический эффект значительно убывает. Так, если осуществить выход на 95 % номинальной скорости за время в пределах, равных девятнадцати периодам главной формы колебаний (а = 0,0055; Тр = 2 с), то динамический коэффициент не превысит единицы (рис. 4, б).

ВЕСТНИК«!

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ^^

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

На рис. 5 показаны скорости в торцевых сечениях, при этом наложение кривых на графике указывает, что сечения вала в данном случае поворачиваются практически тождественно.

О 1 г, с 2

Рис. 5. Угловые скорости торцевых сечений валопровода для времени Т=2 с

Обсуждение (Discussion)

В данной работе исследования выполняются на основе упрощенной математической модели (4). Оценить адекватность принятых допущений для конкретных параметров модели можно сопоставляя значения величин по результатам расчета. Так, например, значения ф и (3 имеют примерно один порядок лишь в течение первого полупериода собственных колебаний, а в дальнейшем ф не превышает 0,1 (3. Однако для более объективной оценки погрешности, которую вносят принятые допущения, необходимы дополнительные теоретические исследования с использованием, например, численных методов расчета, или экспериментальные.

Еще одним упрощением в принятой модели является замена системы ступенчатых валов одним валом с эквивалентной податливостью, однако учитывая большую жесткость системы, что видно из эпюры крутящих моментов (см. рис. 3), данное допущение окажет незначительное влияние на результаты расчетов. Тем не менее при организации системы активного мониторинга работы оборудования для адекватного анализа информации, поступающей с датчиков, необходимо в дальнейшем оценить влияние ступенчатой формы валов на величину усилий и деформаций при динамических нагрузках.

Несмотря на некоторую громоздкость расчетных зависимостей, принятый в работе подход дает универсальный и эффективный метод исследования напряженно-деформированного состояния валопроводов при различных динамических воздействиях, таких как выход винта из воды, попадание в неоднородный набегающий поток, удар винта о лед, резонансные явления при пульсирующей работе двигателя и др. Данную математическую модель также можно дополнить членами, учитывающими внешнее и внутреннее трение, что является актуальным для исследований резонансных явлений и учета потерь от колебательного движения. 2

Заключение (Conclusion)

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

- динамические усилия в валопроводе существенны и могут многократно превышать значения от номинальной нагрузки, если время вывода на номинальную частоту вращения ГД сопоставимо с несколькими периодами главной формы колебаний;

- собственные колебания значительно проявляются лишь в течение короткого промежутка времени (до десяти периодов главной формы колебаний) после динамического воздействия,

N

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

■Р

ЛВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

а затем быстро затухают за счет сил сопротивления на винте, даже без учета потерь на внутреннее и внешнее трение в валопроводе;

- с повышением плавности ввода ГД в рабочий режим (время разгона более двадцати периодов главной формы колебаний) динамический эффект становится незначительным, не создавая дополнительных усилий в валах.

СПИСОК ЛИТЕРАТУРЫ

1. Судовой механик: справ. / Под общ. ред. А. А. Фока. — Одесса: Фешкс, 2008. — Т. 1. — 1036 с.

2. Ларин А. А. Роль исследований крутильных колебаний валопроводов в развитии динамики машин / А. А. Ларин // Питання гстори науки i техшки. — 2009. — № 4. — С. 2-9.

3. Коврижных М. Н. Расчет амплитуд свободных колебаний дискретных многомассовых систем / М. Н. Коврижных, С. С. Глушков // Научные проблемы транспорта Сибири и Дальнего Востока. — 2008. — № 2. — С. 162-164.

4. Soni T. Active vibration control of ship mounted flexible rotor-shaft-bearing system during seakeeping / T. Soni, A. S. Das, J. K. Dutt //Journal of Sound and Vibration. — 2020. — Vol. 467. — Pp. 115046. DOI: 10.1016/ j.jsv.2019.115046.

5. Мартьянов В. В. Оценка угрозы возникновения резонансных колебаний на примере расчета крутильных колебаний судового валопровода пассажирского теплохода пр. Р118 / В. В. Мартьянов // Вестник Государственного университета морского и речного флота имени адмирала С. О. Макарова. — 2020. — Т. 12. — № 2. — С. 359-368. DOI: 10.21821/2309-5180-2020-12-2-359-368.

6. Jee J. Design Improvement of a Viscous-Spring Damper for Controlling Torsional Vibration in a Propulsion Shafting System with an Engine Acceleration Problem / J. Jee, C. Kim, Y. Kim // Journal of Marine Science and Engineering. — 2020. — Vol. 8. — Is. 6. — Pp. 428. DOI: 10.3390/jmse8060428.

7. Murawski L. Shaft line whirling vibrations: effects of numerical assumptions on analysis results / L. Murawski // Marine Technology and SNAME News. — 2005. — Vol. 42. — Is. 2. — Pp. 53-60. DOI: 10.5957/mt1.2005.42.2.53.

8. Halyavkin A. Determination of stiffness coefficient of stern shaft bearing / A. Halyavkin, I. Razov, V. Mamon-tov, G. Kushner // IOP Conference Series: Earth and Environmental Science. — IOP Publishing, 2017. — Vol. 90. — Is. 1. — Pp. 012078. DOI: 10.1088/1755-1315/90/1/012078.

9. Екимов А. С. Измерение крутильных колебаний валопроводов судовых энергетических установок теплоходов проекта 1741А / А. С. Екимов, Н. Г. Петрова, С. С. Глушков // Научные проблемы транспорта Сибири и Дальнего Востока. — 2009. — № 2. — С. 228-230.

10. Тхыонг Н. М. Современное состояние вопроса о динамике судовых валопроводов (опыт Вьетнама) / Н. М. Тхыонг // Судовые энергетические установки: науч.-техн. сб. — 1998. — № 2. — С. 95-99.

11. Мамонтов В. А. Автоматизация методики расчета собственной частоты поперечных колебаний гребного вала / В. А. Мамонтов [и др.] // Вестник Астраханского государственного технического университета. Серия: Морская техника и технология. — 2019. — № 1. — С. 63-70. DOI: 10.24143/2073-1574-2019-1-63-70.

12. Сапига В. В. Совершенствование методов анализа динамики судовых валопроводов / В. В. Сапига, А. Л. Кирюхин, П. В. Черпита // Водный транспорт. — 2014. — № 1 (19). — С. 52-61.

13. Улитин Г. М. Изгибные колебания стержня с переменной жесткостью и распределенной массой / г Г. М. Улитин, С. Н. Царенко // Прикладная математика и механика. — 2015. — Т. 79. — № 6. — С. 817-823.

14. Царенко С. Н. Крутильные колебания стержневых конструкций с осевой неоднородностью гео-Ш метрических характеристик / С. Н. Царенко // Вестник Южно-Уральского государственного университета. а Серия: Математика. Механика. Физика. — 2019. — Т. 11. — № 1. — С. 50-58. DOI: 10.14529/mmph190107.

15. Han H. S. Estimate of the fatigue life of the propulsion shaft from torsional vibration measurement and

CM

the linear damage summation law in ships / H. S. Han, K. H. Lee, S. H. Park // Ocean Engineering. — 2015. — ■ПП6 Vol. 107. — Pp. 212-221. DOI: 10.1016/j.oceaneng.2015.07.023.

16. Huang Q. Numerical modeling and experimental analysis on coupled torsional-longitudinal vibrations of a ship's propeller shaft / Q. Huang, X. Yan, Y. Wang, C. Zhang, Z. Wang // Ocean Engineering. — 2017. — Vol. 136. — Pp. 272-282. DOI: 10.1016/j.oceaneng.2017.03.017.

17. Huang Q. Coupled transverse and torsional vibrations of the marine propeller shaft with multiple impact factors / Q. Huang, X. Yan, C. Zhang, H. Zhu // Ocean Engineering. — 2019. — Vol. 178. — Pp. 48-58. DOI: 10.1016/ j.oceaneng.2019.02.071.

18. Рак А. Н. Особенности расчета параметров вспомогательно-аварийного движения на морских судах с комбинированным дизель-электрическим движителем / А. Н. Рак, С. Н. Царенко // Вестник Государственного университета морского и речного флота имени адмирала С. О. Макарова. — 2021. — Т. 13. — № 1. — С. 115-125. DOI: 10.21821/2309-5180-2021-13-1-115-125.

19. Филиппов А. П. Колебания механических систем / А. П. Филиппов. — Киев: Наукова думка, 1965. —

716 с.

REFERENCES

1. Fock, A. A., ed. Sudovoy mekhanik: sprav. Vol. 1. Odessa: Feniks, 2008.

2. Larin, A. A. "Rol' issledovaniy krutil'nykh kolebaniy valoprovodov v razvitii dinamiki mashin." Pytannya istoriyi nauky i tekhniky 4 (2009): 2-9.

3. Kovrizhnykh, M. N., and S. S. Glushkov. "Raschet amplitud svobodnykh kolebaniy diskretnykh mnogo-massovykh system." Nauchnyyeproblemy transporta Sibiri i Dal'nego Vostoka 2 (2008): 162-164.

4. Soni, Tukesh, A. S. Das, and J. K. Dutt. "Active vibration control of ship mounted flexible rotor-shaft-bearing system during seakeeping." Journal of Sound and Vibration 467 (2020): 115046. DOI: 10.1016/j.jsv.2019.115046.

5. Martianov, Vladimir V. "Assessment of synchronous vibrations risks as exemplified in the shafting torsional vibrations analysis for R118 project passenger ship." Vestnik Gosudarstvennogo universiteta morskogo i rechnogo flota imeni admirala S. O. Makarova 12.2 (2020): 359-368. DOI: 10.21821/2309-5180-2020-12-2-359-368.

6. Jee, Jaehoon, Chongmin Kim, and Yanggon Kim. "Design Improvement of a Viscous-Spring Damper for Controlling Torsional Vibration in a Propulsion Shafting System with an Engine Acceleration Problem." Journal of Marine Science and Engineering 8.6 (2020): 428. DOI: 10.3390/jmse8060428.

7. Murawski, Lech. "Shaft line whirling vibrations: effects of numerical assumptions on analysis results." Marine Technology andSNAME News 42.2 (2005): 53-60. DOI: 10.5957/mt1.2005.42.2.53.

8. Halyavkin, Aleksey, Igor Razov, Victor Mamontov, and Guriy Kushner. "Determination of stiffness coefficient of stern shaft bearing." IOP Conference Series: Earth and Environmental Science. Vol. 90. No. 1. IOP Publishing, 2017. DOI: 10.1088/1755-1315/90/1/012078.

9. Ekimov, A. S., N. G. Petrov, and S. S. Glushkov. "Torsional oscillation in the propeller shafts of ship engine 1741A." Nauchnyye problemy transporta Sibiri i Dal'nego Vostoka 2 (2009): 228-230.

10. Thuong, N. M. "Sovremennoye sostoyaniye voprosa o dinamike sudovykh valoprovodov (opyt V'yetnama)." Sudovyye energeticheskiye ustanovki: nauchn.-tekhn. sb. Vol. 2. 1998. 95-99.

11. Mamontov, Victor Andreevich, Ibragim Musaevich Abacharaev, Vladimir Pavlovich Bulgakov, and Guriy Alekseevich Kushner. "Automating the method of calculating natural frequency of propeller shaft transverse vibrations." Vestnik of Astrakhan State Technical University. Series: Marine Engineering and Technologies 1 (2019): 63-70. DOI: 10.24143/2073-1574-2019-1-63-70.

12. Sapiga, V., A. Kiryukhin, and P. Cherpita. "Perfection methods for analyzing the dynamics of marine shafting." Vodnyy transport 1(19) (2014): 52-61.

13. Ulitin, G.M., and S. N. Tsarenko. "Flexural vibrations of a rod of variable stiffness and distributed mass." Journal of Applied Mathematics and Mechanics 79.6 (2015): 572-576. DOI: 10.1016/j.jappmathmech.2016.04.003.

14. Tsarenko, S. N. "Torsional vibrations of a rod with an axial geometric inhomogeneity." Bulletin of the South Ural State University», series «Mathematics. Mechanics. Physics 11.1 (2019): 50-58. DOI: 10.14529/ mmph190107.

15. Han, HyungSuk, KyungHyun Lee, and SungHo Park. "Estimate of the fatigue life of the propulsion shaft from torsional vibration measurement and the linear damage summation law in ships." Ocean Engineering 107 (2015): 212-221. DOI: 10.1016/j.oceaneng.2015.07.023. |

16. Huang, Qianwen, Xinping Yan, Yikun Wang, Cong Zhang, and Zhihua Wang. "Numerical modeling and -experimental analysis on coupled torsional-longitudinal vibrations of a ship's propeller shaft." Ocean Engineering f 136 (2017): 272-282. DOI: 10.1016/j.oceaneng.2017.03.017. \

17. Huang, Qianwen, Xinping Yan, Cong Zhang, and Hanhua Zhu. "Coupled transverse and torsional vibrations j of the marine propeller shaft with multiple impact factors." Ocean Engineering 178 (2019): 48-58. DOI: 10.1016/j. oceaneng.2019.02.071.

18. Rak, Alexander N., and Sergey N. Tzarenko. "Features of calculating the parameters of auxiliary-emergency traffic on sea-going ships with combined diesel-electric propulsor." Vestnik Gosudarstvennogo universiteta morskogo i rechnogo flota imeni admirala S. O. Makarova 13.1 (2021): 115-125. DOI: 10.21821/23095180-2021-13-1-115-125.

19. Filippov, A. P. Kolebaniya mekhanicheskikh sistem. Kyev: Naukova dumka, 1965.

2 О 2

ЛВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

ИНФОРМАЦИЯ ОБ АВТОРАХ

INFORMATION ABOUT THE AUTHORS

Царенко Сергей Николаевич —

кандидат технических наук, доцент Камчатский государственный технический университет

683003, Российская Федерация, Петропавловск-Камчатский, ул. Ключевская, 35 e-mail: tzarenko@gmail.com Рак Александр Николаевич — кандидат технических наук

Донецкий национальный технический университет

83001, Украина, Донецк, ул. Артема, 58

e-mail: Alexander.Rak@mail.ru

Безлобенко Борис Николаевич —

начальник комплексного судового производства

АО «Петропавловск-Камчатский

судоремонтный завод»

683015, Российская Федерация, Петропавловск-Камчатский, ул. Днепровская, 1 e-mail: bezpsrz@mail.ru

Tzarenko, Sergey N. —

PhD, associate professor Kamchatka State Technical University

35 Klyuchevskaya Str., Petropavlovsk-Kamchatsky,

683003, Russian Federation

e-mail: tzarenko@gmail. com

Rak, Alexander N. —

PhD

Donetsk National Technical University

58 Artyoma Str., Donetsk, 83001, Ukraine

e-mail: Alexander.Rak@mail.ru

Bezlobenko, Boris N. —

Head of the complex marine production

Petropavlovsk-Kamchatsky Shipyard,

JSC

1 Dneprovskaya Str., Petropavlovsk-Kamchatsky, 683015, Russian Federation e-mail: bezpsrz@mail.ru

Статья поступила в редакцию 3 июня 2021 г.

Received: June 3, 2021.

i Надоели баннеры? Вы всегда можете отключить рекламу.