**American Journal of Mechanical Engineering**

## Analysis of Reliability of Modal Parameters Estimation Using High-speed Digital Image Correlation Method

**Martin Hagara**^{1,}, **Róbert Huňady**^{1}, **Pavol Lengvarský**^{1}, **Peter Pavelka**^{1}

^{1}Department of Applied Mechanics and Mechanical Engineering, Technical University of Košice, Faculty of Mechanical Engineering, Košice, Slovakia

Abstract | |

1. | Introduction |

2. | Modan3D |

3. | Experimental Modal Analysis of Rectangular Plate Using Modan3D |

4. | Conclusions |

Acknowledgement | |

References |

### Abstract

This paper deals with description of modification and using high-speed digital image correlation method in experimental modal analysis. In its theoretical part the basic principle of digital image correlation and description of modification of correlation system Q-450 Dantec Dynamics for the purposes of SIMO (Single Input / Multiple Output) assumption of modal parameters is introduced. The practical part – analysis of reliability of program Modan3D, developed for modal parameters estimation – was realized on a simple steel specimen in laboratory conditions. The obtained outputs are present in a form of graphs and tables.

**Keywords:** digital image correlation, experimental modal analysis, reliability, Modan3D

**Copyright**© 2015 Science and Education Publishing. All Rights Reserved.

### Cite this article:

- Martin Hagara, Róbert Huňady, Pavol Lengvarský, Peter Pavelka. Analysis of Reliability of Modal Parameters Estimation Using High-speed Digital Image Correlation Method.
*American Journal of Mechanical Engineering*. Vol. 3, No. 6, 2015, pp 190-194. http://pubs.sciepub.com/ajme/3/6/7

- Hagara, Martin, et al. "Analysis of Reliability of Modal Parameters Estimation Using High-speed Digital Image Correlation Method."
*American Journal of Mechanical Engineering*3.6 (2015): 190-194.

- Hagara, M. , Huňady, R. , Lengvarský, P. , & Pavelka, P. (2015). Analysis of Reliability of Modal Parameters Estimation Using High-speed Digital Image Correlation Method.
*American Journal of Mechanical Engineering*,*3*(6), 190-194.

- Hagara, Martin, Róbert Huňady, Pavol Lengvarský, and Peter Pavelka. "Analysis of Reliability of Modal Parameters Estimation Using High-speed Digital Image Correlation Method."
*American Journal of Mechanical Engineering*3, no. 6 (2015): 190-194.

Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |

### At a glance: Figures

### 1. Introduction

Currently the work of engineers and constructors is focused on designing of reliable structures with comfortable noiseless operation. For a reason that newer materials are increasingly lighter it is necessary to ensure their safety. These requirements suggest that it is even more necessary to pay attention to the dynamic characteristic of the structures.

Dynamic characterization of structures is of significant importance in a wide variety of industries including aerospace, traffic ^{[1]}, petroleum industry ^{[2]}, civil structures, appliances, however also in sport ^{[3]}. In general, for the investigation of the dynamic behavior of structures ^{[4]} a method serving for estimation of modal parameters can be used. Such method, by which natural frequencies, mode shapes and modal damping are estimated, is called modal analysis ^{[5]}. There are two forms of modal analysis – experimental modal analysis (EMA) and operational modal analysis (OMA). While the first one uses for excitation of the structures an additional hardware as impact hammer or shakers, the second one investigates the structures in situ, i.e. during operation.

The responses of the structures are commonly captured using traditional mechanical transducers of accelaration. Several past years ago the non-contact optical devices, based e.g. on laser-Doppler principle, have been used for measuring of structure surface velocity. To the group of optical techniques belongs also a method of digital image correlation (DIC).

DIC is a modern non-contact method serving for 2D or 3D reconstruction of investigated object displacement and deformation fields. The possibility of 3D measurement is dependent on the use of minimally two cameras with CCD (Charge Coupled Devices) or CMOS (Complementary metal–oxide–semiconductor) sensor for sampling of the object during its loading.

**Fig**

**ure**

**1**. Low-speed CCD cameras of correlation system Q-400 Dantec Dynamics

**Fig**

**ure**

**2**. High-speed CMOS cameras of correlation system Q-450 Dantec Dynamics

According to maximal sampling frequency of camera sensor digital image correlation systems can be divided into two groups – low-speed and high-speed systems. While the low-speed cameras (Figure 1) are able to capture maximally several frames per second, the high-speed ones (Figure 2) can take several thousand frames per second.

During correlation digital images are not correlated (compared) as whole, but among small picture elements called facets (Figure 3). Commonly delivered correlation systems allow determine displacement fields in directions x, y and z as well as strain fields , and for each nodal point of virtual grid created on investigated object surface. The spacing between grid nodal points should be equal to 3/4 of facet size.

**Fig**

**ure**

**3**. Investigated object surface divided into facets and virtual grid

