ВЕСТНИК

МОРСКОГО

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

УНИВЕРСИТЕТА

серия

СУДОВОЖДЕНИЕ

Выпуск

Владивосток 2009

УДК 656.61.052(066)

Вестник Морского государственного университета. Вып. 5. Серия: Судовождение. – Владивосток: Мор. гос. ун-т, 2009. – 149 с.

Редакционная коллегия

, д. н. н., проф.(отв. ред.),

, к. т.н., проф. (отв. ред.),

, д. т.н., проф.,

, к. ю.н., доцент.

ISBN -2

© Морской государственный университет имени адмирала ., 2009

СРАВНИТЕЛЬНЫЙ АНАЛИЗ РЕЗУЛЬТАТОВ МОДЕЛИРОВАНИЯ

ДВИЖЕНИЯ СУДОВ НА ЦИРКУЛЯЦИИ В УСЛОВИЯХ МЕЛКОВОДЬЯ

, А.,

МГУ им. адм. , г. Владивосток

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

В постоянно меняющихся эксплуатационных условиях (загрузка судна, его осадка, глубина под килем) точное маневрирование становится чрезвычайно за­труднительным. Судоводитель не имеет инструментария, с помощью которого он мог бы определить, например, на какой траектории окажется судно при по­вороте на новый курс при конкретном угле перекладки руля, тем более в усло­виях мелководья. Имеющаяся в настоящее время на судах информация о ма­невренных элементах, предусмотренная Резолюцией ИМО А.601(15), которая представлена, в большинстве случаев, лишь одной циркуляцией на глубокой воде при максимальном угле перекладки руля, в обычной судоводительской практике планирования и осуществления перехода является практически беспо­лезной и может быть использована лишь в экстренных ситуациях.

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

Моделирование движения судов на циркуляции выполнено на навигацион­ном тренажере NTPro 3000 производства фирмы Тransas Marine. Движение осуществлялось при отсутствии внешних воздействий (ветра, течения, волне­ния). Информация для построения циркуляций получалась следующим об­разом:

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

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

– выбранный для циркуляции модели судна район мелководья выделялся и ограничивался нужной глубиной;

– включались приборы, необходимые для получения данных о движении судна (GPS, эхолот, управление судном переводилось на «автомат»);

– на электронной карте выставлялся буй, вокруг которого планировалось движение судна на циркуляции;

– со станции инструктора запускалось упражнение. Модель судна начинала движение с заданным курсом и скоростью;

– на конец третьей минуты на глубокой воде, когда модель судна выходила на прямой курс, и на конец пятой минут, когда движение на прямом курсе ста­билизировалось, и происходила окончательная просадка модели судна, которая зависела от глубины под килем, на станции инструктора движение останавлива­лось. В момент остановки движения получался «стоп–кадр», т. е. модель судна останавливалась мгновенно, поэтому была возможность снять все необходимые данные, которые в момент остановки находились в «замороженном» состоянии;

– с дисплея снимался курс (с точностью до 0.1о), линейная скорость (с точностью до 0.1 узла), глубина под килем (с точностью до 0.1 м), скорость по­ворота (с точностью до 0.1 град/мин.), управление судном снималось с автома­тического, а руль перекладывался на нужный угол;

– с электронной карты снимались пеленг и расстояние до выставленного буя (пеленг с точностью до 0.1о, расстояние - до 0.1 м);

– производится пуск движения модели судна. В момент пуска модель про­должает движение с теми же параметрами курса и скорости, которые были за­фиксированы в момент остановки. В этот же момент начинается перекладка руля на заданный угол. Модель судна начинает движение по циркуляции с за­данным углом перекладки;

– в дальнейшем модели судна останавливались через каждые 10о измене­ния курса. Снимались те же данные, которые указаны выше, за исключением перекладки руля.

Таким образом, данные для построения циркуляций снимались с достаточ­ной точностью, чтобы можно было бы сделать определенные выводы. Модели­рование движения на циркуляции выполнялось для пяти различных судов с обычными обводами корпуса, основные характеристики которых приведены в таблице 1.

С применением стандартного пакета Excel был выполнен статистический регрессионный анализ всех данных моделирования. Процедура и некоторые результаты этого анализа представлены в работе [1]. В частности, установлено, что зависимость радиуса установившейся циркуляции на мелководье RMуст от отношения глубины к осадке судна для всех судов имеет степенной характер

