DESIGN OF NOISE PROTECTION OF INDUSTRIAL PLANTS- CASE STUDY OF A PLYWOOD FACTORY

Similar documents
Expert paper DESIGNING MAIN FAN NOISE PROTECTION SYSTEM IN MINE JARANDO BALJEVAC

CONFIGURATION OF MODULAR NOISE BARRIERS

NOISE LEVEL CONTROL ON A BUSY INTERSECTION

SOUND INSULATION OF A MECHANICAL WORKSHOP

ACOUSTICAL TREATMENT OF EMERGENCY POWER GENERATOR. Dr. M.Abd El Gawad Saif

The use of layered recycled aggregate concrete barriers in targeting urban noise

A review of the use of different noise prediction models for windfarms and the effects of meteorology

Constructing SoundPLAN Models for Noise Studies That Are Consistent With the HUD Noise Guidebook

ATRIA IN SHOPPING CENTRES, OFFICE BUILDINGS AND HOSPITALS. E.Ph.J. de Ruiter. Adviesbureau Peutz & Associés BV., The Hague, Holland

Safety and Security Glass

Specialist Study on Noise APPENDIX 4 : EXAMPLES OF GOOD OPERATING PRACTICES IN ORDER TO REDUCE NOISE EMISSIONS

Institute of Industrial Science University of Tokyo (Roppongi , Minato-ku, Tokyo 106, Japan)

Applied Noise Control. D. W. Herrin, Ph.D., P.E. University of Kentucky Department of Mechanical Engineering

What is Noise? In daily life, noise means unwanted sound.

THE APPLICATION OF THE NOISE MAPPING TECHNIQUE IN THE COMPLEX INDUSTRIAL ENVIRONMENT

Noise pollution. Sound physics

CONSIDERATIONS REGARDING IMPACT SOUND INSULATION PERFORMANCES OF DIFFERENT FLOOR MATERIALS

Noise Control Case Studies TIM WIENS CONESTOGA-ROVERS & ASSOCIATES

Experimental study of sound insulation performance in the wood-framed building

MCS 020. MCS Planning Standards. For permitted development installations of wind turbines and air source heat pumps on domestic premises

Façade sound insulation of residential houses within Hz

ENVIRONMENTAL NOISE IMPACT ASSESSMENT THE MEADOWS IN HALF MOON BAY PHASE 5

THE REFURBISHMENT OF THE PORTUGUESE PUBLIC SECONDARY SCHOOL NETWORK: AN ACOUSTICAL PERSPECTIVE

Planning. Local. 46 Planning

matteco CONSTRUCTION sound insulation products

Largo & Allegro. Air Compressors. Largo & Allegro

Variations of sound from wind turbines during different weather conditions

Roma to Brisbane Pipeline Dalby Compressor Station Upgrade Environmental Management Plan. Appendix 8. Noise Assessment

Rollair. Air Compressors ROLLAIR 150E AND ROLLAIR 150V-180V-240V

The Clem7 Motorway Tunnel: Mechanical and Electrical Plant Acoustic Design and Performance

ACOUSTIC ENCLOSURES ACOUSTIC ENCLOSURES ACOUSTIC ENCLOSURES

Danish Regulation of Low Frequency Noise from Wind Turbines

Flanking sound transmission in an innovative lightweight clay block building system with an integrated insulation used at multifamily houses

AC-TEC THE COMPANY. Kaplan, Francis, Pelton, Cross Flow and Mini Turbines.

EIS Materials Workshop - Building a Semi-Anechoic Room

Buffer distances for surface roads and elevated highways correlated with pre-existing ambient noise

Wind Turbine Noise, Infrasound and Noise Perception

APPENDIX 6.3-A Noise Modelling Year 2

Ozone Panel 120 Acoustic Assessment and Optimization

Short Communication Infrasound Noise of Natural Sources in the Environment and Infrasound Noise of Wind Turbines

> HGW WATER - WATER AND BRINE - WATER HEAT PUMPS FOR INDOOR INSTALLATION. Available range. Unit description. Options. Accessories

ΗERMES GROUP SERIES AUTOMATIC HORIZONTAL BEAM SAWS. Cutting edge technology since 1977

Stationary Noise Assessment. Car Dealership Strandherd Drive. Ottawa, Ontario

