^{1, 2}Institute of Mechanics and Machine Design Fundamentals, Czestochowa, Poland
^{1}Corresponding author
Journal of Vibroengineering, Vol. 20, Issue 1, 2018, p. 215224.
https://doi.org/10.21595/jve.2017.18733
Received 7 June 2017; received in revised form 30 September 2017; accepted 8 November 2017; published 15 February 2018
JVE Conferences
This paper discusses the formulation and solution for the problem of damped transverse vibrations of the $\mathrm{\Gamma}$ type frame with open cracks. Dissipation of vibration energy in the frame results from the movement constraint in the column and bolt support (constructional damping) and internal damping of viscoelastic material of the frame (rheological model by KelvinVoigt). Presence of a crack impacts local flexibility that has an effect on frame vibration response. The boundary problem for the above system was formulated based on the Hamilton’s principle and solved numerically for the complex eigenvalues ${\omega}^{*}$. The effect of the crack depth and its location on damped vibration was presented. The effect of both types of damping on the degree of amplitude decay was also presented.
Keywords: damped vibration, open crack, amplitude decay factor, Γ frame.
1. Introduction
Cracks are one of the most common causes of structural defects. This kind of damage to structural components may significantly affect machine performance. The dynamic behaviour of a system with a crack on a structural component is changed. Crack in a structure has also a major effect on object vibration. Taking into consideration the constructional and internal damping (and its effect on vibration of the structure) leads to important and interested conclusions. This study was concerned with an effect of the location and depth of cracks on damped vibration of the $\mathrm{\Gamma}$ type frame for selected damping coefficients.
The theoretical, numerical and experimental analysis of vibrations and stability of two types of flat frames i.e. $T$shaped and $\mathrm{\Gamma}$ shaped were presented in studies [1, 2], respectively. The frames were loaded with longitudinal force in relation to the bolt or loaded by a force oriented toward the positive pole. No cracks in the structure were found for any of the frames.
Presence of a crack in an elastic structural component leads to considerable local flexibility due to the strain energy concentration near the crack tip under load. The effect of a transverse open crack on free vibration of a simply supported BernoulliEuler prismatic beam was presented in the study [3]. In order to model the cracks and their continuous flexibility, authors used the displacement field near the crack, found based on fracture mechanics methods. They presented results of two independent cases of evaluation of the problem of crack identification in structures: the continuous cracked beam vibration theory and the lumped cracked beam vibration analysis. The study [4] presented the method to identify cracks based on the assumption of a transverse surface crack. The proposed method is noninvasive and can be easily implemented and used to monitor crack growth. A local flexibility matrix describes a force – displacement relationship. The transverse open singlesided crack is modelled by a rotational spring with calculated stiffness constant. The authors of the study [5] presented an analysis of the effect of two open cracks on frequencies of natural flexural vibrations in a cantilever beam. In this case, the two types of cracks are to be considered: open doublesided and open singlesided. Stability and natural vibration of columns composed of one or several elements with crack in the internal rod were presented in [68]. The study [6] described the effect of cracks in a structure on instability and free vibrations of stepped columns. In papers [7, 8], the vibration frequency, bifurcation load magnitude, bifurcation loadcrack size relationship, and amplitudevibration frequency relationship were studied. The crack was modelled using a rotational spring with a linear characteristic. The results of the numerical calculations presented in [68] indicated that the crack size has a significant effect on bifurcation load magnitude and vibration frequency.
The effect of a local flexibility (crack) on the vibrations of different type of frames was presented in publications [911]. The study [9] discussed the impact of crack localization on transverse vibrations of $\mathrm{\Gamma}$ type frame with and without additional discrete elements. Transverse and longitudinal vibrations of two and threebar frames with cracks were presented in the study [10]. In [9] and [10], local flexibility (crack) was modelled by a local stiffness matrix. A modal analysis of twomember open frame with a transverse crack was the focus of the study [11]. The authors analysed free vibrations of a 3D model of open $L$shape frame for various crack dimensions, locations and values of tip mass for both horizontal and vertical members. The findings of the research presented in [311] showed that natural frequency of the system is lower for location of the crack near fixed end region and increases for farther regions.
In these studies authors focused on the description of stability, free vibration or the effect of cracks on frame vibration, whereas the effect of damping was not taken into consideration. In the study [12], authors formulated and solved the problem of transverse damped vibrations of T type frame with rotational dampers in the points of frame mounting and in the supports.
The focus of the present study is on the analysis of damped vibrations of the $\mathrm{\Gamma}$ type frame with cracks. Previous studies in this field have examined integrated transverse and longitudinal damped vibrations of the $\mathrm{\Gamma}$ type frame without cracks in a structure [13]. In the BernoulliEuler beam model adopted in the study, dissipation of vibration energy derives from the damping. Fig. 1 presents two types of damping: constructional damping in the column support and in fixing bolt frame support and internal damping of the viscoelastic material of the frame. The rotational viscous dampers ${C}_{R}$ were used to model the constructional damping of supports. Internal damping in individual frame beams was denoted by ${E}_{i}^{*}$. The complex eigenvalue of the frame ${\omega}^{*}=Re\left({\omega}^{*}\right)\mp jIm\left({\omega}^{*}\right)$ was given. The real part of the solution described the damped vibration frequency. The imaginary part of the solution described the degree of amplitude decay. The numerical analysis was focused on an effect of the location and depth of cracks on damped vibration frequency and on degree of amplitude decay of the $\mathrm{\Gamma}$ type frame. The results obtained were presented by means of three, and twodimensional graphs. Solutions for additional control of system dynamics and support for early detection of damages are being explored. The problem presented in this work necessitates undertaking actions oriented at the identification. Such actions are: mathematical modelling (included in this paper), experimental research, parameter estimation, model verification and optimization. The all or some of identification stages (actions) are the subject of research for many authors, ex. [1517], which served as an inspiration to undertake this study.
2. Physical and mathematical models of the system
The diagram of the system discussed in the present study is shown in Fig. 1. In this system, cracks are modelled by rotational spring ${K}_{i}$. Rotational cracks compliance was assumed to be dominant in the local flexibility matrix as in [4] in the following form:
where:
and:
where: ${c}_{i}$_{}– compliance, ${E}_{i}$ – Young’s modulus, ${J}_{i}$ – moments of inertia for crosssectional areas of beams, ${I}_{i}$ – dimensionless local compliance function [4], ${a}_{i}$ – crack depth, $h$ – height of the column and bolt frame, $i=$ 1, 2.
The physical and mathematical model of the system takes into account the viscoelasticity of the material characterized by Young’s modulus ${E}_{i}$, viscosity coefficients ${E}_{i}^{*}$ (internal damping) and movement resistance in the supports modelled by rotational viscous dampers (constructional damping).
Fig. 1. The diagram of the $\mathrm{\Gamma}$ type frame with viscous end conditions and cracks in the structure
Equations of motion for individual beams were denoted as:
where: ${W}_{i}\left(x,t\right)$ – transverse displacements, ${E}_{i}$ – Young’s modulus, ${E}_{i}^{*}$ – viscosity coefficients, ${A}_{i}$ – crosssectional areas of beams, ${J}_{i}$ – moments of inertia for crosssectional areas, ${\rho}_{i}$ – beam material densities, $t$ – time, $x$ – spatial coordinate, $i=$ 1, 2, 3, 4.
Solutions of Eq. (4) are in the form:
where ${\omega}^{*}$ – complex eigenvalue of the system, ${\omega}^{*}=Re\left({\omega}^{*}\right)\mp jIm\left({\omega}^{*}\right)$, $j=\sqrt{1}$.
Substitution of Eq. (5) into Eq. (4) leads to:
where:
Geometrical boundary conditions and continuity conditions are given by:
${w}_{3}\left({l}_{3}\right)={w}_{4}\left(0\right),{w}_{4}\left({l}_{4}\right)=0.$
Natural boundary conditions of the system can be presented in the following form:
${{E}_{1}{J}_{1}w}_{1}^{III}\left({l}_{1}\right)={{E}_{2}{J}_{2}w}_{2}^{III}\left(0\right),{{E}_{2}{J}_{2}w}_{2}^{II}\left({l}_{2}\right)={{E}_{3}{J}_{3}w}_{3}^{II}\left(0\right),{{E}_{3}{J}_{3}w}_{3}^{II}\left({l}_{3}\right)={{E}_{4}{J}_{4}w}_{4}^{II}\left(0\right),$
${w}_{3}^{I}\left({l}_{3}\right)+\frac{{E}_{4}{J}_{4}}{{K}_{2}}{w}_{4}^{II}\left(0\right)={w}_{4}^{I}\left(0\right),{{E}_{3}{J}_{3}w}_{3}^{III}\left({l}_{3}\right)={{E}_{4}{J}_{4}w}_{4}^{III}\left(0\right),$
${{E}_{4}{J}_{4}w}_{4}^{II}\left({l}_{4}\right)=j{C}_{R}{\omega}^{*}{w}_{4}^{I}\left({l}_{4}\right).$
The solution for the Eq. (5) is given by:
where:
Substitution of Eq. (9) to Eq. (8ab) yields a homogeneous system of equations with respect to unknown constants ${C}_{ki}$, which, in the matrix form, can be written as:
where:
The system has a nontrivial solution if the determinant of the matrix coefficient (with constant ${C}_{ki}$) equals to zero:
Solution of the above equation (Eq. (13)) yields the real and imaginary parts of the vibration frequency.
3. Results of numerical calculations
Calculations were performed for the variable crack locations in the column (first crack) and in the bolt frame (second crack), selected depths of cracks and selected values of damping coefficients. The data are contained in Table 1.
3.1. Parameters of the system
The nondimensional parameters were used for: internal damping of viscoelastic material of beams – $\eta $, constructional damping in the supports – $\mu $, location of crack in the column and in the bolt – ${L}_{C1}$ and ${L}_{C2}$, respectively:
where:
Table 1. Parameters of the model
Parameter