RMуст = b (Н/d , (1)

где b – константа регрессии, с – коэффициент регрессии.

Таблица 1

Основные характеристики судов, использованных при моделировании

Тип

модели судна

Водоиз-

мещение, т

Длина

наибольшая, м

Ширина,

м

Осадка, м

носом

кормой

Балкер

33089

182.7

22.6

10.1

10.7

Танкер

77100

242.8

32.2

12.5

12.5

Ро-Ро

25400

184.2

30.6

8.2

8.2

Паром

24841

230.9

29.2

8.0

8.0

Контейнеровоз

132540

347.0

42.8

14.0

14.0

В таблице 2 сведены результаты статистического регрессионного анализа зависимости радиуса установившейся циркуляции на мелководье RMуст от отношения глубины к осадке Н/Т. Во всех случаях коэффициент парной корреляции превышает 0,94. Такие же результаты получены для выдвига и тактического диаметра циркуляции.

Таблица 2

Расчетные параметры корреляции между радиусом установившейся циркуляции на мелководье RMуст и отношением глубины к осадке Н/Т для различных судов

Угол перекладки руля (δ)

10°

15°

20°

25°

30°

35°

Константа

регрессии b

Балкер

1211

1087,2

989,51

925,5

892,9

850,45

Контейнеровоз

13885

8238,1

5543,9

4072,1

3216

2837,9

Танкер

2156,4

1748,8

1534,4

1350,1

1169,9

1055,1

Паром

1881,7

1467,8

1209

1053

903,29

818,51

Ро-Ро

5984,2

3981,8

2984,6

2363,4

1962,4

1611,5

Коэффициент регрессии с

Балкер

-0,4216

-0,5115

-0,5547

-0,6293

-0,6997

-0,7591

Контейнеровоз

-1,7724

-1,5432

-1,3904

-1,2881

-1,2325

-1,2355

Танкер

-0,4473

-0,4491

-0,4681

-0,4966

-0,4436

-0,4501

Паром

-0,466

-0,487

-0,506

-0,51

-0,52

-0,5541

Ро-ро

-1,0172

-0,9501

-0,8769

-0,8017

-0,766

-0,7137

Коэффициент парной

корреляции

Балкер

-0,9925

-0,9909

-0,9917

-0,989

-0,9895

-0,9933

Контейнеровоз

-0,9939

-0,9871

-0,9828

-0,9642

-0,9446

-0,9812

Танкер

-0,9658

-0,9729

-0,9752

-0,9642

-0,978

-0,9944

Паром

-0,9553

-0,9462

-0,9533

-0,9642

-0,9445

-0,9649

Ро-ро

-0,9819

-0,971

-0,9687

-0,9733

-0,9446

-0,967

Коэффициент

детерминации R2

Балкер

0,9852

0,9818

0,9835

0,9781

0,9792

0,9866

Контейнеровоз

0,9879

0,9743

0,9704

0,9659

0,9731

0,9627

Танкер

0,9328

0,9466

0,9511

0,9783

0,9742

0,9888

Паром

0,9126

0,8953

0,9089

0,9296

0,8922

0,931

Ро-ро

0,9357

0,9428

0,9383

0,9473

0,9527

0,9389

Для определения параметров циркуляции для любого промежуточного зна­чения угла перекладки руля δ были получены зависимости свободного члена b и коэффициента с регрессии (1), которые имеют, соответственно, логарифмиче­ский и линейный характер.

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

RMуст =bm1(H/d)m2δ =1807,4∙0,8(H/d)0,98δ. (2)

где b, m1 и m2 – константа и коэффициенты регрессии, соответственно.


На рис. 1 показан трехмерный график функции (2), полученный с использованием пакета MatLab 7. Однако для практических целей более удобно использовать двухмерное представление этой функции в виде номограммы, представленной на рис. 2.

Достоверность полученных закономерностей может быть проверена путем сравнения результатов моделирования с фактическими параметрами циркуля­ции реальных судов. Однако такой информации достаточно мало. Так, в работе [2] представлены данные натурных экспериментов по определению параметров циркуляции т/х «Николай Подвойский», имеющего следующие основные размерения: длина наибольшая – 186,2 м; ширина наибольшая – 23,45 м; водоизмещение в полном грузу – 30810 т; осадка в полном грузу – 9,86 м, диаметр винта – 5,9 м. В графическом виде эти данные показаны на рис. 3. вместе с результатами обработки данных циркуляции балкера, основные размерения которого почти совпадают с размерениями т/х «Николай Подвойский». Из рисунка видно, что представленные зависимости практиче­ски идентичны.

Десятипроцентное увеличение радиуса установившейся циркуля­ции т/х «Николай Подвойский» объясняется тем, что его длина на 3 метра больше по сравнению с балкером, и, главное, тем, что его осадка почти на метр превышает осадку балкера.

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

модель Рейнольдса [3] (3)

и модель Гофмана [4] (4)

Однако обеим этим моделям присущ один и тот же недостаток: они не учитывают такие важные параметры, как угол перекладки руля и размеры судна, его осадку. Эти модели устанавливают зависимость только между радиусом установившейся циркуляции и отношением глубины к осадке судна. Хотя один и тот же радиус циркуляции может быть и у судна длиной 50 м, и у судна длиной 250 м, естественно, при разных углах перекладки руля. Тем не менее, с использованием обеих этих моделей были рассчитаны значения радиусов установившейся циркуляции при различных значениях отношения глубины к осадке применительно к двум ситуациям моделирования циркуляции балкера: при угле перекладки 10о, когда радиус установившейся циркуляции был равен 640 м, и при угле перекладки 30о, когда радиус уменьшился до 315 м. Результаты расчетов для двух значений углов перекладки руля в графическом виде представлены на рис. 4.


Расчеты показали, что результаты использования трех методов (моделиро­вание, модель Рейнольдса и модель Гофмана) приближаются друг к другу при уменьшении угла перекладки руля. При этом результаты моделирования зани­мают примерно среднее значение между расчетами по моделям Рейнольдса и Гофмана, когда отношение глубины к осадке становится менее 2 (рис. 4, б). Однако при больших углах перекладки руля в математической модели, исполь­зованной при моделировании, влияние мелководья проявляется гораздо более выражено, чем в обеих моделях, особенно в модели Рейнольдса. По сравне­нию с расчетами по модели Гофмана результаты моделирования отличаются примерно на 10-15%.

Окончательные выводы о степени достоверности любого из рассмотрен­ных методов могут дать лишь обширные натурные эксперименты, однако их проведение практически крайне маловероятно, хотя бы из-за проблем, связан­ных с выборов районов моря с ровным дном и соответствующими глубинами. Тем не менее, учитывая указанный выше основной недостаток моделей Рейнольдса и Гофмана, имеются все основания полагать, что зависимости, по­лученные на основе анализа результатов моделирования, являются вполне рабо­тоспособными и могут быть использованы в судоводительской практике при планировании и осуществлении плавания в стесненных водах.

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

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

2) на мелководье радиус циркуляции заметно увеличивается, причем чем больше угол перекладки руля, тем больше это увеличение. В частности, для балкера при отношении глубины к осадке, равном 1,22, радиус циркуляции увеличивался в 2,8 раза при угле перекладки руля 35° и в 1,8 раза при угле перекладки руля 10°;