Rumster Wind Energy Project

Highway Traffic Noise

A CASE STUDY ON COST/BENEFIT ASSESSMENT OF ROAD TRAFFIC NOISE AMELIORATION WITHIN AND OUTSIDE THE ROAD RESERVE

Acoustic Consulting Australia PtyLtd. I Consultants on Noise and Vibration P0 Box 332

DEVELOPMENT AND STUDY OF ANTI-NOISE BARRIERS BASED ON ALTERNATIVE RAW MATERIAL SOURCES

NOISE IMPACT ASSESSMENT FOR SITE PLAN APPROVAL 47 HAVELOCK STREET OTTAWA, ON. REVISED with ADDENDUM for ROOFTOP OUTDOOR LIVING AREA

Noise Assessment Report

The specimen B W was identified by the client as a Skyfold Classic operable partition, with panels, seals, and clearances in configuration 4E.

Sound Insulation of Residences and Other Buildings Exposed to Aircraft Noise

Building & Ductwork Airtightness Standards

Noise Control Study. for. Avalon Public School Portobello Boulevard at Stormwind Avenue in Ottawa, Ontario. Revision 1

INNOVATIVE AND SUSTAINABLE COATING SYSTEM: PRELIMINARY ACOUSTIC EVALUATION RESULTS

HIGH PERFORMANCE ACOUSTIC INSULATING MATERIAL FOR A QUIETER ENVIRONMENT

ENVIRONMENTAL NOISE CONTROL FOR SEMICONDUCTOR MANUFACTURING FACILITIES

BRITISH COLUMBIA MINISTRY OF TRANSPORTATION AND INFRASTRUCTURE POLICY FOR ASSESSING AND MITIGATING NOISE IMPACTS

COMPARISON OF MEASURED AND PREDICTED SOUND INSULATION FOR A THERMAL RETROFITTED BUILDING

inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE

ACOUSTIC ASPECT OF INNER PARTITION WALLS MATERIALIZATION IN STUDENT HOSTELS ARCHITECTURE

Wind Turbine Noise Assessment in Ontario

SLIDING WEAR BEHAVIOR OF Zn-Al ALLOYS IN CONDITIONS OF BOUNDARY LUBRICATION

SOUND INSULATION DESIGN OF MODULAR CONSTRUCTION HOUSING

VTT Technical Research Centre of Finland VTT GROUP 2017 Place for a photo (no lines around photo) Mikko Nyman VTT Expert Services Ltd

Appendix R Effect of the Green Deck on Local Noise Environment

AN INDIRECT-FIRED ROTARY SUGARCANE BAGASSE TORREFYER: DESIGN AND PERFORMANCE

Environmental Noise Assessment 1020, 1024, 1028, 1032 & 1042 Sixth Line Oakville, ON

Project Report. EFP-06 project Low Frequency Noise from Large Wind Turbines. Results from Sound Power Measurements. Client: Danish Energy Authority

Roadway Traffic Noise Assessment Montreal Road. Ottawa, Ontario

Aluminum Secondary Glazing Specification Guide June 2012

Research of transmission loss properties of thin layered sound reduction systems for building partition elements

Moduline. Acoustic Enclosure System. Rugged Noise Control Structures Using Acoustically Rated & Field Proven Demountable Components

Acoustic Louvers A Complete Range of Certified, High-Performance Acoustic Louvers to Solve Diverse Environmental Noise Pollution Problems

8. Appendix F Acoustic Assessment

August 2008: Issue 1. Non - Residential New Build Built - up Metal Walls

Enhancing Acoustic Performance through the use of Window Shade Fabric White Paper

Noise Feasibility Study EMGO (North Oakville I) Ltd., Town of Oakville, Ontario

Burnt Hill Wind Development

Rectangular sound attenuator

Noise Feasibility Study Proposed Residential Building 103 Dundas Street West Oakville, Ontario

ENGLISH SUMMARY. Rakenteiden Mekaniikka, Vol. 24 No 1, 1991 ON THE EFFECT OF NOISE AND VIBRATION TO MAN. Antero Honkasalo

Stationary Noise Assessment. Orléans Gardens. Ottawa, Ontario

inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE

FLOOR IMPACT SOUND INSULATION PERFORMANCE AND IT S REGULATION OF MULTI-FAMILY DWELLINGS IN JAPAN

