Authors: Kun Zhang, Zhengxian Sun, Jinpeng Su, Xuntao Wei, Mingchao Du, Hongyue Chen
Categories: Article, Engineering, Mechanical engineering
Source: Scientific Reports
Authors: Kun Zhang, Zhengxian Sun, Jinpeng Su, Xuntao Wei, Mingchao Du, Hongyue Chen
A round link chain subject to axial dynamic loads composes a nonlinear viscoelastic system. Unlike the classical pounding problems, the round link chain will not only suffer linear elastic deformation, but also nonlinear plastic or impacting deformation. Based on theoretical formulation and experiments, a new approach is presented in this paper to model and identify the nonlinear dynamic parameters, namely the stiffness and damping for the round link chain. With linear deformation, nonlinear deformation and energy dissipation considered, a modified nonlinear viscoelastic model is developed to describe the vibrational behavior of the chain with numbers of round links. The linear elastic model and impacting model are combined to derive the equivalent nonlinear stiffness, while experiments and the least square fitting method are employed to identify the nonlinear damping according to the modified nonlinear viscoelastic model. The influences of the key parameters such as the length of the chain, elastic module and loading frequency on the dynamic stiffness and damping are investigated. Another test is performed to validate the identification model and good agreements are observed.
As key components of shipping/lifting hoists or bulk material handling machines, round link chains are wildly used in marine, mechanical, mining and civil engineering. Investigations on the dynamic characteristics are of great significance for good performance of the corresponding equipment and machines^1–3^.
Numerous researchers have been devoting to the static and dynamic analysis for kinds of round link chains. Ming et al.^4^ derived an equation to relate the contact area of the round link chains to the static contacting stress. Li et al.^5^ calculated the maximum stress and pressure distribution of the contact area of the ring chain by Hertz theory, and obtained the maximum stress of the meshing contact area of the sprocket chain. Bian et al.^6^ established the mathematical equation of circular chain structure, deduced two static analysis models of chain-ring contact, simulated and analyzed the collision process between circular chains, and revealed the fatigue fracture mechanism and fatigue crack propagation law of circular chain. Li et al.^7^ established a virtual prototype model for dynamics simulation of chain drive system, simulated and analyzed the start-up of ring chain load after shutdown, and obtained the variation law of kinematics and dynamic behavior response of meshing contact. Diao et al.^8^ conducted photo-elasticity experiments and finite element (FE) analysis on the contact stress for the round link chains. Wang et al.^9^ used time-varying dynamic analysis method to obtain the dynamic tension distribution of the chain of heavy scraper conveyor. The finite element analysis was carried out on the three-dimensional contact between the adjacent chains of the straight segment and the bending segment to obtain the three-dimensional stress distribution of the chain.
As reviewed, since the round link chain usually suffer heavy axial load to convey bulk material or to transfer motions, reasonably more attentions are paid to the dynamic characteristics in the axial direction. When the round chain is subjected to harmonic or impacting excitations, impacting, friction and even plastic deformation^10^ can be observed. The round link chain composes a typical viscoelastic system. The damping effects and the nonlinear stiffness should be considered in the dynamic analysis. The Kelvin–Voight model, which was firstly developed for the pounding problems^11^, is widely used to take the energy dissipation into account. However, quite similar to the pounding problems, the contact, friction and the plastic deformation are actually all nonlinear. Thus, the linear model cannot totally consider the nonlinear factors for the damping and stiffness. To address the issues for the pounding problems, kinds of non-linear models^12^, such as the Hertzdamp non-linear model^13^ and nonlinear Hunt–Crossley model^14^ have been proposed. All the nonlinear models are quit efficient and accurate for lots of cases and they are also well employed to model other viscoelastic systems^15–17^.
However, the stiffness and damping coefficients can hardly be determined for complex practical systems. Thereby, considerable work has been done to develop identification methods and experimental techniques to identify or measure the dynamic parameters^18^. Based on the servo motor current and corresponding position deflection, Yang et al.^19^ developed a new identification methodology for the accurate joint stiffness of heavy duty robots. Fan et al.^20^ employed the expected trajectory and external load to identify the stiffness defect of a parallel manipulator. Jin et al.^21^ extended the Hertzdamp model to derive the equivalent stiffness and damping for auto eliminating clearance auxiliary bearing devices. Liu et al.^22^ established a variable motion mapping method to enhance stiffness identification of the remote object. Based on the output properties, Zhao et al.^23^ developed an algorithm to reconstruct simultaneously an anti-damping coefficient for some anti-stable systems, yet the algorithm is not suitable for the Kelvin–Voigt damping and some other spatially varying damping.
With the devolvement of experimental techniques, more and more test-based composite methods are proposed to obtain the dynamic stiffness and damping. Wang et al.^24^ designed a measurement system and derived a tangential contact stiffness model to identify the friction damping and the key parameters of tangential contact stiffness in the transition process. Xu et al.^25–27^ developed a laser-based deformation modeling methods and successfully identified complex deformations for composite structures. Lee et al.^28^ utilized curved-edge diffraction to improve the dimensional measurement technique, and the dynamic parameters for a ball bearing spindle system was obtained. Domenico Lisitano et al.^29^ formulated a so called layer method for damping matrix identification using the experimental receptance- -matrix data together with physical connectivity constraints. The method could overcome the ill-conditions and is stable for various frequency ranges. Budak et al.^30^ put forward a new identification framework for the process parameters in turning and milling. The process damping was obtained from chatter tests and the indentation force coefficient is identified with an energy method. Ben Romdhane et al.^31^ performed experiments to determine the loss factors of a non-obstructive particle system with large number of freedoms. When the damping properties were obtained, they employed an equivalent viscous damping to predict the dynamic characteristics of a beam with the non-obstructive particle system.
The above review indicates that the nonlinear viscoelastic models can be successfully used to study the dynamic characteristics of a wide range of viscoelastic systems, not just limited to the pounding problems. Even though the stiffness and damping coefficients are hardly determined theoretically for complex systems, including the round link chain, the model still can provide physical insights into the viscoelastic problems and can be readily employed to develop experiment-based composite identification methods. Moreover, as presented in the latest literature, the stiffness, linear or nonlinear, can be formulated in a theoretical way, though the expressions of the nonlinear ones are quite complex. However, damping properties for complicated systems are usually obtained using tests, because the energy dissipating mechanism for the systems is comprehensive and sophisticated. Hence, this paper is intended to combine the theoretical modeling and experiments for dynamic parameter identification of a round link chain with numbers of round links. An energy method is utilized to deduce the nonlinear stiffness, including the linear elastic one and the one for the indentations. Based on a modified nonlinear viscoelastic model and testing data, the nonlinear damping for the round link chain is given. The identified results is validated with a new experimental case and the effects of the material and geometrical properties are also examined.
As shown in Fig. 1, when a round link chain is subjected to large dynamic loads in the axial direction, impacting, elastic and plastic deformations will occur. It composes a typical viscoelastic system. In order to accurately simulate stiffness effects and energy dissipation, a nonlinear vibrational system is proposed for a round link chain with n links in this paper. As shown in Fig. 2, the model is based on the non-linear viscoelastic model in Ref12. In the model, the linear stiffness and nonlinear contact stiffness are modeled as an equivalent nonlinear stiffness. The energy dissipation in the contact process is expressed with equivalent damping. The total mass of the n links is simulated with a concentrated mass m.Figure 1Force diagram of the round link chain.Figure 2Equivalent dynamic model for a round link chain with n links.
In this way, the governing equations can be expressed as1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F = F_{m} + F_{k} + F_{c}
\usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{m} = m\ddot{x} $$\end{document}Fm=mx¨ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x $$\end{document}x is the displacement of the chain. The elastic force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{k} = k_{e} \left( x \right)x $$\end{document}Fk=kexx, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ k_{e} \left( x \right) $$\end{document}kex is the equivalent nonlinear stiffness. The dissipated force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{c} = c_{e} \left( x \right)\dot{x} $$\end{document}Fc=cexx˙ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ c_{e} \left( x \right) $$\end{document}cex is the nonlinear equivalent damping. In the model, the total mass *m* is known, while the nonlinear elastic and dissipated forces are needed to be determined. ### Formulation of the equivalent nonlinear elastic force As analyzed above, the equivalent nonlinear elastic force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{k} $$\end{document}Fk can be divided into two parts, namely the linear elastic force and the nonlinear elastic force. The linear elastic force is deduced by the linear elastic deformation of the chain, while the nonlinear elastic force is due to the contact between the adjacent round links. One link of the round link chain is taken and analyzed to derive \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{k} $$\end{document}Fk, as shown in Fig. 3. The round link is assumed to comprise homogeneous and isotropic elastic materials.Figure 3Force diagram of a round link in the chain, (**a**) geometry of the round link and (**b**) simplified force diagram. For the linear elastic deformation of the round link, the assumption of infinitesimal strain theory is employed. Due to the deformation coupling of the curved segment BC in the round link, the application of axial load \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F $$\end{document}F will introduce axial force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ N $$\end{document}N, bending moment \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ M $$\end{document}M and shear force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ Q $$\end{document}Q acting on the straight segment AB. Considering the symmetry of the round link and the axial load, the forces and moment can be defined as2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left\{ \begin{gathered} M = M_{1} \hfill \\ N = \frac{{F_{k} }}{2} \hfill \\ Q = 0 \hfill \\ \end{gathered} \right. $$\end{document}M=M1N=Fk2Q=0 The axial force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ N_{c} $$\end{document}Nc, the moment \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ M_{c} $$\end{document}Mc, and the shear force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ Q_{c} $$\end{document}Qc acting on the curved segments are related to the subtended angle \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \theta $$\end{document}θ and can be expressed as3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left\{ \begin{gathered} N_{C} = \frac{{F_{k} }}{2}\cos \theta \hfill \\ M_{C} = M_{1} - \frac{{F_{k} }}{2}R\left( {1 - \cos \theta } \right) \hfill \\ Q_{C} = - \frac{{F_{k} }}{2}\sin \theta \hfill \\ \end{gathered} \right. $$\end{document}NC=Fk2cosθMC=M1-Fk2R1-cosθQC=-Fk2sinθ The total potential energy of the round link is composed of the potential energy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ U_{AB} $$\end{document}UAB stored in the straight segment and the potential energy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ U_{BC} $$\end{document}UBC for the curved segment. Thus, the total potential energy of the round link chain \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ U $$\end{document}U can be obtained as4\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \begin{gathered} U = U_{AB} + U_{BC} = \int_{0}^{L} {\left( {\frac{{N^{2} }}{{2EA_{S} }} + \frac{{M_{1}^{2} }}{{2EJ_{1} }}} \right)} dx \hfill \\ \int_{0}^{{\frac{\pi }{2}}} {\left( {\frac{{M_{C} \left( \theta \right)^{2} }}{{2EJ_{2} }} + \frac{{N_{C} \left( \theta \right)^{2} }}{{2EA_{S} }} + \frac{{Q_{C} \left( \theta \right)^{2} }}{{2GA_{S} }}} \right)} Rd\theta \hfill \\ \end{gathered} $$\end{document}U=UAB+UBC=∫0LN22EAS+M122EJ1dx∫0π2MCθ22EJ2+NCθ22EAS+QCθ22GASRdθ where, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ E $$\end{document}E is the elastic modulus and the shear modulus \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ G = {E \mathord{\left/ {\vphantom {E {\left( {2\left( {1 + \mu } \right)} \right)}}} \right. \kern-\nulldelimiterspace} {\left( {2\left( {1 + \mu } \right)} \right)}} $$\end{document}G=E/21+μ. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mu $$\end{document}μ is the Poisson’s ratio.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ A_{S} $$\end{document}AS is the cross-sectional area. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ J_{1} $$\end{document}J1 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ J_{2} $$\end{document}J2 are the elastic section modulus of the straight and curved segments, respectively. They can be determined as5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \begin{gathered} J_{2} = J_{1} + \alpha^{2} A_{S} \hfill \\ J_{2} = \frac{{\pi d^{4} }}{64} + \alpha^{2} A_{S} \hfill \\ \alpha = R - d^{2} \left( {8R\left( {1 - \sqrt {1 - \left( \frac{d}{2R} \right)^{2} } } \right)} \right)^{ - 1} \hfill \\ \end{gathered} $$\end{document}J2=J1+α2ASJ2=πd464+α2ASα=R-d28R1-1-d2R2-1 Since the geometry and axial load are assumed to be symmetric, the deformation condition at point A is given as6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \theta_{A} = \frac{\partial U}{{\partial M_{1} \left( t \right)}} = 0 $$\end{document}θA=∂U∂M1t=0 Substituting Eq. (4) into Eq. (6), the moment \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ M $$\end{document}M is obtained as7\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ M_{1} = \frac{{\left( {\pi - 2} \right)R^{2} J_{1} F_{k} }}{{4J_{2} L + 2J_{1} \pi R}} = k_{0} F_{k} $$\end{document}M1=π-2R2J1Fk4J2L+2J1πR=k0Fk where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ k_{0} = {{\left( {\pi - 2} \right)J_{1} R^{2} } \mathord{\left/ {\vphantom {{\left( {\pi - 2} \right)J_{1} R^{2} } {\left( {2LJ_{2} + \pi RJ_{1} } \right)}}} \right. \kern-\nulldelimiterspace} {\left( {2LJ_{2} + \pi RJ_{1} } \right)}} $$\end{document}k0=π-2J1R2/2LJ2+πRJ1. Substituting Eq. (7) into Eq. (4), the total potential energy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ U $$\end{document}U can be obtained as8\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ U = \left( \begin{aligned} \frac{{2k_{0}^{2} L}}{{EJ_{1} }} + \frac{L}{{2EA_{S} }} + \frac{R\pi }{{8EA_{S} }} + \frac{R\pi }{{8GA_{S} }} \hfill \\ + \frac{R}{{2EJ_{2} }}\left( {\frac{\pi }{2}\left( {2k_{0} - R} \right)^{2} + 2R\left( {2k_{0} - R} \right) + \frac{{R^{2} \pi }}{4}} \right) \hfill \\ \end{aligned} \right)N\left( t \right)^{2} $$\end{document}U=2k02LEJ1+L2EAS+Rπ8EAS+Rπ8GAS+R2EJ2π22k0-R2+2R2k0-R+R2π4Nt2 Then, the linear elastic deformation \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Delta L $$\end{document}ΔL of the round link can be expressed as9\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Delta L = \frac{\partial U}{{\partial N}} = \Delta L_{1} + \Delta L_{2} + \Delta L{}_{3} $$\end{document}ΔL=∂U∂N=ΔL1+ΔL2+ΔL3 where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Delta L_{1} $$\end{document}ΔL1,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Delta L_{2} $$\end{document}ΔL2,and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \Delta L{}_{3} $$\end{document}ΔL3 are defined as10\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left\{ \begin{gathered} \Delta L_{1} = \left( {\frac{{2k_{0}^{2} L}}{{EJ{}_{1}}} + \frac{L}{{2EA_{S} }}} \right)F_{k} \hfill \\ \Delta L_{2} = \left( {\frac{R\pi }{{8EA_{S} }} + \frac{R\pi }{{8GA_{S} }}} \right)F_{k} \hfill \\ \Delta L{}_{3} = \frac{R}{{2EJ_{2} }}\left( {\frac{\pi }{2}\left( {2k_{0} - R} \right)^{2} + 2R\left( {2k_{0} - R} \right) + \frac{{R^{2} \pi }}{4}} \right)F_{k} \hfill \\ \end{gathered} \right. $$\end{document}ΔL1=2k02LEJ1+L2EASFkΔL2=Rπ8EAS+Rπ8GASFkΔL3=R2EJ2π22k0-R2+2R2k0-R+R2π4Fk The relationship between the linear elastic deformation and the equivalent elastic force is therefore established. For the nonlinear elastic deformation due to contact, the Hertzian contact theory is used to derive the relationship between the indentation and the equivalent elastic force. As show in Fig. 4, the local indentation is induced between the adjacent round links subject to heavy axial loads. The principal radius of curvatures where the horizontal and vertical loops come into contact are defined as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \rho_{1} $$\end{document}ρ1 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \rho_{1}^{\prime } $$\end{document}ρ1′, respectively. The other principal radius of curvatures are defined as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \rho_{2} $$\end{document}ρ2 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \rho_{2}^{\prime } $$\end{document}ρ2′, respectively.Figure 4Diagram of local indentation deformations in the round link chain. As the vertical and horizontal loops have the same geometry, the sum of the curvatures \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \rho $$\end{document}ρ follows that^12^:11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left\{ \begin{gathered} \rho = \frac{1}{2}\left( {\frac{1}{{\rho_{1} }} + \frac{1}{{\rho_{2} }} + \frac{1}{{\rho_{1}^{\prime } }} + \frac{1}{{\rho_{2}^{\prime } }}} \right) = \frac{4R}{{d\left( {2R + d} \right)}} \hfill \\ \rho_{1} = \rho_{2} = \frac{d}{2} \hfill \\ \rho_{1}^{\prime } = \rho_{2}^{\prime } = - \left( {R + \frac{d}{2}} \right) \hfill \\ \end{gathered} \right. $$\end{document}ρ=121ρ1+1ρ2+1ρ1′+1ρ2′=4Rd2R+dρ1=ρ2=d2ρ1′=ρ2′=-R+d2 The local indentation is assumed as circles, and the radius \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ e $$\end{document}e is much smaller than that of the curvature at the contact point. Hence, according to Hertzian contact theory, the radius \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ e $$\end{document}e and the max stress in the contact area \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sigma_{\max } $$\end{document}σmax can be given as12\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left\{ \begin{gathered} e = \left( {\frac{{3\pi F_{k} \left( {k_{1} + k_{2} } \right)}}{4\rho }} \right)^{\frac{1}{3}} \hfill \\ \sigma_{\max } = \frac{{3F_{k} }}{{2\pi e^{2} }} \hfill \\ \end{gathered} \right. $$\end{document}e=3πFkk1+k24ρ13σmax=3Fk2πe2 where the parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ k_{1} = k_{2} = {{\left( {1 - \mu^{2} } \right)} \mathord{\left/ {\vphantom {{\left( {1 - \mu^{2} } \right)} {\left( {\pi E} \right)}}} \right. \kern-\nulldelimiterspace} {\left( {\pi E} \right)}} $$\end{document}k1=k2=1-μ2/πE. The stress in the area of contact is given as 13\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sigma_{{\left( {x,y} \right)}} = \sigma_{\max } \sqrt {1 - \frac{{x^{2} }}{{e^{2} }} - \frac{{y^{2} }}{{e^{2} }}} $$\end{document}σx,y=σmax1-x2e2-y2e2 The local deformations \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \delta_{1} $$\end{document}δ1 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \delta_{2} $$\end{document}δ2 at the contact area between the horizontal and vertical loops, respectively, are equal on the basis of symmetry. According to the geometrical relationship, the total contact-induced indentation deformation \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \delta $$\end{document}δ is obtained as14\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \delta = \delta_{1} + \delta_{2} = \left( {k_{1} + k_{2} } \right)\iint {\frac{{\sigma_{{\left( {x,y} \right)}} }}{r}ds} = \sqrt[3]{{\frac{{9RF_{k}^{2} \left( {1 - \mu^{2} } \right)^{2} }}{{E^{2} d\left( {2R + d} \right)}}}} $$\end{document}δ=δ1+δ2=k1+k2∬σx,yrds=9RFk21-μ22E2d2R+d3 Assuming the displacement at point A of the round link is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x_{A} $$\end{document}xA, which is a combination of the linear elastic deformation and the contacting indentation, the total displacement of the round link chain with n links \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x $$\end{document}x can be derived as15\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x = nx_{A} = 2n\left( {\Delta L_{1} + \Delta L_{2} + \Delta L_{3} } \right) + \left( {n - 1} \right)\delta $$\end{document}x=nxA=2nΔL1+ΔL2+ΔL3+n-1δ The equivalent nonlinear elastic force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{k} $$\end{document}Fk can then be expressed as a function of the displacement \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x $$\end{document}x. The equivalent nonlinear stiffness can also be obtained as16\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ k_{e} = \left[ {\frac{\partial x}{{\partial F_{k} }}} \right]^{ - 1} $$\end{document}ke=∂x∂Fk-1 ### Experimental identification of the nonlinear damping According to Eq. (1), the equivalent nonlinear damping for the round link chain can be expressed as17\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ c_{e} = \frac{{F - F_{m} - F_{k} }}{{\dot{x}}} $$\end{document}ce=F-Fm-Fkx˙ Since the inertial force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{m} $$\end{document}Fm and equivalent nonlinear elastic force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{k} $$\end{document}Fk have been expanded with respect to the displacement \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x $$\end{document}x, the equivalent nonlinear damping \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ c_{e} $$\end{document}ce can be formulated, once the external force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F $$\end{document}F and displacement \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x $$\end{document}x are tested in the experiment. This paper is focused on the round link chain in the scraper conveyer. In the practical working state, the fluctuating excitation is primarily induced by rotation of the driving system at some speed and initial deformations of the chain can also be observed due to transported coals. Thus, the external excitation can be taken as a simple harmonic motion, defined as18\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x = x_{0} + A\sin \left( {\omega t} \right) $$\end{document}x=x0+Asinωt where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ A $$\end{document}A is the amplitude of the fluctuating displacement, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \omega $$\end{document}ω is the angular velocity, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x_{0} $$\end{document}x0 is the initial deformation due to the pretension force. Thus, the output force corresponding to the excitation can be defined as19\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F = F_{0} + F_{A} \sin \left( {\omega t + \varphi } \right) $$\end{document}F=F0+FAsinωt+φ where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{0} $$\end{document}F0 is the pretension force. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{A} $$\end{document}FA is the amplitude of the pulsatile force response, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varphi $$\end{document}φ is the phase variation between the force and the displacement. The nonlinear dissipated force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{c} $$\end{document}Fc and corresponding equivalent nonlinear damping \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ c_{e} $$\end{document}ce can be obtained as20\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left\{ \begin{gathered} F_{e} = F_{0} + F_{A} \sin \left( {\omega t + \varphi } \right) + m\omega^{2} A\sin \left( {\omega t} \right) - F_{k} \hfill \\ c_{e} = \frac{{F_{c} }}{{\omega A\cos \left( {\omega t} \right)}} \hfill \\ \end{gathered} \right. $$\end{document}Fe=F0+FAsinωt+φ+mω2Asinωt-Fkce=FcωAcosωt It can be found from Eq. (20) the equivalent nonlinear damping is a function of exciting frequency \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \omega $$\end{document}ω and the amplitude of the fluctuating displacement \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ A $$\end{document}A. In this paper, the polynomial approximation method is utilized to approximate the equivalent nonlinear damping \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ c_{e} $$\end{document}ce. Since the highest power of both \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \omega $$\end{document}ω and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ A $$\end{document}A taken in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ c_{e} $$\end{document}ce is 2, as shown in Eq. (20), \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ c_{e} $$\end{document}ce is expanded as21\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ c_{e} = \alpha_{1} \omega^{2} + \alpha_{2} \omega + \alpha_{3} A^{2} + \alpha_{4} A + \alpha_{5} $$\end{document}ce=α1ω2+α2ω+α3A2+α4A+α5 where the coefficients \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \alpha_{1} $$\end{document}α1 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \alpha_{2} $$\end{document}α2 are connected to the square of the velocity. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \alpha_{3} $$\end{document}α3 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \alpha_{4} $$\end{document}α4 are connected to the velocity. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \alpha_{5} $$\end{document}α5 corresponds to the energy dissipation due to the pretension force. With five or more different experimental cases, all the coefficients can be determined and the equivalent nonlinear damping is then obtained. In order to improve the accuracy of the identified results, more than five experimental cases are used and the least square method for multiple variables are employed. Assuming there are *M*(M > 5) pairs of efficient test data and the fitted damping is \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \overline{c}_{e} $$\end{document}c¯e, then the least square sum of the deviations \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varphi $$\end{document}φ can be expressed as22\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \begin{aligned} \varphi & = \sum\limits_{i = 1}^{M} {\left( {c_{e} - \overline{c}_{e} } \right)^{2} } \\ & = \sum\limits_{i = 1}^{M} {\left( {c_{ei} - \alpha_{1} x_{1i} - \alpha_{2} x_{2i} - \alpha_{3} x_{3i} - \alpha_{4} x_{4i} - \alpha_{5} } \right)^{2} } \\ \end{aligned} $$\end{document}φ=∑i=1Mce-c¯e2=∑i=1Mcei-α1x1i-α2x2i-α3x3i-α4x4i-α52 where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x_{1i} $$\end{document}x1i,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x_{2i} $$\end{document}x2i,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x_{3i} $$\end{document}x3i and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x_{4i} $$\end{document}x4i are \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \omega_{i}^{2} $$\end{document}ωi2,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \omega_{i} $$\end{document}ωi,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ A_{i}^{2} $$\end{document}Ai2 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ A_{i} $$\end{document}Ai, respectively. To ensure \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varphi $$\end{document}φ is the minimal, the partial derivatives of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \varphi $$\end{document}φ with respect to the unknown coefficients are taken to be zero. The corresponding system of linear equations in a simple form is obtained as23\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \bf GA = C $$\end{document}GA=C where the elements in the matrix \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ G $$\end{document}G, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ C $$\end{document}C are defined as24\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left\{ \begin{gathered} g_{kq} = \sum\limits_{i = 1}^{M} {x_{ki} x_{qi} } \hfill \\ c_{q} = \sum\limits_{i = 1}^{M} {c_{ei} x_{qi} } \hfill \\ \end{gathered} \right. $$\end{document}gkq=∑i=1Mxkixqicq=∑i=1Mceixqi where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x_{5i} = 1 $$\end{document}x5i=1. The coefficients for the nonlinear damping are then calculated. ## Experimental setup and scheme In order to determine the coefficients for the equivalent nonlinear damping, an electronic universal tensile testing machine was employed to perform the dynamic tensile tests on the 14 × 50 round link chain used in mining applications. The circular chain is extracted from the coal mine of Shandong Energy Group. Figure 5 is the actual working condition of the circular chain used in the coal mine scraper conveyor. The schematic diagram of the experiment system is illustrated in Fig. 6. Harmonic displacement signals are generated by a signal generator controlled by a laptop. The signal is magnified by a power amplifier and drives the exciter to stretch the round link chain. In order to simulate practical loading conditions, a preload is firstly applied to ensure the round link chain is tensile condition in the test. The initial displacement, amplitude and frequency of the fluctuating load are set in the laptop. The total tensile force and displacement are measured by force sensor and displacement sensor, respectively. The measured signals are post processed by the data acquisition equipment and analyzed in the laptop.Figure 5Actual use of underground ring chain in coal mine.Figure 6The schematic diagram of the test system. The round link chain in the test is made of 25MnSi steel (*E* = 210 Gpa) and the mass per unit length is *m* = 4.0 kg/m. The geometrical dimensions of the round link chain used in the tests are shown in Fig. 7. The number of the round links between the fixtures is 16 (*n* = 16), and the effective length of the chain between the fixtures was 800 mm. The ambient temperature during the test is 23 °C.Figure 7Dimensions of the round link. ## Results and discussion ### The equivalent nonlinear stiffness In this part, the proposed theoretical model is employed to formulate the nonlinear stiffness of the 14 × 50 round link chain in the test. As shown in Fig. 8a, the stiffness of the round link chain is nonlinear with respect to the displacement *x*, especially when *x* is in the range from 0 to 3 mm. That is, when the axial deformation of each link is 0–0.3 mm, the nonlinearity of the stiffness is the strongest. With increase in the length of the chain, the stiffness and the nonlinearity are both weakened. When the displacement is large enough, the stiffness of the chain is inversely proportional to the length, which is the case for the linear elastic problems.Figure 8Parameter effects on the nonlinear stiffness, (**a**) effects of the total length and (**b**) effects of the elastic module with *E*~0~ = 2.1 × 10^5^ Mpa and *n* = 16. The influences of the material parameters on the stiffness of the chain are shown in Fig. 8b. A rise in the elastic module seems to reinforce the stiffness of the chain. Nevertheless, the stiffness is not any more proportional to the elastic module. The stiffness tends to be linear, when the material is more elastic. ### The equivalent nonlinear damping In order to obtain the coefficients for the equivalent nonlinear damping, a series of experiments are performed, following the setup described in part 3. One of the tests is presented in detail to demonstrate the detailed identification procedure for the coefficients. In the case, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x_{0} $$\end{document}x0, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ A_{0} $$\end{document}A0, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \omega $$\end{document}ω and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ k_{0} $$\end{document}k0 are taken as 2.2 mm, 0.25 mm, 1.884 rad/s and 0.08973 mm/s, respectively. Thus, the displacement excitation is expressed as25\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x = \left\{ \begin{gathered} 0.08973t{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left( {0 \le t \le t_{0} } \right) \hfill \\ 2.2 + 0.25\sin \left( {1.884t} \right){\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left( {t_{0} < t} \right) \hfill \\ \end{gathered} \right.{\kern 1pt} {\kern 1pt} {\kern 1pt} $$\end{document}x=0.08973t0≤t≤t02.2+0.25sin1.884tt0<t Figure 9a shows time-varying curves for the total displacement of the round link chain, and the dynamic force responses are measured with the force senor, as shown in Fig. 9b. With increment in the deformation of the round link chain, the nonlinearity is weakened. Overall, the frequency of the force responses keeps in agreement with that of the exciting displacement. This phenomenon and the variation of the stiffness to the displacement coincide mutually.Figure 9The loading and responses in time domain. The relationship between the force responses and the displacement excitation in the steady state is obtained, as shown in Fig. 10. The approach track is distinct from the restitution one, which implies that significate energy dissipation is observed in the vibration of the round link chain. Therefore, it is necessary to involve the damping property into the dynamic model of the round chain under axial excitation, especially if the load is relatively large, like the loads in the mining engineering.Figure 10Force–displacement curve. According to the measured data for the dynamic force response, a fitting function in waveform is employed to derive the formula of the force with respect to time *t*, expressed as26\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \begin{aligned} F & = F_{0} + F_{A} \sin \left( {\omega t + \phi } \right) \\ & = 11117.21 + 19522.1\;\sin \left( {1.883t - 0.8350} \right) \\ \end{aligned} $$\end{document}F=F0+FAsinωt+ϕ=11117.21+19522.1sin1.883t-0.8350 It is worth noting that R-square for the fitting is equal to 0.9999 and the fitting is therefore accurate enough. Now that the displacement, the geometrical and material parameters for the tested round link chain are determined, then the equivalent elastic force of the chain can also yield based on Eq. (15). One can calculate the elastic force and then submit it to Eq. (20) at each time step to determine the damping. Alternatively, to simplify the calculation, the formulation of the elastic force can also be fitted in the similar way. According to Eq. (20), the dissipated force and corresponding damping coefficient can be derived by using the property of the trigonometric functions. In this case, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{c} $$\end{document}Fc and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ c_{e} $$\end{document}ce can be expressed as27\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left\{ \begin{gathered} F_{c} = F_{\alpha c} \cos \left( {\omega t} \right) = 1335\cos \left( {1.884t} \right) \hfill \\ c_{e} = \frac{{F_{C} }}{{\omega A\cos \left( {\omega t} \right)}} = {{2834.4\;{\text{N}}\;{\text{s}}} \mathord{\left/ {\vphantom {{2834.4\;{\text{N}}\;{\text{s}}} {{\text{mm}}}}} \right. \kern-\nulldelimiterspace} {{\text{mm}}}} \hfill \\ \end{gathered} \right. $$\end{document}Fc=Fαccosωt=1335cos1.884tce=FCωAcosωt=2834.4Ns/mm As expressed in Eq. (21), the nonlinear damping depends on the amplitude and frequency of the external excitation source. A series of experiments are needed to get the coefficients in Eq. (21). Hence, with the same experimental model (the 14 × 50 round link chain), 9 different tests were implemented based on the aforementioned experimental setup. Corresponding damping properties are determined with the similar framework presented ahead, as listed in Table 1.Table 1Damping of the chain under different loading frequencies.Test number\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \omega $$\end{document}ω/rad/s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ X_{0} $$\end{document}X0/mm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ A $$\end{document}A/mm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{ac} $$\end{document}Fac/N\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ c_{e} $$\end{document}ce/(N s/mm)112.20.2580005095.5242.20.2511003503.2362.20.2513352834.4482.20.2514652332.85102.20.2515421964.3662.150.29852614.1762.250.316432906.9862.30.3520173058.8962.350.423293090.5 A two-variable second-order polynomial equation is defined to fit the equivalent damping of the round link chain, as shown in Eq. (21). According to Eq. (23), the data in Table. 1 is utilized to calculate the coefficients in Eq. (21). The coefficients are finally obtained as28\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left[ {\alpha_{1} {\kern 1pt} {\kern 1pt} \alpha_{2} {\kern 1pt} {\kern 1pt} \alpha_{3} {\kern 1pt} {\kern 1pt} \alpha_{4} {\kern 1pt} {\kern 1pt} \alpha_{5} } \right] = \left[ {28.02{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} - 651.57{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} - 8386.50{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} 7394.84{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} 4378.89} \right] $$\end{document}α1α2α3α4α5=28.02-651.57-8386.507394.844378.89 As shown in Fig. 11, the damping is significant at low frequencies, while increase in frequency leads to decrease in the damping. However, when the frequency rises over about 12.5 Hz, the damping increases rapidly and it tends to be more significant at high frequencies. It may be because that the energy dissipation mainly lie in the plastic deformation of the round link in the low frequency range, while at high frequencies, friction and pounding play a major role in dissipating energy.Figure 11The nonlinear damping properties with respect to vibration amplitude and frequency. ### Validation of the prediction model In order to validate the proposed framework for determination of the dynamic parameters of the round link chain, another test is implemented. Totally different key parameters are taken for the simulated harmonic displacement load to make sure the validation is accurate. The initial displacement \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x_{0} = 2.25\;{\text{mm}} $$\end{document}x0=2.25mm, the loading frequency \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ f = 7\;{\text{Hz}} $$\end{document}f=7Hz, and the amplitude \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ A = 0.32\;{\text{mm}} $$\end{document}A=0.32mm. Then displacement excitation turns to be \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ x = 2.25{ + 0}{\text{.32}}\sin \left( {14\pi t} \right) $$\end{document}x=2.25+0.32sin14πt. According to Eq. (15), the elastic force \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ F_{k} $$\end{document}Fk is obtained, and the damping can be given by submitting \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ A = 0.32\;{\text{mm}} $$\end{document}A=0.32mm and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ f = 7\;{\text{Hz}} $$\end{document}f=7Hz into Eq. (21). According to Eq. (1), the force response is then determined. The curves for the force response to the exciting displacement are employed to analyze the accuracy of the proposed method. The proposed result is compared with the experimentally measured curve, as shown in Fig. 12. Good agreement is observed between the proposed and experimental results.Figure 12Comparison between the measured results and the prediction curve. ## Conclusion This paper presents a new framework to identify the nonlinear stiffness and damping properties of a round link chain under axial dynamic excitation. A modified nonlinear viscoelastic model is developed to study vibrations of the chain with numbers of round links. Based on linear elastic theory and contacting models, a theoretical model is established for the equivalent nonlinear stiffness, including elastic and contacting deformation effects. A series of dynamic experiments are performed and accordingly, the least square fitting method is employed to determine the equivalent damping of the chain. Another test case with totally different loading parameters is designed to validate the proposed framework and the proposed results agree well with the experimental ones. The effects of key material and loading parameters on the nonlinear stiffness and damping are examined. It is found with small dynamic deformations the stiffness is strong nonlinear to the deformation, and the nonlinearity will turn to small if the deformation or the length of the chain increase, or the elasticity of the chain decrease. Due to transformation of the main manners of the energy dissipation, the nonlinear damping is firstly decreased with rise in the loading frequency. However, when the frequency increase over a critical one, the damping turns to increase. That is, the nonlinear damping is significant at both low and high frequencies, while it is relatively small in the middle frequency range. The developed identification framework and proposed results can provide some physical insights into the dynamic properties of the round chain and can help improve the performance of the machines such as the scraper conveyor.