Stroboholography for the analysis of vibration modes of mechanical systems/Mechaniniu sistemu virpesiu komponenciu analize strobo holografijos metodu.

Vasiliauskas, R. ; Palevicius, A. ; Busilas, A. 等

1. Introduction

The holographic interferometry is used in vibration and modal analysis, structural analysis, composite materials and adhesive testing, stress and strain evaluation, and flow, volume/shape, and thermal analysis. All these applications derive from one or more of the four basic methods of the applied holographic interferometry: real-time, double exposure, time-average holography and stroboholography.

The first three holographic methods for investigating static and dynamic processes in various mechanical, electromechanical and microelectromechanical systems, the methodologies of the analysis of obtained interferograms are widely described in the literature, including [1-10]. Effectiveness of one or another method in solving specific problems is preconditioned by several factors--the form of analysed dynamic process (harmonic vibrations, multicomponent vibrations, random vibration processes, impact processes, etc.) the type of the desired information (qualitative or quantitative), conditions of the carried out experiment. In Fig. 1 the examples of interferograms, that are obtained using the methods of holographic interferometry--real-time, double exposure, time-average holography--for analysing and improving operation of the above systems and design solutions are presented.

In Fig. 1, a, b the interferograms of piezoelectric motion transducers obtained by the average-time method for various regimes of operation are presented. The transducers are excited by harmonic signals of different frequencies, so it is important to visualize the nature of the vibrations at different time moments of the vibration period.

The reliability of computer operation largely depends on hard disc housing deformations at different temperature regimes. The visualization of hard disc housing deformations, qualitative and quantitative assessment of the deformations is actual problem at development and design phases. The presented in Fig. 1, c interferogram of hard disc is obtained by the method of double exposure, under the changing temperature regime.

The object shown in Fig. 1, d presents the advances in the design of vibromotors using holographic interference methodologies [1]. The vibromotor itself is interpreted as a mechatronic system consisting of a piezoelectrical actuator and the driven object. The principle of design covers the areas of the application of standing deformation wave energy and propagating wave energy. The specifics of vibromotors are achieving high precision levels in the motion of the driven object [11]. The only way for increasing the accuracy levels of vibromotors is exploiting specific geometric shapes of materials for the actuator and especially in the contacting zones. The method of laser holographic interference is used for determination and control of working characteristics of vibromotors (eigenshapes at appropriate resonance frequencies transfer of waves in the contact zones, etc.). The experimental data is later used as source data for numerical optimization of the design methods of vibromotors. The methods of holography interferometry are used for the extraction of objective experimental characteristics of the working regimes, what later gives the guidelines for optimisation of vibromotors' design. The applied methods of laser holography have to be different for the two working regimes. In the first case, when the output (driven) part is motionlessly fixed, it is possible to use the laser with continuous beam. In that case the method of realtime data processing is used, the motion of objects is averaged in time and time modulation is used (real-time stroboholography). That method can be applied, as the driven part is motionless, and the processes taking place in the contact zone between input and output links do not give rise to longitudinal motion, and do not affect the optical transformation. In case when the driven part is performing longitudinal motion (each of them can differ by its character) the holography method must use the following procedures--impulse holography with compensation of motion of the analysed object

Vibrational dispergator (holographic interferogram of actuator is presented in Fig. 1, e) consists of concentrically located cylinder shaped ultrasound processing chambers [1]. The cylindrical walls of the chambers are formed of piezoelectric actuators fed by the power sources. When the liquid is fed into the ultrasound processing chambers and the piezo-actuators start vibrating at their resonance frequencies, radial standing waves are generated around each actuator. The resonance frequency for the actuators is different and depends on radii of the actuators. Thus the flowing liquid is processed with higher frequency ultrasound every time it passes through the next processing ring. At the same time the high peak constant current actuators radiate high energy impulses to the processed liquid. Thus the whole process of dispergation turns to be more effective. The quality regulation of the dispergation may be performed by control of the level of constant voltage peak. It was proven experimentally that the effectiveness of the process is increased 3 times.

[FIGURE 1 OMITTED]