3) увеличение радиуса циркуляции вследствие влияния мелководья необходимо учитывать при предварительной и исполнительной прокладке при плавании в стесненных районах.

Литература

1. , Лентарёв анализ результатов моде­лирования движения судна на циркуляции. Вестник Морского государственного университета. Вып. 23 Серия: Судовождение. - Владвосток: Мор. гос. ун-т, 2008. - С. 58-65

2. Иванов особенности поворотливости больших танкеров. ЦБНТИ ММФ. Экспресс-информация. Серия «Безопасность мореплавания». М., 1978. - № 4(104). - С. 2-11.

3. Rejnolds J. Ship’s turning characteristics in different water depths. Safety at Sea International. 1976. – № 9

4. Гофман и расчет поворотливости судов внутреннего плавания. Л: Судостроение, 1971.

ПОГРЕШНОСТИ РАСЧЕТА ТОЧЕК ПОВОРОТА ПРИ ПЛАВАНИИ НА МЕЛКОВОДЬЕ

, А.,

МГУ им. адм. , г. Владивосток

Современные условия судоходства, среди прочего, характеризуются тем, что все большая часть ходового времени, особенно для крупных судов, прихо­дится на плавание в стесненных районах. В этих условиях особое значение приобретает проблема надлежащего планирования рейса, что требует уже на стадии предварительной прокладки учета всех возможных факторов, которые влияют на точность плавания, и решения ряда задач, обеспечивающих наиболее точное движение судна по выбранному маршруту. К числу таких факторов и задач относятся влияние мелководья и расчет точек поворота. Помимо того, что такие задачи рекомендуются хорошей морской проктикой, их решение является обязательным условием при выполнении предварительной прокладки в соответствии с требованиями Резолюции ИМО А.893(21) «Руководство по планированию рейса». Эти требования логически вытекают из положений Конвенции ПДМНВ-78/95, где в разделе «Основные принципы наблюдения при несении ходовой вахты» говорится следующее: «Предстоящий рейс должен быть спланирован заранее с учетом всей информации, а любой проложенный курс должен быть проверен до начала рейса».

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

