Comparison of finite element modeling and analytical approach results for
oscillating rod structure with crack
A.V.Cherpakov1, A.N.Soloviev1'2'3, V.V.Gricenko2, S.C.Mohanty4, D.Y.Parshin2, U.I.Butenko5, O.V.Bocharova3
1 Southern Federal University, Rostov-on-Don, Russia 2 Don State Technical University, Rostov-on-Don, Russia 3 Southern Scientific Center of Russian Academy of Sciences, Rostov-on-Don, Russia 4 National institute of Technology, Rourkela, Odisha, India 5 Kazan State University of Architecture and Engineering, Kazan, Russia
Annotation: the paper considers the comparison of natural modes for the first and second transverse oscillation modes obtained by numerical simulation and finite element package of analytical calculation of the core model with an elastic element based on the application of the Euler-Bernoulli hypothesis. The analysis compared waveforms for different values of the defect shows their satisfactory agreement. Keywords: rod design, vibrations, defect, finite element modeling, analytical modeling, waveform
Problems identifying damage in rod structures in a large number of works presented in the review, published not only for a long time, but in recent years [1-12]. Construction of algorithms for solving problems of identification based on selected methods of calculation takes an important place in the study of damaged structures. In recent years, often used a relatively new approach of the development and application of evolutionary algorithms, neural networks and optimization methods [13-14].
An important factor in the identification of defects is the use of sophisticated models as an example of using the finite element approach, and an example of simple Euler-Bernoulli hypothesis [15], and require their comparison
The aim of the work is to compare the natural shapes of the first two modes of transverse vibrations of a rod model. Additional information in this method of identification of damage is a set of natural frequencies obtained from the analytical calculation for the core structure with different sizes of the defect, which is located in a particular place [16].
Model of a cantilever beam with a crack. We consider a cantilever beam of length L of rectangular cross-section (height h, width b) with a transverse notch depth t, located at a distance Lc from the termination (see Fig. 1).
Finite-element full-bodied and a simple analytical model of defective design is presented in [16]. Analytical modeling was based on a composite beam, in which
the incision is replaced by an elastic element (spring) with a constant flexural rigidity.
(a) (b)
Fig. 1. (a) Scheme of the cantilever beam to damage a notch and (b) a composite model of a
cantilever beam with an elastic element
Were calculated and compared waveforms obtained by numerical simulation and finite element package of analytical calculation of the core of the simplified model with an elastic element. Modal analysis was carried out for models with the location of the defect L c = 0.25. Analyzed a model with two values of the defect t = 0.3 and t = 0.7 . When calculated using the numerical expression (6) of [16], the flexural rigidity kt (t) of the elastic element to the analytical model were as
follows: kt (0.7) = 262 Nm / rad and k t (0.3) = 3569 Nm / rad.
Fig. 2 shows the form of the 1st (a, b), 2nd (c, d) modes of transverse vibrations resulting from FE and analytical calculations for the size of the defect (a, c) and (b, d). For comparison, the amplitude of the vibration modes were dimensionless to the amplitude of the oscillation at a point along the length of the bar, located at the free end x = 1.
At the location of the defect has a kink waveform, which is clearly evident in the third and fourth modes of oscillation at the value of the defect t = 0.3 and all selected forms of vibrational modes at a value t = 0.7 for both models. Compared waveforms obtained based on an analytical calculation waveforms obtained by FE at each point along the length of the rod. The relative magnitude of the amplitude waveforms at various points along the length of the rod is calculated as follows:
_ A - A ÄA =-—100%
Ama
(1)
>
1.0 0.8 0.6 0.4 0.2 0
1.0
Amplitude
IM
t = 0.3 .,-■■'
X = 0.25 ..-'■"
..•■' k i -.i, i, FE
' ...... .ii .i >i >. |i,i- 'I- i analit
0
0.2
Amplitude
0.4
0.6
0.8
a)
-0.6 -1.0
II M
t = 0.3
X = 0.25 /
- X
0 0.2 0.4 0.6 /0.8 1
r /
1.0 0.8 0.6 0.4 0.2 0
1.0
0.6
0.2 0 -0.2
-0.6
■ m ■ in I. '¡i I analit
-1.0
Amplitude
0.2
Amplitude
0.4
0.6
0.8
b)
II M
t = 0.7
x = 0.25 /
0 0.2 0.4 0.6 /0.8 1
■ aiia.niT. pacieT
analit
c)
d)
Fig. 2. Forms the 1st (a, b), the 2nd (c, d), transverse mode vibrations resulting FE and analytical calculations for modal quantities crack t = 0.3 (a, b) and t = 0.7 (c,
d)
Analysis comparing the corresponding curves waveforms showed the following. When comparing the amplitudes of the corresponding points along the length of the rod forms of the first oscillation mode the maximum difference is at the value of t = 0.3 Al = 0.2% of the defect and t = 0.7 Al = 1.2%. The
max max
maximum amplitude difference lies in the vicinity of the location of the defect
When comparing the amplitudes of the corresponding points along the length of the rod shape of the second mode oscillation maximum difference is at the value of t = 0.3 AAmax = 1.47% and a defect at t = 0.7 AAmax = 2.61%. The
maximum difference of the amplitudes waveforms observed in inflection waveforms ( x = 0.61). In the vicinity of the location of the defect for the two variants quantities divergence amplitudes of the second waveform does not exceed Al = 0.39%.
X
X
1
0
1
0
Comparison of the amplitudes of the corresponding points along the length of the rod shape of the third oscillation mode shows that the maximum difference of the amplitudes of the defect at the value of t = 0.3 AAmax = 2.7% and t = 0.7
AAmax = 6.5%. The maximum difference of the amplitudes correspond to the points
of the rod in the vicinity of the location of the defect for the two variants quantities.
Comparison of the amplitudes of the corresponding points along the length of the rod forms the fourth oscillation mode shows that the maximum difference is at the value of t = 0.3 AAmax = 1.49% of the defect and t = 0.7 AAmax = 4.7%. Maximum differences correspond to the amplitudes in the neighborhood of the location of the defect.
Conclusions. The analysis compared waveforms for different values of the defect shows that the qualitative characteristics of the curves forms of oscillations in the vicinity of the location of the defect, and along the length of the rod, the same for both models.
The present study was carried out with the partial support of RFBR (grants No 14-38-50933 mol_nr, No 13-08-00516-A)
References
1. Dimarogonas A.D. 1996. Vibration of cracked structures: a state of the art review. Eng. Fract.Mech. 55, pp.831-857. (doi:10.1016/0013-7944(94)00175-8).
2. Friswell M.I. Damage identification using inverse methods // Phil. Trans. R. Soc. A (2007) 365, pp. 393-410. (doi:10.1098/rsta.2006.1930.).
3. Bhat, R. B., and Chakraverty, S. (2007).Numerical analysis in engineering, 2nd Ed., Narosa Publishing House, New Delhi, India. 320 p.
4. Kosenko, E.E., Kosenko, V.V., Cherpakov, A.V. Modelling of rods with defects with different types of fastening // Inzenernyj vestnik Dona (Rus), 2013, №4 URL:ivdon.ru/magazine/archive/n4y2013/2155.
5. Kosenko, E.E., Kosenko, V.V., Cherpakov, A.V. The studying of oscillations of the corpulent cantilever rod model with a defect // Inzenernyj vestnik Dona (Rus), 2013, №4 URL:ivdon.ru/magazine/archive/n4y2013/2153.
6. Denina O.V., Vatulyan A.O. inverse problems for rods. Methods for determination of properties of inhomogeneous elastic rods based acoustic sensing. Saarbrücken, Germany: Lap Lambert Academic Publishing 2011 - 104 p.
7. Burtseva O.A., Chipco S.A., Kaznacheeva O.K., Cherpakov A.V. Vibration control for high-rise constructions. European Journal of Natural History. 2012. № 4. pp. 39-44
8. Akopyan, V. Parameter Estimation of Pre-Destruction State of the Steel Frame Construction Using Vibrodiagnostic Methods [Text]/ V. Akopyan, A. Soloviev, A. Cherpakov // Mechanical Vibrations: Types, Testing and Analysis. - N-Y.: Nova Science Publishers, 2010. - Chapter 4. - pp.147-161.
9. Fan, W., Qiao, P. Vibration-based damage identification methods: A review and comparative study."Struct. Health Monit., 2011. 10(1), pp. 83-111.
10. Mohanty, S.C., I. Ramu. A Review on Free, Forced Vibration Analysis and Dynamic Stability of Ordinary and Functionally Grade Material Plates, Caspian Journal of Applied Sciences Research, 1(13), pp. 57-70, 2012 .
11. Mohanty, S.C., Dash, R.R., Rout, T.,"Static and Dynamic Stability Analysis of Functionally Graded Timoshenko Beam on Winkler Foundation" Journal of Engineering Research and Studies, vol-1, issue-II, pp. 149-165, 2010.
12. Il'gamov M.A., Khakimov A.G. Diagnosis of damage of a cantilever beam with a notch // Russian Journal of Nondestructive Testing. 2009. Vol. 45. №6. pp. 430435.
13. Rao M. Anand, Srinivas, J. and Murthy, B. S. N., Damage detection in vibrating bodies using genetic algorithms,Comput. & Structures, Vol. 82, pp. 963968, 2004.
14. A. A. Krasnoshchekov, B. V. Sobol, A. N. Solov'ev, and A. V. Cherpakov // Identification of Crack Like Defects in Elastic Structural Elements on the Basis of Evolution Algorithms. ISSN 1061_8309, Russian Journal of Nondestructive Testing, 2011, Vol. 47, No. 6, pp. 412-419.
15. Panigrahi, S. K., Chakraverty, S., and Mishra, B. K. Vibrationbased damage detection in a uniform strength beam using genetic algorithm. 2009. Meccanica, 44(6), pp.697-710.
16. Analytical and finite element analysis of vibration parameters in terminal damage/ V. A., Akopian, A. V. Cherpakov, A. N. Soloviev and others, Izvestiya vuzov North Caucasus region. Technical science. 2010. No.5. - pp.21-28.