Влияние проектных характеристик судна на его мореходные и прочностные качества тема диссертации и автореферата по ВАК РФ 05.08.03, кандидат технических наук Суров, Олег Эдуардович

• Суров, Олег Эдуардович
• кандидат технических наук
• 2000

Оглавление диссертации кандидат технических наук Суров, Олег Эдуардович

Введение.

Глава 1. Анализ экспериментальных и теоретических исследований.

1.1. Аварии судов в штормовых условиях.

1.2. История развития теории качки и продольной прочности судов на волнении.

1.3. Экспериментальные исследования продольной качки и прочности.

1.4. Цели и задачи исследования.

Глава 2. Алгоритм и компьютерная программа для расчета продольной качки судна в линейной и нелинейной постановке.

2.1. Алгоритм линейного расчета продольной качки судна.

2.2. Алгоритм нелинейного расчета продольной качки судна.

2.3. Моделирование нетрадиционного волнения.

2.4. Компьютерная программа для расчета продольной качки и прочности судна на волнении.

2.5. Сопоставительные и тестовые расчеты.

Введение диссертации (часть автореферата) на тему «Влияние проектных характеристик судна на его мореходные и прочностные качества»

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

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

При оценке мореходности в первую очередь необходимо знание параметров качки. Основоположником теории продольной качки может по праву считаться А.Н. Крылов, который еще сто лет назад изложил свои исследования в этой области. В дальнейшем на основе теории А.Н. Крылова И.Г. Бубнов разработал практический метод динамического расчета прочности судов при продольной качке на волнении. Последующее развитие теории качки шло по трем основным направлениям: учет возмущений в поле давлений окружающей жидкости, вносимых присутствием и качкой корабля (гидродинамическая теория); учет конечности амплитуд качки и нелинейного характера зависимости действующих сил от кинематических параметров качки (нелинейная теория); учет нерегулярности морского волнения (вероятностная теория качки).

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

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

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

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

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

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

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

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

Для достижения поставленной цели решаются следующие задачи:

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

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

3. Численное теоретическое исследование влияния проектных (форма корпуса, коэффициенты полноты) и эксплуатационных (скорость хода, распределение нагрузки, курсовой угол) характеристик на мореходность и прочность.

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

5. Разработка методики и рекомендаций для проектирования судов, на основе выполненных исследований.

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

Настоящая диссертационная работа состоит из введения, пяти глав, заключения, списка литературы и трех приложений. Объем диссертации 128 страниц машинописного текста, 62 иллюстрации, 26 таблиц, библиография из 147 наименований.

Заключение диссертации по теме «Проектирование и конструкция судов», Суров, Олег Эдуардович

ЗАКЛЮЧЕНИЕ

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

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

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

1. Разработан комплекс компьютерных программ для расчетов продольной качки и волновых нагрузок, действующих на корпус судна при плавании произвольным курсом на регулярном и нерегулярном волнении, в линейной и нелинейной постановке, с различными скоростями хода. В программу включены различные формы волнового профиля, отличающиеся от синусоидального в том числе нерегулярное волнение, с целью исследования влияния нетрадиционного волнения на параметры мореходности и прочности судна. Комплекс позволяет формировать базу данных по теоретическим чертежам корпусов судов с дальнейшим представлением их в графической форме как на экране монитора, так и на бумажном носителе. В удобной форме организован диалог пользователя. Результаты при окончании расчетов записываются в отдельные файлы с возможностью их последующего просмотра, а также выводятся на экран монитора в графической форме. Сопоставительные и тестовые расчеты, выполненные по другим подобным программам, показали хорошее согласование результатов, что дает основания использовать комплекс для дальнейших численных исследований мореходности и прочности судна на волнении. Программа используется в учебном процессе ДВГТУ по дисциплинам "Методы построения теоретического чертежа" и "Динамика корабля", в курсовом и дипломном проектировании.

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

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

