.

Динамика двигателя судна, условия и оценка работы

Динамика двигателя включает в себя обширную информацию про движение КШМ под действием изменяющегося давления и сил инерции, что и рассмотренно ниже.


В процессе работы двигателя его детали находятся под воздействием:

а) переменного давления газов в цилиндре;

б) сил инерции поступательно движущихся масс и моментов от этих сил;

в) центробежных сил вращающихся масс и моментов от этих сил;

г) крутящего момента двигателя и противоположного ему опрокидывающего момента;

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

Суммарное влияние этих факторов определяет “динамику” двигателя — его механическую напряженностьИзменение механической напряженности, характер изменения давления в подшипниках коленчатого вала и условия смазки подшипников, неравномерность вращения вала, вибрацию двигателя и корпуса судна.

Силы, определяющие механическую напряженность двигателя

Строго говоря, механические напряжения в деталях двигателя вызываются не только воздействием перечисленных выше динамических факторов, но и:

  1. Статическими усилиями затяга деталей крепежа при сборке дизеля;
  2. Изменением “внешних” условий — деформацией корпуса судна при его различной осадке или при волнении моря.

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

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

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

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

На кривошипно-шатунный механизм действует 4 рода сил: давления газов Pг, силы инерции Pj, силы веса Pg и силы трения Pтр (рис. 1). Сила от давления газов на поршень равна: Pг = pгFn.

где:

  • pг – давление газов в цилиндре;
  • Fn – площадь поршня.

Эта сила приложена к поршневому пальцу — в тронковых двигателях — или к поперечине крейцкопфа — в крейцкопфных двигателях.

Кривошипно-шатунный механизм
Рис. 1 Силы, действующие в кривошипно-шатунном механизме

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

Pj = Msjs = Gs/g js,         Форм. 1

где:

  • Ms = Gs/g – масса поршня и поступательно движущейся части шатуна;
  • g – ускорение свободного падения;
  • js – ускорение поступательно движущихся масс.

Ускорение js может быть найдено из уравнения перемещения поршня S = f(τ) (обозначения см. рис. 1):

S = OAo  OAi = OAo  OC + CA = R + L  Rcosφ + Lcosβ = R + LR/R  Rcosφ + L/R cosβ.

где:

  • OAo = L + R:
  • L – длина шатуна;
  • R – радиус кривошипа;
  • φ – угол поворота кривошипа от ВМТ;
  • β – угол между шатуном и осью цилиндра.

Обозначим: R/L = λ; тогда:

S = R 1 + 1 / λ  cosφ + 1 / λ cosβ.         Форм. 2

Так как:

Rsinφ = Lsinβ, R/L = sinβ/sinφ = λ; sinβ = λsinφ;

то:

cosβ = 1  Sin2 β = 1  λ2·Sin2φ.

Раскладывая равенство cosβ в биноминальный ряд и взяв первые 2 члена разложения, получим приближенное равенство:

cosβ  1  1/2 λ2 sin2 φ.        Форм. 3

Подставим это равенство в формулу 2. Получим:

S  R 1 + 1/λ  cosφ  1/λ + 1/λ 1/2 λ2 sin2 φ = R 1  cosφ + 1/2 λ sin2 φ.

Так как:

cos2φ = cos2φ  sin2φ = 1  2 sin2φ;

sin2φ = 1/2 1  cos2φ,

то:

S = R 1  cosφ + λ/4 1  cos2φ.        Форм. 4

Взяв производные от пути по времени, найдем законы изменения скорости и ускорения движения поршня:

C = dS/dτ = dS/dφ dφ/dτ = Rω sinφ + λ/2sin2φ,         Форм. 5

J = dC/dτ = dC/dφ dφ/dτ = Rω2 cosφ + λcos2φ,         Форм. 6

где:

  • ω = dφ/dτ – угловая скорость коленчатого вала.

Тогда сила инерции Pj определится зависимостью:

Pj = Ms Rω2 cosφ + λcos2φ.

Первое слагаемое этого равенства называется силой инерции поступательно движущихся масс 1-го порядка:

Pj1 = Ms Rω2 cosφ.        Форм. 7

Второе слагаемое называется силой инерции II-го порядка:

Pj11 = Ms Rω2 λ cos2φ.        Форм. 8

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

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

Характер изменения величин

 Pj1, Pj11

дан на рис. 2. Как видно из рисунка, период колебания величины

 Pj1

равен периоду вращения коленчатого вала, а период колебания

 Pj11

— в 2 раза больше (за 1 оборот коленчатого вала сила

 Pj11

дважды меняет свою величину и на-правление). Сила инерции 1-го порядка имеет максимум в ВМТ и НМТ. При нахождении поршня в этих положениях скорость его равна 0, а ускорение — максимально.

Силы инерции движущихся масс в ДВС
Рис. 2 Характер изменения сил инерции поступательно движущихся масс

В ВМТ силы инерции действуют против направления сил от давления газов, стремятся оторвать поршень вверх. В НМТ, наоборот, — при остановке поршня силы инерции действуют вниз, по направлению действия сил давления газов. При повороте коленчатого вала на 90° от ВМТ сила инерции 1-го порядка равна нулю:

