№22/2018

Norwegian Journal of development of the International Science

ISSN 3453-9875

VOL.1

It was established in November 2016 with support from the Norwegian Academy of Science.

DESCRIPTION
The Scientific journal “Norwegian Journal of development of the International Science” is issued 12 times a year
and is a scientific publication on topical problems of science.

Editor in chief – Karin Kristiansen (University of Oslo, Norway)

The assistant of theeditor in chief – Olof Hansen

 James Smith (University of Birmingham, UK)

 Kristian Nilsen (University Centre in Svalbard, Norway)
 Arne Jensen (Norwegian University of Science and Technology, Norway)
 Sander Svein (University of Tromsø, Norway)
 Lena Meyer (University of Gothenburg, Sweden)
 Hans Rasmussen (University of Southern Denmark, Denmark)
 Chantal Girard (ESC Rennes School of Business, France)
 Ann Claes (University of Groningen, Netherlands)
 Ingrid Karlsen (University of Oslo, Norway)
 Terje Gruterson (Norwegian Institute of Public Health, Norway)
 Sander Langfjord (University Hospital, Norway)
 Fredrik Mardosas (Oslo and Akershus University College, Norway)
 Emil Berger (Ministry of Agriculture and Food, Norway)
 Sofie Olsen (BioFokus, Norway)
 Rolf Ulrich Becker (University of Duisburg-Essen, Germany)
 Lutz Jäncke (University of Zürich, Switzerland)
 Elizabeth Davies (University of Glasgow, UK)
 Chan Jiang(Peking University, China)
and other independent experts

1000 copies
Norwegian Journal of development of the International Science
Iduns gate 4A, 0178, Oslo, Norway
email: [email protected]
site: http://www.njd-iscience.com
 CONTENT
BIOLOGICAL SCIENCES
Loiko S., Romanova Z., Romanov N.
UNMALTED RAW MATERIALS IN
BREWING .................................................................. 3

CHEMICAL SCIENCES
Fedenko Yu.
REGULARITIES OF SORPTION REMOVING PB2+
BY BENTONITE CLAYS OF THE UKRAINIAN
ORIGIN ...................................................................... 8

JURISPRUDENCE
Novikov O.
CURRENT STATUS AND TRENDS OF
INTELLECTUAL PROPERTY CRIMES IN
UKRAINE................................................................. 12

TECHNICAL SCIENCES
Shibaev D., Dorofeev A. Romanov N., Svitoch I., Macebula D.
DEVELOPMENT OF A MOBILE APPLICATION MATHEMATICAL MODELS IDENTIFICATION
FOR TRACKING GEOLOCATION POSITION IN IN CHARACTERISTIC SITUATIONS OF
REAL TIME ............................................................. 19 CONDUCT OF THE BREWERY
Kaganov W. TECHNOLOGICAL COMPLEX ............................. 37
PROTECTION FROM THE ICEOF ELECTRICAL Fyk A.
NETWORKS BY MEANS OF A HIGH- MODELING OF THE PROCESS OF
FREQUENCY ELECTROMAGNETIC WAVE ...... 25 ELECTROMAGNETIC WAVE TRANSFER INTO
Polishchuk M. THE COAXIAL WAVE CABLE THROUGH THE
PRINCIPLES OF DESIGNING MOBILE HIGH-TEMPERATURE SUPERCONDUCTING
ROBOTS ................................................................... 31 ELEMENT FROM THE WAVEGUIDE .................. 44

VETERINARY SCIENCES
Zhukov M., Bosenko T.,
Morgunova V., Kashirina L., Tyurina E.
SABILITY OF BIOCHEMICAL PARAMETERS OF
FECES DURING STORAGE ................................... 55
Norwegian Journal of development of the International Science No 22/2018 3

BIOLOGICAL SCIENCES
UNMALTED RAW MATERIALS IN BREWING

Loiko S.
Technologist, Novocontact Ltd., Berezhanska str.,6a, 04074, Kyiv, Ukraine
Romanova Z.
assistant professor, Candidate of Technical Sciences, Biotechnology of Fermentation products and wine-
making, National University of Food Technologies, Kyiv, Ukraine
Romanov N.
Candidate of Technical Sciences, Integarated Automation of Control Systems, National University of Food
Technologies, Kyiv, Ukraine

Abstract
Purpose: The main raw material for beer production is malt. There is a tendency to use other grain crops in
the production of beer. The using of unmalted materials gives opportunity for brewers and technologist for cost
saving. But replacing of malt makes affects the quality of wort and beer. The purpose of the researches was to
compare quality parameters of 100% malt beers and beer with part of malt replaced by unmalted barley. Meth-
ods: To solve the given problem the methods of the high-performance liquid chromatography, analysis and em-
pirical methods are used. Results: Using a big amount of unmalted grits together with external enzymes gives
brewers opportunity to receive product in the same composition and quality like high-quality product using only
malt. Conclusions: For brewing with a very small amount of malt, barley is a very good choice for replacing
part of it as it has some characteristics of malt. It is possible, economically profitable and technically reasoned to
replace half of the malt in the grind. But it becomes necessary to use exogenous enzymes.

Keywords: barley; beer; cost; enzymes; malt; quality.

1. Introduction tioned activities, including thermostable alpha-

Barley brewing amylase. Products Ceremix are added during mashing.
For brewing with a very small amount of malt, The second option is to use as necessary individ-
barley is a very good choice for replacing part of the ual preparations Ultraflo, Neutrase, Fungamyl and
malt, as it has some its characteristics: Termamyl.
- contains beta-amylase and peptidase The first option is simpler; the second option has
- can be added to the mash tun, since barley the advantage that the brewer himself can more accu-
starch has a low gelatinization temperature rately choose a mixture of enzymes suitable for a spe-
- has a husk, which provides good filtration in the cific composition of the mash. If the part of barley
lauter-tun (with proper crushing) [1-2]. used is very high, then the final degree of fermentation
It is possible to replace half of the malt in the of the wort will be reduced. It can be corrected by
grind. For example, 30% of corn and 70% of malt adding Fungamyl during fermentation or Promozyme
replace 30% of corn grist, 35% of barley and 35% of - during mashing.
malt, and still produce high-quality beer. But it be- Brewing with barley and enzymes gives the
comes necessary to use exogenous enzymes [3]. brewer the following advantages:
For experimental part external enzymes from - the enzyme dosage can be refined to obtain the
Novozymes company were used. required level of activity in the mash
There two options for brewers where found: - less dependence on supply fluctuations, the
The first option is specially prepared mixtures price and quality of the malt
equivalent to malt: Ceremix 2X L, consisting of alpha- - obtaining a wort with a lower viscosity, which
amylase, protease and beta-glucanase; Ceremix 6X means an easier filtration of wort and beer
MG, which also contains cellulase and pentosanase; - cost reduction.
and Ceremix Plus MG, which possesses all the men-

Conditions:
Malt Extract content 75%
1 Barley Extract content 75%
Corn grits Extract content 75%

2 70% of malt + 30% of corn grits: 1000 kg of extract: 700 kg from malt + 300 kg from corn grits

35% of malt + 35% of barley + 30% of corn grits: 1000 kg of extract: 350 kg from malt + 350 kg
3
from barley + 300 kg from corn grits
4 Norwegian Journal of development of the International Science No 22/2018

Enzymes Activities Dosage rate

alpha-a Table 1. Enzymes and dosage rates 1.4-1.8 kg per 1 ton of barley during
Ceremix 2X L
mylase, beta-glucanase, protease mashing
alpha-amylase, beta-glucanase, protease, xy- 0.5 kg per 1 per ton of barley during
Ceremix 2X MG
lanase, cellulase mashing
alpha-amylase, beta-glucanase, xylanase, cellu- 0.7-1.0 kg per 1 ton of barley during
Ceremix Plus MG
lase, protease mashing
0.2-0.5 kg per 1 ton of barley during
Ultraflo L beta-glucanase, cellulase, xylanase
mashing
Neutrase 0.8 L neutral protease 0,6-1,0 kg per 1 ton of grinding
Fungamyl Brew Q maltogenic fungal alpha-amylase 1 kg per 1000 hl in fermentation tank
Promozyme Brew 0.5-1.5 kg per 1 ton of barley during
deamylating amylase
Q mashing
0.5 kg per 1000 hl in fermentation tank
Finizyme 250 L fungal beta-glucanase
or BBT

2. Analysis of the latest research and publica- quality of wort and beer. The purpose of the research-
tions es was to compare quality parameters of 100% malt
The exergetic performance of beer produced by beers and beer with part of malt replaced by unmalted
the conventional malting and brewing process is com- barley.
pared with that of beer produced using an enzyme- 4. Materials and methods of research
assisted process. The aim is to estimate if the use of an Using a big amount of unmalted grits together
exogenous enzyme formulation reduces the environ- with external enzymes gives brewers opportunity to
mental impact of the overall brewing process. [5] receive product in the same composition and quality
Whilst beers have been produced using various like high-quality product using only malt. To solve the
levels of unmalted grains as adjuncts along with malt, given problem the methods of the high-performance
brewing with 100 % unmalted grains in combination liquid chromatography, analysis and empirical meth-
with added mashing enzymes remains mostly un- ods are used.
known [6]. 5. Research results
The conversion of barley into beer represents Barley milling
mankind's oldest and most complex example of ap- As for malt milling and for milling of barley suit-
plied enzymology. Barley is one of the most ancient able as a dry or wet method of milling, you need only
crops, and it has evolved through domestication to a minimal adjustment.
today as a major world crop based on acreage and The wet method of milling is well suited in
production. It has great potential to reclaim some of its which the grains and their husks soften during steep-
prominence as a food grain, largely due to its high ing (but only when working with lauter tun) [8].
nutritional value [7]. When malt is replaced by a small amount of un-
3. Purpose and objectives of the study malted barley, the barley is mixed with the malt before
The main raw material for beer production is milling, but with a high proportion of the barley used
malt. There is a tendency to use other grain crops in it is preferable to mill the barley separately using two-
the production of beer. The using of unmalted materi- rollers mill with a 0.2 mm gap between the rollers
als gives opportunity for brewers and technologist for (normal gap for malt is 0.3-0.5 mm) [9-10].
cost saving. But replacing of malt makes affects the

