Issue 
Mechanics & Industry
Volume 24, 2023



Article Number  11  
Number of page(s)  9  
DOI  https://doi.org/10.1051/meca/2023010  
Published online  17 April 2023 
Regular Article
Damping coefficient calculation method for shoreside container cranes under the seismic conditions
^{1}
School of Mechanical Engineering, Hubei University of Technology, Wuhan 430068, China
^{2}
School of Railway Locomotive and Vehicle, Wuhan Railway Vocational College of Technology, Wuhan
430205, China
^{3}
School of Mechanical Engineering, Yangtze University, Jingzhou 434023, China
^{*} email: 20190134@wru.edu.cn
Received:
27
December
2021
Accepted:
18
March
2023
Shoreside container cranes (SCC) are widely used in ship loading and unloading. However, the SCC is often damaged by the seismic wave shock. In seismic response analysis, the setting of damping coefficients has a large impact on the accuracy of the calculation results. In order to better design the SCC mechanical performance, the improved damping coefficient calculation method is proposed. The results of the simulation and the experiment are compared under the different seismic conditions. The results show that the influence of the damping coefficient calculation method on the SCC seismic response increases with strengthened seismic intensity. When the peak acceleration is small, the simulation results under different damping coefficients are almost the same and are close to the experiment values. When the peak acceleration is large, only the error between the simulation results obtained by the improved method and the experimental values is less than 5%. It is thus concluded that the improved damping coefficient calculation method is accurate and reliable.
Key words: Shoreside container crane / seismic conditions / damping coefficient / experiment
© F. Ye et al., Published by EDP Sciences 2023
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Shoreside container crane (SCC) is a piece of important loading and unloading equipment at ports, and the damage to SCC caused by the earthquake has drawn high attention from many scholars. The designer must master the dynamic performance of the structure under different working conditions, predict the dangerous parts of the structure under different earthquake conditions, and make reasonable improvements to the structure to avoid major accidents when the earthquake comes. However, only similar model tests can be carried out since onsite tests are not available for such large structures. Before the model test, seismic dynamic simulation analysis should be conducted for them, where damping coefficient has always been a difficulty in the research because reasonable setting of damping coefficient may make the simulation results more accurate and reliable.
Although many research achievements have been achieved in seismic tests for large engineering structures, they are targeted at high buildings [1–3], gravity dams [4–6] and bridges [7–9]. In order to obtain the SCC mechanical performance, Jacobs et al. [10] built a 1:20 scale model of SCC to carry out a series of finite element simulation analyses, and they achieved satisfactory results. Kanayama and Kashiwazaki [11] designed a 1:8 test model to research the multirigidbody dynamic behavior of SCC and the damage forms. Japanese researchers [12,13] built a 1:15 test model of pier foundation and the SCC, they carried out a series of seismic tests for the interaction between pier foundation and SCC. Lou et al. [14] discussed the influence of mass proportional damping matrix on the calculation results, and then the determination method of mass proportional damping coefficient was proposed. At present, the research of proportional damping hypothesis has been relatively mature, and there are many cases in the field of large structure as civil engineering and machinery. However, there is no explanation on how to determine the damping coefficient in the numerical calculation of finite elements, and the Rayleigh damping models are often considered as a basis for damping matrix construction commonly utilized in the seismic time history analysis of large structures [15].
Existing studies do not set the corresponding effective damping coefficient for the SCC. In this paper, the improved damping coefficient calculation method is proposed based on the most common four Rayleigh damping coefficient calculation methods, where a simulation analysis model is built for dynamic time history analysis with the SCC model as the experimental subject. In the time history analysis, five groups of damping coefficients obtained by the above methods are used to compare the seismic response of SCC with different damping coefficients. Finally, the shaking table tests are conducted for the model, and the results show that the improved method is more available.
2 Main components of SCC
The research object of this article is the J248 SCC structure, which is mainly welded and assembled by Q345 highstrength steel. The yield limit of the material is 345 MPa, and the elastic modulus is E = 2.06 × 10^{11}N/m^{2}, the Poisson's ratio is 0.3, the density is 7850 kg/m^{3}. The SCC mainly consists of the crane span traveling mechanism, gantry structure (mainly includes landside and seaside pillars, landside and seaside beams, horizontal beams, drawbars, and other components), front and rear girders, drawbar system, machine room, traveling trolley, and slings, etc., as shown in Figure 1.
During operation, the SCC travels along the track that is parallel to the coastal line, while the traveling trolley transfers containers along the track on the beam for loading and unloading operations.
Fig. 1 Main components of SCC structure. 
3 Calculation method for damping coefficient
The existence of damping makes the amplitude of the system in the resonance cannot be infinitely amplified, and it can consume the energy of earthquakes. Unlike the physical quantities such as mass and stiffness, the damping cannot be measured directly on the structure or expressed clearly by mathematical methods. Therefore, the damping matrix [C] is written as the sum of the linear proportions of the mass matrix [M] and the stiffness matrix [K] [16]:
where α and β are damping coefficients.
The setting of damping parameters includes damping coefficient α and β and damping ratio ξ, which needs to be set artificially in the simulation software. The value of α and β can be calculated in the following formula:
where ξ_{m} and ω_{m} are the damping ratio of the number m order vibration mode and its corresponding fixed frequency; ξ_{n} and ω_{n} are the damping ratio of the number n order vibration mode and its corresponding fixed frequency. Upon calculation of α and β, the damping ratio of other order vibration mode may be solved through the following formula:
It can be seen from equation (3) that the damping ratio in time domain calculation is frequency dependent [17]. The selection of appropriate frequency can help obtain the most reasonable α and β. Namely, the data which are closest to the actual situations are obtained through finite element simulation experiments, providing accurate and reliable data reference for the model tests.
At present, in the calculation of the dynamic response of the structure in engineering, the coefficient α and β are usually calculated using the following methods:
Method 1: Assuming that the contributions of mass matrix and stiffness matrix in Rayleigh damping is the same. Then the damping coefficient α and β can be expressed as:
The drawback of this method is that all damping within the frequency range is overestimated, resulting in relatively low dynamic response of structure under the same external dynamic load incentive.
Method 2: The sensitive frequency ranges of the selected object structure is ω_{a} and ω_{b}. The damping ratio at the boundary of frequency may be expressed with the following formula:
This method underestimates the damping between ω_{a} and ω_{b}, it also overestimates the damping beyond the frequency range. Therefore, the structural dynamic response calculated by this method is too high.
Method 3: Based on method 2, relevant scholars improved the calculation method [18], let
where ω = (α/β)^{1/2}, by substituting it into equation (3), the minimum damping ratio ξ_{min} = (αβ)^{1/2} can be obtained. The damping ratio at ω_{a} and ω_{b} may be expressed as:
ξ_{0} = (ξ_{max} + ξ_{min})/2 is defined, where ξ_{max} and ξ_{min} are the maximum and minimum damping ratios, respectively. If the calculated damping ratio ξ < ξ_{0}, then ξ = ξ_{0}. According to equation (7), α and β may be calculated. Then the underestimated damping ratio between ω_{a} and ω_{b} may be partially made up.
Method 4: Based on method 1, the two frequencies of ω_{1} and ω_{2} are used for calculating α and β. ω_{1} is the fundamental frequency of the structure; ω_{2 }= nω_{1}, where n is an odd number which is greater than ω_{e}/ω_{1}, and where ω_{e} is the main frequency of the seismic wave. The frequency characteristics of structure and the characteristics of seismic frequency are considered, but the structure damping within the ranges of ω_{1} and ω_{2} is underestimated. α and β can be expressed as [19]:
Improved method in this paper: The essential difference of the calculation method of damping coefficient mainly lies in the way of frequency selection. Hence, considering the advantages and disadvantages of the above methods, the improved Rayleigh damping coefficient calculation method is proposed. In this method, the ω_{1} as the fundamental frequency of structure in method 4 is changed to the minimum inherent frequency. Here, ω_{1} is the first order inherent frequency of structure vibration under the ground motion excitation. ω_{2 }= nω_{1}, n is an odd number which is greater than ω_{e}/ω_{1}, where, ω_{e} is the main frequency of seismic wave. After two main frequencies are determined, the method 3 is used to calculate damping coefficient α and β. The frequency characteristics of structure and spectrum characteristics of seismic waves are considered, without underestimating or overestimating the damping of structure within the ranges of ω_{1} and ω_{2}.
4 Simulation and experiment
4.1 Model parameters and conditions
3D simulation model of the SCC is established by the software ABAQUS, as shown in Figure 2a. To compare the seismic response of the SCC structure with different damping coefficients and prove the correctness of the improved calculation method, an SCC model is adopted for simulation, the key measurement points of the model are shown in Figure 2b. Among them, the connection between the model leg and the vibrating table is to simulate the real situation, that is, the real contact between the walking wheel at the bottom of the quay crane and the track.
Meanwhile, according to the real wheelrail contact model, the real wheelrail contact model of the experimental model is established, which is a reduced scale (1:15) representation of the J248 SCC structure, as shown in Figure 2c. Similarity relation of the scale model of SCC structure can be seen in Table 1. The main parameters of the shaking table are as follows: The table top of the shaking table is 1500 mm long and 1500 mm wide. The maximum acceleration in the horizontal direction is ± 50 m/S^{2}, the maximum acceleration in the vertical direction is ±30 m/S^{2}, the maximum displacement in the horizontal direction is ±200 mm, the maximum displacement in the vertical direction is ±100 mm, the maximum speed in the horizontal direction is 0.8 m/s, the maximum speed in the vertical direction is 0.8 m/s, and the maximum bearing capacity is 2 t, The frequency range is 0.1–100 Hz. Furthermore, the dynamic response data of each key point of the experiment model under seismic excitation is extracted and analyzed, and then it is compared with the simulation results.
In order to calculate damping coefficient, the frequency of finite element model is first extracted in software ABAQUS. Through hammering modal test, the first 8 order frequencies of SCC are measured, as shown in Table 2.
The scale model of SCC is placed on the shaking table for frequency sweep, thus obtaining the selfvibration frequency of SCC. In the loading setting of shaking table, the range of frequency sweep is set to 10–50 Hz. The acceleration sensors are arranged according to the layout of measuring points. Data recorded by the acceleration sensor are read after the sweep test. The frequency corresponding to the peak value is the selfvibration frequency of SCC structure. Through frequency sweep test, the vibration of the model is the most intense when the excitation frequency is close to 45.97 Hz. Therefore, 45.97 Hz can be determined as the selfvibration frequency (natural frequency) of the model. The frequency of the fourth mode in the finite element modal analysis of the model is 46.60 Hz, which is consistent with the experimental results. The simulation modal analysis results are in good agreement with the swept frequency test results, indicating that the simulation results are comparable with the test values.
Fig. 2 3D simulation model and experiment model. 
Similarity relation of the scale model.
Selfvibration frequency of SCC structure.
4.2 Simulation results
The recorded acceleration of the ELCentro (The name of ground motion) wave and Taft wave (The first seismic wave in the world that successfully records the whole process data) is set as input seismic excitation, and the peak acceleration is adjusted to 0.22 g, 0.4 g and 0.62 g, respectively. Where ELCentro and Taft are two different types of seismic waves. Their dominant frequencies are 27.27 Hz and 34.1 Hz respectively. In the simulation and shaking table experiment, the ELCentro wave and Taft wave in different peak acceleration is input directly. The scope of this study is within the linear range, so the nonlinear factors appearing in the test are constrained. Four damping coefficient calculation methods commonly used in engineering and improved coefficient calculation methods are adopted. The frequency is selected according to the above damping coefficient calculation method, and 5 groups of different damping coefficients are obtained by solving the equations, as shown in Table 3. When method 4 is adopted for calculation, n = 1, namely, ω_{2 }= nω_{1}. At this time, the calculated damping coefficient is the same with that calculated with method 1. Therefore, the finite element simulation results of T1 and T4 are the consistent, they can be merged. In the time history calculation, the measurement points of S3, S6, S11, S12 on the model are monitored.
As shown in Figure 3, it is clear that the trend of the acceleration time domain curve of the measurement point is similar under the different damping parameter settings, and the main difference is the magnitude of the peak acceleration. The maximum peak acceleration is 2721.4 mm/s^{2} under T1 and T4, which is the minimum peak acceleration among several groups of parameters, it is consistent with the calculation characteristics of Method 1. The maximum peak acceleration is 3348.4mm/s^{2} under T2, which is the maximum peak acceleration among several groups of parameters, it is consistent with the calculation characteristics of Method 2. The maximum peak accelerations under T3 and T5 are respectively 3099.8 mm/s^{2} and 2967.8 mm/s^{2}, which is between Method 1 and Method 2.
To clearly show the dynamic effects of damping coefficients on different test points on the model, the acceleration amplification coefficient envelopes of the model on measurement points S3, S6, S11, and S12 are plotted in Figure 4.
From Figure 4, under ElCentro and Taft seismic wave excitation with peak acceleration adjusted to 0.22 g, 0.4 g and 0.62 g, the peak acceleration of measuring points is different under different damping coefficient settings (T1T5), which can be reflected in the acceleration amplification factor. According to the acceleration amplification factor at the measurement points, the deviation between the simulation value of different models and the T5 can be calculated. The deviation is more obvious under the action of ElCentro seismic wave with strength of 0.22 g. The deviation between T2 and T5 at measuring point S3 is 7.81%. Under Taft seismic wave, the deviation is smaller, and the maximum value is 2.29% between T1, T4 and T5. It shows that the damping coefficient has little effect on the SCC when the peak acceleration is adjusted to 0.22 g, and the ElCentro seismic wave model is more sensitive to the damping coefficient than Taft seismic wave. Compared with 0.22 g earthquake excitation, the deviation of acceleration amplification coefficient envelope under 0.4 g earthquake is enlarged. Similarly, the deviation is obvious under ElCentro seismic wave, in which the deviation between T2 and T5 at measuring point S3 reaches 12.9%, and the deviation is smaller under Taft seismic wave, with the maximum value of 2.73% between T1, T4 and T5. When the peak acceleration is adjusted to 0.62 g, the deviation is obvious under ElCentro seismic wave, where the deviation between T2 and T5 reached 17.9% at the measuring point S3. The deviation is small under Taft seismic wave, and the maximum value is 3.76% between T1, T4 and T5.
Damping coefficient of the model.
Fig. 3 The time domain characteristic of S3 point under the condition of EL Centro wave 0.22g (a) T1 and T4; (b) T2; (c) T3; (d) T5. 
Fig. 4 Envelope curve of acceleration amplification coefficient with the seismic wave (0.22 g; 0.4 g; 0.62 g). 
4.3 Experiment results
Whether the calculation method of damping coefficient of SCC is reliable can be verified by shaking table experiment. The dynamic response data of key points under earthquake excitation are extracted and compared with the previous simulation analysis results. If the experiment values are in great agreement with the simulation values obtained by using the improved damping coefficient calculation method, it shows that the improved method is accurate and reliable.
Bolts are used to form a fixed connection with the shaking table at the support of the model legs, which is similar to the boundary conditions in civil engineering. In practice, the boundary condition between the SCC and the port's ground is the contact between the supporting wheel of crane and the port's track. The edited seismic wave is input into the shaking table control system. ElCentro and Taft seismic wave is set to 0.22 g, 0.4 g and 0.62 g.
Figure 5 shows the time domain value of acceleration measured at the points of S3 and S6 of the SCC model under the excitation of 0.22 g ELCentro wave.
The extreme value of acceleration at each point in the test is extracted to calculate the amplification coefficient of its acceleration. Figure 6 is the envelope curve of acceleration amplification coefficient between T1T5 and experiment points.
Under El and Taft waves, the simulation calculation results are compared with the experimental values. In general, the closest finite element calculation result to the test value is the damping parameter obtained with T5 method. That is, the new calculation method proposed in this article. When the peak acceleration is 0.22 g, the deviation between the simulation value (T1–T5) and experimental value is relatively minor, and it is indicated that the five methods are applicable to seismic simulation calculation of the SCC model. However, with the peak acceleration increases, the deviation between simulation value and the experimental value is gradually enlarged. When the peak acceleration is adjusted to 0.62 g, there are measurement point models with deviation exceeds 10% with the calculation method T1T4, while the deviation between the simulation results and experimental value are within 5% with the new calculation method T5. It is clear that the closest simulation results to experimental value are the damping parameters obtained with the T5, or namely, the new calculation method proposed herein. Therefore, the improved damping coefficient calculation method is accurate and reliable.
Fig. 5 Acceleration test value of measurement point under 0.22 g ELCentro wave (a) S3; (b) S6. 
Fig. 6 Envelope curve of acceleration amplification coefficient with the seismic wave (0.22 g; 0.4 g; 0.62 g). 
5 Conclusions
Through the improved damping coefficient calculation method, the simulation results and experimental results of SCC seismic response under different seismic conditions are compared. The main conclusions are summarized as follows:
Four calculation methods of damping coefficient α and β in the dynamic response calculation of engineering structures are summarized. Considering the advantages and disadvantages of the above methods, an improved Rayleigh damping coefficient calculation method is proposed. The calculation method of damping coefficient has a great influence on the SCC mechanical performance analysis with different seismic wave excitation parameters, the deviation of the acceleration amplification coefficient simulation results among the method is up to 17.9%.
The reliability of the calculation method of SCC damping coefficient is verified by shaking table experiment. The dynamic response data of key points with earthquake excitation are extracted and compared with the simulation results. The improved calculation method of damping coefficient is verified, which can be used to predict SCC mechanical performance, as the deviation between the simulation results and experimental value are within 5%.
From the above findings and analysis, the application of the improved damping coefficient calculation method provide an effective way for the study of SCC performance under seismic conditions.
Funding
This study was funded by National Natural Science Foundations of China under Grants 52003078 and the Hubei Provincial Department of Education Q20211406.
Conflict of interest
The authors declare that no conflicts of interest exist in this manuscript.
References
 J.P. Moehle, Seismic analysis, design, and review for tall buildings, Struct. Des. Tall Spe Build. 15, 495–513 (2006) [CrossRef] [Google Scholar]
 D.W. Ko, H.S. Lee, Shaking table tests on a highrise RC building model having torsional eccentricity in soft lower storeys, Earthq. Eng. Struct. Dyn. 35, 1425–1451 (2006) [CrossRef] [Google Scholar]
 H.S. Lee, D.W. Ko, Seismic response characteristics of highrise RC wall buildings having different irregularities in lower stories, Eng. Struct. 29, 3149–3167 (2007) [CrossRef] [Google Scholar]
 A.M. Riyaz, C.A. Taylor, An experimental investigation into earthquakeinduced failure of medium to low height concrete gravity dams, Earthq. Eng. Struct. Dyn. 24, 373–393 (1995) [CrossRef] [Google Scholar]
 R. Tinawi, P. Leger, M. Leclerc, G. Cipolla, Seismic safety of gravity dams: from shake table experiments to numerical analysis, J. Struct. Eng. 126, 518–529 (2000) [CrossRef] [Google Scholar]
 D.G. Zou, B. Xu, X.J. Kong, Study of influence of different methods for calculating Rayleigh damping coefficient on high earthrock dam seismic response, Rock Soil Mech. 32, 797–803 (2011) [Google Scholar]
 M.A. HosamEddin, A.M. AbdelGhaffar, Seismic energy dissipation for cablestayed bridges using passive devices, Earthq. Eng. Struct. Dyn. 23, 877–893 (1994) [CrossRef] [Google Scholar]
 K.R. Karim, F. Yamazaki, Effect of earthquake ground motions on fragility curves of highway bridge piers based on numerical simulation, Earthq. Eng. Struct. Dyn. 30, 1839–1856 (2001) [CrossRef] [Google Scholar]
 B.B. Chen, Z.R. Feng, H.M. You, Seismic response analysis of longspan cablestayed bridge with nonlinear viscous dampers under different periodic seismic waves, Civ. Eng. J. 27, 114–131 (2018) [Google Scholar]
 L.D. Jacobs, D.R. Reginald, T.L. Roberto, Seismic behavior of a Jumbo container crane including uplift, Earthq. Spectra 27, 745–773 (2011) [CrossRef] [Google Scholar]
 T. Kanayama, A. Kashiwazaki, Study on the dynamic behavior of container cranes under strong earthquakes. ASME PVP, Seismic Eng. 364, 276–284 (1998) [Google Scholar]
 S. Yamamoto, T. Sugano, T. Tanabe, The shaking table test of the interaction between the pier type wharf and the container crane during earthquakes. Trans. Jpn. Soc. Mech. Eng. 68, 3209–3216 (2002) [CrossRef] [Google Scholar]
 T. Eto, Model experiment on dynamic interaction of a container crane and a pier in a largescale earthquake, Kuren. 37, 20660–20665 (1999) [Google Scholar]
 M.L. Lou, X.G. Shao, Coefficient selection of damping matrix in explicit calculation of soil layer seismic response, J. Tongji. Univ. 41, 1126–1132 (2013) [Google Scholar]
 Z.Q. Song, C.H. Su, Computation of rayleigh damping coefficients for the seismic analysis of a hydropowerhouse, Shock Vib. 2046345 (2017) [Google Scholar]
 A.K. Chopra, R.K. Goel, A modal pushover analysis procedure for estimating seismic demands for buildings, Earthq. Eng. Struct. Dyn. 31, 561–582 (2002) [CrossRef] [Google Scholar]
 Z.G. Xu, X.L. Du, C.S. Xu, C.Y. Zang, J.W. Jiang, Comparison of determination methods of site Rayleigh damping coefficients in seismic responses analysis of underground structures, Rock Soil Mech. 40, 4838–4847 (2019) [Google Scholar]
 N. Yoshida, S. Kobayashi, I. Suetomi, Equivalent linear method considering frequency dependent characteristics of stiffness and damping, Soil Dyn. Earthq Eng. 22, 205–222 (2002) [CrossRef] [Google Scholar]
 Z. Li, C. Xu, S.Y. Wu, Calculating the Rayleigh damping coefficient under seismic dynamic simulation of a quayside container crane, Aust. J. Mech. Eng. 1–10 (2019) [Google Scholar]
Cite this article as: F. Ye, H. Deng, C. Fang, X. Fu, Z. Li, Damping coefficient calculation method for shoreside container cranes under the seismic conditions, Mechanics & Industry 24, 11 (2023)
All Tables
All Figures
Fig. 1 Main components of SCC structure. 

In the text 
Fig. 2 3D simulation model and experiment model. 

In the text 
Fig. 3 The time domain characteristic of S3 point under the condition of EL Centro wave 0.22g (a) T1 and T4; (b) T2; (c) T3; (d) T5. 

In the text 
Fig. 4 Envelope curve of acceleration amplification coefficient with the seismic wave (0.22 g; 0.4 g; 0.62 g). 

In the text 
Fig. 5 Acceleration test value of measurement point under 0.22 g ELCentro wave (a) S3; (b) S6. 

In the text 
Fig. 6 Envelope curve of acceleration amplification coefficient with the seismic wave (0.22 g; 0.4 g; 0.62 g). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.