Pj1 = 0.

Силы инерции 2-го порядка как в ВМТ, так и в НМТ равны по величине и по направлению и действуют против сил от давления газов на поршень. Максимальные значения силы

 Pj11
, совпадающие с направлением действия сил давления газов, наблюдаются при φ = 90° и 270°.

Причиной появления сил инерции 2-го порядка является конечная длина шатуна. Так, если принять длину шатуна бесконечной: L = ∞, то λ = R / L = 0; соответственно

 Pj11 = 0.

В реальном двигателе из-за конечной длины шатуна при повороте коленчатого вала на 90° после ВМТ поршень проходит более половины своего хода. Эта разница может быть найдена с помощью зависимости 4:

S = Sφ=90  1/2 S = R 1 + 1/4 λ 2  R = 1/2 λ R = R2 / 2L.        Форм. 9

Величина ΔS = R2/2L носит название “поправка профессора Брикса Ф. А. на конечную длину шатуна”. Она указывает, насколько поршень проходит более 1/2 своего хода при повороте кривошипа на 90° от ВМТ. Следовательно, в диапазоне φ = 0 + 90° пкв средняя скорость поршня больше, чем при φ = 90 + 180° пкв. Это обстоятельство оказывает свое влияние на соответствующее изменение ускорения поршня и появление силы инерции 2-го порядка.

Центробежные силы инерции вращательно движущихся масс кривошипа равны:

Pjц = MR jц = MR Rω2,        Форм. 10

где:

  • MR – вращающиеся массы;
  • jц = Rω2 – центробежное ускорение.

Механические напряжения, возникающие в днище головки поршня,определяются только воздействием сил от давления газов. Условия же работы головного соединения определяются уже не только давлением газов, но и силами инерции поршня, его весом и силами трения между поршнем и втулкой цилиндра. Для оценки механической напряженностиИзменение механической напряженности головного, мотылевого соединений и шатуна вводится понятие суммарной движущей силы PΣ, равной алгебраической сумме сил от давления газов Pг, сил инерции поступательно движущихся масс Pj, сил веса поступательно движущихся масс Pg и сил трения Pтр:

PΣ = Pг + Pj + Pg + Pтр.         Форм. 11

Приложив силу PΣ в центре головного подшипника, можно найти нормальное усилие N, действующее на тронк поршня, и усилие S, направленное по шатуну (рис. 3):

N = PΣ tg β  PΣ Sinβ = PΣ λ Sinφ,         Форм. 12

S = PΣ / Cosβ.         Форм. 13

Разложение суммарной движущей силы в КВМ
Рис. 3 Разложение суммарной движущей силы в кривошипно-шатунном механизме

В зависимости (Формула 12) принято допущение: tgβ ≅ Sinβ. При малых углах это допущение не дает большой погрешности. Разложив усилие S на составляющие, найдем касательную T и радиальную R силы, действующие в кривошипношатунном механизме:

T = S Sinφ + β = PΣSinφ + β / Cosβ;        Форм. 14

R = S Cosφ + β = PΣCosφ + β / Cosβ.        Форм. 15

Механическая напряженность кривошипа определяется совместным воздействием сил S и

 Pjц

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

Оценка механической напряженности в условиях эксплуатации

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

Максимальные напряжения в деталях движения и в подшипниках определяются максимальной движущей силой. Величина этой силы может быть принята равной:

PΣmax  Pz  Pjmax,

где:

  • Pz – усилие от максимального давления сгорания.

Можно написать:

Pjmax = Ms Rω2 1 + λ = Ms R πn/302 1 + λ = kj n2

здесь

kj = Ms R π/302 1 + λ.

Тогда:

PΣmax  Pz  kj n2.        Форм. 16

Работоспособность деталей зависит не только от максимальной движущей силы, но и от амплитуды ее изменения, определяющей усталостные разрушения. В 2-тактных двигателях изменение движущей силы может быть принято от 0 (в НМТ) до

 PΣmax

(в районе ВМТ). Поэтому амплитуда силы 2-тактных ДВС равна:

 

PΣ = 1/2 PΣmax  1/2 Pz  kj n2.       Форм. 17

В 4-тактных двигателях диапазон изменения силы — от (-Pjmax) в конце такта выталкивание газов до

 PΣmax

— в начале рабочего хода; амплитуда действующей силы определится равенством:

 

PΣ = 1/2 PΣmax + Pjmax  1/2 Pz.       Форм. 18

Следовательно, 4-тактные двигатели при прочих равных условиях имеют большую амплитуду ΔPΣ, чем 2-тактные, и являются более напряженными.

Максимальные касательные напряжения в шейках коленчатого вала зависят от максимального касательного усилия

 TΣmax

:

 

TΣmax = T + TΣ,

где:

  • TΣ

    – среднее касательное усилие;

  • ΔTΣ – амплитуда отклонения касательного усилия.

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

L = TΣ 2πR.              Форм. 19