Example 1
As an example of obtaining 1000 kg of extract from malt, barley and corn grits, the following variant can be
proposed:
Malt (extract content 76% 475 kg
Mash composition Barley (extract content 68%) 475 kg
Corn grits (extract content 82%) 400 kg
Enzyme for
Termamyl 120 L 0,15 kg
liquefaction:
Enzymes for Ultraflo XL 0,3 kg
mashing: Neutrase 0.8 L 0,3 kg
Norwegian Journal of development of the International Science No 22/2018 5

Pic.1. Mashing process, example 1.

Corn grits is liquefied separately using Termamyl After 60 minutes, the combined mash is heated to
stabilized with a total Ca2 + content of about 100 mg / 76-78 ° C and kept at this temperature until complete
L, with water: additives ratio of about 4: 1. Milled saccharification. Then the wort separation process
malt and barley mashed at the temperature of 52 ° C. follows.
After 30 minutes the boiled part is added to the The wort / beer properties will be similar to the
barley and malt mush to increase its temperature to wort / beer properties obtained using 70% malt and
63-66°C. 30% corn grits.

Example 2
Table 2
Gives some results of lab tests of mashing 50% of malt and 50% of barley using the products of Ceremix.
Parameters, units Control Ceremix 2X L Ceremix 6X MG
Enzymes dosage, w/w of
0 0,14 0,1 0,05 0,025
barley
Extract content, % 19,15 19,65 19,64 19,65 19,71
Total nitrogen, mg /l * 820 1123 1120 1068 945
FAN, mg/l 131 213 208 184 143
Viscosity, mPa*s * 2,84 1,75 1,59 1,63 1,69
Beta-glucan, mg / l * 2203 <30 <30 <30 114
* - for 12% wort

Example 3
The results of full-scale tests before the start of using barley brewing at the breweries in Bulgaria several years
ago.
Table 3
Raw materials and mashing type
Control Variant 1 Variant 2
% kg % kg % kg
Malt 85 5950 60 4200 60 4200
Barley 15 1050 40 2800 40 2800
Ceremix 2X L 0,15 4,2 0,15 4,2
Mashing type decoction infusion decoction
The technical characteristics of the mashing process are in Table 4.
6 Norwegian Journal of development of the International Science No 22/2018
Table 4
Technical characteristics of the mashing process
Characteristics Control Variant 1 Variant 2
pH during mashing 5,50 5,50 5,50
Saccharification time, min 10-15 10-15 10-15
Wort filtration time, min 208 170 200
Speed of filtration, l / m2 / min 13,3 15,7 12,9
Yield, % 68,8 72,3 75,3

6. Discussion of results The extract yield in the experimental brews is

higher than in the control one, and does not depend on
From the table: the amount of used unmalted barley. The yield using
The saccharification time in experimental brews the decoction mashing method is 3% higher than using
does not differ from the saccharification time in the infusion one.
control variant and is within 10-15 minutes. Physico-chemical parameters of the wort quality
For decoction mashing method the filtration are in Table 5.
speed of the experimental wort was almost identical to Table 5
the control. The wort obtained by the infusion method Physico-chemical parameters of the
had a higher speed of filtration.
wort quality
Characteristics Control Variant 1 Variant 2
Extract, % 9,35 9,15 9,13
pH 5,25 5,35 5,40
Color, EBC* 17 17 17
Viscosity, mPa*s * 1,72 1,70 1,64
Beta-glucan, mg / l * 590 167 339
Soluble nitrogen, mg / l * 826 954 937
Alpha-amine nitrogen, mg / l * 139,00 149,00 142,00
Zn2 +, mg / l 0,22 0,21 0,23
* - for 9% wort
It should also be noted that the sugar spectrum that the concentration of individual amino acids is the
determined by HPLC (high-performance liquid chro- same or higher in the experimental brews than in the
matography) were almost identical. As can be seen control one. The only amino acid present in the exper-
from Table 5 total soluble nitrogen and free amine imental brew at a lower concentration is arginine, in
nitrogen are higher in brews using barley and en- the control variant, 9.3 mg / 100 ml of arginine was
zymes. found, compared to 8.6 and 8.4 mg / 100 ml in the
Content analysis of individual amino acids shows experimental brews.
(together with higher values of free amine nitrogen) Finished beer analysis are in Table 6.
Table 6
Finished beer analysis
Characteristics Control Variant 1 Variant 2
First wort extract content, % mass. 9,30 9,14 9,08
Visible extract,% of mass. 2,20 2,10 1,90
Alcohol,% of mass. 2,91 2,89 2,95
Apparent degree of fermentation,% 76,3 77,1 79,1
The difference between the final and
0,2 0,4 0,4
reached degree of fermentation,%
pH 4,4 4,5 4,5
Color, EBC* 14 13 13
Viscosity, mPa*s * 1,65 1,58 1,58
Beta-glucan, mg / l * 211 61 158
Foam stability, Ross and Clark 130 136 136
VDK, mg/l* 0,31 0,24 0,25
Esters, mg / l 22,4 27,1 20,3
Higher alcohols, mg / l 75 79 78
* - for 9% wort
Norwegian Journal of development of the International Science No 22/2018 7
Surprisingly, that VDK content was high in all https://www.sciencedirect.com/science/article/pii/S03
variants - this could be corrected using the Maturex. 60544216312518
Conclusions: 6. Zhuang Sh., Shetty R., Hansen M., Fromberg
Using big amount of unmalted grits together with A., Hansen P.B., Hobley T.J., Brewing with 100 %
external enzymes gives brewers opportunity to receive unmalted grains: barley, wheat, oat and rye, European
wort/beer in the same composition and quality like Food Research and Technology 243(3) · August 2016,
high-quality wort/beer using only malt. Resource use:
https://www.researchgate.net/publication/306306548_
REFERENCES: Brewing_with_100_unmalted_grains_barley_wheat_o
1. Cege P., Kenyan Beer produced with unmalted at_and_rye
Barley, Paper presented at 7th IOB Convention, Nai- 7. Gupta H., Abu-Ghannam N., Gallaghar E.,
robi, 1999, Published in Ferment June/July 1999, No- Barley for Brewing: Characteristic Changes during
vozymes print A 6634 Malting, Brewing and Applications of its By-Products,
2. Graham S., Non-malted Adjuncts to produce April 2010, DOI: 10.1111/j.1541-4337.2010.00112.x,
fermentation Ethanol, Novozymes paper. 2013 Resource use:
3. Aastrup S.,Olsen H.S., Enzymes in Brewing, http://onlinelibrary.wiley.com/doi/10.1111/j.1541-
BIozoom no. 2, 2008 4337.2010.00112.x/full
4. Hannemann W., Brewing trial with specially 8. Kunze W., Technology Brewing and malting
made cheap malt and tailor-made Enzymes, Novo- International Edition, VlB Berlin, 1996
zymes publication 2001-16332-01 9. Houghs J.S., Briggs D.E., Stevens R.,
5. Van Donkelaar L. H.G, Mostert J., Zisopoulos Brewing and malting science, Chapman and Hall,
F.K., Boom R.M., Der Goot A-J, The use of enzymes London, 1971, 658.
for beer brewing: Thermodynamic comparison, 2016, 10. Aschengreen N.H., Brewing technology,
p. 519-527, Resource use: Internal brewing compendium: 1998, update Novo-
zymes 2003
8 Norwegian Journal of development of the International Science No 22/2018

CHEMICAL SCIENCES
REGULARITIES OF SORPTION REMOVING PB2+ BY BENTONITE CLAYS OF THE UKRAINIAN
ORIGIN

Fedenko Yu.
Ph.D., assistant, The Department of Technology of Inorganic Substances, Water Purification and General
Chemical Technology, Igor Sikorsky Kyiv Polytechnic Institute

Abstract
Regularities of sorption removing Pb2+ by bentonite clays of the ukrainian origin had been reviewed. The
initial concentration of Pb2+ on the bentonite met by the highest degree of its removing from water solutions had
been determined. The minimal duration (min) of sorption and sorbent dose keeps to achieve the highest remov-
ing degree of Pb2+ had been determined.

Keywords: sorption, bentonite, heavy metals, pollutant, drinking water, industrial water.

Introduction structural pantry can be expand [3, p. 5]. Such proper-

Rational using of water resources is one of the ty of bentonite clays as the ability to light widely used
key and the most actual task of ecology in conditions in vine industry and is perspective for using in water
of increasing water using in different branches of hu- treatment processes. Possibility of bentonite surface to
man activity. The role of the heavy metals content in be modified opens wide horizons for researches with
this aspect isn’t last. It’s necessary to note that metals the aim to give it properties to remove pollutants of
don’t lost in everything transformations take place in different nature and dispersion.
the environment with different substances. They move Spectra of bentonite clays using at the moment is
from one form to another one with formation of in- quite wide and expending constantly. They used for
cluding water-soluble compounds. This aspect de- removing as non-radioactive heavy metals (Cu2+,
creases the efficiency of water treatment by such Cd2+, Pb2+, Zn2+), and radioactive isotopes 137Cs та
90
widespread methods as coagulation, filtration, bio- Sr for granting lower radioactivity for spent fuel
chemical treatment etc. before disposal [4, p. 17]. However, all these methods
Sorption is the most efficient way keeps to defuse based on using sorbents of foreign origin (Russia in
industrial waste waters and to return purified water particular) can be differ for structure from the Ukrain-
and valuable components to industry. Besides, sorp- ian bentonites and, the most importantly, can be more
tion purification is cheap expressive method keeps to expensive.
remove wide spectra of pollutants with different phase The aim of this work was to find optimal condi-
and dispersion composition including ionic degree of tions of Pb2+ adsorption conditions by bentonite clays
dispersion. of Ukrainian origin for determination of their sorption
Pb2+ is toxic substance of complex action. MPC efficiency and selectivity.
of Pb2+ in drinking tap water (according to STATE Materials and methods
STANDARTS 2.2.4-171-10) is 0,01 mg/L [1, p. 9]. It had been used following reactants for the ex-
Accumulation of Pb2+ in organism by hitting from periment: Pb(NO3)2 (net analysis, STATE
drinking water eventually penalized functions of brain STANDART 4236-77), sodium etylendyamintetraace-
and central nervous system causing coma, convulsions tate (net analysis, STATE STANDART 10652-73),
and even death [2, p. 30]. That’s why scientific foun- acetate buffer mixture (pH 5,6), xylenol orange (net
dation and working of new and improving existing analysis, SPECIFICATIONS 6-09-1509-78), ammo-
methods of water purification from Pb2+ is an actual nia buffer mixture (pH 11), murexide
task. (SPECIFICATIONS 6-09-13-945-94), bentonite.
Sorption can be using in industrial and household Analysis of residual quantities of Pb2+ in water
conditions. solutions had been made by titrimetric method of
Among large variety of sorbents, the preference analysis. Pb2+ model water solutions had been pre-
is usually given for natural materials because it’s us- pared by dissolving of Pb(NO3)2 sample.
ing is the cheapest. Among them clay materials are For sorption conducting bentonite sample (1 g)
stand profitable, especially, saponite and bentonite placed to conical bulbs. Then solutions with concen-
clays with numerous deposits in the territory of trations 1; 0,75; 0,5; 0,2; 0,1 g/L (Pb2+) had been add-
Ukraine (one of the most powerful deposits in the ed. Bulbs closed tightly by plugs and put on mixing
world is situated in the territory of The Cherkassy during 5, 10, 15, 30 min. at room temperature.
region). For building adsorption isotherm of Pb2+ on the
The bentonite (called for the deposite Benton, bentonite solutions with different concentrations
USA) is natural clay mineral, hydro alumina silicate, (0,75; 0,5; 0,2; 0,1 g/L Pb2+) and volume 100 ml had
has the property to swell during hydration. The main been prepared. Bentonite samples (1 g) had been
component of the bentonite (60 – 70 %) montmorillo- placed to conic bulbs. Then prepared solutions had
nite Al2[Si4O10](ОH)2·nH2О is leaf like silicate with been added there. The duration of shaking was 15
Norwegian Journal of development of the International Science No 22/2018 9
minutes. After mixing solutions had been filtrated. moving of Pb2+ bentonite samples (2,5 g) had been
Then residual concentration of heavy metals had been placed into conical bulbs and solution of Pb2+ (0,5 g/L,
determined. 100 ml) had been added there. pH correction had been
For researching of Pb2+ sorption removing ben- carried by HCl and NaOH solutions with concentra-
tonite samples (0,1; 0,2; 1; 2 and 5 g) had been placed tions 0,1 M. pH of model solutions had been estab-
to conical bulbs. Then Pb2+ solution with concentra- lished (2, 4, 6, 8). Bulbs had been placed for mixing
tion 0,5 g/L and volume 100 ml had been added there. during 15 min.
Bulbs had been placed for mixing during 15 min. For Experimental
building of kinetic curves bentonite samples (2,5 g) In figure 1 it’s presented adsorption isotherm of
had been placed into conical bulbs and Pb2+ solution Pb2+ on the bentonite (duration of the process 15 min.,
(0,5 g/L, 100 ml) had been added. Bulbs had been pH 8 (without correction), dose of sorbent – 1 g for
placed for mixing during 5, 10,15, 30 min. For re- 100 ml of the solution).
searching pH influence on effectively of sorption re-

Figure 1 – Sorption isotherm of Pb2+ on the bentonite.

According to sorption isotherm presented in fig- concentration of adsorptive [5, p. 24, 6, p. 17]. Maxi-
ure 1 investigated sorbent on the base of the bentonite mum effectively of sorption removing of Pb 2+ up to 30
clay is mainly microporous, but behavior of adsorb on mg/g watched for equilibrium concentration of Pb 2+
surface, according to Brunnauer’s classification, obvi- 200 mg/L.
ously, subjects to model of Langmuir in the asymptot- In figure 2 it’s presented kinetic curve of sorption
ic approximation of adsorbed Pb2+ quantity to the ca- removing of Pb2+ under similar conditions as for ad-
pacity of monomolecular layer. According to Gills sorption isotherm (figure 1) and primary concentration
model isotherm of adsorption belongs to the class L of Pb2+ 200 mg/L.
type 2 with reaching of the adsorption edge for certain

Figure 2 – Kinetic curve of Pb2+ sorption removing

According to figure 2 in researched period of removing doesn’t change. That’s why, adsorption
time maximum Pb2+ removing percentage reaches 98 charging during more than 30 min doesn’t appropriate.
% for 30 min. For longer duration effectively of Pb 2+ Characteristic of equilibrium installation indicates a
10 Norwegian Journal of development of the International Science No 22/2018
low selectivity of the bentonite sorbent in relation to Pb2+.

Figure 3 – Depending of Pb2+ removing percentage from sorbent dose.

In figure 3 it’s presented depending of Pb 2+ re- adhesion of sorbent particles with a high content of it
moving percentage from sorbent dose (duration 30 in solution. It causes blocking of useful sorption sur-
min, pH 8, primary concentration of Pb 2+ 200 mg/L). face of the bentonite.
According to data presented in figure 3 increas- In figure 4 its illustrated depending of Pb2+ re-
ing sorbent dose from 0,1 to 1 g/100 ml leads to rapid moving percentage from pH of the solution (adsorp-
increasing of removing degree. However, from 1 to tion duration 30 min., sorbent dose 1 g/100 ml, prima-
5 g/100 ml efficiency of sorption removing of Pb2+ ry concentration of Pb2+ 200 mg/L).
almost unchanged. On our opinion, it connects with

Figure 4 – Duration of Pb2+ removing percentage from pH.

Figure 4 data show increasing pH causes increas- ditions from model water solutions are: primary con-
ing removing percentage of Pb2+, moreover, in the centration of Pb2+ 200 mg/L, sorbent dose 1g/100 ml
interval pH 4-6 it is abrupt. The transition of pH from of the solution, duration of sorption 30 min., primary
2 to 8 it caused gradual accession of hydroxyl group to pH 8.
Pb2+. Obviously, in the molecular form (Pb(OH) 2) lead Sorption of Pb2+ on the bentonite happens ac-
sorbs more effectively with prevalence, as we think, cording to Langmuir model and with behavior to im-
adsorption of the physical type. The adruption increas- proving with transition of Pb2+ from ionic to molecu-
ing of the removal percentage (from 50 to 90%) can lar form. It results, as we think, to predominance of
be explained by starting of Pb2+ transition from ionic the mechanism of non-specific physical adsorption.
to molecular form. These data can serve as scientific base for im-
Conclusions plementation to the industrial using cost-effective and
Bentonite clays are effectively domestic cheap effective sorbents of broad spectrum.
sorbents of Pb2+ keeps to clear water to concentrations
normalized for drinking water. Optimal sorption con-
Norwegian Journal of development of the International Science No 22/2018 11
REFERENCES: 4. Кулешова М.Л. Поглощающие свойства бе-
нтонитовых глин как материала для иммобили-
1. ГСанПиН 2.2.4-171-10. Гигиенические тре- зации радионуклидов при захоронении РАО /
бования к воде питьевой, предназначенной для М.Л.Кулешова, В.И.Сергеев, Т.Г.Шимко и др. //
потребления человеком [Текст]. – Взамен ГОСТ Научная конференция ЛОМОНОСОВСКИЕ
2874–82; введ. 2010–05–12. – Киев: Министерство ЧТЕНИЯ МГУ. – 2013. – С. 178-181.
здравоохранения Украины; К.: Изд-во стандартов, 5. Нецкина О.В. Практикум по физической
cop. 2010. – 25 с. химии НГУ. Химическая термодинамика и кине-
2. Адсорбция из растворов на поверхности тика. Адсорбция из растворов на твердой поверх-
твердых тел / Под ред. Г. Парфит и К. Рочестер. ности: метод. Пособие / О.В. Нецкина; Новосиб.
М.: Мир, 1986. – 488 с. гос. ун-т. – Новосибирск: РИЦ НГУ, 2015. – 17 с.
3. Переломов Л.В. Адсорбция свинца натри- 6. Giles C.H. Studies in adsorption: A system of
евым бентонитом и бентонитом, модифицирован- classification of solution adsorption isotherms / C.H.
ным гидроксидом алюминия, в присутствии орга- Giles, T.H. MacEvan, S.N. Nakhwa et.al. // Journal of
нических кислот / Л.В. Переломов, Н.Л. Лагунова, the Chemical Society. – 1960. - № 3. – P. 3973-3993.
И.В. Переломова и др. // Известия ТулГУ . Техни-
ческие науки. – 2013. – Вып. 6. Ч. 2. – С. 237-245.
12 Norwegian Journal of development of the International Science No 22/2018

JURISPRUDENCE
УДК 343.9:343.533
CURRENT STATUS AND TRENDS OF INTELLECTUAL PROPERTY CRIMES IN UKRAINE

Novikov O.
Candidate of Juridical Science,
Assistant of department of criminology and penitentiary law,
Yaroslav Mudryi National Law University
Kharkov, Ukraine

СОВРЕМЕННОЕ СОСТОЯНИЕ И ТЕНДЕНЦИИ РАЗВИТИЯ ПРЕСТУПНОСТИ В СФЕРЕ

ИНТЕЛЛЕКТУАЛЬНОЙ СОБСТВЕННОСТИ В УКРАИНЕ

Новиков О.В.
Кандидат юридических наук, ассистент кафедры
криминологии и уголовно-исполнительного права,
Национальный юридический университет имени Ярослава Мудрого
г. Харьков, Украина

Abstract
The article deals with criminological characteristics of crime in the sphere of intellectual property in
Ukraine during 2007–2017. The author mentions that an essential amount of crimes in the sphere of intellectual
property are annually registered in Ukraine and dynamics of researched crime is specified by some variations.
Following the analysis of quality and quantity indices it has been concluded on negative tendencies related to the
development of crime in the sphere of intellectual property.
Аннотация
В статье представлена криминологическая характеристика преступности в сфере интеллектуальной
собственности в Украине за 2007–2017 гг. Автор отмечает, что в Украине ежегодно регистрируется зна-
чительное количество преступлений в сфере интеллектуальной собственности, а динамика исследуемой
преступности характеризуется некоторыми колебаниями. На основании анализа количественных и каче-
ственных показателей, делается вывод о негативных тенденциях развития преступности в сфере интел-
лектуальной собственности.

Keywords: intellectual property; pirate goods trafficking; crimes against intellectual property; crime in the
sphere of intellectual property; indices of crime measurement in the sphere of intellectual property.
Ключевые слова: интеллектуальная собственность; оборот контрафактной продукции; преступле-
ния против интеллектуальной собственности; преступность в сфере интеллектуальной собственности;
показатели измерения преступности в сфере интеллектуальной собственности.

Постановка проблемы. На пути к обретению По оценкам Украинского альянса по борьбе с

Украиной членства в Европейском Союзе и подделками и пиратством, ежегодные объемы тор-
успешной интеграции в постиндустриальную мо- говли подделками в Украине могут оцениваться в
дель общества, необходимо обеспечить устойчи- 1,3 миллиарда долларов США [2, с. 12]. Значи-
вое развитие экономики государства, а также вы- тельная часть доходов от нарушения прав интел-
полнить ряд международных обязательств. Среди лектуальной собственности поступает в распоря-
которых не последнее место занимает обеспечение жение организованным преступным группам, ко-
надлежащей охраны и защиты права интеллекту- торые в дальнейшем направляют полученные
альной собственности и противодействие пре- средства на финансирование других видов пре-
ступным посягательствам в этой сфере. ступной деятельности, а также на обеспечение
В своём Специальном отчете за 2018 г. Торго- функционирования в Украине «теневой экономи-
вое представительство США отнесло Украину к ки».
странам (Priority Watch List), в которых существу- Единичные преступления, совершаемые в
ют существенные проблемы в сфере охраны прав сфере интеллектуальной собственности, является
интеллектуальной собственности. Констатируется проявлением мощного явления соответствующей
широкое использование нелицензионного про- преступности, угрожающей экономической без-
граммного обеспечения государственными орга- опасности Украины, жизни и здоровью граждан,
нами Украины, а также неспособность правитель- подрывает международный авторитет государства.
ства внедрить эффективные средства борьбы с В связи с этим, преступность в сфере интеллекту-
распространенным в стране интернет-пиратством альной собственности является одной из важней-
[1, с. 56–57]. ших проблем Украины и требует постоянного изу-
чения.
Norwegian Journal of development of the International Science No 22/2018 13
Анализ последних исследований и публи- или использование сведений, составляющих ком-
каций. Криминологический анализ показателей мерческую и банковскую тайну (ст. 231 УК Укра-
измерения преступности в сфере интеллектуаль- ины) и разглашение коммерческой или банковской
ной собственности и отдельных её проявлений тайны (ст. 232 УК Украины). Последние два со-
проводили в своих работах В. Б. Харченко, става относятся к преступлениям против интел-
Л. А. Безуглый, Е. П. Орлюк, А. С. Нерсесян и не- лектуальной собственности только в части неза-
которые другие украинские ученые. конного получения и разглашения сведений, со-
Целью статьи является выяснение современ- ставляющих коммерческую тайну.
ного состояния и тенденций развития преступно- Во-вторых, при осуществлении криминологи-
сти в сфере интеллектуальной собственности в ческого анализа преступности в сфере интеллекту-
Украине. альной собственности мы использовали два источ-
Методология исследования. Наше исследо- ника информации: 1) статистическую отчётность
вание следует начать с некоторых пояснений. Во- департамента информационно-аналитического
первых, под преступностью в сфере интеллекту- обеспечения МВД Украины (с 2013 г. – статисти-
альной собственности мы понимаем совокупность ческую отчётность Управление организационного
преступных посягательств на установленный по- обеспечения Единого реестра досудебных рассле-
рядок владения, пользования и распоряжения объ- дований и информационно-аналитической работы
ектами права интеллектуальной собственности, в Генеральной прокуратуры Украины); и 2) резуль-
результате которых причиняется значительный таты выборочного исследования уголовных дел
ущерб конкретным правообладателям. (уголовных производств) за преступления против
Предложенное понятие является абстрактной интеллектуальной собственности.
моделью реальных преступных проявлений и Эмпирическую базу исследования составили
включает в себя два аспекта: содержательный и 179 уголовных дел (уголовных производств) по
формальный. Содержанием преступности в сфере статьям 176, 177, 229, 231, 232 УК Украины, рас-
интеллектуальной собственности является обще- смотренных с 2007 по 2015 гг. местными судами
ственно опасное явление нарушения права интел- Харьковской, Киевской, Ровенской, Херсонской,
лектуальной собственности в виде противоправно- Львовской, Полтавской областей, АР Крым, г. Ки-
го использования исключительных имуществен- ева и г. Севастополя, а также экспертный опрос 80
ных прав субъектов права интеллектуальной специалистов в сфере интеллектуальной собствен-
собственности, в результате которого причиняется ности Научно-исследовательского института ин-
ущерб в значительном размере. Формой существо- теллектуальной собственности, Государственной
вания такой преступности является относительно службы интеллектуальной собственности, Управ-
массовая деятельность членов общества путем ления Государственной службы борьбы с эконо-
совершение преступлений против права интеллек- мической преступностью и Управления борьбы с
туальной собственности. Иными словами, вся со- киберпреступностью МВД Украины в Харьков-
вокупность преступлений против интеллектуаль- ской области.
ной собственности будет эмпирическим проявле- Изложение основного материала. Традици-
нием преступности в сфере интеллектуальной онно в криминологии для анализа определенного
собственности. вида преступности используют количественные
Следует отметить, что среди криминалистов и (уровень, коэффициенты), качественные (структу-
криминологов нет однозначного мнения относи- ра, география) и количественно-качественные (ди-
тельно круга преступлений, которые непосред- намика) показатели. В связи с этим, считаем целе-
ственно посягают на право интеллектуальной соб- сообразным использовать их для изучения состоя-
ственности. Но большинство ученых соглашаются, ния преступности в сфере интеллектуальной
что к ним относятся: нарушение авторского и собственности.
смежных прав (ст. 176 Уголовного кодекса Украи- Важным количественным показателем пре-
ны (в дальнейшем – УК Украины)); нарушение ступности является ее уровень, то есть зафиксиро-
права на изобретение, полезную модель, промыш- ванное в абсолютных величинах количество заре-
ленный образец, топографию интегральной мик- гистрированных преступлений и лиц, их совер-
росхемы, сорт растений, рационализаторское шивших, в Украине полностью или в отдельных ее
предложение (ст. 177 УК Украины); незаконное регионах за определенный период времени [3,
использование знака для товаров и услуг, фирмен- с. 56]. Уровень преступности в сфере интеллекту-
ного наименования, квалифицированного указания альной собственности за 2007–2017 гг. отражен в
происхождения товара (ст. 229 УК Украины); а табл. 1.
также незаконный сбор с целью использования
14 Norwegian Journal of development of the International Science No 22/2018
Таблица 1
Количество зарегистрированных преступлений против интеллектуальной собственности в Украине за
2007–2017 гг.
Зарегистрированные преступления
(по статьям УК Украины)
Годы
ст. 176 УК ст. 177 УК ст. 229 УК ст. 231 УК ст. 232 УК
Итого
Украины Украины Украины Украины Украины
2007 417 15 77 5 3 517
2008 435 12 42 3 1 493
2009 476 10 56 2 1 545
2010 579 7 78 2 3 669
2011 527 9 153 3 1 693
2012 140 18 91 6 12 267
2013 452 26 327 17 12 834
2014 339 16 250 26 10 641
2015 249 20 131 14 5 419
2016 157 19 168 19 7 370
2017 95 16 121 64 3 299
Проанализировав полученные данные, можно дит постепенное снижение количества зарегистри-
сделать вывод, что ежегодно в Украине регистри- рованных преступлений в этой сфере, что, учиты-
руется значительное количество преступлений в вая экономический кризис, массовое обнищание
сфере интеллектуальной собственности. Всего в населения, девальвацию национальной валюты,
течение 2007–2017 гг. (за одиннадцать лет) в больше свидетельствует о повышении уровня ла-
Украине было зарегистрировано 5747 этих пре- тентности этой группы преступлений и неэффек-
ступлений, то есть в среднем 522 преступления в тивности работы правоохранительных органов,
год. Начиная с 2008 по 2013 гг. суммарное количе- чем об улучшении ситуации с охраной прав ин-
ство всех зарегистрированных преступлений, со- теллектуальной собственности в стране.
вершаемых в сфере интеллектуальной собственно- Количество выявленных лиц, совершивших
сти, постепенно увеличивается (за исключением преступления против интеллектуальной собствен-
2012 г.). После же 2013 г. и по сей день происхо- ности, за 2007–2017 гг. отражено в табл. 2.
Таблица 2
Количество выявленных лиц, совершивших преступления против интеллектуальной собственности в
Украине, за 2007–2017 гг.
Выявлено лиц
(по статьям УК Украины)
Годы
ст. 176 УК ст. 177 УК ст. 229 УК ст. 231 УК ст. 232 УК
Итого
Украины Украины Украины Украины Украины
2007 136 1 0 0 0 137
2008 165 2 0 0 0 167
2009 202 2 0 0 0 204
2010 264 1 47 1 1 314
2011 280 4 59 2 0 345
2012 8* 1* 2* 0* 0* 11*
2013 104 2 35 5 2 148
2014 64 0 36 23 0 123
2015 84 0 12 8 0 104
2016 10 0 15 9 0 34
2017 18 0 12 59 0 89
* По состоянию на 20.11.2012.

