Applications of bispectrum analysis to inspection of fatigue damage in quay cranes

Li Cui1 , Maosen Cao2 , Hao Zhang3 , Shuigen Hu4 , Jie Wang5

1, 2, 3, 4Anhui Provincial International Joint Research Center of Data Diagnosis and Smart Maintenance on Bridge Structures, Chuzhou University, Anhui, People’s Republic of China

1, 2, 3, 4College of Civil and Architecture Engineering, Chuzhou University, Anhui, People’s Republic of China

5Intelligent Transportation and Intelligent Construction Engineering Research Center, Jiangsu Dongjiao Intelligent Control Technology Group Co., Nanjing, China

2Corresponding author

Vibroengineering PROCEDIA, Vol. 44, 2022, p. 28-32. https://doi.org/10.21595/vp.2022.22881
Received 1 August 2022; received in revised form 11 August 2022; accepted 16 August 2022; published 25 August 2022

58th International Conference on Vibroengineering, Ventspils, Latvia, August 25-26, 2022

Copyright © 2022 Li Cui, et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Editor's Pick
Table of Contents Download PDF Acknowledgements References
Cite this article
Views 2
Reads 1
Downloads 70
CrossRef Citations 0
Abstract.

The quay crane is a special mechanical equipment which can lift the container vertically and horizontally through intermittent circulation movement. Fatigue crack is one of the main kinds of damage in quay cranes. The generation and development of fatigue crack are invisible, which makes the identification of fatigue cracks a great challenge. From the dynamic point of view, fatigue cracks have nonlinear dynamic characteristics beyond the linear dynamic category. The traditional damage detection methods based on linear hypothesis are difficult to effectively recognize the nonlinear dynamic characteristics of fatigue cracks, so they cannot be effectively used for structural fatigue crack detection. In order to overcome the limitations of the existing linear damage detection methods, bispectrum analysis, a nonlinear dynamic analysis method, is introduced in this thesis. By drawing on the use of bispectrum analysis to describe the nonlinear dynamic characteristics of rotating structures, a fatigue crack detection method based on nonlinear dynamic characteristics of quay cranes is developed.

Keywords: fatigue damage, bispectrum analysis, nonlinear dynamic characteristics, quay crane.

1. Introduction

Quay Cranes are in a very important and unique position in the increasingly fierce international trade and port competition. Quay Cranes are typical welded metal structures. In addition, when the quay crane is working, the trolley is constantly started and braked, or the suspended load is shuttling between the front and rear main beams at high speed, resulting in alternating stress in the structure and causing dislocation movement accumulation and fatigue cracks. However, the current main health monitoring methods of quay cranes are still based on the visual inspection of experienced personnel and magnetic particle inspection, which obviously cannot meet the huge demand of quay cranes [1].

The Korean scholar Hyawn [2] studied the stress change of the quay cranes under the wind load, and demonstrated that even a very small wind load will cause the stress on the safety device of the quay cranes to increase rapidly, which is related to the increase rate of the wind load. Jacobs [3] evaluated and characterized the critical elastic response and derailment vibration response of portal frame of large quay cranes under earthquake load through the shaking table test of a 1:20-scaling quay crane model. Vöth [4] focused on the safety problems of the quay cranes caused by the high-speed operation of the trolley. Most of the existing studies focus on the safety assessment and reliability research of quay cranes under special loads, especially the earthquake loads. The damage detection of quay cranes under long-term service conditions is rarely considered. The special load will cause great damage to quay cranes, but it is worth noting that on the one hand, the occurrence of the special load is very accidental, and it is very likely that a quay crane will not be affected by such special loads as earthquake loads during the whole service period. On the other hand, most of the faults and damages of quay cranes are caused by long-time and high-speed operation during daily service. Therefore, while paying attention to the damage of the structure caused by special loads, it is also necessary to develop a structural damage detection method based on long-term and daily monitoring data processing of the quay cranes.

2. Bispectrum analysis method

It is assumed that the mean value of real random signal x(n) is zero and k-order is stationary. If the k-order cumulant ckx(τ1,τ2,,τk-1) of x(n) exists, the k-order spectrum is defined as the (k-1)-dimensional Fourier transform of ckx(τ1,τ2,,τk-1):

(1)
Skxω1,ω2,,ωk-1=τ1=-τk-1=-τ1,τ2,,τk-1exp-ij=1k-1ωjτj.

When k= 3, the third-order spectrum S3x is also called bispectrum, expressed in Bx(ω1,ω2).