С другой стороны, можно эту работу выразить через среднее индикаторное давление:

L = pi F 2R i/m,        Форм. 20

где

  • i – число цилиндров;
  • m – коэффициент тактности;
  • 2R = S – ход поршня;
  • F – площадь поршня.

Приравняв правые части зависимостей 19 и 20, можно найти:

TΣ = pi F 2R i/2πR m = Km pi,          Форм. 21

где:

  • Km = F i/(πm) – величина, постоянная для каждого двигателя.

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

TΣ  0,2 Pz,         Форм. 22

С учетом равенства 21 и 22 величина

 TΣmax

выразится зависимостью:

 

TΣmax  Km pi + 0,2 Pz.         Форм. 23

Из приведенных формул видно, что при неизменной частоте вращения величины

 PΣmax, PΣ, TΣ

полностью определится величиной максимального давления в цилиндре (давление сгорания)

Pz, а

 TΣmax

— величина Pz и Pi. Следовательно, в условиях эксплуатации для конкретного скоростного режима вполне обоснованным и наиболее объективным критерием механической напряженности двигателя является максимальное давление цикла Pz и среднее индикаторное давление Pi.

 

Расчетное определение сил, действующих в кривошипно-шатунном механизме

Расчетное определение сил, действующих в кривошипно-шатунном механизме, и последующая расчетная оценка механической напряженности может быть выполнена с помощью совмещенных би-центровой диаграммы проф. Брикса, индикаторной диаграммы и диаграммы сил инерции Толле (рис. 4). Для упрощения расчетов все силы относятся к 1 см2 площади поршня, т. е. имеют размерность давления.

Силы инерции Толле
Рис. 4 Совмещение диаграммы сил инерции Толле, индикаторной диаграммы и би-центровой диаграммы профессора Брикса при определении суммарной движущей силы

Индикаторная диаграмма берется по данным расчета индикаторного процесса или же принимается по данным индицирования реального двигателя. Диаграмма Брикса обычно строится под индикаторной диаграммой, как это показано на рисунке. Поправка Брикса ОО1 на конечную длину шатуна рассчитывается по формуле 9.

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

Pjвмт = Ms/F Rω2 1 + λ = ms Rω2 1 + λ102, МПа.        Форм. 24

Соответственно в НМТ:

Pjнмт = ms Rω2 1  λ102, МПа.        Форм. 25

Здесь:

  • ms = Ms/F – масса поступательно движущихся частей, отнесенная к 1 см2 площади поршня, кг/см2;
  • R – радиус кривошипа, м;
  • ω – угловая скорость, l/с;
  • λ = R/L – отношение радиуса кривошипа к длине шатуна.

Как было отмечено, величина MS, кг, складывается из массы поршня и поступательно движущейся части шатуна: MS = MП + MШS. Для определения величины поступательно движущейся массы шатуна Mшs делается допущение, что масса шатуна сосредоточена в 2-х точках: в центре поршневого пальца – поступательная масса MШS, в центре кривошипа – вращающаяся масса MШR (рис. 5).

Схематическая конструкция шатуна
Рис. 5 Схема шатуна

Величины MШS и MШR находятся из равенств:

MШ = MШS + MШR;

MШR/LS = MШS + LR,        Форм. 26

где:

  • LR и LS – расстояние от центра тяжести шатуна соответственно до центра кривошипа и центра поршневого пальца.

Обычно LR = 0,25 ÷ 0,52 L (верхние значения – для малооборотных ДВС).

Отложив ординаты

 Pjвмт и Pjнмт

на диаграмме Толле и найдя точки A и B, соединяют их прямой и находят точку С. Из точки С восстанавливают перпендикуляр, на котором откладывается величина CD = 3ms2λ. Найденная точка D соединяется с точками A и B, отрезки AD и DB делятся на равное количество частей. Найденные точки 1, 2, 3 и т. д. соединяются прямыми, как это показано на рис. 4. К этим прямым проводится касательная кривая, которая и определяет величину сил инерции в функции от хода поршня (или объема цилиндра).

 
Для нахождения давления газов и сил инерции в функции угла поворота коленчатого вала диаграмма Брикса разбивается из точки O1 лучами через 5 ÷ 15° пкв. Для каждого угла поворота коленчатого вала давление газов и сила инерции находится путем восстановления ординаты из точки пересечения соответствующего луча с полуокружностью радиуса R, проведенной из точки О.

При определении давления газов необходимо учитывать лишь избыточное давление, без учета давления в подпоршневой области (Po или Ps). Силы трения Pтр при расчете PΣ обычно не учитывается – в дальнейшем они учитываются механическим КПД. Силы веса поступательно движущихся частей, отнесенные к 1 см2 площади поршня, откладываются вниз от оси абсцисс на диаграмме Толле. В дальнейшем ординаты Pj + Pg снимаются от линии, проходящей параллельно оси абсцисс на расстоянии Pg.

Pg = ms9.81 102, МПа.

Определение сил динамики ведется в табличной форме в пределах φ = 0 ÷ 360° пкв для 2-тактных ДВС и φ = 0 ÷ 720° пкв для 4-тактных ДВС (табл. 1)