Evaluating Appropriate Residential Ventilation Strategies in Dense Urban Environments and the Challenges for Passive Design

REP Fotenn Bank Phase 1 Feasibility Noise Study with Building Component AmendmentPage 1 of 21

PUMPS/ COMPRESSORS. Mounting, Levelling and Vibration Control of: Pumps. Compressors POWER GENERATION - PUMPS & COMPRESSORS

A CONCEPT FOR USING RECYCLED RUBBER GRANULES IN NOISE REDUCTION CONCRETE S PANELS

Stationary Noise Feasibility Study. The Shops of Tenth Line. Ottawa, Ontario

MAIN STREET PRECISE PLAN ENVIRONMENTAL NOISE STUDY REDWOOD CITY, CALIFORNIA

PUMPS/ COMPRESSORS. Mounting, Levelling and Vibration Control of: Pumps. Compressors POWER GENERATION - PUMPS & COMPRESSORS

Noise measurement and mitigation for urban building foundation excavation

Acoustics/Noise Abatement

NEC Acostix. Sound Attenuation Calculations and Design Software. General Information

Stationary Noise Feasibility Assessment. 315 Chapel Street Ottawa, Ontario

Mansfield University Dormitories Mansfield, Pennsylvania

Identifying impacts of amplified sound in commercial spaces below residences in mixed-use buildings

inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE

Stack Parking Surface Mounted Stackers

Transcription:

DESIGN OF NOISE PROTECTION OF INDUSTRIAL PLANTS- CASE STUDY OF A PLYWOOD FACTORY Branko Radičević, Zoran Petrović, Slobodan Todosijević, Zvonko Petrović University of Kragujevac, Faculty of Mechanical Engineering Kraljevo, Serbia, radicevic.b@mfkv.kg.ac.rs Abstract - The paper presents activities of reduction of noise emitted by plywood factory facilities. Dominant noise sources in the factory, which caused noise levels above limits, were identified. Based on calculation of acoustic fields performed by original software package, a noise map of the endangered area was created, and a noise system was designed. The noise system consists of two sound barriers, insulation of noise sources and changes of organization of the technology process. Application of the noise system led to reduction of specific noise levels for 14 db. Noise levels at a control site in residential areas during day do not overcome limits anymore. Keywords: industrial noise, noise maps, acoustic calculations, noise, noise barriers 1. INTRODUCTION By the example of a chip factory we presented the activities needed to be realized in order to successfully protect the residential area located near the factory. This and similar residential zones naturally sprang up after the factories had been built. Furthermore, the residential zones located in, or closely to the industrial area are difficult to classify into some of the legitimate acoustic zones. This leads to the introduction of a new acoustic zone for the residential areas near an industrial zone. Examples of similar solutions are known in some European countries. In this way the subjects of environment noise would be allowed to mutually cooperate and provide coordination in carrying out measures and decisions. Total amount of noise level in residential area near the chip factory was 76 db (A), which far exceeds the noise indicator values in an environment for the V acoustic zone. There are two causes for the far exceeding. The first is the short distance of the factories, and the second is the frequent city traffic artery between factories and residential blocks. The noise polluted residential zone stretches parallel to the factory halls on the length of around 170 m. Approximate distance between the house facades and the nearest noise sources is around 35 m. Key problems in searching the most efficient method to protect from noise are: an enormous height of some of the factories (up to 25 m), great number of noise sources, distance of the houses and densely populated area. We measured the noise levels of all sources, and then established the dominant noise sources. By measuring the levels of all noise sources, adopting the dominant sources, establishing the type of noise source and weakening of noise through the space by building sound walls, a mathematical method has been formed to establish levels of noise and select the most suitable variant of noise considering the possibilities of technical realization. 2. SOURCES NOISE 2.1 Types of Noise Sources Sources of the noise in factory according to its location can be conditionally divided into three groups: The main noise sources are: the main compressor station and production power processing part which includes three subsystems: a line for sifting raw and dry chips, a line for measuring dust, a line for transport and sifting chip shavings chips preparations (a bridge crane, log transporter, mill) Equipment and machines used in the process of chip production fall into group of stationary sources that generate high levels of noise. In the part of the factory that represents the dominant source of noise, three subsystems can be noticed: a line for sifting raw and dry chips, a line for measuring dust, a line for transport and sifting chip shavings. Fig. 1 Noise sources in the wood processing plant Raw chip is transported by rake and worm transporters from the silo to the dryer of outer layer (DS) and to the dryer of inter layer (MS). From the drier chip is pneumatically conveyed by the ventilator to the cyclone separator. First the coarsest chip fraction is pneumatically conveyed back to the mill (PSKM). In the mill chip is being ground fragmentized and pneumatically conveyed by the ventilator to 71