The estimated bispectrum[5-6]can be expressed as:

(2)
B(ω1,ω2)=X(ω1)X(ω2)X*(ω1+ω2),

where X(k)=n=0N-1x(n)e-i2πNnk(0kN-1) is the Fourier transform determined by the sequence of principal values of x(n); * means complex conjugate.

ω1 and ω2 in B(ω1,ω2) are two mutually independent frequencies, which can represent non minimum phase information. B(ω1,ω2) is a double periodic function with the same period. The period is 2π, and satisfies:

(3)
Bω1,ω2=B*-ω1,-ω2=B-ω1-ω2,ω2==B-ω1-ω2,ω1=Bω2,ω1=B*-ω2,-ω1=Bω2,-ω1-ω2=Bω1,-ω1-ω2.

There are six symmetry lines which are ω1=ω2, 2ω1=-ω2, -ω1=2ω2, ω1=-ω2, ω1=0 and ω2=0 in B(ω1,ω2). Since the bispectrum satisfies such high symmetry, all the bispectrum can pass through the symmetric region ω20, ω1ω2, ω1+ω2π.

3. The finite element simulation

In this paper, an 80 t-65 m quay crane is taken as the prototype to build the finite element model. The relevant parameters are shown in Table 1. The material of the quay crane is Q345 steel, and the material property is: elastic modulus E= 2.06×1011 Pa, Poisson's ratio υ = 0.3, density ρ = 7.85×103 kg/m3. In this paper, the quay crane is set in the working state, that is, the crane is completely locked and in the windless environment. Therefore, the front girder can be regarded as a complex cantilever beam model, and fixed constraints will be set at the bolt holes on the land side. This is to simplify the model and calculation.

Table 1. Main parameters of the model

Front extension (m)
Back extension (m)
Crane gauge (m)
Trolley gauge (m)
Base distance (m)
Load capacity (t)
65
30.5
35
6.3
20.5
Sling
Cliver
65
80

A 50 %- degree breathing crack is set at five positions as shown in Fig. 1. respectively. A simple harmonic excitation Fsinωt is set on the upper surface of the free end near the sea side. F= 50000 N, ω=f1/2 (f1= 0.353 Hz is the first-order natural frequency of the front girder with a breathing crack). The acceleration response of the structure is extracted at the lower surface of the free end near the sea side, and bispectrum analysis is performed. The results are shown in Fig. 2.

When the breathing crack transits from the fixed end near the land side to the free end near the sea side, the second harmonic in the bispectrum gradually changes from obvious to not obvious, but the amplitude of the second harmonic is the smallest in the middle section of the front girder; When the crack is close to the fixed end on the land side and the free end on the sea side, the amplitude of the second harmonic is larger. However, it should be noted that this change rule is not invariable, and may generally change due to the difference between the position where the excitation is applied and the position where the dynamic response is extracted.

Fig. 1. Diagram of damage position on the front girder of the quay bridge

Diagram of damage position on the front girder of the quay bridge

Fig. 2. Three-dimensional bispectrum when the damage is at different position

Three-dimensional bispectrum when the damage is at different position

a) Position 1

Three-dimensional bispectrum when the damage is at different position

b) Position 2

Three-dimensional bispectrum when the damage is at different position

c) Position 3

Three-dimensional bispectrum when the damage is at different position

d) Position 4

Three-dimensional bispectrum when the damage is at different position

e) Position 5

A breathing crack is set at position 1 near the fixed end, and the crack depth p is set as 10 %, 20 %, 30 %, 40 %, and 50 %, respectively to study the change of the bispectral nonlinear dynamic characteristics under different damage degrees. A simple harmonic excitation Fsinωt is set on the upper surface of the free end near the sea side, where the excitation size F is set to 50 000 N and the excitation frequency ω is set to f1/2. The acceleration response of the structure is extracted at the lower surface of the free end near the sea side, and bispectrum analysis is performed. The results are shown in Fig. 3.

Fig. 3. Bispectrum when a) p= 10 %; b) p= 20 %; c) p= 30 %; d) p= 40 %; e) p= 50 %

Bispectrum when a) p= 10 %; b) p= 20 %; c) p= 30 %; d) p= 40 %; e) p= 50 %

a)p= 10 %

Bispectrum when a) p= 10 %; b) p= 20 %; c) p= 30 %; d) p= 40 %; e) p= 50 %

b)p= 20 %

Bispectrum when a) p= 10 %; b) p= 20 %; c) p= 30 %; d) p= 40 %; e) p= 50 %