Holographic interferometry is a powerful tool for the analysis dynamic of microelectromechanical systems [5, 8, 12 - 14]. It is a nondestructive whole field technique capable of registering microoscillations of MEMS components. There exist numerous methods used for the interpretation of patterns of fringes in the holograms of analyzed objects [15 - 17]. Unfortunately, sometimes straightforward application of these motion reconstruction methods (fringe counting technique, etc.) does not produce acceptable results. A typical example is holographic analysis of a micro-electromechanical switch which is described in [5].

The cantilever vibrations of microelectromechanical switch produce nonlinear effects in microstructure. Holographic interferogram of vibrating surface of such cantilever can not be evaluated straightforward using characteristic function for distribution interference fringes in case when time-average method for recording hologram is used. The characteristic function of the distribution of interference fringes of the vibrating surface should be chosen respectively to the nature of nonlinearities taking place in working regimes of microelectromechanical system.

A direct method to reduce the recording time during the pulsed illumination in stroboscopic holographic interferometry is suggested. The technique involves increasing the recording time while the object is in static position and decreasing the recording time during the pulsed illumination accordingly while the object is vibrating. The approach is applied to double-exposure and triple-exposure cases of stroboscopic holographic interferometry, where it is found that considerable reduction in the recording time is possible.

An experimental evaluation of a stroboscopic technique as applied to holographic interferometry is described. This technique is applicable in the analysis of vibrations with high amplitudes and complex wave forms and in the investigation of phase-dependent and nonsymmetrical effects.

2. Theoretical background of the methodology

For the analysis of vibrational processes of the objects, which are excited by harmonic and nonharmonic signals it is very important to carry out the qualitative and quantitative analysis of vibrations distribution on the surface within a single period [18, 19]. In this case, it is reasonable to use the stroboholography method. The paper presents the methodology of strobo pulse formation and synchronization with the corresponding vibration phase, allowing to increase the information generated by holographic interferograms and accuracy of the analysis. In order to ensure optimal hologram exposure time (not more than 1 to 1.5 min) and information in the obtained interferograms it is necessary to select vibrations period filling by strobo pulses of the lowest possible duration given by the expression of [T.sub.1]/k, where [T.sub.1] is period of excitation signal, which depends on the analysed object, and the value of which can be practically from 0.1 to 5 [micro]s k is the number of strobo pulses.

[FIGURE 2 OMITTED]

In Fig. 2 the structural scheme of excitation signal formation, strobo pulse generation and synchronization with the desired vibration phase is presented. The signal from the generator of harmonic signals 8 or nonharmonic signal generator 12, through commutator K is transferred to the rectangular pulse forming scheme 11. The formed rectangular pulse signal the period of which corresponds the period [T.sub.1] of the excitation of signal of the analysed object is fed into the pulse splitting, counting and strobo pulse forming block 7. From this block, the signal is input into strobo pulse commutation and separation block 10 and the excitation signal forming block 6. The separated pulse, the duration of which may vary from 0.1 to 5 [micro]s, is formed with the transfer characteristics ensuring the period of pulse fronts not higher than 10-15% of the minimum strobo pulse duration by strobo pulse forming scheme 9. Strobo pulse amplifier 5 amplifies the strobo pulse up to the voltage needed to control electro--optical modulator.

Strobo pulse synchronization with the corresponding phase of excitation signal is illustrated by the diagram in Fig. 3

[FIGURE 3 OMITTED]

The diagram shows the selection possibility of two strobo pulses: the first strobo pulse corresponds the start time of excitation signal; the second one is shifted in phase with reference to the first one. This methodology of stobo pulse formation has a number of advantages:

* the synchronization between strobo light pulse and the vibration phase of the analysed object is achieved throughout all the hologram exposure time;

* several strobo pulses of the necessary duration can be used for vibration analysis during one period of vibrations;

* strobo light pulse duration can be changed from [T.sub.1]/ up to [T.sub.1]/k.

3. Analysis of the obtained interferograms

The case is analysed where an object, any point of the surface of which in the direction of normal component vibrates in accordance with the law as described by the

equation:

W (x, y, t ) = A ( x, y) sin[[omega].sub.1]t (1)

and holographic interferogram is recorded with the help of short laser light pulses, which are synchronized with the maximum amplitude of the vibrating surface position.