4. Исследовано влияние параметра (хг -х^/Ь, определяющего степень взаимного влияния вертикальной и килевой качки. Практически во всем исследованном диапазоне скоростей хода с ростом параметра (хг -х^/Ь (от отрицательных к положительным значениям) отмечается уменьшение значений

АЧХ продольной качки примерно на 10 %, но АЧХ волновых моментов увеличиваются, причем более значительно - на 20 - 25 %, а в отдельных случаях до 50 %. Такая закономерность прослеживается во всем реальном диапазоне изменения радиуса инерции масс.

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

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

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

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

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

181

Список литературы диссертационного исследования кандидат технических наук Суров, Олег Эдуардович, 2000 год

1. Антоненко C.B. К вопросу о влиянии продольного радиуса инерции масс на волновые изгибающие моменты. НТО Судпрома, вып. 131, Л.: Судостроение, 1969.-С. 81 -90.

2. Антоненко C.B., Суров О.Э. Влияние распределения нагрузки на величину волнового изгибающего момента. // Тр. междунар. конф. Кораблестроение и океанотехника. Проблемы и перспективы: ДВГТУ, Владивосток, 1998. С. 217-221.

3. Антоненко C.B., Суров О.Э. Исследование влияния несимметричного распределения нагрузки на волновые изгибающие моменты. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1999. С. 141 - 145.

4. Антоненко C.B., Суров О.Э. Использование модели ветровых волн для расчета поведения судна в море. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1999.-С. 135 140.

5. Антоненко C.B., Суров О.Э. Расчет мореходности и прочности судна в заданном режиме волнения. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1996. С. 52-56.

6. Антоненко С.В., Суров О.Э. Численное исследование прочности и качки судна на встречном волнении. // Тр. II междунар. конф по судостроению: ЦНИИ им. акад. А.Н. Крылова, СПб., 1998. Т 1. - С. 409 - 416.

7. Ашик В.В Проектирование судов. Л.: Судостроение, 1975. - 352 с. Ю.Бабаев H.H., Лентяков В.Г. Некоторые вопросы общей вибрации судов. - Л.: Судпромгиз, 1961. - 308 с.

8. П.Банди Б. Методы оптимизации. Вводный курс. -М.: Радио и связь.1988,- 128 с.

9. Барабанов Н.В. Конструкция корпуса морских судов. В двух томах. Том 1. Общие вопросы конструирования корпуса судна. СПб.: Судостроение, 1993.-304 с.

10. Барабанов Н.В. Конструкция корпуса морских судов. В двух томах. Том 2. Местная прочность и проектирование отдельных корпусных конструкций судна. СПб.: Судостроение, 1993. - 336 с.

11. М.Барабанов Н.В., Иванов H.A., Новиков В.В., Окишев В.А., Чибиряк И.М. Повреждения судовых конструкций. Л.: Судостроение, 1977. - 400 с.

12. Барабанов Н.В., Иванов H.A., Новиков В.В., Шемендюк Г.П. Повреждения и пути совершенствования судовых конструкций. Л.: Судостроение,1989,-254 с.

13. Барабанов Н.В., Кустов В.Н. Принципы, вызывающие отрыв носовой оконечности судов и ведущие к их гибели. // Сб. докладов междунар. конф., посвященной 100-летию ЦНИИ им. акад. А.Н. Крылова. СПб, 1994.-С. 227-233.

14. Басин A.M. Качка судов. -М.: Транспорт, 1969. 272 с.

15. Бельгова М. А. Изгибающие моменты для судов внутреннего плавания на волнении. Л.: Судостроение, 1966. - 208 с.

16. Беляк Ю.Л. Экспериментальное исследование прочности корпусов судов. Л.: Судостроение, 1964. - 230 с.

17. Благовещенский С. Н., Холодилин А. Н. Справочник по статике и динамике корабля. Л.: Судостроение, 1976. Т - 2, - 176 с.

18. Бойцов Г.В. К вопросу об определении внешних нагрузок при слемин-ге. // Судостроение, 1985. N 5, с.30.

19. Бойцов Г.В. Оптимизация судового корпусах учетом требований снижения его металлоемкости и трудовых затрат// Судостроение, 1984. N 3,

20. Бойцов Г.В., Палий О.М. Прочность и конструкция корпуса судов новых типов. Л.: Судостроение, 1979, 360 с.

21. Болотин В.В. Статистические методы в строительной механике. М.: Строииздат, 1965, 279 с.

22. Бородай И.К., Нецветаев Ю.А. Качка судов на морском волнении. -JL: Судостроение, 1969. 432 с.

23. Бородай И.К., Нецветаев Ю.А. Мореходность судов: Методы оценки. -Л.: Судостроение, 1982. ч

24. Бронников A.B. Особенности проектирования морских транспортных судов. Л.: Судостроение, 1971. - 328 с.

25. Бронский А.И. Корпусные конструкции судов промыслового флота. Л.: Судостроение, 1978, 198 с.

26. Бубнов И.Г. Строительная механика корабля. В 2 частях. ЧЛ, 1912. 4.2,1914.

27. Васильев Ф.П. Численные методы решения экстремальных задач. -М.: Наука, 1980.- 518 с.

28. Васильева Э.Ю. Влияние основных проектных характеристик на качку судна в условиях нерегулярного волнения. Проектирование и конструкции судов: Автореф. дис. канд. техн. наук. Калининград, 1992. - 13 с.

29. Ветер и волны в океанах и морях: Справочные данные. Л.: Транспорт, 1974.-С. 182-192.

30. Внешние силы, действующие на суда: ЦНИИ РУМБ. Л., 1976. - 120 с.

31. Вознесенский А.И. Вероятностный расчет перемещений, скоростей и ускорений различных точек корабля на нерегулярном морском волнении. Доклад к IX науч. техн. конф. НТО Судпрома по теории корабля. Л., 1958.

32. Галахов И.Н., Кандель Ф.Г., Фридлянский А.З. Исследование изгибающих и крутящих моментов в модельном эксперименте на регулярном волнении. Материалы по обмену опытом; - Тр. НТО Судпрома, - JI, 1972. -Вып. 189. - С. 119-130.

33. Герасимов A.B. Энергостатистическая теория нелинейной нерегулярной качки судна. Д.: Судостроение, 1979.

34. Герасимов A.B., Мореншильдт В.А., Нецветаев Ю.А. Разработка инженерной методики расчета продольной качки транспортных судов. В сб.: Экспериментальная гидродинамика судна (Тез. Докладов к VIII науч. - техн. конф.). - Л., НТО Судпрома, 1966. Вып. 80.

35. Екимов. В.В. Методы теории вероятностей в кораблестроении. Л.: Судостроение, 1970. - 272 с.

36. Екимов. В.В. Вероятностные методы в строительной механике корабля. Л.: Судостроение, 1966. - 328 с.

37. Ершов Н.Ф., Попов H.A. Прочность судовых конструкций при локальных динамических нагружениях. Л.: Судостроение, 1989. -200 с.

38. Ершов Н.Ф., Шахверди Г.Г. Метод конечных элементов в задачах гидродинамики и гидроупругости. Л: Судостроение, 1984. - 237 с.

39. Иванов H.A. Экспериментальное исследование экстремальных волновых нагрузок, действующих на судовые конструкции. // Сб. докладов между-нар. конф., посвященной 100-летию ЦНИИ им. акад. А.Н. Крылова. СПб, 1994. С.415-421.

40. Иванов H.A., Каленчук C.B. Анализ условий, сопровождающих катастрофические повреждения судовых конструкций. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1996. С. 38 - 44.

41. Ипатовцев Ю.Н., Короткин Я.И. Строительная механика и прочность корабля. Л.: Судостроение, 1991. - 288 с.

42. Картузова Т.А., Ростовцев Д.М. Оценка волновых и вибрацтонных изгибающих моментов, действующих на корпус судна в услових нерегулярного волнения. // Строительная механика и прочность судовых конструкций: сб. науч. Тр./Л., ЛКИ, 1981.-С. 41 48.

43. Кондриков Д.В., Четыркин Н.В. Применение статистического метода при оценке общей прочности судна по результатам эксперимента. -Тр.//ЦНИИМФ Л.: 1962. - Вып.41. - С.З - 22.

44. Корн Г., Корн Т. Справочник по математике. -М.: Наука, 1968.-720 с.

45. Королев В.И. Экспериментальные исследования величины изгибающих моментов судов на волнении. // Судостроение, 1960. №1.

46. Короткин Я.И., Рабинович О.Н., Ростовцев Д.М. Волновые нагрузки корпуса судна. Л.: Судостроение, 1987. - 236 с.

47. Короткин Я.И., Ростовцев Д.М., Сивере Н.Л. Прочность корабля. Л.: Судостроение, 1974. - 432 с.

48. Короткин. Я.И. Вопросы прочности морских транспортных судов. -Л.: Судостроение, 1965. 390 с.

49. Крылов А.Н. Избранные труды.-Л.: Академия наук СССР, 1958.-804 с.

50. Курдюмов А.А Вибрация корабля. Л.: Судпромгиз, 1961. - 320 с.

51. Лазарев В.Н., Юношева Н.В. Проектирование конструкций судового корпуса и основы прочности судов. - Л.: Судостроение, 1989. - 210 с.

52. Луговский В.В. Гидродинамика нелинейной качки судов. Л.: Судостроение, 1982. - 250 с.

53. Луговский В.В. Динамика корабля. Л.: Судостроение, 1976. - 199 с.

54. Лузянин A.A. Исследование особенностей продольной качки и дополнительного сопротивления морских составных судов. Теория корабля: Авто-реф. дис. канд. техн. наук. С.-Петербург, 1997. - 22 с.

55. Любушин Н.П. Экономическая эффективность проектных решений в судокорпусостроении. Л.: Судостроение, 1982.

56. Максимаджи А. И. Амплитудно-частотные характеристики волновых изгибающих моментов. // Тр. ЦНИИМФ. 1971. - Вып. 134. - С. 11-32.

57. Максимаджи А.И. Экономические факторы при нормировании продольной прочности корпусов судов// Судостроение, N 12, 1990. c.l 1.61 .Максимаджи А.И. Прочность морских транспортных судов. Д.: Судостроение, 1976. - 311 с.

58. Нецветаев Ю. А. Анализ результатов испытаний моделей для исследования волновых моментов в корпусах транспортных судов на встречном волнении. // Тр. ЦНИИ им. А. Н. Крылова. 1968. - Вып.245. - С.77-93.

59. Нецветаев Ю. А. Исследование распределения гидродинамических давлений по корпусу корабля при движении на встречном волнении. ЦНИИ им. А. Н. Крылова. Диссертация, 1958.

60. Ногид JIM. Остойчивость судна и его поведение на взволнованном море. Л.: Судостроение, 1967. - 244 с.

61. Ногид Л.М. Проектирование морских судов. Л.: Судостроение, 1976. - 208 с.

62. Ногид Л.М. Проектирование формы судна и построение теоретического чертежа. Л.: Судпромгиз, 1962.

63. Нормы прочности морских судов. Регистр СССР. 1991, 92 с.

64. Перельмутер A.C. Экспериментальное исследование влияния демпфирования на характеристики килевой качки. ПММ, т. X, 1946.

65. Постнов В.А. Численные методы расчета судовых конструкций. Л.: Судостроение, 1983. - 248 с.

66. Постнов В.А., Тарануха H.A. Оценка напряженно-деформированного состояния корпуса судна методом модуль-элементов. // Судостроение, 1983. №5, С. 5 - 8.

67. Правила классификации и постройки морских судов. Регистр Судоходства, СПб., 1995.

68. Прикладные задачи динамики судов на волнении. Под ред. И.К. Боро-дая Л.: Судостроение, 1989. - 220 с.73 .Продольная качка и слеминг лихтеров лихтеровозной системы. Черноморское ЦПКБ. Руководитель Я.М.Элис.-№ Гр 77030573-0десса, 1977.-20 с.

69. Прочность судов внутреннего плавания. Справочник. М.: Транспорт,1978.

70. Путов Н.Е. Проектирование конструкций корпуса морских судов. Часть 1. Л.: Судостроение, 1976. - 376 с.

71. Путов Н.Е. Проектирование конструкций корпуса морских судов. Часть 2. Л.: Судостроение, 1977. - 424 с.

72. Разработка алгоритмов и программ расчета изгибающих и крутящих моментов, действующих на корпуса судов на волнении. ЛКИ. Руководитель Д.М. Ростовцев. № Гр У 136881 - Л., 1976.- 178 с.

73. Разработка метода расчета сил и моментов, действующих на судно на регулярном волнении. НКИ. Руководитель Ю.В. Ремез. N Гр 75018950 - Николаев, 1976.- 79 с.

74. Разработка теоретических методов расчета волновых нагрузок. ЦНИИМФ Руководитель А.И. Максимаджи. № Гр 76050702 - Л., 1976,- 98 с.

75. Рахманин H.H. Актуальные проблемы мореходности судов. // Судостроение, 1991. N4, С. 3-5.

76. Ремез Ю.В. Качка корабля. Л.: Судостроение, 1983. - 327 с.

77. Ростовцев Д. M. К вопросу определения волнового изгибающего момента. // Тр. ЛКИ. 1960. - Вып.27.

78. Салькаев А.З. Метод определения гидродинамических коэффициентов в уравнениях продольной качки судов на волнении. // Судостроение, 1978.1. N12,-С. 8-13.

79. Семенов-Тян-Шанский В.В., Благовещенский С.Н., Холодилин А.Н. Качка корабля. Л.: Судостроение, 1969. - 392 с.

80. Сердюченко А.Н. Анализ аварий и катастроф судов в море. // Тр. между нар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1999. С. 155 - 160.

81. Сердюченко А.Н. Исследование волновых нагрузок на корпуса судов с учетом групповой структуры морского волнения. Строительная механика корабля: Автореф. дис. канд. техн. наук. Николаев, 1980. - 27 с.

82. Сердюченко А.Н. Нелинейные эффекты в расчетах гидродинамических нагрузок на судовые конструкции. Л.: ЦНИИ "Румб". - 1984. - Вып.40. -С.78-84.

83. Сердюченко А.Н., Федоренко Г.В. Расчет нелинейных гидродинамических нагрузок на судовые корпуса при качке на волнении. // Тр. междунар. конф. Кораблестроение и океанотехника. Проблемы и перспективы: ДВГТУ, Владивосток, 1998. С. 181-186.

84. Сердюченко А.Н. Обеспечение безопасности судна на волнении в катастрофических условиях. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1996. С. 16-21.

85. Слободян С.О. Исследование волновых нагрузок судов на регулярном волнении умеренной крутизны. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1996. С.30-37.

86. Слободян С.О. Оценка влияния некоторых факторов на экстремальные волновые нагрузки. // Тр. междунар. конф. Кораблестроение и океанотехника. Проблемы и перспективы: ДВГТУ, Владивосток, 1998. С. 386-390.

87. Справочник по вероятностным расчетам.-М.: Воениздат, 1970. 536 с.

88. Справочник по строительной механике корабля./ Бойцов Г.В., Палий О.М., Постнов В.А., Чувиковский B.C. В трех томах .- Л.: Судостроение, 1982.

89. Справочник по строительной механики корабля, т 3. Под ред. Ю.А. Шиманского. Л., Судпром гиз, 1960. - 800 с.

90. Справочник по теории корабля: В трех томах. Статика судов. Качка судов. Под ред. Я. И. Войткунского. Л.: Судостроение, 1985.

91. Суров О.Э. Влияние нагрузки судна на продольную качку на регулярном волнении. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1996. С.57-61.

92. Суров О.Э. Исследование волновых изгибающих моментов на регулярном косом волнении.// Материалы II международной конф. "Проблемы транспорта Дальнего Востока", Академия транспорта РФ, Владивосток, 1997. -С.83-84.

93. Суров О.Э., Корепанова Ю.А. Расчетное и экспериментальное исследование волновых изгибающих моментов. // Молодеж и научно-технический прогресс: Материалы конф. ДВГТУ, Владивосток, 1998. - С. 51 - 52.

94. Суслов В.П. О расчетах нагрузок, действующих на суда в экстремальных волновых условиях. // Строительная механика корабля: Сб. науч. тр. -Николаев, НКИ, 1986. С.3-11.

95. Суслов В.П. Об экстремальных волновых нагрузках, действующих на корпуса судов на морском волнении. // Строительная механика корабля: Сб. науч. тр НКИ, 1978. Вып. 136. - С. 15-22.

96. Суслов В.П., Паплаускас А.Р. Алгоритм и ЭВМ программа для расчета нелинейной продольной качки и волновых нагрузок на экстремальном волнении.// Строительная механика корабля: Сб. науч. тр. - Николаев, НКИ, 1990.-С.З-16.

97. Суслов В.П., Суслов C.B. Нелинейный метод расчета волновых нагрузок.// Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1996. С.8-15.

98. Суслов В.П., Суслов C.B. О вероятностно детерминированном методе определения волновых нагрузок, действующих на суда. // Строительная механика корабля: Сб. науч. тр. - Николаев, НКИ, 1994.

99. Ю8.Сухир Э.Л. Амплитудно-частотные характеристики продольной качки и волновых моментов для гибких судов. // Регистр СССР: Сб. науч. тр. Л., Транспорт, 1978. - Вып.8. - С.49- 67.

100. Ю9.Сухир Э.Л. Исследование изгибающих моментов при балластном плавании судов на волнении. Строительная механика и вибрация корабля: Ав-тореф. дис. канд. техн. наук. Одесса, 1967. - 14 с.

101. ПО.Фирсов Г.А. Развитие идей академика А.Н. Крылова в области теории и расчета кораблей на волнении. Тр. НТО Судпрома, 1957. Т. - 7, вып. 2.

102. Хаскинд М.Д. Гидродинамическая теория качки корабля. М.: Наука, 1973.- 327 с.

103. Химмельблау Д. Прикладное нелинейное программирование. М.: Мир, 1975.- 537 с.

104. Холодилин А.Н. Стабилизация судна на волнении. Л.: Судостроение, 1973. - 230 с.

105. Цындаря И.И. Прогиб корабля при нелинейной качке. // Тр. ВМАКВ им. акад. А.Н. Крылова. Л., 1951. Вып. 4, 5.

106. Четыркин Н.В. Об условиях расчета характеристик внешних сил при вероятностной оценке общей прочности корпуса судна. //Тр. ЦНИИМФ Л.: 1962. - Вып.41. - С.23- 42.

107. Пб.Чижиумов С.Д. Исследование динамики слеминга судна с помощью численных моделей. // Тр. междунар. конф. Кораблестроение и океанотехника. Проблемы и перспективы: ДВГТУ, Владивосток, 1998. С. 410-414.

108. Чижиумов С.Д. Численные алгоритмы определения присоединенных масс воды в задачах вибрации и качки судов. Там же С. 415-419.

109. Чувиковский Г.С. Исследование величин внешних сил, действующих на корпус судна в условиях удара о встречную волну. -Тр./НТО Судпрома, 1960,N351.

110. Шаров Я.Ф. Об ударе судна днищем о встречную волну. // Тр. НТО СП, 1960. -Вып.35.

111. Шмырев А.Н., Морендшильдт В.А., Ильина С.Г., Гольдин А.И. Успокоители качки судов. Д.: Судостроение, 1972. - 480 с.

112. Элис Я.М. Исследование гидродинамических характеристик уравнений качки судов. Теория корабля: Автореф. дис. канд. техн. наук. Одесса, 1970.- 22 с.

113. Юрков Н.Н. Некоторые вопросы мореходности судов на встречном волнении. Теория корабля: Автореф. дис. канд. техн. наук. Д., 1974. - 20 с.

114. Antonenko S.V., Korepanova Y.A., Surov O.E. Wave Load Ships Of Various Types On Irregular Waves. // Third International Students' Congress of the Asia-Pacific Region Countries. Congress materials, FESTU, Vladivostok, Russia, 1999.-PartII.-pp. 6-8.

115. Dalzell J.F., Chicco M.J. Wave Load in a Model of the SL-7 Container ship Running at Oblique Headings in Regular Waves. // SSC, 1973.

116. De-Does I. Experimental Determination of Bending Moment for Three Model of Different Fullness in Regular Waves. "Int. Shipbuilding Progress.", 1960. -N68, vol. 7,-pp. 161-180.

117. Fukuda J. Statistic Prediction of Ship Response. Society of Naval Architects of Japan, Second Symposium on Ship Manoeuverability, July 1969.

118. Joosen P., Wahab R., and Woortman J. Vertical Motion and Bending Moment in Regular Waves. "Int. Shipbuilding Progress", 1968. - January, N 16, -pp. 15-31.

119. Kusumoto H. et al. Ship Hull Structure Design by DLA (Dynamic Loading Approach), IHI Engineering Review,Vol.27, No.2, April 1994.

120. Lewis E. V. Ship. Model Test to Determine Bending Moment in Waves. SNAME, Vol. 62, 1954. pp. 431-490.

121. Lotveit M. and Miirer Chr. Further Model Test to Determine Wave Loads on a T-2 Tanker. "European Shipbuilding", 1964. N 3, - pp. 52-73.

122. Lotveit M. and Oth. Wave Loads on a T-2 Tanker Model. "European Shipbuilding", 1961. N 10, - pp. 31-53.

123. Moor D. Longitudinal Bending Moments on Models in Head Seas. TRINA, 1967.-pp. 117-166.

124. Newman J.N. A linearised Theory for the Motion of a Thin Ship in Regular Waves. Journal of Ship Res., 1961. Vol. 4, - pp. 10 - 15.

125. Peters A.S., Stoker J.J. The Motion of a Ship as a Floating Rigid Body in a Seaway. Com. Pure and Appl. Math., 1957. Vol. 10, - pp. 399 - 490.

126. Salvesen N., Tuck E.O. and Faltinsen O. Ship Motions and Sea Load. TSNAME, Vol. 78, 1970.

127. Surov O.E. The Program for Investigation of Vessel Longitudinal Motion and Wave. // Second International Student's Congress of the Asia-Pacific Region Countries. Abstracts, FESTU, Vladivostok, Russia, 1997. pp. 189 - 190.

128. Tomita Y., Tozawa S. Wave Load and High Tensile Steel. Society of Naval Architects of Japan, Symposium on Ship's Strength and Failure Prevention. Oct. 1993.

129. Vossers G. Behavior of Ship in Waves. Harlem, 1962.

130. Vossers G., Swaan W. A. and Rijiken H., Experiments with Series 60 Models in Waves, International Shipbuilding Progress, May, 1961. p. 201.

131. ОПИСАНИЕ ФУНКЦИЙ И ПОДПРОГРАММ

132. DECLARE FUNCTION WINDS (tITLES, LEFT%, RIGHT%, UP%, DOWN%, INK%, PAPER%) DECLARE FUNCTION SAVESCR% (NAMES) DECLARE FUNCTION RESTSCR% (NAMES)

133. DECLARE FUNCTION INTEGRAL.CEB! (X!(), Y!(), N.SHP!, KOL.TOCH!(), SUMMA!)

134. DECLARE SUB SAVE.FILE (Y!(), X!(), KOL.TOCH!(), NAMESS)

135. DECLARE SUB SPLAIN (X!(), Y!(), OSADKA!, N.SHP!, X!, Y!, KOL.TOCH!())

136. DECLARE SUB RED (X!(), Y!(), KOL.TOCH!())

137. DECLARE SUB OPEN.FILE (Y!(), X!(), KOL.TOCH!(), NAMESS)

138. DECLARE SUB RED.TOCH (X!(), Y!(), 1!, N.SHP!, KOL.TOCH!())

139. DECLARE SUB RED.SHP (X!(), Y!(), KOL.TOCH!(), N.SHP!)

140. DECLARE SUB raspred.nagruzki ()

141. DECLARE SUB zalivaemost ()

142. DECLARE SUB runl5 (lv!, SIG!, ka!)

143. DECLARE SUB RUN.NONLINE (ZETA, FI, ZETAD, FID, X)

144. DECLARE SUB FON1 (IN!, ik!, jn!, jk!)

145. DECLARE SUB message (S$, C$)1. DECLARE SUB shapkamenu ()

146. DECLARE SUB RAMKA (XN!, XK!, YN!, YK!, sim!, FOON!)

147. DECLARE SUB Spravka (CS$, DF$)1. DECLARE SUB oshibka (ty!)1. DECLARE SUB sapkatab ()

148. DECLARE SUB star (NAM!, R!)1. DECLARE SUB REDTOTAL ()1. DECLARE SUB neregulyr ()

149. DECLARE SUB prnscreen (wr!)1. DECLARE SUB SAVE ()

150. DECLARE SUB printdat (AS(), NAM!, R!)1. DECLARE SUB openfile ()1. DECLARE SUB Regulyr ()1. DECLARE SUB NEWFILE ()

151. DECLARE SUB message 1 (n$, S$)

152. DECLARE SUB qlstn (X!, Y!, ST$, B!)

153. DECLARE SUB submenu (NAM!, R!)1. DECLARE SUB UPMENU (NAM!)

154. DECLARE SUB MENU (X!, Y!, SH!, n!, S$(), NAM!)1. DECLARE SUB ZASTAV4 ()1. DECLARE SUB ZASTAV2 ()

155. ОПИСАНИЕ ПЕРЕМЕННЫХ И КОНСТАНТ

156. НАЧАЛО ГЛАВНОЙ ПОДПРОГРАММЫ1. H% = 201. CALL ZASTAV2 WIDTH 80, 5010000 CALL FONl(l, 48, 1, 80): PCOPY 0, 110001 CALL star (NAM, R)77 ty = ERR: CALL oshibka(ty) PCOPY 1,0 IF ty = 11 THEN

157. ПОДПРОГРАММЫ *************1. РАСЧЕТ КПМ И КД

158. SUB KPMKD (ka, ТТ(), ВВ(), ВАТАА(), NAP, BP, I, М33(), М330(), L33(), L330(), SIG) DIM BBB(4, 2, 2), A(4, 2, 2) BBB(0, 0

159. M33(I) = MUKI * BB(I)A 2 * 3.141593 * RO/8 + NAP * RO * BP M330(I) = MUI * BB(I) A 2 * 3.141593 * RO/8 + NAP * RO * BP L33(I) = LKYI * BB(I)A 2 * RO * SIGk / 4 L330(I) = LYI * BB(I)A 2 * RO * SIG / 41. END SUB

160. SUB neregulyr 'РАСЧЕТ НА НЕРЕГУЛЯРНОМ ВОЛНЕНИИ

161. BES = 2.5: ALS = .21 * BES

162. AS2 = ALS л 2 BES A 2: BS2 = ALS л 2 + BES A 2

163. MSZ = 0: MZET = 0: MFIM = 0: M.SI.M = 0: M.I.M = 0: MSZG = 0: MAZG = 0 f SIGBN = .5: SIGBKON = 2.5: stepsig = .075 ' DIM S.I.M(IOOO)

164. CATE ynl + 1, xnl + 2: PRINT " ВАРИАНТ I LOCATE ynl + 3, xnl + 2: PRINT "Интенсивность волнения в баллах " LOCATE ynl + 4, xnl + 2: PRINT "Обеспеченность волнения в % " LOCATE ynl + 5, xnl + 2: PRINT "Номер расчетного сечения " xn2=40:хк2 =79

165. CALL RAMKA(xn2, хк2, ynl, ykl, sim, fonn) COLOR sim, fonn

166. CATE ynl + 1, xn2 + 2: PRINT " ВАРИАНТ II

167. CATE ynl + 3, xn2 + 2: PRINT "Высота волны обеспеч. 3%, см "

168. CATE ynl + 4, xn2 + 2: PRINT "Обеспеченность волнения в % "

169. CATE ynl + 5, xn2 + 2: PRINT "Номер расчетного сечения "хпЗ = 20: хкЗ = 60: упЗ = 30: укЗ = 38

170. CALL RAMKA(xn3, хкЗ, упЗ, укЗ, sim, fonn)1. COLOR sim, fonn

171. CATE YN + 1, XN + 2: PRINT " Для расчета продольной качки судна"

172. CATE YN + 2, XN + 2: PRINT "на нерегулярном волнении необходимо"

173. CATE YN + 3, XN + 2: PRINT "ввести дополнительные данные по трем,"

174. CATE YN + 4, XN + 2: PRINT " ниже перечисленным вариантам "

175. XN = 20: ХК = 60: YN = 42: YK = 46

176. CALL RAMKA(XN, ХК, YN, YK, sim + 16, 6)1. COLOR sim, 61444 LOCATE YN + 2, XN + 2: INPUT "ВВЕДИТЕ НОМЕР ВАРИАНТА ВВОДА ДАННЫХ ", variant.dat1. variant.dat <= 0 OR variant.dat >= 4 THEN 1444 IF variant.dat = 1 THEN COLOR sim, fonn

177. CATE ynl + 3, xnl + 34: INPUT"", BAL LOCATE ynl + 4, xnl + 34: INPUT"", obes

178. CATE ynl + 5, xnl + 34: INPUT "", sechenie

179. CATE ynl + 3, xn2 + 35: INPUT "", AH3: LOCATE ynl + 4, xn2 + 35: INPUT "", obes

180. CATE ynl + 5, xn2 + 35: INPUT "", sechenie: AH = AH3 /1.32: TV = 5.9 * (ah /100) л (1 /

181. TV = 3.3 * SQR(AH3 /100): END IF: IF variant.dat = 3 THEN: COLOR sim, fonn

182. CATE yn3 + 3, хпЗ + 35: INPUT "", AH3: LOCATE yn3 + 4, хпЗ + 35: INPUT "", TV

183. CATE yn3 + 5, хпЗ + 35: INPUT "", obes: LOCATE yn3 + 6, хпЗ + 35: INPUT "", sechenie

184. AH = AH3 /1.32: END IF: xnl1 = 1: xkl1 = 17: ynl 1 = 30: ykl 1 = 38

185. CALL RAMKA(xnl 1, xkl 1, ynl 1, ykl 1,14, 4): COLOR sim, 4

186. CATE ynl 1 + 2, xnl 1 + 2: PRINT " 1. СПЕКТР"

187. CATE ynl 1 + 4, xnl 1 + 2: PRINT "РЕКОМЕНДУЕМЫЙ"

188. CATE ynl 1 + 6, xnl 1 + 2: PRINT " 12 МКОБ"xnl2 = 63: xkl2 79: CALL RAMKA(xnl2, xkl2, ynl 1, ykll, 14, 4): COLOR sim, 4

189. CATE ynl 1 + 2, xnl2 + 2: PRINT " 2. СПЕКТР"

190. CATE ynl 1 + 4, xnl2 + 2: PRINT "РЕКОМЕНДУЕМЫЙ"

191. CATE ynl 1 + 6, xnl2 + 2: PRINT " 2 МКПС"1155 XN= 13: XK = 68: YN = 42: YK = 46

192. CATE 45, 20: PRINT "Шаг по частоте";" " ";

193. CATE 46, 18: PRINT "ОТНОСИТЕЛЬНАЯ ЧАСТОТА ОПРЕДЕЛЯЕТСЯ КАК : vL/Lv"

194. CATE 43, 60: INPUT "", SIGBN

195. CATE 44, 60: INPUT "", SIGBKON

196. CATE 45, 3 5: INPUT " ", stepsig1. SCREEN 12: WIDTH 80, 60

197. WINDOW (0, -0)-(1000, 1000)

198. CATE 1,1: PRINT " Спектр "

199. CATE 2, 1: PRINT " волнения "

200. CATE 1,20: PRINT" АМПЛИТУДНО ЧАСТОТНЫЕ ХАРАКТЕРИСТИКИ" LOCATE 2, 20: PRINT

201. CATE 3, 20: PRINT " вертикальной"1.CATE 4,20: PRINT " качки"

202. CATE 3, 40: PRINT " килевой"

203. CATE 4, 40: PRINT " качки"

204. CATE 3,60: PRINT " изгибающего"

205. CATE 4, 60: PRINT " момента"

206. CATE 5, 60: PRINT " для сечения "; sechenie

207. CATE 34, 20: PRINT " СПЕКТРЫ ПРОЦЕССОВ"1.CATE 35,20: PRINT

208. CATE 29,1: PRINT TAB(l); SIGBN; TAB(16); SIGBKON; TAB(22); SIGBN; TAB(36); SIGBKON; TAB(42); SIGBN; TAB(56); SIGBKON; TAB(62); SIGBN; TAB(75); SIGBKON;1.CATE 29, 51: PRINT ""

209. CATE 30, 28: PRINT "Относительная частота vL/Lv"

210. CATE 58, 1: PRINT TAB(22); SIGBN; TAB(36); SIGBKON; TAB(42); SIGBN; TAB(56); SIGBKON; TAB(62); SIGBN; TAB(75); SIGBKON; LOCATE 58, 65: PRINT "";

211. SZ / 300 > MSZ THEN MSZ SZ / 300: SIGB = SIGB + stepsig: WEND: SIGB = SIGBN

212. WHILE SIGB <= SIGBKON: SIG = SIGB * (2 * 3.141593 * 9.81 / LS) Л .5: GOSUB 20000

213. NE (20 + (SIGB stepsig / 5 - SIGBN) * 225 / (SIGBKON - SIGBN), SZ / MSZ + 550)-(20 +

214. SIGB + stepsig / 5 SIGBN) * 225 / (SIGBKON - SIGBN), 550), 15, BF

215. SIGB = SIGB + stepsig: WEND: COLOR 15

216. CATE 7, 5: PRINT "Smax ="; CLNG(MSZ * 30000) / 1001.CATE 8, 5: PRINT" мп.с"

217. COLOR 14: LOCATE 42,1: PRINT "ФОРМА СПЕКТРА"

218. CATE 45, 1: PRINT " ПОДХОДИТ ?": LOCATE 48, 1: INPUT "0-HET: 1-ДА ", QWW IF QWW = 0 THEN SCREEN 0: WIDTH 80, 50: A% = RESTSCR%("NEREG.MJG"): GOTO 1155 OPEN "nereg" FOR OUTPUT AS #4

219. FGG(CH) = ZETAM: FGGl(CH) = FIM: ACH.I.M(CH) = I.M(sechenie) / ZAimm = (GAMMA * LS A 2 * Y(KOL.TOCH(H% \ 2) 2, H% \ 2) * ZA * 2)imm = I.M(sechenie) * 100 / imm: PRINT #4, USING "###.####"; SIGB; ZETAM; FIM; immспектральная плотность изг.момента от (х)

220. S.I.M(CH) = ACH.I.M(CH) А 2 * SZ: S.I.M = S.I.M(CH) + S.I.M: pfim = FIM * ka

221. SZG(CH) = ZETAM Л 2 * SZ: SZG = SZG + SZG(CH)

222. SZ / 300 > MSZ THEN MSZ = SZ / 300

223. ZETAM / 300 > MZET THEN MZET = ZETAM / 300

224. FIM / 300 > MFIM THEN MFIM = FIM / 300

225. S.I.M(CH) / 300 > M.SI.M THEN M.SI.M = S.I.M(CH) / 300

226. ACH.I.M(CH) / 300 > M.I.M THEN M.I.M = ACH.I.M(CH) / 300

227. SZG(CH) / 300 > MSZG THEN MSZG = SZG(CH) / 300

228. SAG(CH) / 300 > MAZG THEN MAZG SAG(CH) / 300

229. CH.= CH + 1: SIGB = SIGB + stepsig: WEND

230. COLOR 15: LOCATE 7, 5: PRINT "Smax ="; CLNG(MSZ * 30000) / 100 LOCATE 8, 5: PRINT " мп.с"

231. CATE 7, 25: PRINT "Zmax ="; CLNG(MZET * 30000) / 100 LOCATE 7, 45: PRINT "Fmax CLNG(MFIM * 30000) / 100 LOCATE 7, 65: PRINT "Mmax ="; CLNG(M.I.M * 30000) / 1001.CATE 8, 65: PRINT" т"

232. CATE 37, 25: PRINT "Szmax ="; CLNG(MSZG * 30000) / 100 LOCATE 38, 25: PRINT " мп-с"

233. CATE 37, 43: PRINT "Sfmax ="; CLNG(MAZG * 30000) / 100

234. CATE 38, 43: PRINT" градус"

235. CATE 37, 61: PRINT "Sm ="; M.SI.M * 300

236. CATE 38, 63: PRINT" (м-т)п.с"

237. SIGB = SIGBN: CH = 1 : 'ВЫВОД ГРАФИКОВ НА ЭКРАН

238. WHILE SIGB <= SIGBKON: SIG = SIGB * (2 * 3.141593 * 9.81 / LS) л .5: GOSUB 20000

239. NE (20 + (SIGB stepsig / 5 - SIGBN) * 225 / (SIGBKON - SIGBN), SZ / MSZ + 550)-(20 +

240. SIGB + stepsig / 5 SIGBN) * 225 / (SIGBKON - SIGBN), 550), 15, BF

241. NE (270 + (SIGB stepsig / 5 - SIGBN) * 225 / (SIGBKON - SIGBN), SZG(CH) / MSZG + 50)270 + (SIGB + stepsig / 5 SIGBN) * 225 / (SIGBKON - SIGBN), 50), 15, BF

242. NE (520 + (SIGB stepsig / 5 - SIGBN) * 225 / (SIGBKON - SIGBN), SAG(CH) / MAZG + 50)520 + (SIGB + stepsig / 5 SIGBN) * 225 / (SIGBKON - SIGBN), 50), 15, BF