Как видим, в течение 2007–2017 гг. количе- преступников на 5747 преступлений. Таким обра-
ство выявленных лиц, совершивших преступления зом, правоохранительными органами выявляется в
в сфере интеллектуальной собственности, харак- 3,43 раза меньше виновных лиц по сравнению с
теризуется определенными незначительными ко- количеством уголовных правонарушений.
лебаниями. В 2012 г. и 2016 г. было выявлено Дополняют сведения о преступности показа-
наименьшее количество таких лиц за последние тели уровня судимости. Уровень судимости – это
одиннадцать лет. Также обращает на себя внима- количество лиц, привлеченных к уголовной ответ-
ние определенный разрыв в соотношении количе- ственности и осужденных к определенным видам
ства зарегистрированных преступлений и числен- наказаний [3, c. 56]. Уровень судимости за пре-
ности выявленных преступников. Так, за исследу- ступления в сфере интеллектуальной собственно-
емый период времени было выявлено всего 1676 сти отображен в табл. 3.
Norwegian Journal of development of the International Science No 22/2018 15
Таблица 3
Уровень судимости за преступления против интеллектуальной собственности в Украине за 2007–2017 гг.
Количество осужденных лиц
(по статьям УК Украины)
Годы
ст. 176 УК ст. 177 УК ст. 229 УК ст. 231 УК ст. 232 УК
Итого
Украины Украины Украины Украины Украины
2007 46 1 5 0 0 52
2008 69 1 11 0 0 81
2009 75 0 11 0 0 86
2010 115 1 8 0 0 124
2011 88 0 7 0 0 95
2012 120 1 29 2 0 152
2013 60 3 16 0 0 79
2014 32 0 12 3 0 47
2015 31 0 13 0 0 44
2016 10 0 6 3 0 19
2017 12 0 7 0 0 19
Итак, согласно официальным данным за удельный вес и соотношение различных видов
одиннадцать лет (2007–2017 гг.) в Украине было преступлений в их общем числе за определенный
осуждено 798 человек за преступления, совершен- период времени на определенной территории [5, c.
ные в сфере интеллектуальной собственности. 105]. Структурное распределение преступлений в
Уровень преступности не всегда дает пра- сфере интеллектуальной собственности в период с
вильное представление о масштабах ее распро- 2007 по 2017 гг. иллюстрирует рис. 1.
странения и интенсивности проявления. В таком Как видим, наибольшее количество преступ-
случае, абсолютные показатели преступности до- лений в структуре преступности в сфере интеллек-
полняются относительными величинами, в частно- туальной собственности – это нарушение автор-
сти: коэффициентом преступной интенсивности ского права и смежных прав (в среднем 67,3% от
и коэффициентом преступной активности насе- всех зарегистрированных преступлений этого ви-
ления [4, c. 158]. По состоянию на 2017 г. коэффи- да); на втором месте – незаконное использование
циент преступной интенсивности преступности в знака для товаров и услуг, фирменного наимено-
сфере интеллектуальной собственности составляет вания, квалифицированного указания происхож-
7 преступлений на 1 млн. населения, а коэффици- дения товара (26%); на третьем месте – нарушение
ент преступной активности 3 преступника на 1 права на изобретение, полезную модель, промыш-
млн. населения достигшего 14-летнего возраста. ленный образец, топографию интегральной мик-
Теперь перейдем к структуре преступности росхемы, сорт растений, рационализаторское
в сфере интеллектуальной собственности в предложение (2,9%); другие преступления в сфере
Украине. Структура преступности, как ее каче- интеллектуальной собственности – 3,8%.
ственная характеристика, представляет собой

Рис 1.
Структура преступности в сфере интеллектуальной собственности в период с 2007 по 2017 гг. в
Украине
16 Norwegian Journal of development of the International Science No 22/2018
Более подробно проанализировать структуру прав сопровождается дополнительной квалифика-
преступности возможно с помощью уголовно- цией по ст. 2031 УК Украины (незаконный оборот
правовых и криминологических показателей, ко- дисков для лазерных систем считывания, матриц,
торые были установлены по результатам выбороч- оборудования и сырья для их производства).
ного исследования 179 уголовных дел, проведен- Встречаются единичные случаи квалификации
ного автором. соответствующей статьи в совокупности ст. 301
Определенное значение для анализа преступ- УК Украины (ввоз, изготовление, сбыт и распро-
ности в сфере интеллектуальной собственности странение порнографических предметов).
имеет способ совершения преступления. Так, Следующим показателем преступности в сфе-
нарушение авторского права и смежных прав (ст. ре интеллектуальной собственности, который мы
176 УК Украины) чаще всего совершается или проанализируем, будет её география. В научной
путем реализации на торговых точках нелицензи- литературе под географией преступности понима-
онных дисков для лазерных систем считывания – ют её распространение по различным регионам
48,1%, или путем незаконного использования объ- страны (административно-территориальное рас-
ектов авторского права и смежных прав в пред- пространение) [6, c. 70]. Наибольшее количество
принимательской деятельности – 19 6% (как пра- выявленных преступлений в сфере интеллектуаль-
вило, использование нелицензионного программ- ной собственности регистрируется в г. Киев,
ного обеспечения при оказании определенных Харьковской, Одесской, Донецкой, Львовской,
услуг). Третью позицию занимает реализация дис- Киевской, Днепропетровской, Херсонской обла-
ков для лазерных систем считывания – 17,3%. До- стях Украины. Таким образом, преступность в
вольно распространены также случаи нарушения сфере интеллектуальной собственности распро-
авторских прав лицами, оказывающими коммерче- странена в густонаселенных промышленных обла-
ские услуги по настройке компьютерной техники. стях, где сосредоточен товарный и денежный обо-
Что касается незаконного использования знака для рот. Существует определенная особенность терри-
товаров и услуг, фирменного наименования, ква- ториального деления преступлений,
лифицированного указания происхождения товара предусмотренных ст. 229 УК Украины. Так,
(ст. 229 УК Украины), то чаще всего встречается наибольшее количество этих преступлений реги-
продажа контрафактных товаров известных произ- стрируется в г. Киеве, Харьковской и Одесской
водителей (78%), а также подделка таких товаров областях, то есть в тех областях, которые являются
(18%). ключевыми транспортными узлами, через которые
Способ совершения преступления определен- зарубежные товары попадают в Украину.
ным образом определяет место совершения пре- Теперь рассмотрим динамику преступности в
ступления. Наиболее распространенными местами сфере интеллектуальной собственности. Динамика
совершения преступлений в сфере интеллектуаль- преступности – это комплексный показатель, ха-
ной собственности являются: магазин (торговая рактеризующий преступность с количественно-
точка) – 46,8%, рынок – 31,2%, улица – 11%, про- качественной стороны. Она отражает изменения
изводственное помещение – 6%, жилое помещение уровня, интенсивности и структуры преступности
– 2%. или отдельного её вида в течение определенного
Что касается преступлений, совершаемых в периода времени [5, c. 169]. Динамика преступно-
совокупности с преступлениями в сфере интеллек- сти в сфере интеллектуальной собственности
туальной собственности, то довольно часто изображена на рис. 2.
(38,5%) нарушение авторского права и смежных

Рис 2.
Динамика преступности в сфере интеллектуальной собственности в период с 2007 по 2017 гг. в
Украине
Norwegian Journal of development of the International Science No 22/2018 17
Как видим, начиная с 2008 по 2011 гг. проис- и по сей день. Вновь созданные или реформиро-
ходит постепенное увеличение количества зареги- ванные подразделения правоохранительных орга-
стрированных преступлений в сфере интеллекту- нов, деятельность которых направлена на проти-
альной собственности (на 28,9% в 2011 г.), вплоть водействие преступности в сфере интеллектуаль-
до 2012 г., когда было зарегистрировано наимень- ной собственности, на сегодняшний день еще не
шее количество преступлений этого вида – 267 (– имеют достаточного опыта для выполнения воз-
61,5%). В 2013 г. произошло резкое увеличение ложенных на них полномочий. Особенно это каса-
зарегистрированных преступлений против интел- ется расследования рассматриваемой группы пре-
лектуальной собственности – 834 (+ 312,4%), что ступлений, совершаемых в глобальной сети Ин-
объясняется следующим. Так, с 01.01.2012 по тернет. А именно там сегодня сосредоточено
20.11.2012 было зарегистрировано только 38 пре- наиболее массовое нарушение прав интеллекту-
ступлений (!) в сфере интеллектуальной собствен- альной собственности в Украине;
ности, а затем за период 21.11.2012 – 31.12.2012 3) с 2014 г. некоторые нетяжкие (и в свою
регистрируется основная доля (229 преступлений). очередь наиболее распространённые) преступле-
Именно после 20.11.2012 ведение уголовной ста- ния в сфере интеллектуальной собственности пе-
тистики перешло от Департамента информацион- решли в разряд дел частного обвинения. В соот-
но-аналитического обеспечения МВД Украины в ветствие со ст. 477 Уголовного процессуального
Управление организационного обеспечения Еди- кодекса Украины, такое производство может быть
ного реестра досудебных расследований и инфор- начато следователем или прокурором лишь на ос-
мационно-аналитической работы Генеральной новании заявления потерпевшего. В связи с про-
прокуратуры Украины. Из-за этого, возникли цессуальными проблемами установления правооб-
определенные сложности с регистрацией преступ- ладателей, или же их нежеланием принимать уча-
лений, совершенных в сфере интеллектуальной стие в официальном расследовании, и возникает
собственности. После 2013 г. уровень зарегистри- снижение возможностей правоохранительных ор-
рованных преступлений в сфере интеллектуальной ганов в борьбе с нарушением прав интеллектуаль-
собственности ежегодно идет на спад. В 2017 г. ной собственности.
было выявлено наименьшее количество преступ- В рамках анализа статистически значимых
лений в этой сфере за последние пять лет (–64,1%). показателей, отражающих социальное бытие пре-
Подобную динамику преступности в сфере интел- ступности в сфере интеллектуальной собственно-
лектуальной собственности в Украине можно объ- сти, особое внимание следует уделить также ис-
яснить следующим: следованию латентности этой группы преступ-
1) с 2013 г. в Украине начался экономический лений. Так, существуют определенные сомнения
кризис, который обострялся торговой войной с относительно полноты выявления и регистрации, а
Российской Федерацией, девальвацией нацио- также объективности тенденций развития указан-
нальной валюты, потерей Автономной Республики ной преступности. Поэтому, приведенные данные
Крым и антитеррористической операцией на во- об уровне, структуре, динамике преступности в
стоке Украины. В связи с этим, Торговое предста- сфере интеллектуальной собственности лишь ча-
вительство США лишило Украину статуса «стра- стично отражают ее реальное состояние. Таким
ны-пирата №1» и одновременно сняло с повестки образом, без установления объемов латентности
дня вопрос о применении экономических санкций, этой преступности невозможно получить полное
которые угрожали стране в случае отсутствия су- представление о ней.
щественных сдвигов по предотвращению преступ- Под латентной преступностью следует пони-
ности в сфере интеллектуальной собственности. мать совокупность фактически совершенных, од-
Поэтому проблема нарушения прав интеллекту- нако не выявленных или таких, что в результате
альной собственности в Украине временно отошла других определенных обстоятельств не стали из-
на второй план. Вышеприведенное, лишний раз вестны правоохранительным и судебным органам,
подтверждает мысль В. Б. Харченко о политиче- преступлений, сведения о которых в связи с этим,
ском характере таких экономических санкций и не находят отражения в официальной уголовно-
ограничений, а также отсутствие их взаимосвязи с правовой статистической отчетности [8]. Результа-
уровнем соблюдения или нарушения прав на объ- ты экспертного опроса уровня латентности пре-
екты интеллектуальной собственности [7, c. 73]; ступности в сфере интеллектуальной собственно-
2) с 2014 г. в Украине началась реформа пра- сти в Украине приведены в табл. 4.
воохранительной системы, которая продолжается
Таблица 4
Результаты экспертного опроса уровня латентности преступности в сфере интеллектуальной собственно-
сти в Украине
Отношение количества незарегистрированных преступлений на одно зареги-
Распределение
стрированное
ответов экс-
1:3 1:4 1:5 1:6 1:7 1:10 >1:10
пертов
1,4 % 1,4 % 2,8 % 1,4 % 1,4 % 80,5 % 11,1 %

Таким образом, по мнению большинства экс- одно из десяти совершенных преступлений в сфе-
пертов (80,5%) в Украине регистрируется только ре интеллектуальной собственности. Определен-
18 Norwegian Journal of development of the International Science No 22/2018
ная часть опрошенных (11,1%) считает, что реги- 3. Кримінологія : підруч. / В. В. Голіна, Б. М.
стрируется более чем одно из десяти преступлений Головкін, М. Ю. Валуйська та ін.; за ред. В. В. Го-
этой группы. Как видим, подавляющее большин- ліни, Б. М. Головкіна. Харків : Право, 2014. 440 с.
ство экспертов считают преступления в сфере ин- 4. Закалюк А. П. Курс сучасної української
теллектуальной собственности высоколатентными кримінології: теорія і практика: у 3 кн. Київ : Ви-
и лишь незначительная часть – среднелатентными. давничий Дім «Ін Юре», 2007. Кн. 1. 424 с.
Выводы. Подытоживая вышесказанное, мож- 5. Криминология : учеб. для вузов / под общ.
но сделать вывод, что преступность в сфере ин- ред. докт. юрид. наук, проф. А. И. Долговой. 3-е
теллектуальной собственности имеет определен- изд., перераб. и доп. Москва : Норма, 2007. 912 с.
ные, свойственные только ей тенденции развития, 6. Іванов Ю. Ф., Джужа О. М. Кримінологія :
а изучение состояние этой преступности на сего- навч. посіб. 2-ге вид., доповн. та перероб. Київ :
дняшний день в Украине остается актуальным для Вид. ПАЛИВОДА А. В., 2008. 292 с.
дальнейших исследований. 7. Харченко В. Б. Протидія інтелектуальному
піратству в Україні: забезпечення охорони прав
СПИСОК ЛИТЕРАТУРЫ: інтелектуальної власності або спосіб впливу на
державну політику та механізм недобросовісної
1. 2018 Special 301 Report – USTR. конкуренції? Теорія і практика інтелектуальної
URL : https://ustr.gov/sites/default/files/files/Press/Re власності. 2015. № 3. С. 66–77. URL:
ports/2018%20Special%20301.pdf. http://nbuv.gov.ua/UJRN/Tpiv_2015_3_10.
2. Ukraine: Розвиток і захист інтелектуальної 8. Оболенцев В. Ф. Латентна злочинність:
власності у Україні. URL : проблеми теорії та практики попередження. Харків
http://sips.gov.ua/i_upload/file/Ukraine_IP_Report_U : Вид. СПД ФО Вапнярчук Н. М., 2005. 128 с.
KRAINE_LR.pdf.
Norwegian Journal of development of the International Science No 22/2018 19

TECHNICAL SCIENCES
DEVELOPMENT OF A MOBILE APPLICATION FOR TRACKING GEOLOCATION POSITION IN
REAL TIME

Shibaev D.
Master student, Department of high technologies, Irkutsk National Research Technical University
Dorofeev A.
Candidate of Technical Sciences, Associate Professor, Department of high technologies,
Irkutsk National Research Technical University

РАЗРАБОТКА МОБИЛЬНОГО ПРИЛОЖЕНИЯ ДЛЯ ТРЕКИНГА ГЕОЛОКАЦИОННОГО

ПОЛОЖЕНИЯ В РЕАЛЬНОМ ВРЕМЕНИ

Шибаев Д.П.
Магистрант, Институт высоких технологий,
Иркутский национальный исследовательский технический университет
Дорофеев А.С.
Кандидат технических наук, доцент, Институт высоких технологий,
Иркутский национальный исследовательский технический университет

Abstract
Today it is difficult to imagine life without smartphones, which are equipped with various functionalities,
including location tracking. The article describes the developed mobile application for tracking geolocation posi-
tion in real time. This application will allow you to track the location of certain groups of people in real time, as
well as to establish communication with a specific group via chat, for example.
Аннотация
Сегодня сложно представить себе жизнь без смартфонов, которые оснащены различными функцио-
нальными возможностями, в т.ч. функцией отслеживания местоположения, геолокации. В статье приво-
дится описание разработанного мобильного приложения для трекинга геолокационного положения в
реальном времени. Данное приложение позволит отслеживать местоположение определенных групп лю-
дей в реальном времени, а также установить связь с конкретной группой, например, с помощью чата.

Keywords: location tracking, geolocation, mobile application, Java, Android, UML diagram.
Ключевые слова: геолокация, отслеживание местоположения, мобильное приложение, Java, An-
droid, диаграмма UML.

На 80 % смартфонов установлена операцион- «Simple Latitude» позволяют отслеживать место-

ная система (ОС) Android, в связи с чем было при- положение пользователя, но не оснащены возмож-
нято решение разрабатывать приложение под дан- ностью общения в этих приложениях. При добав-
ную ОС. Android («Андроид») — операционная лении пользователей в список отслеживаемых, на
система, изначально разработанная компанией карте отображаются все добавленные пользовате-
Android Inc., впоследствии приобретенная компа- ли, что затрудняет визуальное восприятие. Также
нией Google, предназначена для смартфонов, отсутствует возможность отслеживания конкрет-
планшетов, электронных книг, цифровых проиг- ной группы пользователей.
рывателей, наручных часов, игровых приставок, В качестве спутниковых систем рассматрива-
нетбуков, очков Google, телевизоров и других лись системы GPS и ГЛОНАСС. В современных
устройств. В 2015 году появилась поддержка ав- мобильных устройствах имеется совмещенный
томобильных мультимедийных систем и бытовых приемник GPS/ГЛОНАСС, поэтому нет необходи-
роботов [1]. ОС Android основана на ядре Linux и мости выбора конкретной системы: приложение
собственной реализации виртуальной машины автоматически выберет спутник для получения
Java от Google. Android позволяет разрабатывать необходимых данных.
Java-приложения, управляющие устройством через В качестве языка программирования выбран
библиотеки созданные Google[2]. язык Java, который является основой практически
Сравнительный анализ аналогов приложения для всех типов сетевых приложений и всеобщим
показал ряд их достоинств и недостатков. Основ- стандартом для разработки и распространения
ным недостатком приложения «Номера» является встроенных и мобильных приложений, игр, веб-
то, что при регистрации все пользователи сети контента и корпоративного программного обеспе-
видят местоположение пользователя. Кроме того чения [3].
пользователи видят местоположение по номерам Инструменты реализации приложения:
автомобилей, что может привести к повышенному Android SDK (включает библиотеки, документа-
угону автомобилей. Приложения «Найти друзей» , цию и инструменты, которые помогают разраба-
20 Norwegian Journal of development of the International Science No 22/2018
тывать мобильные приложения для платформы приближения или удаления при просмотре карт),
Android), Android Virtual Device (эмулятор, ис- Firebase (база данных для хранения данных в ре-
пользуемый для проектирования, отладки и тести- жиме реального времени с их синхронизацией сра-
рования приложений в реальной среде выполне- зу на всех устройствах).
ния), Android Studio (среда разработки под На рис. 1 представлена диаграмма прецеден-
Android, основанная на IntelliJ IDEA), Google Maps тов разрабатываемой системы с использованием
API (бесплатная служба, позволяющая включить унифицированного языка моделирования (UML).
любую карту из Google Maps на внешнем сайте, Из диаграммы видны пользователи системы и ее
управлять этой картой через JavaScript, например, общая функциональность.
для добавления маркеров географических точек,

Рис. 1. Диаграмма прецедентов. Роли «Авторизованный пользователь» и «Администратор группы»

Интерфейс системы состоит из семи окон (рис. 2).
Norwegian Journal of development of the International Science No 22/2018 21

Рис. 2. Диаграмма классов. Интерфейс приложения

Весь интерфейс приложения реализован по- 1) Окно регистрации (рис.3). При запуске
средством расширяемого языка разметки XML. приложения пользователю необходимо зареги-
Android Studio предоставляет возможность разра- стрироваться в приложении. После регистрации,
ботки интерфейса приложения в тестовом виде зарегистрированный пользователь автоматически
XML, а также в графическом – с помощью разме- попадает на окно карты.
щения объектов на окне способом оперирования
drag and drop.

Рис. 3. Окно регистрации

22 Norwegian Journal of development of the International Science No 22/2018
2) Окно входа (рис. 4). При условии, что поль- для этого из окна регистрации, необходимо перей-
зователь уже зарегистрирован, но вышел из при- ти на окно входа.
ложения, ему необходимо войти в приложение,

Рис. 4. Окно входа приложения

3) Окно карты (рис. 5). Окно карты предо-

ставляет авторизированному пользователю отсле- Рис. 5. Окно карты приложения
живать собственное местоположение, а так же по- 4) Окно чата (рис. 6). В окне чата авторизо-
ложение участников группы, в которой он состоит, ванному пользователю доступна отправка сооб-
и которая выбрана. С окна карты возможен пере- щений, чтение сообщений других участников
ход в окно чата, так же переход в окно групп и группы, а так же выбор конкретной группы из
окно ПДД. В левом верхнем углу карты отобража- списка групп, в которых состоит участник.
ется скорость в км/ч.

Рис. 6. Окно чата приложения

Norwegian Journal of development of the International Science No 22/2018 23
Так как используется время FireBase (рис. 7), менном порядке, несмотря на то, что часы на
все сообщения выстраиваются в правильном вре- устройствах могут быть не синхронизированы.

Рис. 7. База данных FireBase. Хранение сообщений

5) Окно групп (рис. 8). Данное окно предна- другим зарегистрированным пользователем, груп-
значено для создания своей собственной группы, пу.
или для подачи заявки на вступление в созданную,

Рис. 8. Окно групп приложения

6) Окно управления группой (рис. 9). Адми- управления группой. Так же в данном окне суще-
нистратору группы необходимо принимать людей ствует возможность удаления пользователей со-
в группу, так же отклонять заявки неизвестных стоящих в группе, и удаление всей группы цели-
ему зарегистрированных пользователей. Для вы- ком.
полнения данных функций, необходимо окно
24 Norwegian Journal of development of the International Science No 22/2018