The frequency of laser beam pulses, illuminating the object and the hologram at these moments in time is equal to the frequency of the object vibration [f.sub.1] = [[omega].sub.1]/2[pi], while their minimum duration is [[pi].sub.i] = [T.sub.1]/[k.sub.1], where [T.sub.1] is the period of vibrations of the object under investigation; [k.sub.1] is the number of rectangular pulses forming the excitation signal (filling in the vibration period). Each of these pulses can be electrically selected as the strobo pulse. Then the center of strobing pulse will at the time instant [T.sub.1]/4+m [T.sub.1], where m = 0, 1, 2, 3 and the light intensity distribution in the obtained hologram can be described by the following equation [17]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)

where [bar.I](x,y) is intensity of recovered image of point (x, y) against the intensity at the point that during all the period of vibrations was in zero position; [J.sub.0], [J.sub.n] are Bessel function of the first type of order 0 and n accordingly; K is sensitivity vector of holographic interferometer at the point (x, y) [2]; A is amplitude of the object surface vibration at the point (x, y).

With this method of the object and hologram illumination, the reduction of strobo pulse duration significantly reduces the influence of the object displacements on the recovered image quality, but significantly prolongs the duration of the hologram exposure. Assuming [k.sub.1] = [infinity], and then sinn[pi]/[k.sub.1]/n[pi]/[k.sub.1] = 1, light intensity distribution is:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)

Taking into account [(1).sup.p] = cos p[pi] = cos 2 p [pi]/2 and [(-1).sup.p] / i = i sin[(2p - 1) [pi]/2] and using the features of Bessel function [17]: [J.sub.0](z) + 2 [[infinity].[summation over (p=1)] [J.sub.2p] (z)cos 2p [pi]/2 = cos(z) and 2i[[infinity].[summation over (p=1)][J.sub.2p-1](z)sin[(2p - 1)[pi]/2] = i sin z, finally the Eq. (3) can be expressed as:

[bar.I] = [[absolute value of cos(KA) + i sin(KA)].sup.2] = 1 (4)

Thus, recording the hologram with the help of very short light pulses, which are synchronized with the phase of harmonic vibrations the influence of the object's vibrations on the intensity distribution field of the reconstructed image is eliminated. However, with the very short strobo pulse duration, the exposure time of the hologram significantly increases. In order to determine the possibilities to eliminate these shortcomings the calculation of intensity distribution [bar.I] =f(A/[lambda]), at different values of [k.sub.1] is be made according to formula (2), when illumination and observation directions are perpendicular (normal) to the object's surface.

In Fig. 4 the intensity distribution graph of the function [bar.I] = f(A/[lambda]) for various strobo pulse values is presented. It is known that human eye does not distinguish between light-intensity fluctuations, which are equal to or greater than 0.82. Taking this fact into account, it is possible to set such strobo pulse duration that for surface points of the object, vibrating with maximum amplitude the intensity [bar.I] value is not less than 0.82. In this case, visualizing the reconstructed hologram of vibrating object, it will not be different from the hologram of nonvibrating object. From Fig. 3 it can be concluded that for low amplitudes of the tested object vibrations, which are usually investigated by the methods of holographic interferometry the above conditions are fulfilled at very short duration strobo pulses. Considering optimal exposure time of the hologram, informativeness of the interferograms, parameters of electronic part the number of strobo pulses [k.sub.1] is selected to be 20.

[FIGURE 4 OMITTED]

When surface of the object at the same time is excited by several harmonic components, using the method of stroboholography makes it possible to investigate the influence of separate harmonic components of vibrations on the analysis of surface vibrations. By selecting the duration of strobo pulses in such a way that for surface points vibrating with a maximum amplitude [bar.I] = j(A/[lambda]) the value would be not less than 0.82, the image of reconstructed hologram practically will be not different from the image of nonvibrating object.

Under this conclusion the case where surface of the object at the same time vibrates by several harmonic components [8], each corresponding to its three-dimensional form of vibrations will be examined. In this case, position of each surface point of the object during the period can be described by the equation:

W (x, y, t) = [A.sub.1] (x, y) sin [[omega].sub.1]t + [A.sub.2] (x, y)sin ([[omega].sub.2]t + [theta]) (5)

where [A.sub.1] (x,y)sin[[omega].sub.1]t and [A.sub.2](x,y)sin([[omega].sub.2]t + [theta]) are amplitude of vibration components of the point (x, y) the corresponding frequencies of which are [f.sub.1] = [[omega].sub.1]/2[pi] ir [f.sub.2] = = [[omega].sub.2]/2[pi]; [theta] is phase shift of the component.

If interferogram of such a surface is recorded by strobo pulses with the duration of [[tau].sub.i] = [T.sub.1]/[k.sub.1] and the frequency [f.sub.1], and synchronized with the vibrations amplitude position of the selected component [A.sub.1](x, y)sin [[omega].sub.1]t so taking into account the information in Fig. 4 always the strobing pulse duration can be chosen such a one that the indicated vibration component practically will have no influence on the image of reconstructed interferogram. Then the approximate expression of characteristic function, describing the intensity distribution in the reconstructed interferogram will be conditioned only by the vibration component [A.sub.2]sin([[omega].sub.2]t+[theta]):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)

[T.sub.e] is exposure time of the hologram, [phi] is phase of the light wave dissipated by the object's point (x, y) at the hologram's plane. As the sensitivity vector K of holographic interferometer at the point (x, y) is not a function of time the phase of light wave can be expressed as [phi] = [[phi].sub.0] + + K[A.sub.2]sin([[omega].sub.2]t+[theta]), where [theta] is phase shift of the second harmonic component and substituding in to Eq. (6) the following is obtained

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (7)

Expansion of the integrand expression (7) by Besel function gives:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (8)

For the component [A.sub.2]sin([[omega].sub.2]t+[theta]) the phase of strobo pulse will be shifted by the value 5, where 0 < 5 < 1 is a part of the period [T.sub.2] in which the strobo pulse is located. In this case, integral limits at the hologram exposure will be from [delta][T.sub.2] - [[tau].sub.i]/2 up to [delta][T.sub.2] + [[tau].sub.i]/2 and taking into account that the midpoint of strobo pulse in case of the component [A.sub.2]sin([[omega].sub.2]t+[theta]) will correspond the random phase of vibration [[omega].sub.2][delta][T.sub.2] and assuming the strobo pulse duration equivalent to the hologram exposure duration during a single period [[tau].sub.i] = [T.sub.e] the following expression is obtained:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)