There are many possibilities, how to utilize digital image correlation technique. The low-speed correlation systems are commonly used in stress/strain analysis ^{[6]}, or by prediction of locations of strain concentrators ^{[7]}. The high-speed systems serve as tools for motion analysis ^{[8]}, vibration analysis ^{[9]}, crash tests as well as drop tests and can be used also in failure analysis.

For a reason that high-speed digital image correlation systems allow measurements at relatively high sampling frequencies and the results are obtained in every nodal point, it is very convenient to use correlation systems for the purpose of experimental modal analysis.

### 2. Modan3D

Modan3D is a tool, created in Matlab, allowing import and processing of data obtained by Istra4D, what is a software delivered with digital image correlation systems Dantec Dynamics. Istra4D allows export data in a form of AVI, STL as well as HDF5 AVI format. While exporting in AVI format is used for creating of videos, data exported in STL format serves as input for CAx software, HDF5 data format (Figure 4) with hierarchical structure allows relatively easy processing.

As the frequency response function (FRF), which describes the dynamic behavior of structure in frequency spectrum, is given by relation between output and input signal, it was necessary to modify high-speed correlation system Q-450 Dantec Dynamics for acquisition of force signal from impact hammer. For that reason two additional devices were added to measuring string – CCLD (Constant Current Line Drive) amplifier and AC/DC converter NI USB-4431 with four input and one output channel. The block scheme of connection between high-speed correlation system and additional devices can be seen in Figure 5.

**Fig**

**ure**

**4**. Hierarchical structure of HDF5 file exported from Istra4D

**Fig**

**ure**

**5**. Block scheme of modified correlation system Q-450 Dantec Dynamics for the purposes of experimental modal analysis

Program Modan3D allows read the data exported from Istra4D in a form of displacements in particular directions x, y and z as well as the time dependence of force signal obtained from impact hammer. Mentioned data are then transformed into frequency spectrum using Fast Fourier Transform (FFT).

Frequency response matrices in each point of object surface are obtained using following formula:

(1) |

where , ** **and are matrices of object displacements in frequency spectrum and is a force input in frequency spectrum.

Currently Modan3D utilizes for assumption of natural frequencies two functions – Normal Mode Indicator Function (NMIF) and Complex Mode Indicator Function (CMIF) ^{[10]}. The values of NMIF for each spectral line of the frequency spectrum are obtained using:

(2) |

where summation is performed on measured frequency response functions.

The second function, CMIF, is obtained using singular value decomposition from:

(3) |

where symbolizes the left singular matrix corresponding to modal vectors matrix, denotes the right singular matrix corresponding to modal vectors participation matrix and represents a diagonal singular values matrix. Particular singular values are then expressed as:

(4) |

Ideal courses of NMIF and CMIF functions can be seen in Figure 6.

Modan3D can calculate damping ratios, whereby it utilizes the half-power method ^{[11]}. More information about program Modan3D, created at authors department, can be found in publication ^{[12]}.

**Fig**

**ure**

**6**. Ideal courses of CMIF and NMIF functions

### 3. Experimental Modal Analysis of Rectangular Plate Using Modan3D

The estimation of modal parameters in a form of natural frequencies, modal shapes and damping ratios was performed on a steel sheet of thickness 0.8 mm, which shape and dimensions are depicted in Figure 7. The boundary conditions for the specimen were fixation along its narrower edge with free other three edges.

**Fig**

**ure**

**7**. Dimensions and boundary conditions of the specimen

For the excitation an impact hammer Brüel & Kjær 8206. The impact hammer is supplied with three interchangeable impact tips of aluminum, plastic and rubber. The choice of impact tip determines the impulse shape (amplitude and duration) and the bandwidth of the excitation ^{[13]}. The sampling frequency of high-speed cameras was set to 2000 fps, what means that a frequency spectrum from 0 Hz to 1000 Hz was investigated. According to courses depicted in Figure 8 and Figure 9, characterizing the impulse shape and force spectrum of the impact, respectively, the impact hammer with plastic tip for the excitation of the specimen was chosen.

**Fig**

**ure**

**8**. Time responses of impulse shapes for three available hammer tips [13]

**Fig**

**ure**

**9**. Force spectrum of an impact on an aluminum plate

For the analysis of reliability of Modan3D it was necessary to perform the same measurement many times at various conditions such as alternating magnitude of force impacts, locations of impacts etc. The best solution should be that several different people perform the same measurement. As there was a problem to ensure sufficient amount of qualified people, we decided to repeat the measurement ten times (Figure 10).

**Fig**

**ure**

**10**. Excitation of the specimen using impact hammer

The amounts of impact forces measured during experiment varied from minimal value of 12.4 N up to maximal value of 19.8 N are mentioned in Table 1.

The results of each measurement obtained in a form of CMIF and NMIF functions were compared together. As can be seen from Figure 11 and Figure 12, the courses of mentioned courses correspond well together.