Таблица 1. Расчет сил динамики в кривошипно-шатунном механизме
φ° пквPгPj + PgPΣλ SinφNCos φ + βCos βRSin φ + βCos βT1T2T3TiTΣ = ΣTi
123456789101112131415
0
5
10
360

 
Расчет сил PΣ, N, R, Т1 ведется по формулам 11-15. Тригонометрические функции рассчитываются или принимаются по справочникам в зависимости от величины λ. Касательные усилия T1 первого цилиндра (столбец 10) повторяются в последующих столбцах с учетом угла заклинки кривошипов и последовательности работы цилиндров. Так, если в 6-цилиндровом 2-тактном двигателе последовательность работы цилиндров 1-5-3-4-2-6, то для 5-го цилиндра необходимо повторить записи T1 со смещением по фазе на 60° пкв, для 3-го на 120° пкв и т. д.

В многоцилиндровом двигателе коленчатый вал воспринимает и передает суммарный крутящий момент, определяемый суммарным касательным усилием TΣ = ΣTi. Расчетом TΣ для каждого значения φ заканчивается построение таблицы сил динамики. По данным таблицы строятся кривые сил динамики. Примерный вид кривых динамики для 6-цилиндрового 2-тактного дизеля на рис. 6.

Динамика 2-тактного дизеля
Рис. 6 Кривые сил динамики 6-цилиндрового 2-тактного дизеля

Правильность расчета сил динамики проверяется путем определения среднего касательного усилия по всем цилиндрам и нахождения индикаторной мощностиИндикаторная и эффективная мощность двигателя дизеля. Среднее касательное усилие можно найти путем планиметрирования площади под кривой TΣ(φ) или же аналитическим путем – суммированием ординат под кривой TΣ(φ):

TΣ  ΣTx/k,         Форм. 27

где:

  • k = 360m/(Δφ i) – число ординат, приходящихся на один период кривой TΣ(φ);
  • m – коэффициент тактности двигателя;
  • i – число цилиндров;
  • Δφ – шаг по оси абсцисс, с которым определяются ординаты Tx.

По найденному значению

 TΣ

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

 

Mкр = TΣ R F,          Форм. 28

Ni = Mкр 2π n/0.06, кВт.        Форм. 29

Найденная индикаторная мощность сравнивается с мощностью, определенной по индикаторной диаграмме. Считается допустимым расхождение не более 3-4 %.

Рекомендуем к прочтению: Параметры топливоподачи и факторы, их определяющие

Для расчета центробежной силы, действующей в кривошипе, необходимо вращающиеся массы шатуна (MШR), мотылевой шейки (MМШ), щек (MЩ) и противовесов (MПР) привести к одному радиусу по формуле: M′ = Mr/R (из условия

 Pjц = M rω2 = M Rω2

, где M′ – приведенная масса, r – расстояние от центра тяжести массы до оси коленчатого вала, R – радиус приведения). Обычно массы приводятся к радиусу кривошипа R. В случае, когда 2 противовеса расположены противоположно колену, неуравновешенные приведенные массы равны:

 

MR = MШR + MMШ + 2MЩ  2MПР.          Форм. 30

Центробежные силы вращающихся масс определяются зависимостью 10:

Pjц = MR Rω2.

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

Оценка условий работы подшипников коленчатого вала

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

Относительная скорость скольжения в подшипниках судовых дизелей весьма умерена. Она обычно находится в пределах: Cn = 3-12 м/сек (нижние значения — для малооборотных дизелей). Поэтому основным фактором, определяющим надежность работы подшипников судовых ДВС, является величина и характер изменения давления на подшипники.

Давление на головной подшипник полностью определяется величиной и характером изменения суммарной движущей силы PΣ, определенной при расчетах сил динамики (табл. 1). Для оценки давлений на мотылевой и рамовый подшипники также используются данные табл. 1.

Мотылевый подшипник

Сила Pм, действующая на мотылевый подшипник, определяется геометрической суммой векторов (рис. 3):

Pм = S + Pjцш,        Форм. 31

где:

  • S

    – вектор силы, действующей по шатуну;

  • Pjцш

    – вектор центробежной силы вращательно движущихся масс шатуна (

    Pjцш = MШR Rω2

    ).

Для упрощения сложения этих векторов сила 

S

представляется в виде сумме составляющих:

S = R + T,

, найденных при расчете сил динамики (табл. 1). Суммирование векторов удобно производить с помощью вспомогательной векторной диаграммы, как это показано на рис. 7.

Диаграмма давлений в дизеле
Рис. 7 Векторная диаграмма давлений на мотылевую шейку

Откладывая в соответствующем масштабе величины Ri и Ti, получают точки концов векторов 

Si

для каждого расчетного угла φ. Значение этого угла записывается в конце вектора

S

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

 