It is worth mentioning that the period of strobo pulses is [T.sub.s] and in general case [T.sub.s] [not equal to] [eta][T.sub.2] and the fact is that in dependence on the frequency ratio the value 0 < [delta] < 1 can change with each strobo pulse. The interferogram obtained in this way will reconstruct the images captured by each strobo pulse at the same time. Intensity distribution of the reconstructed interferogram in this case will be described as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (10)

where j = 0,1,2 is order number of the strobo pulse. To facilitate the analysis of the obtained expression, acounting that [J.sub.-n] = [(-1).sup.n][J.sub.n] and rearranging the characteristic expression (10) results in:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)

From the expression of characteristic function a conclusion can be made that applying the proposed method of strobo pulse formation and synchronization with the necessary vibration phase, the intensity of the reconstructed image described by this function depends on duration of the strobo light pulse and on the ratio of parameters of different vibration components.

Let us analyse Eq. (12) when frequency ratio of vibration components is irrational number. Recalling the fact--if [f.sub.2] > [f.sub.1] and [f.sub.2]/ [f.sub.1] = q+[delta], where q = 1, 2, 3, ..., 0 < [delta] < 1 the strobo pulses during exposure time will be repeated with the period [T.sub.1] = [T.sub.2](q+[delta]). Evaluating all the facts above and changing summation with respect to [delta] in expression (11) by integration will give:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (12)

[FIGURE 5 OMITTED]

Thus when frequency ratio of the vibration components is irrational number the proposed methodology ensures visualization of time averaged interferogram of vibration form for one of the vibration components of the tested object. Fig. 5 presents an example of the obtained result for a tested plate natural frequencies of which are equal to 379 Hz and 2381.3 Hz.

Fig. 5 presents separation of the plate's vibration forms which correspond two vibration components of different frequencies when the ratio of vibration frequencies is irrational number: a, b - interferograms obtained by time average method representing vibration forms of the plate at [f.sub.1] = 379.0 Hz and [f.sub.2] =2381.3 Hz; c - interferogram of vibration form of the plate obtained by time average method when the plate simultaneously vibrates at the two frequencies [f.sub.1] = 379.0 Hz and [f.sub.2] =2381.3 Hz; d - interferogram of the plate's vibration form corresponding the vibration component the frequency of which is [f.sub.2] when the plate simultaneously vibrates at the two frequencies [f.sub.1] = = 379.0 Hz and [f.sub.2] = 2381.3 Hz obtained by the method of strobo holography.