Value

Total length of the column ${L}_{1}$ [m]

2

Total length of the bolt ${L}_{2}$ [m]

2

Length of the beam ${l}_{1}$ [m]

Variable

Length of the beam ${l}_{2}$ [m]

Variable

Length of the beam ${l}_{3}$ [m]

Variable

Length of the beam ${l}_{4}$ [m]

Variable

Height of the column and bolt $h$ [m]

0.1

Width of the column and bolt $b$ [m]

0.1

Crack depth in the column ${a}_{1}/h$

Variable

Crack depth in the column ${a}_{2}/h$

Variable

Location of crack in the column ${L}_{C1}$

Variable

Location of crack in the column ${L}_{C2}$

Variable

Young modulus of the column and bolt ${E}_{i}$ [Pa]

2.1e11

Internal damping coefficient $\eta $

0.001

Material density in column and bolt ${\rho}_{i}$ [kg/m^{3}]

7860

Value of $\eta $ used in computations is internal damping coefficient of homogeneous structural steel without coating [14]. Selected values of $\mu $ from 0 to 2 describe the main effect of constructional damping coefficient on damped vibration frequency and on the degree of vibration amplitude decay. No effect on $Re\left({\omega}^{*}\right)$ and $Im\left({\omega}^{*}\right)$ for $\mu =$ 0 was found whereas the highest effect on $Im\left({\omega}^{*}\right)$ and $Re\left({\omega}^{*}\right)$ was observed for $\mu =$ 0.4 and $\mu >$ 2, respectively.
3.2. Presentation of research results
The results of the calculations are presented in Figs. 210. The relationships between the first ${\omega}_{1}$ (a) and the second ${\omega}_{2}$ (b) natural frequency of the frame and the depth and location of crack in the column without damping are presented in Fig. 2. The relationships represent reference results. The comparison of the results with those reported in literature demonstrated that the increase in the value of crack depth leads to the decrease in eigenvalues of the undamped system.
Fig. 2. Relationship between the first and the second natural frequencies of the frame and location of: a) the first crack without the second crack, b) without damping
a)
b)
The aim of this paper is to analyse the effect of damping on vibration of a $\mathrm{\Gamma}$ type frame with cracks. In the first case, there was a single crack on the column, whereas in the second case, there were two cracks: one on the column and one on the bolt frame. The relationships of the real $Re\left({\omega}^{*}\right)$ and the imaginary $Im\left({\omega}^{*}\right)$ terms of the first and the second eigenvalue of the frame for different locations of first crack ${L}_{C1}$ and depth ${a}_{1}/h$ are presented in Figs. 3(a), 3(b), 4(a) and 4(b). Calculations were performed for damping coefficients $\eta =$ 0.001 and $\mu =$ 0.4. Figs. 34 show some areas of the curves, with the increase in the value of the crack depth leading to the increase in eigenvalues of the damped system.
Fig. 3. Relationship between the first damped eigenvalue of the frame and location of the first crack without the second crack
a)
b)
Fig. 4. Relationship between the second damped eigenvalue of the frame and location of the first crack without the second crack
a)
b)
Fig. 5. The relationship between the first damped eigenvalue of the frame and location of the first crack for the depth of the second crack ${a}_{2}/h=$ 0.5 and its location ${L}_{C2}=$ 0.5
a)
b)
Fig. 5(a) and 5(b) present the relationships between the real and imaginary parts of the first eigenvalue for the frame and the variable value of depth ${a}_{1}/h$ and location ${L}_{C1}$_{}of the crack in the column for internal damping coefficient $\eta =$ 0.001, constructional damping coefficient $\mu =$ 0.4 and for the depth ${a}_{2}/h=$ 0.5 and location ${L}_{C2}=$ 0.5 of the crack in the frame bolt. Figs. 6(a) and 6(b) present the same relationships for the second eigenvalue for the frame. The comparison of the results shown in Figs. 3, 4 with the results presented in Figs. 5, 6 reveals that, in this second case, changes in the location of the crack on the column lead to the increase in both real and imaginary parts of the first and the second eigenvalue of the system. Some areas of the curves are also nonmonotonic as in Figs. 34.
Fig. 6. The relationship between the second damped eigenvalue of the frame and location of the first crack for the depth of the second crack ${a}_{2}/h=$ 0.5 and its location ${L}_{C2}=$ 0.5
a)
b)
Further part of the study was focused on the determination of the relationship of the first and the second eigenvalue of frame on changes in constructional damping coefficient $\mu $ and the depth of both cracks ${a}_{2}/h={a}_{1}/h$ at ${L}_{C1}={L}_{C2}=$ 0.5 and internal damping coefficient $\eta =$ 0.001. Results of investigations are presented in Figs. 7(a), 7(b), 8(a) and 8(b). A simultaneous change in the depth of both cracks causes insignificant changes in $Re\left({\omega}_{1}*\right)$, $Im\left({\omega}_{1}*\right)$, $Re\left({\omega}_{2}*\right)$ and $Im\left({\omega}_{2}*\right)$ of the system. Considerable changes can be observed in $Im\left(\omega *\right)$ of the first and the second eigenvalue with changes in the value of constructional damping $\mu $.
Fig. 7. The dependence of: a) real $Re\left(\omega *\right)$, b) imaginary $Im\left(\omega *\right)$ parts of the first eigenvalue for the frame on the constructional damping $\mu $ and crack depth ${a}_{1}/h$ for ${a}_{2}/h={a}_{1}/h$ and ${L}_{C1}={L}_{C2}=$ 0.5
a)
b)
Figs. 9(a) and 9(b) present the dependence of the real and the imaginary part of the first eigenvalue of the frame on changes in constructional damping coefficient $\mu $ and changes in cracks location of both cracks ${L}_{C1}={L}_{C2}$ at ${a}_{2}/h={a}_{1}/h=$ 0.5 and for internal damping coefficient $\eta =$ 0.001. The effect of simultaneous changes in crack location ${L}_{C1}={L}_{C2}$ and in constructional damping coefficient $\mu $ on the second eigenvalue of the frame was presented in Figs. 10(a) and 10(b). Substantial changes can be observed in the first and the second eigenvalue of the system for different crack locations. The constructional damping causes similar changes in $Im\left(\omega *\right)$ presented in Figs. 7(b), 8(b) to those in $Im\left(\omega *\right)$ in Figs. 9(b), 10(b).
Fig. 8. The dependence of: a) real $Re\left(\omega *\right)$, b) imaginary $Im\left(\omega *\right)$ parts of the second eigenvalue for the frame on the constructional damping $\mu $ and crack depth ${a}_{1}/h$ for ${a}_{2}/h={a}_{1}/h$ and ${L}_{C1}={L}_{C2}=$ 0.5
a)
b)
Fig. 9. The dependence of: a) real $Re\left(\omega *\right)$, b) imaginary $Im\left(\omega *\right)$ part of the first eigenvalue for the frame on the constructional damping $\mu $ and crack location ${L}_{C1}$ for ${a}_{2}/h={a}_{1}/h=$ 0.5 and ${L}_{C1}={L}_{C2}$
a)
b)
Fig. 10. The dependence of: a) real $Re\left(\omega *\right)$, b) imaginary $Im\left(\omega *\right)$ part of the second eigenvalue for the frame on the constructional damping $\mu $ and crack location ${L}_{C1}$ for ${a}_{2}/h={a}_{1}/h=$ 0.5 and ${L}_{C1}={L}_{C2}$
a)
b)
Important changes can be observed in $Im\left({\omega}_{1}^{*}\right)$ and $Im\left({\omega}_{2}^{*}\right)$ in the case of changes in the constructional damping coefficients (Figs. 7(b), 8(b), 9(b), 10(b)).
Limitation of the scope of the study to the analysis of the two first eigenvalues of the system results from their fundamental importance to engineering practice.
4. Conclusions
This paper analysed two cases of crack location. The effect of damping on vibration of $\mathrm{\Gamma}$ type frame was examined. The findings of the study lead to the following conclusions:
– Taking into account the effect of constructional damping on vibration of the frame causes quantitative and qualitative functional changes, as shown in Figs. 24. These significant changes, between the real part of the first two eigenvalues (Figs. 3(a), 4(a)) and the first two natural frequencies (Fig. 2) resulted from movement constraint of the supports of the frame.
– Location of the crack in the system with the rotational damper in the support has a significant effect on $Re\left({\omega}^{*}\right)$ in both the first and the second eigenvalue of the system. If the crack is located near the support, up to ${L}_{C1}=$ 0.3, the increase in the value of $Re\left({\omega}^{*}\right)$ and its character (Figs. 3(a), 4(a)) is observed compared to the system where dumping was not taken into consideration (Fig. 2.). In the range of ${L}_{C1}$ discussed in the study, the crack depth ${a}_{1}/h$ also has a significant effect on the values of $Re\left({\omega}_{1}^{*}\right)$ and $Re\left({\omega}_{2}^{*}\right)$.
– The method adopted in the study can be successfully used in noninvasive monitoring of the object status. The numerical analysis should assume the presence of natural damping (internal damping and tribological wear or other movement resistance such as structural damping) which occurs in the system studied in order to obtain reference results for the systems with and without cracks. The results of the actual examination of the object, obtained by means of modal analysis should be compared with reference results.
– The adopted methodology needs to be verified experimentally, which will be the aim of future research.
– An increase in $\mu $ results in the increase in $Im\left({\omega}^{*}\right)$ to maximum values, followed by $Im\left({\omega}^{*}\right)\to $ 0 where $\mu \to \infty $ (Figs. 710).
– An increase in $Re\left({\omega}^{*}\right)$(Figs. 710) is caused by the intervention in the conditions of the frame supports, and in extreme cases where $\mu \to \infty $, the supports are changed from joint into rigid mountings.
– Analysis of the results presented in Figs. 710 allows for determination of certain geometrical parameters, such as depth and location of crack and constructional damping parameters of the frame, at which the degree of amplitude decay is the highest i.e. the vibration amplitude decreases the fastest.
– The rotational damper mounted in the support can be considered an additional method to control system dynamics.
References
 Szmidla J. Vibrations and stability of Ttype frame loaded by longitudinal force in relation to its bolt. ThinWalled Structures, Vol. 45, 2007, p. 931935. [Publisher]
 Szmidla J. Free vibrations of a Γ type planar frame loaded by a follower force directed towards the positive pole. Vibrations in Physical Systems, Vol. 24, 2010, p. 405410. [Search CrossRef]
 Chondros T. G., Dimarogonas A. D., Yao J. Continuous cracked beam vibration theory. Journal of Sound and Vibration, Vol. 215, Issue 1, 1998, p. 1734. [Publisher]
 Rizos P. F., Aspragathos N., Dimarogonas A. D. Identification of crack location and magnitude in a cantilever beam from the vibration modes. Journal of Sound and Vibration, Vol. 138, Issue 3, 1990, p. 381388. [Publisher]
 Ostachowicz W. M., Krawczuk M. Analysis of the effect of cracks on the natural frequencies of a cantilever beam. Journal of Sound and Vibration, Vol. 150, Issue 2, 1991, p. 191201. [Publisher]
 Kukla S. Free vibrations and stability of stepped columns with cracks. Journal of Sound and Vibration, Vol. 319, 2009, p. 13011311. [Publisher]
 Sokół K., Uzny S. Instability and vibration of multimember columns subjected to Euler’s load. Archive of Applied Mechanics, Vol. 86, 2016, p. 883905. [Publisher]
 Sokół K. Linear and nonlinear vibrations of a column with an internal crack. Journal of Engineering Mechanics, Vol. 140, Issue 5, 2014, https://doi.org/10.1061/(ASCE)EM.19437889.0000719. [Publisher]
 Ratazzi A. R., Bambill D. V. B., Rossit C. A. Vibrations of a frame structure with a crack. Mecánica Computacional, Vol. 32, 2013, p. 35633574. [Search CrossRef]
 Grossi R. O., Albarracín C. M. Variational approach to vibrations of frames with inclined members, Applied Acoustics, Vol. 74, 2013, p. 325334. [Search CrossRef]
 Mamatha Y. P., Sujith Kumar S., Vivek G., Sumanth D. R., Vinay N. Analysis of crack influence on an open frame Lstructure using FEM. International Journal of Engineering Research and General Science, Vol. 4, Issue 4, 2016, p. 108117. [Search CrossRef]
 Sochacki W., Rosikon P., Topczewska S. Constructional damping mounting influence on T type frame vibrations. Journal of Vibroengineering, Vol. 15, Issue 4, 2013, p. 18661872. [Search CrossRef]
 Sochacki W., Bold M. Transverse and longitudinal damped vibration of the GAMMA type frame. Journal of Applied Mathematics and Computational Mechanics, Vol. 15, Issue 2, 2016, p. 147155. [Publisher]
 Liang J. W., Feeny B. F. Identifying Coulomb and viscous friction from free vibration decrements. Nonlinear Dynamics, Vol. 16, 1998, p. 337347. [Publisher]
 Yu Xue, Jiongming Jiang, Binping Zhao, Tinghuai Ma A selfadaptive articial bee colony algorithm based on global best for global optimization. Soft Computing, 2017, https://doi.org/10.1007/s0050001725471. [Search CrossRef]
 Qi Liu, Weidong Cai, Jian Shen, Zhangjie Fu, Xiaodong Liu, Nigel Linge A speculative approach to spatial‐temporal efficiency with multi‐objective optimization in a heterogeneous cloud environment. Security and Communication Networks, Vol. 9, Issue 17, 2016, p. 40024012. [Publisher]
 Qing Tian, Songcan Chen Crossheterogeneousdatabase age estimation through correlation representation learning. Neurocomputing, Vol. 238, 2017, p. 286295. [Publisher]
Cited By
Journal of Vibroengineering
Marta Bold, Wojciech Sochacki

2019
