Об одной из возможных причин динамической неустойчивости валопроводов судов Текст научной статьи по специальности «Механика и машиностроение»

УДК 629.12.037.001.5

А. И. Миронов, А. Ю. Кукарина

ОБ ОДНОЙ ИЗ ВОЗМОЖНЫХ ПРИЧИН ДИНАМИЧЕСКОЙ НЕУСТОЙЧИВОСТИ ВАЛОПРОВОДОВ СУДОВ

Одной из возможных причин динамической неустойчивости вала, возникающей в процессе эксплуатации валопроводов судов, являются износ дейдвудного подшипника и периодически изменяющиеся под действием поперечных колебаний вала реакции подшипников. Рассматривается кормовой участок гребного вала как имеющий наименьшую собственную частоту. Участок гребного вала моделируется балкой постоянного сечения, опирающейся на смещаемую опору одностороннего действия. Так как в процессе поперечных колебаний балки возможен ее отрыв от опоры, крайняя опора принята защемляющей. В решении учитываются масса вала, а также масса и момент инерции массы гребного винта. При заданных параметрах гребного вала полученное решение позволяет определить собственную частоту колебаний вала, при заданной скорости вращения вала - форму колебаний и суммарные нагрузки на вал от действия статических и динамических нагрузок, а также влияние износа дейдвудного подшипника. Численный эксперимент для валопровода конкретного рыбодо-бывающего судна показал, что поперечные колебания оказывают существенное влияние на характер вращения вала. Начиная с некоторого значения величины износа подшипника его реакция становится периодически отрицательной, т. е. вал будет периодически отрываться от подшипника, и, как следствие, его вращение будет становиться неустойчивым.

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

Введение

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

Динамические нагрузки наиболее опасны. Известно, что небольшие по величине, но периодически воздействующие нагрузки могут вызвать разрушение конструкций, построенных с большим запасом статической прочности.

Одним из механизмов судна, подвергающихся воздействию весьма сложной системы эксплуатационных нагрузок, является валопровод. Любой отказ в работе валопровода чреват серьезными последствиями, вплоть до гибели судна. Известны случаи, когда из-за сильной вибрации вала судно теряло ход, и его дальнейшее движение без буксира становилось невозможным. Отметим, что сильные вибрации начинались без видимых причин. Вал устойчиво работал достаточно длительный отрезок времени, но в какой-то момент начинались сильные вибрации. Известны случаи, когда сильная вибрация начиналась после «центровки» валопровода при запуске двигателя.

Постановка задачи

Причины, приводящие к неустойчивой работе валопровода, могут быть следующими:

— вследствие износа подшипников (особенно дейдвудных) в процессе эксплуатации судна изменяется укладка вала и, как следствие, характер колебаний вала (могут возникнуть параметрические колебания);

— при суммировании реакций от весовых нагрузок на вал (статические реакции при «центровке» вала) и реакций, вызванных колебаниями вала, может произойти отрыв вала от подшипника;

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

Нами рассматривается задача о влиянии износов подшипников валопровода на их суммарные реакции.

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

1 2 q EI 3

i ^ i i i >|r i i ^Г^ч"

F

f

S- B

шт.

\

L

)

М

ГДМ

l

z

Рис. 1. Расчетная схема кормового участка гребного вала: 1 - вал до износа опоры; 2 - вал после износа опоры (пунктирная линия); 3 - гребной винт; f - величина износа опоры; МГдМ - переменный ГДМ; F - вес винта; q - погонный вес гребного вала; EI - жесткость сечения вала на изгиб; l, L - длины участков вала

Так как на вал действуют статические нагрузки от веса вала и гребного винта и переменный гидродинамический момент Мгдм, рассматриваем нагружение вала как сумму статических и динамических нагрузок независимо. При этом, поскольку ГДМ имеет постоянную Мт и переменную Ma sin mt составляющие (Мгдм = Мт + Ma sin т), постоянную составляющую Мт относим к статическим нагрузкам (рис. 2).