243. NE (770 + (SIGB stepsig / 5 - SIGBN) * 225 / (SIGBKON - SIGBN), S.I.M(CH) / M.SI.M +50.-(770 + (SIGB + stepsig / 5 SIGBN) * 225 / (SIGBKON - SIGBN), 50), 15, BFfg(CH) = 270 + (SIGB SIGBN) * 225 / (SIGBKON - SIGBN)

244. FGG(CH) = 550 + FGG(CH) / MZET: FGGl(CH) = 550 + FGGl(CH) / MFIM1. CH = 1 THEN

245. PSET (fg(CH), FGG(CH)), COLLINE

246. PSET (fg(CH) + 250, FGGl(CH)), COLLINE

247. PSET (fg(CH) + 500, 550 + ACH.I.M(CH) / M.I.M), COLLINE1. ELSE

248. NE (fg(CH), FGG(CH))-(fg(CH 1), FGG(CH - 1)), COLLINE

249. NE (fg(CH) + 250, FGG1 (CH))-(fg(CH 1) + 250, FGG1(CH - 1)), COLLINE

250. NE (fg(CH) + 500, 550 + ACH.I.M(CH) / M.I.M)-(fg(CH 1) + 500, 550 + ACH.I.M(CH - 1) /1. M.I.M), COLLINE1. END IF

251. CATE 40, 5: PRINT "СТАНДАРТЫ:"

252. CATE 42, 1: PRINT "Вертикальной качки"

253. CATE 43, 1: PRINT SQR(dszg); "м"

254. CATE 45, 1: PRINT " Килевой качки"

255. CATE 46, 1: PRINT SQR(dsag); "град"

256. CATE 48, 1: PRINT "Изгибающего момента"

257. CATE 49, 1: PRINT SQR(d.i.m); Vm"вероятные амплитуды качкиfip = kd * dsag A .5 'перемещенияzp = kd * dszg Л .5fis = kd * dsas A .5 'скоростиzs = kd * dszs л .5fiu = kd * dsau A .5 'ускоРения zu = kd * dszu A .5

258. ERASE SZG, SAG, fg, FGG, FGG1, szgs, sags, szgu, sagu, SZ, sa COLOR 14: LOCATE 60, 1: PRINT "TAB ПРОДОЛЖИТЬ F9 - ПЕЧАТЬ"; wp$ = "": WHILE wp$ о CHR$(9): wp$ = INKEY$ IF wp$ = CHR$(0) + CHR$(67) THEN , IF (INP(&H379) AND 208) о 208 THEN PCOPY 0, 1

259. BEEP: CALL message("ncmr0T0BTe принтер!", "Продолжить-<Еп1ег> Отменить-<Езс>")wp$ = INKEY$:1. wp$ = CHR$(27) THEN

260. PCOPY 1,0: GOTO 123: END IF: END IF123 END IF

261. WEND: SCREEN 0: WIDTH 80, 50: CALL F0N1(4, 48, 1, 80): PCOPY 0, 1 OPEN "rollingY' + FIL$ FOR INPUT AS #1: CALL printdat(A$(), NAM, rr): CLOSE #1 PRINT #4,"": CLOSE #4: GOTO 90: END 20000 'ВЫБОР СПЕКТРА ВОЛНЕНИЯ IF variant.spektra = 1 THEN '12 МКОБ

262. SZ = (.7795 / SIG A 5) * EXP(-3.11 * 10 A 4 / (AH A 2 * SIG A 4)) 'волновой спектор END IF1. variant.spektra = 2 THEN '2 МКПС

263. SZ = 9.43 * .03575 * (AH3 /100)A 2 / (2 * 3.141593 / TV) * ((.777 * 2 * 3.141593 / TV) / SIG)A6 * EXP(-1.5 * ((.777 * 2 * 3.141593 / TV) / SIG)A 4)

264. SZ = 98.2 * (AH3 /100)A 2 * EXP(-686 / ((TV * SIG)A 4)) / (TV A 4 * SIG A 5)sz = .143 * (ah3 / 100) A 2 * ALS * BS2 * (SIG A 2 + ALS A 2) / (3.14 * (BES A 2 + 2 * ALS A 2) *

265. SIG A 4 + 2 * SIG A 2 * AS2 + BS2 A 2))