Let us analyse equation (12) when frequency ratio of vibration components is rational number. It is assumed that [f.sub.2] > [f.sub.1] and [f.sub.2]/ [f.sub.1] = q, where q = 2,3, .... In this case period [T.sub.2] of the component with the frequency [f.sub.2] will be covered within the period [T.sub.1] so that [T.sub.1]/[T.sub.2] = q is a whole number and in this case strobo pulse duration for the component with the frequency [f.sub.2] per period will be - [[tau].sub.i2] = [T.sub.2]/[k.sup.1] = [T.sub.1]/q[k.sub.1]. Total time of hologram expose in this case will be n[[tau].sub.i2]. Intensity distribution of the recovered hologram in such a case is described by Eq. (8) which after evaluation of strobo pulse duration, total exposure time of the hologram and making the necessary rearrangements will take the same form as Eq. (12). It is obvious that this expression will be correct only I case when the number of strobo pulses during hologram exposure time is significantly greater that the frequency ratio q of the vibration components. Fig. 6 presents an example of the obtained result when a plate natural frequencies of which are equal to 379 Hz and 2381.3 Hz was tested.

[FIGURE 6 OMITTED]

Fig. 6 presents separation of the plate's vibration forms which correspond two vibration components of different frequencies when the ratio of vibration frequencies ([f.sub.2]/[f.sub.1] = 10) is rational number. a, b - Interferograms obtained by time average method representing the plate's vibration forms at [f.sub.1] = 764 Hz and [f.sub.2] = 7640 Hz; c - interferogram of vibration form of the plate obtained by time average method when the plate simultaneously vibrates at the two frequencies [f.sub.1] = 764 Hz and [f.sub.2] = 7640 Hz; d - interferogram of the plate's vibration form corresponding the vibration component the frequency of which is [f.sub.2] when the plate simultaneously vibrates at the two frequencies [f.sub.1] = = 764 Hz and [f.sub.2] = 7640 Hz obtained by the method of strobo holography.

4. Concluding remarks

1. The methodology of strobo light pulses formation allowing one or several light pulses to separate and synchronize with the phase of excitation signal of the tested object in one period of its vibration during the whole duration of the hologram exposure is presented.

2. For the analysis of vibrations several strobo pulses during one period of the vibration can be used and duration of the strobo light pulse can be practically changed from [T.sub.1]/2 up to [T.sub.1]/[k.sub.1].

3. The methodology of component separation of 3D vibration with two different frequency components is theoretically proven for the cases of various frequency ratios of the components.

Received October 06, 2010

Accepted January 28, 2011

References

[1.] Palevicius, A.; Ragulskis, K.; Vasiliauskas, R.; Bubulis, A. 2008. Historical aspects and future perspective of holography in Lithuania. The Society for Experimental Mechanics. Proceedings of the International Symposium to Commemorate the 60th Anniversary of the Invention of Holography. Sheraton Springfield Monach Place Hotel, Springfield, Massachusetts USA, October 27-29: 368-377.

[2.] Palevicius, A.; Vasiliauskas, R. 1998. Application of laser interferometry and holographic methods for analysis and design of vibrational systems. Virpesiu mechanika ir technika. Vibroengineering: monografija. Lietuvos Mokslu akademija. -Vilnius: Lietuvos mokslas, kn. 19: 358-381.

[3.] Vasiliauskas, R.; Palevicius, A.; Ragulskis, K. 1988. Analysis of holographic interferograms by ultrasonic piezoelectric transducers in the investigation of three-dimensional vibrations, Acoustical Physics 34(6): 573-575.

[4.] Patasiene, L.; Vasiliauskas, R.; Fedaravicius, A. 2005. Application of holographic interferometry methods for anal;ysis piezdrives with ring actuators, Ultragarsas 3(56): 17-21.