the cyclone separator from where, by the halt and systems of worm transporters is directed to the sifting grate the sieve. The second and the third fraction are pneumatically conveyed to the drier chip silos of outer layer (DS) and inter layer (MS). The fourth chip fraction is the finest one (the dust) pneumatically conveyed by the compressor to the dust silo. 2.2 The Type of Noise In search for successful model of noise control, evaluation of noise levels and taking adequate measures and methods for reduction of the noise, it is necessary to have the data about the characteristics of noise itself as exact as it is possible. These data are obtained by measuring characteristic noise values in frequency, amplitude, and time range. In order to establish precise noise characteristics, a series of noise measuring have been performed during day and night. Major mechanisms for generating noise in this case can be divided into two groups: mechanical (ventilators, compressors, electro motors, gear transmission, belt transmission, transporters, bearings, and other rotational components), and aerodynamic which are result of pneumatic chip transport. Machine tools of this type generate very high noise levels mainly characterized as tonal, with components ranging from medium to high frequency proportional to rotational speed of mechanical components. Aerodynamic noise is a result of chip motion under the influence of air circulation in pneumatic pipelines which connect silos of great heights (some may even reach 30 m height). Since the shavings in the pneumatic pipelines are moved with high speed through different cross sections and curves, it comes to turbulence, which generates wide belt noise. 2.3 Noise Level The factory works around the clock with no pauses. Table 1 Dominant noise sources in the factory No. Mark Noise level [db(a)] S 1 Rake transporter of inter layer 81 S 2 Gear of rake transporter 81 S 3 Sifting grate sieve 81 S 4 Dry chip Ventilator of inter layer 100 S 5 Dry chip Ventilator of outer layer 97 S 6 Vertical worm transporter 89 S 7 Horizontal rake transporter 89 S 8 Inter layer halts 89 S 9 Ventilator for extra chip drying 94 S 10 Rake transporter Z-silos 86 S 11 Log transporter 83 S 12 Crane (engine) 85 S 13 Reverse pipeline from the hall 85 S 14 Ventilator power 90 kw 87 S 15 Ventilator 81 S 16 Rake transporter 4 kw 80 The first base source is formed by the means of 29 sources divided in three already mentioned groups. Noise levels of particular sources are measured at very close distance (up to 0.5 m), in order to avoid influence of other sources. It should be pointed out that there are different noise sources in the factory which levels are less than 80 db. Only noise levels from 80dB to 100dB are taken into consideration. As a criterion for the choice of dominant sources, noise levels formed by the separate sources are taken, as well as their distance from residential area, the height of the source and their possible shelter from other objects (Table 1). 2.4 Noise Frequency Spectrum In addition to measuring equivalent noise levels for characteristic sources, frequency noise analysis has been conducted. Localization of the noise source: chip preparation and the processing, frequent noise characteristic typical for the factory region is described by the means of 1/3 octave filters. The noise located on the chip preparation has an average equivalent level from 92 db (A) and is wide belted, because it has an approximately even distribution of acoustic energy in a wider frequency spectrum. In the Fig. 2, the description of 1/3 octave frequency spectrum is given. The noise in the processing part of the factory is also wide belted. In the Fig. 3. The frequency characteristic on the platform of the drier is described, where can be seen that noise level is a bit higher regarding the lower and medium frequencies. Fig. 2 Frequency characteristic plywood preparation Fig. 3 Frequency characteristic drier platform Compressor station as third location of noise source isn t separately analyzed in terms of frequent analysis because the compressors are in objects made of bricks. Measuring has shown that is enough to keep windows on the compressor station closed. Under this treatment, the compressor station has no influence on the increase of noise level in the environment. 3. NOISE PROTECTION SYSTEM In order to enable efficient design of noise system, it is necessary to form a mathematical model, which can be used for calculating a noise map for total area. The major advantage of forming the mathematical model is a realization 72

