Математическое моделирование движения траловой системы тема диссертации и автореферата по ВАК РФ 05.18.17, кандидат технических наук Зинченко, Валерий Павлович

• Зинченко, Валерий Павлович
• кандидат технических наук
• 2000

Оглавление диссертации кандидат технических наук Зинченко, Валерий Павлович

Введение.

1. Методы расчета движения траловой системы.

1.1. Общие замечания.

1.2. Метод расчета статики траловой системы В.П.Карпенко.

1.3. Метод расчета статики траловой системы А.Л.Обвинцева.

1.4. Расчет статики и динамики буксируемых систем по НХекгке.

1.5. Двухстержневая модель Б.А.Альтшуля. Маятниковые модели динамики траловой системы.

1.6. Выводы.

2. Постановка задачи расчета движения траловой системы.

2.1. Принципы схематизации траловой системы.

2.1.1 .Моделирование трала.

2.1.2. Моделирование движения распорных досок.

2.1.3. Моделирование действия устройства вертикального раскрытия.

2.1.4. Учет веса и гидродинамического сопротивления ваеров и кабелей.

2.2. Математическая постановка задачи расчета движения траловой системы.

2.3. Численные алгоритмы расчета движения траловой системы.

3. Исследование нестационарного движения разноглубинного трала.

3.1. Численное исследование точности разработанной модели нестационарного движения траловой системы.

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

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

3.4. О возможности расчета коэффициентов модели трала по эмпирическим формулам.

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

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

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

Большое значение для дальнейшего совершенствования тралового лова имеет создание математической модели нестационарного движения траловой системы, обеспечивающей возможность расчета параметров движения с высокой точностью и с учетом настройки системы. Эта модель может использоваться для подготовки судоводителей и мастеров добычи промысловых судов, для чего необходимо создание рыбопромысловых тренажеров, действующих на ее основе. Опыт использования первого отечественного рыбопромыслового тренажера описан в работе [17]. Без подобных моделей невозможно также создание системы автоматического управления скоростью судна, курсом судна и длиной ваеров, обеспечивающей реализацию заданного закона движения трала в водном пространстве.

Отсутствие до настоящего времени математических моделей нестационарного движения траловой системы при различной ее настройке имеет две причины. Первая причина состоит в сложности физических процессов, происходящих при движении трала: гидродинамика нитей, образующих трал, представляет собой самостоятельную научную дисциплину [5]; формоизменяемость и взаимная зависимость геометрической формы трала и приложенных к нему сил известны как главные трудности теории траловой системы [37,21]; движение трала определяется гидродинамическими параметрами его оснастки, получение которых теоретическим путем в настоящее время невозможно. Второй причиной является сложность конструкции трала. Большое количество нитей, составляющих трал и его оснастку, приводит к пропорциональному росту числа уравнений, описывающих движение траловой системы. Даже простые схемы траловой системы, включающие небольшое количество нитей, вызывают существенные математические трудности при расчете статики, не говоря уже о динамике [6,7,16,26]. Именно поэтому расчет статики траловой системы обычно не рассматривается как математическая задача. Такая постановка вопроса препятствует реализации значительно выросших за последние десять лет возможностей современной вычислительной техники в исследованиях траловой системы. Поэтому внедрение в теорию траловой системы математических численных методов является особенно актуальным.

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

Динамика траловой системы в настоящее время рассматривается, в основном, в рамках двухстержневой маятниковой схемы, разработанной Б.А.Альтшулем [2-3]. С помощью этой схемы возможно решение различных задач, в частности, задачи расчета управляющих воздействий для реализации заданного закона движения трала в плоском случае. Однако, для решения задачи расчета движения тралового комплекса при различной его настройке двухстержневая схема не подходит. Для использования двухстержневой схемы необходимо, чтобы изменение углов ориентации распорных досок было малым, что не всегда имеет место. В противном случае изменение гидродинамических сил, приложенных к доскам, не может, вообще говоря, быть учтено в рамках двухстержневой схемы даже в случае плоского движения, что будет показано в гл.З.

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