а) неверный расчет точки поворота (точки начала перекладки руля);

б) неучет влияния мелководья.

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

Сначала рассмотрим ситуацию, когда для конкретных условий плавания выбран правильный радиус циркуляции, но расчет точки поворота выполнен неверно. Иными словами, неверно рассчитано расстояние упреждения (или упреждение), заключенное между точкой поворота и точкой пересечения старого и нового курсов. (В некоторых англоязычных источниках это расстояние называется «расстояние до нового курса» - distance to new course). Помимо обычного промаха в расчетах и графических построениях такая ошибка может быть вызвана применением различных методов определения точки поворота. (Анализ таких методов, оставленный за рамками данной статьи, показывает, что они дают различные результаты). Попытаемся ответить на вопрос о том, насколько судно отклонится от линии намеченного курса (пути), если судоводитель ошибется в выборе точки начала поворота. Эта ситация представлена на рис. 1. Введём следующие обозначения: РРо = ∆р – погрешность в расстоянии упреждения; CA = ∆uо поперечное смещение с линии нового курса; АСо = lо продольное смещение по курсу; ССо = ∆со – ли­нейное смещение от точки выхода на новый курс; < САСо = 90°; < ССоА = К2 К1 = ∆ К – угол поворота.


Очевидно, что поперечное смещение с линии нового курса можно выразить в виде

uо = ∆со∙sin∆ К = ∆р ∙sin∆ К (1)

Тогда продольное смещение будет определяться как

lо = ∆р ∙cos∆ К ( 2)

Поперечное смещение uо возрастает с увеличением угла поворота, дости­гая максимального значения, равного ∆р, при повороте на 90о, затем начинает уменьшаться и становится равным нулю, когда судно ложится на обратный курс. Продольное смещение по курсу lо меняется обратным образом по срав­нению с uо.

Что касается линейного смещения от точки выхода судна на новый курс, то при любом угле поворота его величина равна погрешности в расстоянии до точ­ки начала поворота судна, т. е. ∆со = ∆р, и, кроме того, совпадает с ним по направлению.


Известно, что если глубина под килем становится менее пяти осадок судна, то это начинает сказываться на характеристиках управляемости судна. В частности, увеличивается радиус циркуляции. Следовательно, этот фактор необходимо учитывать при планировании и выполнении рейса, особенно если отношение глубины к осадке становится меньше 2÷3. Рассмотрим погрешно­сти, которые могут возникнуть, если не учитывать влияние мелководья при планировании и реализации поворотов. На рис. 2 пунктирная линия К2 – новый курс, на который нужно лечь судну. Предположим, что для этого на этапе пред­варительной прокладки была рассчитана точка перекладки руля таким образом, чтобы в точке Р судно начало поворот, а в точке С – вышло на новый курс. Од­нако под влиянием мелководья радиус циркуляции увеличивается и, фактиче­ски, судно будет двигаться по циркуляции, показанной сплошной линией. При этом судно выйдет на новый курс в точке СМ. В этом случае АСМ – поперечное смещение с линии курса (причем оно всегда направлено влево от выбранного курса при повороте вправо и наоборот), СА – продольное смещение по курсу, а ССМ – линейное смещение точки выхода на новый курс.