Рис. 9. Окно управления группой. 1) Заявка на вступление в группу. 2) Пользователь, состоящий в группе.

7) Окно правил дорожного движения (ПДД) можность выбрать необходимый пункт ПДД, а
(рис. 10). В окне ПДД пользователь имеет воз- также прочитать данные, выбранного, пункта.
Norwegian Journal of development of the International Science No 22/2018 25

Рис. 10. Окна ПДД. 1) Пункт ПДД. 2)Выбор пункта ПДД

В настоящее время разработанное приложе- СПИСОК ЛИТЕРАТУРЫ:

ние позволяет пользователям определять местопо-
ложения других пользователей, объединять их в 1. Обзор операционной системы Android
группы, общаться при помощи чата. В приложе- [сайт]. [2011]. URL: http://techcave.ru/groups/3 (дата
нии реализован спидометр для отслеживания ско- обращения: 10.09.2018).
рости по данным GPS. Для водителей автотранс- 2. Архитектура Android [сайт]. [2013]. URL:
порта встроен свод правил дорожного движения. В http://studbooks.net/1993712/informatika/opisanie_pr
ближайших планах развития данного приложения edmetnoy_oblasti (дата обращения 14.09.2018).
– разработка его под различные операционные 3. Объектно-ориентированное программиро-
системы, и распространение среди пользователей вание Java «Чат» [сайт]. [2015]. URL:
мобильных устройств. https://knowledge.allbest.ru/programming (дата об-
ращения 15.09.2018).

PROTECTION FROM THE ICEOF ELECTRICAL NETWORKS BY MEANS OF A HIGH-

FREQUENCY ELECTROMAGNETIC WAVE

Kaganov W.
Doctor of techn. Professor, Professor of the RTS Institute,
MIREA-Russian Technological Univers

Abstract
A new way to protect electrical networks from ice is discussed. Theoretical and experimental aspects of this
problem are considered. Information on the establishment of an industrial plant is given
Keywords: ice, electromagnetic waves, generator.

1. The essence of the problem. ing in high-voltage transmission lines, including Rus-
One of the serious causes of accidents in high- sia, Norway, Canada, China, Kazakhstan, Kyrgyzstan
power electric power systems is the formation of a and othe. Only in Russia, because of the ice in the
dense ice precipitate - icing - with the intentions of power systems, there are up to 6-8 serious accidents
supercooled droplets of rain or fog at temperatures per year. The ice, including its "icy rain" variety,
from 0 to -50 C on wires of high-voltage transmission causes enormous damage to the economies of many
lines.The thickness of the ice on transmission lines countries and a page to human casualties
can reach 100 mm or more, substantially weighting Melting is the main method of combating ice.
them and leading to rupture and breakage of supports The melting consists in heating up to 100-1300С of
(Fig. 1). Many countries of the world suffer from ic- short-circuited wires with ice with a constant or alter-
26 Norwegian Journal of development of the International Science No 22/2018
nating current of frequency 50 Hz with a power of sumed, consumers are disconnected for a considerable
more than 200 million volts-amperes and a current of time from the power grid and often unacceptable char-
more than 1000 A. As a result, a lot of energy is con- ring of wires.

Fig.1

2. Brief description of the new method. electromagnetic frequency of about 100 MHz - simul-
Consider the fundamentals of a fundamentally taneously spread along the lines of electrical transmis-
new method of combating ice in high-voltage air pow- sion, figuratively speaking, without interfering with
er networks. The method is protected by the RF patent each other. One wave (frequency 50 Hz) transmits
No. 2356148 dated May 20, 2009 [1] The method is electrical energy through wires, and the other (about
based on the heating of wires with the help of high- 100 MHz) - heats these wires. Short circuit in the net-
frequency electromagnetic waves, which should pre- work and disconnect from it, consumers, as in the
vent the formation of ice on them. melting, produce with this method is not required.
Eliminate the destructive effect of icing on the The transformation of the radiated NF energy in-
power line by heating the wires to a temperature of + to the heat is carried out with the help of special an-
10 -150 ° C . Such heating can be carried out using tenna-heaters with a resistance of 1-2 kOhm, which
high-frequency (HF) electromagnetic waves by con- with a certain step are located on the wires of the
verting their energy into heat. power line (Fig. 2).
To realize this idea, an RF generator is addition-
ally connected to the wires of the power line (Fig. 1).
Both waves - the industrial frequency of 50 Hz and the

Fig.2

Such RF-heating spots, placed with on the wires, V; γ – the conductivity of bodies located within the
heat the latter and ice on them is not formed. volume.
3. The theoretical substantiation of the pro- The first member of sum PG of the equation (1) is
posed meth a power of third-party sources (generators), giving its
There are two possible ways of conversion of mi- energy in the volume V; the second member of sum is
crowave energy into heat: with the help of running a power of electromagnetic energy converted to heat
way and standing wave. The both ways are based on because the active losses in the volume of V; the third
the Poynting equation that determines the time varia- member is the Poynting vector, that determines as
tion of the electromagnetic field energy W in the vol- P=EH and characterizes the energy radiated beyond
ume V : the considered volume V through the surface S.
dW The first option is associated with a traveling
 PG   γ E2 dV   EH dS , (1) electromagnetic wave propagating along a guide struc-
dt V S ture, in this case along a two- or three-wire line. Such
where E – the electric field vector; H – the mag- a wave moves along the elongated structure with ac-
netic field vector; S – a surface bounding the volume tive losses, and gives its energy. The second option is
Norwegian Journal of development of the International Science No 22/2018 27
related to the conversion of energy of electromagnetic i.e. a displacement of a current to the outer shell of
oscillations into heat in the oscillatory volumetric sys- conductor (see Charter 1).
tem in resonance mode, where the standing wave With increasing frequency of the current flowing
mode is established in a resonant cavity. through the wire, the active resistance of the latter
The conversion of energy of the running wave increases due to the so-called skin effect-extrusion of
in the heat. In this case, the power balance equation the current to the outer surface (Fig. 1). As a result,
(3.1) for an electromagnetic wave propagating along the current does not flow through the entire section of
an infinitely long line at a distance x from its origin the wire, but only through a narrow outer layer (Fig. 1,
(Fig. 1) takes the form: whose thickness is determined by formula
P( x)  PG  PR ( x)  PS ( x) , (2)
σ
wherein PG – a wave power supplied to the line δ( F )  0,5 мм , (4)
from an external RF generator; PR – a wave power, μМ F
converting to heat due to of distributed long line of where δ - thickness of the skin layer, , М - rela-
resistance; PS – power of the line radiation. tive magnetic permeability F - frequency in MHz
As amplitude of the incident wave, that propa- The thickness of the skin layer decreases with in-
gates along the line, decreases exponentially, the sec- creasing frequency to several microns, which prede-
ond member of sum in (2): termines a sharp increase, by two orders of magnitude
PR ( x)  PG e2αx ,
or more, of the resistance of the wires compared to the
(3) industrial frequency of 50 Hz, and consequently, a
where  = Rf /  (1/m) – an attenuation constant; more efficient conversion of electromagnetic energy
 – a characteristic impedance of the line; Rf – the into heat
resistance per length unit of one of two identical ca-
bles which value in the microwave range depends on
the surface a of thin layer of δ due to the skin effect,

Fig.3

Skin depth δ with increasing frequency is re- From the formulas (1) and (2) for an aluminum
duced to a few millimeters. This determines the sud- wire, we obtain:
den increase of two orders or more of wire resistance   0,08 1 F , мм , R f R0  6 r F Ом м,
in comparison with the power frequency of 50 Hz, and
hence more efficient conversion of electromagnetic
energy into heat. According to the latest formulas for an aluminum
Due to the skin effect, the layer is thinning for the wire at r = 10mm and a frequency F = 100 MHz, we
current flow (Fig. 1), which increases the conductor's get for the thickness of the skin layer δ = 8 μm and an
specific resistance, the value of which taking into ac- increase in the line resistance of the conductor com-
count (1) at (r / 2δ)> 10: pared to the resistance at a constant current of 600
R f  R0 r 2 δ F  Ом м times.
, (5) Input impedance at a distance l from the end of
where r is the wire radius (mm), F is the frequen- the line:
cy in MHz, R0 (Ohm / m) is the resistivity of the same
conductor to the direct current.
Ux Z n chγλ x   ρ sh γλ x 
 Rin ( x )  jX in ( x )
Z in ( z ) 
   ρ ρ chγλ x   Z n shγλ x 
I x
, (6)

where γ = α + j β is the propagation constant, x = α  R f ρ 1 м , (7)

l / λ is the normalized line length, β = 2π / λ is the
phase constant, λ is the wavelength propagating in the
line, α is the damping constant, equal to: where Rf (Ω / m) is the line resistance of one of
the two identical wires.
28 Norwegian Journal of development of the International Science No 22/2018
An example of calculation according to (3) input resistance with the following initial data: .
graphs of the active and reactive components of the = 600 Ω, λ = 3 m, ZH = 106 Ω are shown in Fig. 4

Fig.4

It follows from the graphs (Fig.2) that as the by converting the radiation energy into heat .To do
length of the line with losses increases with x1 , this, along the main power transmission line (Fig.2).
the active component of the input resistance Rin tends along which the electromagnetic wave propagates,
to the wave resistance ρ, and the reactive component with a certain step it is necessary to place short lengths
Xin - to zero of wire - vertical vibrators (position 2 in Fig.5, a), the
However, as shown by experimental studies, it is radiation power of each of which :
not possible to substantially heat the wires only due to
PV  0,5 R I 2  400 l λ  I 2 ,
2
the skin effect, since the heat transfer of the upper (8)
heated layer of the wire to the surrounding air by con- where  is the wavelength, l is the length of the
vection is great. The power of the convection heat vibrator, RΣ is the radiation resistance of the antenna.
flow QT »QM - the power of the heat flow heating the The realization of such heating of the line is ex-
wire (Fig. 3) plained in Fig.5
A more efficient heating of the wire due to the
traveling electromagnetic wave can be accomplished

Fig.5

At fig.5 a) line 1, with a short vibrator 2 attached When replacing a wire segment with a graffiti-
to it ; b) a resistor 3 of one kΩ, connected to a vibra- ceramic resistor with resistance RRΣ (position 3 in
tor ; c) a resistor attached to the wire . Fig.5 b), most of the signal power will not be emitted
Norwegian Journal of development of the International Science No 22/2018 29
into the surrounding space, but dissipated as heat in RF-heating spots, placed with a certain step on the
this resistor, whose temperature is set to T0. Finally, wires, heat the latter and ice on them is not formed.
by transferring this resistor to the horizontal position When replacing a wire segment with a graffiti-
and leaning it against the main transmission line (Fig. ceramic resistor with resistance RRΣ (position 3 in
1c), it is possible to effect an effective line heating at Fig.5, b), most of the signal power will not be emitted
the expense of radiated power, the total value of which into the surrounding space, but dissipated as heat in
at a distance x is: this resistor whose temperature is set to T0. Finally,
by transferring this resistor to the horizontal position
N and leaning it against the main transmission line (Fig.

PS ( x)  Pi S ,
i 1
(9) 1c), it is possible to effect an effective line heating at
the expense of radiated power, the total value of which
where N is the total number of vibrators in the at a distance x is:
N


main line section
With the help of antennas - vibrators - pieces of PS ( x )  Pi S , (10)
i 1
wire length of 15-20 cm - forcibly increases the radi-
ated RF power, which is then converted into heat by
means of attached resistors to the antenna in one where N is the total number of vibrators in the
kOhm (Fig.5, b). The latter, in turn, are superimposed main line section. .
on the wire, which leads to its heating (Fig. 5c). Such A sketch of the construction of the HF antenna-
heater is shown in Fig. 6

Fig.6

The following designations are accepted in Fig. ry of the firm "Engineering Center ORGRES" in
6: 1 - stranded wire AC-300 (diameter 24 mm) for Khotkovo, Moscow Region .
power lines, 2 - insulating, heat-conducting layer, 3 - As a high-frequency generator with a frequency
graphite layer with resistance 1 kΩ, power dissipation of 81.36 MHz (According to the "Radio Regulations,
up to 20 W, temperature 100-1200С, 4 thermal insula- this frequency is intended for industrial applications)
tor 1.15 kW power was used by the radio transmitter"
4. Experimental verification of the method of Rime ", manufactured by MART (St. Petersburg)
heating the wire by means of a high-frequency elec- (Fig.7)
tromagnetic wave.
The experimental verification of the method of
heating aluminum wires was carried out on the territo-

Figс.7 Fig.8

To match the output resistance of the RF genera- stretched wires type AS-120 - every three meters were
tor Rgen = 50 Ohm with the input resistance of the fixed devices for converting RF radiated energy into
two-wire line Rlin ~ ρ = 600 Ohm, an HF transformer, heat, made on the basis of resistors of the type MLT-2.
matching device was made (Fig. 8). Along both lines - Each transducer consisted of ten resistors with a total
30 Norwegian Journal of development of the International Science No 22/2018
allowable dissipation power of 20 watts. With the RF points of connection of the heaters exceeded the ambi-
generator turned on, the heating of the resistors ent air temperature by 100C
reached 80-1000 C. The graph of temperature variation
along the wire is shown in Fig. 9. At the same time,
the wire temperature at points 1 meter away from the

Fig.9

The main result of experimental studies: - RF generators with a frequency of 81.36 MHz
1. To heat the AC-300 wire per 100С, a specific with a power of 12 kW (1),
power of 20 W / m is required. - matching device for communication of the RF
2. We assume that the distance of 200 m between generator with a three-phase power line (2);
the transmission line supports is 200 m. Therefore, for - block of automatic adjustment of the matching
heating the 400m three-phase transmission line, two device (3),
RF generators with a power of 12 kW are required. - converters-heaters, attached to wires (4);
5. About creation of industrial installation - high-voltage vacuum capacitors for connection
The general scheme of the installation includes to power lines (5) .
(Fig. 10, 11)

It is assumed remote, automatic control system The conclusion.

"Anti-ice" using a satellite system of the type "Tura- 1. Within the framework of the theory of the
ya" or fiber-optic telecommunication system electromagnetic field, the principles of constructing a
system for heating the wires of power lines are con-
Norwegian Journal of development of the International Science No 22/2018 31
sidered in order to prevent the formation of ice on REFERENES:
them. 1. 1.Патент № 2356148 от 15.05.2008, Способ и
2. With the help of the new method it is proposed устройство для борьбы с гололедом на линиях
not to fight ice, but to prevent the formation of ice on электропередачи. МИРЭА и Каганов В.И.
the wires. 2. Jackson, J. Classical Electrodynamics. - Mos-
3. The carried out experimental studies confirm cow: Mir, 1965
the possibility of creating a similar system, the ad- 3. Tamm IE Fundamentals of the theory of elec-
vantage of which are tricity. - Moscow: GITTL, 1956.19.
Prevent the formation of ice on the wires, 4. W / Kaganov, S.Tereshchenko. Microwave
- disconnection of disconnection from the con- electromagnetic wave - promising projects. - Lambert,
sumer network is not performed, 2017.
- absence of burning wires during 5. Kaganov V.I. Radio engineering circuits and
- Saving of electric energy. signals. 4 th edition, -M .: Fo-rum -INFRA, 2018
- the absence of the so-called "dance" wires, 6. Kaganov V.I. Transistor radio transmitters.-M
-Exception of technological short circuit of wires .: Energia, 1976.

PRINCIPLES OF DESIGNING MOBILE ROBOTS

Polishchuk M.
Candidate of Technical Sciences, Associate Professor
National Technical University of Ukraine
“Igor Sikorsky Kyiv Polytechnic Institute”, Kyiv, Ukraine

Abstract
The article proposes three basic principles for the design of mobile robots: the accumulation and conversion
of energy, the integration of motion drives and the use of the generator of reactive pneumatic traction. These
technical solutions allow robots to navigate on surfaces of arbitrary orientation in different coordinate systems.
The implementation of these principles contributes to reducing the total power of the drives and increasing the
reliability of holding robots on the surface of arbitrary orientation in the technological space. The results of
mathematical modelling of constructive and technological parameters of mobile robots are described.

Keywords: mobile robots, walking mechanisms, robots of vertical movement, climbing robots

Introduction the extreme conditions of man-made disasters that are

The evolution of technical systems in the field of dangerous and even unacceptable for human presence.
engineering science, as a rule, led to the emergence of 1. Prerequisites and means for solving the
new means of production. In the twentieth century problem
scientific and technical thought created sufficiently The problem of creating this type of robots is the
reliable means of overcoming the gravitational force lack of a methodology for synthesizing the robot sub-
in the form of flying and reactive equipment. These systems that compensate or overcome the gravitational
funds are used as a reliable transport. However, up to load in order to ensure that the robot is held on an ar-
the present time there are no industrial samples of bitrarily oriented surface during the execution of tech-
equipment for performing contact technological opera- nological operations. Therefore, the purpose of this
tions while simultaneously overcoming the forces of research is to develop principles for the synthesis of
gravity. mobile robots and the introduction of modern comput-
Mobile robots of arbitrary orientation are known er modeling tools for the transition to the design and
in publications as robots of vertical movement or in manufacture of prototypes of mobile robots.
international publications – under the term Climber Theoretical and experimental studies on the crea-
Robot, are a new modification of mobile robots. These tion of robots of arbitrary orientation in the technolog-
robots are equipped with means of fixation on the sur- ical space began in the last decade of the twentieth
face of arbitrary orientation relative to the horizon of century in the countries of Western Europe, Japan, the
the technological space. The creation of this kind of United States, Korea, China and Russia. To date, there
robotics is at the initial stage and is dictated by the are mainly experimental samples of such robots.
need to perform technological operations in such areas Mobile works [2, p. 24; 3, p. 28–32; 4, p. 1; 5, p.
of industrial activity as monitoring of industrial facili- 60] are equipped with devices to fix the robot on sur-
ties, installation and dismantling of building struc- faces of arbitrary orientation, and in studies [6, p. 70–
tures, repair and prevention of their components. In 78; 7, p. 22–28] hybrid drives are proposed that can
the context of the fourth industrial revolution, "Indus- improve the energy efficiency of mobile robots. Tech-
try 4.0" [1, p. 69], with the focus on the use of robotic nical solutions [8, p. 4; 9, p. 6; 10, p. 348; 11, p.54]
systems, information and communication tools, and allow the robot to move on surfaces oriented at differ-
the use of these robots becomes especially urgent in ent angles to the horizon, but only in 2D space, that is,
in a plane. In the general case, the variations of the
32 Norwegian Journal of development of the International Science No 22/2018
constructions of the above-mentioned robots restrict 2. Solution of the problem under consideration
their movement only in the Cartesian coordinate sys- In contrast to the above, the concept of improving
tem. robots of arbitrary orientation based on three funda-
Unlike the aforementioned technical solutions, mental principles of the synthesis of mobile robot de-
the robot model [12, p. 42] allows servicing objects in signs is proposed:
the cylindrical coordinate system, in particular, objects 1) Accumulation of potential energy at each pre-
such as trees, but with a soft porous surface for the vious step of the robot's movement and subsequent
movement of the robot, which limits the technological conversion to kinetic energy at the next step of the
capabilities of mobile robots. At a time when there are movement.
objects facilities that are also close to the cylindrical 2) Integration of drives of longitudinal, vertical
coordinate system, for example, electric grid posts, displacement, and also drives of change of orientation
columns, pipes of thermal power plants and the like. of the robot according to the set route.
In addition, the mobile robot should also work in a 3) The use of traction generators (aerodynamic
system of angular coordinates, which is typical for lift forces) as a means of counteracting the gravita-
humans. The development of systems for connecting tional force to increase the technological load while
the robot to the surface of motion is a technical solu- simultaneously reducing the power of clutch drives
tion [13, p.4], which uses adhesion technology. How- and the movement of robots.
ever, current implementations of this technology are As technical means of implementing these prin-
characterized by an extremely low speed of movement ciples of synthesis, consider the corresponding models
of the robot due to the effect of slow adhesion. This of robots of arbitrary orientation. In Fig. 1 shows a
property still prevents the industrial use of adhesion as mobile robot [14, p. 2] realizing the first of the above
a method of fixing the robot to a surface of arbitrary mentioned principles, namely having the ability to
orientation. Thus, the problem of synthesizing mobile accumulate potential energy at each previous step and
robots capable of performing technological operations converting it into kinetic energy of motion at each
in a space of arbitrary orientation is topical. subsequent step of displacement.

Fig. 1. The robot model of arbitrary orientation with the accumulation and transformation of energy