В диссертационной работе удалось реализовать новую схему траловой системы, изображенную на рис.3. В модели трала представлены не только подборы, но и сборочные, а голые концы, имеющие значительную длину, моделируются отдельно от подбор. Это соответствует особенностям конструкции современных разноглубинных тралов (см.[13], стр.34). На концах подбор и в их середине, а также в точках соединения сборочных смежных пластей трала, приложены силы, моделирующие непрерывное распределение силы натяжения канатно-сетной оболочки трала. Концы подбор соединяются с распорными досками отдельными связями, длина которых может быть различной. Таким образом, впервые появляется возможность моделировать влияние изменения длины голых концов подбор и кабелей на геометрию устья трала.

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

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

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

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

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

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

В диссертационной работе сформулировано количественное условие наступления нестационарного режима движения. Для этого был проделан ряд численных экспериментов. Скорость буксировки изменялась по периодическому закону, период и размах которого выбирались таким образом, чтобы в процессе движения достигалось заданное ускорение буксировки при допустимых значениях скорости. Зависимости координат траловой системы и их производных при различных значениях указанных факторов фиксировались на графиках. Полученные графические зависимости сравнивались между собой. Оказалось, что если ускорение буксировки не превышает 0.03 м/с2, то относительное изменение амплитуды колебаний всех координат системы и их производных до второго порядка включительно при изменении коэффициента демпфирования досок в 2 раза не превышает 25%. Аналогичные утверждения имеют место при изменении в 2 раза коэффициента упругости стержней, а также при изменении массы распорных досок и массы материальных точек модели трала. На основании в работе сделан вывод, что величина ускорения 0.03 м/с2 может использоваться в качестве достаточного условия наступления режима движения, при котором использование разработанной модели обеспечивает приемлемую точность.

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

Научная новизна работы состоит в следующем.

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

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

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

Заключение диссертации по теме «Промышленное рыболовство», Зинченко, Валерий Павлович

Заключение.

В заключение, сформулируем основные результаты выполненных исследований.

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

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

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

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

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

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

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

8. Разработанная математическая модель траловой системы внедрена в рыбопромысловый тренажер, спроектированный в КБ "Вектор" (г.Таганрог) и поставленный в ведущие учебные заведения рыбной отрасли.

Список литературы диссертационного исследования кандидат технических наук Зинченко, Валерий Павлович, 2000 год

1. Альтшуль Б.А., Фридман А.Л. Динамика траловой системы. М.: Агропромиздат, 1990,- 240с.

2. Альтшуль Б.А. Основы теории движения траловой системы // Дисс. на соиск. уч. ст. д.т.н. Калининград: КГТУ, 1993. - 294с.

3. Арушанян О.Б., Залеткин С.Ф. Численное решение обыкновенных дифференциальных уравнений на Фортране. М.: Изд. МГУ, 1990. - 336с.

4. Белов В.А., Короткое В.К., Саврасов В.К., Шимянский С.Л. Буксируемые орудия лова. М.: Агропромиздат, 1987.- 200с.

5. Габрюк В.И. Параметры разноглубинных тралов. М.: Агропроиздат, 1988. - 212с.

6. Габрюк В.И. Компьютерные технологии в промышленном рыболовстве. М.: Колос, 1995. -544с.

7. Данилов Ю.А., Фридман А.Л. Гидродинамическое сопротивление сетей, ниток и тросов (обзор). М.: ВНИРО, 1968. - 64с.

8. Дверник A.B. Совершенствование методики расчета сопротивления рыболовного трала // Дисс. на соиск. уч.ст. к.т.н.- Калининград: КТИРПХ, 1971. 255с.

9. Долгов А.Н., Зинченко В.П., Козлов A.A. Рыбопромысловый тренажер тралового и кошелькового лова //аналит. и реф. инф. ВНИЭРХ, сер. «Промысловая радиоэлектронная аппаратура». М.: ВНИЭРХ, вып.7-8,1998г. - с. 21-24.