Если обозначить АСМ = uМ, СА = lо, ССМ = ∆сМ , R – радиус циркуляции, RМ – радиус циркуляции на мелководье, К2 – К1 = ∆ К – угол поворота, а также учесть, что

ООМ = RМ R,

АC = ОБ = ООМ ∙ sin ∆ К = (RМ R) sin ∆ К,

то после несложных преобразований получим

uМ = CMOMAOM = (RМ R)∙ (1 – cos ∆ К), (3)

(4)

lо = (RМ R) sin ∆ К (5)

Анализ этих зависимостей показывает, что, если при прокладке не учитывать увеличение радиуса циркуляции из-за влияния мелководья, то поперечное смеще­ние с линии нового курса uМ возрастает с увеличением угла поворота, стано­вится равным (RМ R) при повороте на 90о (точка CM90° на рис. 2) и достигает максимального значения, равного 2(RМ R), когда судно ложится на обратный курс. Продольное смещение по курсу достигает максимального значения, равного (RМ R) при повороте на 90о. Линейное смещение точки выхода на но­вый курс ∆сМ с увеличением угла поворота также возрастает, при повороте на 90о достигает значения, равного , при дальнейшем повороте судна про­должает увеличиваться и при выходе судна на обратный курс становится рав­ным как uМ как по величине, так и по направлению, фактически совпадая с по­следним.

Следует полагать, что при совместном действии двух рассмотренных факто­ров – а) неверном расчете точки перекладки руля и б) неучете или неверном уче­те влияния мелководья – общие значения как поперечного (u), так и линейного (∆с) и продольного смещения (l) будут равны сумме обоих составляющих, т. е.

u = ∆uо + ∆uМ , (6)

∆с = ∆со + ∆uМ, (7)

l = ∆lо + ∆lМ , (8)

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

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

Может возникнуть вопрос о правомерности рассмотренных закономерно­стей и достоверности полученных оценок, например, из-за того, что в реально­сти во время маневренного и эволюционного периодов циркуляция имеет пере­менный радиус кривизны. Однако здесь встает проблема расчета мгновенного радиуса кривизны эволюционной части циркуляции, а полученные соотноше­ния остаются справедливыми и в этом случае. В то же время, проблему расчета мгновенного радиуса можно решить, если эволюционную часть циркуляции ап­проксимировать в виде Архимедовой спирали. Тогда радиус кривизны любой точки эволюционной части циркуляции RЭ можно рассчитать следующим образом:

. (9)

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

Таким образом, можно сделать следующие выводы:

1)  при наличии погрешности в расстоянии упреждения расчетное положение судна будет смещено относительно его фактического положения, и величина этого смещения может достигать максимального значения, равного погрешности в упреждении;

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

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

Это обстоятельство необходимо учитывать при планировании маршрута плавания, особенно для крупных судов.

СИСТЕМАТИЗАЦИЯ ВЫРАЖЕНИЙ ДЛЯ РАСЧЕТА И РЕЗУЛЬТАТЫ

РАСЧЕТОВ ОСНОВНЫХ ПАРАМЕТРОВ ГИРОАЗИМУТКОМПАСА «ВЕГА» ПО РЕЗУЛЬТАТАМ ОБРАБОТКИ КРИВЫХ ЗАТУХАЮЩИХ КОЛЕБАНИЙ

, , ёсов,

МГУ им. адм. , г. Владивосток

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

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

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

Среди гирокурсоуказателей отечественный корректируемый гироазимуткомпас (ГАК) «Вега-М» занимает ведущие позиции по количеству установленных комплектов на отечественном транспортном (морском и речном) и рыболовном флотах, а также точности выработки курса.

Принципиально корректируемые ГАК имеют возможность решать следующие задачи [1]:

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

– получение приемлемой точности курсоуказания при маневрировании, для чего нужно увеличить период .

К существенным недостаткам ГАК «Вега-М», которые в основном вытекают из давности разработки самого ГАК, можно отнести:

– аналоговая схема управления переходными процессами ГАК без возможности перенастройки;

– длительное время прихода в меридиан (3–6 часов);

– электромеханическая система выработки корректирующих сигналов.

