МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ РАСПОРНОЙ ТРАЛОВОЙ ДОСКИ ДЛЯ ЗАДАЧ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ ДИНАМИКИ РАЗНОГЛУБИННЫХ ТРАЛОВЫХ СИСТЕМ
Авторы:
1.
Калининградский государственный технический университет
(кандидат технических наук, доцент; заведующий кафедрой промышленного рыболовства)
2.
Калининградский государственный технический университет
(кандидат технических наук; ведущий научный сотрудник)
Тип:
Статья
Страницы:
с 58
по 64
Статус:
Опубликован
Получено:
05.09.2017
Одобрено:
25.09.2017
Опубликовано:
25.09.2017
Классификаторы:
УДК 639.2.081.117.21,519.6
ГРНТИ 34.39 Физиология человека и животных
ГРНТИ 62.13 Биотехнологические процессы и аппараты
ГРНТИ 69.01 Общие вопросы рыбного хозяйства
ГРНТИ 69.25 Аквакультура. Рыбоводство
ГРНТИ 69.31 Промышленное рыболовство
ГРНТИ 69.51 Технология переработки сырья водного происхождения
ГРНТИ 87.19 Загрязнение и охрана вод суши, морей и океанов
ГРНТИ 34.39 Физиология человека и животных
ГРНТИ 62.13 Биотехнологические процессы и аппараты
ГРНТИ 69.01 Общие вопросы рыбного хозяйства
ГРНТИ 69.25 Аквакультура. Рыбоводство
ГРНТИ 69.31 Промышленное рыболовство
ГРНТИ 69.51 Технология переработки сырья водного происхождения
ГРНТИ 87.19 Загрязнение и охрана вод суши, морей и океанов
Язык материала:
русский
Ключевые слова:
разноглубинная траловая система, траловая доска, математическая модель, метод расщепления, точечные массы, графический процессор, компьютерная программа
Аннотация (русский):
Предлагается математическое описание распорной траловой доски разноглубинной траловой системы для численного моделирования методом точечных масс с использованием параллельной архитектуры графического процессора, основанное на представлении зависимости силовых характеристик доски от ее пространственных и кинематических характеристик в табличном виде. Таблицы построены с использованием метода расщепления. Для моделирования динамики разноглубинной траловой системы методом точечных масс необходимо знать все массы, сосредоточенные в узлах, зависимости расстояний между узлами от внутренних сил и значения внешних сил, приложенных к узлам (гидростатических, гидродинамических, сил тяги траулера). Для канатно-сетной части трала связью считается участок каната, троса, нитки или веревки, узлом - место их соединения. Для внедрения траловой доски в общую имитационную модель ее необходимо представить в виде узлов и связей. Для моделирования влияния траловой доски на динамику всей траловой системы по полю скоростей с использованием систем уравнений Навье - Стокса численным методом расщепления вычисляется поле давлений. Поскольку расчет внешних сил из-за большой трудоемкости метода расщепления на каждой итерации алгоритма точечных масс на персональном компьютере в реальном времени не представляется возможным, предлагается более быстрый способ, основанный на композиции методов расщепления для траловых досок и точечных масс для канатно-сетной части траловой системы с использованием табличных преобразований и линейной интерполяции промежуточных значений. Для каждого типа траловой доски производится моделирование методом расщепления для всех возможных комбинаций линейных и угловых скоростей с заданным шагом дискретизации, который определяется путем сравнения с имеющимися данными на основании необходимой точности. По результатам каждого шага вычисляется поле давлений, а затем и внешние силы, приложенные к узлам доски. Полученная шестимерная таблица ставит в соответствие скоростям внешние силы. Для описания конструкции траловой доски предлагается кусочно-линейная аппроксимация профиля всех ее элементов (крыло, предкрылки) в поперечной плоскости сечения с заданной степенью дискретизации, определяемой экспериментально с учетом заданной погрешности вычислений. Модель приемлема для описания большинства применяемых траловых досок.
Предлагается математическое описание распорной траловой доски разноглубинной траловой системы для численного моделирования методом точечных масс с использованием параллельной архитектуры графического процессора, основанное на представлении зависимости силовых характеристик доски от ее пространственных и кинематических характеристик в табличном виде. Таблицы построены с использованием метода расщепления. Для моделирования динамики разноглубинной траловой системы методом точечных масс необходимо знать все массы, сосредоточенные в узлах, зависимости расстояний между узлами от внутренних сил и значения внешних сил, приложенных к узлам (гидростатических, гидродинамических, сил тяги траулера). Для канатно-сетной части трала связью считается участок каната, троса, нитки или веревки, узлом - место их соединения. Для внедрения траловой доски в общую имитационную модель ее необходимо представить в виде узлов и связей. Для моделирования влияния траловой доски на динамику всей траловой системы по полю скоростей с использованием систем уравнений Навье - Стокса численным методом расщепления вычисляется поле давлений. Поскольку расчет внешних сил из-за большой трудоемкости метода расщепления на каждой итерации алгоритма точечных масс на персональном компьютере в реальном времени не представляется возможным, предлагается более быстрый способ, основанный на композиции методов расщепления для траловых досок и точечных масс для канатно-сетной части траловой системы с использованием табличных преобразований и линейной интерполяции промежуточных значений. Для каждого типа траловой доски производится моделирование методом расщепления для всех возможных комбинаций линейных и угловых скоростей с заданным шагом дискретизации, который определяется путем сравнения с имеющимися данными на основании необходимой точности. По результатам каждого шага вычисляется поле давлений, а затем и внешние силы, приложенные к узлам доски. Полученная шестимерная таблица ставит в соответствие скоростям внешние силы. Для описания конструкции траловой доски предлагается кусочно-линейная аппроксимация профиля всех ее элементов (крыло, предкрылки) в поперечной плоскости сечения с заданной степенью дискретизации, определяемой экспериментально с учетом заданной погрешности вычислений. Модель приемлема для описания большинства применяемых траловых досок.
Ключевые слова:
разноглубинная траловая система, траловая доска, математическая модель, метод расщепления, точечные массы, графический процессор, компьютерная программа
разноглубинная траловая система, траловая доска, математическая модель, метод расщепления, точечные массы, графический процессор, компьютерная программа
Введение Постановка задачи. Для моделирования динамики разноглубинной траловой системы методом точечных масс [1, 2] необходимо знать все массы, сосредоточенные в узлах, (точечные массы) с учетом присоединенной массы, зависимости расстояний между всеми связями от внутренних сил (сил натяжения) и значения внешних сил, приложенных ко всем узлам (гидростатических, гидродинамических, сил тяги траулера). Для канатно-веревочного изделия (КВИ), которым является трал, связью считается участок каната, троса, нитки или веревки, а узлом - место их соединения. Для моделирования прогиба каната, троса или веревки в модель точечных масс вводятся дополнительные узлы, расположенные на канате (тросе, веревке …) и делящие его (ее) на части. В траловую систему входят не только КВИ, но и такие элементы конструкции, как распорные траловые доски, которые необходимо внедрить в общую математическую модель. Для этого траловую доску необходимо представить в виде узлов и связей, т. е. создать математическую модель точечных масс для траловой доски. Математическое описание Рассмотрим случай, когда расстояния между точками крепления ваеров и лапок к траловой доске остаются фиксированными в течение всего процесса моделирования (соединение крыльев и предкрылок жесткое). При переходе к математической модели заменим точки крепления ваера (точка 1 на рис. 1) и лапок к траловой доске (точки 2 и 3 на рис. 1) узлами. Введем три инцидентные каждой паре узлов связи - a, b и c (рис. 1). Длины данных связей определены, фиксированы в течение всего процесса моделирования и не зависят от сил сжатия/растяжения (величина, обратная модулю упругости связи при сжатии и растяжении равна нулю). Рис. 1. Распорная траловая доска с обозначением узлов и связей модели точечных масс Для определения массы, сосредоточенной в каждой из трех точек, необходимо знать конструкцию траловой доски и плотность ее материалов. Для определения внешних сил, действующих на узлы траловой доски, необходимо знать вес части траловой доски в воде, сосредоточенный в узле, и поле давлений воды на поверхностях траловой доски. Современные траловые доски имеют достаточно сложную конструкцию, поэтому для точного моделирования их влияния на динамику всей траловой системы поле давлений необходимо вычислять с использованием уравнений Навье - Стокса [3]. Одним из методов численного решения данных уравнений является метод расщепления [4]. Метод достаточно трудоемок с точки зрения вычислений на ЭВМ. Одним из направлений оптимизации вычислений является параллельная (одновременная) обработка множества ячеек трехмерной вычислительной сетки с использованием вычислительных ресурсов графического процессора (графических процессоров) [5]. Нами разрабатывается программный код алгоритма метода расщеплений для параллельных вычислений на персональном компьютере с графическим процессором в операционной системе Windows с использованием программного интерфейса DirectCompute [6]. Из-за большой трудоемкости численного решения уравнений Навье - Стокса методом расщепления, даже с учетом вычислений на графическом процессоре, его применение для расчета внешних сил на каждой итерации алгоритма точечных масс на персональном компьютере в реальном времени не представляется возможным. В связи с этим нами предлагается более быстрый способ, основанный на композиции двух методов: расщепления для траловых досок и точечных масс для канатно-веревочной части траловой системы. Взаимодействие между методами (математическими моделями) предлагается осуществлять через табличные преобразования с линейной интерполяцией промежуточных значений. Опишем подробнее принцип взаимодействия. Траловая доска представляется как абсолютно твердое тело с шестью степенями свободы. Так как разноглубинный трал работает далеко от поверхности дна водоема, то его влиянием на поля скоростей и давлений водной массы можно пренебречь. Можно пренебречь также искажением этих полей, вызванным другими частями траловой системы из-за достаточного их удаления. С учетом вышеупомянутых упрощений для определения поля давлений на поверхностях траловой доски достаточно знать скорость (в векторной форме) перемещения всех точек поверхностей относительно воды в трехмерном пространстве. В свою очередь, для определения скорости любой из точек на поверхностях доски достаточно знать скорость поступательного движения центра масс и три угловые скорости вращения доски вокруг осей декартовой системы координат с началом в центре масс. Для каждого типа траловой доски производится моделирование методом расщепления для всех возможных комбинаций поступательных и угловых скоростей с заданным шагом дискретизации (определяется путем сравнения с имеющимися данными на основании необходимой точности). По результатам каждого шага вычисляется поле давлений, а затем и внешние силы, приложенные к трем узлам доски (в модели точечных масс). В результате мы получаем шестимерную таблицу, ставящую в соответствие скоростям внешние силы. При переходе к математической модели точечных масс с использованием геометрических преобразований достаточно преобразовать мгновенные скорости трех узлов (9 вещественных значений) в поступательную и угловые скорости (6 вещественных значений) и по таблице, с использование линейной интерполяции промежуточных значений, найти значения внешних сил для трех узлов траловой доски. Для математического описания конструкции траловой доски нами предлагается кусочно-линейная аппроксимация [7] профиля всех деталей доски (крыло, предкрылки) в поперечной плоскости сечения с заданной степенью дискретизации, определяемой экспериментально с учетом заданной погрешности вычислений (например, 32). При этом условимся, что профили любого поперечного сечения идентичны. Такая математическая модель приемлема для описания большинства применяемых траловых досок. Поперечный профиль предлагается определять методом кусочно-квадратичной интерполяции [8] по заданным базовым точкам (рис. 2). Рис. 2. Профиль детали в поперечном сечении с базовыми точками Для каждой детали траловой доски необходимо задать профиль, ее геометрические размеры (размах, хорда), положение относительно других деталей и другие характеристики (рис. 3). Необходимо также задать положение точек крепления ваера и лапок. Рис. 3. Характеристики траловой доски, деталей оснастки и точек крепления ваера и лопаток По предложенному алгоритму построения математической модели авторами была разработана компьютерная программа «Конструктор крыловидных траловых досок». Пример трехмерного отображения математической модели траловой доски с крылом и четырьмя предкрылками, созданной в программе, показан на рис. 4 и 5. Рис. 4. Вид сверху поверхностей траловой доски и координатной сетки Рис. 5. Вид снизу поверхностей траловой доски Данная компьютерная программа позволяет создавать трехмерные модели прямокрылых траловых досок с произвольным количеством и положением предкрылок с указанием характеристик крыла и предкрылок, а также точек крепления лопаток и ваера, которые в дальнейшем можно использовать для их моделирования в составе траловой системы. Заключение Для численного моделирования траловых систем методом точечных масс предложено математическое описание распорной траловой доски разноглубинной траловой системы. Данный метод позволяет оптимизировать имитационную модель траловой системы для вычислений с использованием параллельной архитектуры графического процессора. Математическая модель распорной траловой доски, разработанная по предложенному математическому описанию, позволяет представить зависимости силовых характеристик доски от ее пространственных и кинематических характеристик в табличном виде и использовать таблицы в процессе моделирования траловой системы в реальном времени взамен полного вычисления гидродинамики траловой доски по системам уравнений Навье - Стокса на каждом шаге моделирования, значительно уменьшив тем самым трудоемкость вычислений. Для построения таблиц предлагается использовать метод покоординатного расщепления уравнений Навье - Стокса, позволяющий вычислять поле давлений по полю скоростей на поверхностях траловой доски. Для каждого типа траловой доски для всех возможных комбинаций линейных и угловых скоростей с заданным шагом дискретизации вычисляется поле давлений, а затем и внешние силы, приложенные к траловой доске. Полученная шестимерная таблица ставит в соответствие скоростям внешние силы.
1. Недоступ А. А., Ражев А. О., Соколова Е. В., Макаров В. В. Математическое моделирование орудий и процессов рыболовства. Ч. 3: моногр. Калининград: Изд-во КГТУ, 2016. 184 с.
2. Недоступ А. А., Ражев А. О., Соколова Е. В. Моделирование озерных закидных неводов методом конечных элементов // Рыбное хозяйство. 2016. № 2. С. 96-99.
3. Недоступ А. А., Ражев А. О. Расчет сил гидродинамического сопротивления сетных орудий рыболовства в задачах имитационного их моделирования // Изв. Калининград. гос. техн. ун-та. 2016. № 42. С. 185-192.
4. Марчук Г. И. Методы расщепления для решения нестационарных задач // Журнал вычислительной математики и математической физики. 1995. Т. 35, № 6. С. 843-849.
5. Березин С. Б., Пасконов В. М., Сахарных Н. А. Моделирование трехмерных течений методом расщепления с использованием параллельной архитектуры ГПУ // Вычислительные методы и программирование: новые вычислительные технологии. 2012. № 2. С. 75-81.
6. Zink J., Pettineo M., Hoxley J. Practical rendering and computation with Direct3D 11. CRC Press, 2011. 648 p.
7. Антоник В. Г., Срочко В. А. Решение задач оптимального управления на основе методов линеаризации // Журнал вычислительной математики и математической физики. 1992. Т. 32, № 7. С. 979-991.
8. Каханер Д., Моулер К., Нэш С. Численные методы и программное обеспечение. М.: Мир, 1998. 575 с.