10. Доска траловая крыловидная пр.2490. Сборочные чертежи. Калининград: НПО по технике промрыболовства.

11. Изнанкин Ю.А., Лисовой А.П. Учет внутренних свойств рыболовных сетей в устройстве орудий лова. Калининград: КТИРПХ, 1989. - 96с.

12. Зинченко В.П. Численный метод расчета движения траловой системы // рукопись деп. во ВНИЭРХ, реферат опубл. в указателе деп. работ ВИНИТИ N12, 1998г.

13. Зинченко В.П., Долгов А.Н., Козлов A.A. Математическое моделирование нестационарного движения разноглубинного трала //инф. пакет ВНИЭРХ, сер. «Промышленное рыболовство». М.: ВНИЭРХ, вып.1, 2000г. - с. 10-13.

14. Карпенко В.П. Основы теории и расчета устройств раскрытия траловых систем // Дисс. на соиск. уч. ст. д.т.н. Калининград: КГТУ, 1996. - 363с.

15. Когон М.Г., Петелин В.П., Румянцев А.И. Обучение на тренажере разноглубинного лова // Рыбное хозяйство, 1977, N5, с.65-67.

16. Козлов A.A. Методы расчета гидродинамических сил для моделирования движения трала в рыбопромысловом тренажере // Дисс. на соиск. уч.ст. к.т.н. -Калининград: КГТУ, 2000г.

17. Комплексные сравнительные испытания разноглубинных канатных тралов в районе ЮВА // Отчет ВНИРО. М.: ВНИРО, 1985. - 103с.

18. Коротков В.К. Расчет мотенной части с заданными рабочими параметрами. Учебное пособие. Калининград: ВИПК Минрыбхоза СССР, 1990. - 43с.

19. Кулагин В.Д. Статика сетных оболочек. Калининград: КГТУ, 1995. - с. 104.

20. Морские орудия лова Северного бассейна. Рекомендации промысловикам. -Мурманск: ЦПКТБ ВРПО "Севрыба", 1986. 215с.

21. Недоступ A.A. Исследование гидродинамического коэффициента сопротивления тралов // Дисс. на соиск. уч. ст. к.т.н. Калининград: КГТУ, 1999г. - 194с.

22. Обвинцев А.Л. Взаимосвязь между вертикальным и горизонтальным раскрытием трала // Рыбное хозяйство, 1975, N12, с. 39-42.

23. Обвинцев А.Л. Расчет углов ориентации досок разноглубинного трала // Рыбное хозяйство, 1978, N1, с. 50-52.

24. Обвинцев А.Л. Комплексный метод расчета технических характеристик разноглубинных тралов с помощью ЭВМ // Дисс. на соиск. уч. ст. к.т.н. М.: ВНИРО, 1986.- 175с.

25. Разработка методики аналитического расчета конструкций трала // Заключительный отчет о НИР по теме 3134 НПО по технике промрыболовства. Руководитель темы В.А.Белов. Калининград: НПО по технике промрыболовства, 1985. -131с.

26. Рейс БМРТ "Гангут" в районах САХ и ЦВА в декабре 1981- мае 1982 гг. // Рейсовый отчет 2764 НПО по технике промрыболовства. Калининград: НПО по технике промрыболовства, 1981.- 131 с.

27. Розенштейн М.М. Механика орудий промышленного рыболовства. Часть 3. Механика рыболовной сети. Калининград: КГТУ, 1997. - 47с.

28. Розенштейн М.М., Альтшуль Б. А. Механика орудий промышленного рыболовства. 4.1. Механика гибкой нити. Калининград: КГТУ, 1995. - 94с.

29. Розенштейн М.М., Недоступ A.A. Методы расчета сил сопротивления канатно-сетной части трала // обз. инф. ВНИЭРХ, сер. "Промышленное рыболовство". М.: ВНИЭРХ, 1998, вып.2. - с. 1-24.

