CC BY
79
Вестник МГТУ, том 16, №1, 2013 г.
стр. 193-196
УДК 656.618.1 : [629.5.018.712 : 517.958]
Расчет усилий, действующих на объекты буксировки со стороны буксирной связи
Ю.И. Юдин, С.В. Пашенцев, В.В. Каян
Морская академия МГТУ, кафедра судовождения
Аннотация. В статье представлена методика расчета усилий, действующих на объекты буксировки со стороны буксирной связи, используемая в моделировании процесса буксировки одного судна другим. Показано, что применение приведённого алгоритма учета натяжения в гибком растяжимом буксирном тросе под действием собственного веса целесообразно при моделировании движения буксирной системы "буксирующее судно - буксирный трос - буксируемое судно".
Abstract. The paper presents the method of calculation of forces acting on objects engaged in towing operations from towing connection used in modeling of towage process of one vessel by another. It has been shown that the application of the algorithm is expedient for modeling of movement of the towage system "towing vessel -towing line - towed vessel".
Ключевые слова: буксирная система, буксировочный трос, буксируемое судно, буксирующее судно Key words: towage system, towing line, towed vessel, towing vessel
При исследовании процесса буксировки приходится рассматривать три объекта этого процесса -буксировщик, буксируемое судно и буксирный трос, которые в представленной совокупности называются буксирной системой. Математические модели судов, участвующих в процессе буксировки, представляют собой системы дифференциальных уравнений, в левые части которых входят усилия, действующие со стороны буксирного троса. Рассмотрим буксирный трос как объект буксирной системы с целью определения его натяжения при различных способах крепления на судах, участвующих в буксировке.
2. Методика расчета
Рассмотрим геометрию троса в вертикальной плоскости XOZ, которая показана на рис. Предварительно будем считать, что буксирный трос полностью находится в воздухе, а его начальная точка КВ - в кормовой части буксирующего судна, при этом конечная точка буксирного троса NO крепится в носовой части буксируемого судна.
1. Введение
Z
KB
Тг
Z
Рис. Геометрия троса в вертикальной плоскости XOZ
193
Юдин Ю.И. и др. Расчет усилий, действующих на объекты буксировки...
Введём местную систему координат, ось Z которой проходит через точку КВ, при этом начало координат на этой оси не имеет значения для последующих расчётов. Пусть, для определённости, это будет уровень дна или уровень киля буксирующего судна. Ось Х направим в сторону носа буксируемого судна. Придадим начальной и конечной точкам троса координаты KB(xkb, zkb) и NO(xno, zno). Как известно, тяжелый трос, принятый традиционно за гибкую нить, принимает форму цепной линии. Эту форму можно описать одним уравнением цепной линии общего вида, в которое входят три параметра - а, с1 и с2:
z = а х ch (X——) - c2. (1)
a
Уравнение записано во введённой местной системе координат, связанной с точками крепления троса, и при разных комбинациях значений параметров с1, с2 и а описывает все ситуации его расположения. Наша задача состоит в получении зависимости натяжения троса от положения точек КВ и NO.
Параметры с1, с2, а можно определить из двух граничных условий z(xkb) = zkb, z(xno) = zno, приводящих к двум алгебраическим уравнениям:
а х ch(-
а х ch(-
xkb - cl
) - c2 = zkb:
-cl
) - c2 = zn
и одному интегральному для длины l цепной линии:
— - c
l = а х sh(---1)
а
= а х [sh(
-no % - sh(———)].
(2)
(3)
а а
Вычитая одно уравнение системы (2) из другого в любом порядке, устраним константу с2, т.к. она не влияет на форму цепной линии. Вот почему нам не важно, где начало координат на оси Z - в разностном уравнении участвует лишь разность координат z. Теперь три уравнения можно свести к двум уравнениям с неизвестными с1 и а:
а
no
а
а х [ch(-kb Cl ) - ch(Xno Cl )] = za - Zno;
а а (4)
а х [sh(-no - ci ) - sh(-kb - ci )] = l. а а
Это система двух трансцендентных уравнений с гиперболическими функциями, которые возможно решить только численно. Для этого использована вычислительная среда MathCad, которая имеет инструментарий для подобных решений в виде оператора rootfx), x) с указанием функции, корень которой необходимо найти. Приведём решение частной задачи, когда высота крепления буксировочного троса на буксире и объекте совпадают, т.е. разность zkb - zno = 0. Такая частность не будет сказываться на принципиальных выводах из решения задачи. Тогда первое уравнение из (4) приводит к равенству
ch (Xkb-cL) = ch (-no-cL) , а а
что позволяет получить следующее решение:
( Xkb - ci ) = -( Xno - ci )
а а
Это равенство, в свою очередь, позволяет определить значение константы с1 = (xno + xkb)/2. Подстановка найденного значения константы с1 во второе уравнение системы (4) даёт:
sh(———) = l /2а . (5)
2а
Данное трансцендентное уравнение позволяет рассчитать величину параметра а, определяющего горизонтальное натяжение троса. Преобразуем уравнение (5) к виду, наиболее удобному для численного решения:
sh(w) = (l/Ax) х w = ух w,
w = (xno - xkb)/2a = Ax/2a. (6)
Уравнение сведено к простейшей форме sh(w) = у х w, которая допускает простое решение, учитывая, что параметр у = l/Ax > 1 определяет степень провисания троса, т.к. длина буксирного троса l всегда больше или равна горизонтальному расстоянию Ax между точками его крепления на буксирующем и буксируемом судах.
194
Вестник МГТУ, том 16, №1, 2013 г.
стр. 193-196
Численное решение уравнения (6) произведено с помощью оператора MathCad root(sinh(w) - ух w, w), для установленного ряда значений у получены корни, представленные в табл.
Таблица
у 1.01 1.02 1.03 1.04 1.05 1.06 1.07 1.08 1.09 1.10
w 0.271 0.366 0.438 0.488 0.550 0.594 0.642 0.684 0.725 0.765
Тг 6541 4730 3914 3479 3058 2804 2571 2390 2234 2098
Az 30 41 49 54 60 65 70 74 78 82
lw 396 411 418 422 425 427 429 431 433 434
Tr 2.79 2.89 2.94 2.95 2.98 2.99 3.00 3.01 3.02 3.02
Чтобы использовать эти результаты в программных расчётах, нами получена аналитическая аппроксимация значений корней в зависимости от параметра у В указанной зависимости используется полином пятой степени с очень высокой степенью адекватности:
w(y) = (-0.5105 + 2.389у- 4.471/3 + 4.182 у3 - 1.956у4 + 3.657у3)10б. (7)
Горизонтальная составляющая натяжения Тг буксирного троса постоянна по всей его длине и зависит от веса его погонного метра q и константы а, которая вычисляется в ходе решения системы (4). Определение Тг вполне традиционно с помощью уравнения равновесия цепной линии. Оно даёт в ньютонах:
Тг = а х q х g , (8)
где g - ускорение свободного падения, а известный параметр а представлен как размерная величина, имеющая размерность длины.
Вертикальная составляющая может быть найдена из простого равенства Тв = Тг tg(a) (см. рис.), но она не входит в уравнения плоского движения объектов буксировки. При этом tg(a) находится как производная z’x = sh((x - c1)/a). Натяжение буксирного троса можно определить, пользуясь очевидным соотношением между составляющими натяжения
T = д/T2 + T2 = Тг41 + tg2 а =Т 1 + sh2(—-Cl) = Тгch(—-^) = T(z + c2).
V a a a
(.2, x-C_)
a a a
Горизонтальная составляющая натяжения буксирного троса зависит от его длины и от диаметра, определяющего, в свою очередь, величину веса погонного метра буксирного троса q.
В качестве примера расчёта горизонтальной составляющей произведён её расчёт для стального троса ТК-133 ГОСТ 3067-55 диаметром 52.5 мм и длиной 450 м. Именно такой трос выбран для буксируемого судна. Вес погонного метра троса составляет 11.387 кг. Результаты такого расчёта приведены в третьей строке таблице в кН. Эти результаты полезно сопоставить с разрывной прочностью выбранного буксирного троса в 1280, 1600 и 1814 кН при предельной прочности нитей в 120, 150 и 170 кг/мм2 соответственно.
Очень важным фактором при моделировании является провес буксирного троса. Если провес невелик, то трос не касается воды, и его дополнительное сопротивление при движении не учитывается. Если же провес превышает высоту подвеса концов троса, то некоторая часть троса оказывается в воде, и следует учесть её гидродинамическое сопротивление при движении судов. Для решения этих вопросов рассчитаем величины провеса троса. Для этого найдем величину разности координаты z для середины цепной линии, где x = (хкЬ + хпо)/2, и точки подвеса троса, где х = 0 или х = Ах. Тогда получим с учётом (6):
Az = l/2wyx [ch(w/2) - 1)].
(9)
Используя найденные выше значения w при соответствующих величинах у, находим провесы Az. Они приведены в той же таблице в четвертой строке для длины троса l = 450 м. Разумеется, это чисто теоретические расчёты; практически таких провесов при реальной буксировке никогда не допускается.
Трос, попадая в воду, подвергается действию ряда сил, которые уменьшают провес. Это статическая выталкивающая сила, которая действует вертикально, и гидродинамическая подъёмная сила, также имеющая вертикальную составляющую, кроме горизонтальной гидродинамической составляющей. Поэтому, наряду с расчётом провеса, мы рассчитываем длину части троса, находящуюся в воде.
Как и выше (3), длину участка цепной линии находим по интегральной формуле
, ~ w — - c,.
lw = 2a x sh(------1)
= 2a x [sh() - sh(——Cl)] = -2a x shi^—^L),
a
a
a
a
(10)
где хw - координата Х точки пересечения линии троса с уровнем воды (рис.). Для той же точки из уравнения цепной линии имеем:
195
Юдин Ю.И. и др. Расчет усилий, действующих на объекты буксировки...
а х [ch (■
XkjL^~) - ch)] = Zkb - Zw
откуда легко находим
ch (
X C) = ch(X^1)-(Zkb -Zw)/a = ch(Xk^1)-hw/a . a a a
Выражаем гиперболический синус согласно известной формуле для гиперболических функций sh2(x) = ch2(x) - 1 и подставляем в выражение для длины lw (10). Получаем длину в виде:
a
a
lw = 2a. (ch(———) - hw / a)1 -1 . (11)
V a
Для примера рассчитаны длины участка троса в воде при высоте крепления троса над водой hw = 7 м. Эти результаты приведены в таблице пятой строкой. Если использовать для расчёта той же величины формулу, приведённую в работе (Войткунский, 1985), то получим подобные величины, отличающиеся примерно на 2 % в сторону больших значений. Однако сама формула проще и не содержит специальных функций: lw = 2 х (l2 - 8ahw)0 5. Если пойти по пути использования этой формулы, то логично применить формулу справочника (Щетинина, 1975) и для расчёта дополнительного сопротивления троса диаметром d, движущегося в воде со скоростью и. Она такова (кН):
Rtr = 4 х lw х d х и2 х g/1000. (12)
Рассчитанные по этой формуле значения дополнительного сопротивления приведены в таблице в шестой строке для d = 0.045 м и скорости буксировки и = 2 м/с. Хорошо видно, что эти значения незначительны в сравнении с основной силой натяжения троса (строка 3).
Более сложным является вопрос о форме, которую принимает тяжелый трос, движущийся в воде. Она отличается от цепной линии, но при моделировании буксировки форма буксирного троса в воде не учитывается и потому здесь не рассматривается.
Наконец, остается ещё один режим при буксировке, когда расстояние между точками крепления буксировочного троса становится больше, чем первоначальная длина троса. Это означает, что трос получает некоторое упругое удлинение hi, что вызвано силой его натяжения
Ttr = Егх (nd2/4) х (hl/l), (13),
где Etr - модуль упругости троса. Это не просто модуль упругости стали - материала троса, а модуль, который учитывает определённую многопрядную структуру троса, изменяющую величину модуля.
Приведённые выше расчёты позволяют сформулировать следующий алгоритм учёта действия буксировочного троса:
- рассчитывается расстояние hx между точками крепления троса;
- hx сравнивается с первоначальной длиной троса l;
- если hx > 1 (трос не провисает), то сила натяжения троса рассчитывается по формуле (13);
- если hx < 1 (трос провисает), то рассчитывается провес троса hz;
- если hz < hw (трос полностью в воздухе), то рассчитывается натяжение троса по формуле (8);
- если hz > hw (трос частично в воде), то рассчитывается натяжение троса по формуле (8), с умножением результата на множитель 0.875 за счёт его облегчения в воде;
- рассчитывается дополнительное усилие за счёт гидродинамического воздействия на трос по формуле (12) и складывается с основным натяжением;
- полученный результат учитывается при моделировании в уравнениях продольного равновесия участвующих в буксировке судов.
3. Заключение
Предложенный авторами алгоритм учета натяжения в гибком растяжимом буксирном тросе под действием собственного веса позволяет использовать его при моделировании движения буксирной системы "буксирующее судно - буксирный трос - буксируемое судно". Это дает возможность отказаться от простейшей схемы буксировки, когда трос рассматривается как жесткий стержень, и получить тем самым более объективную картину буксирного процесса.
Литература
Войткунский Я.И. Справочник по теории корабля. В 3 т. Под ред. Я.И. Войткунского. Л., Судостроение, т. 3, с.543, 1985.
Щетинина А.И. Управление судном и его техническая эксплуатация. Под ред. А.И. Щетининой. М., Транспорт, 608 c., 1975.
196