c)p= 30 %

Bispectrum when a) p= 10 %; b) p= 20 %; c) p= 30 %; d) p= 40 %; e) p= 50 %

d)p= 40 %

Bispectrum when a) p= 10 %; b) p= 20 %; c) p= 30 %; d) p= 40 %; e) p= 50 %

e)p= 50 %

Fig. 4. Variation trend of the second harmonic amplitude with crack depth

Variation trend of the second harmonic amplitude with crack depth

Fig. 3 shows that the second harmonic component of the bispectrum becomes more and more obvious with the deepening of the damage degree, and the amplitude is also larger and larger. Among them, when the damage degree is 10 % and 20 %, the second harmonic component can be observed, but it is not obvious compared with the case of 50 % damage depth, which may be related to the monitoring position and normalization processing method. Extract the second harmonic amplitude under different crack depths, and the change is shown in Fig. 4. It can be seen from the figure that the amplitude of the second harmonic of the bispectrum monotonically increases with the crack depth. Therefore, the second harmonic component of the bispectrum can be used as a characteristic index to identify the fatigue damage degree.

4. Conclusions

In this paper, the bispectrum dynamic analysis is carried out on the quay crane which is under working condition. The results show that bispectrum analysis can extract the nonlinear dynamic characteristics of the crane structure with fatigue cracks. The bispectral nonlinear dynamic characteristics are positively correlated with the depth of fatigue crack in the structure. The bispectral nonlinear dynamic characteristics no longer show a single change law with the change of crack position. The research results show that the proposed fatigue crack detection method based on bispectrum analysis can capture the nonlinear dynamic characteristics of fatigue cracks in the structure of quay cranes under different working conditions, and provides a new technology for fatigue crack identification of quay cranes. It must be noted that, bispectrum analysis method has shown unprecedented advantages in the fault detection of structure components, especially rotating components, but it is rarely used in the damage detection of large and complex structures. Therefore, this paper only proposes new possibilities for the fatigue crack detection of quay cranes. However, if this method is really to be used in the actual structure monitoring, it must be further studied.

Acknowledgements

This work is supported by Anhui provincial international joint research center of data diagnosis and smart maintenance on bridge structures (No. 2021AHGHYB05), the University scientific research project of Anhui Provincial Department of Education (Nos. KJ2021A1097 and KJ2021A1098), the Key R&D Project of Anhui Science and Technology Department (202004b11020026), and Nanjing International Joint Research and Development Program (No. 202112003).

References

  1. W. Pei, K. Q. Ding, and D. J. Wu, “Method of Safety Health Monitoring and Damage Alarming for Cranes,” Mechanical Engineering and Automation, Vol. 6, pp. 120–122, 2010, https://doi.org/10.3969/j.issn.1672-6413.2010.06.045 [Search CrossRef]
  2. K. D. Hyawn, O. B. Cheol, and C. I. Sik, “Effect of wind velocity variation on Safety factor of container crane tie-downs,” KSCE Journal of Civil and Environmental Engineering Research, Vol. 25, No. 6A, pp. 1127–1131, 2005. [Search CrossRef]
  3. L. D. Jacobs, B. D. Kosbab, R. T. Leon, and R. Desroches, “Seismic behavior of a jumbo container crane including uplift,” Earthquake Spectra, Vol. 27, No. 3, pp. 745–773, Aug. 2011, https://doi.org/10.1193/1.3610238 [Publisher]
  4. S. Vöth, “Safety concepts for container crane hoists,” in Proceedings of the 17th ITI-Symposium, pp. 127–133, 2014. [Search CrossRef]
  5. P. K. Ajmera, N. S. Nehe, D. V. Jadhav, and R. S. Holambe, “Robust feature extraction from spectrum estimated using bispectrum for speaker recognition,” International Journal of Speech Technology, Vol. 15, No. 3, pp. 433–440, Sep. 2012, https://doi.org/10.1007/s10772-012-9153-5 [Publisher]
  6. M. S. Cao, Y. M. Xu, K. Q. Ding, and E. Manoach, “Bispectral dynamics features for characterizing structural fatigue damage,” Journal of Vibroengineering, Vol. 20, No. 5, pp. 2073–2084, Aug. 2018, https://doi.org/10.21595/jve.2018.19963 [Publisher]
JVE Journals is rebranding to Extrica

Inspired by innovations from the previous century and the rapid growth during the last years, we are improving for excellence in your publishing experience

Read to know more
JVE Journals
Extrica