266. END IF: RETURN: 90 END SUB

267. SUB NONLINE (ZATA, FI)' НЕЛИНЕЙНЫЙ РАСЧЕТ КАЧКИ И МОМЕНТОВ

268. DIM Zl(-1 ТО Н% + 1), QQ(-1 ТО Н% + 1), QQI(-1 ТО Н% + 1), ММ(-1 ТО Н% + 1), PIR.S(-1

269. ТО Н% + 1), IZG.M(-1 ТО Н% + 1)

270. SCREEN 12: WIDTH 80, 60: WINDOW (0, -0)-(1000, 1000)

271. SHARED All, Bll, Cll, А12, В12, С12, A21S, B21S, C21S, A22S, B22S, C22S, FF, MM, QN,1. QNX, Fl, F2, Ml, M2

272. OPEN "non" FOR OUTPUT AS #100

273. PRINT #100, "ВРЕМЯ,с";" M(6)(T*M)";" M(10)(T*M)";" М(14)(Т*М)";" ВЕРТИКЛЬНАЯ М";" КИЛЕВАЯ РАД""РЕШЕНИЕ СИСТЕМЫ ДИФ. УРАВНЕНИЙ МЕТОДОМ РУНГЕ-КУТТА-ФЕЛЬБЕРГА"с АВТОМАТИЧЕСКИМ ВЫБОРОМ ШАГА"г**************** ЧИСЛО УРАВНЕНИЙ N

