Issue 
BIO Web Conf.
Volume 10, 2018
Contemporary Research Trends in Agricultural Engineering



Article Number  02004  
Number of page(s)  7  
Section  Engineering and Technology  
DOI  https://doi.org/10.1051/bioconf/20181002004  
Published online  26 March 2018 
Models of adaptation of the milking machines systems
Lviv National Agrarian University, Faculty of Mechanic and Power Engineering, LvivDubliany, Ukraine
^{*} Corresponding author: Dmytriv_V@ukr.net
Two systems of milking machines were considered  biotechnical and vacuum. Methodology for estimation of efficiency of the "machineanimal" biotechnical system was worked out. The dependences of the efficiency parameters of the technical system operation were analyzed. The K_{MMO} operator load factor of the milking machine was proposed. The factor characterizes the technological process of the machine milking. The analytical dependences were worked out for the simulation of the productivity of milking machines and oscillation of the vacuumgage pressure. As to the simulation results, when the vacuum pipeline diameter was increased the oscillation of the vacuumgage pressure decreased. The results of analysis and theoretical researches on technological process of the cow machine milking gave a possibility to define the requirements to the improvement of technological process and technical equipment, which will provide the increase of efficiency of the milking systems. Usage of the developed cyberphysical system of the machine milking of cows will increase the productivity of the milking machine in 1.26…1.85 times. At the vacuum gage pressure oscillation of ΔP_{vp} = 2500 Pa the suction ability of milking machine will be E=4.093 m/s accordingly. The defined index of efficiency of the adapted systems functioning of the milking machine is K_{BTS2} = 55.3.
Key words: milking machine / mathematical model / vacuum / amplitude / load factor.
© The Authors, published by EDP Sciences, 2018
This is an Open Access article distributed under the terms of the Creative Commons Attribution License 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. (http://creativecommons.org/licenses/by/4.0/).
1 Introduction
We will consider the process of cows milking as a biotechnical system of the "personmachineanimal", that does possible the functional ensuring of realization of the genetic potential of the cow productivity with the interaction of objects of this system. Efficiency of functioning of the biotechnical system depends on the parameters of its objects that ensure the quality and efficiency of technological functions.
As a basic performer of the technological process of milking, the milking machine provides adaptation of the technical system to physiology of milk ejection of the cow through vacuum gage pressure with preset parameter and possibility of its adjusting. The milking machine provides also the parameters control and maintenance at the set level during milking of cows. Pulsating air system of milking machine is the peculiarity of unit arrangement that changes the configuration and functional possibilities. Pulsating air system gives an opportunity also to work autonomically, or in composition of the automatic control system of technological process (Industrial Control – IC) of the machine milking. To operate in the composition of IC the milking machine was additionally completed by the measuring device of milk ejection intensity, microprocessor control unit and other elements of electronics [1, 2]. However, the vacuum gage pressure forming is still one of important factors of the parameters forming directly that does possible adaptation of the “machinecow” system. The question of parameters adaptation of the milking machine systems was examined by Dmytriv V.T. and other [35], where influence of the technological and construction parameters on work mode of systems is analyzed. Also, the stability of vacuum gage pressure was estimated directly by the group of researchers consisting of Pazzona A., Murgia L., Zanini L. et al., depending on the method of adjusting in the vacuum pipeline or in a milking machine. They set that stabilizing of vacuum by the regulator of gravitational type is more dynamic, and the permanent to time is twice more lower from the computerbased system [6]. An analogical conclusion is also in accordance with research results of the regulators of pulsating vacuum gage pressure [7]. In particular Reinemann D.J., Schuring N. and Bade1 R.D. investigated the vacuum gage pressure oscillation depending on: a) configuration of the vacuum and milk hose systems (length of pipelines, hose diameter and other parameters influencing the losses of pressure); b) flowrates of milk in the milk pipeline; c) rates of air movement in the vacuum pipeline. They set that the vacuum gage pressure was raised up as in a vacuum and milk pipeline, so in the underteat chamber of the teat cup in the process of diminishing of the intensity of milk ejection [8].
The analysis of research works shows that the grounding of parameters of the effective functioning of the cows milking bioengineering system requires the development of conception and methodology of optimization of the mentioned system parameters and the same for the technical constituent of this system. Optimization of the bioengineering system parameters is also important for the increase of efficiency of milking machines.
2 Materials and methods
2.1 Effectiveness of the "operatormachine animal" bioengineering system
The effectiveness of work of the machine milking operator (MMO) was determined by the productivity for the maximum system load. This parameter depends on the n_{mu} quantity of the milking units and the t_{m} duration of the milking of one cow (with consideration of the preparation and final operations t_{m} = t_{mm} + t_{pf}), the Z_{li} labour inputs:
$${W}_{mm}=f\left({n}_{mu},{t}_{m},{Z}_{li}\right),$$(1)where: W_{mm} – the productivity of the milking machine, cows/h., n_{mu} – quantity of the milking units, pcs., t_{m} – duration of the milking of one cow, h., Z_{li} – labour inputs, MJ.
Effectiveness of functioning of the technical system of milking was estimated by the quality parameter of functioning of the vacuum system, namely the Δp_{VP} amplitude of pressure oscillation of the vacuum system of milking machine. The vacuum gage pressure oscillation depends on the П_{es} parameters of the mentioned technical system elements, the Kf_{es} quality of its interaction and the Z_{TS} energy expense for the system functioning:
$$\Delta {p}_{VP}=f\left({\Pi}_{es},{K}_{fes},{Z}_{TS}\right),$$(2)where Δp_{VP} – the amplitude of vacuum gage pressure oscillation, Pa, Π_{es} – parameters of technical system elements, Kf_{es} – the quality coefficient of the interaction of technical system elements, Z_{TS} – the energy expense for the technical system functioning, MJ.
Effectiveness of the bioengineering system depends on the cow milk ejection that is limited by the pressure oscillation in the vacuum system of the milking machine. The milk yield decrease can be described by the function as a result of the vacuum gage pressure oscillation [9]:
$$\Delta {Q}_{M}={Q}_{P}\cdot \left[\frac{\epsilon}{\epsilon +1}\cdot {\left1\frac{2\cdot {E}_{1}}{3\cdot J}\right}^{3/2}+\frac{1}{\epsilon +1}\cdot \left({\left1\frac{2\cdot {E}_{2}}{3\cdot J}\right}^{3/2}{\left1\frac{2\cdot E}{3\cdot J}\right}^{3/2}\right)\right]$$(3)where ΔQ_{M} – the loss of cow milk yield as a result of the pressure oscillation, MJ, Q_{p} – the potential of the cow milk yield, MJ, ε – the ratio of the pressure oscillation duration to the total milking duration, E_{1},E_{2}  the suction ability of the milking machine with the modes accordingly p_{1}=p+Δp_{VP}/2 i p_{1}=p Δp_{VP}/2, p – the vacuum gage pressure in the under teat chamber of the teat cup, Pa, E – the suction ability of the milking machine for the vacuum gage pressure in the under teat chamber of the teat cup, J=E·2/3 – parameter of the milk ejection.
where K_{f} – the proportion of phases of the milking machine, K_{ε} – flow coefficient, ρ_{M} – milk density, kg/M^{3}.
The suction ability of the milking machine was determined according to the described procedure [9] :
$$E={K}_{f}\cdot {K}_{\epsilon}\cdot {\left(2\cdot p/{\rho}_{M}\right)}^{1/2}/\left({K}_{f}+1\right),$$(4)On the basis of the above mentioned parameter the effectiveness of the milking machine system is expressed by the formula:
$${K}_{BTS}=\frac{1}{{Z}_{BTS}}\sum _{n=1}^{i}{\left({Q}_{p}\Delta {Q}_{M}\right)}^{,}$$(5)where Z_{BTS} – the energy expense for the functioning of the milking machine system, Z_{BTS} = Z_{li} + Z_{TS}, MJ.
The demonstrated formula (5) presents that effectiveness depends of the parameters of the technical system. The mentioned technical system has to ensure the conditions of the milk ejection from the cow udder according to the lactogenesis physiology demand with the reduction of the MMO labour inputs and also total energy expense lowering for the functioning of the milking machine systems.
Taking into account the order of the preparation and final operations of machine milking one can determine the K_{on} load factor of the milking machine operator by the formula [10]:
$${K}_{MNO}=\frac{{t}_{m}+\left({n}_{mu}1\right)\cdot {t}_{f},}{{n}_{mu}\cdot {t}_{pf}}$$(6)where t_{m} – duration of the one cow milking, t_{m} t_{mm}+t_{pf}, s [11], t_{f}− duration of the final operations, s.
If the KMMO < 1 – operator is underloaded, accordingly the milking machine is not overdone on the udder teats, if the K_{MMO} > 1 – the machine milking operator is overloaded and the schedule of preparation and final operations is not kept, or the milking machine is overdone on the udder teats of cow.
From the equation (6) we shall derive the formula for determining the optimal quantity of the milking machines, taking into account that t_{pf} = t_{p}+ t_{f},
$${n}_{mu}={\left[{K}_{MMO}+\frac{{t}_{f}}{{t}_{p}}({K}_{MMO}1)\right]}^{1}\cdot \left[\frac{{t}_{f}}{{t}_{p}}+1\right]$$(7)where t_{p} – the duration of the preparation operations by the MMO, s.
Taking into account the average distribution of the duration of preparation and final operations $\left({t}_{w}={K}_{MMO}\cdot {n}_{mu}\cdot \overline{t}\pm S\left({t}_{f}\right)\right)$ the equation of the productivity of milking machine will have the appearance:
$${W}_{mm}=\frac{3600.{n}_{mu}}{{K}_{MMO}\cdot nmu\cdot},$$(8)where ${\overline{t}}_{p},{\overline{t}}_{f}$ – the average value of the duration of duration preparation and final operations of the cow machine milking accordingly, s, S(t_{p}), S(t_{f}) – the rootmeansquare deviation of the duration of preparation and final operations of the cow machine milking accordingly, s.
2.2 Amplitude of the pressure oscillation of vacuum system of the milking machine
The construction of the milking machine can be presented as the functionally finished elements, which are the constituents of the vacuum system of the milking machine (Fig. 1) [12].
Taking into account the law of mass conservation for the gas in the controlled volume by way of equation of the air motion mechanical energy and work for overcoming of the frictional force, the second airflow during the pumpingout of the vacuum system will be [13]:
$$m={S}_{csa}\cdot \sqrt{\frac{2\cdot g}{1+\xi}\cdot {P}_{1}\cdot {\rho}_{1}\cdot \left({x}^{\frac{2}{n}}{x}^{\frac{n+1}{n}}\right).}$$(9)Differential equation of the V air volume pumping out at the condition of τ air mass by the τ time equals to dM = V·dρ the mass increase [14]:
$$dM={S}_{csa}\cdot \psi \cdot \sqrt{{p}_{i}\cdot}{\rho}_{i}\cdot d\tau ,$$(10)Where $\psi =\sqrt{\frac{2\cdot g}{1+\xi}\cdot \left({x}^{\frac{2}{n}}{x}^{\frac{n+1}{n}}\right)}$ – velocity aspect ratio defines the pressure relation at the air pumpingout, x – ratio of the current value of the vacuum gage pressure to the atmospheric pressure, x = p_{i}p_{A}, p_{i}i ρ_{i}  current value of the pressure and the density in the volume V.
Taking into account that ${\rho}_{i}={\rho}_{A}\cdot {\left(\frac{{p}_{i}}{{p}_{A}}\right)}^{\frac{1}{n},}$ the equation
(10) will have the appearance:
$$\begin{array}{l}dM={S}_{csa}\cdot \psi \cdot \sqrt{{p}_{i}\cdot {\rho}_{A}\cdot {\left(\frac{{p}_{i}}{{p}_{A}}\right)}^{\frac{1}{n}}\cdot d\tau =}\\ ={S}_{csa}\cdot \psi \cdot \sqrt{{p}_{A}\cdot {\rho}_{A}\cdot {\left(\frac{{p}_{i}}{{p}_{A}}\right)}^{\frac{n+1}{n}}\cdot d\tau}\end{array}$$(11)Also taking into account that dM = Vdρ_{i}, we will have:
$$dM=V\cdot {\rho}_{A}\cdot d\left[{\left(\frac{{p}_{i}}{{p}_{A}}\right)}^{\frac{1}{n}}\right]=\frac{V}{n}\cdot {\rho}_{A}\cdot {\left(\frac{{p}_{i}}{{p}_{A}}\right)}^{\frac{1n}{n}}\cdot d{\left(\frac{{p}_{i}}{{p}_{A}}\right)}^{.}$$(12)After the dependences of (11) and (12) were equated, the sign was accounted and formula was reduced by(p_{i}/p_{A})^{1/2}, the differential equation of air pumping out from the V volume will have the appearance:
$$\begin{array}{l}{\rho}_{A}\cdot \frac{V}{n}\cdot {\left(\frac{{p}_{i}}{{p}_{A}}\right)}^{1}\cdot d\left(\frac{{p}_{i}}{{p}_{A}}\right)=\\ ={\mathrm{S}}_{csa}\cdot \psi \cdot \sqrt{{p}_{A}\cdot {\rho}_{A}}\cdot \sqrt{{\left(\frac{{p}_{i}}{{p}_{A}}\right)}^{\frac{n1}{n}}\cdot d\tau}\end{array}$$(13)The (13) formula is transformed and integrated with the limitations from 1 to p_{i}/p_{A} and from 0 to τ:
$$\frac{1}{n}\cdot {\displaystyle \underset{1}{\overset{\frac{{p}_{i}}{{p}_{A}}}{\int}}{\left(\frac{{p}_{i}}{{p}_{A}}\right)}^{\frac{13n}{2n}}}\cdot d\left(\frac{{p}_{i}}{{p}_{A}}\right)={a}_{csa}\cdot \psi \cdot \frac{1}{V}\cdot \sqrt{\frac{{p}_{A}}{{\rho}_{A}}}\cdot {\displaystyle \underset{0}{\overset{t}{\int}}d\tau}\cdot $$(14)The integration results define the duration of air pumping out of the V volume from the pressure p_{A} to the vacuum gage pressure p_{i}, and will be determined in series by the formula:
$$\frac{1}{n}\cdot \frac{{\left(\frac{{p}_{i}}{{p}_{A}}\right)}^{\frac{1n}{2n}}{}_{1}^{\frac{{p}_{i}}{{p}_{A}}}}{\frac{1n}{2n}}={S}_{csa}\cdot \psi \cdot \frac{1}{V}\cdot \sqrt{\frac{{p}_{A}}{{\rho}_{A}}\cdot \tau {}_{0}^{\tau}},$$(15)or
$$\tau =\frac{2}{1n}\cdot \frac{V}{{s}_{csa}\cdot \psi}\cdot \sqrt{\frac{{\tau}_{A}}{{p}_{A}}\cdot \sqrt{{\left(\frac{{p}_{i}}{{p}_{A}}\right)}^{\frac{1n}{n}}1}}.$$(16)To simulate the characteristic of the vacuum gage pressure changes in the vacuum system of the milking machine one has to know the time parameters of the vacuum pressure  the time interval of the air pumping out from the volumes of the variable vacuum gage pressure. From the equation (16) the dependences of pressure changes in time have been derived:
$${p}_{i}={p}_{A}\cdot 1n\sqrt{\frac{{\left(\tau +1\right)}^{2\cdot n}.{\left(1n\right)}^{2\cdot n}\cdot {S}_{csa}^{2\cdot n}\cdot {\psi}^{2\cdot n}\cdot {p}_{A}}{{4}^{n}\cdot {V}^{2\cdot n}\cdot {\rho}^{n}}},$$(17)where v – the volume of chambers of the variable vacuum gage pressure of the milking machine, M^{3}, S_{csa}  the crosssectional area of the pulsator drain port, M^{2}, ψ – velocity aspect ratio defines the pressure relation, m^{1/2}/s, ρ_{A} – the air density in the V volume for the atmospheric pressure, kg/M^{3}, p_{A} – the air atmospheric pressure in the V volume, kg/M^{2}, p_{i} – the value of the vacuum gage pressure in the ith point of time, kg/M^{2}, n – polytropic coefficient, n = 1,41.
Thus, the quantity of pumped out air will be:
$${G}_{i}=\frac{{p}_{i}\cdot V,}{{R}_{n}\cdot \Theta}$$(18)where V – the volume of chambers of the variable pressure, M^{3}, R_{n} – gas constant of the air, Rn = 287,72 J/(kg·°C), Θ – air temperature, °C.
For connection of the i chambers with the given volume one can determine the volume, pressure and the air quantity accordingly [12]:
$$V={\displaystyle \sum _{i=1}^{n}{V}_{i}{}^{,}PV=}\frac{1}{n}{{\displaystyle \sum _{i=1}^{n}{p}_{i}}}^{\text{'}}G={\displaystyle \sum _{\tau =0}^{1}{G}_{i}}$$(19)The amplitude of the vacuum gage pressure oscillation in the vacuum pipeline is determined by the formula:
$$\Delta {p}_{VP}=\frac{{R}_{n}\cdot \Theta \cdot \left({G}_{VP}+G+G\cdot \left({N}_{nmu}^{+}\sum 1\right)\cdot P\left({N}_{nmu}\right)\right),}{{V}_{VP}}$$(20)where G_{Vp}, G – the air quantity in vacuum pipeline and volumes of the variable vacuum gage pressure accordingly, kg, V_{Vp} – the volume of the vacuum pipeline of the specific milking machine, M^{3}, P(N_{nmu}) – the probability of simultaneity and coincidence of the phase and work times of the specific milking machines:
$$P\left({N}_{num}\right)=P\left({N}_{nmu}^{+}\right)\cdot P\left({N}_{nmu}^{}\right),$$(21)Where $P\left({N}_{nmu}^{+}\right)=\frac{{a}_{1}^{{N}_{nmu}}}{{N}_{nmu}^{+}\downarrow}\cdot {\text{e}}^{a1}$ – the probability of simultaneity and coincidence of the phase and work times of the milking machines according to the Poisson limit theorem [15] $P\left({N}_{nmu}^{}\right)=\frac{{a}_{2}^{{N}_{nmu}}}{{N}_{nmu}^{}\downarrow}\cdot {\text{e}}^{a2}$ – the probability of unsimultaneity and uncoincidence of the phase and work times of the milking machines according to the Poisson limit theorem [15],${a}_{1}={n}_{nmu}\cdot p\left(\right),{a}_{2}={n}_{nmu}\cdot p\left(\right),{n}_{mu}$ – the quantity of the simultaneously working milking machines, pcs., N_{NMU}^{+} – the quantity of the simultaneously working milking machines with the coincidence in time of the phase and work times of pulsator, pcs., N_{NMU}^{−} – the quantity of the simultaneously working milking machines with the uncoincidence in time of the phase and work times of pulsator, pcs., p(^{+}_{nmu}) – probability of the coincidence in time of the phases and work times of pulsator for N^{+}_{nmu} – simultaneously working milking machines, p(^{+}_{nmu}) = N^{+}_{nmu}/n_{mnu}, p(^{−}_{nmu}) – probability of the uncoincidence in time of the phases and work times of pulsator for N^{−}_{nmu} simultaneously working milking machines, p(^{−}_{nmu}) = N^{−}_{nmu}/n_{mnu}.
Fig. 1. The structure and functionally scheme of the vacuum system of milking machine: 1 – vacuum pump; 2 – vacuum reservoir; 3 – vacuum hose; 4 – vacuum gage pressure regulator; 5 – pulsator of the milking machine; 5.1 – chamber of the constant vacuum gage pressure; 5.2, 5.3 – chambers of the variable pressure; 5.4 – chamber of the atmospheric pressure; 6 – controller of the milking machine; 7 – interwall chamber of the teet cup; 8 – milk pipeline. 
Results
The productivity of the milking machine (equation 8) depends on the duration of preparation t_{p} and final t_{f} operations, duration of machine milking tmm and the cow productivity. Also the productivity of the milking machine depends of the quantity of the milking units, which are operated by the MMO with the regulated load factor K_{MMO}. The results of the modelling are shown in Fig. 2 and 3.
The analysis of dependence of the productivity of the milking machine from the time tmm of machine milking of cow and ratio of t_{f} /t_{p} (Fig.2) demonstrates when t_{f} → max, t_{p} → max, the productivity of the MMO work is submitted to the specific law. The area on the 3D graph of the 26…34 cows/hour productivity fits with the milking machine for the bucket milking, 34…50 cows/hour – pipeline milking in the barn, from 51 cows/hour – parlor milking
The analysis of dependence of the productivity of the milking machine from the time tmm of machine milking of cow and ratio of t_{f} /t_{p} (Fig.2) demonstrates when t_{f} → min, t_{p} → min, the productivity of the MMO work is submitted to the secondorder equation.
The area on the 3D graph of the 29…36 cows/hour productivity fits with the milking machine for the bucket milking, 37…50 cows/hour with the pipeline milking in the barn, from the 51 cows/hour and high with the parlor milking.
Modelling of the pressure oscillation amplitude in the vacuum pipeline is demonstrated in Fig. 4.
Analysis of results of the vacuum pressure oscillations modelling shows that the maximum oscillations amplitude of vacuum pressure at a frequency of 1 Hz does not exceed Δp_{vp} = 1.4 kPa for vacuum pipeline of D = 40 mm diameter with a probability of P(N_{nmu}) = 0.76212. To vacuum pipeline with diameter of D = 25.4 mm the maximal oscillations amplitude of vacuum pressure is Δp_{vp} = 3.03 kPa with probability of P(N_{nmu}) = 0.628837.
For the number of milking machines which are working synchronous with the maximum probability of coincidence in time cycles and phase of pulsators work the vacuum pressure oscillations does not exceed Δp_{vp} = 2 kPa. For the vacuum pipeline of D = 50 mm diameter the vacuum pressure oscillations are Δp_{vp} = 800…812 Pa with the maximum probability of P(N_{nmu}) = 0.93757…0.95773.
Graphic dependence of the vacuum pressure oscillation amplitude changes Δp_{vp} from the number of simultaneously working milking machines N_{mnu}^{+} with coincidence in time of the cycles and phases of pulsator work with probability P(N_{nmu}) and the vacuum pipeline D diameter was shown in Fig. 5.
The analysis of simulation results shows that with an increase of the number of milking machines the vacuum pressure amplitude is increased. From N_{mnu}^{+} with the n_{mu} = 8 probability of coincidence in time of the pulsator cycles and phases of P(N_{nmu}) = 0.592547 the oscillations amplitude of vacuum pressure will be Δp_{vp} = 4.72 kP. for the diameter of milking machines vacuum pipeline of D = 25.4 mm. For the vacuum pipeline diameter of D = 50 mm and the same number of milking machines the amplitude of vacuum pressure oscillations will be Δp_{vp} = 1.08 kP..
The maximum probability of coincidence in time of the pulsators cycles and phases is P(N_{nmu}) = 0.95773 for N_{mnu}^{+} = 4 machines, then the amplitude of vacuum pressure oscillations will be equal to Δp_{vp} = 3.55 kP. per one cycle of pulsator work (frequency of 1.0 Hz) for milking machines vacuum pipeline diameter of D = 25.4 mm and Δp_{vp} = 0.81 kP. for D = 50 mm.
Fig. 2. Dependence of the productivity of the milking machine from the duration t_{mm} of machine milking of cow and ratio of t_{f}/t_{p}, by the t_{f} → max, t_{p}→ max. 
Fig. 3. Dependence of the productivity of the milking machine from the time t_{mm} of machine milking of cow and ratio of t_{f}/t_{p}, by the t_{f} → max, t_{p}→ max. 
Fig. 4. Graphs of modelling of the vacuum gage pressure oscillation Δp_{vp} in the vacuum pipeline of milking machine in dependence of quantity of the milking units N ^{+}_{nmu} with the coincidence in time of the phase and work times of pulsator with probability of P(N_{nmu}) for diameter of vacuum pipeline D: a – D = 25,4 mm, milking unit with pneumatic pulsator; b – D = 40 mm, milking unit with pneumatic and electric pulsator of the pairwise action; c, d – D = 50 mm, milking unit with pneumatic and electric pulsator of the pairwise action. 
Fig. 5. The dependence of change of the Δp_{vp} vacuum pressure amplitude oscillation from the number of simultaneous working milking machines N_{mnu}^{+} with coincident in time of pulsator cycles and phases and D machine vacuum pipeline diameter. 
3 Discussion and conclusion
The dependences (15) show that efficiency depends on the parameters of the technical system which have to provide the conditions for the milk take out from the cow udder in accordance with the physiology of milk lactogenesis and ejection with reduction of both the MMO labor inputs and energy expense in the process of milking machine system functioning.
The requirements of the machine milking process were fulfilled for the K_{MMO} ≥ 1. The milking duration in the ranging ratio of t_{f} /t_{p} = 0.5…0.8 are reduced due to the optimal number of milking machines which are operated by MMO with the K_{MMO} ≈ 1. At the conditions of KMMO > 1 the MMO has no enough time to fulfil all operations by the requirements of the machine milking process due to unreasonable number of operated milking machines.
It can be noted that the maximum vacuum oscillation in milking machines vacuum pipeline is not more for the permissible parameters (Δp_{vp} ≤ 2.5 kPa) at the milking in the cattle barn with a frequency of 1 Hz. With the diameter of vacuum pipeline increasing the amplitude of vacuum gage pressure oscillation is decreased.
The analysis and the theoretical studies of technological process of the cow machine milking make possible to define the requirements for the improvement of technological process and technical means with the efficiency rise of milking systems in the following way:
 
to provide a flexible relationship amongst the intensity of milk ejection and constructive and technological parameters of milking machine systems, to provide the adaptation of the technical system into the function of milk ejection of each individual cow;
 
to create of cyberphysical milking machine that would conform with the physiological characteristics of cows as to the milk ejection;
 
to optimize of the organizational and technological and structural energetic parameters of technological process and milking machines for adaptation both the technical system and the physiological state of the cow organism.
Using of cyberphysical systems of cows machine milking increases the milking machines productivity by the 1.26…1.85 times.
For the vacuum pressure oscillations of Δp_{vp} = 2500 Pa the suction ability of milking machines for modes will be accordingly of E = 4.093, E_{1} = 4.177, E_{2} = 4.007 m/s. For the milking machine phase proportions were taken the standard value. Accordingly the reducing of milk yields per cow in the energy equivalent is 65.1 MJ (3.08 MJ/kg). The energy equivalent of milk production will be 11450 MJ accordingly. The energy expenses of milking machine for the year in energetic equivalent will be 16556.4 MJ. The efficiency index of conventional milking machine will be K_{BTSS1} = 46.6. According to the results of studies of the technological process of milking, the 450 W (1.62 MJ) reducing of the vacuum machine drive power will enable annual saving of 1773.9 MJ. Accordingly in practice one can use the 1 MMO less with saving in the energy equivalent of 1.971 MJ. Then the efficiency index of adapted systems functioning of the milking machines will be K_{BTS2} = 55.3. Accordingly, the efficiency index of the milking machines systems functioning will increase by 18.7%.
References
 V. T. Dmytriv, New direction of development of the technologies and technical means in dairy cattle husbandry: Materials of the 13th International Symposium on Machine Milking of Farm Animals, 2729 of June, 2006. Milking machine with microprocessor control, 161165 (2006). [Google Scholar]
 V. T. Dmytriv, Mechanization and Electrification of Agriculture, Interdepartmental special scientific coll. Hlevakha, Automated milking machine, 84, 134136 (2010) [Google Scholar]
 V. T. Dmytriv, Mechanization in Agriculture. International Scientific, Scientific Applied and Informational Journal, Sofia, Adaptive machine milking system, 10, 1518 (2015) [Google Scholar]
 V. T. Dmytriv, Mechanics and Technology Bases of the Milking Machine Systems. Theory and Practice: Monograph (2017) [Google Scholar]
 V. Dmytriv, V. Adamchuk, I. Dmytriv, MOTROL. Commission of Motorization and Energetics in Agriculture, Adaptive machine milking with the pneumatic electric magnetic pulsator, 17 (9), 8387 (2015) [Google Scholar]
 A. Pazzona, L. Murgia, L. Zanini, et all. ASAE Annual International Meeting, Riviera Hotel and Convention Center, Las Vegas, Nevada, USA, Dry tests of vacuum stability in milking machines with conventional regulators and adjustable speed vacuum pump controllers (2003) [Google Scholar]
 V. Dmytriv, MOTROL. Commission of Motorization and Energetics in Agriculture, Analytical study of the pressure regulator for pulsating vacuum systems, 18/8, 6976 (2016) [Google Scholar]
 D. J. Reinemann, N. Schuring, R. D. Badel, The Sixth International ASABE Dairy Housing Conference, Minneapolis, Minnesota, USA, Methods for Measuring and Interpreting of the Milking Vacuum (2007) [Google Scholar]
 L.P. Kartashov, Methods for calculating of the biological and technical parameters of the “humanmachineanimal” system, Orenburg: OS AU Publishing Center (2007) [Google Scholar]
 V. T. Dmytriv, Research Journal of National University of Life and Environmental Sciences of Ukraine: Technics and Energy of Agriindustrial Complex, Modeling of the number of milking machines at the optimum workload of the milking machine operator, 212/2, 167173 (2015) [Google Scholar]
 V. T. Dmytriv, Bases of the Theory of the Machinery Usage in the Animal Husbandry, Lviv : Afisha (2008) [Google Scholar]
 V. Dmytriv, Journal of Lviv National Agrarian University : Agrarian Engineer Research, Model of the air flow by the elements of the milking machine, 10, 483488 (2006) [Google Scholar]
 I. V. Dmytriv, V. T. Dmytriv, Research Journal of the Kharkiv Petro Vasylenko National Technical University of Agriculture ”Technical Systems and Technologies of the Animal Husbandry“, Model of air pumping from the chambers of variable vacuum gage pressure of the milking machine, 132, 207212 (2013) [Google Scholar]
 V. Dmytriv, I. Dmytriv, MOTROL. Commission of Motorization and Energetics in Agriculture, Modeling of the time of air outflow from the bounded space, 15/4, 193197 (2013) [Google Scholar]
 P. F. Ovchinnikov, B. M. Lisitsin, V. M. Mikhailenko, Higher Mathematics, Kyiv : Higher school, (1989) [Google Scholar]
All Figures
Fig. 1. The structure and functionally scheme of the vacuum system of milking machine: 1 – vacuum pump; 2 – vacuum reservoir; 3 – vacuum hose; 4 – vacuum gage pressure regulator; 5 – pulsator of the milking machine; 5.1 – chamber of the constant vacuum gage pressure; 5.2, 5.3 – chambers of the variable pressure; 5.4 – chamber of the atmospheric pressure; 6 – controller of the milking machine; 7 – interwall chamber of the teet cup; 8 – milk pipeline. 

In the text 
Fig. 2. Dependence of the productivity of the milking machine from the duration t_{mm} of machine milking of cow and ratio of t_{f}/t_{p}, by the t_{f} → max, t_{p}→ max. 

In the text 
Fig. 3. Dependence of the productivity of the milking machine from the time t_{mm} of machine milking of cow and ratio of t_{f}/t_{p}, by the t_{f} → max, t_{p}→ max. 

In the text 
Fig. 4. Graphs of modelling of the vacuum gage pressure oscillation Δp_{vp} in the vacuum pipeline of milking machine in dependence of quantity of the milking units N ^{+}_{nmu} with the coincidence in time of the phase and work times of pulsator with probability of P(N_{nmu}) for diameter of vacuum pipeline D: a – D = 25,4 mm, milking unit with pneumatic pulsator; b – D = 40 mm, milking unit with pneumatic and electric pulsator of the pairwise action; c, d – D = 50 mm, milking unit with pneumatic and electric pulsator of the pairwise action. 

In the text 
Fig. 5. The dependence of change of the Δp_{vp} vacuum pressure amplitude oscillation from the number of simultaneous working milking machines N_{mnu}^{+} with coincident in time of pulsator cycles and phases and D machine vacuum pipeline diameter. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.