В лаборатории технических средств навигации кафедры технических средств судовождения Морского государственного университета им. адм. была произведена глубокая модернизация ГАК «Вега-М». В результате без существенных изменений осталась только гиросекция (в ней СКВТ курса заменен на цифровой энкодер).

Целью работы является систематизация материалов по основным расчетным формулам для параметров ГАК «Вега-М», обработка на основе этих формул кривых затухающих колебаний модернизированного ГАК «Вега-С», полученных экспериментальным путем, и сравнения полученных результатов с теоретическими для широты г. Владивостока (место установки ГАК «Вега-С»).

Принцип работы ГАК и теоретическая база достаточно полно изложены в работах [2–5].

Для оценки основных параметров ГАК при стендовых испытаниях использовались следующие выражения:

Уравнение для угла отклонения в азимуте [3]

(1)

или в нормальном виде

, (2)

где , .

Уравнение в форме Лапласа записывается в виде

, (3)

где и ,

или . (4)

Уравнения (1–4) описывают переходной процесс инерционного звена второго порядка [6]

, (5)

где Т1 и Т2 – постоянные времени звена. Уравнение (5) можно представить в другом виде:

, (6)

где , . l – коэффициент демпфирования (показатель затухания),

.

Постоянная времени звена равна [6]

. (7)

Переходная функция звена при подаче на вход единичного ступенчатого воздействия (отклонения оси чувствительного элемента ГАК из меридиана) и нулевых начальных значениях имеет вид

, (8)

где , – соответственно вещественная и мнимая части корней (5) или характеристического уравнения (6), h – коэффициент затухания, – частота затухающих колебаний оси чувствительного элемента ГАК. Постоянные времени звена можно представить в виде:

, . (9)

Для вычисления коэффициента затухания и частоты затухающих колебаний по результатам эксперимента можно воспользоваться выражениями:

, , (10)

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

Из уравнений (1–3) можно получить:

, , (11)

, , . (12)

Фактор затухания вычисляется по формулам [2]

. (13)

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

,

, , (14)

.

Аналитическая формула переходной характеристики инерционного звена второго порядка имеет вид [3]

, (15)

которая аналогична формуле (8).

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

В таблице 1 приведены величины основных параметров ГАК «Вега-М» для широты г. Владивостока ().

Таблица 1

Основные параметры ГАК «Вега-М»

Обозначение

Наименование

Значение

Размерность

Н

Кинетический момент гироскопа

0,4

Н·м·с

0,045

0,0171

1/с

0,00077

1/с

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

7,29·10-5

1/с

Tg

Постоянная времени ИГ

60

с

Ay

Модуль управляющего момента

6,84·10-3

Н·м/рад

Az

Модуль демпфирующего момента

3,08·10-4

Н·м/рад

Степень демпфирования (коэффициент затухания)

0,38·10-3

1/с

f

Фактор затухании в φ=0, 60, 80 градусов

3, 6, 30

Т0

Для широты Владивостока

6586/110

с/мин

Тd

Для широты Владивостока

7198/120

с/мин

kf

T0/Td

0,915

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

Рис. 1. Кривая затухающих колебаний

Используя выражения (3), (6) – (14) обработаны кривые затухающих колебаний. Результаты обработки сведены в таблицу 2.

Таблица 2

Результаты эксперимента

Параметры

Дата

k, град.

A1, град.

A2, град.

A3, град.

Td1, сек.

Td2, сек.

11.02.09

90

30,28

10,04

3,56

6130

6330

12.02.09

80

26,4

9,64

2,38

6231

5873

16.02.09

106

36,16

11,14

4,41

6065

6290

17.02.09

107

36,83

11,73

3,82

6150

6060

18.02.09-1

90

29,42

9,85

3,32

6100

6060

18.02.09-2

35

34,33

11,07

3,72

6100

6020

19.02.09-1

47

16,33

5,4

1,85

6035

6040

19.02.09-2

100

34,26

11,24

4,05

6095

6225

Среднее

6113

6112

10.03.09

135

33,28

8,42

2,02

7200

6980

Таблица 3

Результаты расчетов

Параметры

Дата

h, с-1

х10-4

ωd, с-1

х10-3

λ

T1=T, сек

T2, сек

T0, сек.

f

kf

11.02.09