On the robot body 1 are mounted rotary pneumat- The kinetic energy T1 of the robot's foot free from
ic actuators 2 connected through a gear train 10 with adhesion to the moving surface will be
running mechanisms 3. When the elastic members 4 1
T1   (Vx1  Vy1 )dm ,
2 2
are compressed due to the rotation of the legs of the (1)
robot, the potential energy is accumulated during the 2
first half of the step, and in the second half of the step, where: Vx1, Vy1 – the velocity projections on the
these elements convert the potential compression en- coordinate axis 0xyz; m – the mass of the robot.
ergy in the kinetic energy of the robot's motion [15, p. Pedipulator of the robot, whose grips are engaged
117; 16, p. 10]. To implement a plane-parallel motion with the surface, performs a rotary motion with angu-
along the vertical, the robot is equipped with parallelo- lar velocity 1 . The expression for the kinetic energy
grams 5. The latter connect pedipulators with vacuum of a supporting leg of mass m1 and radius R can be
grippers 6 that hold the robot on a surface of arbitrary obtained from expression (1) after integrating it, sub-
orientation. Also, the robot is equipped with a control stituting the velocity of the translational motion of the
unit 7 and a power supply module 8. To overcome robot V = 0 and the angular velocities of the pedipula-
obstacles on the moving surface, the robot is addition-
tors  2  1 :
ally equipped with rolling bearings 9.
Volumes of the potential and kinetic energy of 2
m1 R
T2  ( 1 ) .
2
the robot motion at different stages of displacement (2)
can be determined using Lagrangian equations of the 6
second kind. Thus, the movement of the robot at each second
half of the cycle occurs due to the energy accumulated
Norwegian Journal of development of the International Science No 22/2018 33
at each first half of the travel step. This allows 40% ... The main characteristic of elastic elements is
45% to reduce energy costs for the movement of the their rigidity j – parameter, which determines the force
robot and to direct the resulting energy reserve for the of compression of these elements, and hence the value
execution of technological operations. Pedipulators of the accumulated potential energy in the first half of
work with elastic elements (Fig.1, pos.4), which per- the step of the pedipulator.
form the function of accumulation of potential energy
during the first half of the cycle of displacement.

Fig. 2. The effect of the elastic element j (N/m) on the work performed in the step β1> 45o of the robot movement

In Fig. 2 shows the dependence of the variation tional movements and three rotational, each of which
of work "A" on the stiffness of the elastic element j according to the classical solutions corresponds to an
(N/m) and the forces of weight in the second stage β1> autonomous drive. The method of Fig. 3 – technical
45o of displacement, that is, during the transformation implementation of this principle eliminates the need
of the potential energy into the kinetic energy of the for drives for each of the coordinate axes. To do this,
robot's motion. At the second stage of displacement the robot is equipped with flexible running mecha-
nisms 2 mounted on the body 1. Each pair of legs of
45  1,2  90 drive the of pedipulator is turned
0 0

the robot through the transmissions 3 is provided with

off to save energy resources of the robot, and it moves electric drives 4. The grippers 5 keep the robot on the
only due to kinetic energy. Analysis of these graphs surface moving, and the rotary actuators 6 change the
shows that to increase the kinetic energy of the robot position of the grippers relative to the displacement
movement, it is advisable to increase the rigidity, de- surface. The robot platform has a power supply unit 7,
spite the fact that in this case the counteraction to the a hydraulic or pneumatic valve unit 8 and a gas or
drive increases in the first half of the step, that is, the liquid pressure generator and a controller 9 for con-
efficiency of the drive decreases. However, this nega- trolling the robot. Due to the fact that each foot of the
tive manifestation can be compensated for by an in- robot is made in the form of a compressed set of hem-
crease in the transmission ratio (see item 10, Fig. 1) of ispherical rings inside which corrugated pipes are
pedipulators. placed under different pressures, the robot has the
The second principle, as noted above, involves ability to work in different coordinate systems: rec-
the integration of displacement drives [17, p. 6] with tangular Cartesian, spherical and cylindrical without
the aim of reducing them, and hence of reducing the additional drives on each axis of coordinates.
mass of the robot. It is known that in the Cartesian
space we have six degrees of freedom – three transla-
34 Norwegian Journal of development of the International Science No 22/2018

Fig. 3. Model of the robot with flexible pedipulators