**Fig**

**ure**

**11**. CMIF functions obtained from ten measurements by Modan3D

**Fig**

**ure**

**12**. NMIF functions obtained from ten measurements by Modan3D

The reliability of the modal parameters estimation using program Modan3D was tested by observation of standard deviations, dissipations and average relative deviations of natural frequencies as well as damping ratios amounts of particular modes obtained from mentioned measurements. The obtained data were processed and are present in Table 2 and Table 3.

The abbreviations M1-M10, mentioned in Table 2. and Table 3., denote particular measurements, AA – average amount of the obtained results, SD – standard deviations, D – dissipations and ARD – average relative deviations.

Mode shapes corresponding to mentioned natural frequencies can be found in Figure 13.

**Fig**

**ure**

**13.**Mode shapes of the investigated specimen obtained by Modan3D

### 4. Conclusions

Digital image correlation is a modern technique allowing a wide range of applications. In this article its utilization for the purposes of SIMO experimental modal analysis is described. With regard to classical EMA methodology such performing of modal parameters estimation is very convenient. The main reasons are: it is not necessary to create model of investigated structure, characterize the degrees of freedom and excite the structure (or measure the response of the structure) in all degrees of freedom, what make the measurement process relatively easy and not time-consuming. The data processed from ten measurements, from which average relative deviations lower than 0.1% for natural frequencies estimation and 10% for very small amounts of damping ratios were calculated, approve that the modified high-speed correlation system Q-450 Dantec Dynamics with program Modan3D gives reliable results. The obtained mode shapes of the specimen correspond to the mode shapes computed using finite element method.

### Acknowledgement

This contribution is the result of the projects implementation VEGA 1/0937/12 “Development of nontraditional experimental methods for mechanical and mechatronical systems” and APVV-0091-11 “The use of experimental and numerical methods for the increase of competitiveness and innovation of mechanical and mechatronical systems”.

### References

[1] | Klimenda, F. et al., “Investigation of Vertical Vibration of a Vehicle Model Driving Through a Horizontal Curve”, Manufacturing Technology, Vol. 15, no. 2, 143-148, 2015. | ||

In article | |||

[2] | Zheng, L. et al., “Vibration cause analysis and elimination of reciprocating compressor inlet pipelines”, Engineering Failure Analysis, Vol. 48, no. 2, 272-282, 2015. | ||

In article | |||

[3] | Svoboda, M. et al., “Effect of Impacts on Human Head”, Manufacturing Technology, Vol. 15, no. 2, 226-231, 2015. | ||

In article | |||

[4] | Svoboda, M., Soukup, J., “Dynamic measurement of four-axle railway wagon”, Manufacturing Technology, Vol. 13, no. 4, 2013. | ||

In article | |||

[5] | Nangolo, F.N., Klimenda, F., “Modal Parameter Analysis for Underdamped Mechanical Systems”, Applied Mechanics and Materials, Vol. 732, 247-252, 2015. | ||

In article | View Article | ||

[6] | Zhu, F. et al., “Measurement of true stress-strain curves and evolution of plastic zone of low carbon steel under uniaxial tension using digital image correlation, Optics and Lasers in Engineering, Vol. 65, 81-88, 2014. | ||

In article | View Article | ||

[7] | Fernández, J.Á.P. et al., Measuring Strain Concentrations in Welded Junctions using Digital Image Correlation, in: YPIC 2014 : Young welding Professionals International Conference, Budapest, Hungary, 17-21, 2014. | ||

In article | |||

[8] | Hagara, M. et al, “A determination of the kinematic quantities of a rotating object by digital image correlation method”, American Journal of Mechanical Engineering, Vol. 1, no. 7, 289-292, 2013. | ||

In article | |||

[9] | Huňady, R. et al., The Application of High-speed Digital Image Correlation in Vibration Analysis of a Rotating Object, in: Proceedings of EAN 2012 : 50th annual conference on experimental stress analysis, Tábor, Czech Republic, 1-8, 2012. | ||

In article | |||

[10] | Huňady, R. and Hagara, M., “Experimental Investigation of Mode Shapes of Symmetric Structures”, Acta Mechanica Slovaca, Vol 19, No. 3, 2015, 12-17 p. | ||

In article | |||

[11] | Huňady, R. et al., “Using High-speed Digital Image Correlation to Determine the Damping Ratio”, Procedia Engineering, Vol. 48, 242-249, 2012. | ||

In article | View Article | ||

[12] | Trebuňa, F., Hagara, M., “Experimental modal analysis performed by high-speed digital image correlation system”, Measurement: Journal of the International Measurement Confederation, Vol. 50, no. 1, 78-85, 2014. | ||

In article | |||

[13] | Brüel & Kjær Product Data Sheet: Impact Hammers – types 8206, 8206-001, 8206-002 and 8206-003. | ||

In article | |||