MATHEMATICAL MODELING OF DRIFTING NET OPERATIONS IN THE RIVER DISTURBANCE
Authors:
1.
Kaliningrad State Technical University
(Candidate of Technical Sciences, Assistant Professor; Head of the Department of Industrial Fishery)
2.
Kaliningrad State Technical University
(Candidate of Sciences in Technology; Leading Researcher)
Type:
Article
Pages:
from 45
to 53
Status:
Published
Received:
05.09.2016
Accepted:
25.09.2016
Published:
25.09.2016
Subject area:
UDK 639.2.081.113
GRNTI 34.39 Физиология человека и животных
GRNTI 62.13 Биотехнологические процессы и аппараты
GRNTI 69.01 Общие вопросы рыбного хозяйства
GRNTI 69.25 Аквакультура. Рыбоводство
GRNTI 69.31 Промышленное рыболовство
GRNTI 69.51 Технология переработки сырья водного происхождения
GRNTI 87.19 Загрязнение и охрана вод суши, морей и океанов
GRNTI 34.39 Физиология человека и животных
GRNTI 62.13 Биотехнологические процессы и аппараты
GRNTI 69.01 Общие вопросы рыбного хозяйства
GRNTI 69.25 Аквакультура. Рыбоводство
GRNTI 69.31 Промышленное рыболовство
GRNTI 69.51 Технология переработки сырья водного происхождения
GRNTI 87.19 Загрязнение и охрана вод суши, морей и океанов
Language:
Russian
Keywords:
drifting net, modeling, dynamics, disturbance, material point
Abstract (English):
Drifting net is an enmeshing gear of inland and coastal fisheries. There are three ways of nets operations: float on the surface, mid-water float and bottom float. It is necessary to calculate the parameters of the net so that at different velocities and characteristics of disturbance its form and disclosure meet the conditions of the fishery. The paper proposes a discrete composite mathematical model of the surface single-walled drifting net, allowing, in accordance with the given structure and parameters of the network elements, to determine its power and spatial-temporal characteristics under the river disturbance. The model is based on the composition of Stokes - Relay’s wave model and discrete model of material point. While releasing and drifting of the net, the corresponding hydrodynamic and hydrostatic forces effecting different elements are calculated, the speed of the different elements and their replacements are determined. Based on the proposed mathematical model of the drifting net a numerical experiment that lasted 120 seconds was carried out. During the first 70 seconds the experiment was carried out without disturbance, during the rest 50 seconds the effect of disturbance at ceteris paribus was taken into account. The wave parameters: wave height h = 0.25 m, wave length λ = 5 m. As a result of the experiment the dependence of force and geometrical characteristics of the net on the time of the modeling process was obtained. It is shown that under disturbance release and, therefore, catching capacity are increased, drifting rate is increased as well, but the power of hydrodynamic resistance is decreased.
Drifting net is an enmeshing gear of inland and coastal fisheries. There are three ways of nets operations: float on the surface, mid-water float and bottom float. It is necessary to calculate the parameters of the net so that at different velocities and characteristics of disturbance its form and disclosure meet the conditions of the fishery. The paper proposes a discrete composite mathematical model of the surface single-walled drifting net, allowing, in accordance with the given structure and parameters of the network elements, to determine its power and spatial-temporal characteristics under the river disturbance. The model is based on the composition of Stokes - Relay’s wave model and discrete model of material point. While releasing and drifting of the net, the corresponding hydrodynamic and hydrostatic forces effecting different elements are calculated, the speed of the different elements and their replacements are determined. Based on the proposed mathematical model of the drifting net a numerical experiment that lasted 120 seconds was carried out. During the first 70 seconds the experiment was carried out without disturbance, during the rest 50 seconds the effect of disturbance at ceteris paribus was taken into account. The wave parameters: wave height h = 0.25 m, wave length λ = 5 m. As a result of the experiment the dependence of force and geometrical characteristics of the net on the time of the modeling process was obtained. It is shown that under disturbance release and, therefore, catching capacity are increased, drifting rate is increased as well, but the power of hydrodynamic resistance is decreased.
Keywords:
drifting net, modeling, dynamics, disturbance, material point
drifting net, modeling, dynamics, disturbance, material point
Введение Лов ставными рыболовными сетями в местах с сильным течением, особенно там, где вода несет много посторонних предметов, бывает малоэффективным. Сильное течение и волнение деформируют рыболовную сеть, создавая большое натяжение нитей сетного полотна, что делает сеть малоуловистой. В таких местах эффективнее использовать плавные рыболовные сети, которые свободно сплывают по течению навстречу рыбе, поднимающейся вверх по реке. Лов плавными сетями применяют в низовьях рек Волги и Дона для лова частиковых рыб, на р. Амур для лова лососевых, на р. Обь во время весеннего нерестового хода нельмы, муксуна, сырка, пыжьяна, стерляди, язя и других ценных промысловых рыб. Этот лов применяют и на таких небольших реках, как Кура и Печора, т. к. в этих водоемах он эффективнее, чем лов ставными сетями. Плавная сеть представляет собой объячеивающее орудие внутреннего и прибрежного рыболовства. Различают три способа эксплуатации сетей: плав по поверхности, плав разноглубинный и плав донный [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. Baranov F. I. Teoriya i raschet orudiy rybolovstva. M.: Pischepromizdat, 1948. 436 s.
2. Baranov F. I. Tehnika promyshlennogo rybolovstva. M.: Pischepromizdat, 1960. 696 s.
3. Nedostup A. A., Razhev A. O. Razrabotka matematicheskoy modeli kryla stavnogo podvesnogo nevoda na volnenii i v usloviyah shtorma // Izv. Kaliningrad. gos. tehn. un-ta. 2014. № 35. S. 66-74.
4. Nedostup A. A., Volod'ko D. A., Razhev A. O. Gidrodinamicheskiy raschet rybolovnoy seti // Sb. tr. V Mezhdunar. nauch.-prakt. konf. «Inzhenernye sistemy-2012». M.: RUDN, 2012. S. 27-31.
5. Chebotarev A. I. Obschaya gidrologiya (vody sushi): ucheb. posobie. L.: Gidrometeoizdat, 1975. 530 s.
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. Nedostup A. A. Metody rascheta setnyh passivnyh orudiy vnutrennego i pribrezhnogo rybolovstva: monogr. Kaliningrad: Izd-vo KGTU, 2010. 280 s.
8. Nedostup A. A., Naumov V. A., Razhev A. O., Belyh A. V. Matematicheskoe modelirovanie orudiy i processov rybolovstva: monogr. Kaliningrad: Izd-vo KGTU, 2013. Ch. 1. 253 s.
9. Il'in V. A., Sadovnichiy V. A., Sendov B. H. Matematicheskiy analiz / pod red. A. N. Tihonova. M.: Prospekt, 2004. Ch. 1. 660 s.
10. Patterson David A., Hennessy John L. Computer Organization and Design, Fourth Edition: The Hardware/Software Interface. Morgan Kaufmann, 2011. 914 p.
11. Dorofeev S. N. Geometricheskie preobrazovaniya v primerah i zadachah: ucheb. posobie. Penza: Inform.-izd. centr PGU, 2002. 189 c.