Построенная диаграмма обладает тем свойством, что положительному направлению оси R соответствует сжимающее усилие в щеке (действующее от мотылевой к рамовой шейке), а положительному для оси T — усилие, действующее в сторону направления вращения. Центробежное усилие не изменяется по величине, постоянно действует в сторону, противоположную рамовой шейке, и вызывает растягивающее усилие в щеке. Следовательно, усилие

 Pjцш

можно учесть при расчете

 Pм

, сместив вниз по оси R центр построенной диаграммы на расстояние

 Pjцш

(рис. 7, точка О1. Таким образом вектор

 Pjцш

складывается с вектором R. Тогда сила Pм будет равна расстоянию от точки O1 до конца вектора

 S

.

 

С помощью полярной диаграммы можно найти истинное положение вектора

 Pм

для любого угла поворота кривошипа. Для этого необходимо оси координат вместе с вектором Pм развернуть так, чтобы ось R совпала с направлением кривошипа, повернутым на требуемый угол, и была направлена к рамовой шейке. Тогда новое направление вектора

 Pм

будет истинным.

 
Снимая значение Pм с полярной диаграммы, можно построить развернутую диаграмму Pм(φ) удельных давлений на подшипник (рис. 8).

Диаграмма давлений в дизеле
Рис. 8 Развернутая диаграмма давлений на мотылевую шейку

Эти давления относятся к единице площади поршня. Определив среднее и максимальное значения давлений (

Pмср, Pмmax

), можно рассчитать среднее

 Kмср

и максимальное

 Kмmax

давления на подшипник, отнесенные к единице его площади:

 

 Kмср = Pмср F/DМlМ;

 Kмmax = Pмmax F/DМlМ,         Форм. 32

где:

  • F – площадь поршня;
  • DМ – диаметр подшипника;
  • lМ – длина подшипника.

В судовых дизелях при баббитовой заливке подшипников

 Kмср

достигает величины

 Kмср = 60

бар, а отношение

 Kмmax/Kмср

не превышает значений:

 Kмmax/Kмср = 2÷3.

В современных дизелях с высокими параметрами рабочего процесса средние давления достигли велчины

 Kмср = 7080

бар, отношение

 Kмmax/Kмср = 3,5

, что потребовало более прочного материала для заливки подшипников.

 

Рамовый подшипник

Рамовая шейка нагружена усилиями от 2-х цилиндров, между которыми она расположена. Усилия от каждого цилиндра передается на 2 рамовых подшипника. Поэтому при расчете удельных давлений Pp на шейку, расположенную между i-тым и (i+1)-м цилиндрами, необходимо взять 1/2 геометрической суммы векторов:

Pp = S + Pjцi + S + Pjцi+1/2.        Форм. 33

Центробежное усилие

 Pjц

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

 

 Pjц = mR Rω2.        Форм. 34

Как и при расчете мотылевой шейки, определение равнодействующей Pp удобно выполнять через проекции составляющих сил на осях R и T вспомогательной векторной диаграммы. Если на эти оси откладывать величину сил отстающего цилиндра (примем, что это — i-тый цилиндр), то силы опережающего цилиндра (i+1) будут повернуты в сторону направления вращения на угол Ψ — угол заклинки кривошипа (рис. 9).

Схема действия сил на рамовую шейку
Рис. 9 Схема действия сил во вспомогательной векторной диаграмме при определении давлений на рамовую шейку

Спроектировав силы Ri+1/2 и Ti+1/2 на оси R и T, получим проекции равнодействующей:

Rp = Ri/2 + Ri+1/2 cosψ + Ti+1/2 sinψ;

Tp = Ti/2  Ri+1/2 sinψ + Ti+1/2 cosψ.        Форм. 35

Расчет этих проекций для всех углов φ рекомендуется выполнять в табличной форме (табл. 2):

Таблица 2. Расчет удельных давлений на рамовую шейку
φ° пквRi/2Ri+1 Cosψ2Ti+1 Sinψ2RpTi/2Ri+1 Sinψ2Ti+1 Cosψ2Tp
123456789
0
10
20
360

 

Давление на рамовую шейку
Рис. 10 Векторная диаграмма давления на рамовую шейку

Равнодействующая

 Pp

находится с помощью вспомогательной векторной диаграммы. Силы инерции i-того

 (Pjц i/2)

и (i+1)-ого

 (Pjц i+1/2)

цилиндров учитываются путем смещения центра диаграммы из точки О в точку О′, как это показано на рисунке. По построенной векторной диаграмме строится развернутая диаграмма давлений на рамовую шейку (рис. 11), определяются среднее

 Ppср

и максимальное

 Ppmax

удельные давления на подшипник, отнесенные к единице площади поршня, после чего рассчитывается среднее и максимальное давление на единицу площади подшипника:

 

Kpср = Ppср F/Dplp;

Kpmax = Ppmax F/Dplp,        Форм. 36

где:

  • Dp – диаметр рамового подшипника;
  • lp – длина подшипника.
Давление на рамовую шейку
Рис. 11 Развернутая диаграмма давления на рамовую шейку

Для рамовых подшипников судовых ДВС считается допустимым:

 Kpср = 50

бар, отношение

 Kpmax/Kpср = 3,0.

 

Неравномерность вращения коленчатого вала

Периодичность работы цилиндров приводит к неравномерности воздействия на коленчатый вал суммарного касательного усилия и крутящего момента, что проявляется в неравномерности вращения коленчатого вала. Характер изменения суммарного касательного усилия и угловой скорости вращения дан на рис. 12.

Усилия и угловая скорость вращения коленвала
Рис. 12 Характер изменения суммарного касательного усилия TΣ и угловой скорости ω вращения коленчатого вала 6-цилиндрового дизеля

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

Явления, связанные с периодичностью работы цилиндров, могут быть оценены численно показателем — степенью неравномерности вращения коленчатого вала δ:

δ = ωmax  ωmin/ω,         Форм. 37

где:

  • ω = ωmax + ωmin/2, l/с

    – средняя угловая скорость вращения коленчатого вала.

Рассмотрим, какие факторы определяют величину δ. Для этого на основании теоремы об изменении кинетической энергии системы запишем:

A = I ω2max  ω2min/2,         Форм. 38

где:

  • ΔA – избыток (или недостаток) работы движущих сил по сравнению с работой сил со-противления;
  • I

    – момент инерции вращающихся масс валопровода.

Величина ΔA в соответствующем масштабе определяется площадью ΔF1 = ΔF2 на рис. 12. Преобразуем равенство 38. Запишем:

2A/I = ωmax + ωmin ωmax  ωmin = 2 ω ωmax  ωmin.

Тогда:

δ = ωmax  ωmin/ω = A/I ω2,          Форм. 39

Так как:

ω = πn/30;

I = mr2 = G/g r2,

То:

δ = A g/Gr2 30/πn2  900 A/Gr2n2.         Форм. 40

В этих равенствах:

  • n – частота вращения, об/мин;
  • G – приведенный вес вращающихся частей, кг;
  • r – радиус инерции вращающихся частей, м;
  • g – ускорение свободного падения, м/сек2;

Если взять не радиус, а диаметр инерции вращающихся масс: D = 2r, то равенство 40 можно записать в виде:

δ  3 600 A/GD2n2.         Форм. 41

Величина GD2 – это маховый момент вращающихся масс; он связан с момент инерции зависимостью:

GD2 = M g 2r2 = Mr2 4g = 39,2 I кг м2.        Форм. 42

Как видно из равенств  формул 39 и 41, неравномерность вращения коленчатого вала прямопропорциональна величине избыточной (или недостающей) работы ΔA и обратно пропорциональна моменту инерции вращающихся деталей

 I

(или маховому моменту GD2) и квадрату частоты вращения n2.

 

Избыточная работа может быть найдена с помощью графика TΣ(φ), построенного по итогам расчета сил динамики в кривошипно-шатунном механизме (рис. 12). Определив площадь ΔF наибольшей площадки между кривой TΣ(φ) и линией

 TΣ

, можно рассчитать ΔA:

 

A = F F mf,          Форм. 43

где:

  • F – площадь поршня, см2;
  • mf = mp mφ – масштаб площади диаграммы, кгм/см2мм2;
  • mp – масштаб оси ординат, кг/см2мм;
  • mφ = π 2R/lт – масштаб оси абсцисс, м/мм;
  • R – радиус кривошипа, м;
  • lm – длина оси абсцисс, соответствующая 360° пкв, мм;
  • ΔF – избыточная площадь, мм2.
Момент инерции вращающихся масс определятся суммой моментов инерции всех элементов валопровода
– коленчатого вала
 Iкв
, навешанных механизмов
 Iмех
, маховика
 Iм
, редуктора
 Iр
, валопровода
 Iв
, гребного винта с присоединенной массой воды
 Iгр
:

 

I = Iкв + Iмех + Iр + Iв + Iгр.         Форм. 44

Моменты инерции каждого элемента определяются с помощью зависимости:

Ix = Мx r2,          Форм. 45

где:

  • Mx – масса вращающегося элемента, кг/сек2;
  • r – его радиус инерции, м.

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

r = 0,5d2,

где:

  • d – диаметр вала.

В практике расчетов степени неравномерности вращения коленчатого вала часто пользуются другим способом, отличным от рассмотренного выше, не связанным с необходимостью определения численной величины избыточной работы ΔA. Для этого зависимость 39 преобразуется следующим образом:

δ = A/Iπn/302 Ni/Ni.

Индикаторная мощность Ni может быть выражена через работу Ao, совершаемую двигателем за один оборот коленчатого вала:

Ni = Ao n/60 75 илс.

Тогда:

δ = A Ni/Iπn/30260 75/Ao n = 60 75 302/π2 Ni/In3 A/Ao

Или:

δ = 0,41 106 Ni/In3 A/Ao.

Отношение работы ΔA/Ao численно равно отношению избыточной площади ΔF к площади Fo под линией

 TΣ

(рис. 6) на протяжении 360° пкв независимо от тактности двигателя. Поэтому можно записать:

 

δ = 0,41 106 Ni/In3 F/Fo.         Форм. 46

Если мощность дана в Ni кВт, момент инерции

 I

кг м2, то формула имеет вид:

 

δ = 0б548 107 Ni/In3 F/Fo.         Форм. 47

Обычно степень неравномерности вращения коленчатого вала находится в пределах:

  • δ = 1/200 ÷ 1/300 – для дизель-генераторов переменного тока;
  • δ = 1/100 ÷ 1/150 – для дизель-генераторов постоянного тока;
  • δ = 1/30 ÷ 1/40 – для главных судовых дизелей (без учета вращающихся масс валопровода и винта).

В двигателях малой и средней размерности момент инерции определяется главным образом массой маховика. Поэтому требуемая степень неравномерности вращения коленчатого вала обеспечивается выбором соответствующего маховика. В судовых малооборотных двигателях маховик практически не влияет на неравномерность вращения коленчатого вала — он выполняет роль приводного зубчатого колеса валоповоротного механизма.

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

Крутильные и осевые колебания валопровода

Валопровод любого судна представляет собой колебательную систему, поскольку включает в себя 2 разнородных накопителя энергии: упругость вала (потенциальная энергия) и массу (кинетическая энергия). При изменении внешнего воздействия на валопровод потенциальная энергия “закрученного” вала превращается в кинетическую энергию колеблющихся масс и наоборот — возникают крутильные колебания валопровода. При этом в валопроводе развиваются дополнительные касательные напряжения, иногда приводящие к самой серьезной аварии — поломке коленчатого вала. Характерным признаком разрушения вала из-за крутильных колебаний является расположение плоскости разлома под углом ~45° к оси вала.

Величина касательных напряжений при крутильных колебаниях при прочих равных условиях зависит от 2-х факторов:

  1. От амплитуды изменения крутящего момента (или суммарного касательного усилия);
  2. От совпадения частоты изменения внешнего воздействия и собственной частоты колебаний валопровода. Дизельная силовая установка имеет значительные отклонения Mкр от его среднего значения и поэтому требует особого внимания с точки зрения крутильных колебаний (по сравнению с паротурбинной и газотурбинной силовыми установками, у которых Mкр = const).

По Правилам Регистра перед установкой дизеля на судно выполняется расчет валопроводаРемонт валопроводов и судовых гребных винтов на крутильные колебания — определяются собственная частота колебаний, величина дополнительных касательных напряжений на различных частотах вращения, устанавливаются зоны критических оборотов, опасные для работы дизеля. На головном судне расчеты проверяются экспериментально и уточняются. Материалы расчета на крутильные колебания выдаются на судно. Рассмотрим факторы, определяющие собственную частоту колебаний валопровода. Для этого в простейшем случае для одноузловой колебательной системы (рис. 13) запишем уравнение крутильных колебаний в виде:

Крутильная система дизеля
Рис. 13 Крутильная система с одной массой (а) и ее схематическое изображение (б)

I d2 φ/dτ2 + cφ = 0,          Форм. 48

где:

  • I

    – момент инерции массы M относительно оси вала;

  • φ – угол отклонения массы от равновесного состояния;
  • C – жесткость вала;
  • τ — время.

Жесткость вала определяется зависимостью:

C = CII Wp/l,         Форм. 49

где:

  • GII – модуль упругости материала вала II рода;
  • Wp – полярный момент инерции поперечного сечения вала (для круглового вала Wp = πd4/32);
  • l – длина вала.

Угол закрутки вала φ зависит от величины скручивающего момента Mкр и жесткости вала C:

φ = Mкр/C = Mкрl/GIIWp.

Угловая частота колебаний со связана с жесткостью C и моментом инерции

 I

известной зависимостью:

 

ω = C/I, l/сек.          Форм. 50

Подставив это равенство в формулу 48, получим:

d2φ/dτ2 + ω2φ = 0

Решением этого дифференциального уравнения является уравнение вида:

φ = B sin ωτ + A cos ωτ.         Форм. 51

Приняв начальные уравнения: τ = 0; dφ/dτ = 0; φ = A, получим решение: φ = A cos ωτ.

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

Из зависимостей 49 и 50 можно найти период колебаний T и собственную частоту колебаний ν:

T = 2π/ω = 2π I·l/GII·Wp, с;          Форм. 52

ν = 30/πω = 30/πGII·Wp/I·l, кол/мин.         Форм. 53

Из равенства 53 можно сделать вывод, что собственная частота колебаний снижается при увеличении момента инерции (увеличении массы), увеличении длины валопровода, уменьшении полярного момента (уменьшении диаметра вала) и не зависит от величины крутящего момента. Этот вывод справедлив и для любой сколь угодно сложной крутильной системы.

Схема крутильной 2-массовой системы дизеля
Рис. 14 Схема крутильной 2-массовой 1-узловой колебательной системы

В случае 2-массовой колебательной системы (рис. 14) угловая амплитуда в точке “y” равна 0. Это — узел колебаний. Как видно, узел один, поэтому система называется одноузловой. Частота собственных (свободных) колебаний такой системы определяется по формуле:

ω = I1 + I2·C12/I1·I2,         Форм. 54

где:

  • I1 и I2

    – моменты инерции масс M1 и M2;

  • C1-2 – жесткость вала.

При 3-х массах возможна одно- и двухузловая форма колебаний (рис. 15).

Схема судового валопровода
Рис. 15 Схема валопровода в виде 3-массовой 2-узловой колебательной системы

Частота собственных колебаний определяется зависимостью:

ω = A ± A2  B,         Форм. 55

где:

  • (+) при 2-узловой схеме;
  • (-) при 1-узловой схеме;

A = C12 I1 + I2/2I1 I2 + C23I2 + I3/2I2 I3,         Форм. 56

B = C12 C23 I1 + I2 + I3/I1 I2 I3.        Форм. 57

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

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

В этом случае массы всех цилиндров считаются сосредоточенными в центре тяжести двигателя, что позволяет рассчитать собственную частоту колебаний валопровода наиболее важных форм — одноузловой (по формуле 54) и 2-узловой (по формуле 55).

При расчетах системы валопровода на крутильные колебания гребной и коленчатые валы заменяются общим приведенным (эквивалентным) валом, имеющем по всей длине неизменный диаметр do и обладающий такими же упругими свойствами, что и действительный вал. По условию упругой эквивалентности, должно быть равенство углов скручивания при действии одного и того момента:

φ = Mкр l/GII Wp = Mкр lo/GII Wpo,         Форм. 58

где:

  • lо = l Wpo/Wp – приведенная длина вала;
  • Wpo = πd о4/32 – полярный момент инерции приведенного вала;
  • 2ro = do – диаметр приведенного вала.

Действительные массы заменяются приведенными массами по формуле (из условия

 I = mr2 = moro2

):

 

mo = m r2/ro2.         Форм. 59

Расчет обычно выполняется в относительной форме — приведенная масса одного цилиндра и приведенная длина между 2-мя цилиндрами принимаются равными 1; элементы остальных участков выражаются в долях от принятого за единицу.

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

nрез = νсобств/K,         Форм. 60

где:

  • K = 1, 2, 3, 4… – последовательный ряд целых чисел.

При совпадении частоты вращения коленчатого вала с найденными резонансными частотами в валопроводе развиваются крутильные колебания, увеличивается амплитуда колебаний, появляются дополнительные скручивающие моменты. Однако не все резонансные частоты являются опасными для работы двигателя. Силы трения в колебательной системе (силы “упругого гистерезиса”) препятствуют развитию крутильных колебаний.

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

Наиболее опасными являются частоты вращения, определяемые зависимостью:

n = νсобств/i m,         Форм. 61

где:

  • i – число цилиндров;
  • m – коэффициент тактности двигателя.

При такой частоте вращения частота вспышек работающих цилиндров совпадает с собственной частотой колебаний валопровода.

К примеру, на судах типа “Лисичанск” собственная частота колебаний валопровода в 1-узловой форме равна νсобств= 387 кол/мин. Резонансными являются частоты 387; 193,5; 129; 96,75; 77,4; 64,5; 55,3; 48,4; 43; 38,7;…. Поскольку главный двигатель является 9-цилиндровым 2-тактным, то наиболее опасной является частота вращения n = 387/9 = 43 об/мин.

Расчет на прочность позволил установить, что дополнительные касательные напряжения при этом — более 350 кг/см2, что выше допустимого (рис. 16). Фирма-строитель установила зону запретных частот вращения nкр = 40-47 об/мин. На остальных резонансных частотах в зоне работы двигателя дополнительные касательные напряжения не выходят за пределы 100 кг/см2; эти частоты вращения не являются опасными.

Напряжения в валопроводе т/х “Лисичанск”
Рис. 16 Дополнительные напряжения в валопроводе т/х “Лисичанск” при работе на резонансных частотах

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

  1. Зона критических оборотов смещается путем: а) изменения жесткости валопровода (варьированием длины и диаметра вала); б) изменением массы маховика или противовесов на кривошипах; в) разделением крутильной системы с помощью эластичных муфт;
  2. Амплитуда крутильных колебаний уменьшается путем установки успокоителей — антивибраторов или демпферов. В антивибраторе создаются колебания, противоположные крутильным, что способствует уменьшению амплитуды крутильных колебаний. В демпферах энергия крутильных колебаний срабатывается главным образом за счет трения в упругих элементах.

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

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

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

Схема демпфера
Рис. 17 Схема гидравлического поршневого демпфера осевых колебаний

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

  • вместо момента инерции

     I

    — массу M;

  • вместо крутильной жесткости вала C — осевую жесткость Co;
  • вместо угла закрутки φ — линейную дифференциацию Δl.
Сноски
Sea-Man

Нашли опечатку? Выделите и нажмите CTRL+Enter

Август, 30, 2016 2656 0
Добавить комментарий

Читайте также

Текст скопирован
Пометки
Избранные статьи
Loading

Здесь будут храниться статьи, сохраненные вами в "Избранном". Статьи сохраняются в cookie, поэтому не удаляйте их.

Статья добавлена в избранное! Перезагрузка...