УДК 638.14.015
05.00.00 Технические науки
АДАПТИВНАЯ ТЕХНОЛОГИЯ ЗИМНЕГО ЭЛЕКТРООБОГРЕВА ПЧЕЛ
Оськин Сергей Владимирович д.т.н., профессор, заведующий кафедрой электрических машин и электропривода, БРШ-код: 2746-7547
Потапенко Людмила Владимировна аспирантка кафедры электрических машин и электропривода, БРШ-код: 1244-5113
Овсянников Дмитрий Алексеевич
к.т.н., доцент кафедры электрических машин и
электропривода,
БРШ-код: 2802-4906
Ильченко Яков Андреевич
к.т.н., доцент кафедры электрических машин и
электропривода,
БРШ-код: 8306-5235
Кубанский государственный аграрный университет, Краснодар, Россия
В статье рассматриваются вопрос повышения рентабельности труда в пчеловодстве с помощью внедрения средств механизации и автоматизации основных технологических процессов. Наиболее сложным и ответственным периодом в пчеловодстве, который влияет на дальнейшую продуктивность пчелиной семьи, является зимний период. Особенностью зимнего состояния пчел является их способность к агрегации и поддержанию внутренней температуры на уровне 24...32 °С при широких колебаниях наружной температуры. Для облегчения прохождения зимнего периода применяется электрообогрев, для оптимального режима работы которого необходимо учитывать изменяющуюся геометрию зимнего пчелиного клуба. На параметры микроклимата в улье влияет активность пчел, что следует учитывать при разработке системы управления электрообогревом, которая должна быть адаптивной. Проведенные исследования позволяют получить основные уравнения, описывающие состояние микроклимата в улье с учетом изменяющейся геометрии зимнего клуба пчел. Уравнения связывают основные параметры клуба пчел, характеристики улья, состояние внутреннего воздуха с температурой наружного воздуха. Результаты моделирования могут использоваться при программировании контроллеров управления систем адаптивного электрообогрева, а также в специализированном программном обеспечении для моделирования состояния микроклимата в улье
UDC 638.14.015 Technical sciences
WINTER BEEHIVES ELECTRIC HEATING ADAPTIVE TECHNOLOGY
Oskin Sergey Vladimirovich Doctor of technical sciences, professor, head of department of electrical machines and drives, SPIN-code: 2746-7547
Potapenko Ludmila Vladimirovna, Postgraduate student of department of electrical machines and drives, SPIN-code: 1244-5113
Ovsyannikov Dmitriy Alexeyevich Cand.Tech.Sci., associate professor of electrical machines and drives, SPIN-code: 2802-4906
Ilchenko Iakov Andreevich Cand.Tech.Sci., associate professor of electrical machines and drives, SPIN-code: 8306-5235
Kuban State Agrarian University, Krasnodar, Russia
The article is devoted to increasing beekeeping's earning capacity with rigging and automation of technological processes. Winter is the most complicated and critical period in beekeeping. It influences on further productivity of bee colony. Winter period has specificity - bees are gathering in winter cluster and sustaining temperature in 24...32 ^ range inside. This temperature is depending on fluctuation of outside air temperature. For passing winter without complication electric heating is used. Cluster geometry in winter is changing constantly and for optimal regime of electric heating we need to take it into account. Additionally, bees' activity changing beehive's microclimate and it's important to know when beehive's microclimate adaptive control system is developed. Conducted researches allow us to get basic equations which characterizes microclimate with taking into account of winter cluster's variable geometry. The equations establish linkage basic parameters of beehive, winter cluster, internal air between outdoor air. Simulation results may be used for adaptive heating control system programming. Additionally, it may be used in special-purpose software for beehive's microclimate simulation
Ключевые слова: РЕНТАБЕЛЬНОСТЬ ТРУДА, ПЧЕЛОВОДСТВО, ЭЛЕКТРОТЕХНОЛОГИЯ, МИКРОКЛИМАТ, ЗИМНИЙ КЛУБ, МОДЕЛИРОВАНИЕ, ЭНЕРГЕТИЧЕСКАЯ ХАРАКТЕРИСТИКА
Ро1: 10.21515/1990-4665-132-022
Keywords: EARNING CAPACITY, BEEKEEPING, ELECTRICAL TECHNOLOGY, MICROCLIMATE, WINTER CLUSTER, COMPUTER SIMULATION, ENERGY CHARACTERISTIC
Пчеловодство в России имеет отличительную особенность - при имеющемся большом количестве пчелиных семей очень низкая степень механизации и автоматизации основных технологических процессов. Резко изменить такую ситуацию невозможно, но можно выбрать наиболее нуждающиеся процессы и направить туда научные и инженерные разработки. Для этого нужно провести соответствующие исследования и представить пчеловодам доказательную базу эффективности предлагаемых разработок. Конечно, большая часть пчеловодов ждет отдачу от вложенных средств как можно раньше, то есть за 1 год. Пчеловоды, которые имеют более 100 семей, готовы на более долговременную перспективу. Проведенный анализ рентабельности труда пчеловодов показывает, что сейчас она очень низкая, но ее можно увеличить, если действовать системно [1, 6].
В технологии пчеловождения на юге России очень важной составляющей является зимовка. От успешности прохождения этого периода зависит продуктивность пчелиной семьи, особенно на весенних медосборах. Если после зимовки пчелиная семья сильно ослабла (много подмора), то она не успевает нарастить необходимое количество рабочих пчел к периоду цветения медоносных растений. Публикуется множество исследований, посвященных анализу этого сложного периода состояния пчел. Более тщательный анализ, с использованием современных программных продуктов, проводил Тобоев В. А. и Еськов Е.К. [2, 3, 4]. Особенностью зимнего состояния пчел является их способность к агрегации и поддержании внутренней температуры на уровне 24...32 °С при широких колебаниях наружной температуры. В конце зимы, при появлении расплода пчелы, подни-
мают температуру до 34°С. Поддержание таких температур происходит на основе потребления пчелами меда и постоянной вентиляции зимнего клуба. Не смотря на большое количество исследований в этой области, этот процесс остается до конца не изученным. Необходимость полного понимания жизнедеятельности пчел диктуется задачами повышения медопродуктив-ности отечественного пчеловодства. До сих пор нет твердых рекомендаций по количеству оставляемого меда пчелам на зиму, какая должна быть вентиляция и т.д. Также важно как проводить электрообогрев пчел, в каком режиме. Последнее важно из-за того, что при избыточном обогреве повышается потребление меда и можно спровоцировать пчелиную матку на яйцекладку в холодные дни. Таким образом, нужно продолжать исследования с целью создания оптимальных условий пчелиным семьям в зимний период и не допускать дополнительное потребление корма.
Есть большие различия в зимовке пчелиных семей в зависимости от региона нашей страны. Так, например, отличительной особенностью прохождения такого периода на Северном Кавказе является: поздняя агрегация пчел в клуб, большие колебания температуры окружающего воздуха в осенний и зимний период, раннее появление расплода.
Для анализа физических, процессов протекающих в улье мы приняли в качестве базовой семьи скопление 15000 пчел, расположенных в двенна-дцати-рамочном улье дадановского типа. Известно, что в зимой пчелы собираются в тело шарообразной формы. Мы приняли, что клуб размещается в центре улья и скопление пчел имеет диаметр 220 мм, что соответствует занятию ими 6 улочек. Основные геометрические параметры и пространственное размещения основных элементов представлены на рисунке 1.
Рисунок 1.- Геометрические параметры основных элементов улья и пчелиного зимнего клуба.
Как известно из литературы, плотность клуба пчел может изменяться практически в два раза. Естественно, что максимальная плотность клуба будет при низких температурах, а минимальная - при высоких. Принимаем максимальную плотность клуба при температуре окружающего воздуха -30°С, минимальную плотность - при 0°С. Принимаем линейную зависимость плотности пчел от наружной температуры и тогда уравнение будет выглядеть следующим образом:
рпч = 243 - 8 • Т, кг / мъ
Соответственно коэффициент теплопроводности также будет изменяться линейно согласно уравнению:
1пч = 0,076 - 0,0017 • Т, Вт / м • К
Масса 15 тысячи пчел находится на уровне 1,5 кг. Средняя масса одной пчелы 110-120 г. Так как плотность клуба изменятся, то будет изменяться и радиус клуба в зависимости от температуры окружающего воздуха. Объем клуба пропорционален квадрату радиуса (отдельные цилиндры),
следовательно, при максимальной плотности пчел клуба радиус уменьшиться до 0,71 от первоначального значения. В соответствии с принятыми геометрическими значениями радиусы отдельных улочек будут описываться уравнениями:
Я = 1,07 • Т + 110, мм
1 пч ? ?
Я2 = 0,93 • Т + 95, мм
2 пч ? ?
Я3 = 0,7 • Т + 70, мм
3 пч ? ?
Я4 = 0,3 • Т + 35, мм
4 пч ? ?
где Т - температура окружающего улей воздуха.
Для нашего случая при 15000 пчел энергетическая характеристика будет выглядеть:
Рпчу = -0,3 • Т + 4, Вт
Также энергетику пчел можно связать и с потреблением корма. Известно, что при потреблении меда 0,32 г/ч выделяется мощность 1 Вт. Тогда о расход корма данного клуба в зимний период в зависимости от температуры имеет вид, г/ч
О = -0,1 • Т +1,2
ку ' '
Эти уравнения справедливы для диапазона температур от -5 до +10 °С и когда пчелы не заложили расплод.
Так как для дальнейшего анализа потребуется удельная мощность воспроизводимая пчелами, то разделим полученное выражение на объем занимаемый клубом:
Рпчууд= Рг = -60 • Т + 800, Вт / м3
кл
Аппроксимирующее уравнение теряемой тепловой мощности через вентиляцию будет выглядеть, Вт:
Р = -0,22 -Т + 0,8
венту ' '
Необходимо определить коэффициент теплоотдачи с верхней части поверхности клуба. Известно уравнение для определения мощности необходимой для вентиляции:
Рвент Свозд Р возд Овент
С другой стороны также известно уравнение тепловых потерь с поверхности за счет конвекции:
Рпот = а- 8Пов ■ Л г ,
где а - коэффициент теплоотдачи,
Вт/м-К; Я - площадь теплоот-
дающей поверхности, м2.
В литературе известны следующие границы изменения коэффициента теплоотдачи: 6,0-10-4...3,0-10-3 Вт/(см2-°С).
Приравняем эти мощности для нашего случая, так как пчелы вынуждены отдать эту энергию в окружающее пространство. Тогда получим уравнение для определения коэффициента теплоотдачи:
а- 5 -Лг = с д-р д О -Лг ® а = Свозд'Рвозд' Овент
пов возд г возд вент
Я«
Расход воздуха определяется по формуле:
„ Цн-Опч
О
I -1
вых. в
3
где 1вьа - абсолютное влагосодержание выходящего воздуха, г/м ; 1в -
3
абсолютное влагосодержание входящего воздуха, г/м , дН2Опч - количество
воды, выделившееся в результате окисления корма пчелами,
Чнрпч = 0,68 • , г/ч.
В среднем разница абсолютных влагосодержаний находится на уровне 2,5 г/ м3. Тогда формула для определения расхода воздуха принимает вид:
а- - ^ = 0,27 • а, ч
Подставим выражение расхода корма для улья в зимний период и получим зависимость расхода воздуха от температуры:
а = 0,68 •^ = 0,27 • (_0,1 • т +1,2 ) = -0,027 • Т + 0,324
После соответствующих подстановок получим формулу связывающую коэффициент теплоотдачи и температуру наружного воздуха:
с д • Р о • О с д • р о • (-0,027 • Т + 0,324)
а = возо ' возо ^вент _ возо ' возо I ' ' /
У6 5
/—ч у
Для нашего примера формула принимает вид:
своз0 • Рвоз0 • (-0,027 • Т + 0,324)_ 1007 • 1,29 • (-0,027 • Т + 0,324)
ул
а
У 65 (2 • 3,8 + 2 • 3,3 + 2 • 2,4) • 3,6
£и1 Ул
=19 • (-0,027 • Т + 0,324) = -0,513 • Т + 6,156
Используя полученные уравнения можно приступать к моделированию теплового состояния пчел в различный период зимы. Изменяя температуру окружающего воздуха, можно будет обнаружить или недостаток обогрева или наоборот его избыток.
Для более точных расчетов температурного режима улья нужно учесть изменение энергетической характеристики пчел при температурах близких к нулю. Известно, что в осенне-зимний период энергетическая характеристика зимнего клуба пчел близка к параболе [1] (рис.2).
Рисунок 2 - Зависимости тепловыделений пчелиных семей от температуры окружающего воздуха (для 1000 пчел)
Учитывая это, получено уточненное уравнение аппроксимации энергетической характеристики 15000 пчел зимней агрегации, Вт:
РЬее= 0,016 • Т2 - 0,1 • Т + 4,61
где Т - температура окружающего улей воздуха.
Тогда расход корма в зимний период 0Ьее для 15000 пчел, в зависимости от температуры имеет вид, г/ч:
вЬее = 0,005 • Т2 - 0,032 • Т + 1,475
Так как для дальнейшего анализа потребуется удельная мощность воспроизводимая пчелами, то разделим полученное выражение на объем занимаемый клубом в данном примере [1]:
0
Р*.. 0,016 • Т--0,1 • Т + 4,61
_ Ьее _
Ьее
V
Ьее
0,005
= 3,2 • Т2 - 20 • Т + 922, Вт / м3
где УЬее - объем занимаемый зимним клубом, для данного примера 0,005 м3.
Известно уравнение для определения необходимого расхода воздуха пчелами для удаления образующейся в улье влаги [1]:
Q
¿--ail
А - А
out in
где qHo0 - количество воды, выделавшееся в результате окисления корма, г/ч; Аои{, Ат - абсолютное влагосодержание соответственно выходного и
входного воздуха, г/м .
Количество воды, выделавшееся в результате окисления корма, можно определить по формуле [1]:
qH 2о =0,68 • Gbee
Тогда, учитывая ранее полученные уравнения, выражение для расхода воздуха данного клуба пчел будет иметь вид:
0,0034 • Т2 -0,0216 • Т +1
_ А - А
out Bin
Относительная влажность выходящего из клуба воздуха практически составляет 100%. По данным литературных источников температура выходного воздуха из скопления пчел равна температуре корки клуба и находится на уровне от 10 °С до 15 °С. Тогда абсолютная влажность выходного воздуха находится на уровне 9,4-12,8 г/м . Естественно меньшее значение будет соответствовать минимальной температуре наружного воздуха (30 С), большее значение - максимальной температуре (+15°С). Принимаем, что зависимость абсолютной влажности выходного воздуха от наружной температуры имеет линейный вид. Статистика метеоданных по Краснодарскому краю показывает, что средняя относительная влажность воздуха в зимний период составляет 80%. На основании таблиц перевода относительной влажности в абсолютную (с учетом температуры) получено следующее аппроксимирующее уравнение, м /ч:
0,0034 • Т2 -0,0216 • Т +1 __0,0034 • Т2 - 0,0216 • Т +1_
А.100 - Ап80 " (0,007 • Т +11,5) - (0,0065 • Т2 + 0,3 • Т + 4,03) Следовательно, все основные параметры микроклимата в улье зависят от окружающей температуры. Таким образом, геометрические и теп-лофизические характеристики зимнего клуба пчел можно представить следующей системой уравнений:
R Ьее = 1,07 • Т + 110; R2 bee = 0,93 • Т + 95; R, bee = 0,7 • Т + 70; R, bee= 0,3 • Т + 35
pbee = 243 - 8 • Т;1ee = 0,076 - 0,0017 • Т dar _ 1,25 + 0,025 • Т РЫе = 0,016 • Т2 - 0,1 • Т + 4,61
__0,0034 • Т2 - 0,0216 • Т +1_
Qair _ (0,007 • Т +11,5) - (0,0065 • Т2 + 0,3 • Т + 4,03)
Как видно из данной группы уравнений все основные параметры микроклимата в улье зависят от окружающей температуры. Основные параметры микроклимата внутри клуба пчелы поддерживают за счет изменения его геометрии, внутренней вентиляции, перемещения пчел из более нагретой части к периферии и назад. Необходимость усиленной вентиляции связано с потребностью удаления влаги, образующейся в результате поедания меда. При этом источником энергии является мед. При понижении температуры окружающего воздуха пчелы вынуждены больше есть меда, что приводит к выделению большего количества влаги и, следовательно, им приходится увеличивать вентиляцию, которая выносит и тепло. В зимний период пчелы находятся в пассивном состоянии, что позволяет им резко сократить энергетические затраты. Полученные уравнения можно использовать для моделирования состояния микроклимата в улье как в специализированных программных продуктах, так отдельно разработанных. Результаты моделирования будут использоваться при программировании микроконтроллеров, используемых в системах электрообогрева.
Литература
1. Оськин С.В., Овсянников Д. А. Электротехнологические способы и оборудование для повышения производительности труда в медотоварном пчеловодстве северного Кавказа: монография/ С.В.Оськин, Д.А.Овсянников.- Краснодар: Изд-во «Крон», 2015. -198 с.
2. Еськов Е.К., Тобоев В.А. Сезонная динамика тепловых процессов в межсотовых скоплениях зимующих пчел Apis mtllifera/ Е.К. Еськов, В.А. Тобоев//Зоол. Жур-нал.-2011.Т.90, №3.-С.335-341.
3. Еськов Е.К., Тобоев В. А. Математическое моделирование распределения температурных полей в холодовых агрегациях насекомых/ Е.К. Еськов, В.А. Тобо-ев//Биофизика.-2009.Т.54.Вып. 1.- с.114-119.
4. Тобоев В. А., Толстов М. С. Моделирование тепловых процессов в скоплениях зимующих пчел/В.А. Тобоев, М.С. Толстов//Физические процессы в биологических системах. Июнь. -2014.- с.97-102.
5. Оськин С.В. и др. Теплофизическое моделирование пчелиного улья - составная часть безопасной зимовки/С.В. Оськин, Л.В. Потапенко, Е.В. Пустовойтова, Н.Г. Рассолов// Чрезвычайные ситуации: промышленная и экологическая безопасность; Вып. №29 (1)-Краснодар 2017.-С. 115-121.
6. Оськин С.В., Овсянников Д.А. Внедрение электротехнологий в пчеловодство/ С.В. Оськин, Д.А.Овсянников// Физико-технические проблемы создания новых технологий в агропромышленном комплексе : сборник научных статей по материалам Международной научно-практической конференции. - Ставрополь : СЕКВОЙЯ, 2017. -С.201-206.
References
1. Os'kin S.V., Ovsjannikov D.A. Jelektrotehnologicheskie sposoby i oborudovanie dlja povyshenija proizvoditel'nosti truda v medotovarnom pchelovodstve severnogo Kavkaza: monografija/ S.V.Os'kin, D.A.Ovsjannikov.- Krasnodar: Izd-vo «Kron», 2015. -198 s.
2. Es'kov E.K., Toboev V.A. Sezonnaja dinamika teplovyh processov v mezhsotovyh skoplenijah zimujushhih pchel Apis mtllifera/ E.K. Es'kov, V.A. Toboev//Zool. Zhurnal.-2011.T.90, №3.-s. 335-341.
3. Es'kov E.K., Toboev V.A. Matematicheskoe modelirovanie raspredelenija tem-peraturnyh polej v holodovyh agregacijah nasekomyh/ E.K. Es'kov, V.A. Tobo-ev//Biofizika.-2009.T.54.Vyp.1.- s. 114-119.
4. Toboev V.A., Tolstov M.S. Modelirovanie teplovyh processov v skoplenijah zimujushhih pchel/V.A. Toboev, M.S. Tolstov//Fizicheskie processy v biologicheskih siste-mah. Ijun'. -2014.- s. 97-102.
5. Os'kin S.V. i dr. Teplofizicheskoe modelirovanie pchelinogo ul'ja - sostavnaja chast' bezopasnoj zimovki/S.V. Os'kin, L.V. Potapenko, E.V. Pustovojtova, N.G. Rassolov// Chrez-vychajnye situacii: promyshlennaja i jekologicheskaja bezopasnost'; Vyp. №29 (1)-Krasnodar 2017.- s. 115-121.
6. Os'kin S.V., Ovsjannikov D.A. Vnedrenie jelektrotehnologij v pchelovodstvo/ S.V. Os'kin, D.A.Ovsjannikov// Fiziko-tehnicheskie problemy sozdanija novyh tehnologij v ag-ropromyshlennom komplekse : sbornik nauchnyh statej po materialam Mezh-dunarodnoj nauchno-prakticheskoj konferencii. - Stavropol' : SEKVOJJa, 2017. - s. 201-206.