CC BY
9doi: 10.5862/MCE.66.6
Method of thermotechnical uniformity coefficient evaluation by analyzing thermograms
Метод оценки коэффициента теплотехнической однородности из анализа термограмм
JSC "INSOLAR-INVEST", Moscow, Russia
G.P. Vasilyev, V.A. Lichman, I.A. Yurchenko, M.V. Kolesova,
д-р техн. наук, научный руководитель Г.П. Васильев,
канд. физ.-мат. наук, cтарший научный сотрудник В.А. Личман, ведущий инженер И.А. Юрченко, ведущий инженер-эколог М.В. Колесова,
ОАО "ИНСОЛАР-ИНВЕСТ", г. Москва, Россия
Key words: energy efficiency; the coefficient of thermotechnical uniformity; specific geometrical parameter; thermal resistance; temperature field; temperature; heating; thermogram; building;
Ключевые слова: энергоэффективность; коэффициент теплотехнической однородности; геометрический параметр; сопротивление теплопередаче; ограждающая конструкция здания; температурное поле; теплота; термограмма; здание; строительство
construction
Abstract. One of the primary problems solved for increasing energy efficiency of a building is increasing of the insulating properties of the building envelope. The paper describes the method for determining the coefficient of thermotechnical uniformity by analyzing a thermogram of parts of a building envelope. The method is based on obtaining temperature distribution matrix on the surface of a building envelope fragment by thermography. Two methods for assessing the coefficient of thermotechnical uniformity are proposed. One method for determining the coefficient of thermotechnical uniformity is based on obtaining mean temperature values by numerically processing thermograms, another on numerical integration along the contour lines of temperature curves.
Аннотация. Одной из приоритетных задач, решаемой при повышении энергетической эффективности зданий, является повышение теплозащитных свойств их наружных ограждающих конструкций. В данной работе изложен метод оценки величины коэффициента теплотехнической однородности из анализа термограмм наружных ограждающих частей зданий. Метод основан на получении путем термографирования матрицы распределения температуры на поверхности фрагмента ограждающей конструкции. Предложены два способа оценки величины коэффициента теплотехнической однородности. Один способ оценки величины коэффициента теплотехнической однородности основан на получении усредненных значениях температур путем числовой обработки термограмм, другой - на численном интегрировании вдоль контуров линий температурных кривых.
In Russian Civil Engineering, a dimensionless thermal uniformity coefficient r of the building envelope or its fragment is widely used. This coefficient is determined as r = R0el/R0, where R0el - conventional thermal resistance of a multilayer building envelope having uniformal layers, that is calculated from R0el = Rsi + 2 dk/Àk + Rse (Rsi , Rse - thermal resistances of the internal and the external building envelope layers, m2 • K/W; Ak - thermal conductivity coefficient of the k-th uniformal layer, W/(m • K); dk - thickness of the layer, m ); R0 - reduced thermal resistance of a non-uniformal building envelope, calculated from R0 = (tint - text) • A/Q, where tint, text - internal and external air temperatures, °C; A - surface area of a non-uniformal building envelope or its fragment, m2; Q - total heat flux through this surface, W. The value of the total heat flux Q through the non-uniformal building envelope, containing "thermal bridges" (thermal conductive insertions, gaps, corners, etc.) is determined
Introduction
from numerical computer aided calculation of thermal fields or experimentally. Knowing the value of thermal uniformity coefficient the designer is able to perform quick evaluation of R0 .
Non-uniformal segments of the building envelope (so called "thermal bridges" - joints, dowels, thermal inserts, etc.) take a considerable area of the building envelope that is why they have a considerable impact on insulating properties and durability of the construction [1-3]. Looking for engineering solutions for increasing thermal uniformity of the building envelope is one of the relevant tasks solved by the modern thermal engineering. Both simplified and detailed calculation methods for determining heat transfer (thermal resistance) coefficients of building envelope are well developed and available today [4-6].
Surface thermography is performed according to Russian standards. However, the obtained thermograms due to random impacts of side effects are qualitative. The procedure for thermogram procession has several parameters, providing for alterations in temperature and color match. The core parameters are surface radiation coefficient of the construction - a measure of radiation, emitted by the object, compared to the ideal black body and apparent reflected temperature - a parameter, used to compensate radiation of other bodies [7, 8].
Thermographic analysis is a useful tool for evaluating energy performance a building. It also may be used to verify the presence of air seepage, moisture or water leaks. It is difficult to understand in details all the parameters that influence the thermal scene, that's why infrared thermography is not commonly used for quantitative diagnosis [9, 10].
Most studies involving infrared thermography in building monitoring leads to qualitative results. The technology of infrared cameras has advanced and thermal maps generated by the latest models offer a lot of information in terms of image definition. Better quality of the IR sensors makes the information more accurate. The software for image processing has advanced the most, in terms of graphics management of the pictures: a single shooting makes it possible to make many more evaluations, changing the parameters, to locate homogeneous areas and make very precise calculation of the average values of temperatures. A complete methodology based on the infrared thermovision technique has been developed to reach quantitative data of thermal transmittances of building envelopes [11].
In practice [12] analysis of infrared thermography measurements on a multi-layered wall is a handy tool for thermal monitoring of the walls of a building with tenants. Estimating unknown parameters in the modeling in order to determine the thermal resistance of the wall. Identification is achieved by solving the corresponding inverse problem formulated in the least-squares sense.
In [13] a study is presented to calculate the U-value in building envelopes starting from IR thermography. The absolute deviation between the notional and the measured U-values for IR thermography is found to be at an acceptable level, being in the range of 10 % - 20 %.
The objective of the study presented by [14] is to calculate the sensible heat flux from the exterior surface of buildings using time sequential thermography, based upon an appropriate simplification of the building model. The authors analyzed the thermograms, dividing them into zones by constructing polygons and according to temperature, shape, material and position.
The authors of [15] wished to define a methodology to determine the level of thermal insulation in old buildings through spot measurements of thermal resistance and planar infrared thermography. Currently, the available approaches to evaluate the U-value in existing buildings are: non-destructive methods using heat flux sensors or infrared thermography or destructive methods.
The authors of [16] presented a method to calculate the U-value with IR. To validate their method, they compared IR results with measurements using a heat flux apparatus to obtain the wall heat fluxes and surface temperatures, and U-values were calculated using Fourier's law, taking into account radiation and convection phenomena.
The proposed method is based on the fact that the temperature distribution field, obtained by the thermography, may be corrected ("calibrated") using the values of temperatures, measured by sensors, placed in several specified points. Thus, a qualitative thermogram is transformed into a numerical array of quantative data on temperature distribution across the outer surface of the building envelope.
Methods
Determining the coefficient of thermotechnical uniformity for a fragment of the building envelope
Stage one is selecting (using thermography) of thermally uniformal sections of the building envelope (the reference areas). Next, via contact method, temperature is measured in the selected points of the internal surface of the uniformal fragments of the building envelope, using temperature sensors.
Thermograms are processed using specialized software provided by the thermovisors manufacturers (for example ThermaCAM, FLIR systems). Using the parameters specified above thermograms is calibrated according to the temperature measured by the sensors in the selected points. In this way we achieve equality of temperatures measured by the sensors and obtained from thermography. Then thermograms are digitized and processed using a software tool created in MATLAB. [17-20].
For each building envelope fragment, three types of elements are selected: surface, linear and point. Surface elements are "smooth" parts of the building envelope fragment, having uniform heat flow density. Linear elements - panel joints, window and door esconsons - elements that have one dimension significantly smaller than the other one and insignificant compared to the examined element of the building envelope. Point elements - metallic panel connections, dowels, pins etc. - are elements, with insignificant dimensions of their projection on the surface of the fragment of the building envelope. The coefficient of thermotechnical uniformity of the building envelope fragment containing linear and point thermal non-uniformities is determined from the equation:
r =
i+MZ^z+E^Xk
(1)
, JA J A
where R0/ m2K/W - conditional thermal resistance of the building envelope fragment;
M, = AQj/(tint - text) - specific heat loss through the j-th linear thermal non-uniformity, W/(m • K) (linear heat transfer coefficient);
A Qj - additional heat loss through linear thermal non-uniformity, W/K; tinftext - the temperature of the internal and external air respectively, °C;
Lj/A - ratio of the total length Lj,m of the j-th linear element to the total area of the fragment A/
m2 ;
Xk - specific heat loss through the k-th point thermal non-uniformity, W/K (point or volumetric heat transfer coefficient);
nk = Nk/A - the ratio of average number of point elements Nk to the total area of the fragment A,-
m2.
The first method for assessing the coefficient of thermotechnical uniformity value is based on the average temperature values, obtained by numerical procession of the thermograms:
r =_I__(2)
where,
_mean mean /o\
A _ TQ -TQ (3)
^k mean_t ,
l0 Lext
Tknean _ average temperature of the surface of the building envelope fragment, determined from the thermogram analysis, °C ;
T0iean - average temperature of the uniformal part (reference area) of the external surface of the building envelope, °C ;
bjnean - average width of a linear thermal non-uniformity of the j - kind of the external surface of the building envelope, m;
Amean - average area, with a point thermal non-uniformity of the k - kind, m2.
i
The second method for determining the coefficient of thermotechnical uniformity is based on numerical integration along the contour lines of temperature curves.
The value of V,- W/(m •tf) in the equation (1) is determined from the ratio
_ (4)
(^int-text)
where q^-q^.d/W/m - additional heat loss through linear thermal non-uniformity; qL = /Qdq(x)dx « • - numerical integral of heat flow density along the line, having length of
d,m, perpendicular to the linear non-uniformity; q0,-W/m - value of the heat flow density through the uniformal reference area.
Value of specific heat loss through the point non-uniformity W/m in the equation (1), may be calculated by integrating the volume along the temperature curve contour, obtained from a thermogram:
= AQ = W0raX/2(q)dg = rc-Zsfa (5)
^t (tint-text) (tint-text)
where qt = aext(rext,j - text)/ W/m2 - value of the heat flow density through the external surface of the building envelope fragment; rexM - the temperature in the i-th point of the external surface of the fragment °C, determined from thermograms as well as by direct measurements in several reference points; aext/ - surface coefficient of heat transfer, determined from measurements.
The coefficient of thermotechnical uniformity of the entire building envelope sample is calculated as the weighted average of the fragments:
If beside the temperature sensors a calorimeter for measuring the heat flow is placed in the reference points of the building envelope, it gives an opportunity to determine the surface coefficient of heat transfer of the internal surface allowing determining the value of the thermal resistance of the building envelope fragment.
Results and Discussion.
Let us evaluate the coefficient of thermotechnical uniformity of a panel based on the thermograms shown on Figure 1.
The area of opaque part of the panel equals^ = 5.6 • 2.8 - 1.2 = 14.5 • m2; èi"ean = 1.1 • m; L1 = 5.6 •m; specific geometric parameter is = 0.38 •m-1; average temperatures, obtained from analysis of several dozens of similar areas from different thermograms of the same building: Tmean = -6.5 • oç; Tmean = -7.2 • oç; average external air temperature is text = -7.6 • °C. The coefficient of thermotechnical uniformity of the panels of the examined building equals:
r =-1-= 0 57
r i+i.liL}!.!!-^ °.5/
V14.5/ (-7.2 + 7.6)
Figure 1. Results of processing data from thermograms of a panel building, obtained under field conditions
Let us evaluate the coefficient of thermotechnical uniformity of a balcony based on the thermograms shown on Figure 2. Geometrical parameters of the balcony are A = 3.1 • 2.8 - 2.5 = 6.2 • m2; b™ean = 0.3^m ; L1 = 3.1^m; specific geometrical parameter is L1/A = 0,50 • m-1. Average temperatures, obtained from the thermogram analysis for this building are: °C; T™ = -8.4 • °C; average external air temperature is text = -8.8 • T; the coefficient of
Tmean = -6.9 ■
thermotechnical uniformity equals: r =
- . „ (-6. 9+8. 4)
1+v62/ (-8. 4+8. 8)
= 0.64
Figure 2. Results of processing data from thermograms of a panel building
Let us evaluate the coefficient of thermotechnical uniformity of the place where the window unit fits the panel (Figure 3). Geometrical parameters are = 0.2 • m ; = 4.4 • m;
A = 3 • 2.5 - 1.2 = 6.3 • m2; specific geometrical parameter is = 0.69 • m-1; average external air temperature is text = -6.4 • °C; average temperatures obtained from thermogram analysis are
-2.8 • °C; Tq
—6.0 • T. The coefficient of thermotechnical uniformity equals:
r =
uf44n , (-2.8+6.0) 1+(6:3j-0.1-(-6.o+6.4)
= 0.69
Figure 3. Results of processing data from thermograms of a panel building
Let us evaluate the coefficient of thermotechnical uniformity caused by metallic pins (Figure 4).
is 5; the area of the point thermal non-uniformity is ext = -9.9 • T; average temperature in the
Average number of pins within 1 • m2 pinean = 7.85 •io-3; average external air temperature is t,
vicinity of a point thermal non-uniformity caused by a metal pin is Tlmean = -7.3
of the reference area is T™ean = -9.6 building envelope, containing metal pins equals: r =
1 = -7.3 • T; average temperature T. The coefficient of thermotechnical uniformity of a fragment of a
1+5-7.85-10"
,3 (-7.3+9.6 (-9.6+9.9)
= 0.76
1
T
1
1
1
Figure 4. Results of processing data from thermograms of a panel building
Figure 5 contains thermograms showing linear thermal non-uniformity caused by the gaps between mineral wool panels and point thermal non-uniformity, caused by the metal pins. Values the coefficient of thermotechnical uniformity calculated using equations (5), (4) and (1) equal r = 0.94 (for gaps between panels) and r = 0.81 for metal pins.
Figure 5. Thermograms, showing linear (left) and point (right) thermal non-uniformities of the external surface of samples tested in a climatic chamber
Conclusion
The proposed method is an efficient way to assess the thermal insulation quality of a tenanted building. It allows to:
- assess the thermotechnical uniformity coefficient of a fragment of building envelope under the field conditions;
- locate under the field conditions flawed elements having thermotechnical uniformity coefficient significantly lower than the average value for similar constructions which in turn provides an opportunity to study the insulation quality of said fragments.
The paper presents analytical equations allowing to evaluate the value of the coefficient of thermotechnical uniformity by analyzing thermograms.
A more detailed analysis of design decisions may be carried out based on this method to further improve them and reduce the impact of heat conductive elements on the overall performance of the fragment.
Several examples of processing thermograms obtained in the field conditions using the abovementioned equations are shown.
The described method of determining the coefficient of thermotechnical uniformity and thermal resistance of building envelope fragments tested in a climatic chamber may be used in addition to methods, described in Russian standards, allowing to reduce the testing time. The method is based on combining measuring of the temperature of the uniform building envelope fragments with analyzing temperature fields obtained by thermography.
The coefficient of thermotechnical uniformity is determined by dividing the difference between the internal air temperature and mean temperature of a uniform building envelope fragment by the difference between the internal air temperature and mean temperature of a non-uniform building envelope fragment.
If in addition to the temperature the heat flow density is measured in the chosen spots, it is possible to calculate the surface coefficient of heat transfer for the interior surface. This allows evaluating thermal resistance of the fragment, based on the thermogram.
A software piece for data analysis was developed to determine the value the coefficient of thermotechnical uniformity and thermal resistance of a building envelope fragment based on temperature fields, obtained from thermograms.
Acknowledgement
The Research was conducted by JSC "INSOLAR-INVEST" with the financial support of the Ministry of Education and Science of Russia. Unique identifier of the project RFMEFI57614X0034.
References
Литература
1. Gorshkov A., Vatin N., Nemova D., Tarasova D. The 1. brickwork joints effect on the thermotechnical uniformity of
the exterior walls from gas-concrete blocks. Applied Mechanics and Materials. 2015. No. 725-726. Pp. 3-8.
2. Vatin N.I., Gorshkov A.S., Nemova D.V., Staritcyna A.A., 2. Tarasova D.S. The energy-efficient heat insulation thickness for systems of hinged ventilation facades. Advanced Materials Research. 2014. No. 941-944.
Pp. 905-920.
3. Grinfeldi G.I., Gorshkov A.S., Vatin N.I. Tests results 3. strength and thermophysical properties of aerated concrete block wall samples with the use of polyurethane adhesive. Advanced Materials Research. 2014. No. 941944. Pp. 786-799.
4. Theodosiou T.G., Papadopoulos A.M. The impact of 4. thermal bridges on the energy demand of buildings with double brick wall constructions. Energy and Buildings. 2008. No. 40 Pp. 2083-2089.
5. Renon O. Thermal bridge modeling In energy plus. 5. Building Energy Simulation User News. 2002. Vol. 23.
No. 4.
6. Vatin N., Gorshkov A., Rymkevich P., Nemova D., 6. Tarasova D. Nonstationary thermal conductionthrough the building envelope. Applied Mechanics and Materials. 2014.
No. 670-671. Pp. 365-369.
7. Vatin N., Gamayunova O. Energy efficiency and energy 7. audit: the experience of the Russian Federation and the republic of Belarus. Advanced Materials Research. 2015. Vol. 1065-1069. Pp. 2159-2162.
8. Balaras C. A., Argiriou A. A., Infrared thermography for 8. building diagnostics. Energy and Buildings. 2002. Vol. 34.
No. 2. Pp. 171-183.
9. Ocana S.M., Guerrero I.C., Requena I.G. Thermographic 9. survey of two rural buildings in Spain. Energy and Buildings. 2004. Vol. 36. Pp. 515-523.
10. Bareira E., de Freitas V.P. Evaluation of building materials 10. using infrared thermography. Construction and Building Materials. 2007. Vol. 21. Pp. 218-224.
11. Albatici R., Tonelli A.M. Infrared thermovision technique for 11. the assessment of thermal transmittance value of opaque building elements on site. Energy and Buildings. 2010.
Vol. 42. Pp. 2177-2183.
12. Monchau J., Ibos L., Feuillet V. Diagnosis of insulated 12. building walls using passive infrared thermography and numerical simulations. 7th European Workshop on Structural Health Monitoring. Nantes, France, 2014.
13. Datcu S., Ibos L., Candau Y., Mattel S. Contribution to 13. improvement of building wall surface temperature measurements by infrared thermography. Infrared Physics
and Technology. 2005. Vol. 46. No. 6 Pp. 451-467.
14. Fokaides P., Kalogirou S. Application of infrared thermography for the determination of the overall heat 14.
Gorshkov A., Vatin N., Nemova D., Tarasova D. The brickwork joints effect on the thermotechnical uniformity of the exterior walls from gas-concrete blocks // Applied Mechanics and Materials. 2015. № 725-726. Pp. 3-8. Vatin N.I., Gorshkov A.S., Nemova D.V., Staritcyna A.A., Tarasova D.S. The energy-efficient heat insulation thickness for systems of hinged ventilation facades // Advanced Materials Research. 2014. № 941-944. Pp. 905-920.
Grinfeldi G.I., Gorshkov A.S., Vatin N.I. Tests results strength and thermophysical properties of aerated concrete block wall samples with the use of polyurethane adhesive // Advanced Materials Research. 2014. № 941944. Pp. 786-799.
Theodosiou T.G., Papadopoulos A.M. The impact of thermal bridges on the energy demand of buildings with double brick wall constructions // Energy and Buildings. 2008. № 40 Pp. 2083-2089.
Renon O. Thermal bridge modeling in energy plus // Building Energy Simulation User News. 2002. Vol. 23. № 4.
Vatin N., Gorshkov A., Rymkevich P., Nemova D., Tarasova D. Nonstationary thermal conductionthrough the building envelope // Applied Mechanics and Materials. 2014. № 670-671. Pp. 365-369.
Vatin N., Gamayunova O. Energy efficiency and energy audit: the experience of the Russian Federation and the republic of Belarus // Advanced Materials Research. 2015. Vols. 1065-1069. Pp. 2159-2162.
Balaras C.A., Argiriou A.A. Infrared thermography for building diagnostics // Energy and Buildings. 2002. Vol. 34. № 2. Pp. 171-183.
Ocana S.M., Guerrero I.C., Requena I.G. Thermographic
survey of two rural buildings in Spain // Energy and
Buildings. 2004. Vol. 36. Pp. 515-523.
Bareira E., de Freitas V.P. Evaluation of building materials
using infrared thermography // Construction and Building
Materials. 2007. Vol. 21. Pp. 218-224.
Albatici R., Tonelli A.M. Infrared thermovision technique for
the assessment of thermal transmittance value of opaque
building elements on site // Energy and Buildings. 2010.
Vol. 42. Pp. 2177-2183.
Monchau J., Ibos L., Feuillet V. Diagnosis of insulated building walls using passive infrared thermography and numerical simulations // 7th European Workshop on Structural Health Monitoring. Nantes, France, 2014. Datcu S., Ibos L., Candau Y., Mattel S. Contribution to improvement of building wall surface temperature measurements by infrared thermography // Infrared Physics and Technology. 2005. Vol. 46. № 6. Pp. 451467.
transfer coefficient (U-Value) in building envelopes. Appl. Energy. 2011. Vol. 88. Pp. 4358-4365.
15. Dall'O', Sarto L., Galante A., Pasetti G, Comparison between predicted and actual energy performance for winter heating in high-performance residential buildings in the Lombardy region (Italy). Energy Building. 2012. Vol. 47. Pp. 247-253.
16. Simoes I., Simoes N., Tadeu A. Laboratory assessment of thermal transmittance of homogeneous building elements using infrared thermography. [Electronic resourse]. System requirements: AdobeAcrobatReader. URL: http://dx.doi.org/10.21611/qirt.2014.081 (date of application: 03.10.2016).
17. Vasilyev G.P., Lichman V.A., Peskov N.V., Brodach M.M., Tabunshchikov Y.A., Kolesova M.V. Simulation of heat and moisture transfer in a multiplex structure. Energy & Buildings. 2015. Vol. 86. Pp. 803-807.
18. Vasilyev G.P., Leskov V.A., Gornov V.F., Kolesova M.V., Lichman V.A. Study of heat engineering homogeneity fragments of enclosing structures in the climatic chamber. MATEC Web of Conferences. 2016. Vol. 40. No. of publication 05003
19. Vasilyev G.P., Gornov V.F., Lichman V.A., Yurchenko I.A., A Method of assessing energy consumption of buildings during commissioning. ARPN Journal of Engineering and Applied Sciences. 2015. Vol. 10. No. 15. Pp. 6509-6512.
20. Vasilyev G.P., Leskov V.A., Gornov V.F., Lichman V.A, Kolesova M.V., Dmitriev A.N. Definition of the heat resistance of three-layer ferroconcrete panel. MATEC Web of Conferences. 2016. No. 40. No. of publication 05004
thermography for the determination of the overall heat transfer coefficient (U-Value) in building envelopes // Appl. Energy. 2011. Vol. 88. Pp. 4358-4365.
15. Dall'O', Sarto L., Galante A., Pasetti G. Comparison between predicted and actual energy performance for winter heating in high-performance residential buildings in the Lombardy region (Italy) // Energy Building. 2012. Vol. 47. Pp. 247-253.
16. Simoes I., Simoes N., Tadeu A. Laboratory assessment of thermal transmittance of homogeneous building elements using infrared thermography [Электронный ресурс]. Систем.требования: AdobeAcrobatReader. URL: http://dx.doi.org/10.21611/qirt.2014.081 (дата обращения: 03.10.2016).
17. Vasilyev G.P., Lichman V.A., Peskov N.V., Brodach M.M., Tabunshchikov Y.A., Kolesova M.V. Simulation of heat and moisture transfer in a multiplex structure // Energy & Buildings. 2015. Vol. 86. Pp. 803-807.
18. Vasilyev G.P., Leskov V.A., Gornov V.F., Kolesova M.V., Lichman V.A. Study of heat engineering homogeneity fragments of enclosing structures in the climatic chamber // MATEC Web of Conferences. 2016. Vol. 40. № публикации 05003
19. Vasilyev G.P., Gornov V.F., Lichman V.A., Yurchenko I.A., A Method of assessing energy consumption of buildings during commissioning // ARPN Journal of Engineering and Applied Sciences. 2015. Vol. 10. № 15. Pp. 6509-6512.
20. Vasilyev G.P., Leskov V.A., Gornov V.F., Lichman V.A, Kolesova M.V., Dmitriev A.N. Definition of the heat resistance of three-layer ferroconcrete panel // MATEC Web of Conferences. 2016. № 40. № публикации 05004
Gregory Vasilyev,
+74991440667; [email protected]
Vladimir Lichman, +74991440667; [email protected]
Igor Yurchenko,
+74991440667; [email protected] Marina Kolesova,
+74991440667; [email protected]
Григорий Петрович Васильев, +74991440667; эл. почта: [email protected]
Владимир Александрович Личман, +74991440667; эл. почта: [email protected]
Игорь Андреевич Юрченко, +74991440667; эл. почта: [email protected]
Марина Владимировна Колесова, +74991440667; эл. почта: [email protected]
© Vasilyev G.P., Lichman V.A., Yurchenko I.A., Kolesova M.V., 2016