(

Ma ( I

б

Ma

It

Rac л

A

......ARbc

^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^

i

i

R

ad

A

ШШ

m

mRbd

~Bmrn

F

\

L

)

F

-*- И1

Mm,

M,I

Ми

Ma sin ш

a

l

Рис. 2. Расчетные схемы кормового участка гребного вала: а - при действии постоянных нагрузок; б - при действии переменной составляющей гидродинамического момента; m, M - соответственно погонная масса гребного вала и гребного винта; I - момент инерции винта

В результате получим:

Ra = Rac + Rad,

Rb = Rbc + Rbd, (1)

Ma = Mac + Mad.

Математическое решение задачи

Начнем со статического нагружения (рис. 2, а).

Так как система статически неопределимая, то для определения реакций используем универсальное уравнение изогнутой оси балки: - 0 < г < I:

У = Уо + Фо z +

м о z2 + Qo z3

qz

2EI 6EI 24EI

l < z < L:

У = уо+ф0 z+MoL+Qozi+Rbc (z -1 )3 qz 4

2Е1 6Е1 6Е1 24 Е1 Так как начальные параметры

У0 = 0; фо = 0; Мо = МАС; йо = Яас,

то выражения (2) и (3) принимают следующий вид: — 0 < г < I:

ил 2 о 3 4

У - масг , яасг дг .

У —--+----;

2Е1 6Е1 24Е1

l < z < L:

= MAcz2 , Racz3 , RBC(z -1)3 qz4

у =—AC—+

2EI 6EI

+

6EI

24EI

(2)

(3)

(4)

(5)

Неизвестные реакции ЯАС, ЯВС, МАС, входящие в выражения (4) и (5), находим из следующих граничных условий:

г = I у = - /,

г = Ь Мх = Мт, (6)

г = Ь йу = К

С учетом того, что

d2 у

d3 у

mx = ^. Q = EIdz3

граничные условия (6) принимают вид

M acl2 + RacI3

ql4

2EZ

6EI 24EI

= - f,

L2

mac + racl + rbc (l - l) - ^ = m

2

(7)

rac + rbc - ql = f-

Из системы уравнений (7) находим значения статических реакций ЯАС и Я

вс:

rac =— 2l

f

fo + Mm - F(L -1) + q

L l ^ Ll - — - —

v 2 12 У

Rbc 2l

fo + Mm - F(L - 31) - q

f 1„ L2 l Л — Ll + — + — 3 2 12

где

г -

J0 - 12 ■> .

Решение динамической части задачи (рис. 2, б). Дифференциальное уравнение колебаний балки:

— 0 < г < Ь:

д4 у д2 у

Е1 —у + т —у — 0. (9)

дг4 дГ

Принимаем

у - у(г)ап юг. (10)

Подставляем (10) в уравнение (9). После интегрирования полученного уравнения для форм колебаний балки находим следующие выражения:

— 0 < г < I:

у — у0К1 (аг) + ф0К2 (аг) + -М°-К3 (аг) + К4 (аг); (11)

а а2 EI a3EI

l < z < L:

у = y0K (аz) + K2 (az) + ^K (az) + а a EI

+ K4 (аz) + -Rm^K4 (а(z-1)), (12)

а3 EI а EI

где у0, ф0, M0, Q0 - начальные параметры, определяются из граничных условий; Кь К2, К3, К4 - система фундаментальных функций с единичной матрицей аргументов (az) и (a (z - l)) [1, c. 294].

2

а 4 = ^ (13)

EI

Граничные условия:

z = 0 у = 0 и ф = 0, т. е. уо = 0 и фо = 0; Mo = Mad, Qo = Rad. (14)

z = l у = 0. (15)

z = L Qу = - FHH; Мх = Мин + Ma sin юг, (16)

где Qу и Мх - соответственно поперечная сила и изгибающий момент в сечении на правом конце балки; и Мин - соответственно сила инерции и изгибающий момент от момента инерции I массы гребного винта (рис. 2, б):

F-=- M W •

м =- i д2ф=-I Г ду м " ' dt1 ' Эг I dz

Так как

Э2 у Э3 у

м = EI эу, Q = EI э? ■

граничные условия (16) принимают следующий вид:

- z = L:

EI0 = M 0,

EI —y = -1 —^ I — I + Ma sin mt.

dz д t2 l dz J a

Учитывая условия (14), выражения для форм колебаний (11) и (12) следующие:

0 < z < l:

(17)

(18)

M

y = K3(аz(аz); а2 EI aEI

l < z < L:

(19)

y = Щш. K3 (а z )+4^К4 (а z H^K (а (z -1)). а2EI v ' а3Е1 v ' а3Е1 v '

(20)

Окончательно, используя (19) и (20), после преобразований граничные условия (15), (17) и (18) получаем в следующей форме:

A11mad + A12 rad = °

A21M ad + A22 rad + A23rbd =

A31mad + A32rad + A33rbd = ma.

Коэффициенты Aj (i = 1-3, j = 1-3) имеют следующие значения:

в1

A11 = ^K3(ß1); A12 = K4(ß1); A13 = 0;

A21 = ß[K4(ß) + Yj3K3(ß)]; A22 =[K1(ß) + Y1ßK4(ß)];

A23 =[K1(ß2) + Y1ßK4(ß2)]; A31 = [K1(ß)-Y2ß3K2(ß)];

3

A32 = L[K2(ß)-Y2ß3K3(ß)]; A33 = L[K2(ß2)-Y2ß3K3(ß2)];

М I

в = ®Е; р1 = а1; р2 = а(Е-0; ^ = —; Y2 = —тг.

шЕ тЕ

Используя систему уравнений (21), можно решать две задачи. 1. При заданных параметрах гребного вала определять его собственную частоту. Для этого из условия, что определитель системы уравнений (21) А равен нулю, т. е.

A11 A12 А13

А = A21 A22 A23 =0,

A31 A32 A33

(21)

находим величину р. Затем из выражения (13) вычисляем собственную частоту гребного вала:

ю =

2

vLy

2. При заданной скорости вращения вала ю из системы уравнений (21) определяем динамические составляющие реакций Мао, КАО и Кво.

Численный эксперимент

Оценку влияния износов подшипника f на величину параметров «центровки» валопровода произведем путем численного эксперимента.

Рассмотрим кормовой участок гребного вала со следующими параметрами [2, с. 69]: Ь = 6,11 м; I = 4,385 м; ^ = 73,85 кН; q = 11,05 кН/м; Е1 = 8,064-105 кН- м2; лопастная частота ю = 60,81 1/с, М™ = -138,2 кН/м; МГДМ = -37,69 кН/м; Мт = -87,94 кН/м; Ма = 50,26 кН/м.-

Момент инерции гребного винта, кг ■ м2, вычислялся по формуле [2, с. 97]:

Га ^

2 А

I = 0,02М • Б-

А

А

-+ 3

V Аб

Для рассматриваемого валопровода: - диаметр гребного винта Б = 3,70 м;

А

- дисковое отношение гребного винта-= 0,57.

ао

Для динамического расчета использовалась Программа для ЭВМ [3]. В результате вычислений находим:

- собственная частота гребного вала ю0 = 125 1/с; «отстройка» вала Аю от лопастной частоты ю составляет [2, с. 98]:

Аю =

1

ю

1 -

125

60,8

= 1,056 > 0,2.

Таким образом, условие «отстройки» вала от резонанса выполняется с большим запасом. Вычисляем реакции.

Динамические составляющие реакций вычисляем по выражениям (21), используя Программу для ЭВМ [3].

В результате находим:

МАБ ± 34,3 кН-м; =± 23,63 кН-м; Явв =± 31,68 кН-м.

Статическую составляющую реакции Ево вычисляем, используя выражение (8):

ЯВС = - 2,870-104 кН/м / + 190,4 кН/м. Полная реакция на кормовом дейдвудном подшипнике равна (см. выражение (1)):

яв = ЯвС + Яво =-2,870- 104 кН/м / +190,4 кН/м ± 31,68 кН/м = = -2,870 • 104 кН/м / +252,7 кН/м.

На рис. 3 изображена зависимость реакции на кормовом дейдвудном подшипнике от величины его износа и поперечных колебаний вала.

240 200 160 120 80 40 0 -40

>>

ш ш

1 2 3 4

Рис. 3. Изменение реакции на кормовом дейдвудном подшипнике в зависимости от его износа и поперечных колебаний вала - заштрихованная область; значение статической составляющей реакции - сплошная линия

Согласно рис. 3, поперечные колебания вала оказывают на его движение существенное влияние. Начиная с некоторого значения величины износа подшипника f его суммарная реакция становится периодически отрицательной, тогда как статическая составляющая остается положительной. При одностороннем взаимодействии вала с подшипником это означает, что вал будет периодически отрываться от подшипника, т. е. его вращение становится динамически неустойчивым.

Выводы

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

Следует отметить, что численные значения результатов исследования зависят от конкретного валопровода. Однако качественно они соответствуют полученным в [4].

Таким образом, при оценке динамической устойчивости валопровода недостаточно только проверять его на возможность резонанса. При расчете параметров «центровки» необходимо учитывать поперечные колебания вала и их влияние на величину реакций подшипников.

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

1. Миронов А. И. К исследованию поперечных колебаний гребных валов. Ч. 3. Влияние момента инерции винта на собственную частоту и форму колебаний гребного вала / А. И. Миронов // Вестн. Аст-рахан. гос. техн. ун-та. Сер.: Морская техника и технология. 2014. № 3. С. 21-27.

2. Прочность, устойчивость, колебания: Справочник: в 3 т. / под общ. ред. И. А. Биргера, Я. Г. Пановко. М.: Машиностроение, 1968. Т. 3. 568 с.

3. Миронов А. И. Программа расчета поперечных колебаний гребного вала с учетом момента инерции гребного винта / А. И. Миронов, Е. В. Пономарёва, А. Ю. Кукарина. Свид-во об офиц. регистрации № 2015612797; заяв. № 2014664032; зарег. в реестре программ для ЭВМ 26.02.2015.

4. Комаров В. В. Валопроводы рыбопромысловых судов. Ч. 1. / В. В. Комаров, А. С. Курылев. Астрахань: Изд-во АГТУ, 1997. 166 с.

Статья поступила в редакцию 06.07.2015, в окончательном варианте - 11.07.2015

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

Миронов Альфред Иванович - Россия, 414000, Астрахань; канд. техн. наук, доцент; независимый исследователь; [email protected].

Кукарина Алла Юрьевна - Россия, 414056, Астрахань; Астраханский государственный технический университет; старший преподаватель кафедры «Механика и инженерная графика»; [email protected].

A. I. Mironov, A. Yu. Kukarina

ON ONE OF THE POSSIBLE CAUSES OFDYNAMIC INSTABILITY OF SHIP SHAFTING

Abstract. One of the possible causes of the dynamic instability of shafting, arising during the operation of ship shaftings, is wear of stern-tube bearing and the reactions of stern-tube bearings periodically changing under the influence of transverse vibrations of the shaft. The aftermost part of the propeller shaft is considered as one that has the lowest natural frequency. The part of the propeller shaft is modeled with a beam of constant cross section, which is based on the biased support of the unilateral action. Since in the process of transverse vibrations of the beam the detach-

ment from the support is possible, an extreme support is accepted as a restrained one. The decision takes into account the mass of the shaft, as well as the mass and the moment of inertia of the propeller mass. At the given parameters of the propeller shaft the obtained decision helps determine the natural frequency of the shaft, at the given speed of rotation - waveform and the total load on the shaft from the effects of static and dynamic loads, as well as the impact of stern-tube bearing wear. The numerical experiment for shafting of the particular fishing ship showed that the transverse vibrations have a significant influence on the character of rotation. Starting with a certain value of the bearing wear, its response from time to time becomes negative, i. e., shaft will periodically break away from the bearing and, as a consequence, its rotation will become unstable.

Key words: dynamic instability, shafting, transverse vibrations, static and dynamic loading.

REFERENCES

1. Mironov A. I. K issledovaniiu poperechnykh kolebanii grebnykh valov. Chast' 3. Vliianie momenta in-ertsii vinta na sobstvennuiu chastotu i formu kolebanii grebnogo vala [To the study of transverse vibrations of the propeller shafts]. Vestnik Astrakhanskogo gosudarstvennogo tekhnicheskogo universiteta. Seriia: Morskaia tekhnika i tekhnologiia, 2014, no. 3, pp. 21-27.

2. Prochnost', ustoichivost', kolebaniia: Spravochnik: v 3 t. [Strength, stability, vibrations: Reference: in 3 vol.]. Pod obshchei redaktsiei I. A. Birgera, Ia. G. Panovko. Moscow, Mashinostroenie Publ., 1968. Vol. 3. 568 p.

3. Mironov A. I. Ponomareva E. V., Kukarina A. Iu. Programma rascheta poperechnykh kolebanii grebnogo vala s uchetom momenta inertsii grebnogo vinta [Procedure of the calculation of the transverse vibrations of the propeller shaft taking into account the inertia of the propeller screw]. Svidetel'stvo ob ofitsial'noi registratsii № 2015612797; zaiavka № 2014664032; zaregistrirovano v reestre programm dlia EVM 26.02.2015.

4. Komarov V. V., Kurylev A. S. Valoprovody rybopromyslovykh sudov. Part 1 [Shafting of the fishing vessels]. Astrakhan, AGTU Publ., 1997. 166 p.

The article submitted to the editors 06.07.2015, in the final version - 11.07.2015

INFORMATION ABOUT THE AUTHORS

Mironov Alfred Ivanovich - Russia, 41400, Astrakhan; Candidate of Technical Sciences, Assistant Professor; Independent Researcher; [email protected].

Kukarina Alla Yurievna - Russia, 414056, Astrakhan; Astrakhan State Technical University; Candidate of Technical Sciences; Senior Lecturer of the Department "Mechanics and Engineering Graphics"; [email protected].