of variant acoustic solutions. These solutions, along with other limitations are a foundation for multiple criteria optimization when making choice about specific solutions for noise system. 3.1 Noise at the Checkpoints By the series of measuring in the residential area, the checkpoint has been detected, i.e. the point where the level of noise is the highest in the entire residential zone. As the process of production ceaselessly continues, solution to the noise lies in exploring the system, which at the checkpoint in front of the residential block (in further analysis marked as KT1) provides the level of specific factory noise of 55 db (A) which is a tolerable value of noise indicator in environment during night. Table 2 Noise level at control points Measuring level [dba] [dba] Residual noise Total noise level Tolerable point noise level [dba] day night day night day night КТ 1 63 53 76 73 65 55 If for the existing noise level we adopt the noise level established by measuring noise near the street, which separates the factory and residential area during night, we can calculate a specific noise level of the factory at the measurable checkpoints, i.e. in the residential area. 3.2 Noise Map A necessity for the noise level to be calculated on various points and a need for parallel observation of different sources and barriers, led to the usage of software as a solution to the calculation of noise field. (1) Weakening of the noise level is a result of the correction caused by: C D - divergence of sound waves, C A air absorption, C B barrier, C G terrain absorption, C Z green areas. For sofware application the noise sources (position and level of the noise), position of the barriers, and the position of the checkpoints need to be defined. As a final resul we get noise level in the knots of the web. The step in the web may be taken arbitrary, but usually the web 10x10 m is used as a noise map, Fig. 4. Noise map a desired state. 3.3 Methods for Noise Protection System Project In order to successfully project noise system we need to analyse active and pasive methods of through the effects provoked by the specific variants in the environment. Exact answers in terms of choice of the most efficient and most economic variant of noise, coordinated with technological and building projects, can be only given when analyzed parallelly with technoeconomic analysis of various variants. Fig. 5 Barrier and its base construction Fig. 4 Noise map a desired state The software enables calculation of noise field produced by multiple point noise sources, and their stretching is limited with several line barriers. A model for determination of noise field is based on the spatial model of the field used for the calculation of the noise level, the model of the source, and the model of the barriers. Noise field is described by the noise levels calculated in the points, which are part of the equidistant areas parallel to the basis of the field. An algorithm has been formed to calculate the weakening during the spreading of the noise: C = C D +C A +C B +C G +C Z (1) Fig. 6 Chosen panel and its protective level When deciding about alternative solutions, we considered possibilities of technical realization of suggested variants. As a solution of the problem a combination of active and passive noise has been suggested. Building of the barriers 6 m high and 60 m long and active on the noise sources themselves, with height equal or higher from the barrier height, both represent the most rational solution to protect from the noise coming from the factories in the example mentioned above. 73

example of vertical worm transporter (pict.9) gives lessening of noise level immediatelly behind the panel for 20 db. The same goes for the small barriers made of the same panels and set in front of dominant noise sources. (pict.8) Fig. 7 Isolating power of the panel To build a wall sound we took the panel 80 mm thick (pictures 6, 7), with one perforated surface (32,7%) facing the sources of the noise. Besides its great isolation power, the panel is also fireproof, which is very important here considering the fact that we are talking about chip factory, where the chances for fire are big. Fig. 9 Vertical worm transporter before and after sound 3.4 Protection from Separate Noise Sources Besides the two barriers, 6 m high and 108 m long, noise from separate dominant sources has been carried out. A criterion for establishing priorities in order to realize noise from separate sources primarily was acoustical. The from the noise sources, which produce the highest level of noise, or sources at the shortest distance to the residential area or sources located on the great height, has been carried out. Barriers made of ribbed sheet metal give effect of lessening noise levels (immediatelly behind the barrier) for 8 db. This solution is displayed at the example of inter layer halt of dry chips. Separate sources included by the system are: Vertical worm transporter Ventilator for chip drier with interlayer Inter layer halts Ventilator for chip drier with outlayer Ventilator for extra drying of the chip Some examples of the noise are presented in the figures 8, 9, 10 and 11. Fig. 10 Inter layer halts before and after sound Fig. 8 Ventilator for chip drier with interlayer after sound Fig. 11 Ventilator for extra drying of the chip before and after sound Above mentioned noise sources are protected by a shield, as displayed on the picture which represents vertical worm transporter, or is achieved by putting small barriers. Barriers made of ribbed sheet metal combined with rubber acoustic curtains (immediatelly behind the barrier) give effect of lessening noise levels for 6dB. (fig.11) Protection made by shield of absorptive panels (the same panels as those used for the barrier), what can be seen on the 74