274. CLS : п = 4: DIM YY(n), L(n), W(n), A(n), BB(n), C(n), D(n), E(n), F(n)ka = 2 * 3.141593 / lv: SIG = (ka * 9.81) A .5: CALL runl5(lv, SIG, ka)

275. W(l) = ZETA1: W(2) = FI1: W(3) = ZETA2 * SIGk: W(4) = FI2 * SIGk

276. El = .03 ,*************** ПОГРЕШНОСТЬ ВЫЧИСЛЕНИЙ El

277. H = 2 * 3.141593 / (SIGk * 20) <*************** НАЧАЛЬНЫЙ ШАГ H

278. X = 0 .*************** НАЧАЛЬНОЕ Х0 X

279. FOR J = 1 ТО n: YY(J) = W(J): NEXT J <*************** НАЧАЛЬНОЕ YY0 W(J):1101 E3 = 0: GOSUB 400: DDD = 0: FOR J = 1 TO n

280. A(J) = F(J) * H: YY(J) = W(J) + 2 * A(J) / 9: NEXT J

281. X = X + 2 * H / 9: GOSUB 400: FOR J = 1 TO n

282. BB(J) = H * F(J): YY(J) = W(J) + A(J) /12 + BB(J) / 4: NEXT J X = X + H / 9: GOSUB 400: FOR J = 1 TO n:

283. C(J) = F(J) * H: YY(J) = W(J) + (34.5 * A(J) 121.5 * BB(J) + 135 * C(J)) / 64:

284. NEXT J: X = X + H / 2.4: GOSUB 400: FOR J = 1 TO n

285. D(J) F(J) * H: Q = W(J) - 17 * A(J) /12 + 6.75 * BB(J)

286. YY(J) = Q 5.4 * C(J) + 16 * D(J) /15: NEXT J:

287. X = X + H/4: GOSUB 400: FOR J = 1 TO n:

288. E(J) = F(J) * H: Q = W(J) + 65 * A(J) / 432 .3125 * BB(J)

289. YY(J) = Q + .8125 * C(J) + 4 * D(J) /27 + 5 * E(J) / 144: NEXT J

290. X = X H / 6: GOSUB 400: FOR J = 1 TO n

291. J) = H * F(J): YY(J) = W(J) + A(J) / 9 + .45 * C(J) + 16* D(J) / 45 + E(J) /12 Q = A(J) / 150 .03 * C(J) + 16 * D(J) / 75 + E(J) / 20 E2 = ABS(Q - .24 * L(J)): IF E2 <= El THEN 250 E3 = 1: GOTO 260

292. IF E2 < El / 20 THEN DDD = DDD + 1

293. NEXT J: X = X + H / 6: GOSUB 400: IF E3 = 0 THEN 290

294. X = X H: FOR J = 1 TO n: YY(J) = W(J): NEXT J1. H = H/2: GOTO 1101290 IF DDD = n THEN H = H + H1. X >= 40 THEN1.Ho.l THEN H = . 11. END IF

295. CLS 1: MMZZ 810: MMX = 1000 / LS: MMZ = 170/ X(KOL.TOCH(0) -1,0) FORI = 0 TO H%1. 170 / X(KOL.TOCH(I), I) < MMZ THEN MMZ = 170/ X(KOL.TOCH(I), I) NEXT I

296. DL = LS /H%: alfa# = (ALF / 180) * 3.141592654#: ka 2 * 3.141593 / lv: KA1 = -ka * COS(alfa#) MAX.Q = 0: MAX.S = 0: MAX.M = 0: FOR I = 0 TO H%

297. QQ(I) = -All(I) * F(3) PI 1(1) * YY(3) - CI 1(1) * YY(1) - A12(I) * F(4) - B12(I) * YY(4) - C12(I)

298. YY(2) + FS(I) * COS(SIGk * X) FSS(I) * SIN(SIGk * X) + Q.N(I)

299. QQI(I) = QQI(I 1) + QQ(I -1) + QQ(I): MM® - MM(I -1) + QQI(I -1) + QQI(I) NEXT I1. FOR1 = 0 TO H%

300. QQI(0) 0: MM(0) = 0: QQI(H%) = 0: MM(H%) = 0 PIR.S(I) = QQI(I) * DL / 2: IZG.M(I) = MM(I) * (DL / 2) A 2 IF MAX.Q < ABS(QQ(I)) THEN MAX.Q = ABS(QQ(I)) IF MAX.S < ABS(PIR.S(I)) THEN MAX.S = ABS(PIR.S(I))

301. MAX.M < ABS(IZG.M(I)) THEN MAX.M = ABS(IZG.M(I)): XA(T) = DL * ((H% / 2) -1) Z1(I) (DN - (DN - DK) * I / H% + YY(1) + XA(I) * YY(2) - ZA * COS(KAl * XA(I) + SIGk * X)) * MMZ: XA(I) = LS - DL * I

302. I = 0 THEN: PSET (XA(I) * MMX, Z1(I) + MMZZ): PSET (XA(I) * MMX, X(l, I) + MMZZ)

303. PSET (XA(I) * MMX, MMZ * X(KOL.TOCH(I) -1,1) + MMZZ): ELSE

304. NE (XA(I 1) * MMX, Z1(I - 1) + MMZZ)-(XA(I) * MMX, Z1(I) + MMZZ), 1

305. NE (XA(I 1) * MMX, MMZ * X(l, I - 1) + MMZZ)-(XA(I) * MMX, MMZ * X(l, I) + MMZZ)

306. NE (XA(I 1) * MMX, MMZ * X(KOL.TOCH(I - 1) - 1,1 - 1) + MMZZ)-(XA(I) * MMX, MMZ

307. X(KOL.TOCH(I) 1,1) + MMZZ), 4: END IF

308. NE (XA(I) * MMX, MMZZ 5)-(XA(I) * MMX, MMZ * X(KOL.TOCH(I) - 1,1) + MMZZ) NEXT I: LINE (0, DK * MMZ + MMZZ)-(LS * MMX, DN * MMZ + MMZZ) COLOR 15: LOCATE 13, 1: PRINT "20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 10"

309. FOR Z = 0 ТО 2 * Z.NACH STEP Z.NACH / 5 IF Z = Z.NACH THEN LINE (0, Z)-(MMX * LS, Z), 15 ELSE LINE (0, Z)-(MMX * LS, Z), 3 NEXT Z: FOR I = 0 TO H%

310. NE ((XA(I) DL * .5) * MMX, QQ(I) * MAX.Q + Z.NACH)-((XA(I) + DL * .5) * MMX, QQ(I) * MAX.Q + Z.NACH), 13

311. NE ((XA(I) + DL * .5) * MMX, QQ(I) * MAX.Q + Z.NACH)-((XA(I) + DL * .5) * MMX, Z.NACH), 13

312. NE ((XA(I) DL * .5) * MMX, QQ(I) * MAX.Q + Z.NACH)-((XA(I) - DL * .5) * MMX, Z.NACH), 13

313. NE (XA(I) * MMX, 0)-(XA(I) * MMX, 2 * Z.NACH), 3: NEXT I: FOR I = 0 TO H% 1 LINE (XA(I) * MMX, PIR.S(I) * MAX.S + Z.NACH)-(XA(I + 1) * MMX, PIR.S(I + 1) * MAX.S + Z.NACH), 4

314. NE (XA(I) * MMX, IZG.M(I) * MAX.M + Z.NACH)-(XA(I + 1) * MMX, IZG.M(I + 1) *

315. MAX.M + Z.NACH), 15: NEXT I

316. X >= 40 THEN: PRINT #100, USING "###.###"; X;

317. PRINT #100, USING "#########.##"; IZG.M(6); IZG.M(IO); IZG.M(14);

318. PRINT #100, USING "########.####"; YY(1); YY(2): END IF

319. COLOR 15: LOCATE 15, 15: PRINT " КИНЕМАТИЧЕСКИЕ ХАРАКТЕРИСТИКИ1. ПРОДОЛЬНОЙ КАЧКИ"

320. CATE 16, 5: PRINT "Время T="; USING "###.##"; X; : PRINT "с" LOCATE 16,25: PRINT " Вертикальная"

321. CATE 17, 25: PRINT "Перемещение "; USING "##.####"; YY(1); LOCATE 18, 25: PRINT "Скорость "; USING "##.####"; YY(3); LOCATE 19, 25: PRINT "Ускорение "; USING "##.####"; F(3); LOCATE 16,45: PRINT " Килевая"

322. F(3) = ((B21S * A12/A22S-B11) * YY(3) + (C21S * A12 / A22S CI 1) * YY(1) + (B22S * A12 / A22S - B12) * YY(4) + (C22S * A12 / A22S - C12) * YY(2) - A12 / A22S * MM + FF) / (All -A21S * A12/A22S)

323. F(4) = (-B22S * YY(4) C22S * YY(2) - A21S * F(3) - B21S * YY(3) - C21S * YY(1) + MM) RETURN: END END SUB

324. SUB polnover' ПОЛНОВЕРОЯТНОСТНЫЙ РАСЧЕТ КАЧКИ И МОМЕНТОВ

325. OPEN "polnover" FOR OUTPUT AS #44 sechenie =10: obes = 3

326. SIGBN = .5: SIGBKON = 2: stepsig = .1: AHN = .5: AHK = 12.5: STEP AH = 1: ALFN = 0: ALFK= 180: STEPALF= 15: TNAC = .4: TKON= 1.8: STEPT=.2: STEPLOADB = 1: LOADKB = 4: STEPFRUDB = 1: FRUDKB = 5: KFRUD = .05: NMOMENT = 10: VTV v

327. DM = .025 * (GAMMA * LS л 2 * Y(KOL.TOCH(H% \ 2) 2, H% \ 2) * 2): DZZ = 2.5: DFF = 2.5 'SM(ALFB, TB, HB)

328. DIM SW(30, 10, 13), ACH.M(20, 12), ACH.H(20, 12), ACH.P(20, 12)'АЧХ и сектор волнения

329. DIM AH.M(20, 12, 5), AH.H(20, 12, 5), AH.P(20, 12, 5)

330. DIM ACHM(20, 12, 5), АСНЩ20, 12, 5), ACHP(20,12, 5)

331. DIM SM(12, 10, 13), SH(12, 10, 13), SP(12, 10, 13) 'стандарты

332. DIM PALF(12), PAH(3, 13), PTE(8), PLOAD(5)

333. DIM PUM(IO), PTU(13), PRAION(3), QWM(IO), QHT(10, 4), SPEED(5)

334. DIM PUH(IO), PTUHH(13), QWH(10), QHTHH(10, 4)

335. DIM PUP(10), PTUPP(13), QWP(10), QHTPP(10, 4): NER= I

336. ВЕРОЯТНОСТИ РАЙОНОВ 'ВЕРОЯТНОСТИ НАГРУЗКИ 'ВЕРОЯТНОСТИ ПЕРИОДОВ В ДОЛЯХ ОТ СРЕДНЕГО

337. PRAION(l) PRAION(2) PRAION(3)3333 .3333 .3333

338. PLOAD(O) = .2: PTE(0) = .075 PLOAD(l) = .2: PTE(1) = .124 PTE(2) = .191 PTE(3) = .219 PTE(4) = .191

339. PLOAD(2) = .2 PLOAD(3) = .2: PLOAD(4) = .2:

340. ВСТРЕЧОЕ '15 '30 '45 '60 '75 '90 '105 '120 '135 '150 '165 '180

341. THEN: 'ПРИ УМЕРЕННОМ ВОЛНЕНИИ '0 ВСТРЕЧОЕ '15 '30 '45 '60 '75 '90 '105 '1201. КУРС.УГЛЫ1. PALF(3) =1/24: '1351. PALF(2) = 1 / 24: '1501. PALF(l) =1/24: '1651. PALF(O) =1/48: '180

