Text
                    В. И. ЕРМАКОВ
Использование
программируемых
ми крокал ькуляторов
в судовых условиях

Внедрение малых вычислительных средств в практику работы судоводителей облегчает их труд и способствует повышению надеж- ности решения задач, связанных с контролем за плавучестью и остойчивостью судна, с расчетами инерционно-тормозных харак- теристик судна, с навигационными и эксплуатационными вопросами. На судах флота рыбной промышленности внедрены программи- руемые микрокалькуляторы типов „Электроника - БЗ-34, МК-52, -54, -56", имеющие примерно одинаковые технические данные. Опыт использования этих микрокалькуляторов в учебном процессе пока- зал, что для их эффективной эксплуатации необходимы удачно сос- тавленные программы и хорошее знание лексики и грамматики языка микрокалькулятора. В настоящем пособии сделана попытка помочь слушателям и пользователям квалифицированно освоить имеющиеся программы. ЛЕНИНГРАД „СУДОСТРОЕНИЕ" 1990
Ермаков В. И. Е72 Использование программируемых микрокалькуляторов в судо- вых условиях: Справочник. — Л.: Судостроение, 1990. — 104 с.: ил. ISBN 5-7355-0224-7 Кандидат технических наук с большим стажем преподавательской работы составил практическое руководство по эксплуатации программируемых микро- калькуляторов отечественного производства для решения задач, связанных с судовождением, с контролем остойчивости судна, определением его инерци- онно-тормозных характеристик. В пособии обобщен опыт использования таких приборов на судах флота и в учебных заведениях. При разработке программ учитывались по возможности следующие требования: удобство ввода исходных данных и получения конечных результатов; краткость и лаконичность про- грамм; высокая степень автоматизации способа расчета часовых углов солнца и звезд и склонения солнца; полнота сведений по оптимальному использованию микрокалькулятора. Все программы прошли детальную проверку и одобрены экипажами судов рыбопромыслового флота. Для судоводителей. Полезна также преподавателям и курсантам'морских учебных заведений. „ 2705140300-055 ,, „„ Е 048(01)-90 31 90 ЬГ.К 32.844 ISBN 5-7355-0224 ® Ермаков В. И.. 1990
ПРЕДИСЛОВИЕ В справочнике содержатся сведения о реализации возможностей программируемого микрокалькулятора „Электроника МК-52” для реше- ния навигационных, астронавигационных и эксплуатационных задач в судовых условиях. Каждая программа имеет инструкцию по ее использованию, а также сопровождается примерами и комментариями. В справочнике дана инфор- мация по применению блоков расширения памяти микрокалькулятора. В книге приводятся также сведения об использовании микрокомпьютера „Электроника МК-85”, служащего для выполнения научных, инженерных и статистических расчетов и прикладных задач с помощью программ, составленных на языке Бейсик. Книга предназначена для судоводителей транспортного и промысло- вого флотов. Однако она может быть полезна рреподавателям и учащим- ся судоводительской специальности. Для тех, кто желает пополнить свои знания, рекомендуется обратиться к пособиям [4,5, 7, 9]. Большую помощь в составлении программ и ценные замечания по улучшению содержания книги внес профессор Роман Андреевич Скубко, к сожалению, недавно ушедший из жизни. Отзывы и предложения можно присылать по адресу: 191065, Ленин- град, ул. Гоголя, 8, издательство „Судостроение”.
OCIIOIIHblf OliO JHA'U НИМ Геометрия и остойчивость судна L - длина судна, м В - ширина судна, м <1 - осадка судна (средняя) , м Cg --.коэффициент общей плотности судна Д - массовое водоизмещение, т S7 - объемное водоизмещение, м3 XQ - абсцисса центра массы судна a G - аппликата центра массы судна 1g - плечо статической остойчивости, м 1Р- !М ~ плечи остойчивости формы относительно основной, произвольного полюса, метацентра соответственно, м h - начальная поперечная метацентрическая высота, м в - угол крена с - инерционный коэффициент, с • м1/2 т - период собственных бортовых колебаний судна, с Астронавигация <Р - широта наблюдателя X - долгота наблюдателя 6 — склонение светила, девиация компаса гм - местный часовой угол h - высота небесного светила А — азимут светила г’ - часовой угол точки Овна г*- звездное дополнение S - расстояние, мили К - курс судна, град Д'Г - интервал времени, ч Кинематика и динамика судна v - Скороть судна К - гидродинамический коэффициент S - тормозной путь, м Д - расстояние, измеренное радиолокатором, мили m - масса судна, кг F - сила упора винта, м t — интервал времени, с
Глава 1 ИСПОЛЬЗОВАНИЕ МИКРОКАЛЬКУЛЯТОРА В НЕПРОГРАММИРУЕМОМ ВАРИАНТЕ 1.1. Назначение клавиш Пульт управления микрокалькулятора содержит 32 клавиши с двой- ным или тройным управлением. В непрограммируемом варианте исполь- зуются следующие клавиши (рис. 1): F, К - префиксные; 0 ... 9 цифровые; □ - отделение целой части от дробной; В1 - разделение вводимых чисел и передвижение информации в регистрах стека; , г, X — арифметические операции; — операции обмена содержимым между регистрами X и Y; Сх - сброс содержимого регистра X (индикато- ра); /—/- - операция смены знака мантиссы или порядка числа; ВП - подготовительная клавиша для ввода порядка числа; X -* П - запись числа в регистры памяти; П ->Х - вызов из регистра памяти. Префиксные клавиши F и К используются для реализации так назы- ваемых совмещенных функций. Это означает, что при нажатии, например, клавиш F и 8 реализуется функция косинуса угла, а К и S - функция перевода угла, выраженного в градусах и минутах, в градусы и доли градуса. Подробно об этом указано в руководстве по эксплуатации мик- рокалькулятора. 1.2. Выполнение одноместных операций Одноместными называются операции, в которых используется одно число: например, отыскание значения тригонометрической функции, Рис. 1. Внешний вид микрокалькулятора „Электроника МК-52”
логарифма или числа по логарифму, вызов числа я и др. Рассмотрим изложенное на примере нахождения синуса угла. Отметим, что при последовательном нажатии клавиши F и другой реализуется операция, указанная над соответствующей клавишей. В этом случае принято писать то обозначение, которое написано над клавишей. Для удобства клавиша F и надклавишное обозначение имеют одинаковый цвет - желтый. Так, например, при нажатии клавиш F и 7 реализуется функция синуса. Для нахождения синуса 20° необходимо: переключатель Г/Р/ГРД перевести в положение Г (крайнее левое положение); набрать на индикаторе число 20; нажать клавиши F и 7 и ждать результата. Через несколько секунд на индикаторе высветится число 3,4202014—01, равное 0,34202014. Следует помнить, что если число по абсолютному значению меньше единицы, то оно представляется в форме плавающей запятой. Так как любое число можно представить в виде N = т • 10", где т — мантисса, а п - порядок числа на основании 10, то число 0,34202014 будет записано как 3,42020-Ю-1. Однако на индикаторе высвечивается только порядок. Рассмотрим другой пример. Пусть на индикаторе находится число N = 3,4202014 - 01 и нужно найти его arcsin. Для этого необходимо нажать клавиши F и 4. Тогда после завершения цикла на индикаторе высветится число 20. Поскольку переключатель находится в положении Г, то, следовательно, мы получили 20°. Приведем еще несколько комбинаций клавиш с пояснениями: F В I — восстановление (на индикаторе) предыдущего результата. Так, для пре- дыдущего примера после нажатия клавиши на индикаторе получим число 3,42202014 — 01; F □ — кольцевое перемещение чисел в стековых реги- страх; F С F — сброс действия клавиши F. 1.3. Занесение чисел в регистры памяти и извлечение из них Микрокалькулятор имеет 15 регистров памяти. Запись числа в реги- стры памяти производится с индикатора следующим образом: набирается число (если оно не является результатом предыдущей операции); нажима- ются клавиши X ->П и символ адресуемого регистра (числовые О ... 9, буквенные а ... е). При этом копия числа с индикатора будет перене- сена в соответствующий регистр памяти, а само число останется на инди- каторе. Извлечение числа из любого регистра памяти на индикатор произво- дится нажатием клавиш П -> X и символа регистра. Для контроля (в на- чальный период работы с микрокалькулятором) рекомендуется после выполнения операции Запись числа в регистр памяти произвести операцию извлечения числа из того же регистра памяти. Дело в том, что иногда вместо клавиши X -* П (при записи числа в регистр памяти) ошибочно нажимают клавишу П -> X. В этом случае ошибку также можно обнару- жить сличением ошибочно вызванного числа (нажата клавиша П -* X) 6
с числом, вводимым в регистр памяти. Однако, как показывает опыт учебной работы, первый способ более надежен. Пример. Требуется занести число 1 2 3 в 4-й регистр памяти. Решение. Последовательно нажимаем клавиши 1 2 3 - на индекаторе высветится число 123, Затем нажимаем клавиши X— П 4. Число 123 будет записано в 4-й регистр памяти. Для контроля следует нажать клавиши П-*Х 4. Если число, извлеченное из 4-го регистра памяти, будет отличаться от исходного, то необходимо исправить ошибку одним из двух способов: 1) нажать клавишу ♦», т. е. вызвать заданное число из регистра Y, так как после нажатия клавиш П-*Х4 число, находившееся в регистре X, переместится в регистр Y; 2) набрать заново, после чего нажать клавиши X -> П 4. Примечание. При наборе чисел необходимо четко нажимать клавиши и задерживать клавишу на долю секунды в нижнем положении. Для вычисления тригонометрических функций микрокалькулятору требуется около 2 с, поэтому необходимо выдержать паузу перед очередной операцией (на- пример, перед операцией умножения синуса угла на число и т. д.). 1.4. Выполнение двухместных операций Арифметическая операция выполняется следующим образом: наби- рается первое число; нажимается клавиша В t (число, находящееся на ин- дикаторе, т. е. в регистре X, переводится в регистр Y); набирается второе число; нажимается клавиша соответствующей арифметической операции; с индикатора (регистр X) снимается результат. Пример. Требуется найти сумму двух чисел: Л\ = 5, IV, =7. Решение. Нажимаемые клавиши...... 5 В t 7 + Индикация................ 5 5 7 12 Если хотя бы одно из чисел вводится вызовом из памяти или получается как результат одноместной или предыдущей арифметической операции, то клавишу В 1 нажимать не надо. Возведение числа X в степень У, т. е. нахождение числа N= X?, производится так: на индикатор вводится показатель степени У; нажимается клавиша В t (число переводится в регистр Y); на индикатор число X; нажимаются клавиши F X Y. Пример. Дано: У = 2/3; Х= 2500. Требуется найти N = 2500 - 2/3. Решение. Нажимаемые клавиши ..... 2 В t 3 2500 F XY Индикация................... 2 2 3 6.666..-01 2500 184.20157 Ответ. N= 184.20157. Если при выполнении сложной арифметической операции необходимо запом- нить промежуточные результаты, их надо занести в регистры памяти. Использование стековой памяти. Стековая память микрокалькуля- тора содержит 4 регистра: Т, Z, Y, X, из которых X и Y являются опера- ционными, т. е. в арифметической операции участвуют числа, находящие- ся в регистрах X и Y (отметим, что делитель и вычитаемое находятся в регистре X). Дополнительно к ним имеется регистр предыдущего ре- зультата — X].
Ввод числа нажатием цифровых клавиш или вызовом из соответст- вующего регистра памяти, как указано ранее, всегда производится в ре- гистр X. После ввода очередного числа информация в стеках перемещается на один регистр вверх. Если нажимается только клавиша В t, то число, находившееся в регистре X, сохранится, а его копия перейдет в регистр Y. Схематически перемещение чисел в стеках будет выглядеть так: Клавиши Схема перемещения Цифровые.............Х-> Y-+ Z -*Т-> исчезает Bt...................Х-> Y-> Z-*Т-► исчезает После выполнения арифметической операции числа перемещаются на один регистр вниз. Схематически это выглядит так: Т ->-Z -> Y ->|*| ->Х^ -»х1. Знак * здесь означает результат арифметической операции, произве- денной с содержимым регистров X и Y. При нажатии клавиш происходит кольцевое перемещение чисел: X->T->Z ->Y->X. В качестве примера работы стека рассмотрим перемещение информа- ции в регистрах стека. Нажимаемая 1 В t 2 Bt 3 Bt 4 F 5 клавиша Отметим, что умелое использование возможностей стековых регист- ров позволит избежать лишних действий при составлении программы вычислений. Решение вспомогательных задач. Микрокалькуляторы „Электрони- ка МК-52 и МК-61” имеют встроенные подпрограммы для перевода вре- менных величин, выраженных в часах, минутах, секундах и долях секунд, в часы и обратно, а также угловых величин, выраженных в градусах, минутах и долях минут, в градусы и доли градуса и обратно. Рассмотрим несколько примеров. Пример. Дано: а = 120° 34.7'. Требуется найти значение угла, выраженное в градусах и долях градуса. Решение. Нажимаемые клавиши Индикация 1 2 3 120. И 120. 3 4 7 120.347 К+ 120.57832 Ответ: а = 120.57832°. Пример. Дано а= 36.12345°. Требуется найти значение угла в градусах, мину- тах и долях минуты. 8
Решение. Нажимаемые клавиши 36 □ 1 2 345 Кб Ответ: а. = 36° 07.407'. Индикация 36. 36. 36.12345 36.07.407 Пример. Дано: Ггр = 12 ч 37 м 42 с. Требуется найти значение Гринвичского времени в часах и долях часа. Решение. Нажимаемые клавиши Индикация 1 2 12. □ 12. 3 74 2 12.3742 Ко 12.628332 Ответ: Тг„ = 12.628332 ч. Пример. Дано: ТГр = 13,48971 ч. Требуется найти значение Тур в часах, минутах и секундах. Решение. Нажимаемые клавиши Индикация 489 7 1 13.48971 К 3 13.292296 Ответ: Тгр = 13 ч 29 м 23 с. Глава 2 ИСПОЛЬЗОВАНИЕ МИКРОКАЛЬКУЛЯТОРА В РЕЖИМЕ „ВЫЧИСЛЕНИЕ ПО ПРОГРАММАМ” 2.1. Общие положения В режиме „Вычисление по программам” дополнительно использу- ются следующие клавиши, реализующие соответствующие команды: F ПРГ - перевод в режим программирования; F АВТ — перевод в режим автоматических вы- числений (вывод из режима про- граммирования) ; F (X < 0, X = 0, X > 0, X т6 0) - ввод операций условных переходов; БП - ввод операции безусловного пере- хода; ПП — ввод подпрограмм (и потактовое про- t хождение программы); ШГ - смещение программы на один шаг „вперёд”; 9
ШГ - смещение программы на один шаг „назад”; В/0 - установка программы по нулево.му адресу и ввод опера- ции возврата из подпрограммы; С/П - операция Останов и пуск вычислений по программе; F L (0 1 2 3) - ввод операций для организации циклов; К НОП - ввод команды НЕТ ОПЕРАЦИИ (исключение ошибочной команды); К - реализация косвенной адресации (совместно с клавиша- миП->ХиХ-*Пи номерами адресуемых регистров). Программа представляет собой последовательность кодов операций (команд), которые микрокалькулятор выполняет автоматически в про- цессе вычисления. При этом имеется ввиду, что исходные данные занесе- ны в соответствующие регистры микрокалькулятора. 2.2. Программирование задачи Рассмотрим простейшую задачу, которая, однако, имеет практическое значение при решении астрономических задач на микрокалькуляторе „Электроника МК-56”. Речь идет о переводе угла, выраженного градуса- ми и минутами, в градусы. Формула в этом случае имеет вид а0 = а° +------а2. 60 (2.1) где а° — число градусов; а'2 — число минут с десятыми долями. Прежде всего проанализируем формулу и выявим константы и пере- менные. Из формулы видно, что в исходных данных содержится одна константа - число 60 и две переменные, одна из которых (at) выражена в градусах, а другая — в минутах. После этого необходимо определить регистры, в которых будут находиться исходные данные перед началом решения задачи, а также ре- гистр, где будет храниться результат (в любом случае он прежде всего попадает в регистр X). Для этой цели составим таблицу. Из табл. 2.1 видно, что число 60 занесено в регистр 0, это позволит использовать его многократно. Часть угла, выраженного в градусах, Параметр Обозначения Символ регистра Константа Часть угла, выраженная в градусах Часть угла, выраженная в дуговых минутах Угол, выраженный в градусах и долях градуса 10
будет находиться в регистре Y, другая часть, выраженная в минутах, — в регистре X. Полагаем, что результат должен оставаться в регистре X. Если бы возникла необходимость занести полученный результат в реги- стр 1, то вместо X было бы написано число 1. Теперь, используя формулу (2.1) и табл. 2.Г, составим программу для перевода углов в градусы. Программа 1 00 01 02 03 П -+Х0 v + С/П 60 13 10 50 Здесь первая строка — адреса команд, вторая — нажимаемые клави- ши, третья - коды команд. По адресу 03 записана команда ОСТАНОВ. Затем следует составить инструкцию по использованию программы, т. е. подробный перечень последовательных действий, необходимых для получения результата. Инструкция по использованию Программы 1 Нажимаемые клавиши F ПРГ F АВТ 60 X В/0 С/П 1. Включить микрокалькулятор 2. Перейти в режим „Программирование” 3. Занести программу 4. Перейти в режим „Автоматическая ра- 5. Занести исходные данные в регистры: 60 в регистр 0 О; в регистр Y аа в регистр X 6. По трем параметрам, занесенным в ре- гистры в качестве исходных данных, опреде- лить угол, выраженный в градусах Примечание. В действительности рассматриваемая задача является вспо- могательной и для ее решения инструкцию не составляют. Более того, в микрокаль- куляторах „Электроника МК-54”, „Электроника МК-52”, „Электроника МК-61” имеются встроенные подпрограммы, служащие для решения такой задачи (см. примеры на с. 18), Однако в качестве иллюстрации этапов вычисления по програм- мам она удобна. 2.3. Ввод и редактирование программы Этот этап работы выполняется в соответствии с Инструкцией по использованию программы. Перед включением микрокалькулятора необ- ходимо убедиться в том, что переключатель Г/Р/ГРД находится в левом положении, т. е. в положении Г. После выполнения п. 2 Инструкции на индикаторе высветится на- чальный адрес (в данном случае 00), с которого будет записываться программа. При выполнении п. 3 Инструкции на индикаторе последовательно будут появляться коды 60, 13, 10, 50, которые соответствуют адресам 11
00, 01, 02, 03. После появления кода 50 код 60 исчезает, т. е. максималь- ное количество одновременно индицируемых кодов равно трем. Таким образом, индикатор после ввода программы будет выглядеть так: 50 10 13 04. Если при вводе программы допущена ошибка и замечена сразу (по сличению кодов в Программе 1 и на индикаторе), то нажимают клавишу ШГ (перемещают программу „назад” на один шаг) и далее нужную клави- шу. Если же ошибка находится далеко от текущего адреса, то данный способ (перемещение по одному шагу) будет неэффективным. В этом случае поступают следующим образом: переводят микрокалькулятор в режим „Автоматическая работа” нажатием клавиш F АВТ; устанавливают счетчик адресов на номер ошибочной команды нажа- тием клавиш БП N, где N - адрес ошибочной команды; переводят микрокалькулятор в режим „Программирование” нажа- тием клавиш F ПРГ; исправляют ошибку. При этом, если имеется лишняя команда, ее устраняют нажатием клавиш К ИОП, а если пропущена команда, то набирают часть программы заново, начиная с пропущенной команды. Вернувшись к нашему примеру, предположим, что при наборе предпоследней команды вместо клавиши ES была нажата клавиша EL Тогда при сличении кодов вместо числа 10 на индикаторе микрокалькулятора будет обнаружено число 11. Для ис- правления ошибки дважды нажимают клавишу ШП, а затем клавишу И. Для отработки навыка можно нажать последовательно следующие клави- ши: F АВТ; БП 02; F ПРГ; +. После выполнения указанных действий ошибка будет устранена. 2.4. Занесение исходных данных и решение задачи Отладка программы заключается в том, чтобы решить аналогичную задачу с известным ответом. Выполнение этого этапа позволяет устано- вить, правильно ли составлена и занесена программа, исправен ли микро- калькулятор. Для проверки можно решить следующую задачу. Известно, что 20о30р = 20,5°, Поэтому, задав а = 20°30' и выполнив п. 5.6 Инструк- ции, т. е., нажав последовательно клавиши 20 В t 30 В/0 С/П, получим 20,5°. Это будет свидетельствовать о правильности решения контроль- ного примера и об исправности микрокалькулятора. Последний этап мало отличается от предыдущего по составу операций. Разница заключается в том, что ответ неизвестен и мы считываем его с индикатора. При этом константа 60 заново не вводится. Теперь решим примеры на микрокаль- куляторе МК-52. 12
Пример. Дано: <рс = 38O38’N; 6 = 16° 19*9'3; гм = 3° 11,8'W, Требуется выра- зить исходные данные в градусной мере и изменить знак 6. Решение. Нажимаемые клавиши Индикация 38.38 К + 38.633333 16.19 9 К + /—/ -16.331667 3.11,8 К + 3.1966667 2.5. Оптимизация программы командами переходов и организации циклов 2.5.1. Использование команды условного перехода. Рассмотрим составление более сложной программы, в которой используется команда условного перехода, - программы для расчета сличимой высоты и ази- мута светила. Исходные формулы [2]: tg Ас = sinlM/ (tg 6 • cos </>с - cos tM • sin<pc); cos hc = sinlM cos 6 cpsec Ac, где Ac - счиспимый азимут светила в полукруговом счете; hc - счисли- мая высота светила; <рс — счислимая широта наблюдателя; 8 — склонение светила; tM — практический часовой угол светила. Отметим, что азимут получается в полукруговом счете, а счислимая высота светила вычисляется через азимут. При этом синус (косеканс) азимута определяется через функцию его тангенса, а здесь содержится „ловушка”. Дело в том, что при Ас > 90° тангенс азимута имеет знак минус. Следовательно, мы получим вместо азимута его дополнение до 180 со знаком минус. Последующее вычисление синуса азимута, где в качестве аргумента будет использоваться полученный угол, приведет к ошибочному резуль- тату. Поясним это на примере. Нажимаемые клавиши Индикация 1. Включаем микрокалькулятор, переклю- чатель ставим в положение Г 2. 120 120. 3. F9 -1.7320508 4. F6 -59,99999 5. F7 -8,6602544 - 01 Таким образом, после выполнения операции нахождения синуса угла мы получили ошибочный ответ, который приведет к ошибочному результату при вычислении счислимой высоты. Если же для вычисления sin^c использовать аргумент Ас = 120°, то sin/c = 0,86602544. 13
Следовательно, для исправления указанной ошибки надо ввести в программу такое условие: если азимут получится отрицательным, то к нему надо прибавить 180°. Забегая вперед, укажем, что в программе (см. далее Программу 2) это условие реализовано по адресам 20 ... 25. Фрагмент программы выглядит так: 20 21 22 23 24 FXC 0 24 П-»Хс + Х-П 5 5С 24 6С 10 45 На языке микрокалькулятора это означает: если полученный угол (после выполнения команды по адресу 19) будет меньше нуля (т. е. иметь отрицательное значение), то условие выполняется. Следовательно, микрокалькулятор пропускает адрес перехода (команду, записанную по адресу 21) и начинает выполнять программу с 22-го адреса: извле- кает содержимое из регистра С (180°), складывает его с полученным значением угла (т. е. вычисляет азимут) и заносит его копию в регистр 5. Оставшееся на индикаторе число (азимут) используется дальше по назна- чению. Если же после выполнения команды по адресу 19 значение угла бу- дет больше нуля (положительное значение), то микрокалькулятор после проверки условия (условие не выполняется) перейдет на использование команды по адресу 21, т. е. занесет азимут в регистр 5. Изложенное долж- но насторожить пользователей, желающих самостоятельно составить про- грамму для решения задачи с использованием тригонометрических функ- ций и их комбинаций. Теперь перейдем к написанию программы, предварительно выделив для констант, переменных, а также полученных результатов соответствую- щие регистры памяти. Затем на основе исходной формулы и табл. 2.2 составляем програм- му с учетом необходимости использования условного перехода и целе- сообразности наличия программы для перевода углов в градусную меру, о чем говорилось ранее. Для того чтобы безошибочно составить програм- му, рекомендуется использовать „развертку” стековых регистров так, как это показано ниже в схеме передвижения информации в стековых регистрах: т о о о о о Z О 0 0 0 о Y 0 0 cos ipc cos >рс 0 X cos ‘Pc s tg S cos ¥>c tg 5 Клавиша П->Х1 F8 П-»Х2 F9 X Из фрагмента схемы видно, как перемещается информация в регист- рах стековой памяти в ходе реализации программы. Поскольку при сос- тавлении программы, как правило, возникают „опечатки”, рекомендует- ся пользоваться в записи мягкими карандашами. После завершения написания „развертки” на основе данных послед- ней строки (нажимаемых клавиш) оформляют программу, т. е. „стыку- ют” адреса и записывают коды. 14
Таблица 2.2. Исходные данные и результаты вычислений счнслимой высоты и азимута Параметр Обозначение (значение) Символ регистра Константа 60 0 Константа 180 с Широта счислимая •Рс Склонение светила 2 Местный часовой угол светила 'м 3 Азимут светила 5 Высота счислимая Х/6 Для оперативного контроля при записи программы полезно запом- нить формулу наиболее употребительных кодов. Так, например, формула кода записи числа в какой-либо числовой регистр памяти имеет вид NK = 4NP, где NK.~ номер кода; 4 - цифра, означающая запись в регистр памяти; 7Vp — номер числового регистра. Так, например, код записи числа в 5-й регистр памяти соответствует числу 45. Код извлечения числа из числового регистра памяти соответ- ствует формуле NK = 6Np, где NK - номер кода; 6 - цифра, означающая извлечение числа из регистра памяти; Np — номер числового регистра. Арифметические коды запоминаются легко: + — 10; X - 12; ----11; -г - 13. Коды записи и извлечения из буквенных регистров памяти соответст- венно имеют вид Далее представлена программа расчета счислимой высоты и азимута светила. Программа 2 00 01 02 03 04 05 06 07 08 п-хо + с/п П-Х 1 F8 П-Х2 F9 X 60 13 10 50 61 1Г 62 1Е 12 09 10 11 12 13 14 15 16 17 п-хз F8 П-Х1 F7 X — П-Х 3 F7 4- 63 1Г 61 1С 12 11 63 1С 13 18 19 20 21 22 23 24 25 26 F т F6 F Х< 0 24 П-Х с 4- Х-П5 П-ХЗ F7 23 1L 5С 24 6С 10 45 63 1С 27 28 29 30 31 32 33 34 -35 П-Х2 F8 X П-Х5 F7 — F5 Х-П 6 С/П 62 1Г 12 65 1С 13 1 — 46 50 15
Здесь содержатся две программы: первая (адреса 00 . . . 03) служит для перевода углов в градусную меру, а вторая (адреса 04 . . . 35) для расчета счислимой высоты и азимута. При этом для вычисления азимута используются команды с адреса 03 по адрес 24, т. е. когда азимут записы- вается в регистр 5. Для вычисления счислимой высоты используется ос- тавшаяся часть программы с адреса 25. Инструкция по использованию Программы 2 Нажимаемые клавиши 1. Включить микрокалькулятор, пере- ключатель установить в положении Г (граду- дусы) 2. Перейти в режим „Программирова- F ПРГ 3. Занести в программу 4, Перейти в режим „Автоматическая pa- F АВТ 60 X 180Х БП 04 С/П 5. Занести исходные данные для реше- ния контрольного примера: для микрокалькулятора „Электроника МК-56”: 60 в регистр 0 180 в регистр С ¥>с в регистр 1 6 в регистр 2 гм в регистр 3 6> Вычислить счислимую высоту и ази- мут светила 7. Сличить полученные данные с конт- рольными значениями и, если они совпадут, перейти к решению задачи, начиная с третьей позиции п. 5. Примечание. 1.Для перевода дробной части высоты необходимо вычесть целое число градусов, а оставшуюся часть умножить на 60; число 60 вызывается из регистра 0. 2. Широта считается всегда положительной. Склонение считается положитель- ным, если оно одноименно с широтой <р, и отрицательным, если разноименно с ¥>. Пример. Дано: = 38° 38,0’N; 6 = 16° 19,9'N; rM = 3° 11,8'W. Требуется вы- Решение. Предположим, что выполнены п. 1 -4. Инструкции. Поэтому пере- ходим в соответствующие регистры памяти (60 X -» П 0 для МК-56): 38.38 К + (38,633333) Х-П1 16.199К+ (16,331667) Х->П2 3.188 К + (3,1966667) Х->ПЗ После ввода исходных данных нажимаем клавиши БП 04 С/П и ждем заверше- ния цикла вычислений (примерно 25 с). На индикаторе высвечивается значение счислимой высоты: Лс = 67,52286°. Для отделения десятых долей градуса и перевода в дуговые минуты поступаем следующим образом: нажимаем клавиши К 6; инди- кация - 67.31 37, Следовательно, Лс = 67° 31,4'. Для вызова азимута из регистра 5 на индикатор нажимаем клавиши П--Х5, тогда на индикаторе высветится число 171,96377. Следовательно, Ас = 172°. По- скольку первая буква азимута полукругового счета одноименна с широтой, а вто- рая - с наименованием полукругового часового угла, то окончательно получим Ас = = 172° NW. 16
Теперь рассмотрим программу, в которой вместо обсервованной высоты вычисляется перенос высотной линии положения. Так как эта программа незначительно отличается от предыдущей, то ограничимся лишь описанием различий в программах. Размещение исходных данных и конечных результатов показано ниже: 60 о 180 ¥>с йс 6 6 2 П X Здесь в первой строке указаны исходные данные и конечные результаты; во второй — регистры. Программа вычислений переноса высотной линии положения будет отличаться наличием дополнительных команд, а именно: 35 П-Х4 64 38 П-ХО 60 39 40 X С/П 12 50 Видно, что программа увеличилась на 5 адресов, но обрела иное ка- чество удобство пользования, которое уменьшает вероятность возник- новения ошибки при расчете переноса высотной линии положения обыч- ным способом. В инструкции по использованию этой (модифицированной) програм- мы в п. 5 добавится одна строка (для ввода йо) — й0 - в регистр 4: йо • йо К + X -> П 4. После завершения цикла вычислений на индикаторе высветится перенос высотной линии положения, выраженный в дуговых минутах. Если для предыдущего примера примем й0 = 67° 33,4', то после вы- полнения п. 6 Инструкции получим: П = +2,0г. 2.5.2. Использование команды безусловного перехода. Команда безусловного перехода реализуется клавишами БП и адреса перехода. При выполнении команды безусловного перехода прерывается последо- вательность выполнения команд в возрастающем порядке адресов и микрокалькулятор переходит на выполнение команды по указанному в программе адресу перехода. Наиболее часто безусловный переход ис- пользуется для облегчения работы пользователя исключением .лишних” нажатий клавиш при решении задачи. Так, например, если в программе 1 после адреса 03 написать команду безусловного перехода БП 00, то при решении второго и последующих примеров клавишу В/0 нажимать не надо. 2.5.3. Использование команды организации цикла. Команда органи- зации цикла используется при многократном выполнении одних и тех же операций. Она реализуется с помощью клавиш F L (0 ... 3), как было указано ранее, и адресов перехода. Команда, записанная по адресу перехода, является началом цикла, а команда F L (0 ... 3) вместе с адресом перехода — концом цикла. При каждом обращении к концу цикла из числа, находящегося в регистре 0 ... 3 (смотря по тому, кото- рый из них использован), вычитается единица, и если разность не равна 17 664-2
нулю, то процесс вычислений начинается снова с выполнении ••*«••011.1. записанной в адрес перехода. Количество циклов равно Ко/Ш'и »> еди- ниц, содержащемуся в упомянутом регистре. После этого мю- )«>ь.1льку- лятор, игнорируя адрес перехода, перейдет к выполнению команды, находящейся за адресом перехода. Более подробно эта процедура в комментариях к программе „Расчет средней квадратической погреш- ности”. Глава 3 РЕШЕНИЕ ОСНОВНЫХ ЗАДАЧ, СВЯЗАННЫХ С КОНТРОЛЕМ ПЛАВУЧЕСТИ И ОСТОЙЧИВОСТИ СУДНА 3.1. Расчет координат центра массы судна Первый вариант. Исходные формулы [4] имеют вид До*со+ ^Mxi у* xG = --------------------------------= ------; (3.1) П д д До + Sp/ 1 AozGo+Sp/x,- Щ. 1 о Mz zg =--------------------=------------= — ; (3-2) n Д Д До + S pt 1 Д=Д0 + 2р/, (3.2) где xG,zG - абсцисса и ордината центра массы судна в данном случае нагрузки, м; xGq, zGq - абсцисса и ордината центра массы судна порож- нем, м; До, Д — массовое водоизмещение судна порожнем И в данном случае нагрузки соответственно, т; р/ - массы статей нагрузки судна, т; Xj.Zj — координаты центра массы соответствующих статей нагрузки, м; п - количество статей нагрузки; М Л| ciawici кин момс«и i и статьи нагрузки относительно миделя, тс • м; Л/х статический момент i -й 18
статьи нагрузки относительно основной, м; Мх — статический момент массы судна относительно миделя, тс • м; Mz - статиче.ский момент мас- сы судна относительно основной, т • м. В данном варианте программы (см. Программу 3) ввод исходных данных (р/, х/, z/) производится пользователем в адресуемые регистры памяти (1, 2, 3 соответственно). По сравнению со вторым вариантом (см. Программу 4) при решении задачи следует в каждом цикле вычис- лений выполнять на 4 действия больше, т. е. производить больше нажатий на клавиши. Однако, как показала практика обучения работе на микро- калькуляторе, этот вариант программы, поскольку он проще, воспри- нимается пользователями на начальной стадии обучения лучше. В табл. 3.1 представлено распределение исходных данных, промежуточных и конеч- ных результатов. Таблица 3.1. Распределение исходных данных, промежуточных и конечных результатов по соответствующим регистрам памяти Параметр Обозначение Регистр Масса i -й статьи нагрузки, м Pi 1 Абсцисса i -й статьи нагрузки, м Xi 2 Аппликата 1 -й статьи нагрузки, м 3 Статический момент массы i -й статьи нагрузки MXj 4 относительно миделя, т • м Статический момент массы i -й нагрузки относи- MZi 5 тельно основной, т м Статический момент массы судна относительно Mx 6 миделя, т • м Статический момент массы судна относительно Mz 7 основной, т • м Абсцисса центра массы судна, м XG 8 Аппликата центра массы судна относительно ос- *G 9 новной,м Массовое водоизмещение судна, т Д a Расчет координат xG и zG центра массы судна приведен в Програм- ме 3. Программа 3 00 01 02 03 04 05 06 07 08 П--Х 1 П-Х2 X Х-П4 П-Х6 + Х-П6 П->Х1 П-Х 3 61 62 12 44 66 10 46. 61 63 09 10 12 13 14 15 16 17 X Х-»П5 П-Х7 + Х-П7 П-Х 1 П-Ха + Х-Па 12 45 67 10 47 61 6- 10 4 18 19 20 21 22 23 24 25 26 С/П П-Х6 П-Ха 4- Х-П8 П-*Х7 П-Ха 4 Х-»П9 50 27 С/П 50 66 6- 13 48 67 6- 13 49 19
Инструкция по использованию первого варианта Программ* 1 1. Включить микрокалькулятор 2. Перевести микрокалькулятор в режим „Про- граммирование” 3, Занести в память микрокалькулятора про- грамму 4. Перевести микрокалькулятор в режим „Ав- томатическая работа” 5. Произвести очистку регистров 6, 7, а (для при- обретения навыка выполнения этой операции в необ- ходимых случаях, например при повторном решении задачи и т.д.) 6. Произвести первый (контрольный) цикл вы- числений Ввести исходные данные в соответствующие регистры памяти Произвести цикл вычислений Проверить содержимое 4-го регистра памяти и сличить со значением Мх Проверить содержимое 5-го регистра памяти и сравнить его со значением Mz Проверить содержимое 6-го регистра памяти и сличить его со значением МХо Проверить содержимое 7-го регистра памяти и сличить его со значением Mz Произвести цикл вычислений координат и сли- чить полученное значение со значением zq Извлечь содержимое 8-го регистра памяти и сличить его со значением xq 7. Произвести второй цикл вычислений Ввести исходные данные в соответствующие регистры памяти П Вычислить моменты и массу Извлечь из соответствующих регистров памяти Мх и Mz для записи в таблицу расчетов (по форме, принятой в Информации об остойчи- вости судов) 8. Выполнить, как указано в п. 7, операции для всех статей нагрузки, имея в виду, что статья „Судно порожнем” уже введена 9. Извлечь значения Мх и Mz для записи в таб- 10. Произвести расчет xqwzq 11. Извлечь значения xq, z q и Д । павиши I HIT I- Л ВТ rf/o с/п БП 19 С/П В/0 С/П -Х4 (Мх ) + Х5 (Mz ) В/0 С/П П-Х 4 (МХ1) П-Х5 (Л/.,) П - X 6 (Мх) П -X 7 (М.) БП 19 С/П П-Х» (Д) Примечание. При снятии груза (или для исключения ошибочного цикла вычислений) масса груза вводится со знаком минус, я координаты пводягся со сво- ими знаками. Затем производится цикл вычислений нажпгигм клавиш В/0 С/П. В результате описанная выше операция будет выполнена. 20
2. При пересчете координат центра массы судна после расходования части груза, а также в случае пополнения запасов рекомендуется следующий способ: при выполнении п. 6 Инструкции использовать данные судна в предыдущем случае нагрузки (Д, xq, zq)\ снять тс статьи нагрузки (с их прежними данными), которые подвергнуты изменению; ввести снова те же статьи нагрузки сданными по состоянию на момент расчета; вычислить координаты центра массы судна после коррекции нагрузки. Второй вариант. Рассмотрим другой вариант программы для расчета координат центра массы судна. Программа 4 00 01 02 03 04 05 06 07 08 09 Х-ПЗ - Х-П2 - FX=0 09 F Х->П2 F Х-П1 43 14 42 11 5 Е 09 25 42 25 41 10 11 12 13 14 15 16 17 18 19 П-ХЗ X Х-*П5 П-Х7 + Х-П7 П->Х1 П-»Х2 X Х-П4 63 12 45 67 10 47 61 62 12 44 20 21 22 23 24 25 26 27 28 29 П-Х6 + Х-П6 П-»Х1 П-Ха + X-Па П-Х4 С/П П-Х7 66 10 46 61 6- ’ 10 4- 64 50 67 30 31 32 33 34 35 36 37 П->Ха -г- Х-»П9 П->Х6 П->Ха 4 Х-»П8 С/П 6- 13 49 66 6- 13 48 50 Таблица 3.2. Распределение исходных данных и промежуточных результатов Параметр Обозначение Регистр Масса i -й статьи нагрузки, т Pi Z/! Абсцисса i -й статьи нагрузки, м *i Y/2 Аппликата 1 -й статьи нагрузки, м zt X/3 Статический момент массы i -й нагрузки относи- тельно миделя, т • м 4 Статический момент массы i -й нагрузки относи- тельно основной, т • м Mz/ 5 Статический момент массы судна относительно основной, т • м 6 Статический момент массы судна относительно основной, т м M2 7 Абсцисса центра массы судна в грузу,м XG 8 Аппликата центра массы судна в грузу, м ZG 9 Массовое водоизмещение судна, т Д a Инструкция по использованию программы 4 Нажимаемая клавиша 1. Включить микрокалькулятор (переключатель в любое положение) 2. Перейти в режим „Программирование” F ПРГ 3, Записать программу в память микрокалькуля- 4. Перейти в режим „Автоматическая работа” F АВТ 21
5, Произвести очистку регистров 6, 7, а 6. Произвести первый (он же контрольный цикл вычислений Ввести исходные данные в соответствующие стековые регистры памяти калькулятора Произвести пуск программы с адреса 00 Проверить содержание регистров памяти 1 ... 3 (в них должны содержаться соответственно Ao.*Go.zGo) Проверить содержимое регистров 4 и 5 (в них должны содержаться МХо иМ,о соответсгвен- Проверить содержимое регистров 6 и 7 Произвести пуск программы с адреса 29 Снять с индикатора значение xq^ Извлечь из регистра 9 значение zq^ 7. Произвести второй и последующие циклы вы- числений Ввести исходные данные в стековые регистры памяти Произвести пуск программы с адреса 00 Снять с индикатора значение My Извлечь содержимое регистра 5 8. Извлечь содержимое регистров 6 и 7 9. Произвести расчет xq и zq 10. Снять с индикатора значение xq 11. Извлечь из регистра 9 значение zq 12. Извлечь из регистра а значение Д й. II 1 I,. И ' В/0 €711 II • X 1 II • X 2 П-ХЗ П-Х4 П-»Х5 П-Х6 (Мх) П-»Х7 (Mz ) С/П (БП 29 С°/П) П-»Х9 Pi В t Xj В t Zi В/0 С/П П-Х5 (Mz.) П-+Х6 (MXJ П-Х7 (Мг) С/П П-*Х9 П-Ха Краткое пояснение к программе. Исходные данные (р^, х, Zf) вводят- ся в стековые регистры (X, Y, Z) и сохраняются в регистрах 1, 2, 3 соот- ветственно до следующего цикла вычислений. Это сделано с той целью, чтобы можно было осуществить контроль после вычисления статических моментов. Следует иметь в виду, что очистка регистра X перед вводом данных не требуется, а если после ввода г,- будет ошибочно нажата кла- виша В t, то сработает блокировка и ошибка будет устранена без вмеша- тельства пользователя. Сущность блокировки сводится к следующему. Если содержимое регистра X (команда по адресу 03) будет равно нулю, то. следовательно, в регистрах X и Y содержались одинаковые числа — Zj. В этом случае командой FQ (адрес 06) будет перемещен в регистр Т, а в регистр X переместится х,, в регистр Y - Z/, т. е. ошибка будет исклю- чена, процесс вычислений не прервется. Если же исходные данные будут введены правильно, то разность будет отлична от нуля (за исключением случая, когда х/ = z(-, что практически исключено), следовательно, коман- ды по адресам 06 ... 08 выполняться не будут. Для исключения ошибочного цикла вычислений необходимо ввести массу р[ со знаком минус, а координаты со своими знаками и произ- вести пуск программы с адреса 00. Эти же операции выполняются при снятии груза. 22
Для отработки навыков по обоим вариантам программы рекомен- дуется использовать типовые варианты расчета координат центра массы, имеющиеся в информации об остойчивости судна, т. е. примеры с готовы- ми ответами. 3.2. Расчет координат центра массы судна при приеме (снятии) и перемещении грузов При решении задачи, связанной с определением координат центра массы судна при различных вариантах нагрузки, можно использовать любую из приведенных программ. Поскольку в процессе производствен- ной деятельности не все статьи нагрузки подвергаются изменению, ре- комендуется следующее: за исходные координаты центра массы и массовое водоизмещение брать результаты, полученные при предыдущем расчете; предыдущие значения масс статей нагрузок, подвергшихся измене- нию (например, топлива, воды), вписать в таблицу со знаками минус, а их координаты без изменения знаков; текущие значения масс указанньк статей и их координаты записать в таблицу по фактическим данным; произвести расчеты обычным путем, вводя исходные данные с их знаками. Этим обеспечивается универсальность решения задачи. Отметим, что подобную задачу приходится решать при балластиров- ке судна, комбинируя различные варианты, обеспечивающие и остойчи- вость, и надлежащую посадку судна. Особенно удобно для этой цели использовать микрокалькулятор „Электроника МК-52”, который обес- печивает сохранность программы при отключенном питании. 3.3. Расчет плеч статической остойчивости судна Плечо статической остойчивости судна нормируется Правилами Регистра СССР. Рассмотрим задачу расчета плеча применительно к исполь- зованию микрокалькулятора. В качестве исходной выберем следующую формулу [6]: /в zGsine, (3.4) где 1д — плечо статической остойчивости, м; Ir — плечо остойчивости фор- мы относительно основной, м; zG — аппликата центра массы, м; 0 — угол крена, град. Исходную формулу представим в следующем виде: 1в =1к~ zG sin (О,,- ! + 80), (3.5) где 0п _ ! — предыдущее значение угла крена, град; 80 - приращение угла крена, град. 23
Ниже представлены табл. 3.3 и Программа *> пыЧМСЛвНМ* ГШ»Ч стати- ческой остойчивости судна. Таблица 3.3. Размещение исходных данных при расчете статической остойчивости судна Параметр Обозначение регистра Плечо остойчивости формы ‘К X Аппликата центра массы судна 0 Приращение угла крена 60 । Угол крена Плечо статической остойчивости судна 'в X Программа 5 00 01 02 03 04 05 06 07 08 П-»Х0 П->Х1 П-»Х 2 + Х-П2 F7 X - С/П 60 61 62 10 42 1С 12 11 50 Инструкция для работы с Программой 5 1. Включить микрокалькулятор, переключатель поставить в положение Г (градусы) 2. Перейти в режим „Программирование" 3. Занести программу 4, Перейти в режим „Автоматическая работа” 5. Занести константы zG и 60 в регистры 0 и 1 соответственно 6. Произвести очистку регистра 2 (в целях отра- ботки навыка выполнения этой операции в необходи- мых случаях) 7. Занести в регистр X (на индикатор) плечо ос- тойчивости формы для в = 10° 8. Произвести первый цикл вычислений 9. После завершения цикла вычислений снять с индикатора значение плеча статической остойчивости для 0 = 10° 10. Произвести второй и последующий циклы вы- числений, выполнив п. 7 ... 9 Инструкции Нажимаемые клавиши F ПРГ F АВТ zGX->n0 60 Х -П 1 СХХ-*П2 В/0С/П Так как программой предусмотрен автоматический способ вычисле- ний угла крена с приращением 80, то повторное нажатие клавиш В/0 С/П приводит к увеличению угла 0, Следовательно, если нажатие произведено до ввода очередного значения /^, то получится ошибочный результат. В целях проверки правильности результата (значение 1g) для произ- вольного угла в необходимо занести в регистр 2 угол, ранный 0 - 80, в регистр X - 1К, соответствующее углу 0, и произвести пуск программы нажатием клавиш В/0 С/П. После завершения цикла вычислений на инди- каторе высветится искомый результат. 24
Для расчета плеч статической остойчивости при данной аппликате центра массы судна, но при других значениях объемного водоизмещения необходимо выполнить п. 6 ... 10 Инструкции. Если же аппликата также будет другой, необходимо занести ее значение в регистр 0. Пример. Дано: судно типа БАТ „Спрут”, Плечи остойчивости формы даны отно- сительно основной, учтен коэффициент выступающих частей; объемное водоизме- щение 57 = 5500 м3; аппликата центра массы судна z = 7,12 м; углы крена 10 ... 80°. Требуется рассчитать плечи статической остойчивости. Решение, Будем полагать, что программа набрана и отлажена, и микрокаль- кулятор находится в режиме „Автоматическая работа”. Поэтому переходим к вы- полнению п. 5 Инструкции. Нажимаемые клавиши 7.12Х-П0 10 Х-П1 Индикация СХ Х-П2 1 . 34 В/0 С/П 2.72 В/0 С/П 4 . 14 В/0 С/П 5.46 В/0 С/П 6.5 В/0 С/П 7.22 В/0 С/П В/0 С/П 7.6 В/0 С/П 1.34 1.036 - 01 2.72 2.848 - 01 4.14 7.6 9.094 - 01 7.6 5.881 - 01 Плечи статической остойчивости БАТ „Спрут” при V= 5500 м3, zg = = 7,12 м имеют следующие значения: е°.................... 10 20 30 40 50 60 70 80 1К. м................ 1,34 2,72 4,14 5,46 6,50 7,22 7,60 7,60 /fl.M................ 0,10 0,28 0,58 0,88 1,04 1,05 0,91 0,59 3.4. Расчеты, связанные с построением универсальной диаграммы статической остойчивости Регистр СССР рекомендует иметь в информации об остойчивости судна универсальную диаграмму статической остойчивости для определе- ния значений нормируемых параметров статической остойчивости в кри- тической ситуации, т, е. когда нагрузка судна необычна, а остойчивость на пределе. „Необычность” нагрузки по Регистру следует понимать в том смысле, что критическая (допустимая) аппликата центра массы судна в данной ситуации становится недействительной, например, при сущест- венном увеличении площади парусности вследствие принятого палубного груза, или учете массы льда при обледении судна по нормам, отличаю- щимся от принятых Регистром СССР, и т. д. 25
Достоинством универсальной диаграммы статической •>. ошчивости (УДСО) является то, что она позволяет определить значения параметров статической остойчивости графическим способом по следунипнм аргу- ментам: аппликата центра массы судна (или начальная поперечная мета- центрическая высота судна), водоизмещение (массовое или объемное) и угол крена. В настоящее время наиболее употребительными являются два варианта УДСО: ср шкалой аппликат центра массы судна z <;; со шка- лой начальных поперечных метацентрических высот. В связи с этим ограничимся рассмотрением упомянутых вариантов диаграммы. Прежде всего отметим, что особенностью любо'го варианта УДСО является синусная шкала углов крена, т. е. шкала, точка которой рассчи- тана в соответствии с уравнением x=xosinft (3.6) где х — абсцисса точки синусной шкалы, мм; х0 - база шкалы углов крена, мм; в - угол крена, град. Ниже представлена программа для расчета абсцисс точек синусной шкалы УДСО. Программа 6 О 01 02 03 04 05 06 07 08 09 П-Х0 П-Х1 П-Х2 + Х-П2 7 X С/П Б/П 00 60 61 62 10 42 1 12 50 51 00 Краткая инструкция по работе с программой. После набора програм- мы и перехода в режим „Автоматическая работа” необходимо записать: х0 - в регистр 0; 86 - в регистр 1; 0 - в регистр 2. После этого произ- вести первый цикл вычислений нажатием клавиш В/0 С/П. Второй и последующий циклы произвести нажатием клавиши С/П. Пример. Дано: х0 = 200 мм; 6в = 10°; 0=0- 90°. Требуется вычислить значения а для указанного диапазона углов. Решение. Будем полагать, что микрокалькулятор готов к выполнению первого цикла вычислений (переключатель в положении Г). Нажимаемые клавиши Индикация В/0 С/П 34.73 (мм) С/П 68.40 С/П 100. С/П 12836 С/П 153.21 С/П 173.20 С/П 187.94 С/П 196.96 С/П 200. Ниже представлены значения абсцисс синусной шкалы углов крена для соответствующих углов крена;. 0°........ 10 20 30 40 50 60 70 80 90 х,мм..., 34,7 68,4 100 128.6 153.2 173.2 187.9 197 200 Теперь рассмотрим каждый вариант УДСО отдельно. 3.4.1. УДСО со шкалой аппликат центра масс судла. Схематически 26
этот вариант диаграммы выглядит так, как показано на рис. 2. Здесь в - синусная шкала углов крена; 1Р — шкала плеч остойчивости формы относительно произвольно- го полюса; zq — шкала аппликат центра массы судна. Отметим, что шкалы 1р и zq имеют одинаковый масштаб, а точке пересения шкалы 0 со шкалой zg соответствует ап- пликата zg = zp< где zp — аппли- ката произвольного полюса. В ос- нове этого варианта УДСО формула 1о=1к~ (ZG - Z?) sin0, Рис. 2, Универсальная диаграмма стати- ческой остойчивости (3-7) где 1g — плечо статической остойчивости судна; 1К — плечо остойчивости формы относительно основной; zp — аппликата произвольного полюса. Эта формула формально может быть получена из -исходной (3.4) следующей комбинацией: 1в = (lK - zPsinty—(zQ — zPsin0), где lK- Ip sinfl = 1Р. (3.8) Переход от исходной к формуле (3.7) необходим для того, чтобы исключить нерабочую часть диаграммы, а, следовательно, сделать ее ком- пактной. При выборе zp рекомендуется руководствоваться следующим правилом: исходя из типовых случаев расчета zG, имеющихся в информа- ции об остойчивости судна, найти среднее значение, округленное до деся- тых долей метра, и принять его за zp. Для этой цели можно также исполь- зовать среднее значение допустимой аппликаты центра массы судна, найденное из информации об остойчивости судна (из графиков или таб- лиц). После этого необходимо произвести пересчет плеч остойчивости формы, т. е. перейти от 1К к 1Р. Поскольку формула пересчета аналогична формуле расчета плеч статической остойчивости, для этой цели можно использовать рассмотренную ранее программу. Ниже приведены значения 1к и 1р БАТ „Спрут” при zp = 7,0 м для различных значений объемного водоизмещения. О 10 20 30 40 50 60 70 80 Х7,м3 0 1,42 2,86 4,26 5,46 6,48 7,24 7,62 7,62 4500 0 0,20 0,47 0,76 0,96 1,12 1,18 1,04 0,73 0 1,34 2,72 4,14 5,46 6,50 7,22 7,60 7,60 5500 0 0,12 0,33 0,64 0,96 1,14 1,16 1,02 0,71 0 1,32 2,66 4,08 5,44 6,48 7,18 7,52 7,56 6500 0 0,10 0,27 0,58 0,94 1,12 1,12 0,94 0,67 0 1,30 2,64 4,06 5,40 6,44 7,10 7,44 7,49 7500 0 0,08 0,25 0,56 0,90 1,08 1,04 0,86 0,60 27
Сопоставительный анализ показывает, что после пересчета плечи уменьшились в 5-6 раз. Как указывалось, это позволяет выбрать при тех же размерах диаграммы больший масштаб. После выполнения процеду- ры пересчета плеч переходят к построению диаграммы (см. рис. 2). Снача- ла, пользуясь вычисленными значениями абсцисс, разбивают синусную шкалу. Затем, восстановив перпендикуляр в точке, соответствующей нулю синусной шкалы, разбивают (оцифровывают) шкалу 1р. После этого выбранные для заданного водоизмещения и различных углов крена значения 1р наносят на диаграмму, соединяют плавной кривой, надписы- вают значение водоизмещения. После нанесения всех кривых 1р произ- водят привязку шкалы zq следующим образом. В точке 0 = 90 восста- навливают перпендикуляр к синусной шкале углов крена и разбивают его в масштабе шкалы 1р. Поскольку в точке пересечения шкал zg и Ozq — zp = 0, точке пересечения двух шкал соответствует аппликата zG = Zp, что и фиксируется на шкале zq. В нашем примере zp = 7,0 м. Отметим, что вниз от шкалы в zq < Zp, а вверх — Zq > zp. Как видно из описания, построение УДСО со шкалой аппликат центра масс не пред- ставляет трудности. Примечания. 1. Определение 1д с помощью УДСО заключается в следующем (рис. 3): на шкале zg откладывают аппликату центра массы судна, вычисленную обычным способом; полученную точку а соеди- Zc няют прямой с началом синусной __________________шкалы; в точке b синусной шкалы, соответствующей заданному углу Zp крена 0, восстанавливают перпендику- Рис. 3. Схема определения плеча ста- тической остойчивости ляр и продолжают его до пересечения с кривой плеч 1р (для заданного водоизмещения) и вспомогательной прямой Оа\ снимают отрезок cd перпендикуляра от вспомогательной прямой (она называется прямой плеч остойчивости массы) до кривой /Р; с помощью шкалы 1р переводят в метры, получают плечо статической остойчивости 1д. 2. Если участок кривой плеч IР окажется ниже прямой плеч остой- чивости массы, это будет означать, что имеет место отрицательная остой- чивость и, следовательно, судно получит статический крен, угол которого соответствует точке пересечения прямой и кривой плеч. Действительно, в этом случае 10 = 0, тогда как левее этого угла - 10 < 0. 3. Определение максимального плеча статической остойчивости, угла заката, метацентрической высоты и их регулирования являются частными вариантами первого примера, поэтому здесь не рассматриваются. Однако отметим, что аппликату начального поперечного метацентра определяют проведением касательной к начальному участку кривой /г и продолже- нием его до пересечения со шкалой zq. В точке пересечения z(- = zm. Это легко поясняется следующим примером. МетаиеН1риче< кая формула 28
остойчивости, как известно, имеет вид lg = /zsin0. Так как в случае касательной 1в = 0 при в & 0, то, следовательно, h = 0. Получаем, что 2Ш = ZG- Таким образом, с помощью УДСО вследствие простоты способа решения задач можно (дополнительно) успешно проводи!., техническую учебу. 3.4.2. Универсальная диаграмма статической остойчивости со шка- лой h. Этот вариант диаграммы базируется на формуле 1д = 1м + Л sinft (3.9) где Iм - плечо остойчивости формы относительно начального поперечного метацентра; h — начальная поперечная метацентрическая высота; 0 — угол крена. Для пересчета плеч остойчивости формы служит формула 1м= 1 *к ~ zmsinO, (3.10) где 1К - плечо остойчивости формы относительно основной; zm - аппли- ката метацентра (водоизмещение, м3), выбираемая из гидростатических кривых по аргументу. Ниже приведены исходные данные значений плеч остойчивбсти формы относительно начального поперечного метацентра Б АТ „Спрут” для объем- ного водоизмещения \7 = 4500 м3 и zm = 8,18 м и результаты. е°.............. о 10 1к. м........... 0 1,42 /дт. м.......... 0 0,00 20 30 40 50 2,86 4,26 5,46 6,48 0,07 0,17 0,21 0,22 Примечания. 1. Перерасчет 60 70 80 7,24 7,62 7,62 0,16 -0,05 -0,43 плеч для другого водоизмещения производится аналогично. При этом zm выбирается, естественно, для этого водоизмещения. 2. Синусная шкала углов крена может быть выбрана такого же раз- мера, как и в первом случае построения УДСО. 3. Нанесение кривых плеч остойчивости формы производится анало- гично. 4. Привязка шкалы h к диаграмме производится на основании прави- ла: в точке пересечения синусной шкалы углов крена (О = 90°) значение h равно нулю; вверх значения h имеют знак минус, вниз - плюс. 5. Для лучшей ориентации диаграммы в некоторых случаях произ- водят пересчет плеч по формуле lMi =1м+Ь sine, (З.П) где lMl - новое значение плеча остойчивости формы; b - отрезок, на который точка нуль шкалы h поднимается вверх относительно первона- чального положения. Так, если принять b = 0,6 м, то получим соответст- венно следующие значения: 29
е°............ О 10 20 30 40 50 60 70 80 1м,-м......... О 0,10 0,28 0,47 0,60 0,68 0.68 0.51 0.16 Следует отметить, что УДСО со шкалой zG имеет некоторые преиму- щества перед диаграммой со шкалой h с точки зрения как построения (zp = const), так и использования аргумента zG. Однако суда ГДР, напри- мер, имеют второй вариант УДСО. 3.5. Расчет плеч динамической остойчивости Диаграмма динамической остойчивости представляет собой график работы восстанавливающего момента или плеч динамической остойчи- вости в зависимости от угла крена. При вычислении, работы восстанавли- вающего момента используют формулу Ав = 0,560 (Mo + 2Af, + ... +2М„_ , +Мп), (3.12) где 50 — интервал угла на ДСО между смежными моментами, выражен- ный в радианной мере; Мо, .., Мп - восстанавливающие моменты для соответствующих углов крена. Разделив обе части уравнения (3.12) на Д, получим /^ = Дв/Д= 0,560 (/0 + 2/1 + ... 2/„_ ,+/„), (3.13) где ld — плечо динамической остойчивости, м. Если принять 60 = 10° и выражение в скобках обозначить через S (интегральная сумма), то окончательно получим ld = 0,0873 2/. (3-14) Сначала рассмотрим способ расчета плеч динамической остойчивости судна с помощью калькулятора в непрограммируемом варианте. Предположим, что уже найдены с помощью УДСО. Дано: судно типа БМРТ пр. 394; Д = 3100 т (30411кН); Vo =3000 м3; zG = 5,67 м; О = 0 ... 90°; 60 = 10°. Расчет ведем с использованием микрокалькулятора „Электроника Мк-56”, результаты которого приведены ниже: Нажимаемые клавиши 1, Включить МК 2, 5 В t 57,3 х->по 3. 0,09 В t 5. 0,23 + 6. F Вх + 0,47 + 7. F Вх + 0,78 + Индикация 0. 8.726 ... -02 1.8 4,1 2J6 -02 -01 0,1 30
Нажимаемые клавиши 8. F Вх + 1,12 + 9. F Вх + 1,40 + 10. F Вх+ 1,43 + 11. FBX+1.35 + 12. FBx + l,13 + 4,26 6.78 9.61 12,39 14,87 После расчета S' переходам к вычислению ld. Нажимаемые клавиши 1.0,09 П-ХОХ 2.0,41 П-ХОХ 3.0,11 П-ХОХ 4, 2,36 П-ХОХ Индикация 7,853 ... -03 3,57 ... -02 9,68 ... -02 2,059 -01 Аналогично выполняют п. 5 ... 10. 10. 14,87 П-ХОХ 1,297 Значения плеч ниже: статической и динамической остойчивости судна даны 0°............ 0 10 20 30 40 50 60 70 80 90 'в.м......... 0 0,09 0,23 0,47 0,78 1,12 1,40 1,43 1,35 1,13 S.M.......... 0 0,09 0,41 1,11 2,36 4,26 6,78 9,61 12,39 14,87 ....... О 0,01 0,04 0,10 0,21 0,37 0,59 0,84 1,08 1,30 По данным в и ld (см. Программу 3) построена диаграмма динами- ческой остойчивости (ДНО). Теперь рассмотрим способ определения динамического угла крена по ДЦО (рис. 4). Динамический угол крена, как известно, определяется исходя из равенства Ад — Ad — Md0d, (3.15) Разделив обе части уравнения (3.15) на Д, получим 31
57,3° Рис. 5. Определение динамического угла крена (а) и1с (6) Нетрудно заметить, что график (/dKp) представляет собой прямую, проходящую через начало координат. Следовательно, чтобы нанести ее на поле диаграммы, необходимо найти координаты второй точки. Для этой цели надо, например, задаться значением угла и вычислить/dKp. Если принять 0d = 1 рад, то получим /ЛР-МЛ/Д. (3.17) Следовательно, координатами второй точки будут 6 = 1 рад; ldKp = = Md^. Таким образом, определение динамического угла крена сводится к следующему. На шкале углов крена (рис. 5, а) откладывают вправо от нуля угол, равный одному радиану (точка а). В полученной точке а восстанавлива- ют перпендикуляр, а по нему измерителем откладывают отрезок ab, равный в масштабе ldKp для 0d = 1 рад. Через точки 0 и b проводят отрезок прямой. Из точки пересечения прямой с ДДО с опускают перпендикуляр на шкалу углов крена 00. В точке d снимают значение динамического угла крена. Пример. Дано: Д= 3100 т (30411 кН); zq = 5,67 м; Md = 592,0 т м (5807,5 кН • м), Требуется определить 0d, Решение. Находим ldKp = 592/3100 = 0,19 м. Произведя необходимые построения на ДДО, получим 9d = 30°. Определение момента Md по динамическому углу крена производится в обратном порядке. Рассмотрим табл. 3, 4 и Программу 6 расчета плеч статической и ди- намической остойчивости, 32
Таблица 3.4. Размещение исходных данных и результата при расчете плеч статической и динамической остойчивости Параметр Обозначение Символ регистра Плечо остойчивости формы ‘К X Аппликата центра массы судна ZG 0 Приращение угла крена 60 1 Угол крена в 2 Половина приращения угла крена в радианной 6в /2 3 мере Плечо статической остойчивости ‘в 4 Плечо динамической остойчивости ‘d X Программа 6 00 01 02 03 04 05 06 07 08 09 П->Х 1 2 4- F 9 Х-ПЗ F П-Х 1 П-Х2 + Х-П2 61 02 13 1 Е 43 25 61 62 10 42 10 12 13 14 15 16 17 18 19 F 7 П^ХО X — П-Х4 -• Х-»П4 + П->Х5 + 1 С 60 12 64 14 44 10 65 10 20 21 22 23 24 25 Х-П5 П-»Х 3 X С/П БП 06 45 63 12 50 51 06 Краткая инструкция по работе с программой. После набора програм- мы и перехода в режим „Автоматическая работа” выполняют следующие операции: заносят константы zq и 80 в регистры 0 и 1 соответственно; производят очистку регистров 2, 4 и 5 (хотя после включения микро- калькулятора эти регистры памяти нет необходимости очищать, все же эту операцию необходимо выполнить в целях приобретения навыка); при повторном решении примера или решении другого примера очистка ре- гистров 2, 4 и 5 обязательна; заносят в регистр X (на индикатор) значение плеча остойчивости формы для угла крена. 10°; производят первый цикл вычислений путем нажатия клавиш В/0 С/П; после завершения цикла вычислений с индикатора снимают значение плеча динамической остойчивости ld, а из регистра 4 извлекают значение плеча статической остойчивости для 0= 10°. Второй и последующий циклы вычислений производят нажатием клавиши С/П (в программе предусмотрена команда пуска БП 06) после ввода на индикатор 1К. Примечание. Если по каким-либо причинам не удалось снять значение ld, то рекомендуется повторить вычисление части программы: произвести очистку регистра X нажатием клавиши Сх; нажать клавиши БП 18 С/П. После завершения цикла вычислений на индикаторе высветится значение Zd. При этом содержимое регистров 2,4, 5 не изменится. 33 664 —3
При решении второго примера, как указано выше, необходимо про- извести очистку регистров 2,4, 5. Для отработки навыка по использованию программы рекомендуется использовать данные плеч статической и динамической остойчивости судна, приведенные ранее. 3.6. Расчеты, связанные со снятием судна с мели (камня) При посадке .судна на мель уменьшается осадка судна, вследствие чего возникает реакция грунта, равная по массе потерянной плавучести (части массового водоизмещения): 2V= MX)(dab-d)q, (3-18) где N — реакция грунта, кН; dab — средняя осадка судна после посадки на мель, м; d — средняя осадка судна до посадки на мель; q — число килоньютонов на один сантиметр осадки при осадке, равной dab. Если судно получило пробоину, то реакция грунта увеличивается на значение массы влившейся воды. При посадке судна на риф (камень) и сохранении герметичности корпуса реакцию грунта определяют по формуле jV= - pS8d, грр N- реакция грунта, кН; bd - изменение средней осадки в район центра тяжести площади действующей ватерлинии, м; S — площадь действующей ватерлинии, м2; р — плотность воды, кг/м3. Значение тягового усилия, необходимого для снятия судна с мели, определяют по формуле F = fN, (3.19) где N — реакция грунта, кН; f — коэффициент трения скольжения корпуса судна о грунт. Коэффициент трения скольжения о грунт зависит от материала корпу- са судна, степени его шероховатости, наличия деформаций и структуры грунта. Ниже приведены ориентировочные данные в зависимости от упомянутых параметров: Коэффициент трения Характер грунта Жидкая глина................. Мягкая глина................. Глина с песком............... Песок........................ Галька....................... Каменная плита............... Камень....................... Корпус деревянный 0,15 0,20 0,30 0,35 0,38 0,60 0,70 Корпус стальной 0,20 0,30 0,35 0,40 0,42 0,35 ... 0,40 0,40 ... 0,50 34
Иногда перед снятием судна с мели производят расчеты по уменьше- нию реакции грунта путем предполагаемого снятия или продольного пере- мещения груза. Далее следует иметь в виду изменение глубины в районе посадки судна вследствие прилива и отлива моря. При увеличении глуби- ны акватории реакция грунта будет уменьшаться, при уменьшении, наобо- рот, — возрастать. Пример 1. Судно, имеющее массовое водоизмещение Д = 36 297 кН (3700 т), село на мель. Расчеты после замера осадок носом и кормой позволили получить значения: Sd = -0,3 м; q = 87,4 кН/см (8,9 т/см), Требуется рассчитать.реакцию грунта. Решение. № -100 • 87 4 • 0,3 =-2620 кН (267т-с). Пример 2. Рассчитать значение тягового усилия, необходимого для снятия судна, с мели по данным предыдущего примера, если грунт представляет собой глину с песком, а корпус судна стальной. Решение. F = /№ 0,35 • 2620 кН = 917 кН (93,45 • с). Пример 3. Определить уменьшение метацентрической высоты судна, используя данные примера 1. Решение. Уменьшение метацентрической высоты судна h вычисляют по формуле где 6й - изменение метацентрической высоты, м; d осадка судна до посадки на мель, м; h - метацентрическая высота судна до посадки судна, м; /V - реакция грунта, кН; Д - массовое водоизмещение судна до посадки на мель, кН. Пример 4. Дано: Д = 36 297 кН (3700 т); d = 5,75 м, h = 0,25 м, 6d = -0,3 м; N = -2620 кН (267 т). Требуется рассчитать 6/1. Решение. Подставив в формулу (3.20) исходные данные, получим -2620 6Л --------------- 36297 - 2620 = -0,08 м. Следовательно, метацентрическая высота уменьшится на 8 см. 3.7. Расчет опрокидывающего момента по методике Регистра СССР Опрокидывающий момент Мс является максимально возможным кренящим моментом, приложенным динамически, который с учетом энергии качки наклоняет судно,-не опрокидывая, а при наличии угла заливания наклоняет судно на угол, меньший угла заливания Ситуация. Судно находится лагом к волне, без хода, с выключенны- ми успокоителями качки (если таковые имеются). В момент наклонения 35
судна наветренным бортом на угол, равный вг (амплитуда качки с обеспе- ченностью около 2%, т. е, максимальная из 50 колебаний), на него нале- тает шквал. Эта ситуация для судна является наихудшей, так как и кренящий и восстанавливающий моменты будут действовать в одном направлении. Определение Мс при отсутствии угла заливания. На ДДО (рис. 5, б) вправо от нуля откладывают угол, равный вг (точка а'). В полученной точке восстанавливают перпендикуляр к шкале углов крена и продол- жают его до пересечения с ДЦО (точка А1). Через полученную точку А‘ проводят прямую, параллельную шкале углов крена, а по ней влево от нуля откладывают угол, равный вг, полу- чают так называемую исходную точку А. Из этой точки проводят каса- тельную АС к ДДО, которая представляет собой график плеч опрокиды- вающего момента для данной ситуации. Далее по прямой, проходящей через точки АА', откладывают в масш- табе отрезок АВ, равный 1 рад (57,3°). В точке В восстанавливают пер- пендикуляр и продолжают его до пересечения с „асательнойЛС (точка £). Снимают отрезок BE, переводят с помощью масштаба плеч в метры, получают 1С. Опрокидывающий момент вычисляют по формуле Мс = Ыс. (3.21) Определение опрокидывающего момента при наличии угла заливания. В этом случае вместо касательной АС производится секущая AF, так как диаграмма обрывается на угле заливания В остальном задача решается аналогично (см. рис. 5, б) Рекомендуется использовать аналитический способ расчета 1С. Д1я пояснения обратимся к рис. 5, б, в частности, к прямоугольным подобным треугольникам AFf и АЕ1В, произведем элементарные выкладки, после чего получим . , . 57,3° ВЬ‘ -l„~ (1„ - а'Л') . (3.22) иг + иг где ld - плечо динамической остойчивости судна, соответствующее углу заливания fy-; а'Аг — плечо динамической остойчивости, соответствую- щее амплитуде качки вг. Пример. Дано; Д = 3100 т; вг = 20°; 0р= 50. Требуется рассчитать 1С. Решение. Соответственно находим: а'А*' — 0,04 м; ld = 0,37 м. Подставив исходные данные в формулу, получим: 57,3° 57,3° 1С = (0,37 - 0,04) ---- = 0,33 ----- = 0,27 м. 70° 70° Ответ: 1С = 0,27 м. 36
3.8. Расчет поперечной метацентрической высоты по измеренному периоду собственных бортовых колебаний судна Сущность способа, рекомендованного ИМКО (ныне ИМО) в 1965 г. для судов длиной до 70 м, заключается в следующем. Соответствующим способом на тихой воде или на волнении измеряется серия (не менее 5) периодов собственных бортовых колебаний судна, затем вычисляется среднее значение. После этого вычисляют h по формуле h = (сВ/тв)2, где h — начальная поперечная метацентрическая высота, м; с - инерцион- ный коэффициент, с • м-0’5 (определяется при креновании судна); В - ширина судна, м; тд — среднее значение измеренного периода собствен- ных бортовых колебаний судна, с. Пример. Дано: В = 9,8 м; с = 0,79 с-м-1/2; тд = 10,3 с. Требуется опреде- лить h. Решение, h = (0,79 • 9,8/10,3)3 = 0,56 м. В целях упрощения расчетов используют различного рода графики, в том числе номограмму, предложенную ИМКО- Удобными для исполь- зования следует считать графики, построенные на основе формулы (3.22) . Для расчета точек кривой h = /(те) при фиксированных значениях Вкс рекомендуется использовать программируемый микрокалькулятор. Ниже приведены Программа 7 расчета метацентрической высоты по из- меренному периоду собственных бортовых колебаний судна и пример этого расчета. Программа 7 00 01 02 03 04 05 06 07 08 П-Х1 П'->Х 2 + Х-»П2 П-»Х0 -> + FX С/П 61 62 10 42 60 14 13 22 50 Распределение исходных данных и результата по соответствующим регистрам памяти приведено ниже: Параметры.............с • В бтд то Регистры................ 0 1 2 Примечание, бтд — приращение периода тд Пример. Дано: судно типа БМРТ; В = 14,0 м; с = 0,76 с • м1/2; тд =10 4-30 с; St о = 0,5 с. Требуется вычислить значения Л в заданном интервале тд. Решение представлено в виде последовательных операций. Нажимаемые клавиши 1. Включить МК 2. Перевести МК в режим „Программиро- F ПРГ (В/0 F ПРГ) вание” 37
3. Записать программу 4. Перевести МК в режим „Автоматичес- кая работа” 5. Вычислить произведение чисел № С X X В и записать в регистр О 6. Занести приращение периода (0,5 с) в регистр 1 7. Число 9.5 занести в регистр 2 8. Произвести необходимое количество циклов вычислений Программа 7, строка 2 F АВТ Х-»П0 0.5Х->П 1 9.5 Х-*-П2 В/0 С/П Ниже приведен фрагмент результатов вычислений т0...........10,0 10,5 11,0 11,5 12,0 12,5 13,0 13,5 14,0 ... 30 Л, м........ 1,13 1,02 0,94 0,86 0,79 0,72 0,67 0,62 0,58...0,13 Примечания, 1. При расчете метацентрических высот для друго- го значения с необходимо предварительно выполнить действия, изложен- ные в пп. 5 и 7 примера, а затем перейти к п. 8. 2, При определении метацентрической высоты на волнении следует выбирать такие эпизоды, когда собственные колебания судна существен- но преобладают над вынужденными. 3. В целях повышения точности измерения периода собственных борто- вых колебаний необходим индикатор вертикали на базе простейшего гиромаятника с бесконтактным способом записи колебаний на бумаге (например, электроискровым). Исследования автора показали, что при наличии регистратора качки можно просто определить период собствен- ных колебаний судна, 4. Так как инерционный коэффициент зависит как от момента инер- ции массы судна, так и от массового водоизмещения, то необходимо иметь значения с для некоторых (характерных) значений водоизмещения с учетом способа размещения статей нагрузки относительно диаметраль- ной плоскости. 3.9. Построение универсальной диаграммы динамической остойчивости Как указывалось, Регистр СССР рекомендует иметь в „Информации об остойчивости судна” или в приложении к ней справочные материалы „для более точного определения всех показателей остойчивости и посадки судна при необычных случаях нагрузки судна, когда остойчивость его может оказаться на пределе”. К таким материалам с полным основанием можно отнести универсальную диаграмму динамической остойчийости (УДДО). Как известно, одним из нормируемых параметров является критерий погоды К = Мс/М„, (3.23) 38
где Мс — опрокидывающий момент, определяемый по методике Регистра СССР; Mv — момент шквала, определяемый также по методике Регистра СССР. Отметим, что Мс можно определить как по диаграмме статической остойчивости (глазомерной оценкой равенства двух площадей, одна из которых имеет, как правило, сложную конфигурацию), так и с помощью диаграммы ДНО. Второй способ является более точным, так как исклю- чаются ошибки, возникающие от неточности глазомерной оценки равенст- ва площадей, которые имеют место в первом способе. Однако он гораздо сложнее, поскольку в классическом варианте, как мы убедились, требу- ется сначала найти плечи статической остойчивости, потом рассчитать плечи динамической остойчивости, При использовании УДЦО плечи дина- мической остойчивости находятся графическим способом по аргументам: водоизмещение (весовое или объемное); аппликата центра величины для прямого положения судна zc0; аппликата центра тяжести судна для данного случая нагрузки zG; углы крена в/. УДЦО строится на основе формулы [6] ~ - zc0)cos9 - (1 -cos9) (zG - zc0), (3.24) где Ус, z с — координаты центра величины, соответствующие углу крена в и фиксированному значению угла дифферента. Уравнение (3.24) представим в следующем виде: ld = ldl~^d2, (3.25) тогда Idi =J>csin9 - (zc- zc0)cos9 - (1 - cos 9)zP, (3.26) ld2 = (1 - cos 9) (zG -Zp- zc(i), где Zp — число, соответствующее точке пересечения шкалы (zG — zc0) со шкалой углов крена и служащее для ориентировки диаграммы (рис. 6). Zs-Za, Рис, 6. Схема универсальной диаг- Рис. 7. Универсальная диаграмма ди- раммы динамической остойчиво- намической остойчивости 39
Из уравнения (3.26) видно, что значения ldt можно вычислить еще на стадии проектирования судна, так же как и плечи остойчивости формы, и представить их в виде кривых/^ =/( Vo, 0), где Vo — теоретическое объемное водоизмещение. Что касается Zj2, то для его нахождения графи- ческим способом необходимо, по аналогии с универсальной диаграммой статической остойчивости, иметь две шкалы: горизонтальную, пропорцио- нальную значениям 1 - cos 0, вертикальную, проходящую через точку 0 = 90° для значений zG — zp — zcQ. Отметим, что zp имеет постоянное значение и служит для устранения нерабочей части, а также ориентировки диаграммы. Поэтому целесообразно в качестве входного аргумента шка- лы принять разность ZG - ze0. (3.28) Для выбора zр рекомендуется использовать уравнение Zp — — Zco)cp> где zGср - среднее значение допустимой (критической) аппликаты центра масс судна, взятое из всех возможных значений водоизмещения и условий плавания; zcOcp — среднее значение аппликаты центра величины для пря- мого положения судна, взятое из возможных значений водоизмещения судна. На рис. 6 схематически показаны УДДО и принцип определения ld. Д1Я больщей наглядности приведена только одна кривая плеч динамичес- кой остойчивости. Из изложенного видно, что построение УДДО аналогично построению УДСО. Теперь рассмотрим второй способ построения УДДО [6, с. 35]. Сущ- ность его заключается в том, что сначала вычисляют ld для нескольких значений водоизмещения при одной и той же аппликате центра массы судна, а потом „вмещают” их в каркас диаграммы. В качестве иллюстра- ции приведем пример, ограничившись двумя значениями водоизмещения (рис. 7). Пример. Дано: судно типа БМРТ пр. 394; 70] = 3000 м3; V„ = 3600 м3; ZG~ 5,67 м. Требуется построить две ветви УДДО. Решение. 1. Принимаем х0 = 200 мм и по формуле х= 200 • (1 - cos в) = 400 • sin1 (0/2) вычисляем х,-; заносим эти значения в таблицу, которая представлена ниже. 2. С помощью аргументов \7„ = 3000 м3; V, = 3600 м3 и zG = 5,67 м на- ходим сначала плечи статической остойчивости, а затем плечи динамической остой- чивости и заносим их в третью и четвертую строки таблицы соответственно. 0°.........0 10 20 х, мм. ... 0 3,0 12,1 ld,M ... . 0 0,01 0,05 Id- м . . . . 0 0,01 0,04 30 40 50 60 26,8 46,8 71,4 100 0,12 0,23 0,40 0,61 0,10 0,21 0,37 0,59 70 80 90 V„,m3 131,6 165,3 200 0,84 1,06 1,25 3600 0,84 1,08 2,30 3000 40
3. По аргументам Vo = 3000 м’ центра величины: zCo = 2,7i м; zCOi = 4. Находим разности: zq - zCo = -3,15 = 2,52 м. 5. Принимая zc„c = 2,8 м, a ZGcp : 6. По данным таблицы строим шкалу углов крена диаграммы. 7. Через точки 0 и 90 шкалы углов крена проводим две вертикальные шкалы, одна из которых (левая) представляет собой шкалу , а вторая - шкалу zq — — zp - zco. Общие шкалы имеют одинаковый масштаб. При этом левая шкала оцифровывается от нуля: положительные значения - вверх, отрицательные — вниз от шкалы углов крена. Правая шкала оцифровывается следующим образом. Точка пересечения со шка- лой углов крена соответствует значению zp. Вверх от шкалы 0 числа будут возрас- тать, вниз - уменьшаться. При этом шкала оцифровывается в значениях zq — zCo. На этом построение каркаса диаграммы заканчивается. 8. Откладываем разность zq — zco = 2,95 м по правой шкале и проводим прямую Оя (см. рис. 6). От прямой 0а откладываем (в масштабе) вверх, перпенди- кулярно шкале углов крена, значения (отрезок cd) для соответствующих углов крена (третья строка таблицы). Полученные точки соединяем плавной кривой, надписываем значение водоизмещения (см. рис. 7): Vo = 3000 м’. 9. Аналогично производится построение второй кривой для Vo = 3600 м’. 3600 м’ находим аппликаты Использование УДДО рассмотрим на примере. Пример. Дано: судно типа БМРТ пр. 394; Д= 3093 т; KV = 1006; р = = 1,025 т/м’; zq = 5,97 м; 0 = 0* 60° . Требуется найти 1^. Решение. Находим теоретическое объемное водоизмещение 3093 V» = -------------= 3000 м’. 1,006 -1-.025 2, По аргументу V» = 3000 м’ находим zCo = 2,72 м. 3. Вычисляем разность zq - zCa = 5,97 - 2,72 = 3,20 м. 4. Откладываема (у— zCo = 3,20 по правой шкале, а затем проводим прямую, соединяющую начало координат с полученной точкой. 5. -С помощью измерителя снимаем lj (Д = 3093 т; zq = 5,97 м) : 10 20 0,02 30 0,06 40 50 60 может быть Примечание, i. Второй способ построения УДДО прост и реализован даже в учебном процессе. 2. Учет влияния дифферента на плечи динамической остойчивости может быть реализован в виде поправок. 3. Наличие УДДО в информации об остойчивости судна и пособиях позволит повысить эффективность технической учебы по наиболее сложному разделу динами- ческой остойчивости судна. Для удобства расчета точек шкалы углов крена приводим программу для микрокалькулятора „Электроника МК-56”. Константы и переменные по регистрам памяти распределятся следую- щим образом: 41
Параметры............... х„ 60 О Регистры.................... О 2 3 Далее дана программа для расчета . Программа 8 Г 01 02 03 04 05 06 07 08 09 10 11 12 2 П-ХО П—X 2 П-»ХЗ + Х-ПЗ F8 П-Х1 1 *> - X С/П 3 60 62 63 10 43 1Г 61 01 14 11 12 50 В ней предусмотрено автоматическое вычисление угла (как и в про- грамме, служащей для пересчета плеч остойчивости формы). Результат вычислений высвечивается на индикаторе. Если принять хо = 150 мм, то получим следующие значения: в".............. 10 20 30 40 50 60 70 80 90 х,мм..._,...... 2,3 9,0 20,1 35,1 53,1 75 98,7 124,0 150 Глава 4 РАСЧЕТ ИНЕРЦИОННО-ТОРМОЗНЫХ ХАРАКТЕРИСТИК СУДНА В УСЛОВИЯХ ПРОМЫСЛА Общие сведения. В правиле 6 МППСС-72 указывается, что „каждое судно должно всегда следовать с безопасной скоростью с тем, чтобы оно могло предпринять надлежащее и эффективное действие для преду- преждения столкновения и могло быть остановлено в пределах расстоя- ния, требуемого при существующих обстоятельствах и условиях”. При выборе безопасной скорости учитывают следующие факторы: маневрен- ные возможности судна и в особенности расстояние, необходимое для полной остановки судна; поворотливость судна в преобладающих усло- виях. Таким образом, знание инерционно-тормозных характеристик край- не необходимо судоводителю при всех обстоятельствах плавания (ветер, волнение, степень обрастаний корпуса, состояние нагрузки и посадки судна), особенно при плавании в условиях большой плотности судопо- тока. Если суда морского флота практически ежемесячно могут произ- водить испытания на мерной миле, то суда рыбопромыслового флота, находящиеся в море длительное время, лишены этой возможности. Следо- вательно, наша задача состоит в том, чтобы дать судоводителю промыс- лового флота метод, который позволил бы ему самостоятельно опреде- лять инерционно-тормозные характеристики судна опытно-расчетным путем в открытом море при различных загрузке и степени обрастания корпуса. 42
Как известно, существуют методы определения маневренных элемен- тов судна, основанные на „абсолютной” привязке судна относительно береговых ориентиров, однако они непригодны для применения в откры- том море. В основном судоводители пользуются данными о маневрен- ных элементах, полученными при заводских испытаниях серийных (голов- ных) судов. Однако, как показывает сравнение их с результатами спе- циальных наблюдений, в большинстве случаев они не более точны, чем результаты судовых наблюдений, выполненных в ходе эксплуатации. В связи с этим необходимость определения инерционно-тормозных харак- теристик в условиях плавания очевидна. В практике управления судном существует два способа торможения судна: активное и пассивное. В первом случае торможение судна производится работой движителя на задний ход, во втором — инерция судна гасится сопротивлением воды движению судна. Активное торможение, в свою очередь, делят на два этапа: 1) предварительный — с момента подачи команды (машинным теле- графом) до момента времени, когда движитель станет работать на задний ход (для упрощения задачи будем полагать, что переходный процесс от- сутствует) . Этот этап, по существу, представляет собой период пассивного торможения; 2) период активного торможения — от момента начала работы винта на задний ход до полной остановки судна. Для упрощения вопроса и сохранения основных закономерностей движения судна, имея в виду цель — предложить судоводителю наиболее простой и удобный математический аппарат для использования непосред- ственно в судовых условиях, предположим, что погода штилевая, судно не имеет ни крена, ни дифферента, движение судна в период торможения прямолинейное; сила упора винта (фиксированного шага) при заданном режиме движения судна на заднем ходу остается практически постоянной. 4.1. Пассивное торможение Исходные уравнения: v = voe<-~K/m)s.; (4.1) (4.2) Приводим их следующему виду: г/1>о =e<~K/m>s-, (4.3) vot = w(e<K/m>s - I)/К, (4.4) 43
где v0 _ начальная скорость судна, м/с;р - мгновенная скорость судна; е — основание натуральных логарифмов; S тормозной путь судна, м; т — масса судна с учетом присоединенной массы воды, кг; К — гидродина- мический коэффициент, Н • м-1. Для учебных целей коэффициент К рекомендуется определять по следующей формуле: K = &l3g\ (4.5) где А - массовое водоизмещение судна, т; g - ускорение свободного падения. Так, для А = 3200 т получим К = 217,1534 9,81 = 2130,3 Н м-1. Определение коэффициента К в натурных условиях сводится к замеру скорости судна при пассивном торможении через интервал времени, когда скорость уменьшается примерно наполовину, и расчету К по формуле . (4.6) t \ V Vo / Здесь К — гидродинамический коэффициент, Н • м; t — интервал времени торможения, соответствующий уменьшению скорости с v0 до V, м/с. Пример. Дано: т = 3520 • 103 кг; v0 =10 уз = 5,14 м/с; v = 5,5 уз = 2,83 м/с; t = 266 с. Требуется рассчитать К. 3520 • 10э / 5,14 - 2,83 \ Решение. К =---------------I -----------) = 2101,5 Н • м"1. 5,14-266 \ 2,83 / Ниже представлены программа для расчета инерционно-тормозных характеристик и распределение исходных данных (табл. 4.1). Таблица 4.1. Размещение исходных данных и результата при расчете инерционно-тормозных характеристик судна при пассивном торможении Параметр Обозначение (значение) Символ регистра Константа для перевода скорости в узлы 1,94 0 Гидродинамический коэффициент К 1 Масса судна с учетом присоединенной массы т 2 воды Приращение тормозного пути SS 3 Текущее значение тормозного пути S 4 Отношение гидродинамического коэффициента К/т 5 к массе судна с учетом просоединенной мае- сы воды Отношение текущей скорости к начальной v/vt ff Произведение начальной скорости на интервал Х/7 времени торможения 44
Программа 9 00 - 01 02 03 04 05 06 07 08 П-Х1 П-Х2 Х-П5 П-Х5 П-ХЗ П^Х4 =• Х-П4 61 62 13 45 65 63 64 10 44 09 10 11 12 13 14 15 16 X /-/ F1 Х-П6 П->Х0 П-»Х5 4 П->Х6 12 0L 16 46 60 65 13 66 17 18 19 20 21 22 23 24 F4 1 - X Х-П7 С/П БП 04 . 23 01 11 12 47 50 51 04 Инструкция по работе с программой 9 1. Включить МК 2. Перейти в режим „Программирование” 3. Занести программу 4. Перейти в режим „Автоматическая работа” 5. Занести исходные данные в Соответствующие Нажимаемые клавиши F ПРГ F АВТ регистры памяти: 1,94 - в регистр 0 К - в регистр 1 m - в регистр 2 85 - в регистр 3 6. Произвести очистку регистра 4 (если до этого МК был в работе) 7. Произвести первый цикл вычислений 8. Занести с индикатора в таблицу полученное значение v/v0 9. Извлечь из регистра 6 на индикатор получен- ное значение v0 г 10. Занести полученное значение в таблицу 11. Произвести второй и последующие циклы вы- числений 1.94Х-П0 КХ->П 1 Х->П2 Х->ПЗ СХХ->П4 В/0 С/П П-Х6 С/П Пример. Дано: Д = 3200 т; m = 3520 • 103 кг; К = 2101,5 Н м"‘; 5 = 0- - 3500 м. Требуется рассчитать v01 и v/ v0. Решение. 1. Принимаем 85 = 500 м. 2. Выполняем операции, указанные в пп. 1 ... 4 Инструкции, регламентирую- щей работу с этой программой. 3. Заносим исходные данные в соответствующие регистры памяти и производим очистку регистра 4. Нажимаемые клавиши 1,94 Х->П0 2101,5 Х-П 1 3520 ВП 03 Х -П 2 500Х-*ПЗ СХ Х-П4 4. Производим первый цикл вычислений. Нажимаемые клавиши В/0 С/П Индикация 1,94 2101,5 3520000 500. 0. Индикация ИЗО. (число округлено до целых секунд) 45
5. Извлекаем из регистра 6 значение v/ v„ для первого цикла вычислений. Нажимаемые клавиши Индикация П—Х6 7,4 ... -0,1 (округление производит- ся до сотых долей) 6. Производим последующие циклы вычислений. Результаты инерционно-тор- мозных характеристик судна при пассивном торможении следующие: 5,м...............О 500 v/v„..............1 0,74 vot, уз • с......0 ИЗО 1000 1500 2000 0,55 0,40 0,30 2654 4707 7475 2500 3000 3500 0,22 0,17 0,12 11206 16234 23010 Полученные данные используются при оценке ситуации. 4.2. Активное торможение Расчетные формулы va = \Jf[k tg=[arctg (y/Kvo!F) ~ 5 т 1 1Г1п (у20 + f/k \ U + ) где га — мгновенное значение скорости движения судна, м/с; р0 — ско- рость движения судна в начальный момент активного торможения, м/с; Та - интервал времени активного торможения, отсчитываемый от началь- ного момента активного торможения, с; 5а — тормозной путь, м; F — сила упора винта на заднем ходу, Н. Расчет инерционно-тормозных характеристик судна при активном торможении представлен в Программе 10. 00 01 02 П-Х2 П-Х0 П-»Х1 62 60 61 10 11 12 П->Х0 X F — 60 12 21 20 21 22 П-Х4- П-Х5 + 64 65 10 30 31 32 Fx n-»Xd Fx 22 6Г 22 40 41 42 F3 П-»ХЗ X 18 63 12 47 48 49 Х-П9 С/П БП 49 50 51 Программа 10 FBx П-Х2 10 0 62 43 44 45 2 =• П-»Х1 02 13 61 50 18 18 06 07 08 09 + F6 Х-»ПсП-»Х1 13 IL 4С 61 16 17 18 19 4- Х-Пв П-»Хс П->Хв 13 4L 6С 6L 26 27 28 29 F9 n->Xd X Х-П6 IE 6Г 12 46 36 37 38 39 Fx + «* ч- 22 10 14 13 46 13 46
Размещение исходных данных и результатов вычислений приведены ниже: Параметры....... F К v, т ra 6ta va тг/180 За Регистры........ О 1 2 3 4 5 6 7 X/ 9 Инструкция по работе с программой 10 Нажимаемые клавиши 1. Включить МК — 2. Перейти в режим „Программирование" F ПРГ 3. Занести программу 4. Перейти в режим „Автоматическая работа” F АВТ 5. Занести исходные данные в соответствующие регистры памяти: F - в регистр 0 FX-ПО К- ” 1 КХ-П1 v, - ” 2 v0 X —П2 m - ” 3 тХ —П 3 6ta 5 61аХ-П5 п/180 — в регистр 7 Firl80-i-X-n7 6. Произвести очистку регистра 4 (если МК был СхХ-П4 в работе) 7. Произвести первый цикл вычислений В/0 С/П 8. Снять с индикатора значение Sa для первого — цикла 9. Произвести индикацию va первого цикла П-Х6 10. Произвести второй и последующие циклы вы- С/П числений 11. Снять значение За 12. Произвести индикацию va, как указано в п. 9 П -> X 6 П р и м е ч а н и сДПсрвый цикл вычислений длится 18 ... 20 с, второй - 9 ... 11 с. Если за указанный период вычисления не завершатся, необходимо нажать клавишу С/П и найти ошибку, допущенную при записи программы или при занесении исходных данных. 2. При отрицательном значении va необходимо: занести в регистр 4 прежнее значение Га, т. с. на 6га меньшее последнего значе- ния, при котором получен отрицательный результат; выбрать новое значение 6га1, в несколько раз меньшее 6га, и занести его в ре- продолжить вычисления. Пример. Дано: А = 3200 т; К = 2101,5 Н • м-1; v„ = 10 уз = 5,14 м/с; F = = 13 200 кге = 129492 Н; m = 3520 • 10’кг. Требуется рассчитать va и За до полной остановки судна. Решение. Принимаем 6га = 60 с. Выполняем операции, указанные в пп. 1 ... 4 Инструкции. Заносим исходные данные в соответствующие регистры памяти: *> Нажимаемые клавиши Индикация 129492 X —ПО 129492 2101,5 Х-П 1 2101,5 5.14Х-П2 5,14 3520 ВП 0 3 3520000 Х-П 3 60Х-П5 60. Fir 180-гХ —П 7 1,7453 -02 СхХ —П4 0. 47
123 Производим очистку регистра 4 и первый цикл вычислений: В/0 С/П Снимаем значение S3 Производим индикацию va Производим второй цикл вычислений: С/П На индикаторе Sa, м Производим индикацию va, м/с 298 Производим третий цикл вычислений, получаем отрицательное значение va. По- этому задаемся 6ta, =1 с и нажимаем следующие клавиши: Нажимаемые клавиши Индикация 120Х-П4 120. БП18С/П 298. 1,29 -01 Ниже представлены результаты вычислений, которые контрольных чисел при отладке использованы в качестве могут быть программы: 60 120 122 0,056 298,1 121 0,092 298 123 0,02 ...... 298,2 Примечания. 1. Тормозной путь Sa до полной остановки судна (скорость равна 0,02 м/с) составляет 298 м. 2. Задаваясь исходными данными для других режимов, получим соответствую- щие значения va и га. 223 298 4.3. Пассивное торможение судна при движении с тралом Расчетные формулы: К"т + vot S = £т1п , V где Vo — начальная скорость движения судна с тралом, м/с; v — текущее значение скорости судна, м/с; Кт — инерционный коэффициент, характе- ризующий интенсивность торможения судна (определяется при натурных испытаниях); t — интервал времени торможения, с; S - тормозной путь, пройденный судном за промежуток времени t, м. Расчет тормозных характеристик судна с тралом при пассивном тор- можении представлен в Программе 11. Программа 11 00 01 02 03 04 05 06 07 08 09 П-»Х0 П-»Х1 X X-*П 6 П->Х 2 П-* X 3 + Х-ПЗП-Х1 X 60 61 12 46 62 63 10 43 61 12 10 11 12 13 14 15 16 17 18 19 48
П-»ХО + П-Х6 -» -=• Х-»П4П-*Х1 -* т F 3 60 10 66 14 13 44 61 14 - 18 20 21 22 23 24 25 26 27 28 29 П->Х0 X Х-»П5П-*Х4 1 • 9 4 X Х-П4 60 12 45 64 01 0- 09 04 12 44 30 31 32 С/П БП 04 50 51 04 Размещение исходных данных и результатов дано ниже; Параметры.........Кт v, 6t t v S Регистры........... 0 1 2 3 4/X 5/ Y Инструкция no использованию программы. После набора программы и перевода микрокалькулятора в режим „Автоматическая работа” необ- ходимо исходные данные занести в регистры 0, 1 и 2 в любой последо- вательности, а затем произвести очистку регистра 3 нажатием клавиш Сх X -* П 3. После этого перейти к вычислению искомых параметров v и 5: первый цикл - нажатие клавиш В/П и С/П; второй и последующие циклы — нажатие клавиши С/П. После каждого цикла вычислений на индикаторе высвечивается значение текущей скорости судна в узлах (множитель 1,94 записан по ад- ресам 24 ... 27, а операция умножения — по адресу 28), а в регистре Y (и в регистре 5) хранится значение тормозного пути в метрах. Отметим, что значение текущей скорости в узлах хранится в регистре 4. Пример. Дано: к0 =4,5 уз = 2,32 м/с; Кт — 68,52 м; Sf = 20 с. Требуется рас- считать значения v и S для Г = 0 ... 240 с. Решение. Будем полагать, что калькулятор готов к решению задачи, поэто- му приступим к вводу данных и решению задачи. Операция Нажимаемые клавиши Индикация Ввод К? в регистр 0 68.52 Х^ПО 68.52 Ввод в регистр 1 2.32 Х-»П 1 2.32 Ввод St в регистр 2 20Х-*П2 20. Очистка регистра 3 схх-пз' 0. Первый цикл вычислений В/0 С/П 2.68 Индикация S FQ 35.43 Второй цикл вычислений С/П 1.91 Индикация S Ff; 58.67 П р и м е ч а н и е. Остальные циклы вычислений производятся аналогично. Если необходимо уточнить текущее значение г, следует нажать клавиши П -»Х 3. Ниже даны значения тормозных характеристик судна с тралом при пассивном торможении (Кг = 68,52 i; v0 =4,5 уз); г, с...................... 20 40 60 v, уз.................... 2,7 1,9 1,5 S, м..................... 35 59 7,6 80 100 120 140 160 180 1,2 1,0 0,9 0,8 0,7 0,6 90 101 111 120 127 134 664 - 4 49
Глава 5 Рис. 8. Схема определения обсерво- ванных координат по двум радио- навигационным линиям положения ИСПОЛЬЗОВАНИЕ МИКРОКАЛЬКУЛЯТОРА ДЛЯ РЕШЕНИЯ ОСНОВНЫХ ЗАДАЧ СУДОВОЖДЕНИЯ 5.1. Определение места судна по импульсно-фазовым РНС Для повышения точности получения обсервованных координат и уп- рощения процедуры графических построений предлагается программа расчета обсервованных координат по- средством координат точек, выбран- ных из таблиц. Принято считать широ- ту точки С наибольшей, а широту точки А меньше широты точки С. Схема определения обсервован- ных координат представлена на рис. 8, где АВ — первая линия положения, CD — вторая; Мо — обсервованное место (точка пересечения линий по- ложения) . Отметим, что обсервованные ко- ординаты определяются относительно координат точки С. Из схемы на рис. 8 видно, что разность широт (РШ) равна отрезку CG, а отношествие (ОТШ) — отрезку M0G, т. е. РШ = M0Ccos а2; <ро =tPc — MoG = рс ~ РШ; ( РШ Ч ОТШ = Л/оС51па2; ДА = ОТШ-sec 1у>с ——2/’ где а2 - угол пересечения второй линии положения с меридианом; М^С = = М0Е + СЕ. Отрезок М0Е определяется из треугольника M0£L4. Учитывая различные ситуации, для определения отрезков АЕ и М0Е можем записать следующие формулы: АЕ = AF ± CFtg а2; MqE = Jinaj/cos (a2 ±aj). Ниже приведена вспомогательная таблица, служащая для определения знаков: 50
Кл~кВ *C—KD “i CF AX Примечание. Долготы (разность долгот) имеют знаки: „+” - восточная (к востоку); западная (к западу). В необходимых случаях (например, при плавании в восточном полу- шарии, или в южных широтах) судоводителю следует составить таблицу, аналогичную табл. 5.1, исходя из возможной комбинации пересечения линий положения. Таблица 5,1, Размещение исходных данных и промежуточных результатов для расчета обсервованных координат Регистры Исходные данные Промежуточные ре- зультаты Примечание. Расположение точек соответствует занесению исходных дан- ных косвенным методом адресования, начиная с регистра d. Для этой цели в регистр 1 заносится число 14 (см. Инструкцию по использованию Программы 12). При же- лании занесение исходных данных производят обычном способом, однако в первом случае занесение исходных данных удобнее, так как нажимаются одни и те же кла- виши - В/О С/П. Ниже представлены формулы, использованные при составлении Про- граммы 12 и программа расчета обсервованных координат места судна: tga2 = I Хд - Xcl (ч>с + Pd \ ~--------------cosl------------1; I Pc - V>dI X 2 / АЕ= AF ± CFtg a2; М0Е = АЕ sina-i/cos (a2±aj); М0С= MqE + СЕ = M0C+ |<^с — |cos a2; РШ = Af0Ccos a2; / РШ \ ОТШ = AfoCsina2; Да = ОТШ • sec 1<рс - 5 Ро ~Рс — |РШ|; Х0=Хс + ДХ. Программа 12 00 01 02 03 04 05 06 07 08 09 П-Х0 4 + С/П КХ-П1 С/П П-ХбП-ХЬ - Bt 60 13 10 50 1 50 6Г 6 11 0Е 10 11 12 13 14 1.5 16 17 18 19 2 4 П-ХЬ + Г- 8 П-ХСП-Ха - X 4 02 13 6 10 1Г 6С 6- 11 12 13 20 21 22 23 24 25 26 27 28 29 51
F6 *Х-»П1 1 41 30 31 F8 П-»Х8 1Г 68 40 41 Х-ПЬ П-ХС 4 6С 50 51 ВТ П-ХЬ 0Е 6 60 61 + П-Х 1 10 61 70 71 П-Х6 + 66 10 80 81 F 8 £ П-Х9 69 31 П-Х6 66 42 П-Ха 6- 52 F 8 1Г 62 F 7 1С 72 В t ОЕ 82 П-ХЬ 6 П-Х7 - ВТ 67 11 ОЕ 33 34 35 X F X 11 12 22 43 44 45 - П-Xd F8 11 6Г 1Г 53 54 55 v Х-П6 F 13 46 25 63 64 65 X П-ХЬП-Х 12 6L 61 73 74 75 П-ХЬ F8 X 6 1Г 12 83 84 85 F7 X Х-П 1С 12 43 2 -г П-Х7 02 13 67 36 ‘ 37 38 21 14 13 46 47 48 X П-Х9П-Х8 12 69 6Г 56 57 58 П-ХЬ F9 X 6 1Е 12 66 67 68 I С/П + F 8 50 10 1Г 76 77 78 П-Х9 69 14 11 86 87 88 ; с/п п-Х8 *+ 50 68 14 10 39 F 6 1L 49 11 59 С/П 50 69 13 79 Х-П 2 42 89 10 Инструкция по использованию Программы 12 Адреса команд 00 ... 03 04 ... 05 06 ... 59 60 ... 66 67 ... 79 80 ... 86 87 ... 90 Пояснение к структуре программу Операции Перевод дуговых минут в доли градуса и суммирование с градуса- ми угла Запись чисел в адресуемые регистры косвенным способом с исполь- зованием регистра 1. В данном случае в регистр 1 заносится число 14 и запись координат производится в регистры d, с ..., т. е. в сторону уменьшения нумерации Расчет углов а, и а2, отрезка EF и запись их в регистры 1, 6 и Ь соотвественно. При этом значения долготы точки D и широты точки В аннулируются, так как в дальнейшем они не нужны. Это следует учитывать при повторном решении задачи для этих же данных). На индикаторе высвечивается значение EF, а в регистре Y находится значение AF. Для того чтобы вычислить отрезок АЕ, необходимо определить знак EF: при АЕ = AF — EF необхо- димо нажать клавишу /—/ и продолжить вычисления по программе Подготовительная работа к вычислению отрезка М „Е. На индика- торе высвечивается значение угла а,. Если по схеме требуется най- ти сумму углов, то знак угла а, остается без изменений, если разность, то нажимается клавиша /—/ и продолжают вычисления по программе Вычисление обсервованной широты и запись ее в регистр 2 Вычисление разности долгот и запись ее в регистр 3. На индикаторе высвечивается разность долгот ДХ (абсолютное значение). Так как обсервованная долгота вычисляется относительно долготы точки С, то необходимо соответственно определить ее знак. Если он минус, то нажимают клавишу /—/ Вычисляют обсервованную долготу 1. Включить микрокалькулятор, переключатель поставить в положение Г 2. Перевести микрокалькулятор в режим „Программиро- вание” 3. Записать программу в память микрокалькулятора 4. Перевести микрокалькулятор в режим „Автоматическая работа” 5. Занести исходные данные в адресуемые регистры памяти: 5.1. Константу 60 в регистр 0 5.2. Константу 14 в регистр 1 5.3, Координаты точек А, В, С, D последовательно перевес- ти в градусы, затем перенести в регистры d, с, ... 6 соответственно Сначала обрабатывается ширата соответствующей точ- ки, а затем долгота 6. Произвести пуск программы на счет с адреса 06 7. Определить знак отрезка АЕ, если знак отрицательный, то откорректировать и произвести пуск программы с адреса 60 8. Определить знак угла at; если знак отрицательный, то откор- ректировать, произвести пуск программы с адреса 67 9. Определить знак разности долгот; если знак отрицатель- ный, то откорректировать, произвести пуск программы с адреса 87 10. Перевести обсервованную долготу в градусы и минуты Нажимаемые клавиши F ПРГ F АВТ 60 X С/П С/П с/п С/П с/п 11. Извлечь из регистра 2 обсервованную широту и перевести <р° - ее в градусы и минуты П — X О X Пример. Находясь в Тс = 12 ч 20 м по счислению в точке с координатами у>с = 27 50.0 N,XC = 65°55.0*W, измерили интервалы времени по поверхностным сигналам от двух пар станций 9930 - W и 9930 - Y. Получили соответственно Т = = 13 944,0 и Г = 53 179,0 мкс. После предварительной обработки [ 7, с. 128] следую- щие результаты: Точка Широта N 27° 56.4' 27’40.5' 28’00.0’ 27’00.0' Долгота W 66’00.0’1 65° 00.0* Г 65’56.4'1 65° 10.1'J № линии положения I II Требуется рассчитать обсервованные координаты места судна на момент наблю- дений. Решение. Будем полагать, что микрокалькулятор готов к приему исходных данных, а переключатель стоит в положении Г (градусы). Ниже приведен порядок решения задачи с сохранением нумерации пунктов Инструкции. Операции 5.1. Ввод числа 60 в регистр 0 (для МК-56) 5.2. Ввод числа 14 в регистр 1 5.3. Ввод координат точек А, В, С, D в регистры d, с, ... 6 соответственно Нажимаемые Индикация клавиши 60 Х-П О 60. 14Х-П1 14. 53
5.3.1. Ввод широты точки А в ре- гистр d Ввод долготы точки А в регистр С 5.3.2, Ввод широты точки В в ре- Ввод долготы точки В в ре- гистр а 5,3.3. Ввод широты точки С в ре- гистр 9 Ввод долготы точки С в ре- гистр 8 5.3.4. Ввод широты точки D в ре- Ввод долготы точки D в ре- гистр 6 6. Пуск программы на счет с адреса Об 7. Определение знака отрезка АЕ и пуск программы с адреса 60 8. Определение знака угла а, и пуск программы с адреса 67 9, Определение знака разности долгот и расчет обсервованной долготы 10. Перевод обсервованной долготы в градусы и минуты 11. Индикация обсервованной широты и перевод в градусы и минуты Ответ: у, =27°55.5’N; Хо =65°52.4’. 27.564 С/П 66 С/П 27. 40 5 С/П 65 С/П 28. С/П 65.56 4 К + С/П 27 С/П 65. 10 1 К + С/П С/П С/П С/П С/П Кб П~*Х2 Кб 27.94 27.94 66. 66. 27.675 27.675 65. 65. 28. 28. 65.94 65.94 27. 27. 65.168333 65.168333 4.1068622 16.678087 6.6689849-02 -6.6689979-02 65.8,7331 65.52 39842 27.55 47 ... 5.2. Ведение промыслового счисления При ведении промысла возникает необходимость частой смены кур- сов. При этом судоводитель выполняет большой объем графических работ на планшете или карте крупного масштаба. Использование микро- калькулятора существенно упрощает работу судоводителя. Расчетные формулы: •Рс2 - 'Pci + РШ; -\:2 -\:1 + РД; PUI = Д7^с cos ПУ + Д7гт cos Кт; РД = AT'O’cSinny + V? sinKT) secipCp, где<рС1, ХС1 — счислимые координаты предыдущей точки; <рс2, Хс2 — счислимые координаты последующей (текущей) точки; РШ, РД — прира- щение координат; рс, рт — скорость судна по лагу и скорость течения, уз; 54
AT — интервал времени плавания между точками 1 и 2, ч; Кт - направ- ление (курс) течения, град; ПУ — путь при дрейфе, град; (/зср - средняя широта, град, Распределение исходных данных и промежуточных результатов по регистрам памяти и Программа 13 расчета счислимых координат приве- дены ниже. Параметры Регистры . 'Pci дг ПУ Z/4 6,0 00 01 02 03 Программа 13 04 05 06 07 08 09 42 10 14 14 20 25 66 43 12 22 25 44 14 63 16 65 19 23 10 24 62 25 12 26 67 27 13 28 0Е 29 60 30 10 31 40 32 02 34 35 12 40 64 1С 42 43 65 44 57 45 36 41 41 46 37 23 38 63 39 10 50 62 13 10 67 53 БП 66 47 -П 1 41 12 49 1С 48 С/П 50 12 52 С/П 50 Инструкция по использованию программы 13 Нажимаемые клавиши 54 49 49 61 10 67 1. Включить микрокалькулятор, переключатель по- ставить в положение Г 2, Перевести микрокалькулятор в режим „Про- граммирование" 3. Записать программу в память микрокальку- лятора 4. Перевести микрокалькулятор в режим „Автома- тическая работа” 5. Записать константу 60 в регистр 7 (для МК-56) 6. Перевести счислимую широту в градусы и занес- ти в регистр 0 (если широта южная, то после перевода в градусы изменить знак, а затем занести в регистр 1) 7. Перевести счислимую долготу в градусы (если она западная, то изменить знак) и записать ее в регистр 1 8. Записать значение скорости течения (в узлах) в регистр 5 9. Записать направление течения Кт й градусах в ре- гистр 6 10. Перед сменой курса занести в стековые регистры последовательно: путевой угол, скорость судна в узлах и время следования этим курсом в долях часа 11. Пустить программу на счет с адреса 00 12. Снять с индикатора значение долготы в градусах и долях градуса, перевести в градусы и минуты F ПРГ F АВГ 60Х ПУ В t vc В t ДТ В/0 С/П 55
13. Извлечь из регистра 0 широту в градусах и пере- П -♦ X О вести дробную часть в минуты Пример, а) Тс = 11 ч 40 м; <рс = 43°20.0*N; Хс = 65°10.0’W. Судно легло на курс Кс = 65°, скорость судна по лагу vc = 12,3 уз, течение не учитываем; б) Тс = 17 ч 50 м, Кс = 68°, течение учитываем; vT = 1,3 уз, Кт = 150°. Требуется определить счислимые координаты па Тс = 17 ч 50 мин. Решение. Будем полагать, что микрокалькулятор готов к решению задачи, поэтому приступаем к вводу данных, сохраняя нумерацию операций, приведенных в Инструкции по использованию Программы 13. Операции Нажимаемые клавиши Индикация 5. Ввод константы (для МК-56) 60Х-П 7 60.0 6. Ввод <₽с в регистр 0 43.20 К + 44.333333 Х->П0 7. Ввод Хс в регистр 1 65. 10 К + 65.166667 1-1 Х-П 1 -65.166667 8, Ввод i>T = 0 в регистр 5 СхХ->П5 0 9. Пропускаем, так как vT = 0 J0. Ввод в стековые регистры ПУ, vc, ДТЧ 65 В t 12,3 В t 6.17 11. Пуск программы с адреса 00 В/0 С/П -63.58368 12. Перевод Хс в градусы и ми- Кб -63.3502 нуты 13. Извлечение <рс из регистра 0 п->хо 43.867882 и перевод в градусы и минуты Кб 43'52.07 Ответ: на Тс = 17 ч 50 м у>с = 43°51.1*N; Хс = 63° 35,0’W. Пример. По данным предыдущего примера рассчитать счислимые координаты на Тс= 19 ч 30 мин с учетом течения и нового курса Кс = 68° . Решение. Операции Нажимаемые клавиши Индикация 1. Ввод vT в регистр 5 1,3 Х-П5 1.3 2. Ввод Кт в pei-истр 6 150Х-П6 150. 3. Ввод ПУ, vc, Д7Ч в регистры 68 В1 68. стека 12,3 В t 12,3 1,67 1,67 4. Пуск программы с адреса 00 В/0 С/П 63.117918 5. Перевод долготы в градусы и Кб -63.0707 минуты 6. Извлечение широты из регист- П-ХО 43.964793 ра 0 7. Перевод широты в градусы и Кб 43.57.887 минуты Ответ: на Тс = 17 ч 50 мин v>c = 43°57.9’N; Хс = 63°07.1’W. Краткие пояснения по структуре Программы 13 Адрес Операция 00 ... 04 Запись исходных данных ДТ, vc, ПУ из стековых регистров в регистры памяти 2, 3, 4 соответственно 05 ... 07 Вычисление разности широт без учета течения 08 ... 13 Расчет разности широт от влияния течения. Если течение нс учи- тывается, то срабатывает команда условного перехода и коман-. ды 10 ... 13 выполняться не будут 56
14 ... 18 19 ... 22 23 ... 28 29 ... 33 34 ... 39 Расчет алгебраической суммы разности широт в минутах (ми- лях) и перевод в градусы и доли градуса Расчет широты у>сз и запись ее в регистр 0 Расчет секанса средней широты Расчет отшествия в милях без учета течения Расчет отшествия в милях от влияния течения. Если течение не 40 ... 44 45 ... 47 49 ... 51 53 ... 54 учитывается, то команды по адресам 37 ... 39 не выполняются, так как срабатывает команда условного перехода Расчет алгебраической суммы разности долгот в долях градуса Расчет долготы ХСз в градусах и долях градуса Фрагмент программы для перевода угла, выраженного в граду- сах и минутах, в градусы и доли градуса Команда безусловного перехода на адрес 49. Она позволяет упростить процедуру перевода долготы и доли градуса. 5.3. Плавание по дуге большого круга Порядок расчета начального и промежуточных курсов. С помощью МК вычисляют: ортодромическое расстояние между начальной и конеч- ной точками S, а также начальное ортодромическое направление А в полу- круговом счете; ортодромическую поправку, служащую для перевода начального ортодромического направления (полукругового счета) в лок- содромическое S Ф1 =~^ tg<pxsin^opT1; локсодромическое направление и соответственно курс по формулам Ллок 1 = ЛорТ1 + ; К. = Ллок (NE) = 360° - Ллок i (NW); ортодромическое направление для следующего курса: Лортг = Ллок 1 + Ф>; широту первой точки поворота 'Pi = Ра +5’cos^nOKi. где - широта начальной точки; S - длина хорды, град.; ЛпОк1 - лок- содромическое направление от начальной точки. Ортодромическая поправка для расчета локсодромического на- правления от первой точки поворота S Фг =------tg р2 Sin А орт2; 57
локсодромическое направление и курс при плавании по второй хорде •‘4лок2 — ^орт2 + 'Pit К2 = Ллок2 (NE) = 360° - Ллок2 (NW). Примечания. 1. Расчет локсодромического направления и кур- сов для последующих хорд осуществляется аналогично. 2. Прокладка хорд производится следующим образом. От исходной точки прокладывают курс, по нему - локсодромическое расстояние S. Из конца данного отрезка прокладывают следующий курс, а по нему - локсодромическое расстояние S и т. д. З.При наличии большой невязки, обнаруженной в процессе плавания, необходимо повторить расчеты. Если же невязка будет незначительная (3 ... 5 миль), то новый курс следует проложить от обсервованной точки параллельно старого курса. Расчетные формулы определения ортодромического расстояния и начального ортодромического направления. Ортодромическое расстояние для плавания на сфере вычисляют по формуле cos 5орт = sin sin<PB + cos <рА cos i/^cos РД, где 5орт — расстояние по ортодромии между начальным и конечным пунктами плавания, град; Фв ~ широты начального и конечного пунктов плавания, град; РД — разность долгот, град. Следует иметь в виду, что широта >рА начального пункта считается всегда положительной, а широта конечного пункта \рв может быть как положительной, (если <рв одноименна с <рА). так и отрицательной (если разноименна с tpA). Начальное ортодромическое направление (полукругового счета) вычисляют по формуле cos Лорт = (sin<pB - sin</u cos5OPt)/cos >рА sinSopT, где Лорт - начальное ортодромическое направление в полукруговом счете. Первая буква Лорт всегда одноименна с <рА, а вторая — с разностью долгот. Программирование задачи и вычисление по программе. Распределим исходные данные и константы по регистрам памяти: число 60 — в регистр 0; <Ра ~ 1; 'Рв - 2; РД - 3; 5-6, где 5 - плавание по хорде, выражен- ное в градусах дуги; Промежуточные результаты занесем в следующие регистры памяти: 5 - в регистр 4; Лорт1 - в регистр 5. Вычисление текущей ортодромической поправки (|Д- = 0,5 5° tg щ sin -40рт|) оформим в виде подпрограммы, которую запишем с адреса 37. 58
Программу вычисления 5орт и Лорт запишем по адресам 04 ... 35. Напомним, что программа по адресам 00 ... 03 служит для перевода угла, выраженного в градусах и минутах, в градусы. Приведем Програм- му расчета элементов дуги большого круга. Программа 14 Инструкция по работе с Программой 14 1. Включить МК 2. Передать в режим „Программирование” F ПРГ 3, Занести программу 4. Передать в режим „Автоматическая работа” F АВТ 5. Занести исходные данные в регистры памяти: 5.1. Число 60 - в регистр 0 (для МК-5.6) 5.2. ” 1 ’ 5.3. ” 2 5.4. РД ”3 5.5,5° ”6 6. Вычислить 5орт и Аорт БП04С/П 7. Вычислить Алок для первой хорды БП 37 С/П 8. Определить значение ср, П -»Х 1 9. Вычислить Алок для второй хорды БП 37 С/П 10. Определить значение П-»Х1 Примечание. Последующие значения f, и АЛ0К|- вычисляют аналогично, т, е. выполнением операций, указанных в пп. 9 ч 10. Пример. Дано: = 48° 10*14; Хл = 6°20*W; ~рВ = 22°50*N; Кв = 82°25'N. Требуется выполнить расчеты для плавания по дуге большого круга способом хорд, если длина хорды 5 равна 120 милям. Решение. Выполняем действия в соответствии с Инструкцией. Операции Нажим'аемые клавиши 1. Включаем МК 2. Переходим в режим „Программирование" F ПРГ 59
5. Заносим исходные данные и константы в соот- ветствующие регистры: 5.1. Число 60 - в регистр О 5.2. -РА " 1 5.3. <₽в ” 2 5.4. РД ” 3 5.4.1. Вычисляем РД 60 Х-П о 48. 10 К + Х-П 1 22. 50 К + Х-П 2 82. 25 К + 6. 20 К + 5.4,2. РД - в регистр 3 5.5.5° ” 6 6. Вычисляем 5орт и ЛорТ1 7. Вычисляем ЛЛОК1 для пер- вой х.орды 8. Определяем значение Х-П 3 2Х-П6 БП 04 С/П БП 37 С/П 9. Вычисляем Ллок для второй хорды 10. Определяем 11. Определяем А лок для треть- ей хорды 12. Находим БП 37 С/П П-Х1 Кб БП 37 С/П П-Х1 Кб 60. 48.166667 22.833333 82.416667 6.3333333 76.083334 76.083334 84.02 ..... 85.13...... 48.33631° 48.202’ 87.37 ..... 48.428197° 48.25' 89.62 48.441631° 48.22 9* Примечания. 1. Аналогичным способом рассчитывают искомые значения для последующих хорд. Имеет смысл рассчитывать не более 3-4 хорд, так как в случае невязок, неизбежно возникающих при плава- нии, особенно в океане, потребуются новые расчеты. Если при вычислении очередного локсодромического направления произошел сбой, т. е. допущена ошибка, то необходимо занести >рА в ре- гистр 1, а Аорт — в регистр 5. После этого повторить весь цикл вычисле- ний, начиная с п. 7, т. е. нажать клавиши БП 37 С/П и цолучить Алок для очередной хорды, 2. Для перевода локсодромических направлений в курсы используют соотношения: К = Ane; К = 360° — Anw Таким образом, получим: Ллок, град....... 85,1 87,4 89,6 91,9 94,1 96,4 К................ 274,9 272,5 270,4 268,1 265,9 263,6 Для индикации нажмем клавиши П -> X 4, после чего получим 64,088797°. Для перевода в мили нажмем клавиши 60 X, тогда S = 3845,3 мили. Для определения выигрыша воспользуемся табл. 23 МТ-75. Получим AS = 29 миль. После расчета истинных курсов прокладывают их, начиная с точки Я, а по ним откладывают расстояния S = 120 миль. 60
Для расчета выигрыша в расстоянии (Д5) можно составить програм- му, используя какую-либо приближенную формулу. Однако учитывая ограниченную память МК, с одной стороны,а с другой,— простоту решения этой приближенной задачи с помощью табл. 23 МТ-75, которая всегда под рукой, целесообразнее вбспользоваться приведенной рекомендацией. 5.4. Расчет девиации магнитного компаса Расчетные формулы: 5 = A + 5sinK + CcosK + Dsin2K + £'cos 2К, где 5 — девиация магнитного компаса, град; А, В, С, D, Е — коэффициен- ты, вычисляемые по формулам: А = 1/8 1 8,-; В = 0,25 [(53 - 8Т) + 0,71 (62 + S4 — 56 — S8)]; С = 0,25 [(8, - 65) + 0,71 (62 - 84 - 86 +68)]; D = 0,25- (82 — 34 + 86 — 8g); Е = 0,25- (61 - 83 +6s -67); 61,82, . .., 88 — значения девиации, полученной из наблюдений на курсах, с интервалом 45°, начиная от 0. Будем полагать, что коэффициенты целесообразно вычислять в обыч- ном режиме. Что касается девиации, то она вычисляется по Программе 15, размещение исходных данных приведено ниже: Параметры........... Регистры............ ABODE ДК К 14 6 d с Ь а 9 8 7 0 Х/1 Программа 15 00 01 02 03 04 05 06 КХ-»П0 С/П П-»Х8 П-*Х7 + Х-П7 ВТ L0 50 68 67 10 47 0Е 10 11 12 13 14 15 16 F8 П->ХЬ X + П->Х7 2 14 1Г 6 12 10 67 02 20 21 22 23 24 25 26 П-Ха X f F8 П-»Х9 X + 6- 12 14 1Г 69 12 10 30 31 32 33 34 35 36 Х->П1 С/П БП 02 41 50 51 02 07 08 09 F7 П->Хс X 1С 6 12 17 18 19 X В t F 7 12 0Е 1С 27 28 29 + П->Хё + 10 6Г 10 37 38 39 61
Инструкция по использованию Программы 15 1. Включить микрокалькулятор, переключатель по- ставить в положение Г (градусы) 2. Перейти в режим „Программирование” 3. Записать программу в память микрокалькулятора 4. Перейти в режим „Автоматическая работа” 5. Занести число 14 в регистр памяти О 6. Занести коэффициенты А, В, С, D, Е в регистры d ... 9 соответсгвенно 7. Занести ДК = 10° в регистр 8 8. Произвести очистку регистра 7 9. Вычислить значение девиации для диапазона К = = 10 - 360° 'Нажимаемые клавиши F ПРГ F АВТ 14Х-П0 А В/0 С/П Е В/0 С/П 10 В/0 С/П СХ В/0 С/П БП 02 С/П Примечание. После завершения очередного цикла вычислений на индика- торе высвечивается значение девиации компаса в градусах, оно же хранится в регист- pc 1. Значение текущей величины К хранится в регистре 7. Пример. Из наблюдений получены следующие значения девиации: К°.......... 0 45 90 135 180 225 270 315 6 .......... 1,0 0,8 0,0 -1,0 -0,8 0,6 1,4 1,2- Требуется рассчитать таблицу девиации компаса. Решение. Вычислим значения коэффициентов А ... Е: А = 1/8 • (1,0 + 0,8 - 1,0 - 0,8 + 0,6 + 1,4 + 1,2) = 0,4° ; В = 1/4 = 1-1,4 + 0,71 (0,8- 1,0- 1,2)1 = -0,7°; С = 1/4 • [ (1,0 + 0,8) + 0,71 (0,8 + 1,0 - 0,6 + 1,2)) = 0,9° ; D = 1/4 • (0,8 + 1,0 + 0,6 - 1,2) = 0,3° ; £=1/4» (1,0-0,8- 1,4) = -0,3°. 2. Выполним действия, указанные в пп. 1 ... 8 Инструкции. 3. Выполним п. 9 Инструкции: Нажимаемые клавиши Индикация 3.1 С/П (К= 10°) 9.85 -01 3.2 с/п (К =20°) 9.69 -01 3.3 с/п (К = 30°) 9.39 -01 3.4 с/п (К = 40°) 8.82 -01 3.10 с/п (К = 100°) -2.66 -01 Примечания. 1. Аналогично выполняются все остальные циклы вычислений. 2. Приращение ДК выбирается самим пользователем. При ДК = 45° можно произвести контроль правильности расчета, так как полученные значения девиации должны быть равны соответствующим значениям наблюдаемых величин. В Программе 16 дан расчет коэффициентов девиации магнитного компаса. 62
00 <->по 10 01 С/П 50 02 Сх ог 03 с-по го Программа 16 04 05 10 5L 15 06 03 03 07 08 08 09 6С 20 10 22 68 23 24 10 30 67 32 04 34 66 25 К-П0П L0 35 26 61 27 28 6L 29 40 66 50 кх-поп L0 10 42 61 43 6Г 45 53 54 68 55 60 6Г 61 С/П 50 6L 69 10 6С 36 10 46 56 67 37 38 68 39 48 04 49 66 58 04 59 В Параметры Регистры . Примечание. Многоточие означает пропуск символов 6. соответственно. Инструкция по использованию программы. После ввода программы и перехода в режим „Автоматическая работа” выполняют следующие операции: заносят число 14 в регистр 0; заносят девиации 6] . .. 68 в регистры d ... 6 последовательно внесением значения девиации в регистр X (на индикатор) и нажатием клавиш В/0 С/П; заносят снова число 14 в регистр 0; записывают число 8 в регистр 2; записывают число 0,71 в регистр 1; вычисляют коэффициенты А, В, С, D, Е нажатием клавиши С/П (пуск программы в адрес 02). После завершения вычислений на индикаторе высвечивается коэф- фициент Е. Для индикации коэффициента D необходимо нажать клави- ши FQ. Коэффициенты А, В, С извлекаются из регистров памяти 5, 4, 3 соответственно. Ниже приведены результаты решения примера: К°........ 0 45 90 135 180 225 270 315 Л В С D Е 6 ........ 1,0 0,8 0 -1,0 -0,8 0,6 1,4 1,2 0,4 -0,7 0,9 0,3 -0,3 63
5.5. Расчет элементов движения встречного судна по данным радиолокатора В правиле 19 (d) МППСС-72 указывается, что если при ограничен- ной видимости на экране радиолокатора (РЛС) появился сигнал другого судна и существует опасность столкновения, то необходимо предпринять своевременные действия для выхода из создавшейся ситуации. В предла- гаемой методике принимается, что допустимым сближением является расстояние Дзад, равное 2—3 милям. Программой предусмотрена такая последовательность расчета: сначала по данным двух наблюдений за ин- тервал времени 3 мин вычисляется ожидаемое кратчайшее расстояние сближения Дкр. Если Дкр меньше 2-3 миль, то вычисления продолжа- ются вплоть до расчета курса, на который должно лечь судно спустя 3 мин после вторых наблюдений. Расчетные формулы: П? =«П°2 0° = ----------- ; Дкр = Д1 sin/3; Д1/Д2-1 Д1 - д2 /?О=Г11=КС; Д| _ j = , Д] _ з = 2-Д1_2‘, cos (3 Г В = (Д1 - 2 + Ус - 2 • Д1 _ 2VCCOS0)1/2; Д1 — 2 sine? = sin0o--- 5 = Дзад/Дз; ^в Дз = Д1 -Д1_3;'Кв = <7±Кс±18О°; К = П2 + (П, - П3) +<7р+ 180°; sin0 = sin (К - Кв)---: •'с Кс1 = П2 + (П1 -П2)+0 + <7р. Здесь П1, П2, Д], Д2 - пеленги и дистанции в моменты времени первых и вторых наблюдений соответственно, град, мили; Дкр — кратчайшее расстояние сближения между судами, мили; Д] _ 2 — относительная скорость сближения, уз/3 мин.; Д] _ 3 - расстояние по линии относитель- но движения (ЛОД) от первой точки до точки изменения ЛОД, мили; ус - скорость судна, уз/3 мин; vB - скорость встречного судна, уз/3 мин; Дзад — заданное кратчайшее расстояние сближения между судами; Кв - курс встречного судна, град; КС1 — курс судна, обеспечивающий рас- хождение с встречным судном на заданном расстоянии Дзад; углыД,0о, q показаны на рис. 9, а, а угол q - на рис. 9, б. 64
При расчете курса встречного судна Кв необходимо учесть следующие правила. Судно обнаружено справа (П2 >-Кс): —180° при «у + Кс > 180°; t 180° при q + Кс < 180°. Судно обнаружено слева (П2 < Кс): — 180° при |<7- Кс| < 180°; + 180° при |q--KJ >180°. Если получено отрицательное значение, нужно изменить знак числа нажатием клавиши /—/. Для расчета элементов движения встречного судна и курса расхожде- ния с ним составлена Программа 17. Размещение исходных данных приве- дено в табл. 5.2. Таблица 5.2. Размещение исходных данных и промежуточных результатов для расчета курса расхождения со встречным судном Параметр Обозначение Символ (значение) регистра Первая дистанция Первый пеленг Вторая дистанция Второй пеленг Курс своего судна Скорость своего судна Относительная скорость Модуль разности пеленгов Вспомогательный угол Константа Кратчайшее расстояние Скорость встречного судна Вспомогательный угол Разность пеленгов Вспомогательный угол Сумма углов 664 —5 65
Краткое пояснение по структуре Программ^ 17 Адрес команд 00 ... 03 04 ... 06 07 ... 14 Операция Вычисление разности пеленгов Вычисление модуля разности пеленгов и запись в регистр 7 Вычисление вспомогательного угла (3 и занесение его в ре- 15 ... 19 20 ... 24 24 ... 46 гистр 8; попутно вычисляется разность Д, - Д2, которая зано- сится в регистр 6 Вычисление Дкр в милях и запись его в регистр а Вычисление относительной скорости (мили/3 мин) Д1 _ 2 и запись в регистр 6 Вычисление скорости встречного судна (мили/3 мин) и зане- сение значения в регистр. Попутно (команды по адресам 29 ... 36) вычисляется вспомогательный угол /30 и его значе- ние записывается в регистр 8, где до этого находилось значение 47 ... 54 55 ... 61 62 ... 66 67 ... 69 угла (3 Вычисляется угол q' и его значение записывается в регистр О Вычисляется угол qp и записывается в регистр а, в котором до этого находилось значение Дкр Вычисляется угол q + П2 + (П2 - 11,) и записывается в ре- гистр 7, в котором до этого находилось значение угла П2 ~ П2 Подготовительная операция для определения значения q ± Кс: из регистра С извлекается q, из регистра 4 — Кс, после чего микрокалькулятор прекращает вычисления (см. п. 2 приме- чания) Вычисление курса расхождения КС1 и его индикация 81 66
Примечания. 1. Команды по адресам 04, 05 и 32, 33 служат для определе- ния модуля разности пеленгов и курсового угла обнаруженного судна. 2. После выполнения команды по адресу 69 микрокалькулятор прекращает вычисления, на индикаторе высвечивается значение Кс. Теперь необходимо вручную вычислить курс встречного судна (см. п. 9 Инструкции по использованию Програм- мы 17) и после этого перейти к выполнению последней части программы. Инструкция по использованию Программы 17 1. Включить микрокалькулятор, переключатель по- ставить в положение Г 2. Перевести микрокалькулятор в режим „Автома- тическая работа” 3. Занести программу в памятымикрокалькулято- ра 4. Перевести микрокалькулятор в режим „Автома- тическая работа” 5, Занести исходные данные в соответствующие ре- гистры памяти микрокалькулятора: 5.1. Курс своего судна в регистр 4 5.2. Скорость своего судна, мили/ч (деленные на 20) в регистр 5 5.3. Первую дистанцию (в милях) в регистр 0 5.4. Первый пеленг в регистр 1 5.5. Вторую дистанцию (в милях) в регистр 2 5.6. Второй пеленг в регистр 3 6. Пуск программы с адреса 00 7. Снять значение Дкр, Если Дкр меньше 3 миль, выполнить п. 8 Инструкции 8. Пуск программы с адреса 20 9. Убедиться в том, что после завершения вычисле- ний на индикаторе высветилось значение курса своего судна (Кс), и приступить к вычислению курса встречно- го судна, используя следующие правила: 9.1. Встречное судно справа Всгречное судно слева 9.2. Полученное значение меньше 180 Полученное значение больше 180 Полученное значение по модулю меньше 180 Полученное значение по модулю больше 180 9.3. Если число (Кв) отрицательное, то изменить Нажимаемые клавиши F ПРГ F АВТ С/П /-/ 10. Пуск программы с адреса 70 С/П 11. Снять с индикатора значение курса расхождения со встречным судном 12. Извлечь из регистра скорость встречного судна и П - X Ь перевести в узлы 20 X Пример. Судно следует курсом Кс = 48°, скорость vc = 15,4 уз. В Тс = = 9 ч 15 мин на экране радиолокатора обнаружена цель. Произведение измерения: 9 ч 18° мин 9 ч 21 мин Дистанция, мили Пеленг, град 16,6 17,6 67
Требуется рассчитать элементы движения встречного судна и, если необходимо, рассчитать курс судна, на который следует лечь в 9 ч 24 мин. с тем чтобы кратчайшее расстояние сближения между судами Дкр было не менее 3 миль. Решение. Будем полагать, что программа занесена в память микрокальку- лятора, который находится в режиме „Автоматическая работа”, переключатель в положении Г. Далее предположим, что некоторые исходные данные уже занесены в соответствующие регистры памяти: курс судна Кс = 48° - в регистр 4; скорость судна vc = 0,75 уз/3 мин - в регистр 5; константа 180 - в регистр 9. Введем наблюденные данные в соответствующие регистры памяти: Параметр Д. П1 Регистр 0 Нажимаемые клавиши 9,1 Х-»П 0 16,6 Х-П 1 8Х->П2 17,6 Х-ПЗ После ввода данных в регистры памяти 0 ... 3 производим пуск программы на счет с адреса 00 (нажатием клавиш В/0 С/П). В результате выполнения части программы до адреса 19 вычисляется кратчайшее расстояние Дкр в милях и индици- руется в регистре X. В нашем примере Дкр = Г,15 мили. Так как Дкр < 3 миль, продолжаем вычисления'нажатием клавиши С/П. После завершения выполнения части программы с адреса 20 по адрес 69 на индикаторе высветится значение Кс, а в регистре Y будет храниться значение угла q. Руководствуясь Инструкцией по использованию Программы 17, рассчитаем курс встречного судна (цели). Нажимаемые клавиши - (так как цель слева) П-Х9 - (q- Кс< 180) /-/ Следовательно, Кв = 147,5°. После этого рассчитываем КС1: С/П Ответ: КС1 =110,2°. Индикация 32.5 180. -147.48 147.48 110.2 Таким образом, чтобы разойтись с данной целью на расстоянии Дзад = 3 мили, судно должно к моменту времени Тс = 9 ч 24 мин лечь на курс 110,2°. Примечание. Для индикации скорости цели необходимо извлечь число из регистра b (нажать клавиши П -► X Ь), а затем умножить его на 20. В результате получим скорость в узлах. В нашем примере П-ХЬ 7.025 20 X 14.05 (уз) Теперь рассмотрим пример, когда цель расположена справа. Дано: Кс = 27°; vc = 14,5 уз; Г, = 10 ч 26 мин; II, =38°; Д; =11,5 уз; Г, =10 ч 29 мин; П, =36° ; Д2 = 10,0 мили. Требуется рассчитать Дкр, а при необходимости КС1. 68
Решение. Предположим, что, как и в предыдущем случае, микрокалькуля- тор готов к решению задачи. Нажимаемые клавиши П 4 (Кс) Индикация 20 i) 7.25 - 01 7.25-01 П 9 (если до этого в регистре 9 находилось 180. другое число) В/0 С/П (ДКр) (так как Д < 3, продолжим вычисления) С/П + (цель справа) 38. 10. 36.0 2.65 + (цель справа) С/П (КС1) 70.06 180. 250.06 34.99 Ответ: Дкр = 2,65 мили; КС1 =35°; vB = 18,6 уз. Примечания. 1. Чтобы разойтись с целью на заданном расстоянии (Дзад = = 3 мили), необходимо, чтобы к моменту времени Тс = 10 ч 32 м судно легло на курс КС1 =35°. 2. Если скорость относительного движения судна и цели определяется за 6 мин, то при вводе гс в регистр 5 необходимо скорость в узлах разделить на 10. Соответ- ственно при переводе ve в узлы необходимо число, извлеченное из регистра В, умно- жить на 10. 3. Для повторною вычисления курса расхождения со встречным судном нет необходимости выполнять всю программу. Достаточно извлечь на индикатор угол О из регистра А, нажать клавиши БП 62 С/П, на индикаторе высветится значение Кс. После этого необходимо вычислить значение Кв, а затем нажать клавишу С/П. После завершения вычислений получим значение КС1. 5.6. Расчет средней квадратической погрешности результатов наблюдений Общие положения. В соответствии с ГОСТ 16263-70 „Термины и опре- деления” погрешностью измерения считается отклонение результата измерения от истинного значения измеряемой величины. Все погрешности делят на случайные, систематические и грубые. Случайная погрешность — это составляющая погрешности измерения, изменяющаяся случайным образом при повторых измерениях одной и той же величины. Систематическая погрешность - составляющая погрешности измере- ния, остающаяся постоянной или закономерной, изменяющаяся при повторых измерениях одной и той же величины. Грубая погрешность - погрешность измерений, существенно превы- шающая ожидаемую погрешность при данных условиях наблюдений. 69
Инструментальная погрешность измерений — составляющая погреш- ности измерения, зависящая от погрешности применяемых средств изме- рений. Погрешность метода измерений — составляющая погрешности изме- рений, происходящая от несовершенства метода измерений. В дальнейшем будем рассматривать только случайные погрешности измерений, подчиненные нормальному закону распределения. Моделирование погрешностей нормального распределения, Для прове- дения качественного эксперимента в учебных целях [3,с. 76—77] необхо- димо уметь воспроизводить случайные погрешности нормального распре- деления. Для этого рекомендуется следующий прием: сначала получить числа равномерного распределения и посредством их числа (погрешности) нормального распределения [2]. Исходные формулы: >'«+ 1 = О(1Ьп + я); Д„ = [21n (l/vn + j)]°-ssin2jrv„, (5.1) (5.2) где Vn, v„+ 1 — предыдущие и последующие положительные числа равно- мерного распределения в интервале 0 . . . 1; D ( ) - дробная часть числа, заключенного в скобках; Д„ — случайная погрешность, подчинен- ная нормальному закону распределения. Отметим, что с помощью формулы (52) генерируется серия чисел с характеристиками о = 1, До =0, где До - погрешность арифметического среднего (математическое ожидание). Расчет чисел нормального распределения представлен в Программе 18. Программа 18 00 01 02 03 04 05 06 07 08 09 п-хо 1 х Г Я + х ->П<1 КП-Xd —ь П-Xd 60 01 01 12 20 10 4Г ГГ 14 6Г 10 11 12 13 14 15 16 17 18 19 ы п-хо 2 X г я X F 7 —> Х-ПО F 1/Х 60 02 12 20 12 1С 14 40 23 20 21 22 23 24 25 F In 2 X F— X С/П 18 02 12 21 12 50 Инструкция по использованию программы 18 Нажимаемые клавиши 1. Включить МК, переключатель поставить в положе- ние Р 2. Перевести МК в режим „Программирование” 3. Записать программу 4. Перевести МК в режим „Автоматическая работа” 5. Занести случайное (положительное) число, мень- шее 1,в регистр 0 6. Пустить программу на счет с адреса 00 F ПРГ F АВТ 0.1234 Х-ПО В/0 С/П 70
Примечания. 1. После завершения первого цикла вычислений на индика- торе появится число 8.25 ... -01,т. е. \ =0.825. 2. Последующие циклы вычислений после нажатия клавиш дадут следующие результаты: 0.006; -1,654; 0.0628; -0.14267; 0.2624; -0.408. 3. Цикл вычислений длится около 12 с. Если на этот период вычисления не завершатся, то необходимо нажать клавишу С/П, что приведет к аварийной останов- ке процесса вычислений. После этого необходимо найти ошибку в программе и во- зобновить вычисления. 4. Предложенная программа позволяет рассчитать 6000 неповторяющихся чисел равномерного распределения, а на их основе - соответственно числа нормаль- ного распределения. Расчет средней квадратической погрешности результатов наблюдений (1-й вариант). Исходные формулы: 1 п 10 =----- S Vi -l0; п 1 т0 = [ S Vil(n- I)]0,5; т0 = rn/y/n, где lj — измеренные величины; п — количество измеренных величин, содержащихся в данной серии; т — средняя квадратическая погрешность одиночного измерения; т0 - средняя квадратическая погрешность ариф- метического среднего; г,- - отклонения от арифметического среднего. Ниже представлены табл. 5.3 с исходными данными и Программа 19 для расчета средней квадратической погрешности. Таблица 5.3. Размещение исходных данных и результатов Параметр Обозначение Символ регистра Измеренная величина ‘i d, с ... 5 Количество измеренных величин п 1,2 Константа 14 0 Сумма измеренных величин Slj 3 Сумма квадратов отклонений S17 4 Арифметическое среднее 1ч X Средняя квадратическая погрешность арифмети- ческого среднего то Y Средняя квадратическая погрешность одиночно- го измерения т Z 71
Программа 19 Краткие пояснения к программе Диапазон команд Операция 00 ... 01 Косвенный способ записи измеренных величин в адресуемые регистры 02 ... 07 Расчет суммы измеренных величин посредством команд орга- низации циклов, запись полученной суммы в регистр 3 08 ... 09 Запись числа п в регистр 1 10 ... 11 Расчет среднего арифметического измеренных величин, запись в регистр 3 12 ... 14 Запись числа 14 в регистр С 15 ... 17 Расчет отклонений от арифметического среднего с использова- нием чисел из адресуемых регистров 18 ... 23 Расчет суммы квадратов отклонений от арифметического сред- него и запись ее в регистр 4 с использованием команды органи- зации цикла 24 ... 29 Расчет средней квадратической погрешности одиночного изме- рения, запись копии в регистр 4 ,29 ... 32 Расчет средней квадратической погрешности арифметического среднего 33 Извлечение арифметического среднего из регистра 3 для его очистки и обработки последующей серии измеренных величин. При этом арифметическое среднее хранится в стековом регистре 34 ... 36 Перезарядка регистра 0 37 ... 39 Очистка индикатора и регистров 3, 4 40 Перемещение чисел в регистрах стека, с тем чтобы избавиться от нуля, находящегося в регистре X 41 ... 42 Перезарядка регистра 1, т. е. запись числа п в регистр 1, при этом нуль из регистра Т исчезнет. Вернувшись к операции по адресу 40, отметим, что при ее невыполнении была бы потеряна полезная информация - число m После решения задачи измеренные величины мых регистрах, как указано ниже: сохраняются в адресуе- Регистрь................... d с Ь а 9 8 7 6 5 i.......................... I 2 3 4 5 6 7 8 9 72
Информация может служить „ключом” при проверке правильности занесения значений измеренных величин в адресуемые регистры памяти. Так, например, для того чтобы проверить правильность занесения /5, не- обходимо нажать клавиши П -*Х 9. Если в серии измеренных величин содержится более 9, необходимо использовать 2-й вариант программы. Расчет средней квадратической погрешности без ограничения коли- чества измеренных величин (2-й вариант). Расчетные формулы: ли ' 1л Е/,-(Е//)2/л 10= — Ъ It, 11 л 1 л - 1 т0 = т!\/п. Программа 20 39 40 тС/П Инструкция по использованию Программы 20 1. Включить микрокалькулятор (переключатель по- ставить в положение Г) 2. Перевести микрокалькулятор в режим „Про- граммирование” 3. Записать программу в память микрокалькуля- 4. Перейти в режим „Автоматическая работа” 5. Установить счетчик адресов на адрес 00 6. Поочередно (одно за другим) ввести измеренные величины, произвести пуск программы 7. После ввода всех измеренных величин (отсчет на индикаторе равен л) произвести пуск программы с адре- са 19 Нажимаемые клавиши F ПРГ F АВТ В/0 I/ С/п шгс/п (БП 19 С/П) FQ 8. После завершения цикла вычислений снять с ин- дикатора значение т, 9. Перевести из регистра Y в регистр X (на индика- тор) т 10. Перевести из регистра У в регистр X 664 — 6 73
Примечания. 1. Измеренные величины второй и последующих серий (или повторное решение данного примера) обрабатываются в соответствии с Инструк- цией, начиная с п. 5. 2. После решения очередной задачи исходные данные не сохраняются, так как это происходит в 1-м варианте Программы. 3. Клавиша ШГ (см. п. 7 Инструкции) нажимается для того, чтобы „обезгла- вить” команду безусловного перехода. В этом случае адрес перехода не будет влиять на правильность выполнения последующих команд программы. 5.7. Расчет коэффициента точности счисления Общие положения. Безопасность мореплавания и бесконфликтное ведение промысла напрямую зависят от точного и аккуратного счисления пути судна. Следует отметить, что в соответствии с требованиями ИМО точность счисления должна быть увязана с точностью обсерваций и макси- мально допустимым интервалом между ними с учетом ожидаемых по- грешностей от неточного учета приливо-отливных и других течений, погрешности в определении угла дрейфа. При этом считается, что счисле- ние ведется лагом и гирокомпасом, удовлетворяющим конвенционным требованиям. Для оценки точности текущего места используется допус- тимая радиальная погрешность места Мс (с вероятностью 95%) в зависи- мости от дистанции до опасности [3]. Оценка точности текущего места производится с учетом погрешности последней обсервации и погрешности счисления Для определения Л/с (Г) применяют два способа: статистический и априорный. Статистический способ определения Мс (Г) основан на совместном учете всех факторов, влияющих на точность счисления пути судна. Этот способ базируется на использовании теории случайных функций. Доказа- но, что средняя квадратическая погрешность счисления Мс (7) зависит от продолжительности плавания между надежными обсервациями и может быть определена по формулам: Л/с (Г) = 0,71Кс -г при t <2 ч; (5.3) Мс (t) = Ксу/Г при t > 2 ч. (5.4) В дальнейшем будем пользоваться полученными соотношениями. Кс вычисляют по'формуле
где Кс — коэффициент точности счисления; с,- — невязки счиспимых мест при продолжительности плавания 2 <5 ч; - продолжительность плавания на i-м курсе, ч; 1,13 - коэффициент, соответствующий вероят- ности 68,3%. Для того чтобы средняя квадратическая погрешность определения Кс не превышала определенного значения, необходимо определить коли- чество невязок, не менее чем и > 0,3 К2/т2кс, где п - максимальное коли- чество невязок; тК(, - СКП коэффициента счисления. Так, например, если Кс = 1,3, ткс = 0,1, то п > 0,3 • 1,3/0,01 = 51. Г/ = л f = и 11/2 I S Д?/ Е Г,-(и- 1)1 , где п — число невязок; Д,- — разность между соответствующими невяз- ками, вычисленными по формуле (5.4) и снятыми с карты. Относитель- ная погрешность [3, с. 102] Кс = 0,08 0,1 считается допустимой. Вполне понятно, что для получения значения коэффициента точности счисления с большей надежностью необходима систематизация и обра- ботка большого объема статистического материала как для данного судна, так и для группы однотипных судов. Расчет коэффициента точности счисления Кс представлен в Програм- ме 21, исходные данные - в табл. 5.4. Таблица 5.4, Распределение исходных данных и промежуточных результатов при расчете коэффициента точности счисления Параметр Обозначение Продолжительность плавания, ч Невязка, мили Суммарная продолжительность плавания, ч Сумма произведения чисел Константа Коэффициент точности счисления Х/1 Y/2 Примечание. Невязка (в милях) и продолжительность плавания на каж- дом курсе вводятся в регистры X и Y, а хранятся (до очередного цикла вычислений) в регистрах 1 и 2 соответственно. Это позволяет в случае сомнения в правильности ипода исходных данных произвести контроль после выполнения цикла вычислений. Программа 21 00 01 02 03 04 05 06 07 08 09 Х-П 1 fQ X—П2 П-Х 1 I- L- П-Х 1 П-Х 3 + Х-П 3 Ю 41 25 42 61 21 61 63 10 43 25 10 12 13 14 15 17 18 19 X П-Х4 + X—П4 С/П БП 00 П-Х4 П-ХЗ а. 64 10 44 50 51 00 64 63 13 20 21 22 П-Х 5 X С/П 65 12 50 75
Инструкция по использованию программы 21 1. Включить микрокалькулятор (положение пере- ключателя любое) 2. Перевести микрокалькулятор в режим „Програм- мирование” 3. Записать программу в память микрокалькулятора 4, Перевести микрокалькулятор в режим „Автома- тическая работа” 5. Произвести очистку регистров 3 и 4 (только пе- рец первым циклом вычислений и при повторном реше- нии задачи) 6. Записать константу 1,13 в регистр 5 7. Установить счетчик адресов по адресу 00 8. Произвести л циклов вычислений (расчет сумм, хранимых в регистрах 3 и 4) 9. Вычислить коэффициент точности счисления Примечания, 1, После выполнения вычислений Нажимаемые клавиши F ПРГ F АВТ СХХ-ПЗ Х-П4 1.13Х-П5 В/0 С/П шгс/п по п. 9 на индикаторе высветится значение Кс. 2. Операция ШГ нужна для того, чтобы „обезглавить” команду безусловного перехода, записанную по адресам 15 и 16. Вместо этой операции можно использовать другую: БП 17. Пример. Используя данные, приведенные ниже, рассчитать коэффициент точ- ности счисления Кс: Продолжительность плава- ния, ч.................. 2,5 2.9 3.2 2.0 2.0 Невязка, мили...........1,3 0,9 1,7 1.0 1.2 Решение. Будем полагать, что микрокалькулятор готов к решению задачи, поэтому описание примера приведено с п. 5 Инструкции. z Нажимаемые клавиши Индикация 5. Очистка регистров 3 и 4 СхХ-П 3 Х-П4 0. 6. Запись константы 1,13 в регистр 5 1.13Х-»П5 1.13 7. Установка счетчика адресов на ад- В/0 1.13 рес 00 1.3 8. Первый цикл вычислений 1.13В t 2.5 2.5 С/П 2.05 ... Второй 0.9 В t 6.9 2.9 2.9 С/П 3.588 ... Третий 1.7 В t 1.7 3.2 3.2 С/П 6.629 ... Четвертый 1 В t 2 2. С/П 8.04 Пятый 1.2 В t 1.2 2 2. С/П 9.74 9. Вычисление коэффициента точности С/П 8.735 -01 счисления Ответ: Кс = 0,87. 76
Примечание. Приведенный контрольный пример,ввиду небольшого коли- чества исходных данных (всего 5!) носит иллюстративный характер и может быть использован только при отладке программы. 5.8. Расчет поправки к высоте светила на наклонение зрительного луча При исправлении высоты светила, измеренной относительно частей корпуса другого („вспомогательного”) судна, вместо поправки за накло- нение видимого горизонта вводят поправку за наклонение зрительного луча, которая вычисляется по формуле Ыг=- [0,425 + 1,856(е - H)/S] , где ДЛ — поправка к высоте за наклонение зрительного луча, дуг.мин; 5 - расстояние до „вспомогательного” судна, мили; е — возвышение глаза наблюдателя над уровнем моря, м; Н - возвышение „вспомога- тельной точки” над уровнем моря, м. Ниже представлены табл. 5.5 с исходными данными и Программа 22 расчета поправок высоты, в которой предусмотрен расчет дистанции через приращение, задаваемое пользователем. Таблица 5.5. Размещение исходных данных и промежуточных результатов Параметр Обозначение Символ регистра Констант* 0.42 0 Начальное значение дистанции до вспомогательною S, 1 судна, мили Приращение дистанции, мили 6S 2 Текущее значение дистанции, мили 1 Константа 1.856 3 Разность возвышений, м е-Н 4 Поправка высоты за наклонение зрительною луча. ЛИ X дуг. мин Программа 22 00 01 02 03 04 05 06 07 08 09 1I-X2 П-Х1 + Х-П1 П—Х0 X П-ХЗ П-Х4 х П-Х1 62 61 10 41 60 12 63 64 12' 61 10 11 12 13 14 15 + + /-/ С/П БП 00 13 10 0L 50 51 00 77
Инструкция по работе с программой 22 1. Включить микрокалькулятор 2. Перейти в режим „Программирование” 3, Занести программу 4. Перейти в режим „Автоматическая работа” 5. Занести константы в соответствующие регистры Нажимаемые клавиши F ПРГ F АВТ памяти: 5.1. 0,42 - в регистр 0 5.2. Число, равное 6S, - в регистр 2 5.3. 1,856 - в регистр 3 6. Выбрать начальное значение расстояния So и за- нести число S„ - 6S в регистр 1 7. Выбрать первое значение е — Н и занести в реги- стр 4 (если число меньше нуля, то нажать клавишу /-/) 8, Установить счетчик адресов на нулевой адрес 9. Вычислить значения АЛ для всего диапазона рас- стояний 10. Повторить п. 6 11. Выбрать второе значение е- Ни занести в реги- стр 4 (если е — Н < 0, нажать клавишу /—/) 12. Вычислить значения АЛ для всего диапазона рас- стояний: 0,42 X-ПО 6SX-П2 1,856 Х-ПЗ So -6SX —П 1 (е-Н) (/-/) Х-П4 В/0 С/П So -8S Х-П 1 (е-„)(/-/) Х-П 4 остальные циклы В/0 С/П С/П Примечания, 1. Подробно этот вопрос изложен в работе ( 31. 2. Выполняя пп. 9 и 10, вычислить значения поправки высоты для всего диапа- зона разности возвышений и составить таблицу. Пример. Дано: е - И = ± 0,5 м; S = 1,5 * 6,0 миль; SS = 0,5 мили. Требуется рассчитать значения поправки высоты. Ниже (табл. 5.6) представлен фрагмент таблицы поправок высоты. Таблица 5.6, Поправка высоты за наклонение зрительного луча +0,5 -1,25 -1,30-1,42 -1,56—1,73—1,91 -2,10 - 2,29 - 2,48 -2,67 -2,87 0,0 -0,63 -0,84 -1,05 -1,26 -1,47 -1,68 -1,89 -2,10 -2,31 -2,52 -2,73 -0,5 -0,11 -0,38-0,68 -0,95 -1,2 -1,45 -1,67 ........... Теперь, приняв шаг изменения аргумента е - Н равным 0,5 м, выполним пп. 8 . . . 11 Инструкции: а) занесем е — Н = 0,5 в регистр 4; б) занесем So - SS = = 1.0 мили в регистр 1; в) произведем вычисления, нажимая клавиши В/0 С/П в каждом цикле. 78
5.9. Расчет элементов высотной линии положения с использованием суточных эфемерид Солнца и точки Овна 5.9.1. Общие положения. Известно, что в пределах суток Гринвичский часовой угол и склонение Солнца изменяются пропорционально интерва- лу времени. Часовой угол точки Овна изменяется равномерно за любой промежуток времени. Поэтому целесообразно при использовании про- граммируемого микрокалькулятора вычислять часовые углы Солнца и точки Овна, а также склонение Солнца автоматически по формулам t°M = t?P + Тгр Д/г®₽ ± Му; = ^fp + Тгр Д//р ± Х^.; 5й = 5®+ ГГрД6°; где /„, /„, tu - местные часовые углы Солнца, точки Овна и звезды соответственно; /®р> /f₽ - Гринвичские часовые углы Солнца и точки Овна соответственно на начало суток по Гринвичскому времени; Д/°р, Д/гр — приращения Гринвичского часового угла Солнца и точки Овна соответственно в градусах за один час; Ггр - Гринвичское время наблю- дения, ч; 5'г —, табличное значение склонения Солнца в градусах на Ггр = = 0 ч; Д6® - приращение склонения Солнца в градусах за один час Грин- вичского времени; X - счислимая долгота места наблюдателя, град. Для расчета Гринвичского времени в часах на момент наблюдения использована формула где ДГ?р, ДТ^р, ДТгр — доли Гринвичского времени наблюдений, выра- женные в часах, минутах и секундах соответственно. Расчет приращений координат. Для решения задачи в соответствии с методикой, изложенной выше, необходимо иметь значения приращений Гринвичского часового угла и склонения Солнца за один час Гринвичс- кого времени. В общем случае приращения вычисляют по формулам = ----- (,t?P2 - /?р>. + 360°); 24 Д5° = ------ (61 -6Т), 24 79
где /‘гр?, *rpi - Гринвичские часовые углы Солнца на 0 ч Гринвичского времени на последующие и текущие сутки соответственно, выбранные из МАЕ; 5®, 6® - склонение Солнца на 0 ч Гринвичского времени на по- следующие сутки соответственно. Для расчета приращений координат можно использовать подпрограм- му (адреса 00 ... 03 включительно) Программы 23, служащей для пере- вода угла, выраженного в минутах и градусах, в градусную меру. Однако исследования показали, что решение этой задачи можно упростить на ос- нове следующих предпосылок. Так как между значениями часовых углов Т®р2 и г?р| разность невелика, то Д1®р можно представить в следующем виде: Д^р = 24 -------+ 360е 60 где Д/'2, Д/\ - части часовых углов в дуговых минутах за вычетом гра- дусов. Так, например, для вычисления приращения часового угла на 01.05.88 г. будем иметь (см. Приложение 1) Д/’2 = 44,9'; Д/2 = 43,2'. Приращения Д6° вычисляем по формуле где Д62, Д5\ — части склонения, выраженные в дуговых минутах, за вы- четом целого числа градусов. Однако, если Д62 < ДЙ!, то к Д62 добавляется 60'. Так, например, для определения приращения Д5 на 04.05.88 г. (см. Приложение 1) будем иметь: Д62 = 75,4'; Д6] = 58,1'. Ниже представлена программа для расчета приращений координат часового угла и склонения. Программа 23 00 01 02 03 04 05 П-Х0 * П-Х1 + 2 11 60 13 61 10 02 06 07 08 4 т С/П 04 13 50 Размещение констант и переменных, а также результата показано ниже: Константы и переменные ..... 60 360(0) Дг'г г) Дг^Дб,) Регистры ................... 0 1 Y X Примечание. Параметры, приведенные в скобках, 80
Инструкция по расчету приращений координат 1. Включить МК. 2. Перевести МК в режим „Программирование” нажатием клавиш F ПРГ. 3. Набрать программу 23. 4. Перевести МК в режим „Автоматическая работа” нажатием клавиш F АВТ. 5. Занести константы 60 и 360 в регистры 0 и 1 соответственно. 6. Произвести отладку программы (решают контрольный пример), используя данные, представленные ниже: Нажимаемые клавиши Индикация 15.00.1319 В/0 С/П Если результат отличается от приведенного (табличного), приступают к поиску и исправлению ошибки. 7. Рассчитать Дг®р так, как это указано в и. 6. При этом имеется в виду, что значения t®p (на Ггр = 0 ч) на необходимое количество суток предварительно выбраны из МАЕ и занесены в таблицу, аналогично таблице в приложении 1. 8. После расчета Д?®р переходят к вычислению Д6?Предварительно производят очистку регистра 1 нажатием клавиш Сх X -> П 1. В остальном задача решается так же, как и при расчете Дг?р. Примечание. Знаки 6 и Д8 зависят от соотношения <р и 6s: если у и 8 раз- ноименны, то 8 имеет знак минус; если и 8 одноименны, то 6 имеет знак плюс; если б2 > 61, то знак Д6 совпадает со знаком 6®; еслиб, < 8то наоборот. В приложении 1 представлены координаты Солнца (îРи 5®) на 1988 г. на период с 1 по 21 мая. Рекомендуется в период предрейсовой подготовки заранее составить таблицу приращений координат. Расчет элементов высотной линии положения Солнца. В Программе 24 дня расчета элементов высотной линии положения использованы фор- мулы: sin/1 с — 'Pc sin5 + cos i/)c cos 5 cos tM; sinS — sin<pcsinAc cos A'c = ------------------- cos i/>c cosЛс где hc - счислимая высота; Ac - счислимый азимут; >/>c широта; 5 - склонение светила; tM - местный часовой угол. счислимая Программа 24 00 01 02 03 04 05 06 07 08 09 В t п-»ха X П ->ХЪ + П-*Х9 + х->п: П^Х 0Е 6Г 12 10 69 10 43 25 6С 10 11 12 13 14 15 16 17 18 19 X П-Ха + Х-’П 2 F 7 П -*Х 1 X П -X 1 F 8 12 6- 10 42 1С 61 1С 12 61 1С 20 21 22 23 24 25 26 27 28 29 81
X с/п 12 50 Пример. 12.05.87 г. наблюдали Солнце.м>с = 53° 10.0'N; Хс = 35°40.0’W; Ггр = - 14 ч 50 м 40 с; h0 = 54° 23.0'. Требуется рассчитать элементы высотной линии по- ложения. Решение. 1. Включить микрокалькулятор 2. Перейти в режим „Программирование” 3, Записать Программу 24 в оперативную память 4. Перейти в режим „Автоматическая работа” 5. Записать Дг в регистр d 6. Записать Д8 в регастр с Нажимаемые клавиши 7. Записать Тгр (на Тгр = 0 ч) в регистр В 8. Записать 8 (на Ггр = 0 ч) в регистр а 9. Записать ~рс в регистр 1 10. Записать Хс в регистр 9 11. Обсервованную высоту Ло записать в регистр 7 F ПРГ F АВТ 15.000208 Х-П d 1.048611 ВП02 /-/Х-П с 180-547 К+ Х-П Ь Х-П а (17.946667) 53.10 К+ Х-П 1 (53.166667) 35.40 К+ /-/Х-П9 (-35.666667) 54.23 К+ X - П 7 (54.383332) 12. Гринвичское время наблюдений перевести в 14.5040 К — в регистр 3 (14.84442) 13. Пустить программу на счет В/0 С/П (-8.74 -01) Примечания. 1. После завершения цикла вычислений на индикаторе высве- тится перенос, равный -0,87', т. е. -0,9*. Остальные величины хранятся в регистрах, указанных ниже: Параметры...................... 6 Регистры....................... 2 2. Для удобства решения задачи необходимо иметь трафаретку для записи значений исходных данных. Расчет элементов высотной линии положения звезды. Пример. Дано: 16.09.87 г. <рс = 22°27.0*^ Хс = 54°06.0'W; Тгр Измерена высота звезды Капелла (табл. 5.7). = 8 ч 36 мин. Исходные данные распределены следующим образом: Параметры..................... Дг Д6 Ггр 8 Регистры....................... d С b f 5? 82
Таблица 5.7. Исходные данные для расчета элементов высотной линии положения звезд на 16.09.87 г. г?р = 354° 24.6' (на Ггр = 0-ч) Звезда т* 6* л. Ггр Капелла 28 Г 03.3' 45° 59.2’ 65“ 1 l.O* 8 ч 34 мин 25 с Процион 245 20.3 5 15.7 43 16.5 8 35 49 Хамаль 328 22.6 23 24.4 54 15.7 8 36 51 Ригель 281 30.8 8 12.6 58 19.5 8 38 50 Дополнительное условие: микрокалькулятор использовался для обработки дневных наблюдений, т. е. для расчета элементов высотной, линии положения Солн- ца и расчета азимута для определения поправки компаса. 1. Дг ’ занести в регистр d 2. Произвести очистку регистра с 3. гт (на Тгр = 0 ч) занести в регистр 8 4. Счислимую широту записать в регистр 1 Нажимаемые клавиши 15.041042 Х-П d 354.246 К + 22.27 К На этом предварительная работа заканчивается и элементов высотной линии положения первой звезды. 5. Вычислить Гринвичский часовой угол звезды Ка- пелла на Тгр = 0 ч и записать в регистр Ь можно перейти к расчету 281.033 К + (281.055) П-Х8 + (635.465) Х-П b 45.592 К + (45.986666) Х-П а 65.11 К + (65.183332) Х-П 7 8.3425 К-В (8.573611) В/0 С/П (6.7 - 02) П-Х5 (16.08) П-ХЗ (710.32) Поскольку гм = 710.32° , т. е. 9,68° Е, вторая буква азимута будет иметь наиме- нование Е. 6. Склонение звезды записать в регистр а 7. Обсервованную высоту записать в регистр 7 8. Гринвичское время наблюдений (в часах и долях часа) записать на индикатор 9. Пустить программу на счет с адреса 00. (После завершения цикла вычислений высвечивается перенос) 10. Извлечь из регистра 5 счислимый азимут 11. Извлечь из регистра 3 местный часовой угол Теперь приступают к расчету второй и последующих линий положе- ния. Для этой цели выполняют действия, изложенные в пп. 6 ... 12. Пример. Используя данные табл. 5.7 и предушего примера, рассчитать элементы высотной линии положения звезды Процион. Решение. Нажимаемые клавиши 1. Вычислить Гринвичский часовой угол звезды Про- 245.203 К + цион на 0 ч Гринвичского времени (в градусах и долях (245.33832) П—Х8 + градуса) и записать в регистр В (599.74832) X — П b 83
2. Склонение звезды записать в регистр а 3. Обсервованную высоту (в градусах и долях граду- са) записать в регистр 7 4. Гринвичское время наблюдения (в часах и долях часа) записать в регистр X (индикатор) 5. Пустить программу на счет с адреса 00 (после за- вершения цикла вычислений на индикаторе высвечивает- ся перенос в милях) 6. Извлечь из регистра 5 Ас 7. Извлечь из регистра 3 местный часовой угол (в градусах и долях градуса). Так как практический часовой восточный гм = 45° ,0Е, следовательно, вторая буква ази- мута будет Е. Ответ: п = 1,1'; Лс = 104,6° NE. 5.157 К + (5.2616666) Х-Па 43.165 К + (43.275) Х-П 7 8.3549 К (8.5969442) В/0 С/П (1.06) П-Х5 (104.62) П-ХЗ (647.95 ...) Пример. По данным табл. 5.7 и предыдущей задачи рассчитать элементы высот- ной линии положения звезды Хамаль. 1. Вычислить Гринвичский часовой угол звезды Ха- маль на Ггр = 0 ч в градусах и долях градуса и занести в регистр Ь 2. Склонение звезды (в градусах и долях градуса) за- писать в регистр а Нажимаемые клавиши 328.226 К + (328.37666) П-Х8 + (682.78666) X - П b 23.244 К + (23.406666) 54.517 К + (54.861666) 8.3651 К - (8.6141666) В/0 С/П (8.45 -02) (80.8) 3. Обсервованную высоту (в градусах и долях гра- дуса) записать в регистр 7 4. Гранвичское время наблюдений (в часах и долях часа) записать в регистр X (индикатор) 5. Пустить программу на счет с адреса 00 (после за- вершения вычислений на индикаторе высвечивается пе- ренос высотной линии положения в милях) 6. Извлечь из регистра 5 счислимый азимут (в граду- сах и долях градуса) 7. Извлечь из регистра 5 счислимый местный часовой угол светила. Так как практический часовой угол запад- ный гм = 38.25°, то вторая буква азимута имеет наиме- нование W. Ответ: п = 0,1'; Ас = 80,8° NW. Пример. По данным предыдущего примера рассчитать элементы высотной линии положения звезды Ригель. (758.25) 1. Вычислить Гринвичский часовой угол звезды Ри- гель на 0 ч Гринвичского времени (в градусах и долях градуса) и записать в регистр Ь 2. Склонение звезды (в градусах и долях градуса) записать в регистр а Так как склонение разноименно с широтой, изме- нить знак на минус Нажимаемые клавиши 281.308 К+ (281.51332) П-Х8 + П-ХВ 8.126 К + /-/ (-8.21) Х-П а 84
3. Обсервованную высоту (в градусах и долях граду- 58.195 К + са) занести в регистр 7 (48.325) X — П 7 4. Гринвичское время наблюдений (в часах и долях 8.3850 К -• часа) записать в регистр X (8.647222) 5. Пустить программу на счет с адреса 00. (После за- В/0 С/П вершения цикла вычислений на индикаторе высвечивает- (-7.9 -02) ся перенос в милях) 6. Извлечь из регистра 5 счислимый азимут (в граду- П — X 5 ах и долях градуса) (164.57) 7. Извлечь из регистра 3 местный часовой угол све- П — X 3 тила. Так как практически часовой угол восточный гм = (711.886) = 8,1° Е, то вторая буква азимута будет иметь наименова- ние Е. Ответ: п = -0,1’; Ас= 164.6° NE. Для приведения линий положения к одному зениту предлагается дополнительная программа, составленная на основе формулы Дй'г = vAT/Cos (Aj - ПУ), где Дй'г - поправка переноса в минутах дуги (милях); v — скорость суд- на, уз;. АГ, - интервал времени приведения в долях часа (Д7) = То - 7)); А/ - азимут светила в круговом счете; ПУ — путевой угол судна; i = = 1,2 (и — 1), где п — количество линий положения; То — момент време- ни, к которому приводятся линии положения. Порядок размещения исходных данных, результат и Программа 25 расчета поправки приведены ниже: Параметры..................... v ПУ Г, А/ Т/ Дйг Внкстры....................... 2 3 4 Y X X Программа 25 52 53 55 56 57 58 59 60 61 П "Х4 Х-Пс F П-ХЗ F8 П-Хе 14 11 4Е 25 63 11 1Г 6Е 62 63 64 65 X П-Х 2 X С/П 12 62 12 50 Пример. Дано; v = 4,5 уз; ПУ = 223°. Требуется привести к одному зениту пинии положения предыдущего примера. Решение. После набора и отладки программы занесем исходные данные для первой линии положения. 1. Скорость движения судна в узлах 2. Путевой угол в градусах - в регистр 3 3. Го в часах - в регистр 4 4. Азимут (в круговом счете) первой линии положе- нии записывают в регистр Y Нажимаемые клавиши 4.5 Х-П2 223 Х-ПЗ 8.3850 К — (8.6472222) Х — П4 16.1 В 1 85
5. 7', в долях часа - в регистр X 6. Пустить программу на счет с адреса 52 Ответ: bhz = —6,3'. 8.3425 К « (8.573611) БП 52С1П (-0.30) После этого переходят к расчету поправки для второй линии положения. Теперь необходимо выполнить только пп. 4 ... 6. 1. Азимут (в круговом счете) второй линии положе- ния записать в регистр Y 2. Момент наблюдений (в долях часа) записать в ре- гистр X 3. Пустить программу на счет с адреса 52 Ответ: ДЛ2 = 0,1'. Нажимаемые клавиши 104. В1 8.3549 К -* (8.5969442) БП52С/П (-*0,11) И, наконец, рассчитать поправку для третьей линии положения. Нажимаемые клавиши 1. Азимут (в круговом счете) записать в регистр Н 279.1 В t 2. Время наблюдений (в долях часа) записать в реги- 8.3651 К-» стр X (8.6141666) 3. Пустить программу на счет с адреса 52 БП 52 С/П Ответ: Ahz = -0,1'. 5.10. Расчет обсервованных координат по наблюдениям двух светил Определение обсервованных координат по двум высотным линиям положения возможно как при разновременных наблюдениях (например, по высотным линиям положения Солнца), так и при одновременных (по наблюдениям звезд, с учетом приведения высот к одному зениту). Рассмотрим эту задачу только применительно к наблюдениям Солнца. В любом случае для расчета переноса и азимута используется одна и та же программа. Расчет обсервованных координат по двум высотным линиям положе- ния Солнца. При определении места судна по Солнцу вторые наблюдения производят после того, как азймут Солнца изменится на значение не ме- нее 30°. К этому времени, в общем случае, координаты судна значительно изменяются и, следовательно, для расчета элементов второй линии положе- ния необходимо в регистр 1 занести новое значение широты, а в регистр 9 — новое значение счислимой долготы. После расчета элементов второй высотной линии положения можно приступить к расчету обсервованных координат. Сущность способа заключается в следующем. На основе обще- известных формул вычисляются поправки к счислимым координатам, а затем сами обсервованные координаты на момент которых наблюдений. Для этой цели используются часть программы с адреса 61 (Программа 26) и табл. 5.8. 86
Таблица 5.8. Распределение исходных данных при расчете поправок координат по элементам двух высотных линий положения Солнца Параметр Единица измерения Символ регистра Перенос первой высотной линии положения мили 2 Азимут первой высотной линии положения в круго- градусы 3 вом счете Перенос второй высотной линии положения мили 4 Азимут второй высотной линии положения градусы 5 Разность широт мили X Разность долгот минуты дуги X Примечание. После получения поправок обсервованные координаты широт считается положительной, если она к северу, а разность долгот — к востоку. Программа 26 — 61 62 •63 64 65 66 67 68 69 70 П-Х5 П-Х 3 — F 7 Х-П8 П-Х2 П-Х 5 F 7 X П-Х4 65 62 11 1С 48 62 65 1С 12 64 71 72 73 74 75 76 77 78 79' 80 П-Х 3 F 7 X — П-Х8 - 4- С/П П-Х4 П-Х 3 F 8 63 1С 12 68 13 50 64 63 1Г 81 82 '83 84 85 86 87 88 89 90 1 Х П-Х2 П-Х5 F 8 X — П-Х8 4 П-Х 1 F 8 12 62 65 1Г 12 1 1 68 13 61 1Г 91 92 С/П 13 50 Инструкция по использованию программы 26 После расчета элементов второй высотной линии положения необходимо: Нажимаемые клавиши 1. Занести перенос второй высотной линии положе- Д6, X—П4 ния в регистр 4 2. Перевести азимут второй линии положения в кру- П — X 5 говой счет и записать в регистр 5 , (360 — - если лс NW) X —П4 3. Занести перенос первой высотной линии положе- Ыц X —П2 ния в регистр 2 4. Занести азимут первой линии положения (в круге- X — П 3 вом счете) в регистр 3 5. Рассчитать поправку широты ' С/П (БП 61 С/П) 6. Рассчитать поправку долготы С/П (БП 78 С/П) 87
Пример. 12.05.87 г. Судно следует курсом ГКК = 298° (+1,3°), скорость по ла- гу гл = 18 уз, отсчет лага ОЛ, = 25,40, поправка лага ДЛ = +2,6%, течение и дрейф не учитываем. В Ггр =114 10 мин 47 с, ОЛ = 85.18, ДЛ =+2,6% рассчитали счислимые координаты ¥>с= 53°46,9'N; Ас = 31° 15.9*W. Находясь по счислению в точке с коорди- натами у>с = 53°16.0'N; Х<; = 35°46.0"W в Тгр = 8 ч 21 мин 31 с получили обсерво- ванную высоту Солнца Ло = 14е38.9'. Используя микрокалькулятор и Программу 23, рассчитали элемент высотной линии положения: ДЛ, = -1,7'; Лс, =79.2°NE. В Тгр = И ч 10 мин 47 с, находясь по счислению в точке с координатами у>с = = 53°46.0*N; = 37° 15.9'W, получили вторую обсервованную высоту = = 37°58.0’ и рассчитали элементы второй высотной линии положения: ДЛ3 = -1,6 ; АС1 =115.1° NE. Требуется рассчитать поправки к координатам и обсервованные координаты на момент вторых наблюдений. Решение. Будем полагать, что часть программы, служащая для расчета поправок широты и долготы, набрана в память микрокалькулятора, на индикаторе находится значение переноса второй высотной линии положения (счетчик адресов указывает адрес следующей команды - 61), а азимут второй высотной линии поло- жения находится в регистре 5. Отметим, что так как оба азимута имеют наименова- ние NE, то их значения соответствуют азимуту в круговом счете. Следовательно, никакого перевода делать не надо. Ниже приводится последовательность решения задачи: Операции Н а жи маемые Индикация 1. Запись переноса второй линии положе- клавиши (-1,6) -1,6 ния в регистр 4 2. Перевод азимута второй линии положе- Х-П4 (П-Х5 (А, ния в круговой счет и запись в регистр 5 (в 360 360 данном примере эта операция не нужна) Х-+П5) 360 -А,) 3. Запись переноса первой высотной ли- 1,7 1.7 нии положения в регистр 2 Х-П2 -1.7 4. Запись азимута в круговом счете пер- 79.2 79.2 вой линии положения в регистр 3 х->п 3 79.2 5. Пуск программы (расчет поправки широты) с/п 5.49 -02 6. Расчет минут обсервованной широты 46 46.0 7. Пуск программы (расчет поправки С/П 46.05 -2.94 ... долготы) 8. Расчет минут обсервованной долготы 15.9/-/ -15,9 Ответ: = 53° 46.0* N; Хо = 37° 18.8* W. Примечания. 1. Счислимые координаты на момент -18.84 ... вторых наблюдений могут быть сняты с карты или вычислены аналитически. 2. Для решения рассматриваемой задачи целесообразно использовать програм- му, содержащуюся в блоке расширения памяти. 5.11. Работа с энергонезависимой памятью Общие сведения. У микрокалькулятора „Электроника МК-52” имеет- ся запоминающее устройство ППЗУ, служащее для хранения (использо- вания) записанных программ в объеме до 512 команд. Обращение к ППЗУ производится с помощью клавиш A t и t Ф по адресу, набираемому на индикаторе [8]. 88
Адрес обращения к ППЗУ (число) состоит из незначащей информации W (одна любая цифра, кроме нуля), начального адреса (0000 . , . 1023); f числа шагов; незначащей информации. Так, например, если необходимо обратиться к начальному адресу И 0023 для использования 20 шагов программы, то адрес обращения может иметь вид: 40000207 10"4. Здесь цифра 4 является незначащей, но обяза- И тельной; число 7 10~4 представляет собой незначащую, необязательную информацию. Обращение к ППЗУ осуществляется в следующем порядке: устанавли- вают режимы работы ППЗУ переключателем С/3,СЧ (сгирание, запись, считывание); набирают на индикаторе адрес (не менее семи цифр); последовательно нажимают клавиши A t и t 4. При нажатии этих кла- ,1 виш индицируется признак обращения к ППЗУ (знак минус во всех разрядах), сигнализирующий о том, что в это время нажимать другие В клавиши нельзя. Использование ППЗУ в режиме записи. Информация считывается I из памяти микрокалькулятора в ППЗУ последовательными пачками по семь шагов и размещается там так, что первый считанный шаг в первой пачке записывается в ППЗУ на 13-м и 14-м местах, считая от начального адреса, в адресе обращения к ППЗУ. Аналогичный способ укладки инфор- мации в накопитель ППЗУ используется для последующих доз информа- ции. После записи информации в ППЗУ программная память микрокаль- кулятора становится чистой. При записи информации в ППЗУ число шагов в адресе обращения к ППЗУ должно соответствовать числу шагов записываемой программы и fti.itb кратным семи. JL'iи «описи программы в ППЗУ необходимо: BK ihxmih микрокалькулятор; nrpciiiH н режим ..llpoi раммирование*'; юнее г и программу в программную память микрокалькулятора; не ре ключа! ель С/З/СЧ установить в положение 3 (запись); переключатель Д/П установить в положение П (запись из програм- мной памяти); перейги в режим „Автоматическая работа”; набрать на индикаторе адрес обращения (число, состоящее из набора I цифр нс менее семи); нажать клавишу АI, затем 11; переключатель С/З/СЧ установить в положение СЧ (во избежание Порчи информации); выключить микрокалькулятор. Рассмотрим подготовительную работу по записи программы в ППЗУ на примере программы, содержащей 21 шаг. Будем полагать, что запись В программы производится с адреса 00. Тогда адрес обращения к ППЗУ ж будет 1000021. Для определения конечного адреса ^'воспользуемся " формулой Nt = 2N - 1. Получим М = 2-21 - 1 = 41. Следует помнить, 89
что в данном случае участок ППЗУ с адреса 00 по 41 должен быть очищен, так как в противном случае новая информация наложится на старую, что недопустимо. Использование ППЗУ в режиме стирания. Как отмечалось, этот режим используется в том случае, когда на место старой информации нужно записать новую. Информация в ППЗУ стирается построчно. Начало и ко- нец стираемой информации определяется адресом обращения к ППЗУ. Для стирания информации по всей строке достаточно в адресе обращения к ППЗУ указать хотя бы один адрес четырехразрядного слова, располо- женного на данной строке. Максимальное количество стираемой информа- ции за один прием составляет 14 строк. Отметим, что для очистки ППЗУ с адреса 00 по 41 необходимо: установить микрокалькулятор в режим „Автоматическая работа”; установить переключатель С/З/СЧ в положение С (стирание); установить переключатель Д/П в положение П; набрать на индикаторе адрес обращения к ППЗУ, т. е. число 1000021; нажать клавишу At, затем, после прекращения признака обращения к ППЗУ, клавишу 11. На этом работа по чистке указанного объема инфор- мации будет завершена. Для очистки всего накопителя ППЗУ необходимо использовать новый начальный адрес, т. е. один из адресов на следующий неочищенной строке, в каждом последующем цикле. Подробные сведения об этом приведены в руководстве по эксплуатации микрокалькулятора. Использование ППЗУ в режиме считывания. Для считывания програм- мы из ППЗУ необходимо: включить микрокалькулятор; установить переключатель С/З/СЧ/ в положение СЧ (считывание); . установить переключатель Д/П в положение П (считывание в програм- мную память микрокалькулятора); набрать на индикаторе адрес обращения к ППЗУ (для нашего примера число 1000021); нажать клавишу A t, а после прекращения йндицирования признака обращения к ППЗУ — клавишу t J-; перейти в режим „Программирование” (после прекращения индици- рования признака обращения к ППЗУ) и сличить программу, записан- ную в программную память микрокалькулятора с табличной, имеющейся в пособии. Считанная программа, как это делается всегда, используется для ре- шения задач в режиме „Автоматическая работа” с соответствующими исходными данными. Ниже приведены сводные данные для записи программ в ППЗУ с на- чального адреса 0000: 90
Длина про- граммы (шаги) Требуемая длина Начальный адрес ячейки Адрес обра- щения к ППЗУ Конечный ячейки № строки памяти 62 63 0000 1000063 0125 8 82 84 0128 1012584 0292 19 55 56 0304 1030456 0415 26 91 91 0417 1041791 Длина записываемой в ППЗУ (и считываемой из ППЗУ) программы в шагах должна быть кратна семи. Начальный адрес ячейки последующей программы определяется по конечному адресу ячейки предыдущей про- граммы; он соответствует номеру первой ячейки на последующей строке памяти накопителя. Так, например, конечный адрес-первой программы равен 0125 и расположен н строке 8. Следовательно, начальный адрес программы 2 равен 0128 (номер первой ячейки строки 9 накопителя). При другой комбинации гирлянды программ начальные адреса ячеек второй и последующих программ будут иными. Пример. Требуется записать в оперативную память микрокалькулятора „Элект- роника МК-52" Программу 3. Р е in е н и с. Устанавливаем переключатель С/З/СЧ/ в положение СЧ (считыва- ние). Устанавливаем переключатель Д/П в положение П (считывание в программную память), набираем на индикаторе адрес обращения к ППЗУ - 1030456; нажимаем клавиши AI и 11. После прекращения признака обращения к ППЗУ (тире во всех разрядах микрокалькулятора) программа будет переписана в ОЗУ для дальнейшего использования. 5.12. Краткие сведения о блоке расширения памяти БРП-2 Блок расширения памяти (БРП) микрокалькуляторов „Электрони- ка МК-52” и „Электроника-Астро” предназначен для хранения программ астронавигационных задач, значений средних экваториальных координат 35 ярких звезд, а также программы навигационных задач. Имеется спе- циальная программа для проверки исправности микрокалькулятора. Информация, хранящаяся в БРП-2, считывается в ОЗУ. Соединение БРП-2 с микрокалькулятором производится при отключенном питании (габа- риты блока не более 44X58,5 X 18,5, масса блока не более 0,025 кг). В перечне астронавигационных программ содержатся следующие: вычис- ление Гринвичского часового угла точки Овна и экваториальных коорди- нат Солнца; вычисление счислимых высот и азимутов; вычисление эле- ментов высотной линии положения; вычисление обсервованных коорди- нат корабля по двум и более линиям положения и оценка его точности; исправление измеренных высот светил и другие программы. В перечне навигационных задач содержатся программы: аналитическое счисление 1 (вычисление координат места);. аналитическое счисление II (вычисление курса и расстояния) и другие программы. 91
Решение задачи происходит в 3 этапа: считывание программы задачи из ПЗУ БРП-2; подготовка к решению задачи; решение задачи. Считывание программы задачи производится с помощью адреса обращения, так же как из ППЗУ. В данном случае переключатель БРП-2 устанавливается в положение 1 (астронавигационные задачи) или в поло- жение II (координаты звезд, навигационные задачи). При этом, если считывается программа, то переключатель микрокалькулятора Д/П, как известно, ставится в положение II, а если считываются данные - то в положение Д. Подготовка к решению задачи выполняется при последовательном нажатии клавиш В/0 С/П. Решение задачи происходит в диалоговом режиме: после появления на индикаторе диалогового символа (число от 1 до 22) необходимо ввести на индикатор числовое значение запрашиваемого параметра и нажать клавишу С/П. Переключатель микрокалькулятора Р/ГРД/Г должен находиться в положении Г. В руководстве по эксплуатации БРП-2 [2] приведены подробные све- дения о правилах эксплуатации БРП-2 совместно с микрокалькулятором „Электроника МК-52”, а также контрольные примеры на весь перечень задач. В качестве иллюстрации рассмотрим пример расчета Гринвичского часового угла точки Овна и экваториальных координат Солнца на момент 7р = 23 ч 59 мин 59 с 01.05, 87 г Операции Нажимаемые Индикация клавиши 1. Включи гь МК 0. 2. Ввод начального адреса задачи 100С098 -1000098 3. Считывание программы из ПЗУ БРП-2 В/0 С/П 1 4. Ввод даты на запрос (1) 010587 С/П 2. 5. Ввод Тгр на запрос (цифра 2) 23.5959 С/П -1019698 6. Промежуточные операции В/0 С/П -1039298 7. Расчет rj.p В/0 С/П 219° 22.3' 8. Промежуточные операции С/П -1058898 9. Промежуточные операции В/0 С/П -1078498 .10. Расчет 6 е В/0 С/П 15° 09.2'N 11. Расчет а ° С/П 38° 38.7' 12. Расчет г гр С/П 180° 43.6' 1 3. Промежуточные операции -1098056. 14. Расчет А ’ t I 15.897 ... В/0 С/П (15,9') 15. Изменение знака R°, если измерения /-/ -15.897 производились по верхнему краю 16. Расчет адреса перехода на задачу (вы- С/П 3180677 числение счислимых высот и азимутов светил) Ответ: ггр = 180° 43.6'; 6 = 15° 09.2’N. 92
Полученные результаты соответствуют значениям, выбранным из МЛЕ, на Ггр = 0 ч 02.05.87 г. Детальные сведения об использовании блока расширения памяти приведены в пособии [2]. 5.13. Краткие сведения о блоке расширения • памяти БРП-3 БРП-3 предназначен дня постоянного хранения пакета 60 прикладных программ для использования с помощью микрокалькулятора „Электро- ника-52” [3]. В их числе: вычисление корней квадратного и кубического уровней; вычисление определителей 3-го и 4-го порядков; решение системы (2-4) линейных алгебраических уравнений; выполнение опера- ций с матрицами; численное решение обыкновенного дифференциального уравнения (I ... 4 порядков) и другие. Габаритные размеры блока нс более 58.5 X 42.3 X 10. Масса блока нс более 0,035 кг. Информация в постоянное запоминающее устройство (ПЗУ) записана в процессе изго- товления. Питание ПЗУ и управление его работой осуществляются от микрокалькулятора Порядок решения задачи такой же, как и при исполь- зовании БРП-2 И качестве иллюстрации рассмотрим пример вычисления корней кцад|штног(1 уравнении ил*2 + hx + c - 0. I. Набран. анк мим к ПЗУ-1000035, нажать клавиши А I и t I. » Паалть клавиши В/0 С/П. 2 и I pel метра 5 II р к м г ч л и и и. I. Время решения - около 7 с. .’.После завершения процесса вычислений на дисплее высвечивается значение дискриминанта. Если оно положительно, то, следовательно, корни действительные. 1лли Же дискриминант имеет отрицательное значение, то в регистрах 5 и 6 будут па ходиться действительная и мнимая части комплексных корней соответственно. Подробные сведения об использовании БРП-3 приведены в пособии | 3]. Таким об- разом, наличие блоков расширения памяти существенно повышает возможности микрокалькулятора „Электроника МК-52”.
Приложение I по Гринвичскому времени на май 1988 г. Дата Точка Овна (рр Солнце Г ®р 01.05 219° 08.2' 180°43.2' 15°04.8'N 02.05 220 07.4 44.9 22.8 03.05 221 06.5 46.5 40.6 04.05 222 05.6 48.0 58.1 05.05 223 04.8 49.4 16 15.4 06.05 224 03.9 50.6 32.3 07.05 225 03.1 51.6 50.0 08.05 226 02.2 52.6 17 05.5 09.05 227 01.3 ’ 53.4 21,6 10.05 228 00.5 54.0 37.5 11.05 228 59,6 54.6 53.0 12.05 229 58.8 180 54.9 18 08.3 13.05 230 57.9 55.2 23.2 14.05 231 57.0 55.3 37.8 15.05 232 56.2 55.2 52.2 16.05 233 55.3 55.0 19 06.2 17.05 234 54.5 54.7 19.9 18.05 235 53.6 54.2 33.2 19.05 236 52.7 53.6 46.2 20.05 237 51.9 52.8 58.9 21.05 238 51.0 51.9 20 11.3 Приложение 2 Плечи остойчивости формы относительно основной судна типа БМРТ пр. 394 Vo d 0° 10° 20° 30° 40° so° 60° 70° 80° 90° zm zc. 2000 3,60 0 1,08 2,20 3,30 4,34 5,28 6,21 6,83 7,06 6,94 6,22 1,98 2200 3,82 0 1,08 2,19 3,30 4,35 5,32 6,24 6,85 7,05 6,92 6,14 2,12 2400 4,15 0 1,08 2,18 3,29 4,38 5,38 5,38 6,85 7,02 7,00 6,10 2,28 2600 4,41 0 1,08 2,17 3,29 4,40 5,42 6,28 6,84 6,99 6,87 6,12 2,44 2800 4,63 0 1,08 2,17 3.29 4,42 5,45 6,29 6,80 6,95 6,85 6,14 2,60 3000 4,91 0 1,08 2,17 3,30 4,44 5,47 6,30 6,76 6,93 6,80 6,18 2,72 3200 5,22 0 1,10 2,18 3,30 4,45 5,47 6,28 6,73 6,88 6,76 6,22 2,88 3400 5,41 0 1.12 2,20 3,33 4.46 5,47 6,24 6,67 6,84 6,73 6,28 3,02 3600 5,62 0 1,12 2,22 3,34 4,47 5,47 6,20 6,62 6,78 6,67 6,33 3,15 3800 5,91 0 1,12 2,23 3,35 4,47 5,47 6,17 6,56 6,74 6,63 6,36 3,29 4000 6,09 0 1,13 2,23 3,36 4,48 5,47 6,14 6,53 6,70 6,58 6,42 3,42 94
Приложение 2 продолжение Плечи остойчивости формы относительно основной судна типа СРТМ-5О2Э V„ d zc. zm 10° 20° 30° 40° 50° 60° 70° 80° 90° 750 2,97 1.77 4,83 0,846 1,687 2,472 3,141 3.732 4.172 4.437 4,475 4,320 800 3,10 1,86 4,80 0,843 1,684 2,465 3,142 3,731 4,156 4,407 4,450 4,303 850 3,23 1,92 4,76 0,840 1,680 2,458 3,142 3.729 4.140 4,373 4.426 4.284 900 3,37 2,0 4,75 0,836 1,676 2,450 3,143 3,728 4.122 4,337 4,404 4,269 950 3,54 2,10 4,73 0,836 1,672 2,445 3,144 3,722 4,104 4,313 4,369 4,252 1000 3,67 2,18 4.76 0,836 1,668 2,439 3,146 3,716 4,289 4,335 4,335 4,235 1050 3,80 2.26 4,74 0,836 1,662 2,434 3,148 3,710 4,072 4,263 4,300 4.220 1100 3,93 2,33 4,77 0,838 1,657 2,430 3,145 3,709 4,054 4,239 4,263 4,202 1150 4,04 2,40 4,78 0,840 1,652 2,426 3,142 3,707 4,036 4,219 4,234 4,184 1200 4,17 4,79 0,842 1,647 2,421 3,138 3,706 4.023 4.204 4,204 4,204 4,166 Плечи остойчивости формы относительно основной судна типа СРТМ-502М V. zm 'CO 10° 20° 30° 40° 50° 60° 70° 80° 550 2.54 4,55 .53 0.306 1.599 2,343 2,904 3,314 3.664 3.951 4,026 600 2.71 4.53 .63 0,797 1.588 2.329 2,877 3.296 3,652 3,928 4.015 650 2,82 4.52 ,72 0,790 1,581 2,311 2.853 3,285 3,646 3,903 4,001 700 3,05 4.49 .82 0,786 1.576 2,290 2.831 3,275 3,640 3,880 3,984 750 3.20 4.48 0,783 1.573 2,266 2,812 3,268 3,634 3,861 3,964 800 3.34 4,46 .0 0.783 1.570 2.242 2,795 3,263 3,627 3,844 3.938 850 3,45 4.42 .05 0.784 1.563 2,219 2,780 3,259 3,620 3,829 3.916 900 3,77 4,47 .21 0,789 1,5 38 2,177 2,755 3,251 3,602 3,806 3,882 1000 3.88 4,47 0,792 2,159 2,744 3,244 3,593 3.796 3,872 1050 4.08 4,49 .39 0.795 1,505 2,142 2,734 3,236 3,583 3,787 3.864 1100 4.17 4.55 .45 0,785 1,490 2,127 2,725 3,225 3,572 3,778 3,857 Приложение 3 Средние квадратические погрешности N ---------------------------------------------------------- /0 т„ т 1 2 3 4 5 6 7 8 I 14,3 13,5 12,9 14,0 13,8 14,2 13,6 13,9 13.78 0,16 0,45 2 104,5 105,4 105,5 104.8 106.1 106,8 106,3 105,4 105,6 0,27 0,77 3 4.3 4.7 3,9 4,0 4,3 4,7 3,6 3,8 4,16 0.14 0,41 4 12,0 12,7 12,3 ‘ 12,5 13,0 12.4 12,2 12,1 12.4 0,12 0,33 5 4 4 5,0 4,0 4.8 4.6 3,9 4,6 4,5 4,55 0.11 0,32 6 12.4 12,0 13.0 12.7 12.1 11,9 12,7 12,6 12,42 0.14 0,39 20,5 20,3 20,8 20.9 20,4 20,7 20.9 20.1 20,58 0,10 0,30 8 19,7 20,1 18,9 19.4 19,0 20,1 20,2 20,7 19,76 0,22 0,63 9 220,1 221,6 222.0 220.9 220,4 220,9 221,1 220,8 220.98 0.22 0,61 10 14,2 13,9 14,3 13.8 14,4 14,0 14,4 13,8 14.10 0,09 0,26 11 12.8 13,4 12,6 14,0 12,1 12.4 12.7 11.9 12.74 0.24 0,68 12 4.3 2,9 3,8 2 6 3.5 3.9 2,7 4.0 03.46 0,23 0,65 13 3.6 3,7 3.9 3.4 3.1 3,0 3,1 3.6 03,42 0.12 0,33 95
Продолжение прил. 3 Приложение 4 Персональный микрокомпьютер „Электроника МК-85” Общие сведения. „Электроника МК-85” [8] представляет собой кар- манный персональный компьютер (ПК), используемый в режиме обыч- ного калькулятора, так и в режиме „Вычисление по программам”. Причем в отличие от микрокалькулятора „Электроника-52” в нем используется одна из версий языка Бейсик (рис. 10). Алфавит Бейсика включает в себя: прописные и строчные буквы латинского алфавита от А до Z; прописные и строчные буквы русского алфавита; знаки арифметических операций: +. —, X, / - сложения, вычитания, умножения, деления соответственно и другие специальные знаки. Рис, 10. Внешний вид компьютера „Электроника МК-85” 96
На панели ПК расположены 54 клавиши и жидкокристаллический дисплей, на который можно выводить одновременно до 12 символов. Дня реализации совмещенных функций имеются префиксные кла- виши S и F, а также клавиша MODE, которая обеспечивает выбор режима работы ПК. Коды режимов работы ПК указаны на передней панели, в их числе: пуск программы; запись программы в память, про- верка и редактирование ее; ввод в режим построчной отладки программы; вывод из режима построчной отладки программы и др. ПК оперирует угловыми величинами, выраженными в градусах, радианах или градах. Для этой цели используются клавиши DEG, RAD, GRA. Приведем некоторые основные характеристики ПК МК-85 Число разрядов мантиссы Число разрядов порядка........................ -, Число адресуемых регистров памяти . . (26 с возможностью расширения до 1 78) 10 10 (Р0 ... Р9) Объем энергонезависимой памяти, байт Число одновременно хранимых про- грамм (файлов)..................... Максимальное число шагов программы 1221 (около 150 строк программы на языке Бейсик) Бейсик 4 элемента СЦ 0,18 или блок пита- Язык программирования Потребляемая мощность, Вт ...... . Диапазон рабочих температур, ° С . . . . । абаритные ра тмсры. мм........... Мам и. Г. не более................. не более 0.02 5 ... 40 13 X 165 X 73 150 Испош. тование принципа „бегущей строки позволяет записывать строки программы длиной до 63 символов. При помощи клавиш перемещения курсора ->и «- возможен просмотр всей строки, ПК имеет примерно такие же встроенные функции, как и микрокалькулятор „Электроника МК-52". Точность вычислений состав- ляет ± 1 в десятом знаке. ПК в режиме калькулятора (прямой режим) позволяет выполнять арифметические операции следующим образом: с помощью клавиатуры вводятся числа и арифметические операции (в по- рядке их следования); нажимается клавиша EXE; снимается результат. Пример. Дано выражение N= 42 - 1 1 + 4 X 7/5. Требуется найти чис- ло N. Решение. 1. Нажимаем последовательно клавиши MODE и ф. На дисплее появится надпись RUN. 2. Вводим правую часть выражения 42-11 +4X7/5. 3. Нажимаем клавишу EXE. 97
4. Снимаем результат: N= 36.6. Использование ПК в режиме ,,Программирование”. Как указывалось, ПК имеет 1221 шаг энергонезависимой памяти. Для определения коли- чества шагов программы следует пользоваться правилом: номер строки программы занимает 2 шага программы; программируемая команда, функциональная команда, символ, нажатие клавиши EXE (аналог коман- ды - возврат каретки) соответствует одному шагу. Например, строка вида l<p0FORl = lTON занимает 8 шагов в памяти. Так как макси- мальная длина строки составляет 63 символа, то при записи в нее несколь- ких операторов можно сэкономить память на к шагов, где к - число сэкономленных адресов. ПК оперирует как с числовыми, так и строковыми переменными длиной до семи символов (букв,цифр, знаков). Бейсик ПК МК-85 содер- жит 21 оператор. Рассмотрим некоторые из них. Оператор ввода INPUT предназначен для ввода переменных с клавиа- туры во время выполнения программы. Встретив оператор INPUT, Бейсик прекращает выполнение программы и печатает вопрос, а если в программе записано наименование переменной, то напечатает и наименование. После ввода "мелового значения переменной и нажатия клавиш ЕХЕ Бейсик продолжит вычисления. Если программой предусмотрен ввод числовых данных для нескольких переменных, то ввод значений каждой переменной сопровождается нажатием клавиши ЕХЕ. После ввода значений для после- , дней переменной начнется процесс вычислений. Пример. 2 ф INPITT Z Z. Оператор IF THEN предназначен для передачи управления по условию. Если условие выполняется, то происходит переход на строку, номер которой записан после слова THEN. Если же условие не выполняется, то команда перехода игнорируется и Бейсик выполняет команду, записан- ную в следующей строке. Зф INPUT L 4 0IFL = 0 THEN 11 ф 5ф ................ 11 0 END В данном примере, если L = 0, Бейсик после строки 3 ф прекратит естественный порядок выполнения программы и перейдет к выполнению команды, записанной в строке 110, если же условие не будет выполнять- ся, то Бейсик будет выполнять команду по адресу 5 0. Если управление необходимо передать подпрограмме, записанной в другом файле, то вместо адреса перехода записывается знак =$= и номер файла перехода. Оператор GOSUB используется для безусловного перехода к под- программе в данном файле или в другом файле. Пример. 1 0 GOSUB 4 0 7 0 GOSUB #Р1 Здесь 4 0 - первая строка подпрограммы, на которую следует перейти, а Р1 — файл, в котором находится подпрограмма. 98
Оператор RETURN записывается в конце подпрограммы. Он обеспе- чивает возврат к строке, следующей за той, в которой записан оператор перехода к подпрограмме. Оператор окончания программы END используется для обозначе- ния конца программы. Оператор выполнения программы RUN предназначен для пуска про- граммы на счет. Если пуск программы производится со строки, то после набора оператора RUN набирается номер строки и нажимается клави- ша ЕХЕ Операторы стирания. Оператор CLEAR предназначен для стирания программы только в одном файле, номер которого указан перед опера- тором CLEAR. Для уничтожения всех программ используется команда CLEARA. Ввод и использование программы для решения задач. Вводы програм- мы в компьютер (ПК) осуществляются клавишами в следующем по- рядке: I. ПК переводится в режим записи WRT нажатием клавиш MODE 1. Выбираете» нужный файл нажатием клавиш S и цифровой от ф до 9 (для выбо|Ш файла <> I’ ф до I’ 9) 3. Построчно в поди । с и программа. Переход на следующую строку прои 1но/иися нажатием клавиши ЕХЕ. При этом на индикаторе высвечи- вается количество оставшихся шагов памяти, что очень удобно для конт- 4 После МПнеи П|х>|раммы ПК переводится в режим вычислений нажатием клавиш MODE Ф (.чедует помнить, что если при выборе файла вместо номера на дисп- лее 11ОН11И1СП нюбраженне ромба, это означает, что данный файл занят, । с птм laitHcaira программа. Если решено вместо нее записать другую, г» старую программу нужно стереть. Для этого нужно набрать: номер файла, «тем слово ( LI AR и. наконец, нажать клавишу ЕХЕ. Выполнение программы рассмотрим исходя из двух ситуаций: а) про- грамма только что записана в какой-либо файл и ПК находится в режиме записи; б) ПК находится в режиме вычислений. В первом случае поступают так: производят очистку индикатора нажатием клавиши АС; набирают номер нужного файла нажатием кла- виш S и номера файла; производят пуск программы нажатием клавиш SRUN (или набирают слово RUN, а затем нажимают клавишу ЕХЕ). Во втором случае, т.е. когда ПК находится в режиме вычислений, поступают следующим образом: набирают номер нужного файла, нажимают клавишу ЕХЕ. Если про- грамма выполняется не с начального адреса, то поступают так: вводят команду RUN и номер адреса; нажимают клавишу ЕХЕ. Если при отладке программы (решении контрольного примера) будет обнаружена ошибка (ошибки), то имеется возможность отредакти- ровать программу. Для просмотра программы используется команда 99
LIST. Однако редактирование программы возможно только лишь в режи- ме записи программы (WRT). Даже коаткое описание возможностей микрокомпьютера позволяет заключить, что он обладает большими возможностями по сравнению с микрокалькуляторами, а наличие энергонезависимой памяти, позволяю- щей записать более 1000 команд, в сочетании с алгоритмическим языком Бейсик и малыми габаритами микрокомпьютера делают его надежным помощником судового специалиста.
СПИСОК ЛИТЕРАТУРЫ \. Астронавигационный альманах на 1986 -1990 гг. Л.: ГУНиО МО, 1987. 2. Блок расширения памяти БРП-2 „Электроника-астро”.: Руководство по эксплуатации, 277 с. 3. Блок расширения памяти „Электроника БРП-3”: Руководство по эксплуа- тации. 196 с, 4. Дьяконов В. П. Справочник по расчетам на микрокалькуляторах. 2-е изд., испр. М.: Наука, 1986. 224 с. 5. Брмакон В. И. Методы оперативною контроля остойчивости судов флота рыбной промышленности: Учсбн. пособ. Калинин! рад ВИНК Минрыбхоза СССР, 1987. М -.. п.//'.чакон В. И. Промысловое судовождение расчеты. М.: АгрОпромиздат. 1984.14’1- I. Mi I'>->цI:и обработки паинишионных измерений с оценкой точности. Л.- ГУИНОМО, 1985.52 с. И А1к«. ptuoi,1ьн у, hi гор ..Эчектроннм МК-52": Руководство по эксплуатации. Him I. I IHIIГИ, 198 Ч. Г..х..нти1гг1>л по iip.iK 1Н-1Г, кому Применению программируемых микрокаль-
ОГЛАВЛЕНИЕ Предисловие.............................................................. 3 Основные обозначения..................................................... 4 Глава 1. Использование микрокалькулятора в непрограммируемом варианте . . 5 1.1. Назначение клавиш....................................... 5 1.2. Выполнение одноместных операций......................... 5 1.3. Занесение чисел в регистры памяти и извлечение из них... 6 1.4. Выполнение двухместных операций......................... 2 Глава 2. Использование микрокалькулятора в режиме „Вычисление по про- граммам!'............................................................. 9 2.1. Общие положения........................................ 9 2.2. Программирование задачи................................. Ю 2.3. Вводи редактирование программы.........................1' 2.4. Занесение исходных данных и решение задачи.............12 2.5. Оптимизация программы командами переходов и организации циклов....................................................... 12 Глава 3. Решение основных задач, связанных с контролем плавучести и остой- чивости судна........................................................ 18 3.1. Расчет координат центра массы судна.................... 18 3.2. Расчет координат центра массы судна при приеме (снятии) и пере- мещении грузов................................................23 3.3. Расчет плеч статической остойчивости судна..............23 3.4. Расчеты, связанные с построением универсальной диаграммы ста- тической остойчивости ....................................... 25 3.5. Расчет плеч динамической остойчивости...................30 3.6. Расчеты, связанные со снятием судна с мели (камня)......34 3.7. Расчет опрокидывающего момента по методике Регистра СССР ... 25 3.8. Расчет поперечной метацентрической высоты по измеренному периоду собственных бортовых колебаний судна................37 3.9. Построение универсальной диаграммы динамической остойчи- вости ..........................................................38 Глава 4. Расчет инерционно-тормозных характеристик судна в условиях про- мысла .....................................................42 4.1. Пассивное торможение..................................43 4.2. Активное торможение...................................46 4.3. Пассивное торможение судна при движении с тралом......4g 102
Глава 5. Использование микрокалькулятора для решения основных задач судовождения.................................................50 5.1. Определение места судна по импульсно-фазовым РНС........50 5.2. Ведение промыслового счисления......................... 54 5.3. Плавание по дуге большого круга......................... 5.4. Расчет девиации магнитного компаса.................... g.j 5.5. Расчет элементов движения встречного судна по данным радио- локатора .................................................... 64 5.6. Расчет средней квадратической погрешности результатов наблю-' деннй...........................................;.......... . 59 5.7. Расчет коэффициента точности счисления................. 74 5.8. Расчет поправки к высоте светила на наклонение зрительного луча......................................................... 77 5.9. Расчет элементов высотной линии положения с использованием ' суточных эфемерид Солнца и точки Овна...................... 7q 5.10. Расчет обсервованных координат по наблюдениям двух светил . . §6 5.11. Работа с энергонезависимой памятью..................... 88 5.12. Краткие сведения о блоке расширения памяти БРП-2....... 91 5.13. Краткие сведения о блоке расширения памяти БРП-3....... 93 Приложение 1. Значение экваториальных координат на начало суток по Гринвич- скому времени на май 1988 г.. .......................................... 94 Приложение 2. Плечи остойчивости формы относительно основной судна типа БМРТ пр. 394 ........................................................... 94 Приложение 3. Средние квадратические погрешности...................... 95 Приложение 4. Персональный микрокомпьютер „Электроника МК-85”......... 96 < пнеок литературы......................................................101