4. EFFECTS OF NOISE PROTECTION SYSTEM Carrying out the measures, relevant noise level of the factory has been realized and it does not exceed tolerable limit of 65 db during daily working conditions. Further measurs are to be taken in accordance with recommendations given in this paper, in order to reach tolerable limit of 55 db for night. Effects of lessening of noise levels on the checkpoint are displayed in the following table: Table 3 Effect of reductionof total noise level in the factory after realization of projected measures of noise Tipe of Protection Two bariers (6x108 m) Five dominant noise sources Noise level db(a) Before After Reduction of Total Noise Level db(a) Separately Cumulatively 76 69 7 7 69 66 3 10 Total noise level in the residential area on the checkpoint is lessened for 10 db (A), starting from the beginning of realization of projected measures for noise. This means reduction of the total noise level of the factory for 3.16 times, i.e. 316.2%. Table 4 Effect of reductionof specific noise level in the factory after realization of projected measures of noise Tipe of Protection Two bariers (6x108 m) Five dominant noise sources Noise level db(a) Before After Reduction of Specific Noise Level db(a) Separately Cumulatively 76 69 8 8 69 66 6 14 The specific noise level on the checkpoint has been lessen for 14 db(a), starting from the beginning of realization of projected measures for noise. This means reduction of the specific noise level of the factory for 5 times, i.e. 501.2%. 5. CONSLUSION AND DIRECTION OF FURTHER PROTECTION ACTIVITIES By measures of noise, a specific factory noise level of 62 db has been realized. Reducing the noise level to 14 db has considerably improved acoustic comfort in the residential zone near the chip factory. This level of noise is in a range of tolerable values for noise indicator in the environment, refferring to time period during a day. Since the factory works around the clock, it is necessary to take extra measures in order to protect from noise and thus reach tolerable noise level during night. The directions of further actions is difficult to realize at this point, because there is a large number of point noise sources with levels around 80 db and they are located on great heights. Moreover, there are several rake transporters, which practically represent line noise source and they need to be isolated by entire length. Some of the possibilities are connecting elastic pipelines to the basis of huge electromotors and fans. Above presented example confirms that by design the noise system great effects can be achomplished in terms of reduction of noise level in the environment. 6. REFERENCES [1] Petrović, Z., Radičević, B.,Šoškić Z., Bjelić, M. The from the noise of industrial plants - Noise and vibrations, XXI Conference with international participation, Tara 2008. [2] Petrović, Z., Radičević, B., Bjelić, M., The proposal for the solution for from noise ŠPIK IVERICA FANTONI GROUP - Ivanjica, University of Kragujevac, Faculty of Mechanical Engineering Kraljevo, 2007. [3] M. Kolarević, Z. Šoškić, Z. Petrović, B. Radičević: Noise Protection In Urban Environment- Description Of A Project, Mechanics, Transport, Communications, Academic journal, Todor Kableshkov University of Transport, Sofia, ISSN 1312-3823, issue 3/2011, pp.iv- 69-IV-78. [4] Praščević, M., Cvetković, D., Environmental Noise, University of Niš, Faculity of occupational Safety Niš, 2005 [5] B. Kotzen, C. English: Environmental Noise Barriers, A guide to their acoustic and visual design, Second edition, Taylor & Francis e-library, London and New York, 2009. [6] Uzunović, R., Protetction from Noise and Vibration, LOLA Institute, Belgrade,1997 [7] Veličković, D., Noise and Vibration 2, University of Niš, Faculity of occupational Safety, Niš, 1990 [8] Standards: SRPS ISO 1996-1 and SRPS ISO 1996-2 Acknowledgement: The paper is a part of the research performed within the project TR37020. The authors would like to thank to the Ministry of Education and Science of the Republic of Serbia for supporting this research. 75