Влияние проектных характеристик судна на его мореходные и прочностные качества
Диссертация
При оценке мореходности в первую очередь необходимо знание параметров качки. Основоположником теории продольной качки может по праву считаться А. Н. Крылов, который еще сто лет назад изложил свои исследования в этой области. В дальнейшем на основе теории А. Н. Крылова И.Г. Бубнов разработал практический метод динамического расчета прочности судов при продольной качке на волнении. Последующее… Читать ещё >
Список литературы
- Антоненко C.B. К вопросу о влиянии продольного радиуса инерции масс на волновые изгибающие моменты. НТО Судпрома, вып. 131, Л.: Судостроение, 1969.-С. 81 -90.
- Антоненко C.B., Суров О. Э. Влияние распределения нагрузки на величину волнового изгибающего момента. // Тр. междунар. конф. Кораблестроение и океанотехника. Проблемы и перспективы: ДВГТУ, Владивосток, 1998. С. 217−221.
- Антоненко C.B., Суров О. Э. Исследование влияния несимметричного распределения нагрузки на волновые изгибающие моменты. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1999. С. 141 — 145.
- Антоненко C.B., Суров О. Э. Использование модели ветровых волн для расчета поведения судна в море. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1999.-С. 135 140.
- Антоненко C.B., Суров О. Э. Расчет мореходности и прочности судна в заданном режиме волнения. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1996. С. 52−56.
- Антоненко С.В., Суров О. Э. Численное исследование прочности и качки судна на встречном волнении. // Тр. II междунар. конф по судостроению: ЦНИИ им. акад. А. Н. Крылова, СПб., 1998. Т 1. — С. 409 — 416.
- Ашик В. В Проектирование судов. Л.: Судостроение, 1975. — 352 с. Ю. Бабаев H.H., Лентяков В. Г. Некоторые вопросы общей вибрации судов. — Л.: Судпромгиз, 1961. — 308 с.
- П.Банди Б. Методы оптимизации. Вводный курс. -М.: Радио и связь.1988, — 128 с.
- Барабанов Н.В. Конструкция корпуса морских судов. В двух томах. Том 1. Общие вопросы конструирования корпуса судна. СПб.: Судостроение, 1993.-304 с.
- Барабанов Н.В. Конструкция корпуса морских судов. В двух томах. Том 2. Местная прочность и проектирование отдельных корпусных конструкций судна. СПб.: Судостроение, 1993. — 336 с.
- М.Барабанов Н. В., Иванов H.A., Новиков В. В., Окишев В. А., Чибиряк И. М. Повреждения судовых конструкций. Л.: Судостроение, 1977. — 400 с.
- Барабанов Н.В., Иванов H.A., Новиков В. В., Шемендюк Г. П. Повреждения и пути совершенствования судовых конструкций. Л.: Судостроение, 1989,-254 с.
- Барабанов Н.В., Кустов В. Н. Принципы, вызывающие отрыв носовой оконечности судов и ведущие к их гибели. // Сб. докладов междунар. конф., посвященной 100-летию ЦНИИ им. акад. А. Н. Крылова. СПб, 1994.-С. 227−233.
- Басин A.M. Качка судов. -М.: Транспорт, 1969. 272 с.
- Бельгова М. А. Изгибающие моменты для судов внутреннего плавания на волнении. Л.: Судостроение, 1966. — 208 с.
- Беляк Ю.Л. Экспериментальное исследование прочности корпусов судов. Л.: Судостроение, 1964. — 230 с.
- Благовещенский С. Н., Холодилин А. Н. Справочник по статике и динамике корабля. Л.: Судостроение, 1976. Т — 2, — 176 с.
- Бойцов Г. В. К вопросу об определении внешних нагрузок при слемин-ге. // Судостроение, 1985. N 5, с. 30.
- Бойцов Г. В. Оптимизация судового корпусах учетом требований снижения его металлоемкости и трудовых затрат// Судостроение, 1984. N 3,
- Бойцов Г. В., Палий О. М. Прочность и конструкция корпуса судов новых типов. Л.: Судостроение, 1979, 360 с.
- Болотин В.В. Статистические методы в строительной механике. М.: Строииздат, 1965, 279 с.
- Бородай И.К., Нецветаев Ю. А. Качка судов на морском волнении. -JL: Судостроение, 1969. 432 с.
- Бородай И.К., Нецветаев Ю. А. Мореходность судов: Методы оценки. -Л.: Судостроение, 1982. ч
- Бронников A.B. Особенности проектирования морских транспортных судов. Л.: Судостроение, 1971. — 328 с.
- Бронский А.И. Корпусные конструкции судов промыслового флота. Л.: Судостроение, 1978, 198 с.
- Бубнов И.Г. Строительная механика корабля. В 2 частях. ЧЛ, 1912. 4.2,1914.
- Васильев Ф.П. Численные методы решения экстремальных задач. -М.: Наука, 1980.- 518 с.
- Васильева Э.Ю. Влияние основных проектных характеристик на качку судна в условиях нерегулярного волнения. Проектирование и конструкции судов: Автореф. дис. канд. техн. наук. Калининград, 1992. — 13 с.
- Ветер и волны в океанах и морях: Справочные данные. Л.: Транспорт, 1974.-С. 182−192.
- Внешние силы, действующие на суда: ЦНИИ РУМБ. Л., 1976. — 120 с.
- Вознесенский А.И. Вероятностный расчет перемещений, скоростей и ускорений различных точек корабля на нерегулярном морском волнении. Доклад к IX науч. техн. конф. НТО Судпрома по теории корабля. Л., 1958.
- Галахов И.Н., Кандель Ф. Г., Фридлянский А. З. Исследование изгибающих и крутящих моментов в модельном эксперименте на регулярном волнении. Материалы по обмену опытом- - Тр. НТО Судпрома, — JI, 1972. -Вып. 189. — С. 119−130.
- Герасимов A.B. Энергостатистическая теория нелинейной нерегулярной качки судна. Д.: Судостроение, 1979.
- Герасимов A.B., Мореншильдт В. А., Нецветаев Ю. А. Разработка инженерной методики расчета продольной качки транспортных судов. В сб.: Экспериментальная гидродинамика судна (Тез. Докладов к VIII науч. — техн. конф.). — Л., НТО Судпрома, 1966. Вып. 80.
- Екимов. В.В. Методы теории вероятностей в кораблестроении. Л.: Судостроение, 1970. — 272 с.
- Екимов. В.В. Вероятностные методы в строительной механике корабля. Л.: Судостроение, 1966. — 328 с.
- Ершов Н.Ф., Попов H.A. Прочность судовых конструкций при локальных динамических нагружениях. Л.: Судостроение, 1989. -200 с.
- Ершов Н.Ф., Шахверди Г. Г. Метод конечных элементов в задачах гидродинамики и гидроупругости. Л: Судостроение, 1984. — 237 с.
- Иванов H.A. Экспериментальное исследование экстремальных волновых нагрузок, действующих на судовые конструкции. // Сб. докладов между-нар. конф., посвященной 100-летию ЦНИИ им. акад. А. Н. Крылова. СПб, 1994. С.415−421.
- Иванов H.A., Каленчук C.B. Анализ условий, сопровождающих катастрофические повреждения судовых конструкций. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1996. С. 38 — 44.
- Ипатовцев Ю.Н., Короткин Я. И. Строительная механика и прочность корабля. Л.: Судостроение, 1991. — 288 с.
- Картузова Т.А., Ростовцев Д. М. Оценка волновых и вибрацтонных изгибающих моментов, действующих на корпус судна в услових нерегулярного волнения. // Строительная механика и прочность судовых конструкций: сб. науч. Тр./Л., ЛКИ, 1981.-С. 41 48.
- Кондриков Д.В., Четыркин Н. В. Применение статистического метода при оценке общей прочности судна по результатам эксперимента. -Тр.//ЦНИИМФ Л.: 1962. — Вып.41. — С. З — 22.
- Корн Г., Корн Т. Справочник по математике. -М.: Наука, 1968.-720 с.
- Королев В.И. Экспериментальные исследования величины изгибающих моментов судов на волнении. // Судостроение, 1960. № 1.
- Короткин Я.И., Рабинович О. Н., Ростовцев Д. М. Волновые нагрузки корпуса судна. Л.: Судостроение, 1987. — 236 с.
- Короткин Я.И., Ростовцев Д. М., Сивере Н. Л. Прочность корабля. Л.: Судостроение, 1974. — 432 с.
- Короткин. Я.И. Вопросы прочности морских транспортных судов. -Л.: Судостроение, 1965. 390 с.
- Крылов А.Н. Избранные труды.-Л.: Академия наук СССР, 1958.-804 с.
- Курдюмов А. А Вибрация корабля. Л.: Судпромгиз, 1961. — 320 с.
- Лазарев В.Н., Юношева Н. В. Проектирование конструкций судового корпуса и основы прочности судов. — Л.: Судостроение, 1989. — 210 с.
- Луговский В.В. Гидродинамика нелинейной качки судов. Л.: Судостроение, 1982. — 250 с.
- Луговский В.В. Динамика корабля. Л.: Судостроение, 1976. — 199 с.
- Лузянин A.A. Исследование особенностей продольной качки и дополнительного сопротивления морских составных судов. Теория корабля: Авто-реф. дис. канд. техн. наук. С.-Петербург, 1997. — 22 с.
- Любушин Н.П. Экономическая эффективность проектных решений в судокорпусостроении. Л.: Судостроение, 1982.
- Максимаджи А. И. Амплитудно-частотные характеристики волновых изгибающих моментов. // Тр. ЦНИИМФ. 1971. — Вып. 134. — С. 11−32.
- Максимаджи А.И. Экономические факторы при нормировании продольной прочности корпусов судов// Судостроение, N 12, 1990. c. l 1.61 .Максимаджи А. И. Прочность морских транспортных судов. Д.: Судостроение, 1976. — 311 с.
- Нецветаев Ю. А. Анализ результатов испытаний моделей для исследования волновых моментов в корпусах транспортных судов на встречном волнении. // Тр. ЦНИИ им. А. Н. Крылова. 1968. — Вып.245. — С.77−93.
- Нецветаев Ю. А. Исследование распределения гидродинамических давлений по корпусу корабля при движении на встречном волнении. ЦНИИ им. А. Н. Крылова. Диссертация, 1958.
- Ногид JIM. Остойчивость судна и его поведение на взволнованном море. Л.: Судостроение, 1967. — 244 с.
- Ногид Л.М. Проектирование морских судов. Л.: Судостроение, 1976. — 208 с.
- Ногид Л.М. Проектирование формы судна и построение теоретического чертежа. Л.: Судпромгиз, 1962.
- Нормы прочности морских судов. Регистр СССР. 1991, 92 с.
- Перельмутер A.C. Экспериментальное исследование влияния демпфирования на характеристики килевой качки. ПММ, т. X, 1946.
- Постнов В.А. Численные методы расчета судовых конструкций. Л.: Судостроение, 1983. — 248 с.
- Постнов В.А., Тарануха H.A. Оценка напряженно-деформированного состояния корпуса судна методом модуль-элементов. // Судостроение, 1983. № 5, С. 5 — 8.
- Правила классификации и постройки морских судов. Регистр Судоходства, СПб., 1995.
- Прикладные задачи динамики судов на волнении. Под ред. И.К. Боро-дая Л.: Судостроение, 1989. — 220 с. 73 .Продольная качка и слеминг лихтеров лихтеровозной системы. Черноморское ЦПКБ. Руководитель Я. М. Элис.-№ Гр 77 030 573−0десса, 1977.-20 с.
- Прочность судов внутреннего плавания. Справочник. М.: Транспорт, 1978.
- Путов Н.Е. Проектирование конструкций корпуса морских судов. Часть 1. Л.: Судостроение, 1976. — 376 с.
- Путов Н.Е. Проектирование конструкций корпуса морских судов. Часть 2. Л.: Судостроение, 1977. — 424 с.
- Разработка алгоритмов и программ расчета изгибающих и крутящих моментов, действующих на корпуса судов на волнении. ЛКИ. Руководитель Д. М. Ростовцев. № Гр У 136 881 — Л., 1976.- 178 с.
- Разработка метода расчета сил и моментов, действующих на судно на регулярном волнении. НКИ. Руководитель Ю. В. Ремез. N Гр 75 018 950 — Николаев, 1976.- 79 с.
- Разработка теоретических методов расчета волновых нагрузок. ЦНИИМФ Руководитель А. И. Максимаджи. № Гр 76 050 702 — Л., 1976, — 98 с.
- Рахманин H.H. Актуальные проблемы мореходности судов. // Судостроение, 1991. N4, С. 3−5.
- Ремез Ю.В. Качка корабля. Л.: Судостроение, 1983. — 327 с.
- Ростовцев Д. M. К вопросу определения волнового изгибающего момента. // Тр. ЛКИ. 1960. — Вып.27.
- Салькаев А.З. Метод определения гидродинамических коэффициентов в уравнениях продольной качки судов на волнении. // Судостроение, 1978.1. N12,-С. 8−13.
- Семенов-Тян-Шанский В.В., Благовещенский С. Н., Холодилин А. Н. Качка корабля. Л.: Судостроение, 1969. — 392 с.
- Сердюченко А.Н. Анализ аварий и катастроф судов в море. // Тр. между нар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1999. С. 155 — 160.
- Сердюченко А.Н. Исследование волновых нагрузок на корпуса судов с учетом групповой структуры морского волнения. Строительная механика корабля: Автореф. дис. канд. техн. наук. Николаев, 1980. — 27 с.
- Сердюченко А.Н. Нелинейные эффекты в расчетах гидродинамических нагрузок на судовые конструкции. Л.: ЦНИИ «Румб». — 1984. — Вып.40. -С.78−84.
- Сердюченко А.Н., Федоренко Г. В. Расчет нелинейных гидродинамических нагрузок на судовые корпуса при качке на волнении. // Тр. междунар. конф. Кораблестроение и океанотехника. Проблемы и перспективы: ДВГТУ, Владивосток, 1998. С. 181−186.
- Сердюченко А.Н. Обеспечение безопасности судна на волнении в катастрофических условиях. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1996. С. 16−21.
- Слободян С.О. Исследование волновых нагрузок судов на регулярном волнении умеренной крутизны. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1996. С.30−37.
- Слободян С.О. Оценка влияния некоторых факторов на экстремальные волновые нагрузки. // Тр. междунар. конф. Кораблестроение и океанотехника. Проблемы и перспективы: ДВГТУ, Владивосток, 1998. С. 386−390.
- Справочник по вероятностным расчетам.-М.: Воениздат, 1970. 536 с.
- Справочник по строительной механике корабля./ Бойцов Г. В., Палий О. М., Постнов В. А., Чувиковский B.C. В трех томах .- Л.: Судостроение, 1982.
- Справочник по строительной механики корабля, т 3. Под ред. Ю. А. Шиманского. Л., Судпром гиз, 1960. — 800 с.
- Справочник по теории корабля: В трех томах. Статика судов. Качка судов. Под ред. Я. И. Войткунского. Л.: Судостроение, 1985.
- Суров О.Э. Влияние нагрузки судна на продольную качку на регулярном волнении. // Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1996. С.57−61.
- Суров О.Э. Исследование волновых изгибающих моментов на регулярном косом волнении.// Материалы II международной конф. «Проблемы транспорта Дальнего Востока», Академия транспорта РФ, Владивосток, 1997. -С.83−84.
- Суров О.Э., Корепанова Ю. А. Расчетное и экспериментальное исследование волновых изгибающих моментов. // Молодеж и научно-технический прогресс: Материалы конф. ДВГТУ, Владивосток, 1998. — С. 51 — 52.
- Суслов В.П. О расчетах нагрузок, действующих на суда в экстремальных волновых условиях. // Строительная механика корабля: Сб. науч. тр. -Николаев, НКИ, 1986. С.3−11.
- Суслов В.П. Об экстремальных волновых нагрузках, действующих на корпуса судов на морском волнении. // Строительная механика корабля: Сб. науч. тр НКИ, 1978. Вып. 136. — С. 15−22.
- Суслов В.П., Паплаускас А. Р. Алгоритм и ЭВМ программа для расчета нелинейной продольной качки и волновых нагрузок на экстремальном волнении.// Строительная механика корабля: Сб. науч. тр. — Николаев, НКИ, 1990.-С.З-16.
- Суслов В.П., Суслов C.B. Нелинейный метод расчета волновых нагрузок.// Тр. междунар. конф. Проблемы прочности и эксплуатационной надежности судов: ДВГТУ, Владивосток, 1996. С.8−15.
- Суслов В.П., Суслов C.B. О вероятностно детерминированном методе определения волновых нагрузок, действующих на суда. // Строительная механика корабля: Сб. науч. тр. — Николаев, НКИ, 1994.
- Ю8.Сухир Э. Л. Амплитудно-частотные характеристики продольной качки и волновых моментов для гибких судов. // Регистр СССР: Сб. науч. тр. Л., Транспорт, 1978. — Вып.8. — С.49- 67.
- Ю9.Сухир Э. Л. Исследование изгибающих моментов при балластном плавании судов на волнении. Строительная механика и вибрация корабля: Ав-тореф. дис. канд. техн. наук. Одесса, 1967. — 14 с.
- ПО.Фирсов Г. А. Развитие идей академика А. Н. Крылова в области теории и расчета кораблей на волнении. Тр. НТО Судпрома, 1957. Т. — 7, вып. 2.
- Хаскинд М.Д. Гидродинамическая теория качки корабля. М.: Наука, 1973.- 327 с.
- Химмельблау Д. Прикладное нелинейное программирование. М.: Мир, 1975.- 537 с.
- Холодилин А.Н. Стабилизация судна на волнении. Л.: Судостроение, 1973. — 230 с.
- Цындаря И.И. Прогиб корабля при нелинейной качке. // Тр. ВМАКВ им. акад. А. Н. Крылова. Л., 1951. Вып. 4, 5.
- Четыркин Н.В. Об условиях расчета характеристик внешних сил при вероятностной оценке общей прочности корпуса судна. //Тр. ЦНИИМФ Л.: 1962. — Вып.41. — С.23- 42.
- Пб.Чижиумов С. Д. Исследование динамики слеминга судна с помощью численных моделей. // Тр. междунар. конф. Кораблестроение и океанотехника. Проблемы и перспективы: ДВГТУ, Владивосток, 1998. С. 410−414.
- Чижиумов С.Д. Численные алгоритмы определения присоединенных масс воды в задачах вибрации и качки судов. Там же С. 415−419.
- Чувиковский Г. С. Исследование величин внешних сил, действующих на корпус судна в условиях удара о встречную волну. -Тр./НТО Судпрома, 1960, N351.
- Шаров Я.Ф. Об ударе судна днищем о встречную волну. // Тр. НТО СП, 1960. -Вып.35.
- Шмырев А.Н., Морендшильдт В. А., Ильина С. Г., Гольдин А. И. Успокоители качки судов. Д.: Судостроение, 1972. — 480 с.
- Элис Я.М. Исследование гидродинамических характеристик уравнений качки судов. Теория корабля: Автореф. дис. канд. техн. наук. Одесса, 1970.- 22 с.
- Юрков Н.Н. Некоторые вопросы мореходности судов на встречном волнении. Теория корабля: Автореф. дис. канд. техн. наук. Д., 1974. — 20 с.
- 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.
- 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.
- 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.
- Fukuda J. Statistic Prediction of Ship Response. Society of Naval Architects of Japan, Second Symposium on Ship Manoeuverability, July 1969.
- 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.
- Kusumoto H. et al. Ship Hull Structure Design by DLA (Dynamic Loading Approach), IHI Engineering Review, Vol.27, No.2, April 1994.
- Lewis E. V. Ship. Model Test to Determine Bending Moment in Waves. SNAME, Vol. 62, 1954. pp. 431−490.
- 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.
- Lotveit M. and Oth. Wave Loads on a T-2 Tanker Model. «European Shipbuilding», 1961. N 10, — pp. 31−53.
- Moor D. Longitudinal Bending Moments on Models in Head Seas. TRINA, 1967.-pp. 117−166.
- 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.
- 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.
- Salvesen N., Tuck E.O. and Faltinsen O. Ship Motions and Sea Load. TSNAME, Vol. 78, 1970.
- 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.
- 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.
- Vossers G. Behavior of Ship in Waves. Harlem, 1962.
- Vossers G., Swaan W. A. and Rijiken H., Experiments with Series 60 Models in Waves, International Shipbuilding Progress, May, 1961. p. 201.
- ОПИСАНИЕ ФУНКЦИЙ И ПОДПРОГРАММ
- DECLARE FUNCTION WINDS (tITLES, LEFT%, RIGHT%, UP%, DOWN%, INK%, PAPER%) DECLARE FUNCTION SAVESCR% (NAMES) DECLARE FUNCTION RESTSCR% (NAMES)
- DECLARE FUNCTION INTEGRAL. CEB! (X!(), Y!(), N. SHP!, KOL. TOCH!(), SUMMA!)
- DECLARE SUB SAVE. FILE (Y!(), X!(), KOL. TOCH!(), NAMESS)
- DECLARE SUB SPLAIN (X!(), Y!(), OSADKA!, N. SHP!, X!, Y!, KOL. TOCH!())
- DECLARE SUB RED (X!(), Y!(), KOL. TOCH!())
- DECLARE SUB OPEN. FILE (Y!(), X!(), KOL. TOCH!(), NAMESS)
- DECLARE SUB RED. TOCH (X!(), Y!(), 1!, N. SHP!, KOL. TOCH!())
- DECLARE SUB RED. SHP (X!(), Y!(), KOL. TOCH!(), N. SHP!)
- DECLARE SUB raspred. nagruzki ()
- DECLARE SUB zalivaemost ()
- DECLARE SUB runl5 (lv!, SIG!, ka!)
- DECLARE SUB RUN. NONLINE (ZETA, FI, ZETAD, FID, X)
- DECLARE SUB FON1 (IN!, ik!, jn!, jk!)
- DECLARE SUB message (S$, C$)1. DECLARE SUB shapkamenu ()
- DECLARE SUB RAMKA (XN!, XK!, YN!, YK!, sim!, FOON!)
- DECLARE SUB Spravka (CS$, DF$)1. DECLARE SUB oshibka (ty!)1. DECLARE SUB sapkatab ()
- DECLARE SUB star (NAM!, R!)1. DECLARE SUB REDTOTAL ()1. DECLARE SUB neregulyr ()
- DECLARE SUB prnscreen (wr!)1. DECLARE SUB SAVE ()
- DECLARE SUB printdat (AS (), NAM!, R!)1. DECLARE SUB openfile ()1. DECLARE SUB Regulyr ()1. DECLARE SUB NEWFILE ()
- DECLARE SUB message 1 (n$, S$)
- DECLARE SUB qlstn (X!, Y!, ST$, B!)
- DECLARE SUB submenu (NAM!, R!)1. DECLARE SUB UPMENU (NAM!)
- DECLARE SUB MENU (X!, Y!, SH!, n!, S$(), NAM!)1. DECLARE SUB ZASTAV4 ()1. DECLARE SUB ZASTAV2 ()
- ОПИСАНИЕ ПЕРЕМЕННЫХ И КОНСТАНТ
- НАЧАЛО ГЛАВНОЙ ПОДПРОГРАММЫ1. H% = 201. CALL ZASTAV2 WIDTH 80, 5 010 000 CALL FONl (l, 48, 1, 80): PCOPY 0, 110 001 CALL star (NAM, R)77 ty = ERR: CALL oshibka (ty) PCOPY 1,0 IF ty = 11 THEN
- ПОДПРОГРАММЫ *************1. РАСЧЕТ КПМ И КД
- SUB KPMKD (ka, ТТ (), ВВ (), ВАТАА (), NAP, BP, I, М33(), М330(), L33(), L330(), SIG) DIM BBB (4, 2, 2), A (4, 2, 2) BBB (0, 0
- M33(I) = MUKI * BB (I)A 2 * 3.141 593 * RO/8 + NAP * RO * BP M330(I) = MUI * BB (I) A 2 * 3.141 593 * 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
- SUB neregulyr 'РАСЧЕТ НА НЕРЕГУЛЯРНОМ ВОЛНЕНИИ
- BES = 2.5: ALS = .21 * BES
- AS2 = ALS л 2 BES A 2: BS2 = ALS л 2 + BES A 2
- 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)
- 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
- CALL RAMKA (xn2, хк2, ynl, ykl, sim, fonn) COLOR sim, fonn
- CATE ynl + 1, xn2 + 2: PRINT «ВАРИАНТ II
- CATE ynl + 3, xn2 + 2: PRINT «Высота волны обеспеч. 3%, см «
- CATE ynl + 4, xn2 + 2: PRINT «Обеспеченность волнения в % «
- CATE ynl + 5, xn2 + 2: PRINT «Номер расчетного сечения «хпЗ = 20: хкЗ = 60: упЗ = 30: укЗ = 38
- CALL RAMKA (xn3, хкЗ, упЗ, укЗ, sim, fonn)1. COLOR sim, fonn
- CATE YN + 1, XN + 2: PRINT «Для расчета продольной качки судна»
- CATE YN + 2, XN + 2: PRINT «на нерегулярном волнении необходимо»
- CATE YN + 3, XN + 2: PRINT «ввести дополнительные данные по трем,»
- CATE YN + 4, XN + 2: PRINT «ниже перечисленным вариантам «
- XN = 20: ХК = 60: YN = 42: YK = 46
- CALL RAMKA (XN, ХК, YN, YK, sim + 16, 6)1. COLOR sim, 61 444 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
- CATE ynl + 3, xnl + 34: INPUT»», BAL LOCATE ynl + 4, xnl + 34: INPUT»», obes
- CATE ynl + 5, xnl + 34: INPUT ««, sechenie
- CATE ynl + 3, xn2 + 35: INPUT ««, AH3: LOCATE ynl + 4, xn2 + 35: INPUT ««, obes
- CATE ynl + 5, xn2 + 35: INPUT ««, sechenie: AH = AH3 /1.32: TV = 5.9 * (ah /100) л (1 /
- TV = 3.3 * SQR (AH3 /100): END IF: IF variant. dat = 3 THEN: COLOR sim, fonn
- CATE yn3 + 3, хпЗ + 35: INPUT ««, AH3: LOCATE yn3 + 4, хпЗ + 35: INPUT ««, TV
- CATE yn3 + 5, хпЗ + 35: INPUT ««, obes: LOCATE yn3 + 6, хпЗ + 35: INPUT ««, sechenie
- AH = AH3 /1.32: END IF: xnl1 = 1: xkl1 = 17: ynl 1 = 30: ykl 1 = 38
- CALL RAMKA (xnl 1, xkl 1, ynl 1, ykl 1,14, 4): COLOR sim, 4
- CATE ynl 1 + 2, xnl 1 + 2: PRINT «1. СПЕКТР»
- CATE ynl 1 + 4, xnl 1 + 2: PRINT «РЕКОМЕНДУЕМЫЙ»
- 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
- CATE ynl 1 + 2, xnl2 + 2: PRINT «2. СПЕКТР»
- CATE ynl 1 + 4, xnl2 + 2: PRINT «РЕКОМЕНДУЕМЫЙ»
- CATE ynl 1 + 6, xnl2 + 2: PRINT «2 МКПС"1155 XN= 13: XK = 68: YN = 42: YK = 46
- CATE 45, 20: PRINT «Шаг по частоте"-» ««-
- CATE 46, 18: PRINT «ОТНОСИТЕЛЬНАЯ ЧАСТОТА ОПРЕДЕЛЯЕТСЯ КАК: vL/Lv»
- CATE 43, 60: INPUT ««, SIGBN
- CATE 44, 60: INPUT ««, SIGBKON
- CATE 45, 3 5: INPUT ««, stepsig1. SCREEN 12: WIDTH 80, 60
- WINDOW (0, -0)-(1000, 1000)
- CATE 1,1: PRINT «Спектр «
- CATE 2, 1: PRINT «волнения «
- CATE 1,20: PRINT» АМПЛИТУДНО ЧАСТОТНЫЕ ХАРАКТЕРИСТИКИ» LOCATE 2, 20: PRINT
- CATE 3, 20: PRINT «вертикальной"1.CATE 4,20: PRINT «качки»
- CATE 3, 40: PRINT «килевой»
- CATE 4, 40: PRINT «качки»
- CATE 3,60: PRINT «изгибающего»
- CATE 4, 60: PRINT «момента»
- CATE 5, 60: PRINT «для сечения «- sechenie
- CATE 34, 20: PRINT «СПЕКТРЫ ПРОЦЕССОВ"1.CATE 35,20: PRINT
- 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 ««
- CATE 30, 28: PRINT «Относительная частота vL/Lv»
- 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 ««-
- SZ / 300 > MSZ THEN MSZ SZ / 300: SIGB = SIGB + stepsig: WEND: SIGB = SIGBN
- WHILE SIGB ≤ SIGBKON: SIG = SIGB * (2 * 3.141 593 * 9.81 / LS) Л .5: GOSUB 20 000
- NE (20 + (SIGB stepsig / 5 — SIGBN) * 225 / (SIGBKON — SIGBN), SZ / MSZ + 550)-(20 +
- SIGB + stepsig / 5 SIGBN) * 225 / (SIGBKON — SIGBN), 550), 15, BF
- SIGB = SIGB + stepsig: WEND: COLOR 15
- CATE 7, 5: PRINT «Smax =" — CLNG (MSZ * 30 000) / 1001. CATE 8, 5: PRINT» мп. с»
- COLOR 14: LOCATE 42,1: PRINT «ФОРМА СПЕКТРА»
- 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
- 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спектральная плотность изг. момента от (х)
- S.I.M (CH) = ACH.I.M (CH) А 2 * SZ: S.I.M = S.I.M (CH) + S.I.M: pfim = FIM * ka
- SZG (CH) = ZETAM Л 2 * SZ: SZG = SZG + SZG (CH)
- SZ / 300 > MSZ THEN MSZ = SZ / 300
- ZETAM / 300 > MZET THEN MZET = ZETAM / 300
- FIM / 300 > MFIM THEN MFIM = FIM / 300
- S.I.M (CH) / 300 > M.SI.M THEN M.SI.M = S.I.M (CH) / 300
- ACH.I.M (CH) / 300 > M.I.M THEN M.I.M = ACH.I.M (CH) / 300
- SZG (CH) / 300 > MSZG THEN MSZG = SZG (CH) / 300
- SAG (CH) / 300 > MAZG THEN MAZG SAG (CH) / 300
- CH.= CH + 1: SIGB = SIGB + stepsig: WEND
- COLOR 15: LOCATE 7, 5: PRINT «Smax =" — CLNG (MSZ * 30 000) / 100 LOCATE 8, 5: PRINT «мп.с»
- CATE 7, 25: PRINT «Zmax =" — CLNG (MZET * 30 000) / 100 LOCATE 7, 45: PRINT «Fmax CLNG (MFIM * 30 000) / 100 LOCATE 7, 65: PRINT «Mmax =" — CLNG (M.I.M * 30 000) / 1001. CATE 8, 65: PRINT» т»
- CATE 37, 25: PRINT «Szmax =" — CLNG (MSZG * 30 000) / 100 LOCATE 38, 25: PRINT «мп-с»
- CATE 37, 43: PRINT «Sfmax =" — CLNG (MAZG * 30 000) / 100
- CATE 38, 43: PRINT» градус»
- CATE 37, 61: PRINT «Sm =" — M.SI.M * 300
- CATE 38, 63: PRINT» (м-т)п.с»
- SIGB = SIGBN: CH = 1: 'ВЫВОД ГРАФИКОВ НА ЭКРАН
- WHILE SIGB ≤ SIGBKON: SIG = SIGB * (2 * 3.141 593 * 9.81 / LS) л .5: GOSUB 20 000
- NE (20 + (SIGB stepsig / 5 — SIGBN) * 225 / (SIGBKON — SIGBN), SZ / MSZ + 550)-(20 +
- SIGB + stepsig / 5 SIGBN) * 225 / (SIGBKON — SIGBN), 550), 15, BF
- NE (270 + (SIGB stepsig / 5 — SIGBN) * 225 / (SIGBKON — SIGBN), SZG (CH) / MSZG + 50)270 + (SIGB + stepsig / 5 SIGBN) * 225 / (SIGBKON — SIGBN), 50), 15, BF
- NE (520 + (SIGB stepsig / 5 — SIGBN) * 225 / (SIGBKON — SIGBN), SAG (CH) / MAZG + 50)520 + (SIGB + stepsig / 5 SIGBN) * 225 / (SIGBKON — SIGBN), 50), 15, BF
- 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)
- FGG (CH) = 550 + FGG (CH) / MZET: FGGl (CH) = 550 + FGGl (CH) / MFIM1. CH = 1 THEN
- PSET (fg (CH), FGG (CH)), COLLINE
- PSET (fg (CH) + 250, FGGl (CH)), COLLINE
- PSET (fg (CH) + 500, 550 + ACH.I.M (CH) / M.I.M), COLLINE1. ELSE
- NE (fg (CH), FGG (CH))-(fg (CH 1), FGG (CH — 1)), COLLINE
- NE (fg (CH) + 250, FGG1 (CH))-(fg (CH 1) + 250, FGG1(CH — 1)), COLLINE
- 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
- CATE 40, 5: PRINT «СТАНДАРТЫ:»
- CATE 42, 1: PRINT «Вертикальной качки»
- CATE 43, 1: PRINT SQR (dszg) — «м»
- CATE 45, 1: PRINT «Килевой качки»
- CATE 46, 1: PRINT SQR (dsag) — «град»
- CATE 48, 1: PRINT «Изгибающего момента»
- 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
- 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
- BEEP: CALL message («ncmr0T0BTe принтер!», «Продолжить-<�Еп1ег> Отменить-<�Езс>»)wp$ = INKEY$:1. wp$ = CHR$(27) THEN
- PCOPY 1,0: GOTO 123: END IF: END IF123 END IF
- 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 20 000 'ВЫБОР СПЕКТРА ВОЛНЕНИЯ IF variant. spektra = 1 THEN '12 МКОБ
- SZ = (.7795 / SIG A 5) * EXP (-3.11 * 10 A 4 / (AH A 2 * SIG A 4)) 'волновой спектор END IF1. variant. spektra = 2 THEN '2 МКПС
- SZ = 9.43 * .3 575 * (AH3 /100)A 2 / (2 * 3.141 593 / TV) * ((.777 * 2 * 3.141 593 / TV) / SIG) A6 * EXP (-1.5 * ((.777 * 2 * 3.141 593 / TV) / SIG) A 4)
- 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) *
- SIG A 4 + 2 * SIG A 2 * AS2 + BS2 A 2))
- END IF: RETURN: 90 END SUB
- SUB NONLINE (ZATA, FI)' НЕЛИНЕЙНЫЙ РАСЧЕТ КАЧКИ И МОМЕНТОВ
- DIM Zl (-1 ТО Н% + 1), QQ (-1 ТО Н% + 1), QQI (-1 ТО Н% + 1), ММ (-1 ТО Н% + 1), PIR. S (-1
- ТО Н% + 1), IZG. M (-1 ТО Н% + 1)
- SCREEN 12: WIDTH 80, 60: WINDOW (0, -0)-(1000, 1000)
- SHARED All, Bll, Cll, А12, В12, С12, A21S, B21S, C21S, A22S, B22S, C22S, FF, MM, QN, 1. QNX, Fl, F2, Ml, M2
- OPEN «non» FOR OUTPUT AS #100
- PRINT #100, «ВРЕМЯ, с"-» M (6)(T*M)"-» M (10)(T*M)"-» М (14)(Т*М)"-» ВЕРТИКЛЬНАЯ М"-» КИЛЕВАЯ РАД""РЕШЕНИЕ СИСТЕМЫ ДИФ. УРАВНЕНИЙ МЕТОДОМ РУНГЕ-КУТТА-ФЕЛЬБЕРГА"с АВТОМАТИЧЕСКИМ ВЫБОРОМ ШАГА"г**************** ЧИСЛО УРАВНЕНИЙ N
- CLS: п = 4: DIM YY (n), L (n), W (n), A (n), BB (n), C (n), D (n), E (n), F (n)ka = 2 * 3.141 593 / lv: SIG = (ka * 9.81) A .5: CALL runl5(lv, SIG, ka)
- W (l) = ZETA1: W (2) = FI1: W (3) = ZETA2 * SIGk: W (4) = FI2 * SIGk
- El = .03 ,*************** ПОГРЕШНОСТЬ ВЫЧИСЛЕНИЙ El
- H = 2 * 3.141 593 / (SIGk * 20) <*************** НАЧАЛЬНЫЙ ШАГ H
- X = 0 .*************** НАЧАЛЬНОЕ Х0 X
- 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
- A (J) = F (J) * H: YY (J) = W (J) + 2 * A (J) / 9: NEXT J
- X = X + 2 * H / 9: GOSUB 400: FOR J = 1 TO n
- 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:
- C (J) = F (J) * H: YY (J) = W (J) + (34.5 * A (J) 121.5 * BB (J) + 135 * C (J)) / 64:
- NEXT J: X = X + H / 2.4: GOSUB 400: FOR J = 1 TO n
- D (J) F (J) * H: Q = W (J) — 17 * A (J) /12 + 6.75 * BB (J)
- YY (J) = Q 5.4 * C (J) + 16 * D (J) /15: NEXT J:
- X = X + H/4: GOSUB 400: FOR J = 1 TO n:
- E (J) = F (J) * H: Q = W (J) + 65 * A (J) / 432 .3125 * BB (J)
- YY (J) = Q + .8125 * C (J) + 4 * D (J) /27 + 5 * E (J) / 144: NEXT J
- X = X H / 6: GOSUB 400: FOR J = 1 TO n
- 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
- IF E2 < El / 20 THEN DDD = DDD + 1
- NEXT J: X = X + H / 6: GOSUB 400: IF E3 = 0 THEN 290
- X = X H: FOR J = 1 TO n: YY (J) = W (J): NEXT J1. H = H/2: GOTO 1 101 290 IF DDD = n THEN H = H + H1. X ≥ 40 THEN1.Ho.l THEN H =. 11. END IF
- 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
- DL = LS /H%: alfa# = (ALF / 180) * 3.141 592 654#: ka 2 * 3.141 593 / lv: KA1 = -ka * COS (alfa#) MAX. Q = 0: MAX. S = 0: MAX. M = 0: FOR I = 0 TO H%
- 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)
- YY (2) + FS (I) * COS (SIGk * X) FSS (I) * SIN (SIGk * X) + Q. N (I)
- 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%
- 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))
- 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
- I = 0 THEN: PSET (XA (I) * MMX, Z1(I) + MMZZ): PSET (XA (I) * MMX, X (l, I) + MMZZ)
- PSET (XA (I) * MMX, MMZ * X (KOL.TOCH (I) -1,1) + MMZZ): ELSE
- NE (XA (I 1) * MMX, Z1(I — 1) + MMZZ)-(XA (I) * MMX, Z1(I) + MMZZ), 1
- NE (XA (I 1) * MMX, MMZ * X (l, I — 1) + MMZZ)-(XA (I) * MMX, MMZ * X (l, I) + MMZZ)
- NE (XA (I 1) * MMX, MMZ * X (KOL.TOCH (I — 1) — 1,1 — 1) + MMZZ)-(XA (I) * MMX, MMZ
- X (KOL.TOCH (I) 1,1) + MMZZ), 4: END IF
- 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»
- 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%
- NE ((XA (I) DL * .5) * MMX, QQ (I) * MAX. Q + Z. NACH)-((XA (I) + DL * .5) * MMX, QQ (I) * MAX. Q + Z. NACH), 13
- NE ((XA (I) + DL * .5) * MMX, QQ (I) * MAX. Q + Z. NACH)-((XA (I) + DL * .5) * MMX, Z. NACH), 13
- NE ((XA (I) DL * .5) * MMX, QQ (I) * MAX. Q + Z. NACH)-((XA (I) — DL * .5) * MMX, Z. NACH), 13
- 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
- NE (XA (I) * MMX, IZG. M (I) * MAX. M + Z. NACH)-(XA (I + 1) * MMX, IZG. M (I + 1) *
- MAX.M + Z. NACH), 15: NEXT I
- X ≥ 40 THEN: PRINT #100, USING «###.###" — X-
- PRINT #100, USING «#########.##" — IZG. M (6) — IZG. M (IO) — IZG. M (14) —
- PRINT #100, USING «########.####" — YY (1) — YY (2): END IF
- COLOR 15: LOCATE 15, 15: PRINT «КИНЕМАТИЧЕСКИЕ ХАРАКТЕРИСТИКИ1. ПРОДОЛЬНОЙ КАЧКИ»
- CATE 16, 5: PRINT «Время T=" — USING «###.##" — X-: PRINT «с» LOCATE 16,25: PRINT «Вертикальная»
- 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 «Килевая»
- 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)
- F (4) = (-B22S * YY (4) C22S * YY (2) — A21S * F (3) — B21S * YY (3) — C21S * YY (1) + MM) RETURN: END END SUB
- SUB polnover' ПОЛНОВЕРОЯТНОСТНЫЙ РАСЧЕТ КАЧКИ И МОМЕНТОВ
- OPEN «polnover» FOR OUTPUT AS #44 sechenie =10: obes = 3
- 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
- DM = .025 * (GAMMA * LS л 2 * Y (KOL.TOCH (H% 2) 2, H% 2) * 2): DZZ = 2.5: DFF = 2.5 'SM (ALFB, TB, HB)
- DIM SW (30, 10, 13), ACH. M (20, 12), ACH. H (20, 12), ACH. P (20, 12)'АЧХ и сектор волнения
- DIM AH. M (20, 12, 5), AH. H (20, 12, 5), AH. P (20, 12, 5)
- DIM ACHM (20, 12, 5), АСНЩ20, 12, 5), ACHP (20,12, 5)
- DIM SM (12, 10, 13), SH (12, 10, 13), SP (12, 10, 13) 'стандарты
- DIM PALF (12), PAH (3, 13), PTE (8), PLOAD (5)
- DIM PUM (IO), PTU (13), PRAION (3), QWM (IO), QHT (10, 4), SPEED (5)
- DIM PUH (IO), PTUHH (13), QWH (10), QHTHH (10, 4)
- DIM PUP (10), PTUPP (13), QWP (10), QHTPP (10, 4): NER= I
- ВЕРОЯТНОСТИ РАЙОНОВ 'ВЕРОЯТНОСТИ НАГРУЗКИ 'ВЕРОЯТНОСТИ ПЕРИОДОВ В ДОЛЯХ ОТ СРЕДНЕГО
- PRAION (l) PRAION (2) PRAION (3)3333 .3333 .3333
- PLOAD (O) = .2: PTE (0) = .075 PLOAD (l) = .2: PTE (1) = .124 PTE (2) = .191 PTE (3) = .219 PTE (4) = .191
- PLOAD (2) = .2 PLOAD (3) = .2: PLOAD (4) = .2:
- ВСТРЕЧОЕ '15 '30 '45 '60 '75 '90 '105 '120 '135 '150 '165 '180
- 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
- END IF: 'ВЫБОР СЕЗОНА: tITLE$ = «СЕЗОН» NAM = 1: n = 4: XX = 28: YY = 9: DIM CS$(n)
- CS$(1) = «ЗИМА «:CS$(2) = «ВЕСНА «: CS$(3) = «ЛЕТО
- CS$(4) = «ОСЕНЬ «: A = v. menu (tITLE$, CS$(), n, XX, YY, NAM, sim%, FON%)
- SEZON = NAM ZAGOLOVOK$ = tITLE$ + «-» + RTRIM$(CS$(NAM)) + «» + ZAGOLOVOK$
- ERASE CS$: 'ВЫБОР МОРЯ: tITLE$ = «MOPE»: NAM = 1: n = 6: XX = 46: YY = 12:1. DIM CS$(n)
- CATE 27, 30: COLOR 1, 4: PRINT «--------------------»
- ACHP (OM, ALFB, FRUDB) = ACHP (OM, ALFB, FRUDB) + FIM * ka * 57.29 578 * PLOAD (load) 'АЧХ килевой качки NEXT LOADB OM = OM +1 NEXT SIGB
- ALFB = ALFB + 1 NEXT ALF LOW = FRUDB
- FOR FRUDB = 0 ТО FRUDKB STEP STEPFRUDB 'цикл по скорости STEPFRUDB=1 FRUDKB=5
- 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ом = о
- 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)))
- 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)))
- FOR ALF = ALFN TO ALFK STEP STEP ALF alfa# = (ALF / 180) * 3.141 592 654# SM (ALFB, ТВ, HB) = 0: SH (ALFB, ТВ, HB) 0: SP (ALFB, ТВ, HB) — 0ом = о
- FOR SIGB = SIGBN TO SIGBKON STEP stepsig '===============dsig = stepsig * (2 * 3.141 593 * 9.81 / LS) A .5 '|dsigk = dsig dsig л 2 * v * COS (alfa#) / 9.8066 «j
- 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
- CATE 27, 30 + LO: COLOR 14, 4: PRINT «-»: NEXT АНЗ i******************* РАСЧЕТ ПОЛНОЙ ВЕРОЯТНОСТИ ******************** FOR MOMENT = 0 ТО NMOMENT 'ЦИКЛ ПО МОМЕНТАМ
- МО = DM * MOMENT: ZZO = DZZ * MOMENT: FFO = DFF * MOMENT QWM (MOMENT) = 0: QWH (MOMENT) = 0: QWP (MOMENT) = 0 FOR RAION = 1 TO NRAION 'ЦИКЛ ПО РАЙОНАМ ПЛАВАНИЯ
- QHT (MOMENT, RAION) = 0: QHTHH (MOMENT, RAION) = 0: QHTPP (MOMENT, RAION) = 0 HB = 1 'НАЧАЛО ЦИКЛА ПО АНЗ%
- 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
- TV (TB * STEPT + TNAC) * TS PUM (TB) = 0: PUH (TB) = 0: PUP (TB) = 0
- ALFB = 0 'НАЧАЛО ЦИКЛА ПО ALF
- PRINT #44, «PTU (HB)», AH3, PTU (HB)
- 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 АНЗ
- 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 'ЦИКЛ ПО МОМЕНТАМ
- МО = 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
- NE (100 + -LOG (QWM (MOMENT)) / LOG (IO) * XMM, M0 * ZMM + 100)-(100 +
- G (QWM (MOMENT + 1)) / LOG (IO) * XMM, Ml * ZMM + 100), 14
- NE (100 + -LOG (QWH (MOMENT)) / LOG (IO) * XMM, ZZ0 * ZMHP + 100)-(100 +
- G (QWH (MOMENT + 1)) / LOG (IO) * XMM, ZZ1 * ZMHP + 100), 15
- NE (100 + -LOG (QWP (MOMENT)) / LOG (IO) * XMM, FF0 * ZMHP + 100)-(100 +
- G (QWP (MOMENT + 1)) / LOG (IO) * XMM, FF1 * ZMHP + 100), 131. NEXT MOMENT14 ТРАФИК
- COLOR 14: LOCATE 1,1: PRINT «ОБЕСПЕЧЕННОСТЬ ВОЛНОВЫХ
- ВЕРТИКАЛЬНЫХ ИЗГИБАЮЩИХ МОМЕНТОВ»
- CATE 3,40 LEN (ZAGOLOVOK$) 2: PRINT ZAGOLOVOK$
- CATE 56, 9: PRINT «0 1 2 3 4 5 6 7 8 9 10 -Log Q" —
- CATE 2, 2: PRINT «Mb»: LOCATE 3, 2: PRINT «—»: LOCATE 4, 2: PRINT «pgBLa»
- CATE 6, 4: PRINT «0.250»: LOCATE 11,4: PRINT «0.225»:
- CATE 16, 4: PRINT «0.200»: LOCATE 21, 4: PRINT «0.175»
- CATE 26, 4: PRINT «0.150»: LOCATE 30, 4: PRINT «0.125»
- CATE 35, 4: PRINT «0.100»: LOCATE 40, 4: PRINT «0.075»
- CATE 45, 4: PRINT «0.050»: LOCATE 50, 4: PRINT «0.025»
- CATE 60, 10: PRINT «ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ «-: DO: LOOP WHILE INKEY$ о CHR$(27)
- SCREEN 0: WIDTH 80, 50: OPEN «rollingV + FIL$ FOR INPUT AS #1: CALL printdat (A$(), NAM, rr): CLOSE #1
- 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
- TV = (TB * STEPT + TNAC) * TS
- ALFB = 0 'НАЧАЛО ЦИКЛА ПО ALF
- 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
- PRINT #44, «PUM (TB)»: 'FOR ТВ = 0 TO 7 STEP 1 'НАЧАЛО ЦИКЛА ПО TV
- PRINT #44, ТВ, PUM (TB): 'NEXT ТВ
- 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 'ЦИКЛ ПО МОМЕНТАМ
- FOR RAION = 1 TO NRAION 'ЦИКЛ ПО РАЙОНАМ ПЛАВАНИЯ
- PRINT #44, MOMENT, RAION, QHT (MOMENT, RAION): NEXT RAION NEXT MOMENT: END SUB
- SUB REGULYR «РАСЧЕТ НА РЕГУЛЯРНОМ ВОЛНЕНИИ В ЛИНЕЙНОЙ ПОСТАНОВКЕ
- DIM Z1(H% + 1), ZMP (-1 ТО Н% + 1), ZMA (-1 ТО Н% + 1)
- Z1(I) = (DN (DN — DK) * I / H% + ZETAG (IE) + XA (I) * FIG (IE) / 57.29 578 — ZA * COS (KAl * XA (I) + SIGk * IE)) * MMZ: XA (I) = LS — DL * I IF I = H% OR I = 0 THEN
- 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
- NE (XA (I) * MMX, MMZ * X (l, I) + MMZZ)-(XA (I) * MMX, MMZ * X (KOL.TOCH (I) -1,1) + MMZZ), 4: ELSE
- NE (0, DK * MMZ + MMZZ)-(LS * MMX, DN * MMZ + MMZZ), 3
- CATE IE * 6 + 6, 59: PRINT USING «#" — IE: LOCATE IE * 6 + 6, 60: PRINT «•П/4»
- CATE IE * 6 + 6, 66: PRINT USING «###.###" — ZETAG (IE)
- CATE IE * 6 + 6, 74: PRINT USING «###.###" — FIG (IE)
- NE (710, MMZZ)-(1000, MMZZ): NEXT IE
- CATE 3, 58: PRINT «Положен.Вертик. Измен.»
- CATE 4, 58: PRINT «волны перем, м. угла, гр»
- NE (XA (I) * MMX, 5 + MMZ * X (l, I))-(XA (I) * MMX, 5 + MMZ * X (KOL.TOCH (I) 1,1)), 4
- ELSE: LINE (XA (I 1) * MMX, 5 + MMZ * X (l, I — 1))-(XA (I) * MMX, 5 + MMZ * X (l, I)), 4
- NE (XA (I 1) * MMX, 5 + MMZ * X (KO.L.TOCH (I -1) -1,1 — 1))-(XA (I) * MMX, 5 + MMZ *
- X (KOL.TOCH (I) 1,1)), 4: END IF: LINE (XA (I) * MMX, 10)-(XA (I) * MMX, 5)
- NE (0, 5 + DK * MMZ)-(LS * MMX, 5 + DN * MMZ), 3: IF I = 0 THEN 3
- NE (XA (I) * MMX, 5 + ZMP (I) * MMZ)-(XA (I -1) * MMX, 5 + ZMP (I -1) * MMZ)
- 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
- X (KOL.TOCH (I) 1,1) + MMZZ), 4: END IF
- 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»
- ОТНОСИТЕЛЬНЫЙ МОМЕНТ Mw. = Mw/(LA2 * В * Rw* gamma)
- COLOR 15: LOCATE 19, 15: PRINT «ХАРАКТЕРИСТИКИ НАГРУЗКИ»
- COLOR 13: LOCATE 20, 5: PRINT «Максим, погонная нагрузка на шпацию t/m Р =" — MAX. Q
- COLOR 11: LOCATE 21,5: PRINT «Амплитуда макс, перерезывающей силы, t R =" — MAX. S-1. P. S (10)
- 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
- 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
- NE ((XA (I) DL * .5) * MMX, PP * MAX. Q + Z. NACH)-((XA (I) + DL * .5) * MMX, PP * MAX. Q + Z. NACH), 13
- NE ((XA (I) + DL * .5) * MMX, PP * MAX. Q + Z. NACH)-((XA (I) + DL * .5) * MMX, Z. NACH), 13
- 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
- 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
- NE (XA (I) * MMX, I. M (I) * ABS (MAX) + Z. NACH)-(XA (I + 1) * MMX, I. M (I + 1) * ABS (MAX) + Z. NACH), 15
- SCREEN 0: WIDTH 80, 50: A% = RESTSCR%(«OPENF.MJG»): GGG = 0 END SUB
- SUB RUN. NONLINE (ZETA, FI, ZETAD, FID, IE)' НЕЛИНЕЙНЫЙ РАСЧЕТ КОФ. УР-НИЙ КАЧКИ
- DIM QQ (H% + 1), QQI (H°/o + 1), MM (H% + 1), HP (H% + 1), BP (H% + 1), NACH. B (H% + 1)
- DIM KSI11#(H% + 1), KSI21#(H% + 1), KSI51#(H% + 1), KSI12#(H% + 1), KSI22#(H% + 1),
- DIM sil 1#(H% + 1, 12), si21#(H% + 1, 12), si51#(H% + 1, 12), sil2#(H% + 1, 12), si22#(H% + 1)
- DIM KS11#(H% + 1), KS21#(H% + 1), KS51#(H% + 1), TT (H% + 1), BB (H% + 1), GG (H% + 1)
- DIM DZ (H% + 1), DZM3(H% + 1), DZM30(H% + 1), DZL3(H% + 1), DZL30(H% + 1)
- DIM M33K (H% + 1), M330K (H% + 1), L33K (H% + 1), L330K (H% + 1)
- DIM M33N (H% + 1), M330N (H% + 1), L33N (H% + 1), L330N (H% + 1)
- РАСЧЕТ ВСПОМОГАТЕЛЬНЫХ ВЕЛИЧИНka = 2 * 3.141 593 / lv: SIG = (ka * 9.81) A .5alfa# = (ALF/ 180) * 3.141 592 654#
- DL = LS / H%: KA1 = -ka * COS (alfa#): KA2 = ka * SIN (alfa#): I = 0: WHILE I ≤ H% XA (I) = DL * ((H% / 2) -1)
- NAP (I) = TT (I) + X (1,1) HP (I) TT (I) — HP (I) — X (l, I)
- 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)
- NAP (I) = TT (I) + X (l, I) HP (I) TT (I) = HP (I) — X (l, I)
- NAP (I) ≥ .4 * BP (I) THEN NAP (I) = .4 * BP (I) END IF
- NAP (I) = TT (I) + X (l, I) HP (I) TT (I) = HP (I)-X (1,I)
- NAP (I) ≥ .4 * BP (I) THEN NAP® = .4 * BP (I) END IF
- 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)
- CALL SPLAIN (X (), Y (), TT (I), I, TT (I) + X (l, I), Y, KOL. TOCH ()):
- BB (I) = Y * 2: NACH. B (I) = BB (I)
- BB (I) ≤ 0 THEN BB (I) = .0001
- BATAA (I) = GG (I) / (TT (I) * BB (I))1. PACHETKnMHK?
- CALL KPMKD (ka, TT (), BB (), BATAA (), NAP (I), BP (I), I, M33(), M330(), L33(), L330(), SIG) 55:1 = 1 + 1: WEND: D = IN * GAMMA
- DM33(0) = (M33(0) M33(l)) * 2: DM330(0) = (M330(0) — M330(l)) * 2
- DM33(H%) = (M33(H% 1) — M33(H%)) * 2: DM330(H%) = (M330(H% - 1) — M330(H%)) * 2
- DL33(0) = (L33(0) L33(l)) * 2: DL330(0) = (L330(0) — L330(l)) * 2
- DL33(H%) = (L33(H% 1) — L33(H%)) * 2: DL330(H%) — (L330(H% - 1) — L330(H%)) * 21. FOR I = 1 TO H% -1
- 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
- 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)
- 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)
- 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
- КАРА = (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))
- KAPA1 = ksiil * SIN (alfa#) + ksii5 * COS (alfa#) л 2
- КАРА > 1 THEN КАРА = 1: IF KAPA1 > 1 THEN KAPA1 = 166: El = КАРА * GAMMA * BB (I): E2 = L330(I)
- НЗ = KAPAl * SIG л 2 * M330(I): E4 = KAPA1 * v / 2 * DL330(I) / (DL * 2) E5 = v / 2 * DM330(I) / (DL * 2)
- F11 = ZA * (El E3 — E4): F22 = ZA * KAPA1 * SIG * (E2 — E5)
- Ml 1 = Fll * X: M22 = F22 * X
- CKX = C0S (KA1 * X): SKX = SIN (KA1 * X)
- 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
- Q.N (I) = -.5 * (ZXFVFZ * (DZM3(I) * ZXFVF DZM30(I) * ZVD) + ZXFZ * (DZL3(I) * ZXF -DZL30(I) * ZVD))
- ZXFVFZ < 0 THEN Q. N (I) = -.5 * (ZXFZ * (DZL3(I) * ZXF DZL30(I) * ZVD))
- I = 0 OR I = H% THEN Q. N (I) = Q. N (I) / 2: QN = QN + Q. N (I): QNX = QNX + Q. N (I) * X
- NEXT I: FF = F1* СOS (SIGk * IE) F2 * SIN (SIGk * IE) + QN
- MM = Ml * COS (SIGk * IE) M2 * SIN (SIGk * IE) + QNX
- ERASE МЗЗ, M330, L33, L330, DM33, DM330, DL33, DL330
- ERASE KSI11#, KSI21#, KSI51#, KSI12#, KSI22#, KSI52#
- ERASE sil 1#, si21#, si51#, sil2#, si22#, si52#, KS11#, KS21#, KS51#
- ERASE TT, BB, GG, BATAA, NAP1. END SUB
- SUB run 15 (Iv, SIG, ка)» РАСЧЕТ КОЭФФИЦИЕНТОВ УРАВНЕНИЙ КАЧКИ
- SHARED A11, B11, CI 1, A12, B12, C12, A21S, B21S, C21S, A22S, B22S, C22S, FF, MM, QN, 1. QNX, F1, F2,M1,M21. ПО ВЫВЕДЕННЫМ ФОРМУЛАМ
- 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
- 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)
- DIM QQ (H% + 1), QQI (H% + 1), MM (H% + 1), GS (H% + 1)
- DIM sil 1#(12), si21#(12), si51#(12), sil2#(12), si22#(12), si52#(12)
- RO = GAMMA / 9.81: alfa# = (ALF / 180) * 3.141 592 654#: DL = LS / H%
- SIGk = SIG ((SIG A 2 / 9.81) * v * COS (alfa#)): KA1 = -ka * COS (alfa#)
- KA2 = ka * SIN (alfa#): ZA = .11 * lv A .715 'полувысота волны
- ZA = 4.67:1 = 0: WHILE I ≤ H%
- CALL KPMKD (ka, T (), B (), BATA (), 0, 0,1, M33(), M330(), L33(), L330(), SIG)
- 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
- DM33(I) = 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
- DM33(0) = (M33(0) M33(l)) * 2: DM330(0) = (M330(0) — M330(l)) * 2
- DM33(H%) = (M33(H% -1) M33(H%)) * 2: DM330(H%) = (M330(H% - 1) — M330(H%)) * 2
- DL33(0) = (L33(0) L33(l)) * 2: DL330(0) = (L330(0) — L330(l)) * 2
- DL33(H%) = (L33(H% -1) L33(H%)) * 2: DL330(H%) = (L330(H% -1) — L330(H%)) * 2
- CALL raspred. nagruzki: FOR I = 0 TO H%: XA (I) = DL * ((H% / 2) -1): IF I = 0 OR I = H%
- THEN: A11(1) = (2 * PI (I) / (9.81 * DL) + M33(I)) / 2: B11(1) = (L33(I) v / 2 * DM33(I) / (2 *
- DL)) / 2: CI 1(1) = (GAMMA * B (I) v / 2 * DL33(I) / (2 * DL)) / 2: A12(I) = All (I) * XA (I)
- КАРА = (М / (Т (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))
- KAPA1 = ksiil * SIN (alfa#) + ksii5 * COS (alfa#) A 2
- 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)
- F11 = ZA * (El E3 — E4): F22 = ZA * KAPA1 * SIG * (E2 — E5)
- Mil — F11 * X: M22 = F22 * X
- CKX = C0S (KA1 * X): SKX = SIN (KA1 * X)
- FIM = ((FI1 A 2 + FI2 A 2) A .5) / alfO ERASE A, letb, C, letg, letx
- J = 1: WHILE J ≤ 2: 'определение производных FI, ZETA
- РЕШЕНИЕ СИСТЕМЫ FID, ZETAD
- FID = FI1 * COS (OMEGA#) + FI2 * SIN (OMEGA#)
- ZETAD = ZETA1 * COS (OMEGA#) + ZETA2 * SIN (OMEGA#)
- FIM = ((FI1 A 2 + FI2 л 2) л .5) / alffl1. FOR I = 0 TO H%1. I = 0 OR I = H% THEN
- 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
- QQ (I) = -A11(1) * ZETAD2 В11 (I) * ZETAD 1 — CI 1(1) * ZETAD — A12(I) * FID2 — B12(I) *
- FID1 С 12(1) * FID + FS (I) * COS (OMEGA#) — FSS (I) * SIN (OMEGA#)
- QQI (I) = QQI (I 1) + QQ (I -1) + QQ (I): MM® = MM (I -1) + QQI (I -1) + QQI (I): END IF
- PIR.SILA (I, J) = QQI (I) * DL/2: IZG. MOMENT (I, J) = MM (I) * (DL / 2) л 2: NEXT I: J = J + 1
- OMEGA# = 3,141 592 654# / 2: WEND: MAX. S = 0: MAX. M = 0: FOR I = 0 TO H%
- P.S (I) = SQR (PIR.SILA (I, 1) л 2 + PIR. SILA (I, 2) л 2)
- M (I) = SQR (IZG.MOMENT (I, 1) л 2 + IZG. MOMENT (I, 2) л 2)
- MAX.S < ABS (P.S (I)) THEN MAX. S = ABS (P.S (I))
- MAX.M < ABS (I.M (I)) THEN MAX. M = ABS (I.M (I)): NEXT I
- ERASE МЗЗ, M330, L33, L330, DM33, DM330, DL33, DL330
- ERASE sill#, si21#, si51#, sil2#, si22#, si52#1. ERASE QQI, MM1. END SUB
- ПОДПРОГРАММА ИСХОДНЫХ ДАННЫХ
- OPEN «rollingV + NEMA$ FOR OUTPUT AS #1
- PRINT #1, TAB (28)-» ИСХОДНЫЕ ДАННЫЕ ДЛЯ «
- PRINT #1, TAB (24) — «РАСЧЕТА ПРОДОЛЬНОЙ КАЧКИ СУДНА"1. PRINT #1,»»
- PRINT #1, «Количество расчетных шпаций Н = «- Н%
- PRINT #1, «Водоизмещение судна, т. D = «- D
- PRINT #1, «Длина судна по ГВЛ. м. L = «- LS
- PRINT #1, «Осадка судна на носовом перпен-ре м Тн= «- DN
- PRINT #1, «Осадка судна на кормовом перпен-ре м Тк= «- DK
- PRINT #1, «Скорость хода судна, м/с. V = «- v
- PRINT #1, «Угол встречи с волной, градусы ALFA = «- ALF
- PRINT #1, «Абсцисса центра тяжести судна, м. Xg= «- Xg
- PRINT #1, «Момент инерции относит. ЦТ. м. А4 Iyy = «- iyy
- PRINT #1, «Момент инерции площади ГВЛ. м. л4 ffi = «- Ifi
- PRINT #1, «Длина волны, м. LV =" — lv
- PRINT # 1, «Удельный вес воды, т/куб.м GAMMA = «- GAMMA1. PRINT #1,»»
- SUB star (NAM, R) «ОПРЕДЕЛЯЮЩАЯ ПОДПРОГРАММА
- 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
- ON KEY (9) GOSUB 1090: 'KEY (9) ON: ON KEY (IO) GOSUB 1000: KEY (10) ON
- NAM = 1: R = 1: Ql = 1: WHILE Q1 = 1: CALL submenu (NAM, R) файл
- SELECT CASE NAM CASE 1: SELECT CASER CASE 1: CALL NEWFILE CASE 2: CALL openfile CASE 3:1. CALL SAVE CASE 4:
- CALL messagel (n$, «СОХРАНИТЬ КАК.»): CALL SAVE END SELECT 'редактирование CASE 2: SELECT CASE R
- CASE 1: CALL REDTOTAL CASE 2: CALL RED (X (), Y (), KOL. TOCH ()) END SELECT 'расчет
- CASE 1: CALL prnscreen (47) END SELECT CASE 7: SELECT CASER
- CASE 2: CALL qlstn (22, 30, «ВЫХОД ?», F) IF F = 1 THEN
- CALL qlstn (22, 30, «СОХРАНИТЬ ДАННЫЕ ?», Ql) IFQ1 = 1 THEN CALL SAVE END IF
- Ql = 0: END IF: END SELECT: END SELECT: WEND: KILL «*'.MJG»: END END SUB
- Calculation of Ship Motion and Bending Moment and Comparison ofthe Programs
- Oleg E. Surov, Far -Eastern State Technical University. II-Woong Kim, University of Ulsan1998 Ulsan
- GENERAL INFORMATION GENERAL PARTICULARS
- THIS BOOKLET IS PREPARED TO ENABLE THE TRIM & STABILITY AND LONGITUDINAL STRENGTH CALCULATION TO BE EASILY PERFORMED.
- 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
- LENGTH BETWEEN PERPENDICULARS3. BREADTH (MOULDED)4. J3EPTH (MOULDED)5. DESIGN DRAUGHT (MOULDED)
- M/TUKONG NAVIGATOR» 315.00 M 57.20 M 30.40 M 20.45 M 20.45 M
- 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 5 154 651 I 83 050 166 919.1 / 156 995.9
- KOREAN REGISTER OF SHIPPING +KRS1, +KRM1, UMA-2
- AMERICAN BEREAU OF SHIPPING
- A1(E), «OIL CARRIER», +AMS, +ACCUii. tank plan1. PBOft-E