342. END IF: 'ВЫБОР СЕЗОНА: tITLE$ = "СЕЗОН" NAM = 1: n = 4: XX = 28: YY = 9: DIM CS$(n)

343. CS$(1) = " ЗИМА ":CS$(2) = " ВЕСНА ": CS$(3) = " ЛЕТО

344. CS$(4) = " ОСЕНЬ ": A = v.menu(tITLE$, CS$(), n, XX, YY, NAM, sim%, FON%)

345. SEZON = NAM ZAGOLOVOK$ = tITLE$ + " -" + RTRIM$(CS$(NAM)) + " " + ZAGOLOVOK$

346. ERASE CS$: 'ВЫБОР МОРЯ: tITLE$ = "MOPE": NAM = 1: n = 6: XX = 46: YY = 12:1. DIM CS$(n)

347. CATE 27, 30: COLOR 1, 4: PRINT "--------------------"

348. ACHP(OM, ALFB, FRUDB) = ACHP(OM, ALFB, FRUDB) + FIM * ka * 57.29578 * PLOAD(load) 'АЧХ килевой качки NEXT LOADB OM = OM +1 NEXT SIGB

349. ALFB = ALFB + 1 NEXT ALF LOW = FRUDB

350. FOR FRUDB = 0 ТО FRUDKB STEP STEPFRUDB 'цикл по скорости STEPFRUDB=1 FRUDKB=5

351. FRUD FRUDB * KFRUD 'KFRUD=0.05: SPEED(FRUDB) = FRUD * (9.81 * LS) Л .5 IF SPEEDW <= SPEED(FRUDB) THEN: ALFB = 0 FOR ALF = ALFN TO ALFK STEP STEP ALFом = о

352. FOR SIGB = SIGBN TO SIGBKON STEP stepsig ACH.M(OM, ALFB) = ACHM(OM, ALFB, FRUDB -1) + (ACHM(OM, ALFB, FRUDB) -ACHM(OM, ALFB, FRUDB -1) * (SPEEDW SPEED(FRUDB -1)) / (SPEED(FRUDB) -SPEED(FRUDB - 1)))

353. ACH.H(OM, ALFB) = ACHH(OM, ALFB, FRUDB 1) + (ACHH(OM, ALFB, FRUDB) -ACHH(OM, ALFB, FRUDB - 1) * (SPEEDW - SPEED(FRUDB - 1)) / (SPEED(FRUDB) -SPEED(FRUDB -1)))

354. FOR ALF = ALFN TO ALFK STEP STEP ALF alfa# = (ALF / 180) * 3.141592654# SM(ALFB, ТВ, HB) = 0: SH(ALFB, ТВ, HB) 0: SP(ALFB, ТВ, HB) - 0ом = о