In each leg of the robot there are four corrugated to the displacement surface and the permissible tech-
pipes. Two pipelines with pressure p1 p2 are placed in nological load to ensure the reliability of its industrial
a vertical plane and two other pipelines – in a horizon- operation. Having formulated the system of equilibri-
tal plane with pressure p3 and p4. Due to the action of um equations on the basis of the Herman–Euler–
these pressures, forces appear that, while bending the D`Alambert principle (here we omit the record for
robot's leg, orient the robot in the technological space: brevity), we find the corresponding reaction forces
d2 d2 d2  d 2 (3) N1,2 and the frictional forces Q1,2 by the robot foot to
F1  p1 ; F2  p2 ; F3  p3 ; F4  p4 , the displacement surface and then compare them with
4 4 4 4
where d – the internal diameter of the corrugated the technological load N in depending on the angle α
tubes. Since the axes of the pipelines are offset from of the robot inclination to the horizon.
the pedipulator axis by the eccentricity e, there are mo- The reaction forces of N2 and the frictional forces
ments M1 and M2 that flex the robot's leg: Q2y of the robot supports with the displacement sur-
face are determined as follows (the designation of the
d d
2 2
M1  ( p1  p2 )e; M 2 
( p3  p4 )e (4) parameters, see Fig.3):
4 4 N2  Q2  a3G  b3 N ; Q2 y  d3G  h3 N ,(5)
where: e – eccentricity of placement of corrugat-
where for the compactness of the incoming val-
ed pipelines in the plane of the coordinate system. ues of variables is denoted:
To develop a robot, it is necessary to establish a
connection between the forces of adhesion of its legs
a3  d1b2  b1d 2  ; b3  b2 h1  b1h2  ; d1  yc cos   zc sin  ;
b2  x2 cos(   ); d3  d2 a1  d1a2  ; h3  h2 a1  a2 h  ;
a1  y2 cos(   )  z2 sin(   ); a2   x2 sin(   );
d2  xc sin  ; b1  y2 sin(   )  z2 cos(   ); h1  yc cos  zc sin ; h2  xc sin ;   a1b2  b1a2 ;
x2, y2, z2 – coordinates of the contact point of the the same system of equilibrium equations, we find the
second leg of the robot with the displacement surface; remaining unknown reactions N1 and the frictional
xc, yc, zc are the coordinates of the center of gravity of forces Q1y:
the robot; α, φ – angles of inclination to the horizon of N1  Q1  Ga4  Nh4 ; Q1 y  Gh5  Nh6 , (6)
surfaces on which the robot's legs rest; ψ is the angle
of inclination of the central axis of the robot passing where also for the compactness of writing varia-
through its center of gravity G (see Fig.3). Then, from bles is defined:
a4  cos   a3 cos(   )  d3 sin(   ); h4  cos  b3 cos(   )  h3 sin(   );
h5  sin   d3 cos(   )  a3 sin(   ); h6  sin  h3 cos(   )  b3 sin(   ).
For the stability of the robot, the frictional forces where μ – the coefficient of friction of the grip of
of each of its legs must not exceed the boundary val- the robot's leg with the surface along which the robot
ues: moves. Substituting in expression (7) the expressions
Q1y   N1; Q2 y   N2 ; N1  0; N2  0, (7) above the found reactions of forces (5) and (6), we
find limitations for the technological load of the robot
taking into account the forces acting on it:
Norwegian Journal of development of the International Science No 22/2018 35
Q1  Ga4 Q  Ga3
N1  0  N  ; N2  0  N  2 ;
h4 b3
(8)
Q2  G( a3  d3 ) Q1  G(  a4  h5 )
N ; N .
h3  b3 h6   h4
Among the values of the reaction N of the tech- legs of the robot, and curves 3 and 4 are the beginning
nological load calculated in accordance with condi- of slippage of the said robot legs, respectively.
tions (8), we choose the largest, which simultaneously The robot of arbitrary orientation can be on such
satisfies all inequalities, which allows to determine the surfaces as a floor, a wall or a ceiling. Accordingly, if
maximum technological load of the robot, for exam- the robot is on the floor, then in the above dependenc-
ple, for drilling, rivets, dowels, etc. technological op- es it is necessary to substitute the value
erations.   0о ,   0о ; if on a vertical wall, then
As a result of the simulation, the limiting values
of the technological load are obtained (Fig.4): curves 1   90о ,   90о and if on the ceiling, then
and 2, respectively, determine the separation states   180о ,   180о and so on.
from the displacement surface of the first and second

Fig. 4. Graphs of the boundary values of the technological force N of the robot as a function of the angle α of the
inclination to the horizon

And, finally, the third principle – the use of trac- in Fig. 5. Like the previous one, it also has flexible
tion generators as a means of counteracting the gravi- pedipulators 2, grippers 3, gear 4 and electric drives 5
tational force is realized by the robot [18, p.6], shown on the body 1.

Fig. 5. Mobile robot with pneumatic generator of aerodynamic lift

The main difference of this robot is the installa- force G, regardless of the position of the robot in the
tion in the center of its masses suspension Cardan 6 XYZ space. This principle allows us to differentiate
with three degrees of freedom and a pneumatic gener- the approach to regulating the lifting force of the ro-
ator of traction 7. The location of the thrust generator bot, depending on its orientation in space.
on the Cardan suspension allows the thrust generator Similarly to the previous case, from the system of
to maintain the coincidence of the lines of action of equilibrium equations of the robot determine the forc-
opposing forces: the rise of G1 and the gravitational
36 Norwegian Journal of development of the International Science No 22/2018
es of normal reactions Ni of the legs of the robot and The frictional force Q1y and the normal reaction N1
the corresponding frictional forces Qi (see Fig.5): are determined as
N 2  b12 (G  G1 ) , (9) 1 1
Q1 y  (G  G1 )( b12  sin  ); N1  (G  G1 )(b12  cos  ).
where: G – the weight of the robot; G1 – traction 2 2
force; α - angle of inclination of the plane of move-
ment of the robot to the hori- As can be seen from the graphs of Fig. 6 with
zon b12  cos  yc  sin  zc 2 y2 ; yc, zc – coordinates positive reactions N1 and N2 with the angle of inclina-
tion of the displacement surface α ≤ 540, the weight of
of the center of gravity of the robot; y2 – coordinate of the robot increases the technological load.
contact with the moving surface of the second part of
the robot.

Fig. 6. Graphs of the dependence of the normal reactions N1, N2 of the robot on the angle α of its inclination to
the horizon

This means that the inclusion of the jet engine is sion of technological capabilities of mobile robots of
more suitable for the values of the angle of inclination arbitrary orientation in the technological space.
of the robot to the horizon at α ≥ 540. Of course, the 4. Conclusion
critical angle of inclination depends on the other cen- The concept of the synthesis of mobile robots is
trifugal characteristics of the robot. However, using based on three principles of constructing their designs.
the generator of aerodynamic forces ensures reliable These principles can be used autonomously and in
retention of the robot on surfaces of arbitrary orienta- combination, depending on the technological purpose
tion with any structural and technological parameters of the robot and its profitability, which, in turn, is de-
of the robots. To calculate the critical angle of inclina- termined by the area of industrial operation. The pro-
tion to the horizon and, accordingly, to regulate the posed approach to the synthesis of mobile robots can
force of the reactive thrust, that is, the aerodynamic reduce the weight of structures by reducing the total
force of the pneumatic generator, it is necessary to number of drives. This increases the energy resource
take into account other characteristics of the robot. for improving the efficiency of both transport and
3. Results and discussion technological operations performed by the robot in
Creation of means for accumulation of potential various areas of industry. The next stage of the re-
energy of drives with the subsequent transformation search is the experimental approbation of this concept
into kinetic energy of movement of the robot, and also of the synthesis of mobile robots.
integration of drives allows reducing considerably the
general power that provides reduction of gravitational REFERENCES:
load.
The use of a pneumatic motor as a means of 1. Kuznetsov Yu.N. Challenges of the fourth in-
counteracting the gravitational force makes it possible, dustrial revolution "Industry 4.0" in front of scientists of
when adjusting the traction force, depending on the Ukraine. Vestnik KhNTU, Kherson. 2 (61), 2017, p. 67-
angle of the robot's inclination to the horizon, to in- 75.
crease the reliability of holding the robot on the sur- 2. Raju D.D., Jaju S.B. Developments in wall
face of an arbitrary orientation. climbing robots: a review. International journal of engi-
A fundamentally new realization of the legs of neering research and general science. 2 (3), 2014, 37p.
mobile robots in the form of a set of hemispherical 3. Gradetskii V.G., Veshnikov V.B., Kalinichen-
rings and corrugated pipelines at different gas or liq- ko S.V. Controlled movement of mobile robots over
uid pressures makes it possible to achieve an arbitrary arbitrarily oriented surfaces in space. M: Nauka, 2001,
orientation of the mobile robot in different working 360p.
spaces: a rectangular Cartesian, spherical and cylindri- 4. Chernous'ko F.L. Mobile robots: research, de-
cal coordinate system. This effect ensures the expan- velopment, prospects.
http://www.ras.ru/news/shownews.aspx?id=f5c75bcf-
Norwegian Journal of development of the International Science No 22/2018 37
2fa5-40e6-b067-4492c5ab22&print=1. (Chernous'ko 12. Tin Lun Lam., Yangsheng Xu Tree Climbing
F.L. Bolotnik N.N. and Gradetskii V.G.) Robot: Design, Kinematics and Motion Planning. New
5. Polishchuk M.N. Modernization of the vacuum York: Springer Heidelberg. 2012, p. 37–46.
gripper of the mobile walking robot. Mechanics and 13. Yehya M., Wasim A., Jahanzaib M., Abdalla
Advanced Technologies, vol. 2 (80), 2017, p. 59–64. H. A. Unipolar Electroadhesion Pad Technology for
6. Kuznjecov Ju.M., Shybecjkyj V.Ju. Robotic Adhesion Mechanism of Wall Climbing Robot.
systems complexes of pharmaceutical and biotechno- International Journal of Robotics and Mechatronics, vol. 2 (1),
logical productions. Kiev, GhNOZIS, 2012, 335 p. 2014, p. 1–10.
7. Pavlenko I.I. Industrial work: calculation ba- 14. Yampolskiy L.S. Method and device for
sics and design. KNTU. Kirovoghrad, 2007, 420p. movement of pedipulators of walking robot. State
8. Sinev A.V., Gradetskii V.G., Kravchuk L.N. The Register of Patents of Ukraine, Pat. 111021, 2016,
vehicle for moving along the surfaces arbitrarily ori- 11p. (Yampolskiy L.S., Polishchuk M.N., Persikov V.
ented in space. State Register of Patents of Russia, K.).
Pat. № 2262461, 2005, 6p. 15. Polishchuk M., Oliinyk V. Mobile climbing
9. Baranov V. E. Vehicle for moving on inclined robot with elastic energy accumulators. Mechanics
and vertical surfaces. State Register of Patents of Rus- and Advanced Technologies, vol. 1 (82), 2018, p. 116
sia Pat. № 2267434, 2006, 8p. (Baranov V. E., Vagin – 122.
S.G., Zatsepin V. A., Kokorin A. G.) 16. Kuznetsov Yu.M., Polishchuk M.M. Walking
10. Chashchukhin V.G. Investigation of the mobile robot Kuznetsova – Polishchuka. Ukrainian
movement parameters of a robot with a sliding seal.– Institute of Intellectual Property, Ukraine, application
Vestnik Nizhegorodskogo universiteta im. N.I. for a patent № а2018079762018, 16p.
Lobachevskogo, vol 4 (2), 2011, pp. 347–349. 17. Polishchuk M., Opashnianskyi M., Suyazov
11. Shamim Hasan. Smartphone Controlled Spy N. Walking Mobile Robot of Arbitrary Orienta-
Robot with Video Transmission and Object Collector, tion, International Journal of Engineering and Manu-
International Journal of Engineering and Manufactur- facturing (IJEM), vol.8, No.3, 2018, pp.1–11.
ing (IJEM), Vol.7, No.6, 2017, p.50-58. (Shamim Ha- 18. Polishchuk M.N. Antigravity mobile robot
san, Khalid Hossain Jewel, Niaz Mostakim, Nabil Polishchuk. Ukrainian Institute of Intellectual Property,
Hossain Bhuiyan, M.K. Rahman, Sheikh Dobir Hoss- application for a patent № а201805661, 2018, p. 8–10.
ain and Khalid Hossain).

MATHEMATICAL MODELS IDENTIFICATION IN CHARACTERISTIC SITUATIONS OF CONDUCT

OF THE BREWERY TECHNOLOGICAL COMPLEX

Romanov N.
PhD, National University of Food Technologies, Kiev
Svitoch I.
PhD student, Igor Sikorsky Kyiv Polytechnic Institute, Kiev
Macebula D.
PhD student National University of Food Technologies, Kiev

ИДЕНТИФИКАЦИЯ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ В ХАРАКТЕРНЫХ СИТУАЦИЯХ

ПОВЕДЕНИЯ КОМПЛЕКСА ПИВЗАВОДА

Романов Н.
Кандидат технических наук, старший преподаватель,
Национальный университет пищевых технологий, г. Киев
Свиточ И.
Аспирант, Киевский политехнический институт имени Игоря Сикорського, г. Киев
Мацебула Д.
Аспирант, Национальный университет пищевых технологий, г. Киев

Abstract
Goal. The purpose is to create mathematical models that allow to describe such basic qualitative
characteristics of brewing process as quality of production at all stages of the production process, productivity
and losses.
Method. To pass the passive experiment in the production conditions of the brewery according to the results
of which (data of measurements of technological parameters and data of laboratory analyzes) to carry out the
identification of mathematical models. In the Matlab, using the least squares method, evaluate the competitive
models of certain structures and the parameters of these models based on the criteria: determination, Akaike
information criterion, Hawking and Fisher statistics.
Results. The use of the above criteria of adequacy allowed us to select the optimal model structures for each
situational-significant area in terms of their accuracy and complexity.
38 Norwegian Journal of development of the International Science No 22/2018
Scientific novelty. On the basis of the conducted research, models were developed describing the
dependencies of the output process variables, which were the main qualitative indicators that can be quickly
measured by objective methods and definitely characterize the quality of semi-finished products and finished
products (malt extraction, concentration of hopped wort, stage of fermentation), equipment productivity, specific
losses of resources at each stage of production depending on the regime parameters and input products. Models
were determined with fuzzy values of input variables.
Practical significance. The created models can be used to create new or to improve the existing systems of
automation of brewing processes.
Аннотация
Цель. Цель работы заключается в создании математических моделей, позволяющих описать следу-
ющие основные качественные характеристики процесса пивоварения как качество продукции на всех
стадиях производственного процесса, производительность и потери.
Методика. Провести пассивный эксперимент в производственных условиях пивзавода по результа-
там которого (данные измерений технологических параметров и данные лабораторных анализов) осуще-
ствить идентификацию математических моделей. В среде Matlab методом наименьших квадратов, оце-
нить модели-конкуренты определенных структур и параметры этих моделей на основе критериев: детер-
минации, информационного критерия Акаике, статистик Хокинга и Фишера.
Результаты. Использование приведенных критериев адекватности позволило осуществить выбор оп-
тимальных структур моделей для каждой ситуационно-значимой зоны с точки зрения их точности и
сложности.
Научная новизна. На основании проведенных исследований созданы модели, описывающие зависи-
мости выходных переменных процесса, которыми были основные качественные показатели, которые
могут оперативно измеряться объективными методами и изначально характеризовать качество полуфаб-
рикатов и готовой продукции (экстрактивность солода, концентрация охмеленного сусла, степень сбра-
живания), производительность оборудования , удельные потери ресурсов на каждой стадии производства
в зависимости от режимных параметров и входной продукции. Модели определялись при нечетких зна-
чениях входных переменных.
Практическая значимость. Созданные модели могут быть использованы для создания новых или со-
вершенствования существующих систем автоматизации технологических процессов пивоварения.

Keywords: brewing, fuzzy values of variables, identification of mathematical models, situational-

significant zones.
Ключевые слова: пивоварение, нечеткие значения переменных, идентификация математических
моделей, ситуационно-значимые зоны.

Пивоваренная промышленность является од- Несмотря на большое качественное и количе-

ной из ведущих отраслей пищевой промышленно- ственное разнообразие отдельных технологиче-
сти Украины. Основными направлениями повы- ских процессов, их различные масштабы, режимы
шения эффективности пивоваренного производ- и условия реализации, все они имеют общие свой-
ства являются: улучшение качества продукции, ства: они структурно сложны, состоят из отдель-
повышение производительности технологического ных частей или элементов; перерабатывают веще-
оборудования и рациональное использование ство, энергию или информацию; связаны с други-
энергетических и материальных ресурсов. Техно- ми технологическими или сопутствующими
логические процессы в пивоварении характеризу- процессами.
ются, как правило, большой сложностью и разно- Анализ работ по технологии и автоматизации
образием. Общей характерной чертой всех суще- пивоваренных производств показал, что созданы
ствующих технологических процессов является то, необходимые технологические условия для обес-
что для преобразований исходного сырья, матери- печения эффективного автоматизированного
ала или полуфабриката в определенный конечный управления производственными процессами пиво-
продукт, необходимо большое количество функ- варения. Среди таких работ можно выделить рабо-
ционально отличных степеней переработки. ты ученых В. Кунце, Л. Нарцисса, В.А. Домарец-
Для целенаправленного течения этих кого. В области автоматизации производственных
процессов на отдельных ступенях или на их процессов неоценимым вкладом являются работы
конкретных стадиях необходимы различные виды ученых: П. Телиса, Л.Третяк, Е. Егорова. Сетецен-
энергии, вспомогательных веществ и информации. трический подход в области управления развит в
Условия и режимы реализации отдельных стадий работах С.А. Юдицкого, Ю.С. Затуливетра, С.А.
технологических процессов могут также Куджа [2,3]. Использование сценарного подхода
значительно отличаться (щелочные, кислотные, при управлении сложными объектами обосновано
маслянистые и другие среды; высокие в работах В.В. Кульбы, С.А. Косяченко,
температуры или низкие, высокое давление и Д.А.Кононова [4]. Обобщая, можно отметить, что
вакуум; статическая нагрузка или вибрация; отсутствует комплексный подход к проблемам
мощные электромагнитные поля и т.д.) [1]. автоматизированного управления технологиче-
ским комплексом пивзавода в целом. На основе
Norwegian Journal of development of the International Science No 22/2018 39
анализа принципов управления сложными техно- продукции. Модели определялись при нечетких
логическими комплексами и достижений совре- значениях входных переменных.
менной теории и практики автоматизированного Обработка экспериментальных данных осу-
управления определены актуальные проблемы ществлялась методом наименьших квадратов, ко-
систем управления пивзаводом. торый позволяет по экспериментальным данным
Необходимость планирования на предприятии подобрать такую аналитическую функцию, кото-
в условиях высокой изменчивости окружающей рая подходит настолько близко к эксперименталь-
среды связана с неопределенностью, риском, а ным точкам, насколько это возможно [5].
значит и с появлением проблемных ситуаций на Идея метода наименьших квадратов заключа-
предприятии. Причиной появления проблемных ется в том, что функцию:
или кризисных ситуаций может быть не только Y  f ( x, a0 , a1 ,...ak ) (1)
объект, но и субъект управления.
Возникновение различных ситуаций, опасных необходимо подобрать таким образом, чтобы
для производства и в целом для предприятия, объ- сумма квадратов отклонений измеренных значе-
ективно характеризуется большим количеством ний уі от расчетных Yi была наименьшей:
n
взаимосвязанных факторов, повышающих слож-
ность управления. Практика показывает, что про- 
S  [ y i  f ( xi , a0 , a1 ,...ak )]2  min (2)
i 1
блемные ситуации возникают и протекают по- Задача сводится к определению коэффициен-
разному. Ни одна производственная система не тов а из условия (2). Для ее решения необходимо
і
функционирует в вакууме, ее окружает очень раз- составить систему уравнений:
нородная среда. Все, что находится вне предприя-
тия, в той или иной степени постоянно влияет на S 
 0,
него. Это влияние протекает во времени неравно- a0
мерно или периодически с разной силой и прояв- 
S
ляется неодинаково. Предприятие должно хорошо  0, 

знать не только среду, в которой находится, но и a1  (3)
уметь адекватно реагировать на каждое из внеш- ... 
них воздействий. S 
Для получения математических моделей для  0.
оценки критериев оптимальности в ситуационно- ak 

значимых зонах в производственных условиях Дальнейший расчет выполняется на основе
пивзавода ПАО «Оболонь» был проведен пассив- алгоритма вычисления среднеквадратичного от-
ный эксперимент, в ходе которого осуществля- клонения, а также других критериев.
лась, с интервалом обеспечивающим статистиче- При постановке задачи предполагается, что
скую независимость измерений, фиксация контро- структура математической модели известна, по
лируемых переменных, характеризующих теоретическим основам и основным гипотезам о
технологические процессы приготовления пива механизмах явлений, которые происходят в техно-
(общее количество переменных 31). логических процессах пивоваренной промышлен-
По результатам эксперимента (данные изме- ности [6,7].
рений технологических параметров и данные ла- Метод наименьших квадратов реализован с
бораторных анализов) осуществлялась идентифи- помощью программного средства Matlab
кация математических моделей методом наимень- Optimization Toolbox. В зависимости от того как
ших квадратов в среде Matlab, где оценивались параметры входят в аналитическую зависимость
модели-конкуренты определенных структур и па- (4), мы получаем или линейную, или нелинейную
раметры этих моделей на основе критериев: де- функцию.
терминации, информационного критерия Акаике, Полином имеет следующий вид:
статистик Хокинга и Фишера. Использование при- Y  AX1  BX 2  CX1 X 2  DX1  EX 2  F (4)
2 2

веденных критериев адекватности позволило осу- де A, B, C, D, E, F – искомые коэффициенты

ществить выбор оптимальных структур моделей полинома.
для каждой ситуационно-значимой зоны с точки В качестве примера, рассмотрим построение
зрения их точности и сложности. модели для некоторых основных контролируемых
Модели описывали зависимости выходных переменных процессов пивоварения.
переменных процесса, коими были основные каче- Модель качественного показателя процесса
ственные показатели, которые могут оперативно получения товарного солода.
измеряться объективными методами и изначально Данная модель соответствует области изме-
характеризовать качество полуфабрикатов и гото- нения параметра нечеткого понятия "Экстрактив-
вой продукции (экстрактивность солода Е, кон- ность" "высокая". При этом, число Кольбаха "вы-
центрация охмеленного сусла Cохм , степень сбра- сокое", содержание редуцирующих сахаров "нор-
живания S збр ), производительность оборудования, ма", время осахаривания "норма" влагосодержание
солода "норма", содержание азотистых веществ
удельные потери ресурсов на каждой стадии про- "норма", температура сушения "высокая". Полу-
изводства от режимных параметров и входной чили поленом следующего вида:
40 Norwegian Journal of development of the International Science No 22/2018
E ( K , С ред ,  оц , t суш , Wсол )  0,629 K  0,040083 С ред  0,2004417  оц  0,675639 t суш 
 0,040861Wсол  0,01366 K 2  0,00174 KС ред  0,0087 K оц  0,00055 С ред оц 
 0,00138 оц 2  0,0000556 С ред2  0,4249 Kt суш  0,0257 KWсол  0,02708 С редt суш  (5)
 0,00163 С редWсол  0,1354 t сушWсол  0,00819 оцWсол  0,00156 t суш  0,0000057 Wсол 
2 2

 0,00019 t сушWсол  5.383

где E - экстрактивность, %; K - число Коль- Решение уравнения регрессии приведены на
баха, %; С ред - содержание редуцирующих саха- рис. 1, из которого видно, что при росте темпера-
туры сушения и числа Кольбаха экстрактивность
ров, г/100г сухого вещества;  оц - время осахари- увеличивается и соответственно выход экстракта
вания, мин.; t суш - температура сушения, °С; Wсол из солода увеличится, что приводит к повышению
- влагосодержание солода, %; производительности участка.

Рис. 1. Изменение экстрактивности товарного солода при нечетких значениях переменных число
Кольбаха "высокое", содержание редуцирующих сахаров "норма", время осахаривания "норма"
влагосодержание солода "норма", содержание азотистых веществ "норма", температура сушения
"высокая".

Модель производительности процесса полу- та "меньше нормы", температура сушения "нор-

чения товарного солода. ма", общая кислотность "выше нормы", содержа-
Изменение параметра нечеткого понятия вы- ние редуцирующих сахаров "выше нормы" проис-
ход экстракта "высокий" при содержание экстрак- ходит по модели вида:
Вихекст р( Вм екст р, t суш , К заг , С ред )  0,355638 Вм екст р  1,942362 t суш  0,12312 К заг 
 4,75976 С ред  0,01066 Вм екст р2  0,05824 t суш 2  0,01803 Вм екст рt суш 
(6)
 0,000508 Вм екст рК заг  0,019646 Вм екст рС ред  0,002776 t суш К заг  0,1073 t сушС ред 
 0,00099 К заг2  0,03825 С ред2  0,00193 К загС ред  63,24
где Вихекстр- выход экстракта, %; Вмекстр- Решение данного полинома показано на рис.2,
из которого очевидно, что при содержанию экс-
содержание экстракта, %; t суш - температура су-
тракта "меньше нормы", и значениях температуры
шения, °С; К заг - общая кислотность, см3 рас- сушения "норма", общей кислотности "выше нор-
3 мы", содержания редуцирующих сахаров "выше
твора NaOH конц. 1 моль/дм в 100cm3 сусла; нормы", мы сможем получить выход экстракта в
С ред - содержание редуцирующих сахаров, г/100г пределах нормы.
сухого вещества;
Norwegian Journal of development of the International Science No 22/2018 41

Рис. 2. Изменение выхода экстракта при значениях нечетких переменных содержание экстракта
"меньше нормы", температура сушения "норма", общая кислотность "выше нормы", содержание
редуцирующих сахаров "выше нормы"

Модель потерь процесса получения товарно- получена модель, которая соответствует области
го солода. изменения параметра нечеткого понятия "Потери
Для ситуационно-значимой зоны, в которой экстрактивных веществ при производстве солода"
потери редуцирующих сахаров "низкие", темпера- "низкие" вида:
тура сушки "низкая", число Кольбаха "высокое"
Втр е. р. ( Втр ред , t суш , K )  4,304 Втр ред  0,6902 t суш  0,389 K  2,89 Втр ред2 
 0,127096 Втр редt суш  0,0716 Втр ред K  0,0175 t суш  0,0098 K  0,0056 t суш 2  (7)
 0,0128 K 2  21,33
Втр е. р. - потери экстрактивных веществ при рующих сахаров, %; t суш - температура сушения,
производстве солода, %; Втр ред - потери редуци- °С; K- число Кольбаха, %.
Решение данного полинома приведено на
рис.3.
42 Norwegian Journal of development of the International Science No 22/2018

Рис. 3. Изменение потерь экстрактивных веществ при производстве солода при нечетких значениях
входных факторов потери редуцирующих сахаров "низкие", температура сушения "низкая", число
Кольбаха "высокое".

Разработанные математические модели для квадратов оценка коэффициента 1 примет вид

различных ситуационно-значимых зон основных
Cov( X 1 , X 2 )
технологических процессов пивоварения (затира- 1  1   2 , где присутствует
ние, варка сусла с хмелем, брожение и дображива- Var( X 1 )
ние). Cov( X 1 , X 2 )
При построении математической модели, ко- смещение  2
торая адекватно описывает технологический про- Var( X 1 ) .
цесс, очень важен анализ правильности ее специ- При включении несущественной переменной
фикации. в модель смещение оценок коэффициентов не
Негативно на свойствах модели обозначается происходит, но появляется другой недостаток -
как отсутствие значимой переменной, так и чрез- растут стандартные ошибки коэффициентов.
мерное присутствие незначимой. Оценки становятся статистически незначимыми.
В случае, когда в модель не включена суще- Если точная спецификация модели неизвестна (что
ственная переменная, наблюдаются следующие практически всегда и бывает), то используют кри-
последствия: терии, позволяющие выбрать из некоторого мно-
1. Исчезает возможность правильной оценки и жества моделей лучшую.
интерпретации уравнений. Наиболее распространенными критериями яв-
2. Оставшиеся коэффициенты при перемен- ляются критерий Шварца (Schwarz) и критерий
ных становятся смещенными. Акаике (Akaike), которые позволяют выбрать
3. Стандартные ошибки коэффициентов и t- лучшую модель из множества различных специ-
статистики некорректны и не могут быть исполь- фикаций. Критерии численно построены так, что-
зованы для оценки качества созданной модели. бы учесть влияние на качество подгонки модели
Например, если модель имеет вид двух противоположных тенденций. При увеличе-
Yi    1 X1i  2 X 2i   i и из нее изъять пере- нии числа переменных в модели адекватность мо-
дели в общем случае увеличивается. Причем, чис-
менную X 2 , то в новой спецификации фактически ло регрессоров должно быть разумным, чтобы не
будем рассматривать модель Yi    X1i  ui , де вызвать "искусственной подгонки" зависимой пе-
ui   2 X 2i   i . Если переменные X 1 и X 2 кор- ременной поясняющими. С другой стороны, недо-
статочное включение переменных в модель дает
релированы, то нарушается предпосылка теоремы большую стандартную ошибку, и качество под-
Гаусса-Маркова о некоррелированности случайно- гонки снижается.
го члена и регрессоров, так как в этом случае меж- Формулы для расчета критериев Akaike и
ду X 1 и u существует ненулевая корреля- Schwarz записываются:
ция.Оценки, полученные по методу наименьших 2s
квадратов для данной модели, уже не эффективны AIC  ln(ˆ 2 )  (8)
n
среди линейных оценок. Оценки являются даже
смещенными, поскольку для метода наименьших
Norwegian Journal of development of the International Science No 22/2018 43
s ln n 2. Скорректированный коэффициент детер-
SC  ln(ˆ 2 )  (9) минизации:
n
n
где ̂ - выборочная дисперсия; s – число
2 Rk2  1  (1  R 2 ) (11)
ограничений ступеней свободы, n – длина выборки ns
данных. 3. Критерий «предсказания» квадратичной
Значение s в данном случае равно числу неза- погрешности:
висимых переменных в уравнении, включая сво- 1 2
PSE( s)  RSS  ˆ 2 (13)
бодный член. Первое слагаемое в обеих формулах n n
- штраф за большую дисперсию, второй - штраф за 4. Статистика Хокинга:
использование дополнительных переменных. Кри- RSS
SS  (14)
терии рассчитываются для каждой спецификации, (n  s)(n  s  1)
которая рассматривается, и предпочтение отдается 5. AEF критерий средней оцененной диспер-
той модели, где значения критериев меньше. сии:
Кроме названных выше критериев структур- RSS
ной идентификации (структурной адекватности) IS  (15)
регрессионных моделей используют также [8]: n
n  s(ln( )  1) ln 3
1. Коэффициент детерминизации свидетель- s
ствует о плотности взаимодействия между факто- где 3 – «уровень доверия», который показыва-
рами и результирующей переменной и определя- ет с вероятностью 1-3 значения I s дает оценку
ется: среднего риска.
RSS 6. Статистика Фишера указывает на силу
R2  1 (10)
SST взаимосвязи между зависимой переменной и неза-
где висимыми:
n n RSS
RSS   ( y  yˆ i ) 2 F (S ) 
n  s || y  yˆ || 2
(16)
i 1
n Была проведена оценка разработанных моде-
SST   ( yi  yi ) 2 лей-конкурентов определенных структур и пара-
i 1 метров этих моделей на основе критериев: детер-
yi - данные измерений; ŷi - данные, полу- минации, информационного критерия Акаике, ста-
ченные в результате прогнозирования с помощью тистик Хокинга и Фишера (табл.2.6).
модели; yi - усредненное значение данных.
Таблица 1
Критерии адекватности математических моделей.
Созданная модель Критерии
R 2 AIC SS F (s)
Модель (1) 0,31 5,43 0,27 4,48
Модель (2) 0,29 4,36 0,24 4,26
Модель (3) 0,35 5,32 0,31 4,52

Относительная приведенная погрешность нова, В. С. Зубченко, М. В. Карпутіна, М. С. Рома-

прогноза по разработанным моделям не нов // Харчова і переробна промисловість. – 2005.
превышала 5,5%. – № 11. – С. 21-22.
Выводы 2. Юдицкий С.А. Временной анализ систем с
Проведен пассивный эксперимент на всех сетецентрическим управлением на основе непре-
стадиях приготовления пива и осуществлена иден- рывно-дискретной модели / Юдицкий С.А., Точ
тификация сетевых математических моделей в Д.С. // Управление большими системами. – 2010. –
ситуационно-значимых зонах. С. 67-78.
Полученные результаты могут быть 3. Затуливетер Ю. С. Графодинамические си-
использованы в процессе проектирования, стемы с сетецентрическим управлением в матема-
разработки и внедрения новых и модернизации тически однородном поле компьютерной инфор-
существующих систем автоматизации мации / Затуливетер Ю. С. , Фищенко Е. А. //
технологических объектов в пищевой Управление большими системами. – 2010. – 30.1
промышленности, а также при проектировании – С. 567-604.
систем диагностирования. 4. Кульба В. В. Сценарное исследование
сложных систем: анализ методов группового
СПИСОК ЛИТЕРАТУРЫ: управления / Кульба В. В., Кононов Д. А., Чернов
1. Романова З. М. Дослідження впливу И. В., Рощин П. Е., Шулигина О. А // Управление
постійного магнітного поля на інтенсифікацію большими системами. – 2010. – 30.1 – С. 154-186.
процесу приготування пивного сусла / З. М. Рома- 5. Кулинич А. А. Методология когнитивного
44 Norwegian Journal of development of the International Science No 22/2018
моделирования сложных плохо определенных си- 7. Томашевський В. М. Моделювання систем /
туаций / А. А. Кулинич // Избранные труды второй В. М. Томашевський. – К. : Видавнича група BNV,
международной конференции по проблемам 2007. – 352 с.: іл.
управления, 17-19 июня 2003 г. – М. : ИПУ РАН, 8. Дилигенский Н. В. Нечеткое моделирование
2003 г. – С. 219-226. и многокритериальная оптимизация производст-
6. Алтунин А. Е. Модели и алгоритмы приня- венных систем в условиях неопределенности: тех-
тия решений в нечетких условиях: Монография / нология, экономика, экология: Монография / Н. В.
А. Е. Алтунин, М. В. Семухин. – Тюмень : Издате- Дилигенский, Л. Г. Дымова, П. В. Севастьянов. –
льство Тюменского государственного университе- М. : Машиностроение-1, 2004. – 336 с.
та, 2000. – 352 с.

MODELING OF THE PROCESS OF ELECTROMAGNETIC WAVE TRANSFER INTO THE COAXIAL

WAVE CABLE THROUGH THE HIGH-TEMPERATURE SUPERCONDUCTING ELEMENT FROM
THE WAVEGUIDE

Fyk A.
Сandidate of engineering sciences, associate professor, National Academy National Guardia of Ukraine

МОДЕЛИРОВАНИЕ ПРОЦЕССА ПЕРЕДАЧИ ЭЛЕКТРОМАГНИТНЫХ ВОЛН ИЗ

ПРЯМОУГОЛЬНОГО ВОЛНОВОДА В КОАКСИАЛЬНЫЙ КАБЕЛЬ ЧЕРЕЗ
ВЫСОКОТЕМПЕРАТУРНЫЙ СВЕРХПРОВОДЯЩИЙ ЭЛЕМЕНТ

Фык А.И.
кандидат технических наук, доцент, доцент Национальной Академии
Национальной гвардии Украины
г.Харьков, Украина

Abstract
The paper presents the results of the simulation of the penetration process into a coaxial cable of an elec-
tromagnetic wave through a high-temperature superconductor, which is located in the center of a wide wall of a
rectangular waveguide.
Аннотация
В статье представлены результаты моделирования процесса проникновения в коаксиальный кабель
электромагнитной волны через высокотемпературный сверхпроводник, который расположен в центре
широкой стенки прямоугольного волновода.

Keywords: high temperature superconductor (HTS); waveguide.

Ключевые слова: високотемпературний сверхпроводник (ВТСП); волновод.

Постановка проблемы. ника необходимо разработать математический

В современных радиотехнических систем аппарат расчета ее компонент, используя теорию
(РТС) все чаще стали использовать сверхчувстви- электромагнитного поля и теорию сверхпроводи-
тельные полупроводниковые и сверхпроводящие мости[4]
элементы (антенны, элементы памяти), которые Обзор последних исследований и публика-
изготавливаются с использованием методов инте- ций.
гральной технологии, что позволяет существенно В настоящее время разработано множество
уменьшить габариты, вес и стоимость радиоэле- математических моделей (в том числе программ
ментов устройств [1,2]. А использование сверх- (CSC, MW Office и другие), описывающих с раз-
проводящих материалов для СВЧ-систем (волно- ной степенью точности электромагнитные поля в
воды, резонаторы, антенны и т.д.) позволило: различных волноводных системах [5]. Однако, они
уменьшить потери энергии, передаваемой (за счет не позволяют учитывать использование микропо-
уменьшении на три-четыре порядка сопротивле- лосковых элементов (МПЭ) произвольной струк-
ния по сравнению с другими материалам), увели- туры из высокотемпературного сверхпроводящего
чить добротность резонансных элементов, полу- материала, потому что, все электромагнитные тео-
чить возможности управления в скоростных си- рии высокотемпературных сверхпроводников яв-
стемах перенастройки (фазовращатели, ляется эмпирическими. Эти теории высокотемпе-
коммутаторы, электромагнитные защитные ратурной сверхпроводимости используются при
устройства и прочее.). Поэтому значительный ин- известных ограничениях (тип сверхпроводника
терес к разработке математических моделей таких (Монокристалы или пленка разной толщины), па-
сверхпроводящих элементов по настоящее время раметры электромагнитной волны, передаваемой
не утихает [3]. Для изучения распространения (частота, мощность), наличие системы охлаждения
электромагнитной волны в волновом с учетом и прочее.
включения высокотемпературного сверхпровод-
Norwegian Journal of development of the International Science No 22/2018 45
Поэтому для постановки задачи моделирова- Формулировка задачи исследования. Цель
ния процесса передачи электромагнитной волны с статьи - моделирование процесса передачи элек-
прямоугольного волновода к коаксиального кабеля тромагнитной волны из прямоугольного волново-
через высокотемпературную сверхпроводниковых да в коаксиальный кабель через высокотемпера-
пленку нанесенную на цилиндрический стержень турную сверхпроводниковых пленку нанесенную
необходимо разработать математический аппарат на цилиндрический стержень [8].
расчета ее электромагнитных компонент (электри- Изложение основного материала исследо-
ческие и магнитные поля), используя в совокупно- вания
сти теорию электромагнитного поля и теорию 1. Особенности применения уравнений
сверхпроводимости. В статье поясняется переход Максвелла и эквивалентных схем при наличии
от уравнений Максвелла к эквивалентным уравне- в СВЧ-тракте высокотемпературного сверх-
ний поля приведенных в статье [5]. Кроме того, проводящего элемента.
автором ранее [4] было обоснован выбор эмпири- Согласно теории электромагнитного поля со-
ческих моделей сверхпроводимости и это будет отношение между напряжённостью электрическо-
учтено без доказательств при моделировании про- го и магнитного полей в любом СВЧ тракте опре-
цесса передачи электромагнитной волны из волно- деляется уравнениями Максвелла [7]
вода в коаксиальный кабель.
H 
rotE    
t 

E 
rotH  
t 
. (1)
Решение (1) для волновода в прямоугольной системе координат при условии его заполнения диэлек-
триком представляет собой систему из шести скалярных уравнений:
E z E x H y H z H x E y 
   ,   
x z t x z t 
. (2)
E y E z H x H y H z E 
   ,   x
z y t z y t 
E x E y H z H x H y E 
   ,   z
y x t y x t 
Пример структуры волны H Ez  0, Hz  0 приведен рис. 1.
Ey
Z
Hx
Y

Рис.1–Структура волны типа Н в волноводе

Предположим, что электромагнитная гармоническая волна распространяется вдоль оси Z. Тогда ре-
шение (2) можно рассматривать в виде:
E x, y, z, t   E x, y e jt  z 


H x, y, z, t   H x, y e jt  z 


. (3)
46 Norwegian Journal of development of the International Science No 22/2018
jt
При этом в выражении (3) множитель e описывает гармоническое изменение поля во времени с
частотой  , а  – коэффициент распространения ЭМВ. Тогда уравнение (2), учитывая, что поле одно-
родно вдоль оси Y, можно записать в виде:
E z 
 E x   jH y 
x 
 E y   jH x 

E y 
   jH z 
x

H z
 H x  jE y  . (4)
x 

 H y  jE x 

H y

  jE z
x 
Подставив выражения для Ey и H y в уравнения (4), получим следующую систему уравнений:
H z 
  2 H x   
x 
 j H z 
Ey   2 
 x 

E

  2 E x   z
dx , (5)

j  j E z 

H y   Ex   2
  x 

где  2  2   2 - волновое число.

Учитывая, что в прямоугольном волноводе основным типом волны является H волна
E z  0, H z  0 , перепишем выражение (5) в виде [7,9]:
 H z 
H x   2 
 x 

 j H z 
Ey   2
 x 
. (6)
Решение системы уравнений (6) можно записать следующим образом [5,7]:

E y   e z F x, y e jt 



H x   h z G x, y e jt 

 
H z   q z g x, y e jt  ,

(7)
 
где e  z, h  z, q  z - функции распределения по координате t, которые именуют скаляр-
ными комплексными амплитудами, F (x,y); G  (x,y); g  (x,y) - собственные функции, характеризую-
щие распределение составляющих поля по координатам поперечного сечения для волн Н-типа,  - по-

рядковый номер волны, распространяющейся в волноводе   1 H10 [5,7]. При включении в

волноводный тракт ВТСП-элемента, состоящего из диэлектрической подложки цилиндрической формы
и нанесенной на её поверхность сверхпроводящей плёнки (ВТСП-стержень) (рис 2) [5].
Norwegian Journal of development of the International Science No 22/2018 47

I II Z

z0

n
z z+∆z
X

Рис. 2.–Схематическое отображение области с высокотемпературным-стерженем

Предполагая, что E,H – это поля, наводящиеся на сверхпроводящем элементе, а E1 , H1 – век-

торы независимого электромагнитного поля при отсутствии ВТСП-стержня, перепишем выражение (4)
для обеих волн в виде [7]:
rotE   jH
 

rotH  jE  J 

rotE 1   jH
1 
  jE  J 
rotH 1 1 1 . (8)
Умножая векторно слева третье уравнение (8) на H  , четвёртое на E
 , второе на E , первое на H
1 1
соответственно и просуммировав почленно полученные равенства, получим:
H rotE  E rotH  H
 rotE  E rotH
  J E  J E .(9)
1 1 1 1 1 1
Проинтегрировав выражение (9) по объёму и используя теорему Гаусса-Остроградского, с учётом,
что divab  brota  arotb , получим [5]:
  E 1  H  E  H 1dS   J1E  JE 1 dV
S V . (10)
Анализируя распределения электромагнитного поля в волноводе в окрестности вокруг сверхпрово-
 
дящего стержня z, z  z , поверхностный интеграл в (10) представим как сумму трёх поверхност-
ных интегралов: интегралов по поперечным сечениям I и II и интеграла по боковой поверхности слоя
вокруг стержня:

 A d S   A d S   A d S   A d S
S I II Sбок

  
, (11)
где   E  H
A   E  H
 .
1 1 
Поскольку в рассматриваемых интегралах векторы dS ориентированы по внешним нормалям, то,
o
вводя единичные векторы – продольный z и орт внутренней нормали n , выражение (11) запишем в
виде[5]:

 A dS   A z dS   A z dS   A ndS
o o

S I II Sбок . (12)

Интеграл от произведения ( An ) по боковой поверхности можно преобразовать к виду[5]:
48 Norwegian Journal of development of the International Science No 22/2018

  ndS 
A  E 1  H  E  H 1ndS   n  E 1H  n  E H 1dS . (13)

Sбок Sбок

Для E 1 , представляющего собой поле в волноводе без источника (стержня), справедливо граничное
условие  
n  E  0 . Тогда выражение (13) можно записать следующим образом:
1
 A ndS    n  E H 1dS
Sбок Sбок . (14)
Обозначая интеграл от  n по поперечному сечению z как:
A
 z odS  E  H
Qz    A  1   E  H 1z dS
o ,(15)

Sп Sп
где Sп -площадь поперечного сечения волновода,запишем соотношение (10) в виде[5]:
Qz  z   Qz     J E 1dV   n  E H 1dS
.(16)

V Sбок
Соотношение (16) показывает, что различия между интегралами по поперечным сечениям
Qz  z  и Qz  определяются эквивалентным источником (ВТСП-стержнем). В дальнейшем для
расчётов предполагаем, что в прямоугольном волноводе распространяется волна H10 . Поэтому напря-
жённости в области стержня и вне его согласно [9] можно записать соответственно:
E  E y , H
 H
 H
x
 

z
E  E , H
 H
 H 
1 1y 1 1x 1z 
. (17)
Учитывая выражение (7), можно записать:
E  ez F x, y 
 

H  hz G x, y   q z g x, y  

E1  e 1 z F x, y  



H 1  h 1 z G x, y   e 1 z g x, y 
 1
(18)
 
Подставляя выражение (18) в уравнение (15), получим, что слагаемые с продольными векторами да-
o
дут ноль (т.к в смешанное произведение входит продольный орт z ). Меняя порядок интегрирования и
суммирования, запишем (15) в виде:
Qz    
E 1  H 
  E  H 
 z odS  e h F Gz odS e h
 
o . (19)
1 1  1 FG z dS
Sп Sп Sп
Учитывая, что F  G z o  F G  z o   F F , и F  G   FF , выражения (19), следует за-
писать в виде:


Qz   e 1h  e h 1 F FdS
. Sп (20)
Из условий нормировки, рассмотренных в работе [5,7,8], получим   1, т.е. возможно только рас-
пространение Н-волны через рассматриваемую область z, z  z  [5]. Тогда выражение(20) можно
записать в виде:

Qz   e1z h z   e1z h 1z  , (21)


1 z  e 1 z .
где e 
Разделив обе части равенства (16) на z и переходя к пределу при z  0 , записываем:
Norwegian Journal of development of the International Science No 22/2018 49

 JE 1dS   n  E H 1dl

dQ

dz
Sп C . (22)
Учитывая выражение (18) и (21), перепишем выражение (22) в виде:
 
 
 
d
    
e1z h z   e z h1z   e1z   J FdS 
  1


n  E gdl  
dz  
 C 
 
 
 
 S

 h1z   J fdS  n  E Gdl ,

 п C  (23)
где C – рассматриваемый замкнутый контур  z, z  z
Для упрощения (23) введём вспомогательные функции:

 nE gdl

z    
1
J FdS 

Sп  C

 
 z    JfdS  n  E Gdl 


Sп C . (24)
Тогда выражение (23) можно переписать в виде:

e1z h z   e z h 1z   e1z  z   h 1z  z 

d
dz . (25)
Предположим, что e1 z  e , h 1 z  e 
z  z 1
соответствует прямой волне, а

e1 z  e z и h 1 z  e z1 соответствует обратной волне. Тогда, подставив в выражение (25)

и взяв полусумму, получим[5]:
dh 
 Y0e   z 
dz 
, (26)
de
  Z0 h   z 
dz 

где Y0  , Z 0   .

Волновод c источниками тока и напряжения можно представить в виде эквивалентной длинной ли-
нии. Тогда, учитывая (26), уравнение волны перепишем в виде:
dI    z 
 Y0 U 
dz
 . (27)

  Z0I   z  
dU
dz 
Для волновода со сверхпроводящим стержнем систему уравнений (27) можно записать в виде:
dI    z  .
 Y0 U (28)
dz
Тогда анализ волноводного тракта со сверхпроводником можно проводить по схеме, представлен-
ной рис .3.
50 Norwegian Journal of development of the International Science No 22/2018

I I  dI

  dU

Y0dz 
 (z)dz U

dz
Рис. 3. Эквивалентная схема волновода с ВТСП-элементом
Из уравнения (26) можно определить h  z , а из (18) определить поперечные составляю-
e  z и
щие поля. Из выражения (4) можно определить выражение для продольной составляющей поля H z . Но
выражение (26) фактически есть описание волноводного тракта как длинной линии, и поэтому необхо-
димо в приведенных терминах определить продольную составляющую поля  H z  . Дифференцируя
второе уравнение в системе (26) и подставляя в него первое уравнение с учётом (24), получим:
 
d 2e d   d  

 Z0  Y0e   z     JfdS    n  E Gdl  ,

  (29)
 dz  C 
2 dz 
dz
 Sп  
  e z  .
где e
В данном уравнении дифференцирование интегралов даёт ноль, что позволяет переписать (29) в ви-
де:

 n  E gdl
d 2e
  2e   . (30)
2
dz
Sп
В работе [10] доказано, что:
 
 
1  
q z    e z   n  E gdl  , (31)
 2 
 Sп 
тогда

e z   q 
 2  n  Egdl . (32)

Sп
Подставив e  z в уравнение (30) и учитывая (32), получим:

 
     
d 2q d 
   2  2
  n  E gdl    q  n  E gdl    n  E gdl .
2
(33)
2 2 2
dz  dz   
 Sп  Sп Sп
или:

d 2q   2 2 d 
2
 q  2
   n  Egdl .  (34)
dz 2  2  dz 2 
Sп
Norwegian Journal of development of the International Science No 22/2018 51

Обозначая  02   2   2 , перепишем (34) в виде:


d2H   2 d 2  
z 2
  Hz   0  gUz , (35)
dz 2  э
 2
dz 

  1 d x 0 , y0  
2
 z gx, y  ,

1
где n  E gdl  g nEdl  gUz  , э  
 
2
, Hz  q

Sп
x0 , y0 - координаты проводящего стержня,  02   2 ,U(z)(7,8)
Таким образом, систему для определения поля из (18), (26) и (35) преобразуем к виду[5]:

dH 
x  G x, y  Y0 e   z  
dz 

dE y
dz

 Fx, y   Z 0 h   z  
. (36)

2
d Hz   2 d  2 
2
  Hz     gUz 
dz 2  э  0 dz 2  

Приведенное распределение верно, если ВТСП-стержень находится в сверхпроводящем состоянии.
Определив приведенным методом эквивалентных длинных линий токи и напряжения (27), из (36) можно
получить электромагнитное поле в исследуемой области.
Очевидно, что передача энергии поля из волновода в коаксиальный кабель производится по сверх-
проводяшему стержню. Тогда магнитное поле, которое наводится в нем, можно определить из уравнени-
ий Лондона[4]:
J  J n  Js 
J n  n E 


B ,
1
rotJ   (37)
 0 2L 
1 
Js  E 
 0 2L 
где J n , Js -плотность тока нормальных и сверхпроводящих носителей заряда, 0 - магнитная
проницаемость вакуума, L - лондоновская глубина проникновения магнитного поля в сверхпровод-
ник.
2. Результаты моделирования процесса распространения электромагнитной волны из волно-
вода в коаксиальный кабель.
Выражение для векторов поля бегущей волны типа Н10 в прямоугольном волноводе без потерь со-
гласно (36)имеет вид:

Es  A1F( x, y)e j(t  z) 


j(t  z ) 
Hs  A1G ( x, y)e 

H z  A1g( x, y)e j( wt  z) 
 , (38)
52 Norwegian Journal of development of the International Science No 22/2018
j(t z)
где e - волновой множитель, который описывает закон распределения поля вдоль коор-
динаты z для различных моментов времени; A1 – постоянный множитель, определяющий интенсивность
поля,C-const[7-9].
Известно, что перенос энергии через поперечное сечение волновода производится каждым типом
волны независимо от другого типа, а этот факт позволяет применить здесь понятие ортогональности[7].

a b

  
2
F ds  C2 sin 2 xdxdy  1
a
x 0 y 0 . (39)
Из выражения (39) определяем константу С:

2
С
ab . (40)
Теперь учтем наличие ВТСП-стержня и запишем выражение для напряженности электрического по-
ля (36):

E s  e (z)F( x, y)e jt , (41)

.
где e( z ) - функция, характеризующая распределение составляющей поля Es по координате z.
Из формул (38 - 41) можно получить выражение для напряженности электрического поля [5]:
 
h д0 sin x1 sin x (sin 2 l1  j sin 2l1)e j(t z) , (42)
2 1
E y  I 0
ab a a 2
a
где I0 – ток в стержне, x1  - расстояние между узкой стенкой и стержнем,l1 – расстояние от
2
2
стержня до закороченной стенки волновода,   - коэффициент фазы. Схематично на рис 4а,б по-
в
казан процесс переноса энергии (изменение направления вектора Умова-Пойтинга) из волновода в коак-
сиальный кабель. На данном рисунке также отображена динамика наведения тока в стержне. На рис. 4
в,г изображено пространственное распределение поля для стандартного прямоугольного волновода

i
. Ey

а в
Norwegian Journal of development of the International Science No 22/2018 53
Ey t
2

t2>
t t1
1

0 z

б г
Рис.4-Результаты моделирования процесса наведения токов напряжённостью поля електромагнитной
волной в волноводе:
а,б- Процесса передачи энергии из волновода в коаксиальный кабель
в,г-Напряжённость поля (Ex) на сверхпроводящем елементе на моменты времени t 1,t2;
Определение напряжённости магнитного поля производится согласно выражению :

1
Hx  Ey
0 (43)
При выводе формул (3.16 – 3.21) предполагалось, что стенки волновода являются идеальными про-
водниками. Однако реально они обладают не бесконечно большой, а конечной проводимостью (Т=77К).
Поэтому некоторая доля энергии, передаваемой по волноводу, затрачивается на нагревание стенок, по
которым протекают высокочастотные токи. Поскольку для различных типов волн распределение и зна-
чение токов разные, то и потери в стенках волновода будут зависеть от типа волны, которая распростра-
няется в волноводе. Величина потерь для волны Hmn характеризуется коэффициентом затухания [50]:

1  2
2 a m n  b n m
2
2 a n  b m 
m, n  (1  1 )  1 , (44)
2 2 2 2
2 
120(1  1 )  a n b m ab 

 m2 n 2 1
где 1   ;  - удельная проводимость стенок (медь =5,8107 (Ом  м) ,
2 a 2 b2

1  т, n  0
 n , m  
2  n, m  0 . (45)
Из формул (44) и (45) можем получить величину затухания для волны Н10:

1  2a 2 
10    
 2 
120(4a 2  2 )  b a b  (46)
В выражение (42) и (43) добавляется множитель e
 z и учитываются таким образом потери в
стенках волновода..
Передаваемая мощность из волновода в коаксиальный кабель рассчитывается следующим образом
[7]:
54 Norwegian Journal of development of the International Science No 22/2018

1
P  Re(I 2э Z
 э)
2 . (47)
Подставляя выражения для Iэ и Z э [5], получим выражение для расчета мощности, передаваемой
на волне типа Н10.

(I 0m h д ) 2 2 
P  0 sin ( x1) sin 2 (l1)
ab a . (48)
Максимальная длина волны макс, при которой может распространяться волна Н10 в прямоугольном
волноводе, равна: макспр. Н10 =2а. Минимальная длина волны min определяется из условия min
  H 20 =а. Опеределим частотный диапазон КВП:
cc
f max  f крН20  
крН20 a 

c c
f min  f крН10  
крН10 2a 
. (49)
Вывод 5. Гупта К., Гардж Р., Чадха Р. Машинное
В результате электродинамического модели- проектирование СВЧ устройств. /Пер. с англ. -М.:
рования процесса передачи электромагнитной Радио и связь. 1987. - 432с.
волны из пространства волновода в коаксиальный 6. Фык, А.И. Моделирование электромагнит-
кабель через высокотемпературный сверхпрово- ного поля в прямоугольном волноводе с сверхпро-
дящий стержень получении аналитические выра- водящим стержнем // Международный Научный
жения для определения магнитной, электрической Институт "Educatio" Ежемесячный научный жур-
составляющей поля в переходе волновод- нал, / –Новосибирск: 2015. – Вып № 6 (13). Ч..2–
сверхпроводник-коаксиал с учетом использования С. 17–21
эмпирических моделей сверхпроводимости. Пред- 7. Ширман Я.Д. Радиоволноводы и объёмные
ложенный математический аппарат расчета с ис- резонаторы М.: Сов. радио, 1959.-370с.
пользованием теории электромагнитного поля и 8. Фык А.И Experimental research of microstrip
эмпирических моделей сверхпроводимости позво- aerial as protective device of receiver from electromag-
ляет описать электродинамические процессы в netic defeat / D.B. Кucher, R.О. Gonchar // EUREKA:
СВЧ-системах где стандартные программы элек- Physics and Engineering" (Tallin, Estonia),No 5 (2017)
тродинамического расчета (AWR, CST STUDIO, – P.54-53
Proteus и т.д) использовать затруднительно. 9. Гольдештей Л.Д. Электромагнитные поля и
волны, изд.2.,-М.: Сов. радио. 1971.201с.
СПИСОК ЛИТЕРАТУРЫ: 10. Ковалёв И.С. Основы теории и расчёта
устройств СВЧ.-Минск:Наука и техника,-1972.-
1. Щербаков В.И. Эффективность плоских ан- 200с
тенных решеток для систем непосредственного 11. Кучер Д.Б, Використання тонких плівок
спутникового телевизионного вещания. XXVI високотемпературних надпровідників для захисту
НТК «Теория и техника антенн». АО «Радиофизи- елементів радіоелектронної апаратури від
ка».-М., 1994,С.158-161. електроерозійного руйнування при впливі потуж-
2. Huang, “Microstrip reflectarray antenna for the них електромагнітних випромінювань / Фик О.І.
SCANSCAT radar application”, JPL Publication No. В.О.Колпаков В.Д.Калугін //Системи обробки
90 - 45, November 15, 1990. інформації Збірник наукових праць -2002-
3. Кучер Д. Б. Мощные электромагнитные из- Вип.5(21)-Харків: ХВУ, – С. 36-41
лучения и сверхпроводящие защитные устройства. 12. London,H, Proc.Roy. Soc.A,1940,176, no.
- Монографія / Севастополь: Ахтиар, 1997. - 188 с. 964, pp.522-533
4. Fyk А.І . Model of process of influence of the 13. Гинзбург В.Л. Замечания о высокотем-
electromagnetic field on high temperature supercon- пературной сверхпроводимости // ВТСП, 1990. NI.
ducting microstrip antenna Magyar Tudományos 0.7-26.
Journal (Budapest, Hungary) VOL 1, No 9 (1) (2017)
– P.34-44
Norwegian Journal of development of the International Science No 22/2018 55

VETERINARY SCIENCES
SABILITY OF BIOCHEMICAL PARAMETERS OF FECES DURING STORAGE

Zhukov M.
Candidate of Veterinary Sciences, Senior Researcher, Department of Experimental Therapy, State Scientific
Institution All-Russian Veterinary Research Institute of Pathology, Pharmacology and Therapy of the Russian
Academy of Agricultural Sciences, Russian Federation, Voronezh
Bosenko T.
Graduate student, State Scientific Institution All-Russian Veterinary Research Institute of Pathology, Phar-
macology and Therapy of the Russian Academy of Agricultural Sciences, Russian Federation, Voronezh
Morgunova V.
Candidate of veterinary sciences, leading researcher of the laboratory of biochemistry of blood, State Sci-
entific Institution All-Russian Veterinary Research Institute of Pathology, Pharmacology and Therapy of the
Russian Academy of Agricultural Sciences, Russian Federation, Voronezh
Kashirina L.
Junior Researcher of the Laboratory of Biochemistry of Blood, State Scientific Institution All-Russian Vet-
erinary Research Institute of Pathology, Pharmacology and Therapy of the Russian Academy of Agricultural
Sciences, Russian Federation, Voronezh
Tyurina E.
Junior Researcher of the Laboratory of Biochemistry of Blood, State Scientific Institution All-Russian
Veterinary Research Institute of Pathology, Pharmacology and Therapy of the Russian Academy of Agricul-
tural Sciences, Russian Federation, Voronezh

СТАБИЛЬНОСТЬ БИОХИМИЧЕСКИХ ПОКАЗАТЕЛЕЙ КАЛА ВО ВРЕМЯ ЕГО

ХРАНЕНИЯ

Жуков М.С.
Кандидат ветеринарных наук, старший научный сотрудник отдела экспериментальной терапии,
ГНУ Всероссийский научно-исследовательский ветеринарный институт патологии, фармакологии и
терапии Россельхозакадемии, Российская Федерация, г. Воронеж
Босенко Т.Д.
Аспирант, ГНУ Всероссийский научно-исследовательский ветеринарный институт патологии,
фармакологии и терапии Россельхозакадемии, Российская Федерация, г. Воронеж
Моргунова В.И.
Кандидат ветеринарных наук, ведущий научный сотрудник лаборатории биохимии крови, ГНУ Все-
российский научно-исследовательский ветеринарный институт патологии, фармакологии и терапии
Россельхозакадемии, Российская Федерация, г. Воронеж
Каширина Л.Н.
Младший научный сотрудник лаборатории биохимии крови, ГНУ Всероссийский научно-
исследовательский ветеринарный институт патологии, фармакологии и терапии Россельхозакадемии,
Российская Федерация, г. Воронеж
Тюрина Е.В.
Кандидат ветеринарных наук, ведущий научный сотрудник лаборатории биохимии крови, ГНУ Все-
российский научно-исследовательский ветеринарный институт патологии, фармакологии и терапии
Россельхозакадемии, Российская Федерация, г. Воронеж

Abstract
The comparatively high risk of an analytic error associated with changes in the composition of stool during
storage limits the ability to perform its analysis. Studies have been carried out to study the variability of the feces
of cattle and pigs in the different term of storage. The absence of veracious interspecific differences in the stabil-
ity of feces during storage is shown. It was revealed that during the storage of stool specimens at a temperature
of -18 ℃, changes in its composition occur, which must be taken into account when carrying out coprological
studies. The features of changes in feces of content glucose, albumins, total lipids and triglycerides were studied.
Maximum storage stability was noted in albumins, but minimal in triglycerides. In order to reduce the risk of an
analytical error, it is recommended to conduct feces analysis in the first hours after the selection, but if neces-
sary, it is allowed to freeze for a term to 24 hours.
Аннотация
Сравнительно высокий риск аналитической ошибки, связанный с изменением состава кала во время
хранения, ограничивает возможности выполнение его анализа. Проведены исследования по изучению
вариабельности показателей кала крупного рогатого скота и свиней при разном сроке его хранения. По-
казано отсутствие достоверных межвидовых различий в стабильности показателей кала при хранении.
56 Norwegian Journal of development of the International Science No 22/2018
Выявлено, что в процессе хранения проб кала при температуре -18 ℃, происходят изменения его состава,
что необходимо учитывать при проведении копрологических исследований. Изучены особенности изме-
нения в кале содержания глюкозы, альбуминов, общих липидов и триглицеридов. Максимальная устой-
чивость при хранении отмечена у альбуминов, но минимальная у триглицеридов. С целью снижения
риска аналитической ошибки рекомендуется проводить исследования кала в первые часы после отбора,
но при необходимости допускается его заморозка на срок до 24 часов.

Keywords: cattle, pigs, feces, stability of indices during storage, pre-analysis stage of study, permissible er-
ror limits
Ключевые слова: крупный рогатый скот, свиньи, кал, стабильность показателей при хранении,
преданалитический этап исследования, допустимые пределы погрешности.

Введение. Заболевания желудочно- Материалы и методы исследования. В

кишечного тракта у животных входят в число условиях предприятия специализирующегося на
наиболее актуальных в плане распространённости откорме молодняка крупного рогатого скота на
и экономической значимости. Они могут быть как основании результатов клинического обследова-
самостоятельной патологией незаразной природы, ния были выбраны 7 телят в возрасте 6 месяцев
так и компонентом симптомокомплекса инфекций больные энтероколитом. На промышленном сви-
[4]. Постоянный контакт пищеварительного тракта новодческом комплексе с аналогичным диагнозом
с внешней средой, тесная физиологическая связь с выделили 7 поросят в возрасте 45 суток. У всех,
другими органами и системами организма, обу- задействованных в опыте животных утром за 1,5-2
славливают сложную этиологическую структуру часа до кормления из прямой кишки отобрали
желудочно-кишечных болезней со сравнительно пробы кала. Материал в количестве 50 г/гол в спе-
высоким риском их развития, а многогранный па- циализированном термостабильном контейнере,
тогенез затрудняет их диагностику и терапию. По- оснащённом хладореагентами, доставили в лабо-
этому, несмотря на достижения науки и практики, раторию в течение 1,5 часа. В условиях лаборато-
количество не решённых и спорных проблем в рии кал от каждого животного гомогенизировали,
современной гастроэнтерологии постоянно увели- а затем делили на 4 образца по 10 г: №1, №2, №3 и
чивается [10]. На большинстве животноводческих №4. К каждому из которых добавили 20 г дистил-
предприятиях регистрируются массовые желудоч- лированной воды. Полученную смесь тщательно
но-кишечные болезни телят [7] и поросят [1], диа- перемешивали стеклянной палочкой в течении 1
гностика которых проводится комплексно с ис- минуты после чего она была оставлена в тёмном
пользованием клинических и лабораторных мето- месте на 30 минут при комнатной температуре.
дов. Из числа лабораторных тестов преимущество Далее полученную взвесь центрифугировали 15
используются микробиологические исследования мин при 3000 об/мин и отбирали супернатант.
тканей органов и крови, так же часто проводятся Исследование супернатанта образца №1 было
биохимический анализ крови, но относительно проведено через 4-6 часов после отбора, а осталь-
реже анализируется кал, что обусловлено особен- ные образцы были заморожены в холодильной
ностями данного материала. Так, сложный состав камере при температуре «-18,0 °С». Образец №2
кала, включающий большое количество микроор- анализировали через 24 часа после отбора, №3 и
ганизмов, свободных метаболитов и биологически №4 соотвественно через 48 и 144 часов. При этом
активных соединений (ферменты, гормоны и др.) определяли уровень глюкозы, альбуминов, общих
повышают риск вариабельности показателей и в липидов и триглицеридов по общепринятым мето-
первую очередь аналитических погрешностей, дикам с помощью коммерческих наборов реакти-
связанных с процессами отбора, транспортировки вов.
и преданалитической подготовки. В одних норма- Математико-статистическую обработку полу-
тивных документах допускается хранение кала ченных данных проводили с помощью приклад-
при температуре от 0 до 4 °С не более суток [6], в ных программ Statistica v6.1 и Microsoft Excel. Рас-
других до 12 часов при температуре 3-5 ℃ [5] при считывали среднюю арифметическую и её ошиб-
этом все авторы не рекомендуют замораживать ку, а так же достоверность разницы (р) по
пробы. Однако в реальных условиях животновод- критерию Стьюдента.
ства и организации лабораторных исследованиях Результаты исследований и обсуждения.
выполнить эти требования невозможно, необхо- Данные анализа кала представлены в таблице 1, из
димы более приемлемые требования. Поэтому бы- которой видно, что у телят и поросят количество
ли проведены исследования по изучению вариа- глюкозы в сравнении с исходным уровнем умень-
бельности показателей кала при разном сроке его шилось в течение 24, 48 и 144 часов соотвественно
хранения. на 1,5 и 2,0 %; 17,0 и 16,5 %, 34,5 и 33,8 %.
Norwegian Journal of development of the International Science No 22/2018 57
Таблица – 1
Изменение биохимических показателей проб кала телят и поросят во время их хранения
Продолжительность хранения проб (час)
Показатели
4-6 24 48 144
2,32±0,130 2,28±0,132 1,93±0,110* 1,52±0,101***
Глюкоза, ммоль/л
4,50±0,370 4,41±0,119 3,76±0,204* 2,98±0,133***
11,5±0,58 11,3±0,71 10,8±0,49 9,5±0,60*
Альбумин, г/л
7,90±0,301 7,79±0,251 7,44±0,270 6,58±0,301**
3,20±0,091 3,12±0,100 2,72±0,091*** 1,98±0,10***
Общие липиды, г/л
4,60±0,240 4,51±0,182 3,89±0,220* 2,72±0,089***
1,10±0,057 1,03±0,050 0,85±0,032*** 0,62±0,027***
Триглицериды, ммоль/л
2,80±0,095 2,64±0,088 2,16±0,042*** 1,60±0,073***
Примечание: в верхнем сроке показатели телят, в нижней – поросят; * - р≤0,05; ** - р≤0,01 и *** -
р≤0,001 в сравнении с данными предыдущего периода хранения.

Выбор глюкозы для оценки вариабельности стическое значение [3, 9]. Однако кал обладает
состава кала был обусловлен тем, что это наиболее сравнительно высокой протеолитической активно-
востребованный метаболит в процессах полостно- стью, и в случае сохранения её значимого уровня
го и симбионтного пищеварения в кишечнике [2], белки будут разрушаться, поэтому динамика аль-
поэтому, уменьшение её количества отражает сте- бумина отражает степень торможения функций
пень активности свободных (полостных) фермен- ферментов во время хранения. Результаты наших
тов и жизнедеятельности микрофлоры. Получен- исследований показали, что замораживание оказы-
ные нами результаты показали, что замораживание вает на белки более выраженный стабилизирую-
снижает, но не останавливает потребление глюко- щий эффект чем на углеводы и жиры.
зы в исследуемом материале. Таким образом, в процессе хранения проис-
Липидограмма кала имеет высокую диагно- ходят изменения состава кала, выраженность ко-
стическую ценность, позволяя оценивать функции торых не имеет существенного межвидового раз-
поджелудочной железы и кишечника. У здоровых личия. Для оценки аналитической значимости вы-
животных в кишечнике всасывается до 95-98 % явленных изменений мы использовали два
жира, поступающего с кормом, поэтому его со- подхода. Во-первых, была рассчитана достовер-
держание в кале минимально, но при патологии ность различия по критерию t Стьюдента. В меди-
оно увеличивается [3, 9]. При этом потери липидов ко-биологических исследованиях достоверными
во время хранения могут нивелировать диагности- являются различия при t=2, что соответствует Р =
ческое значение липидограммы. Так же следует 95% [8]. Из данных таблицы 1 видно, что досто-
отметить, что нейтральные жиры (триглицериды) верные изменения произошли в содержании глю-
легко плавятся и они легче воды. При этом, высок козы через 48 часов, альбуминов – 144 ч, общих
риск их потерь на фоне происходящих структур- липидов – 48 ч и триглицеридов – на через 48 ч.
ных метаморфоз во время замораживания и размо- Однако, помимо статистической достоверно-
раживания, что так же увеличивает риск аналити- сти различий, имеется допустимый предел по-
ческой ошибки. Во время хранения образцов от грешности значений, определяемый вариабельно-
телят содержание общих липидов и триглицеридов стью метода исследования и метаболического
в сравнении с исходным уровнем сократилось в профиля организма. Поэтому вторым вариантом
течение 24 часов на 2,5 и 6,0 %, 48 часов – на 15,0 оценки значимости выявленного нами изменения
и 23,0 % и 144 часов – на 41,0 и 43,5 %. У поросят показателей кала во время хранения, была оценка
так же отмечено уменьшение указанных показате- соответствия наших результатов, допустимым
лей соотвественно на 2,0 и 5,8%, 15,4 и 22,7%, 40,8 пределам аналитической вариации, принятым ре-
и 43,0%. Представленные данные подтвердили ферентной группой экспертов по научной разра-
прогнозируемый риск потерь липидов в кале во ботке и унификации методов и средств клиниче-
время его хранения. Так, в сравнении с первыми ской лабораторной диагностики [5]. В соответ-
24 часами хранения во время вторых суток интен- ствии с решением этой группы допускаются
сивность потерь общих липидов возраста у телят и пределы вариации показателей альбуминов – 3%,
поросят соотвественно в 6,2 и 7,7 раза, а триглице- триглицеридов – 7%, глюкозы и общих липидов -
ридов в 3,8 и 3,9 раза. 5%. Полученные нами данные показали, что в те-
Содержание альбумина в течение 24, 48 и 144 чение 24 часов хранении изменения всех изучае-
часов хранения кала телят и поросят снизилось мых показателей не превышает допустимые пре-
соотвественно на 1,6 и 1,4 %, 6,0 и 5,8 %, 17,0 и делы вариации, но при увеличении этого срока
16,7 %. Известно, что увеличение содержания аль- уровень аналитической вариабельности становится
бумина в кале указывает на наличие воспалитель- не приемлемым (диаграмма 1).
ного процесса, что определяет важное его диагно-
58 Norwegian Journal of development of the International Science No 22/2018

Диаграмма – 1. Динамика биохимических показателей проб кала крупного рогатого скота и свиней в
процессе их хранения в замороженном состоянии (среднее процентное отклонение от исходного уровня)

Заключение. Копрологические исследования 2. Бергнер Х. Научные основы питания сель-

существенно расширяют методологические воз- скохозяйственных животных / Х. Бергнер, Х-А.
можности ветеринарной гастроэнтерологии, по- Кетц. – М.: Колос, 1973. – 597 с.
вышая объективность диагноза и обоснованность 3. Битти А.Д. Диагностические тесты в гастро-
назначаемой терапии. Однако высокий риск ана- энтерологии // А.Д. Бити. Пер. с англ. – М.: Меди-
литической ошибки связанный с изменением со- цина, 1995. – 224 с.
става кала во время хранения ограничивают воз- 4. Калюжный И.И. Клиническая гастроэнтеро-
можности проведения его анализа. Полученные логия животных / под ред. И.И. Калюжного. –
нами данные показали отсутствие достоверных СПб.: Изд-во Лань, 2015. – 448 с.
межвидовых различий в стабильности показателей 5. Лабораторные методы исследования в кли-
кала при хранении. Выявлено, что в процессе хра- нике: Справочник // В.В. Меньшиков, Л.Н. Делек-
нения проб кала при температуре -18°С, происхо- торская, Р.П. Золотницкая и др. – М.: Медицина,
дят изменения его состава, что необходимо учиты- 1987. – 368 с.
вать при проведении копрологических исследова- 6. Лабораторная диагностика гельминтозов и
ний. Изучены особенности изменений в кале протозоозов: Методические указания. 2-е изд. –
содержания глюкозы, альбуминов, общих липидов М.: ФБУЗ «Федеральный центр гигиены и эпиде-
и триглицеридов. Максимальная устойчивость при миологии» Роспотребнадзора, 2014. – 118 с.
хранении отмечена у альбуминов, но минимальная 7. Николаева О.Н. Этиология и профилактика
у триглицеридов. С целью снижения риска анали- желудочно-кишечных болезней телят / О.Н. Нико-
тической ошибки рекомендуется проводить иссле- лаева // Практик. – 2010. – №1. – С. 26 – 31.
дования кала в первые часы после отбора, но при 8. Реброва О.Ю. Статистический анализ меди-
необходимости допускается его заморозка на срок цинских данных // О.Ю. Реброва. – МедиаСфера,
до 24 часов. 2002. – 312 с.
9. Шевченко И.А. Лабораторные методы ис-
СПИСОК ЛИТЕРАТУРЫ: следования при заболеваниях органов пищеваре-
ния // И.А. Шевченко. – Л.: Медицина, 1982. – 200
1. Алтухов Н. Пути профилактики желудочно- с.
кишечных болезней поросят в период отъёма / Н. 10. Циммерман Я.С. Нерешенные и
Алтухов, Ю. Бригадиров, А. Шамардина // Свино- спорные проблемы современной гастроэнтероло-
водство. – 2005. – №6. – С. 21-22. гии / Я.С. Циммерман. – М.: МЕДпрессинформ,
2013. – 224 с.
 

№22/2018

Norwegian Journal of development of the International Science

ISSN 3453-9875

VOL.1

It was established in November 2016 with support from the Norwegian Academy of Science.

DESCRIPTION
The Scientific journal “Norwegian Journal of development of the International Science” is issued 12 times a year
and is a scientific publication on topical problems of science.

Editor in chief – Karin Kristiansen (University of Oslo, Norway)

The assistant of theeditor in chief – Olof Hansen

 James Smith (University of Birmingham, UK)

 Kristian Nilsen (University Centre in Svalbard, Norway)
 Arne Jensen (Norwegian University of Science and Technology, Norway)
 Sander Svein (University of Tromsø, Norway)
 Lena Meyer (University of Gothenburg, Sweden)
 Hans Rasmussen (University of Southern Denmark, Denmark)
 Chantal Girard (ESC Rennes School of Business, France)
 Ann Claes (University of Groningen, Netherlands)
 Ingrid Karlsen (University of Oslo, Norway)
 Terje Gruterson (Norwegian Institute of Public Health, Norway)
 Sander Langfjord (University Hospital, Norway)
 Fredrik Mardosas (Oslo and Akershus University College, Norway)
 Emil Berger (Ministry of Agriculture and Food, Norway)
 Sofie Olsen (BioFokus, Norway)
 Rolf Ulrich Becker (University of Duisburg-Essen, Germany)
 Lutz Jäncke (University of Zürich, Switzerland)
 Elizabeth Davies (University of Glasgow, UK)
 Chan Jiang(Peking University, China)
and other independent experts

1000 copies
Norwegian Journal of development of the International Science
Iduns gate 4A, 0178, Oslo, Norway
email: [email protected]
site: http://www.njd-iscience.com
–