Технико-технологическое обеспечение подсистемы выращивания нетели
О.В. Ужик, к.т.н., Белгородская ГСХА
Один из критериев оценки целевой функции подсистемы выращивания нетели — биомасса молочной железы. Особо важным является период второй половины стельности животного. Начиная с 5—6 мес. происходит естественная активизация развития вымени, а также имеет место дополнительное наращивание его биомассы в результате возмущающего воздействия — массажа.
На рисунке приведена схема подсистемы выращивания нетели, состоящей из нескольких звеньев: главное звено — оператор; звено выращивания нетелей до 6-месячной стельности; звено выращивания нетелей с 6-месячной стельности.
Переходное состояние в данном случае — изменение физиологического состояния животного
под воздействием внешних факторов и внутренних процессов, происходящих в организме.
Действенным сигналом подсистемы выращивания нетели во второй половине стельности является массаж вымени животных, который стимулирует его рост (наращивание биомассы), что в конечном счёте должно привести к росту молочной продуктивности коровы. Это и будет служить оценкой эффективности выращивания нетели.
Предложенная нами подсистема формирования вымени нетели снабжена обратной связью, обеспечивающей оптимизацию управлений действием исполнительных машин и механизмов технологического процесса в текущий момент времени, оцениваемой по реакции как объекта воздействия, так и внешнего окружения.
При этом полагаем, что внешняя среда — субстрат, представляющий собой базовые эле-
менты технологии содержания нетели с соответствующим её физиологии режимами поения, кормления, микроклимата и т.д.
Однако, как показано на рисунке, в подсистеме формирования вымени нетели присутствуют и случайные факторы. Следует отметить, что оптимизация управлений, формируемая на опыте выращивания текущих животных, более значимой оказывается для последующих групп животных.
Каждое звено характеризуется своим выходным сигналом. Для первого звена выходным сигналом является масса вымени, формируемая в результате естественных биологических процессов, протекающих в организме нетели в соответствии с её физиологическим состоянием. Для второго звена масса вымени, формируемая за счёт механического воздействия на молочную железу (массажа), наряду с естественной активизацией тканеобразующих процессов сообразно периоду стельности животного.
Очевидно, что деятельность как каждого звена, так и подсистемы в целом должна оцениваться своей целевой функцией:
У1(та) - целевая функция звена выращивания нетелей до 6 мес. стельности;
/ С)
/2(тр) - целевая функция звена выращивания нетелей с 6 мес . стельности до отёла,
/Управляющее воздействие^
(1)
Ввод осеменённой тёлочки в нетельную группу
Контроль параметров вымени
Нетель в возрасте с 1 по 6 мес. стельности
Контроль параметров вымени
I—
Нетель в возрасте с 6 по 8 мес. стельности
Случайное воздействие
Направленное воздействие
Контроль параметров вымени
Л
Управляющее воздействие >
СУБСТРАТ
где /({) — целевая функция управления подсистемой выращивания нетелей.
Конечная цель функционирования подсистемы — получение максимальной массы вымени нетели. Известно, что математической моделью данной подсистемы является интегратор [1, 4], который с учётом смены режима работы в момент времени можно представить в виде:
т(г, т)
т1(г),г < т
у20) (т )Чг (г, т) +\' 42 (г, и) Х2 (и, т )йп, г > т,
V т
(2)
где ) — текущее значение массы вымени в период времени ;
ш^) — конечное значение первого периода выращивания нетели;
^(0)( ) — исходное значение выходной переменной второго периода выращивания нетели;
) и #2(?,и) — функции роста вымени в период времени и и,и 6( ,?); — х2(и, ) — воздействие.
Так как технологический процесс формирования вымени непрерывен, то:
у(20)( ) = ш().
(3)
Отсюда
т(г, т)<
(4)
Рис. - Схема подсистемы выращивания нетели
| т1 (г), г < т
Ц(т) + У2(г -т), г >т,
где ш1( ) — выходная переменная при 1 ; у2(1 — ) — прирост вымени в результате внешнего воздействия, осуществляемого в интервале времени и 6 ( , {). Математическое ожидание вектора выходных переменных подсистемы выращивания нетели может быть представлено суммой двух слагаемых, соответствующих возможному режиму функционирования:
т(г, т) = ух + у2, (5)
где 71 и у — математическое ожидание выходных переменных первого и второго режимов. Опуская промежуточные выкладки [1, 5], мы можем привести уравнение, полностью определяющее положение в пространстве подсистемы выращивания нетели:
в к т(г, т) = X Рг (та ) тг + X Р? (Ик )т?, (6)
¿=1 г=1
где г 1 — множество внешних воздействий;
Р — вероятности принятия значений массы ш;
— множество моментов времени нахождения нетелей в первом режиме; ик — множество моментов времени нахождения нетелей во втором режиме; а также целевую функцию управления технологическим процессом:
Л(га+и ) = т(та ) + т? (Ик ). (7)
Получение максимальной массы (объёма) вымени — главная задача выращивания нетели. Поэтому условие оптимального управления технологическим процессом будет иметь вид:
max ЛОа+н) = max m(Ta) + max m"(uk)• (8)
Анализ уравнения (8) показывает, что максимум целевой функции f (t + и) может быть достигнут при максимуме функции первого режима, а также при выборе оптимальных управлений, обеспечивающих максимум прироста вымени, во втором режиме.
Известны многие математические выражения, характеризующие динамику развития биологических объектов в результате преобразования энергии [6]. Если эта система находится в равновесии, т.е. не получает из внешней среды и, наоборот, ничего не теряет, то для неё будет справедливо выражение:
dm dS
-=--• (9)
dt dt
Тогда
dm dS d .
-+ — = — (m + S) = 0.
dt dt dt
Из уравнения (10) следует, что
(10)
dm = kmS, dt
(12)
где к — коэффициент пропорциональности, который в свою очередь определяется как:
k =
m
(13)
f
где ^ — удельный темп роста вымени нетели.
Так как из (11) следует, что при Sf = 0 (условие завершения преобразования энергии среды обитания):
S = Const = mf - m,
(14)
то уравнение (12) можно представить в виде:
dm dt
а с учётом (13): dm д dt
■ km(mf - m),
=-m(mf - m) = ]Xm
m
f
m m
f
В таком случае:
> m0
mf - m
v J
1
+ — m
\
-dm =1 ttdt.
0
(15)
(16)
(17)
/
Проинтегрировав (17), получим уравнение:
m =
m0mf-er
mf + m0e -1)
(18)
m + S = Const = m0 + S0 = mf + Sf = C, (11)
где m0 и S0 — исходные значения массы вымени m и энергии S источника преобразования в момент времени t = 0;
mf и Sf — значения, к которым приближаются значения m и S при t .
В таком случае наиболее приемлемым для описания характера преобразования энергии питательной среды в конечный продукт будет уравнение логистического роста, при использовании которого делается двоякое допущение: энергия роста пропорциональна массе вымени m; механизм роста работает со скоростью, пропорциональной ресурсу источника энергии S; процесс роста необратим. Тогда:
из которого следует, что рост вымени нетели осуществляется по экспоненте и носит асимптотический характер с замедлением темпа роста при достижении mf.
Исходя из условия формирования вымени нетели путём проведения массажа во второй половине стельности представляется весьма важным установление времени, с которого необходимо начинать применять эту технологическую операцию. Наиболее вероятный момент — точка перегиба кривой роста вымени. Определить её можно, продифференцировав уравнение (16) и приравняв его правую часть нулю. Из его решения следует, что перегиб кривой роста имеет место при достижении массы вымени:
1
m = — m
■f>
(19)
в момент времени 1мас, который можно определить путём подстановки (19) в уравнение (18):
= iln
m.
■m
m
(20)
Рост вымени нетели под воздействием массажа, реализуемого во второй половине стельности, достаточно подробно рассмотрен в работе В.Ф. Ужика [6]. Полученное им уравнение роста молочной железы имеет вид:
m,.
m + -
Pn
10
r
(1 - e-").
(21)
где m — суммарная масса вымени, кг;
т — масса вымени, определяемая полученным нами уравнением роста молочной железы (18), кг;
Р — мощность процесса массажа вымени нетели, Вт;
п10 — удельный прирост массы вымени нетели на единицу сообщаемой мощности в процессе массажа в период времени I =0, кг/Вт;
t
mac
г — коэффициент, учитывающий изменение влияния на прирост вымени в процессе массажа по мере старения организма нетели.
После подстановки (18) в (21) последнее приобретет вид:
+
ms =
m0 m^er
т^ + m0
+ -
Pn
10
(ev-1) (1 - e-")
(22)
X
a^L
aAL
1 W Amax 2 a
m1 + m2 2 X 2(m1 + m2)
Sß\/ Amax a 2
ax,y m1 + m2 2 4(m1 + m2)
Sß\/ Amax m2 a2 4m22 J
(23)
где т1 — масса вымени, кг;
т2 — масса массажного устройства, кг; а — коэффициент пропорциональности, кг/с;
£ — площадь поперечного сечения выме-
2
ни, м ;
— коэффициент пропорциональности деформации вымени, Н/м7/2; Атах — амплитуда колебаний системы вымя — массажное устройство, м. Подставив данное выражение в уравнение (22):
т„
m0m^e
т^ + m0
vt
/ -x+(l-e"rt)
(ev-1) 1 '
10 nr
m
1 + m2
2 (m1 + m2 ]
X
X
Таким образом, нами получено уточнённое уравнение, характеризующее динамику увеличения массы вымени нетели при наличии возмущающего воздействия на него.
После отёла из отелившихся нетелей формируют группы раздоя, где их оценивают по динамике роста молочной продуктивности в сравнении со средним показателем по предприятию. Не соответствующих критериям оценки коров и первотёлок по завершении периода раздоя выбраковывают.
Одним из механических приспособлений, обеспечивающих пневмомеханическое воздействие на рецепторные зоны молочной железы, может быть предложенное нами устройство [7]. Принцип его действия основан на вызове колебаний системы массажное устройство — вымя под воздействием пневмовибратора.
Мощность Рмас процесса массажа вымени нетели, развиваемая предлагаемым устройством, равна:
aALa
SßV Amax a 2
m1 + m2 4 (m1 + m2 )2
Ы Amax a2
m2 4m2
(24)
мы получим математическую модель роста вымени нетели при внешнем воздействии на молочную железу путём массажа с учётом конструктивно-режимных параметров массажного устройства.
Экспериментальными исследованиями массажного устройства установлено, что амплитуда колебаний вымени под его воздействием зависит от массы устройства, массы груза пневмови-братора, частоты колебаний, а также диаметра мембраны. Так, для обеспечения амплитуды колебаний вымени, равной 18 мм, при массе молочной железы 3 кг, устройство для массажа вымени должно обладать следующими параметрами: масса 0,5 кг, масса груза пневмовибратора
0.13.кг, частота пульсаций — 1,7 Гц, диаметр мембраны — 0,12 м.
Обработкой результатов производственных испытаний массажного устройства подтверждена адекватность математических и эмпирических моделей роста вымени нетели при наличии возмущающего воздействия — массажа молочной железы. Установлено, что первотёлки, которым массаж вымени в нетельный период осуществляли экспериментальным устройством, превосходили своих аналогов в группе с массажем устройством АПМ-Ф-1 на 78 кг, а животных контрольной группы — на 156,7 кг молока. Это указывает на целесообразность его использования в производстве при формировании вымени высокопродуктивных коров.
Литература
1. Ужик В.Ф., Ужик О.В., "Укик Я.В. Теория технологий и технических средств в животноводстве: монография. Белгород: Изд-во БелГСХА, 2009. 198 с.
2. Скляревич A. E, Скляревич Ф. К. Вероятностные модели объектов с возможными изменениями. Рига: Зинатне, 1989. 366 с.
3. Поляк Б.Т. Введение в оптимизацию. М.: Наука, 1983. 384 с.
4. Вентцель Е.С., Овчаров Л.А. Теория вероятностей и её инженерные приложения. М.: Наука, 1988. 480 с.
5. Ужик В.Ф. Совершенствование средств механизации для формирования вымени высокопродуктивных коров: автореф. дисс.... докт. техн. наук. Оренбург, 1994. 38 с.
6. Франс Дж., Торнли Дж.Х.М. Математические модели в сельском хозяйстве / Пер. с англ. А.С. Каменского; под. ред. Ф.И. Ерешко. Предисл. Ф.И. Ерешко, А.С. Каменского. М.: Агропромиздат, 1987. 400 с.
7. Патент №2368133 RU, С1 МПК А 01 J 7/00 Устройство для массажа вымени нетелей / Ужик О.В. (RU). №2008140357/12; Заявлено 10.10.2008; Опубл. 27.09.2009, Бюл. №27.
2