[5.] Ostasevicius, V.; Palevicius, A.; Daugela, A; Ragulskis, M.; Palevicius, R. 2004. Holographic imaging technique for characterization of MEMS switch dynamics. Proc. SPIE Vol. 5389, Smart Structures and Materials 2004: Smart Electronics, MEMS, BioMEMS, and Nanotechnology; Vijay K. Varadan; Ed.: 73-84.

[6.] Patasiene, L.; Vasiliauskas, R.; Fedaravicius, A. 2007. Holgraphic interferometry method for determination of layer piezostack parameters, Ultragarsas 1(6)2: 23-25.

[7.] Palevicius, A.; Busilas, A.; Vasiliauskas, R. 2008. Analysis holographic interferogram of links of vibrating systems, Journal of Vibroengineering (Vibromechanika), 2008 June. -Vilnius, Lithuania, volume 10, ISSUE 2, numbers of publications from 348 to 370: 245-250.

[8.] Vasiliauskas, R.; Palevicius, A. 2008. Analysis of holographic interferometry of microelectromechanical system with non harmonic cycles excitation, Journal of Vibroengineering (Vibromechanika), January/March. Vilnius, Lithuania, vol. 10 ISSUE 1, numbers of publications from 326 to 347: 42-46.

[9.] Patasien?, L.; Federavicius, A.; Vasiliauskas, R. 2008. Holographic interferometry methods for analysis and design of rotary converter, Journal of Vibroengineering (Vibromechanika), January/March. -Vilnius, Lithuania, vol. 10, ISSUE 1, numbers of publications from 326 to 347: 33-36.

[10.] Ragulskis, K.; Dabkevicius, A.; Kibirkstis, E.; Bivainis, V.; Miliunas, V.; Ragulskis, L. 2009. Investigation of vibrations of a multilayered polymeric film, Mechanika 6(80): 30-36.

[11.] Mazeika, D.; Bansevicius, R. 2009. Study of resonant vibrations shapes of the beam type piezoelectric actuator with preloaded mass, Mechanika 2(76): 33-37.

[12.] Tonin, R.; Bies, D.A. 1978. Time-average holography for the study of three-dimentional vibrations, Journal of sound and vibrations 52: 315.

[13.] Vikram, C.S. 1974. Stroboscopic holographic interferometry of vibration simultaneously in two sinusoidal modes, Optics communications, v.II, Nr.4: 360-364.

[14.] Shajenko, P.; Johnson, C.D. 1968. Stroboscopic holographic interferometry, Appl. Phys. Lett. 13(1): 4446.

[15.] Harris, J.S.; Fusek, R.L.; Marcheski, J.S. 1979. Stroboscopic interferometer, Appl. Opt. 18: 2368-2371.

[16.] Zarucki, MA.; Gurevich, V.S. 1986. On research of vibration characteristics of the parts with stress concentrators by the method of stroboholography. In textbook: Stress Concentration in the Elements of Aviation Engines. Inter university research materials. Ufa aviation Institute of S. Ordzonikidze: 97-100 (in Russian).

[17.] Vest, Ch. 1982. Holographic interferometry. Editor in chief J.S. Ostrovski. Moscow: Mir. 504p. (in Russian)

[18.] Lee, Y-S.; Kim, S-H.; Skibniewaki, M.J. 2009. Analytical and experimental approach for assessing vibration serviceability of highway bridges due to heavy vehicle traffic, The Baltic Journal of Road and Bridge Engineering 4(3): 123-133.

[19.] Bayraktar, A.; Birinci, F.; Altunisik, A.C.; Turker, T.; Sevim, B. 2009. Finite element model updating of Senyuva historical arch bridge using ambient vibration tests, The Baltic Journal of Road and Bridge Engineering 4(4): 177-185.

R. Vasiliauskas, Mykolas Romeris University, Kaunas Faculty of Public Security, V. Putvinskio 70, 44211 Kaunas, Lithuania, E-mail: [email protected]

A. Palevicius, Kaunas University of Technology, Mickeviciaus 37, 44244 Kaunas, Lithuania, E-mail: [email protected]

A. Busilas, Vilnius Gediminas Technical University, Sauletekio al. 11, 10223 Vilnius, Lithuania, E-mail: [email protected]

K. Pilkauskas, Kaunas University of Technology, Mickeviciaus 37, 44244 Kaunas, Lithuania, E-mail: [email protected]