30. Соловьев A.A. Переходные режимы движения судна при буксировке разноглубинного трала // дисс. на соиск. уч.ст. к.т.н. Ленинград: ЛВИМУ, 1989. - 159с.

31. Фонарев А.Л. Гидромеханика // Методические разработки по теме "Сопротивление сетной части орудий лова". Калининград: КГТУ, 1995. - 93с.

32. Фридман А.Л. Теория и проектирование орудий промышленного рыболовства. М.: Легкая и пищевая промышленность, 1981. - 328с.

33. Leitzke Н. Theorie der Zugsysteme ein neues Berechnungsverfahren zur Loesung fischereitechnischer Aufgabenstellungen // Fischerei Forschung. - Rostock, В19, N3, 1981, S65-68.

34. Leitzke H. Berechnung von im Wasser beschleunigt geschleppten, komplizierten diskretisierten Systemen in endlich grossen Zeitschritten // Fischerei Forschung. Rostock, B25, N1, 1987, S46-52.

35. Moderhak W. Opor hydrodynamiczny czesci sieciowej wloka pelagicznego// Stud. Mater. Mor. Inst. Ryb. Gdynia, Ser. B, N64, 1993. - 119 s.

36. Положение ц.м. указано для доски с дугой без грузов. Масса доски с дугой без грузов тс1 =1.346 т

37. Масса одного груза тг = 0.06 тprogram model;uses Forms,

38. Opengl in 'd:\program files\borland\delphi4\source\rtl\win\opengl.pas'; {$R *.RES} begin1. Application.Initialize;

39. Application.CreateForm(Tmodelf, modelf); Application.Run; end.unit umodel;interfaceuses

40. Windows, Messages, SysUtils, Classes, Graphics, Controls, Forms, Dialogs, ExtCtrls,StdCtrls, Buttons,ComCtrls,math,umath,umcon,ugcon,uvcon,uvtypes, uplotter, uvessel, upagecontolf,ueyef; type

41. ГФ коорд. точек подбор трала и их производныеrp0ip,ix.1.:=rp0[ip,ix][i]+(vp[ip,ix][i]+vfu[i])*dt;rpip,ix. 1.:=rp0[ip,ix][i]-rs [i] ;

42. ГФ скорость точек подбор отн. верхнего течения vpip,ix.1.:=vp[ip,ix][i]+vpt[ip,ix][i]*dt; end;for ip:=0 to 7 do begin ip1:=nextp(ip); rpOip,1.:=rp0[ipl,-1]; rp[ip,1]:=rp[ipl,-1]; vp[ip,1]:=vp[ipl,-1]; end;for i:=l to 3 do begin

43. ГФ координаты центра трала отн. земли rcO1.:=rc0i.+(vc[i]+vfu[i])*dt; rc[i]:=rc0[i]-rs[i];

44. ГФ скорость ц.т. отн. верхнего течения vc1.:=vci.+vct[i]*dt; end; end;end; // if trlnev=l• modeltime:=modeltime+dt;***конец расчета нового состояния***end; end.unit udoors; interfaceuses math,umath,umcon,uvcon,utrawl; type

45. Матрица перехода в/из систему, повернутую вокруг трех осей // из исходной ГФ системы.

46. Директор конструкторского морской электроники к.т.н.1. Долгов А.Н.

47. Внедрения результатов диссертационной работы Зинченко В.П.

48. Вид внедренных результатов: математическая модель нестационарного движения траловой системы.

49. Характеристика масштаба внедрения: два вида промысловых тренажеров, установленных в МГТУ и БГАРФ.

50. Новизна результатов: создан метод математического моделирования нестационарного движения траловой системы с учетом ее настройки.

51. Социальный и научно-технический эффект: внедрение результатов позволило создать тренажер с принципиально новыми потребительскими свойствами (подводная визуализация).1. Сыротюк А.П.1. Кудрявцев H.H.s I —

52. Десятерик М.Н. , "А/ ^fw-1