355. FOR SIGB = SIGBN TO SIGBKON STEP stepsig '===============dsig = stepsig * (2 * 3.141593 * 9.81 / LS) A .5 '|dsigk = dsig dsig л 2 * v * COS(alfa#) / 9.8066 "j

356. SM(ALFB, TB, HB) = SQR(SM(ALFB, TB, HB) * dsigk) SH(ALFB, TB, HB) SQR(SH(ALFB, TB, HB) * dsigk) SP(ALFB, TB, HB) = SQR(SP(ALFB, TB, HB) * dsigk) ALFB = ALFB + 1 NEXT ALF NEXT TB HB = HB + 1 LO = LOW + HB \ 2

357. CATE 27, 30 + LO: COLOR 14, 4: PRINT "-": NEXT АНЗ i******************* РАСЧЕТ ПОЛНОЙ ВЕРОЯТНОСТИ ******************** FOR MOMENT = 0 ТО NMOMENT 'ЦИКЛ ПО МОМЕНТАМ

358. МО = DM * MOMENT: ZZO = DZZ * MOMENT: FFO = DFF * MOMENT QWM(MOMENT) = 0: QWH(MOMENT) = 0: QWP(MOMENT) = 0 FOR RAION = 1 TO NRAION 'ЦИКЛ ПО РАЙОНАМ ПЛАВАНИЯ

359. QHT(MOMENT,RAION) = 0:QHTHH(MOMENT,RAION) = 0: QHTPP(MOMENT, RAION) = 0 HB = 1 'НАЧАЛО ЦИКЛА ПО АНЗ%

360. FOR АНЗ = AHN ТО АНК STEP STEP АН TS = 3.3 * SQR(AH3) PTU(HB) = 0: PTUHH(HB) = 0: PTUPP(HB) = 0 FOR TB = 0 TO 7 STEP 1 'НАЧАЛО ЦИКЛА ПО TV

361. TV (TB * STEPT + TNAC) * TS PUM(TB) = 0: PUH(TB) = 0: PUP(TB) = 0

362. ALFB = 0 'НАЧАЛО ЦИКЛА ПО ALF

363. PRINT #44, "PTU(HB)", AH3, PTU(HB)

364. QHT(MOMENT, RAION) = QHT(MOMENT, RAION) + PTU(HB) * PAH(RAION, HB) QHTHH(MOMENT,RAION)=QHTHH(MOMENT,RAION)+PTUHH(HB) * PAH(RAlON, HB) QHTPP(MOMENT, RAION) = QHTPP(MOMENT, RAION) + PTUPP(HB) * PAH(RAION, HB) HB = HB + 1: NEXT АНЗ

365. SCREEN 12: WIDTH 80, 60: WINDOW (0, -0)-(Ю00,1000): ZMM = 800 / (DM * NMOMENT) ZMHP=800/(DZZ*NMOMENT): XMM = 800 / 10'XMM: A = setka(100, 900, 10, 100, 900, 10) FOR MOMENT = 0 TO NMOMENT -1 'ЦИКЛ ПО МОМЕНТАМ

366. МО = DM * MOMENT: ZZ0 = DZZ * MOMENT: FF0 = DFF * MOMENT Ml = DM * (MOMENT + 1): ZZ1 = DZZ * (MOMENT + 1): FF1 = DFF * (MOMENT + 1) IF -LOG(QWM(MOMENT + 1)) / LOG(IO) > 10.5 THEN 14

367. NE (100 + -LOG(QWM(MOMENT)) / LOG(IO) * XMM, M0 * ZMM + 100)-(100 +

368. G(QWM(MOMENT + 1)) / LOG(IO) * XMM, Ml * ZMM + 100), 14

369. NE (100 + -LOG(QWH(MOMENT)) / LOG(IO) * XMM, ZZ0 * ZMHP + 100)-(100 +

370. G(QWH(MOMENT + 1)) / LOG(IO) * XMM, ZZ1 * ZMHP + 100), 15

371. NE (100 + -LOG(QWP(MOMENT)) / LOG(IO) * XMM, FF0 * ZMHP + 100)-(100 +

372. G(QWP(MOMENT + 1)) / LOG(IO) * XMM, FF1 * ZMHP + 100), 131. NEXT MOMENT14 ТРАФИК

373. COLOR 14: LOCATE 1,1: PRINT " ОБЕСПЕЧЕННОСТЬ ВОЛНОВЫХ

374. ВЕРТИКАЛЬНЫХ ИЗГИБАЮЩИХ МОМЕНТОВ"

375. CATE 3,40 LEN(ZAGOLOVOK$) \ 2: PRINT ZAGOLOVOK$

376. CATE 56, 9: PRINT "0 1 2 3 4 5 6 7 8 9 10 -Log Q";

377. CATE 2, 2: PRINT " Mb": LOCATE 3, 2: PRINT "—": LOCATE 4, 2: PRINT "pgBLa"

378. CATE 6, 4: PRINT "0.250": LOCATE 11,4: PRINT "0.225":

379. CATE 16, 4: PRINT "0.200": LOCATE 21, 4: PRINT "0.175"

380. CATE 26, 4: PRINT "0.150": LOCATE 30, 4: PRINT "0.125"

381. CATE 35, 4: PRINT "0.100": LOCATE 40, 4: PRINT "0.075"

382. CATE 45, 4: PRINT "0.050": LOCATE 50, 4: PRINT "0.025"

383. CATE 60, 10: PRINT "ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ <ESC>";: DO: LOOP WHILE INKEY$ о CHR$(27)

384. SCREEN 0: WIDTH 80, 50: OPEN "rollingV + FIL$ FOR INPUT AS #1: CALL printdat(A$(), NAM, rr): CLOSE #1

385. OPEN "polnover" FOR OUTPUT AS #44: HB = 1: PRINT #44, "СТАНДАРТЫ МОМЕНТОВ" PRINT #44," ALF, TV, AH3, SM(ALFB, TB, HB)": FOR AH3 = AHN TO AHK STEP STEP AH TS = 3.3 * SQR(AH3): FOR ТВ = 0 TO 7 STEP 1 'НАЧАЛО ЦИКЛА ПО TV

386. TV = (TB * STEPT + TNAC) * TS

387. ALFB = 0 'НАЧАЛО ЦИКЛА ПО ALF

388. FOR ALF = ALFN TO ALFK STEP STEP ALF PRINT #44, ALF, TV, AH3, SM(ALFB, TB, HB) ALFB ALFB + 1 NEXT ALF NEXT TB HB = HB + 1 NEXT AH3

389. PRINT #44, "PUM(TB)": 'FOR ТВ = 0 TO 7 STEP 1 'НАЧАЛО ЦИКЛА ПО TV

390. PRINT #44, ТВ, PUM(TB): 'NEXT ТВ

391. PRINT #44, "PTU(HB)": 'HB = 1: 'FOR АНЗ = AHN TO AHK STEP STEP AH 'PRINT #44, АНЗ, PTU(HB): 'HB = HB + 1: 'NEXT АНЗ PRINT #44, "MOMENT, RAION, QHT(MOMENT, RAION)" FOR MOMENT = 0 TO NMOMENT 'ЦИКЛ ПО МОМЕНТАМ

392. FOR RAION = 1 TO NRAION 'ЦИКЛ ПО РАЙОНАМ ПЛАВАНИЯ

393. PRINT #44, MOMENT, RAION, QHT(MOMENT, RAION): NEXT RAION NEXT MOMENT: END SUB

394. SUB REGULYR "РАСЧЕТ НА РЕГУЛЯРНОМ ВОЛНЕНИИ В ЛИНЕЙНОЙ ПОСТАНОВКЕ

395. DIM Z1(H% + 1), ZMP(-1 ТО Н% + 1), ZMA(-1 ТО Н% + 1)

396. Z1(I) = (DN (DN - DK) * I / H% + ZETAG(IE) + XA(I) * FIG(IE) / 57.29578 - ZA * COS(KAl * XA(I) + SIGk * IE)) * MMZ: XA(I) = LS - DL * I IF I = H% OR I = 0 THEN

397. NE (XA(I) * MMX, MMZ * X(l, I) + MMZZ)-(XA(I) * MMX, MMZ * X(KOL.TOCH(I) 1,1) + MMZZ), 4: END IF: IF 1 = 0 THEN

398. NE (XA(I) * MMX, MMZ * X(l, I) + MMZZ)-(XA(I) * MMX, MMZ * X(KOL.TOCH(I) -1,1) + MMZZ), 4: ELSE

399. NE (0, DK * MMZ + MMZZ)-(LS * MMX, DN * MMZ + MMZZ), 3

400. CATE IE * 6 + 6, 59: PRINT USING "#"; IE: LOCATE IE * 6 + 6, 60: PRINT "•П/4"

401. CATE IE * 6 + 6, 66: PRINT USING "###.###"; ZETAG(IE)

402. CATE IE * 6 + 6, 74: PRINT USING "###.###"; FIG(IE)

403. NE (710, MMZZ)-(1000, MMZZ): NEXT IE

404. CATE 3, 58: PRINT "Положен.Вертик. Измен."

405. CATE 4, 58: PRINT " волны перем,м.угла,гр"

406. NE (XA(I) * MMX, 5 + MMZ * X(l, I))-(XA(I) * MMX, 5 + MMZ * X(KOL.TOCH(I) 1,1)), 4

407. ELSE: LINE (XA(I 1) * MMX, 5 + MMZ * X(l, I - 1))-(XA(I) * MMX, 5 + MMZ * X(l, I)), 4

408. NE (XA(I 1) * MMX, 5 + MMZ * X(KO.L.TOCH(I -1) -1,1 - 1))-(XA(I) * MMX, 5 + MMZ *

409. X(KOL.TOCH(I) 1,1)), 4: END IF: LINE (XA(I) * MMX, 10)-(XA(I) * MMX, 5)

410. NE (0, 5 + DK * MMZ)-(LS * MMX, 5 + DN * MMZ), 3: IF I = 0 THEN 3

411. NE (XA(I) * MMX, 5 + ZMP(I) * MMZ)-(XA(I -1) * MMX, 5 + ZMP(I -1) * MMZ)

412. NE (XA(I 1) * MMX, MMZ * X(l, I - 1) + MMZZ)-(XA(I) * MMX, MMZ * X(l, I) + MMZZ) LINE (XA(I - 1) * MMX, MMZ * X(KOL.TOCH(I -1) - 1,1 -1) + MMZZ)-(XA(I) * MMX, MMZ

413. X(KOL.TOCH(I) 1,1) + MMZZ), 4: END IF

414. NE (XA(I) * MMX, MMZZ 5)-(XA(I) * MMX, MMZ * X(KOL.TOCH(I) - 1,1) + MMZZ) NEXT I: LINE (0, DK * MMZ + MMZZ)-(LS * MMX, DN * MMZ + MMZZ) COLOR 15: LOCATE 13,1: PRINT "20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 0"

415. ОТНОСИТЕЛЬНЫЙ МОМЕНТ Mw. = Mw/(LA2 * В * Rw* gamma)

416. COLOR 15: LOCATE 19, 15: PRINT " ХАРАКТЕРИСТИКИ НАГРУЗКИ"

417. COLOR 13: LOCATE 20, 5: PRINT "Максим, погонная нагрузка на шпацию t/m Р ="; MAX.Q

418. COLOR 11: LOCATE 21,5: PRINT "Амплитуда макс, перерезывающей силы, t R ="; MAX.S;1. P.S(10)

419. MAX = 0: IF ABS(MAX) < ABS(MAX.M) THEN MAX = MAX.M IF ABS(MAX) < ABS(MAX.MTV) THEN MAX = MAX.MTV IF ABS(MAX) < ABS(M.SUM) THEN MAX = M.SUM MAX = Z.NACH / MAX 'МАСШТАБ ДЛЯ MOMEHTOBB FOR Z = 0 TO 2 * Z.NACH STEP Z.NACH / 5

420. Z = Z.NACH THEN LINE (0, Z)-(MMX * LS, Z), 15 ELSE LINE (0, Z)-(MMX * LS, Z), 3 NEXT Z: FOR I = 0 TO H%: PP = PI(I) / DL

421. NE ((XA(I) DL * .5) * MMX, PP * MAX.Q + Z.NACH)-((XA(I) + DL * .5) * MMX, PP * MAX.Q + Z.NACH), 13

422. NE ((XA(I) + DL * .5) * MMX, PP * MAX.Q + Z.NACH)-((XA(I) + DL * .5) * MMX, Z.NACH), 13

423. NE ((XA(I) DL * . 5) * MMX,PP * MAX. Q+Z.N ACH)-((XA(I) - DL * .5) * MMX, Z.NACH), 13 NEXT I: FOR I = 0 TO H% -1

424. NE (XA(I) * MMX, P.S(I) * ABS(MAX.S) + Z.NACH)-(XA(I + 1) * MMX, P.S(I + 1) * ABS(MAX.S) + Z.NACH), 11

425. NE (XA(I) * MMX, I.M(I) * ABS(MAX) + Z.NACH)-(XA(I + 1) * MMX, I.M(I + 1) * ABS(MAX) + Z.NACH), 15

426. SCREEN 0: WIDTH 80, 50: A% = RESTSCR%("OPENF.MJG"): GGG = 0 END SUB

427. SUB RUN.NONLINE (ZETA, FI, ZETAD, FID, IE)' НЕЛИНЕЙНЫЙ РАСЧЕТ КОФ. УР-НИЙ КАЧКИ

428. DIM QQ(H% + 1), QQI(H°/o + 1), MM(H% + 1), HP(H% + 1), BP(H% + 1), NACH.B(H% + 1)

429. DIM KSI11#(H% + 1), KSI21#(H% + 1), KSI51#(H% + 1), KSI12#(H% + 1), KSI22#(H% + 1),

430. DIM sil 1#(H% + 1, 12), si21#(H% + 1, 12), si51#(H% + 1, 12), sil2#(H% + 1, 12), si22#(H% + 1)

431. DIM KS11#(H% + 1), KS21#(H% + 1), KS51#(H% + 1), TT(H% + 1), BB(H% + 1), GG(H% + 1)

432. DIM DZ(H% + 1), DZM3(H% + 1), DZM30(H% + 1), DZL3(H% + 1), DZL30(H% + 1)

433. DIM M33K(H% + 1), M330K(H% + 1), L33K(H% + 1), L330K(H% + 1)

434. DIM M33N(H% + 1), M330N(H% + 1), L33N(H% + 1), L330N(H% + 1)

435. РАСЧЕТ ВСПОМОГАТЕЛЬНЫХ ВЕЛИЧИНka = 2 * 3.141593 / lv: SIG = (ka * 9.81) A .5alfa# = (ALF/ 180) * 3.141592654#

436. DL = LS / H%: KA1 = -ka * COS(alfa#): KA2 = ka * SIN(alfa#): I = 0: WHILE I <= H% XA(I) = DL * ((H% / 2) -1)

437. NAP (I) = TT(I) + X( 1,1) HP(I) TT(I) - HP(I) - X(l, I)

438. BP(I) = Y(KOL.TOCH(I) 2,1) * 2: HP(I) = X(KOL.TOCH(I) - 2,1) IF TT(I) + X(l, I) >= HP(I) THEN BB(I) = BP(I)

439. NAP(I) = TT(I) + X(l, I) HP(I) TT(I) = HP(I) - X(l, I)

440. NAP (I) >= .4 * BP(I) THEN NAP(I) = .4 * BP(I) END IF

441. NAP(I) = TT(I) + X(l, I) HP(I) TT(I) = HP(I)-X(1,I)

442. NAP (I) >= .4 * BP(I) THEN NAP® = .4 * BP(I) END IF

443. GG(I) = 2 * INTEGRAL. CEB!(X(), Y(), I, KOL.TOCH(), GG(I), TT(I) + X(l, I))1. = INETG.TRAP(LS, GG(I), H%, I, IN)

444. CALL SPLAIN(X(), Y(), TT(I), I, TT(I) + X(l, I), Y, KOL.TOCH()):

445. BB(I) = Y * 2: NACH.B(I) = BB(I)

446. BB(I) <= 0 THEN BB(I) = .0001

447. BATAA(I) = GG(I) / (TT(I) * BB(I))1. PACHETKnMHK£

448. CALL KPMKD(ka, TT(), BB(), BATAA(), NAP(I), BP(I), I, M33(), M330(), L33(), L330(), SIG) 55:1 = 1 + 1: WEND: D = IN * GAMMA

449. DM33(0) = (M33(0) M33(l)) * 2: DM330(0) = (M330(0) - M330(l)) * 2

450. DM33(H%) = (M33(H% 1) - M33(H%)) * 2: DM330(H%) = (M330(H% - 1) - M330(H%)) * 2

451. DL33(0) = (L33(0) L33(l)) * 2: DL330(0) = (L330(0) - L330(l)) * 2

452. DL33(H%) = (L33(H% 1) - L33(H%)) * 2: DL330(H%) - (L330(H% - 1) - L330(H%)) * 21. FOR I = 1 TO H% -1

453. DM33(1) = M33(I -1) M33(I + 1): DM330(I) = M330(I -1) - M330(I + 1) DL33(I) = L33(I -1) - L33(I + 1): DL330(I) = L330(I -1) - L330(I + 1) NEXT I

454. B12(I) = (L33(I) * XA(I) 2 * v * M33(I) - v / 2 * DM33(1) * XA(I) / (2 * DL)) 12 CI2(1) = CI 1(1) * XA(I) - (2 * v * L33(I) + v A 2 / 2 * DM33(I) / (2 * DL)) / 2: ELSE A11(1) = PI(I) / (9.81 * DL) + M33(I) B11(1) - L33(I) - v / 2 * DM33(1) / (2 * DL)

455. Z1 = (jl / qla) * TT(I) * (1 -1 / qla) wl = EXP(-ka * Zl)wwwl# = KA2 * .5 * BB(I) * (1 (Zl / TT(I)) л M) sill#(I,jl) = wl * SIN(wwwl#) si21#(I, jl) = wl * COS(wwwl#) * Zl л (M -1) si51#(I, jl) = Zl * si21#(I,jl)

456. KSI11 #(I) = KSI11 #(I) + si 11 #(I, j 1) KSI21#(I) = KSI21#(I) + si21#(I, jl) KSI51 #(I) = KSI51 #(I) + si51 #(I, j 1) jl =jl + 1: WEND

457. КАРА = (M / (TT(I) л (M 1) * qla * Q2A)) * ((Q2A - Q2A / qla) * KS21#(I) + (KSI22#(I) -(si22#(I, 0) + si22#(I, Q2A)) / 2))ksii5 = ((M + 1) / ((TT(I) л M) * qla * Q2A)) * ((Q2A Q2A / qla) * KS51#(I) + (KSI52#(I) si52#(I, 0) + si52#(I, Q2A)) / 2))

458. KAPA1 = ksiil * SIN(alfa#) + ksii5 * COS(alfa#) л 2

459. КАРА > 1 THEN КАРА = 1: IF KAPA1 > 1 THEN KAPA1 = 166: El = КАРА * GAMMA * BB(I): E2 = L330(I)

460. НЗ = KAPAl * SIG л 2 * M330(I): E4 = KAPA1 * v / 2 * DL330(I) / (DL * 2) E5 = v / 2 * DM330(I) / (DL * 2)

461. F11 = ZA * (El E3 - E4): F22 = ZA * KAPA1 * SIG * (E2 - E5)

462. Ml 1 = Fll * X: M22 = F22 * X

463. CKX = C0S(KA1 * X): SKX = SIN(KA1 * X)

464. ZV = ZA * COS(KAl * X + SIGk * IE): ZVD -ZA * SIGk * SIN(KA1 * X + SIGk * IE) ZXFVFZ = ZETAD + X * FID - v * FI - ZVD: ZXFVF = ZXFVFZ + ZVD ZXFZ = ZETA + X * FI - ZV: ZXF = ZXFZ + ZV

465. Q.N(I) = -.5 * (ZXFVFZ * (DZM3(I) * ZXFVF DZM30(I) * ZVD) + ZXFZ * (DZL3(I) * ZXF -DZL30(I) * ZVD))

466. ZXFVFZ < 0 THEN Q.N(I) = -.5 * (ZXFZ * (DZL3(I) * ZXF DZL30(I) * ZVD))

467. I = 0 OR I = H% THEN Q.N(I) = Q.N(I) / 2: QN = QN + Q.N(I): QNX = QNX + Q.N(I) * X

468. NEXT I: FF = F1* СOS(SIGk * IE) F2 * SIN(SIGk * IE) + QN

469. MM = Ml * COS(SIGk * IE) M2 * SIN(SIGk * IE) + QNX

470. ERASE МЗЗ, M330, L33, L330, DM33, DM330, DL33, DL330

471. ERASE KSI11#, KSI21#, KSI51#, KSI12#, KSI22#, KSI52#

472. ERASE sil 1#, si21#, si51#, sil2#, si22#, si52#, KS11#, KS21#, KS51#

473. ERASE TT, BB, GG, BATAA, NAP1. END SUB

474. SUB run 15 (Iv, SIG, ка)" РАСЧЕТ КОЭФФИЦИЕНТОВ УРАВНЕНИЙ КАЧКИ

475. SHARED A11, B11, CI 1, A12, B12, C12, A21S, B21S, C21S, A22S, B22S, C22S, FF, MM, QN,1. QNX, F1,F2,M1,M21. ПО ВЫВЕДЕННЫМ ФОРМУЛАМ

476. A11 = 0: B11 = 0: CI 1 = 0: A12 = 0: B12 = 0: C12 0: A21S = 0: B21S = 0: C21S = 0: A22S = 0: B22S = 0: C22S = 0: В = 0: Fl - 0: F2 = 0: Ml = 0: M2 = 0

477. DIM M33(H% + 1), M330(H% + 1), L33(H% + 1), L330(H% + 1), DM33(H% + 1), DM330(H% +1., DL33(H% + 1), DL330(H% + 1)

478. DIM QQ(H% + 1), QQI(H% + 1), MM(H% + 1), GS(H% + 1)

479. DIM sil 1#(12), si21#(12), si51#(12), sil2#(12), si22#(12), si52#(12)

480. RO = GAMMA / 9.81: alfa# = (ALF / 180) * 3.141592654#: DL = LS / H%

481. SIGk = SIG ((SIG A 2 / 9.81) * v * COS(alfa#)): KA1 = -ka * COS(alfa#)

482. KA2 = ka * SIN(alfa#): ZA = .11 * lv A .715 'полувысота волны

483. ZA = 4.67:1 = 0: WHILE I <= H%

484. CALL KPMKD(ka, T(), B(), BATA(), 0, 0,1, M33(), M330(), L33(), L330(), SIG)

485. DMN3(I) = M33(I): DMN30(I) = M330(I): DLN3(I) = L33(I): DLN30(I) = L330(I): TTN(I) = T(I)1 = 1+1: WEND: FOR I = 1 TO H% 1

486. DM33(I) = M33(I 1) - M33(I + 1): DM330(I) = M330(I -1) - M330(I + 1)

487. DL33(I) = L33(I 1) - L33(I + 1): DL330(I) = L330(I -1) - L330(I + 1): NEXT I

488. DM33(0) = (M33(0) M33(l)) * 2: DM330(0) = (M330(0) - M330(l)) * 2

489. DM33(H%) = (M33(H% -1) M33(H%)) * 2: DM330(H%) = (M330(H% - 1) - M330(H%)) * 2

490. DL33(0) = (L33(0) L33(l)) * 2: DL330(0) = (L330(0) - L330(l)) * 2

491. DL33(H%) = (L33(H% -1) L33(H%)) * 2: DL330(H%) = (L330(H% -1) - L330(H%)) * 2

492. CALL raspred.nagruzki: FOR I = 0 TO H%: XA(I) = DL * ((H% / 2) -1): IF I = 0 OR I = H%

493. THEN: A11(1) = (2 * PI(I) / (9.81 * DL) + M33(I)) / 2: B11(1) = (L33(I) v / 2 * DM33(I) / (2 *

494. DL)) / 2: CI 1(1) = (GAMMA * B(I) v / 2 * DL33(I) / (2 * DL)) / 2: A12(I) = All(I) * XA(I)

495. КАРА = (М / (Т(1) л (М 1) * qla * Q2A)) * ((Q2A - Q2A / qla) * KS21# + (KSI22# - (si22#(0) + si22#(Q2A)) / 2))ksii5 = ((M + 1) / ((T(I) л M) * qla * Q2A)) * ((Q2A Q2A / qla) * KS51# + (KSI52# - (si52#(0) + si52#(Q2A)) / 2))

496. KAPA1 = ksiil * SIN(alfa#) + ksii5 * COS(alfa#) A 2

497. G(I) = GS(I): IF КАРА > 1 THEN КАРА = 1: IF KAPA1 > 1 THEN KAPA1 = 1 , El = КАРА * GAMMA * B(I): E3 = KAPA1 * SIG A 2 * M330(I) E4 = KAPA1 * v / 2 * DL330(I) / (DL * 2): E2 = L330(I) E5 = v / 2 * DM3 30(1) / (DL * 2)

498. F11 = ZA * (El E3 - E4): F22 = ZA * KAPA1 * SIG * (E2 - E5)

499. Mil — F11 * X: M22 = F22 * X

500. CKX = C0S(KA1 * X): SKX = SIN(KA1 * X)

501. FIM = ((FI1 A 2 + FI2 A 2) A .5) / alfO ERASE A, letb, C, letg, letx

502. J = 1: WHILE J <= 2: 'определение производных FI, ZETA

503. РЕШЕНИЕ СИСТЕМЫ FID, ZETAD

504. FID = FI1 * COS(OMEGA#) + FI2 * SIN(OMEGA#)

505. ZETAD = ZETA1 * COS(OMEGA#) + ZETA2 * SIN(OMEGA#)

506. FIM = ((FI1 A 2 + FI2 л 2) л .5) / alffl1. FOR I = 0 TO H%1. I = 0 OR I = H% THEN

507. QQ(I) = (-A11(1) * ZETAD2 В11 (I) * ZETAD1 - CI 1(1) * ZETAD - A12(I) * FID2 - B12(I) * FID1 - С 12(1) * FID + FS(I) * COS(OMEGA#) - FSS(I) * SIN(OMEGA#)) * 2 QQI(I) = 0: MM(I) = 0: ELSE

508. QQ(I) = -A11(1) * ZETAD2 В11 (I) * ZETAD 1 - CI 1(1) * ZETAD - A12(I) * FID2 - B12(I) *

509. FID1 С 12(1) * FID + FS(I) * COS(OMEGA#) - FSS(I) * SIN(OMEGA#)

510. QQI(I) = QQI(I 1) + QQ(I -1) + QQ(I): MM® = MM(I -1) + QQI(I -1) + QQI(I): END IF

511. PIR.SILA(I, J) = QQI(I) * DL/2: IZG.MOMENT(I, J) = MM(I) * (DL / 2) л 2: NEXT I: J = J + 1

512. OMEGA# = 3,141592654# / 2: WEND: MAX.S = 0: MAX.M = 0: FOR I = 0 TO H%

513. P.S(I) = SQR(PIR.SILA(I, 1) л 2 + PIR.SILA(I, 2) л 2)

514. M(I) = SQR(IZG.MOMENT(I, 1) л 2 + IZG.MOMENT(I, 2) л 2)

515. MAX.S < ABS(P.S(I)) THEN MAX.S = ABS(P.S(I))

516. MAX.M < ABS(I.M(I)) THEN MAX.M = ABS(I.M(I)): NEXT I

517. ERASE МЗЗ, M330, L33, L330, DM33, DM330, DL33, DL330

518. ERASE sill#, si21#, si51#, sil2#, si22#, si52#1. ERASE QQI, MM1. END SUB

519. ПОДПРОГРАММА ИСХОДНЫХ ДАННЫХ

520. OPEN "rollingV + NEMA$ FOR OUTPUT AS #1

521. PRINT #1, TAB(28);" ИСХОДНЫЕ ДАННЫЕ ДЛЯ "

522. PRINT #1, TAB(24); "РАСЧЕТА ПРОДОЛЬНОЙ КАЧКИ СУДНА"1. PRINT #1,""

523. PRINT #1, "Количество расчетных шпаций Н = "; Н%

524. PRINT #1, "Водоизмещение судна, т. D = "; D

525. PRINT #1, "Длина судна по ГВЛ. м. L = "; LS

526. PRINT #1, "Осадка судна на носовом перпен-ре м Тн= "; DN

527. PRINT #1, "Осадка судна на кормовом перпен-ре м Тк= "; DK

528. PRINT #1, "Скорость хода судна, м/с. V = "; v

529. PRINT #1, "Угол встречи с волной, градусы ALFA = "; ALF

530. PRINT #1, "Абсцисса центра тяжести судна, м. Xg= "; Xg

531. PRINT #1, "Момент инерции относит. ЦТ. м.А4 Iyy = "; iyy

532. PRINT #1, "Момент инерции площади ГВЛ. м.л4 ffi = "; Ifi

533. PRINT #1, "Длина волны, м. LV ="; lv

534. PRINT # 1, "Удельный вес воды, т/куб.м GAMMA = "; GAMMA1. PRINT #1,""

535. SUB star (NAM, R) " ОПРЕДЕЛЯЮЩАЯ ПОДПРОГРАММА

536. ON KEY(l) GOSUB 1010: 'KEY(l) ON: 'ON ERROR GOTO 77: 'ON KEY(2) GOSUB 1020 'KEY(2) ON: 'ON KEY(3) GOSUB 1030: 'KEY(3) ON: 'ON KEY(4) GOSUB 1040: KEY(4) ON 'ON KEY(5) GOSUB 1050: 'KEY(5) ON: 'ON KEY(6) GOSUB 1060: 'KEY(6) ON

537. ON KEY(9) GOSUB 1090: 'KEY(9) ON: ON KEY(IO) GOSUB 1000: KEY(10) ON

538. NAM = 1: R = 1: Ql = 1: WHILE Q1 = 1: CALL submenu(NAM, R)файл

539. SELECT CASE NAM CASE 1: SELECT CASER CASE 1: CALL NEWFILE CASE 2: CALL openfile CASE 3:1. CALL SAVE CASE 4:

540. CALL messagel(n$, "СОХРАНИТЬ КАК."): CALL SAVE END SELECT 'редактирование CASE 2: SELECT CASE R

541. CASE 1: CALL REDTOTAL CASE 2: CALL RED(X(), Y(), KOL.TOCH()) END SELECT 'расчет

542. CASE 1: CALL prnscreen(47) END SELECT CASE 7: SELECT CASER

543. CASE 2: CALL qlstn(22, 30, "ВЫХОД ?", F) IF F = 1 THEN

544. CALL qlstn(22, 30, "СОХРАНИТЬ ДАННЫЕ ?", Ql) IFQ1 = 1 THEN CALL SAVE END IF

545. Ql = 0: END IF: END SELECT: END SELECT: WEND: KILL "*'.MJG": END END SUB

546. Calculation of Ship Motion and Bending Moment and Comparison ofthe Programs

547. Oleg E.Surov, Far -Eastern State Technical University. II-Woong Kim, University of Ulsan1998 Ulsan

548. GENERAL INFORMATION GENERAL PARTICULARS

549. THIS BOOKLET IS PREPARED TO ENABLE THE TRIM & STABILITY AND LONGITUDINAL STRENGTH CALCULATION TO BE EASILY PERFORMED.

550. THE LONGITUDINAL STRENGTH CALCULATION IS MADE TO CHECK THAT THE BENDING MOMENTS AND HULL-GIRDER SHEAR FORCES ARE WITHIN THE PERMISSIBLE LIMITS.1. PARTICULARS OF THE VESSEL1. NAME OF THE VESSEL

551. LENGTH BETWEEN PERPENDICULARS3. BREADTH (MOULDED)4. J3EPTH (MOULDED)5. DESIGN DRAUGHT (MOULDED)

552. M/T ^UKONG NAVIGATOR" 315.00 M 57.20 M 30.40 M 20.45 M 20.45 M

553. B&W LICENSEE 7S80MC 31,670 BHP x 72.2 RPM 28,500 BHP x 69.7 RPM1. YUK0N6- LINE LIMITED1. PANAMA1. PANAMA3 F J U 5154651 I 83050166919.1 / 156995.9

554. KOREAN REGISTER OF SHIPPING +KRS1, +KRM1, UMA-2

555. AMERICAN BEREAU OF SHIPPING

556. A1(E), "OIL CARRIER", +AMS, +ACCUii. tank plan1. PBOft-E