3,492

1,025

0,323

923,5

595,7

5802

2,92

0,947

12.02.09

3,862

1,008

0,358

926,1

662,4

5808

3,39

0,932

16.02.09

3,469

1,036

0,318

915,3

581,3

5748

2,89

0,948

17.02.09

3,685

1,022

0,339

920,8

624,8

5784

3,1

0,941

18.02.09-1

3,577

1,03

0,328

917,0

601,7

5760

2,98

0,945

18.02.09-2

3,643

1,03

0,333

915,3

610,4

5748

3,04

0,943

19.02.09-1

3,609

1,041

0,327

907,5

594,4

5700

2,97

0,945

19.02.09-2

3,503

1,031

0,322

918,5

591,1

5772

2,0

0,947

Среднее

3,605

1,0279

0,331

918,0

607,7

5766

2,91

0,944

10.03.09

3,981

0,8727

0,407

1047

852,5

6576

4,06

0,913

Таблица 4

Результаты расчетов

Параметры

Дата

Ay, Н∙м/рад

х10-3

Az, Н∙м/рад

х10-4

Az/Ay

Td, сек

Среднее

8,923

2,884

0,032

1,068

0,767

6120

1,176

1,142

10.03.09

6,86

3,185

0,046

0,967

0,997

7090

1,015

1,01

Усреднённые вычисленные параметры модернизированного ГАК «Вега» не соответствуют параметрам серийного ГАК «Вега-М». После пересчета коэффициентов передачи по горизонтальному и вертикальному каналам получены экспериментальные данные за 10 марта 2009 г.

Вид экспериментальной кривой затухающих колебаний за 10.03.09 г. и рассчитанных по формулам (8) и (15) переходных процессов показан на рис. 2. Несовпадение кривых на начальном участке объясняется следствием того, что после режима «Горизонтирования» ГАК был выведен ускоренно из меридиана, а затем переведен в режим «Работа». Индикатор горизонта длительное время находился на упоре (t =25¸75 мин.), что исказило переходной процесс в начальный период.

Из таблицы 3 (эксперимент за 10.03.09) видно, что отклонения экспериментальных результатов от теоретических не превышает 3,3% по параметру Az и примерно 1% по остальным нормируемым параметрам, что вполне приемлемо для инженерных расчетов.

Рис.2. Вид экспериментальной КЗК и рассчитанных переходных процессов ГАК

Список использованных источников

1.  Новые технические средства судовождения [Текст] / и др. Под общ. Ред. . – М.: Транспорт, 1973. – 264 с.

2.  Коган, гироазимуткомпас «Вега» [Текст] / , . – М. : Транспорт, 1983. – 200 с.

3.  Смирнов, средства судовождения: Теория: Учебник для вузов [Текст] / , , / Под ред. . – М.: Транспорт, 1988. – 376 с.

4.  Блинов, приборы: Учебник для вузов ММФ [Текст] / , , . – 4-е изд., перераб. и доп. – М.: Транспорт, 1980. – 448 с.

5.  Перминов, параметров двухрежимного корректируемого гирокомпаса [Текст] / Приборостроение, сб. научн. тр. № 000. Пермь: 1974. – С. 30–34.

6.  Власенко, электроавтоматика [Текст] / , . М. : Транспорт, 1983. – 258 с.

ГИРОАЗИМУТКОМПАС «ВЕГА»

С АПЕРИОДИЧЕСКИМ ХАРАКТЕРОМ ПРИХОДА В МЕРИДИАН

, , ёсов,

МГУ им. адм. , г. Владивосток

Приводятся экспериментальные данные апериодического процесса прихода в меридиан корректируемого гироазимуткомпаса «Вега-С» с цифровой системой управления. Модули маятникового и демпфирующего моментов рассчитаны по заданным фактору затухания и времени прихода в меридиан.

В лаборатории технических средств навигации кафедры технических средств судовождения Морского государственного университета им. адм. была произведена глубокая модернизация гироазимуткомпаса (ГАК) «Вега-М». Разработана цифровая система управления гироблоком. Без существенных изменений осталась только гиросекция (в ней СКВТ курса заменен на цифровой энкодер).

Из за большого объема эта статья размещена на нескольких страницах:
1 2 3 4 5 6 7