МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАБОТЫ ПЛАВНОЙ СЕТИ В УСЛОВИЯХ ВОЛНЕНИЯ РЕКИ
Авторы:
1.
Калининградский государственный технический университет
(канд. техн. наук, доцент; зав. кафедрой промышленного рыболовства )
2.
Калининградский государственный технический университет
(младший научный сотрудник УНИД )
Тип:
Статья
Страницы:
с 45
по 53
Статус:
Опубликован
Получено:
05.09.2016
Одобрено:
18.01.2023
Опубликовано:
25.09.2016
Классификаторы:
УДК 639.2.081.113
ГРНТИ 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 Загрязнение и охрана вод суши, морей и океанов
Язык материала:
русский
Ключевые слова:
плавная сеть, моделирование, динамика, волнение, точечная масса
Аннотация (русский):
Плавная сеть представляет собой объячеивающее орудие внутреннего и прибрежного рыболовства. Различают три способа эксплуатации сетей: плав по поверхности, плав разноглубинный и плав донный. Необходимо рассчитать параметры сети так, чтобы при различных значениях скорости течения и характеристиках волнения ее форма и раскрытие удовлетворяли условиям промысла. Предлагается дискретная композитная математическая модель поверхностной одностенной плавной сети, позволяющая по заданной конструкции и параметрам элементов сети определить ее силовые и пространственно-временные характеристики с учетом волнения. Модель основана на композиции волновой модели Стокса - Релея и дискретной модели точечных масс. В процессе раскрытия и сплава сети рассчитываются соответствующие гидродинамические и гидростатические силы, действующие на различные элементы, определяются скорость различных элементов и их перемещение. На основе предложенной математической модели плавной сети был проведен численный эксперимент продолжительностью 120 секунд. В течение первых 70 секунд эксперимент проводился в отсутствие волнения, в течение остальных 50 секунд учитывалось влияние волнения при прочих равных условиях. Параметры волны: высота волны h = 0,25 м; длина волны λ = 5 м. В результате эксперимента были получены графики зависимостей силовых и геометрических характеристик сети от времени процесса моделирования. Показано, что при волнении раскрытие и, следовательно, уловистость сети увеличиваются, увеличивается также скорость плава, а сила гидродинамического сопротивления уменьшается.
Плавная сеть представляет собой объячеивающее орудие внутреннего и прибрежного рыболовства. Различают три способа эксплуатации сетей: плав по поверхности, плав разноглубинный и плав донный. Необходимо рассчитать параметры сети так, чтобы при различных значениях скорости течения и характеристиках волнения ее форма и раскрытие удовлетворяли условиям промысла. Предлагается дискретная композитная математическая модель поверхностной одностенной плавной сети, позволяющая по заданной конструкции и параметрам элементов сети определить ее силовые и пространственно-временные характеристики с учетом волнения. Модель основана на композиции волновой модели Стокса - Релея и дискретной модели точечных масс. В процессе раскрытия и сплава сети рассчитываются соответствующие гидродинамические и гидростатические силы, действующие на различные элементы, определяются скорость различных элементов и их перемещение. На основе предложенной математической модели плавной сети был проведен численный эксперимент продолжительностью 120 секунд. В течение первых 70 секунд эксперимент проводился в отсутствие волнения, в течение остальных 50 секунд учитывалось влияние волнения при прочих равных условиях. Параметры волны: высота волны h = 0,25 м; длина волны λ = 5 м. В результате эксперимента были получены графики зависимостей силовых и геометрических характеристик сети от времени процесса моделирования. Показано, что при волнении раскрытие и, следовательно, уловистость сети увеличиваются, увеличивается также скорость плава, а сила гидродинамического сопротивления уменьшается.
Ключевые слова:
плавная сеть, моделирование, динамика, волнение, точечная масса
плавная сеть, моделирование, динамика, волнение, точечная масса
Введение Лов ставными рыболовными сетями в местах с сильным течением, особенно там, где вода несет много посторонних предметов, бывает малоэффективным. Сильное течение и волнение деформируют рыболовную сеть, создавая большое натяжение нитей сетного полотна, что делает сеть малоуловистой. В таких местах эффективнее использовать плавные рыболовные сети, которые свободно сплывают по течению навстречу рыбе, поднимающейся вверх по реке. Лов плавными сетями применяют в низовьях рек Волги и Дона для лова частиковых рыб, на р. Амур для лова лососевых, на р. Обь во время весеннего нерестового хода нельмы, муксуна, сырка, пыжьяна, стерляди, язя и других ценных промысловых рыб. Этот лов применяют и на таких небольших реках, как Кура и Печора, т. к. в этих водоемах он эффективнее, чем лов ставными сетями. Плавная сеть представляет собой объячеивающее орудие внутреннего и прибрежного рыболовства. Различают три способа эксплуатации сетей: плав по поверхности, плав разноглубинный и плав донный [1, 2]. Целью исследования была разработка математической модели донной одностенной плавной сети, позволяющей по заданной конструкции и параметрам элементов сети определить ее силовые и пространственно-временные характеристики с учетом волнения. Математическая модель Донная плавная сеть состоит из сетной оболочки, посадочных нитей, верхней, нижней и боковых подбор, оснастки верхней подборы (плава), оснастки нижней подборы (загрузки), угловых огонов сети, оттяжек, поводцов и вешек. Вид донной плавной сети в плоскости OXZ показан на рис. 1. В процессе раскрытия и сплава сети рассчитываются соответствующие гидродинамические и гидростатические силы, действующие на различные элементы, определяется скорость различных элементов и их перемещение. Необходимо рассчитать параметры сети так, чтобы при различных значениях скорости течения и характеристиках волнения ее форма и раскрытие удовлетворяли условиям промысла. Рис. 1. Вид донной плавной сети в плоскости OXZ: Q - верхняя подбора; Qv - вешка; G - нижняя подбора; K - точка соединения оттяжки и поводцов; Т - сила натяжения Рассмотрим экстремальный случай, когда фронт волны распространяется противоположно направлению сплава сети. В указанных условиях задачу можно свести к двумерному случаю, рассматривая процесс в плоскости, перпендикулярной течению и фронту волны. В данной плоскости, согласно волновой модели Стокса - Релея, описанной в [3], частицы воды у поверхности движутся по эллиптическим траекториям (рис. 2). Рис. 2. Профиль волны, соответствующий модели Стокса - Релея Вектор скорости воды (u, w) относительно дна водоема в точке (x0, z0) является суммой трех составляющих: , где v - скорость вынужденного течения (речного); uв - скорость волнового течения; uw и ww - горизонтальная и вертикальная проекции скорости вращения воды, вызванного волной, на оси декартовой системы координат. Проекции скорости вращения воды, вызванного волной, на оси декартовой системы координат определяются системой параметрических уравнений [4]: где h - высота волны; c - скорость распространения волны; g - ускорение свободного падения; k = 2π/λ; λ - длина волны; H - глубина места лова; z0 - высота от дна водоема; tw - промежуточная величина системы параметрических уравнений. Скорость распространения волны определяется зависимостью , R = 1/k. Скорость волнового течения [5]: , где uвп = 0,25h2k2c - скорость ветрового течения на поверхности. Скорость вынужденного течения на высоте от дна водоема z0 определяется зависимостью , где vп - скорость вынужденного течения на поверхности; C - коэффициент формулы Шези, изменяющийся для рек в пределах 10-60 м1/2/с. При создании математической модели необходимо учесть, что элементы плавной сети могут находиться на поверхности воды. В связи с этим необходимо знать зависимость профиля волны от времени в декартовой системе координат (x, z). Данная зависимость определяется системой параметрических уравнений: где t - модельное время; a0 и b0 - горизонтальная и вертикальная полуоси эллипса (рис. 2); d0 = 0,25πh2/λ - поднятие центра орбит, вызванное несимметричностью волны относительно линии, определяющей положение воды в покое [5]. Согласно волновой модели Стокса - Релея, размер полуосей эллипса определяется по следующим формулам: При расчете сил гидродинамического сопротивления необходимо учитывать скорость движения элементов плавной сети относительно воды. Результирующий вектор мгновенной скорости в произвольной точке невода определяется как разность вектора скорости движения этой точки и вектора движения воды относительно дна водоема (u, w). Рассмотрим процесс дискретного моделирования плавной сети. Разобьем сеть на множество дискретных элементов, обрабатываемых по одному алгоритму. Представим плавную сеть в виде пяти групп элементов: сетное полотно; верхняя подбора; нижняя подбора; оттяжки и поводцы; вешки. Выведем зависимости, связывающие силовые и геометрические характеристики элементов сети по каждой группе [3]. Для элементов сетного полотна, оттяжек и поводцов в модели точечных масс [6] проекции сил гидродинамического сопротивления на оси декартовой системы координат, в которой ось OX направлена по скорости движения элемента сети относительно воды, имеют следующие эвристические зависимости [7]: , , (1) где Rx, Rz - проекции сил гидродинамического сопротивления; c0 - гидродинамический коэффициент для сети, расположенной параллельно потоку воды; c90 - гидродинамический коэффициент для сети, расположенной перпендикулярно потоку воды; r - длина связи между узлами модели точечных масс; rx, rz - длина проекции элемента на оси прямоугольной декартовой системы координат, ось OX которой направлена в направлении скорости перемещения узла относительно воды; v - абсолютное значение скорости перемещения узла относительно воды; ρ - плотность воды; S - площадь участка сети; Fo - сплошность участка сети. Гидродинамические коэффициенты [8] находятся отдельно для сетного полотна и оттяжек с поводцами: , , , (2) где Re - число Рейнольдса для сети; d - диаметр ниток сети; ϑ - кинематическая вязкость воды; b - коэффициент, учитывающий тип материала сети (b = 0,165 для капроновой сети, b = 0,16 - для полиэтиленовой сети). Для вычисления гидродинамических коэффициентов используется математическая операция возведения в степень. Данная операция достаточно трудоемка при вычислениях на ЭВМ. Рассмотрим метод оптимизации вычисления сил гидродинамического сопротивления, основанный на разложении их функции от скорости в ряд Тейлора [9] и применении SIMD-инструкций [10] центрального и графического процессоров. Подставив (2) в (1) можно преобразовать (1) к виду , (3) где R - одна из проекций сил гидродинамического сопротивления; k0, k1 - постоянные для участка сети коэффициенты, зависящие от параметров посадки и материала участка сети; χ - безразмерный параметр, определяющий ориентацию участка в сети относительно вектора относительной скорости v. Коэффициент k90 одинаков для каждой проекции R. Коэффициент k0 различен для каждой проекции R. Данные коэффициенты постоянны в течение всего процесса моделирования и вычисляются перед началом работы основного цикла алгоритма. Параметр χ является переменным и вычисляется на каждой итерации алгоритма моделирования для каждой проекции R. В выражении (3) математические операции возведения переменной v в степень запишем в виде , (4) где p = 2 - b для k0 и p = 1,72 для k90. Разложим выражение (4) в ряд Тейлора четвертого порядка [9]. Общий вид разложения задается выражением , где f (i)(a) - значение производной i-го порядка разлагаемой функции в точке a. При a = 1 при разложении четвертого порядка получим (5) В выражении (5) используются только операции умножения, сложения и вычитания. Данные операции на современных процессорах выполняются аппаратно и быстро. Для преобразования длин проекций участка сети из системы координат, ось OX которой параллельна дну водоема, в систему координат, ось OX которой совпадает с направлением движения элемента, а также для обратного преобразования проекций сил гидродинамического сопротивления элементов сети используется аффинное преобразование [11]. Для подбор и вешек при расчете гидродинамических сил сопротивления используем следующие зависимости (как для цилиндра): , , , , , , где RG - сила гидродинамического сопротивления нижней подборы; dG - диаметр нижней подборы; RQ - сила гидродинамического сопротивления участка верхней подборы; dQ - диаметр верхней подборы; Rv - сила гидродинамического сопротивления вешки; dv - диаметр пенопласта вешки; hv - высота пенопласта вешки; bv = 3,0 - коэффициент для гладкого цилиндра; b1 = 4,2 - коэффициент для шестипрядных полиэтиленовых канатов [4]; k1G, k1Q, kv - поправочные коэффициенты, определяющие удлинение цилиндра; L - длина сети; ν - кинематическая вязкость воды. При расчете внешних сил, действующих на узлы сети в модели точечных масс, необходимо также учесть гидростатические гравитационные силы, действующие на ее элементы и силы (силы натяжения) взаимодействия между элементами сети. Суммарная сила натяжения на произвольном участке сети определяется из уравнения где A - сумма площадей поперечного сечения элементов на участке плоскостью, проходящей между его узлами параллельно плоскости OXZ; r0 - длина участка при T = 0 (без растяжения); E - модуль упругости материала. Составим дифференциальные уравнения движения и определим начальные условия для каждого узла: где vi - скорость узла i; ri0 - начальные координаты узла i; t - время процесса моделирования; mi - сосредоточенная масса в узле i с учетом присоединенной массы; wi - узловое ускорение; F - внешняя сила, действующая со стороны среды на участок сети; ri - координаты узла i; rj - координаты узла j; Fi - равнодействующая внешних сил, действующая со стороны среды на узел i; Ni - множество всех узлов, смежных узлу i по ребру (i, j). Найдем массы, сосредоточенные в узле, по каждой группе элементов сети: , где ρij - усредненная плотность ниток участка (i, j); Vij - объем ниток участка (i, j); m0i - собственная масса узла i (с учетом присоединенной) без участков. Для узлов сетного полотна примем m0i = 0. Для нижней подборы , где G - вес нижней подборы; ρG - усредненная плотность нижней подборы с загрузкой. Для верхней подборы и вешки , где Q - плавучесть верхней подборы или вешки; ρQ - усредненная плотность верхней подборы или вешки. На основе предложенной математической модели плавной сети был проведен численный эксперимент, который длился 120 секунд: в первые 70 секунд эксперимент проводился без волнения, в течение остальных 50 секунд - с учетом влияния волнения при прочих равных условиях. Параметры волны: высота волны h = 0,25 м; длина волны λ = 5 м. Графики зависимости различных величин от времени процесса моделирования, построенные по результатам эксперимента, приведены на рис. 3-5. Рис. 3. Динамика плава сети: Qv - вешка; К - точка соединения поводцов с оттяжкой вешки; G - нижняя подбора с грузом; Q - верхняя подбора с плавами Рис. 4. Динамика величины раскрытия сети и вертикальных координат ее элементов: H(Q - G) - высота раскрытия сети, остальные обозначения соответствуют рис. 3 Рис. 5. Динамика силы гидродинамического сопротивления и сил натяжения сети: Tv - натяжение в оттяжке вешки; Tн - натяжение в нижнем поводце; Tв - натяжение в верхнем поводце; R - сила гидродинамического сопротивления сети По графикам можно сделать вывод, что при волнении раскрытие и, соответственно, уловистость сети увеличиваются, скорость плава увеличивается, а сила гидродинамического сопротивления уменьшается. Происходящее обусловлено «подбрасыванием» сети на вершинах волны вверх, что приводит к уменьшению силы трения нижней подборы о грунт водоема. В свою очередь, уменьшение силы трения (уменьшение торможения сети) приводит к увеличению скорости плава сети. Заключение Таким образом, в ходе исследования получены следующие результаты. 1. Разработана математическая модель динамических процессов, происходящих в плавной сети в условиях волнения, основанная на композиции двух математических моделей: модели волнения Стокса - Релея и модели точечных масс. 2. Предложен метод оптимизации вычисления сил гидродинамического сопротивления на ЭВМ, основанный на разложении в ряд Тейлора. 3. Проведен численный эксперимент, по результатам которого построены графики зависимости различных величин от времени процесса моделирования. 4. Установлено, что при волнении раскрытие и, соответственно, уловистость сети увеличиваются, скорость плава увеличивается, а сила гидродинамического сопротивления уменьшается.
1. Баранов Ф. И. Теория и расчет орудий рыболовства. М.: Пищепромиздат, 1948. 436 с.
2. Баранов Ф. И. Техника промышленного рыболовства. М.: Пищепромиздат, 1960. 696 с.
3. Недоступ А. А., Ражев А. О. Разработка математической модели крыла ставного подвесного невода на волнении и в условиях шторма // Изв. Калининград. гос. техн. ун-та. 2014. № 35. С. 66-74.
4. Недоступ А. А., Володько Д. А., Ражев А. О. Гидродинамический расчет рыболовной сети // Сб. тр. V Междунар. науч.-практ. конф. «Инженерные системы-2012». М.: РУДН, 2012. С. 27-31.
5. Чеботарев А. И. Общая гидрология (воды суши): учеб. пособие. Л.: Гидрометеоиздат, 1975. 530 с.
6. Nedostup A. A., Razhev A. O. A discrete model of gill nets for static and dynamic problems // 11th International workshop - Contributions on the theory of fishing gears and related marine systems DEMAT 2013. Germany, Rostock, 2013. Vol. 8. P. 13-22.
7. Недоступ А. А. Методы расчета сетных пассивных орудий внутреннего и прибрежного рыболовства: моногр. Калининград: Изд-во КГТУ, 2010. 280 с.
8. Недоступ А. А., Наумов В. А., Ражев А. О., Белых А. В. Математическое моделирование орудий и процессов рыболовства: моногр. Калининград: Изд-во КГТУ, 2013. Ч. 1. 253 с.
9. Ильин В. А., Садовничий В. А., Сендов Б. Х. Математический анализ / под ред. А. Н. Тихонова. М.: Проспект, 2004. Ч. 1. 660 с.
10. Patterson David A., Hennessy John L. Computer Organization and Design, Fourth Edition: The Hardware/Software Interface. Morgan Kaufmann, 2011. 914 p.
11. Дорофеев С. Н. Геометрические преобразования в примерах и задачах: учеб. пособие. Пенза: Информ.-изд. центр ПГУ, 2002. 189 c.
