Nonlinear parametric vibration analysis of radial gates

Abstract: In this study, based on the dynamic characteristics of radial gates, a nonlinear differential equation of motion for a radial gate arm is derived. The excitation conditions of both principal and subharmonic parametric resonance cases are obtained by using multi-scale and numerical methods. The parameter analysis reveals the following: In comparison to the analytical calculation method (e.g., Multi-scale) of dynamic instability region division, this method effectively analyzes the parametric vibration of radial gates in terms of the end moment, vibration duration, and amplitude. The vibration amplitudes of the arm increase with an increase in its length and excitation amplitude as well as a decrease in the arm inclined angle. Moreover, the parametric resonance excitation is facilitated, and the resonance region becomes larger with an increase in the initial end moment. Comparing with damping effect, vibration response of the arm is mainly affected by the nonlinear behavior in system. Therefore, the energy transfer method (e.g. Add tuned mass dampers) should be taken into account in vibration control.

Keywords: Radial gate; parametric vibration; multi-scale method; parameter analysis

1. Introduction

A radial gate is a main regulating structure employed in hydraulic engineering, and its stability is directly related to the safe operation of a reservoir. However, the dynamic instability mechanism of radial gates has not been elucidated thus far, and as a result, the existing structures that are constructed without appropriate anti-vibration design are vulnerable. Even if such structures are designed for anti-vibration performance, their design may not comply with current design code [1] and thus may exhibit deficient anti-vibration performance. In recent years, several engineering practices, such as Wachi Dam (1967) in Japan and the Folsom Dam (1995) in the United States, have shown that radial gates will produce different degrees of vibration and even damage structure considerably to flow-induced excitation when the gates were partially opened [2-3]. The lightness, slenderness, and flexibility of the arm render the structure to vibration and dynamic instability easily, when subjected to water loads. Moreover, the increase of water level will reduce the structural frequency, this characteristic indicates that the radial gate should be highly vulnerable under high water head [4]. Consequently, the dynamic stability of the radial gate arm has become a crucial issue with many influence factors that needs to be resolved in both the design and operation stages [5]. According to excitations such as cavity flow, Kármán vortices and special hydrodynamic load, the vibration patterns of the gates can be categorized as follows: forced vibration, self-excited vibration, and parametric vibration. Based on the elongated nature of radial gate arms as mentioned above, among the various types of gates, the parametric vibration only appears in radial gates. The problems of parametric vibration are encountered in many branches of physics and engineering, have been studied extensively, such as in cable structures [6], beam structures [7], plate and shell structures [8], and coupled models of multi-member interaction [9]. In 1980s, Zhang, J.G. [10] observed the phenomenon of parametric vibration during gate operation in dynamic water, considered it as one of the important causes of dynamic instability in the arm through engineering failure examples, and proposed a concise analytical calculation method. Subsequently, there are thought-provoking researches that have been carried out on the division of the dynamic instability region [11]. Although this region can be used to determine whether parametric resonance occurs or not, it does not yield the response amplitude at an arbitrary time. In general, the radial gate is subjected to flow-induced load for a relatively short period of time, and the arm amplitudes continue to fluctuate within a narrow range without the magnification factor tending to infinity. Hence, the method of dynamic instability region may not be an accurate criterion for defining the dynamic instability without considering the time factor. Recently, Zhang, J. and Xie, Z.X. [12] investigated the unstable behavior of nonlinear systems in the radial gate arm under certain dynamic loads and proved the necessity to analyze the stability of radial gate in terms of the excitation amplitude and time duration.

In this study, according to the mechanical characteristics of the radial gate, with by completely considering the geometric nonlinear factors caused by the bending moment and large displacement at the high head, the nonlinear differential equation of motion for an arm is derived. The stability of parametric vibration is obtained using multi-scale and numerical methods. Moreover, the effects of excitation factors (water head, frequency, time, and amplitude) and structural factors (inclined angle, length, and damping) on the vibration response of the arm are quantitatively analyzed.

2. Derivation of the nonlinear differential equation of motion

Three dimensional overview of the radial gate is shown in Figure 1, the gate contains a face plate to retain water, two major crossbeams, four inclined arms with box-sectional area, major vertical beams and minor beams [13]. The anchor of radial gate arm is attached to piers, which can be regard as a semi-rigid connection. In accordance with Chinese design code SL74-2013 [1], the structure can be simplified as an extended beam, both ends constrained by water seal in y direction, and assembled with arms by welding. The water pressure can pass to piers through face plane, main beams and arms, while, the load acting on main beam can be considered as uniform hydrostatic pressure q by equivalent calculation (see Figure 2). Then on the base of these, further analyze the arm vibration under local coordinate, the model as shown in Figure 3, with the length of the arm L, the transverse displacement v(x,t), flow-induced excitation frequency ω_d, the tangential hydrostatic force T₀= q (L_b + 2c)/ 2, and the additional tangential hydrodynamic force τ. Due to parametric vibration is induced by harmonic load, so τ can be assumed as the function of harmonic displacement of flow excitation D_d = A_dcosω_dt.

Figure 1. Three dimensional overview of the radial gate.	Figure 2. The main frame of in-plane radial gate.

The following basic hypotheses are made before establishing the equation of motion for radial gate arms:

(1) Material properties: based on previous hydraulic experiments on full-scale or model gates, the arm, which are considered homogeneous and isotropic, and always behave elastically under design loads, even fails in a predominantly elastic mode (sway-type failures). For convenience and safety, most study of hydraulic experiments, suppose the whole structure work in elastic state.

(2) Element type: due to the ends of arm are restricted by major crossbeams and piers, respectively, and its length is much larger than dimensions of cross-section, so the arm can be regarded as a Bernoulli-Euler beam, the torsion and shear effect are negligible.

Recently, Zhang, X.C. et al [14] proposed a non-moment design method for eliminating end bending moment of vertical main frame (x-z plane). This method reduces out-of-plane eccentric compression on the arm and improve the overall stability bearing capacity. But the study of in-plane (x-y plane) layout scheme about end bending moment has not been reported thus far, while Hu, M.S. et al. [15] examined that the effect of in-plane initial end bending moment M₀ cannot be neglected from prototype observation. According to the above researches, in-plane vibration of arm is discussed in this work. The relation between the deflection and moment can be written as:

$EI{\eta }^{''}=-M_{0}$

(1)

where, according to design code (SL74-2013), in order to ensure the safety of the calculation results, the major crossbeam is assumed as a double cantilever beam and M₀ = qc²/2.

In order to keep derivations in this study clear, according to the relatively simple boundary conditions without loss of generality (Figure 2), The sag of the arm is expressed as follows:

$\eta ={\frac {M_{0}x^{2}}{2EI}}{\mbox{+}}{\frac {M_{0}x}{R}}$

(2)

where EI denotes the flexural rigidity and R, the rotational stiffness of semi-rigid connection.

Figure 3. Vibration model of a radial gate arm.

The equation of motion for the arm is as follows:

${\frac {\partial }{\partial s}}\left[\left(T_{0}+\tau \right)\left({\frac {{\mbox{d}}\eta }{{\mbox{d}}s}}+\right.\right.$$\left.\left.{\frac {{\mbox{d}}v}{{\mbox{d}}s}}\right)\right]{\mbox{=}}m{\frac {{\partial }^{2}v}{\partial t^{2}}}+$$EI\left({\frac {{\partial }^{4}\eta }{\partial s_{}^{4}}}+\right.$$\left.{\frac {{\partial }^{4}v}{\partial s_{}^{4}}}\right)+$$\mu {\frac {\partial v}{\partial t}}$

(3)

where s represents the lagrange coordinate by the arc length, m is the mass per unit length of the arm and the damping coefficient μ. When α = 90°, the arm is upright, same as layout of Figure 2.

According to the geometric deformation and stress state of the arm model, the following results are obtained:

${\frac {\tau }{h}}={\frac {T_{0}}{H_{0}}}={\frac {{\mbox{d}}s}{{\mbox{d}}x}}=$${\left[1+{\left({\frac {{\mbox{d}}\eta }{{\mbox{d}}x}}\right)}^{2}\right]}^{\frac {1}{2}}$

(4)

where H₀ and h represent the axial force and additional axial hydrodynamic force, respectively.

By Taylor formula, Equation (5) is expanded at dη/dx = 0:

${\left[1+{\left({\frac {{\mbox{d}}\eta }{{\mbox{d}}x}}\right)}^{2}\right]}^{\frac {1}{2}}=$$1+{\frac {1}{2}}{\left({\frac {M_{0}x}{EI}}{\mbox{+}}{\frac {M_{0}}{R}}\right)}^{2}$

(5)

According to static force equilibrium equation of radial gate arm:

$H_{0}{\frac {{\mbox{d}}^{2}\eta }{{\mbox{d}}s_{}^{2}}}{\mbox{=}}EI{\frac {{\partial }^{4}\eta }{\partial s_{}^{4}}}$

(6)

Substituting Equations (4), (5) and (6) into Equation (3), the following equation is obtained:

$\left[H_{0}{\frac {{\mbox{d}}^{2}v}{{\mbox{d}}s_{}^{2}}}{\mbox{+}}h\left({\frac {{\mbox{d}}^{2}\eta }{{\mbox{d}}s_{}^{2}}}+\right.\right.$$\left.\left.{\frac {{\mbox{d}}^{2}v}{{\mbox{d}}s_{}^{2}}}\right)\right]\left[1+\right.$$\left.{\frac {1}{2}}{\left({\frac {M_{0}x}{EI}}{\mbox{+}}{\frac {M_{0}}{R}}\right)}^{2}\right]{\mbox{=}}m{\frac {{\partial }^{2}v}{\partial t^{2}}}+$$EI{\frac {{\partial }^{4}v}{\partial s_{}^{4}}}+\mu {\frac {\partial v}{\partial t}}$

(7)

The length variation of micro-segment ds' is assumed as follows:

${\left({\mbox{d}}{s^{'}}_{}\right)}^{2}={\left({\mbox{d}}x+{\frac {\partial u}{\partial x}}{\mbox{d}}x\right)}^{2}+$${\left({\mbox{d}}\eta +{\frac {\partial u}{\partial \eta }}{\mbox{d}}\eta \right)}^{2}$

(8)

where u (x, t) is the longitudinal displacement.

Substituting Equation (8) into Equation (4) and neglecting the quadratic term of longitudinal strain gives the following equation:

$\epsilon ={\frac {{\mbox{d}}s^{'}-{\mbox{d}}s}{{\mbox{d}}s}}\approx {\frac {{\left({\mbox{d}}s^{'}\right)}^{2}-{\left({\mbox{d}}s\right)}^{2}}{2{\left({\mbox{d}}s\right)}^{2}}}{\mbox{=}}{\frac {\partial u}{\partial s}}{\frac {{\mbox{d}}x}{{\mbox{d}}s}}+$${\frac {\partial v}{\partial s}}{\frac {{\mbox{d}}\eta }{{\mbox{d}}s}}+$${\frac {1}{2}}{\left({\frac {\partial v}{\partial s}}\right)}^{2}$

(9)

As ds ≈ dx, Equation (9) can be rewritten as follows:

$\epsilon ={\frac {\partial u}{\partial x}}+{\frac {\partial v}{\partial x}}{\frac {{\mbox{d}}\eta }{{\mbox{d}}x}}+$${\frac {1}{2}}{\left({\frac {\partial v}{\partial x}}\right)}^{2}$

(10)

The boundary conditions for generalized displacements vector (Figure 2) are written as:

$v\left(0,t\right)=0$ , $v\left(L,t\right)=0$ , $u\left(0,t\right)=D_{\mbox{d}}cos\alpha$ , $u\left(L,t\right)=0$ .

(11)

Satisfying the boundary conditions (11), by separation of variables, the vibration mode of the arm can be expressed as:

$v\left(x,t\right)=V\left(t\right)\sum _{n=1}^{N}sin{\frac {n\pi x}{L}}$

(12)

where V(t) is the displacement amplitude. According to early analysis, Liu, J.L. et al [11] pointed out that, for radial gate arm vibration induced by flow excitation, the first mode is dominated. Hence, We can use the single mode approximation with neglecting modal interaction.

The axial hydrodynamic force of radial gate arm can be expressed as follows:

$h={\frac {EA}{L}}{\int }_{0}^{L}\epsilon {\mbox{d}}x={\frac {EA}{L}}\left(u_{\mbox{d}}{\mbox{+}}{\frac {2L^{2}M_{0}V}{{\pi }^{3}EI}}+\right.$$\left.{\frac {{\pi }^{2}V_{}^{2}}{4L}}\right)$

(13)

where $u_{\mbox{d}}=U_{\mbox{d}}cos\left({\omega }_{\mbox{d}}t\right)=$$D_{\mbox{d}}cos\left(\alpha \right)$ .

Substituting Equations (2), (12), and (13) into Equation (7) yields the following equation:

${\begin{array}{c}\left[H_{0}\left(-V{\frac {{\pi }^{2}}{L_{}^{2}}}sin{\frac {\pi x}{L}}\right)+{\frac {EA}{L}}\left(u_{\mbox{d}}{\mbox{+}}{\frac {2L^{2}M_{0}V}{{\pi }^{3}EI}}+{\frac {{\pi }^{2}V_{}^{2}}{4L}}\right)\left(-{\frac {M_{0}}{EI}}-V{\frac {{\pi }^{2}}{L_{}^{2}}}sin{\frac {\pi x}{L}}\right)\right]\left[1+{\frac {1}{2}}{\left({\frac {M_{0}x}{EI}}{\mbox{+}}{\frac {M_{0}}{R}}\right)}^{2}\right]\\=m{\ddot {V}}sin{\frac {\pi x}{L}}+EIV{\frac {{\pi }^{4}}{L_{}^{4}}}sin{\frac {\pi x}{L}}+\mu {\dot {V}}sin{\frac {\pi x}{L}}\end{array}}$

(14)

Then, according to the Galerkin method, the equation of motion is obtained as follows:

${\ddot {V}}+2\xi {\omega }_{\mbox{a}}{\dot {V}}+\left({\omega }_{\mbox{a}}^{2}+\right.$$\left.au_{\mbox{d}}\right)V+bV_{}^{2}+cV_{}^{3}+du_{\mbox{d}}=$$0$

(15)

where

$2\xi {\omega }_{\mbox{a}}={\frac {\mu }{m}}$ , ${\omega }_{\mbox{a}}^{2}={\omega }_{0}^{2}\left(1+\lambda \right)$ , ${\omega }_{0}^{2}={\frac {{\pi }^{2}H_{0}}{mL_{}^{2}}}+{\frac {EI{\pi }^{4}}{mL_{}^{4}}}$ ,

$\lambda =\left[-{\frac {{\pi }^{2}H_{0}}{2L}}\left(1{\mbox{+}}{\frac {M_{0}^{2}}{2R^{2}}}\right)-\right.$$\left.{\frac {{\pi }^{2}H_{0}M_{0}^{2}}{2E^{2}I^{2}L^{2}}}\left({\frac {L^{3}}{6}}-\right.\right.$$\left.\left.{\frac {L^{3}}{4{\pi }^{2}}}\right)-{\frac {M_{0}^{2}H_{0}{\pi }^{2}}{4REI}}-\right.$$\left.{\frac {2L^{3}M_{0}^{2}A}{{\pi }^{4}EI^{2}}}\right]/\left[\left({\frac {{\pi }^{2}H_{0}}{mL_{}^{2}}}+\right.\right.$$\left.\left.{\frac {EI{\pi }^{4}}{mL_{}^{4}}}\right)\left(-\right.\right.$$\left.\left.{\frac {mL}{2}}\right)\right]$ ,

$a=-{\frac {{\pi }^{2}EA}{L^{3}}}\left[{\frac {L}{2}}\left(1{\mbox{+}}{\frac {M_{0}^{2}}{2R^{2}}}\right)+\right.$$\left.{\frac {M_{0}^{2}}{2E^{2}I^{2}}}\left({\frac {L^{3}}{6}}-{\frac {L^{3}}{4{\pi }^{4}}}\right)+\right.$$\left.{\frac {M_{0}^{2}L^{2}}{4REI}}\right]/\left(-{\frac {mL}{2}}\right)$ ,

$b=\left\{{\begin{array}{c}\\\end{array}}\left({\frac {-2{\pi }^{2}M_{0}A}{{\pi }^{3}IL}}\right)\left[{\frac {L}{2}}\left(1{\mbox{+}}{\frac {M_{0}^{2}}{2R^{2}}}\right)+\right.\right.$$\left.\left.{\frac {M_{0}^{2}}{2E^{2}I^{2}}}\left({\frac {L^{3}}{6}}-\right.\right.\right.$$\left.\left.\left.{\frac {L^{3}}{4{\pi }^{2}}}\right)+\right.\right.$$\left.\left.{\frac {M_{0}^{2}L^{2}}{4REI}}\right]-{\frac {M_{0}{\pi }^{2}A}{4IL^{2}}}\left[{\frac {L^{2}}{\pi }}\left(1{\mbox{+}}{\frac {M_{0}^{2}}{2R^{2}}}\right)+\right.\right.$$\left.\left.{\frac {M_{0}^{2}}{2E^{2}I^{2}}}\left({\frac {L^{4}}{\pi }}-\right.\right.\right.$$\left.\left.\left.{\frac {6L^{4}}{{\pi }^{3}}}\right)+\right.\right.$$\left.\left.{\frac {M_{0}^{2}L^{2}}{REI\pi }}\right]\right\}/\left(-\right.$$\left.{\frac {mL}{2}}\right)$ ,

$c=\left({\frac {-EA{\pi }^{4}}{4L^{4}}}\right)\left[{\frac {L}{2}}\left(1{\mbox{+}}{\frac {M_{0}^{2}}{2R^{2}}}\right)+\right.$$\left.{\frac {M_{0}^{2}}{2E^{2}I^{2}}}\left({\frac {L^{3}}{6}}-{\frac {L^{3}}{4{\pi }^{2}}}\right)+\right.$$\left.{\frac {M_{0}^{2}L^{2}}{4REI}}\right]/\left(-{\frac {mL}{2}}\right)$ ,

$d=\left(-{\frac {AM_{0}}{LI}}\right)\left[{\frac {L^{2}}{\pi }}\left(1{\mbox{+}}{\frac {M_{0}^{2}}{2R^{2}}}\right)+\right.$$\left.{\frac {M_{0}^{2}}{2E^{2}I^{2}}}\left({\frac {L^{4}}{\pi }}-\right.\right.$$\left.\left.{\frac {6L^{4}}{{\pi }^{3}}}\right)+{\frac {M_{0}^{2}L^{2}}{REI\pi }}\right]/\left(-\right.$$\left.{\frac {mL}{2}}\right)$ .

where ξ is the damping ratio, ω₀ corresponds to the first-order natural frequency and λ is the arm parameter of the end moment.

In comparison to the classical parametric vibration models [2, 11-12], this model is highly comprehensive and it reflects the characteristics of the parametric vibration of radial gates with greater accuracy.

3. Stability analysis

In the present investigation, as both the stretching of the arm and the flow excitation are small, the nonlinear terms and the parametric excitation terms are smaller than the linear terms. Therefore, the method of multiple scales is applied to Equation (15) to asymptotically describe the slow dynamics of the model [16-17]. To weaken the nonlinear terms, the following substitutions are made:

${\epsilon }^{2}\chi =-2\xi {\omega }_{\mbox{a}}$ , ${\epsilon }^{2}f_{1}=-a$ , ${\epsilon }^{2}f_{2}=-d$

(16)

where ε is a small bookkeeping parameter. Then, Equation (15) can be rewritten as:

${\ddot {V}}+{\omega }_{\mbox{a}}^{2}V={\epsilon }^{2}\chi {\dot {V}}+$${\epsilon }^{2}f_{1}u_{\mbox{d}}V-bV_{}^{2}-cV_{}^{3}+{\epsilon }^{2}f_{2}u_{\mbox{d}}$

(17)

Supposing the solution to Equation (17) can be expressed as:

$V=\epsilon V_{11}\left(T_{0},T_{1},T_{2}\right)+{\epsilon }^{2}V_{12}\left(T_{0},T_{1},T_{2}\right)+$${\epsilon }^{3}V_{13}\left(T_{0},T_{1},T_{2}\right)$

(18)

where T₀ = t, T₁ = tε, T₂ = tε². Substituting Equation (18) into Equation (17) and letting the coefficients ε, ε² and ε³ equal zero, the following equations are obtained:

$D_{0}^{2}V_{11}+{\omega }_{\mbox{a}}^{2}V_{11}=0$

(19)

$D_{0}^{2}V_{12}+{\omega }_{1}^{2}V_{12}=-2D_{0}D_{1}V_{11}-bV_{11}^{2}+$$f_{2}U_{\mbox{d}}sin\left({\omega }_{\mbox{d}}t\right)$

(20)

$D_{0}^{2}V_{13}+{\omega }_{\mbox{a}}^{2}V_{13}=-2D_{0}D_{1}V_{12}-$$\left(D_{1}^{2}+2D_{0}D_{2}\right)V_{11}+\chi D_{0}V_{11}+f_{1}U_{\mbox{d}}sin\left({\omega }_{\mbox{d}}t\right)V_{11}-$$2bV_{11}V_{12}-cV_{11}^{3}$

(21)

The complex solution to Equation (19) can be assumed as:

$V_{11}=A\left(T_{1},T_{2}\right)e^{i{\omega }_{\mbox{a}}T_{0}}+$${\overline {A}}\left(T_{1},T_{2}\right)e^{-i{\omega }_{\mbox{a}}T_{0}}$

(22)

where ${\overline {A}}$ is the complex conjugate of A.

Substituting Equation (22) into Equation (20), one gets:

$D_{0}^{2}V_{12}+{\omega }_{a}^{2}V_{12}=-2i{\omega }_{a}D_{1}Ae^{i{\omega }_{\mbox{a}}T_{0}}-$$bA^{2}e^{i2{\omega }_{\mbox{a}}T_{0}}+i{\frac {f_{2}}{2}}U_{\mbox{d}}e^{i{\omega }_{\mbox{d}}T_{0}}+$$bA{\overline {A}}+cc$

(23)

where cc indicates the complex conjugate part. Eliminating the secular terms yields:

$V_{12}={\frac {bA^{2}}{2{\omega }_{\mbox{a}}^{2}}}e^{i2{\omega }_{\mbox{a}}T_{0}}-$${\frac {if_{2}U_{\mbox{d}}}{2\left({\omega }_{\mbox{d}}^{2}-{\omega }_{\mbox{a}}^{2}\right)}}e^{i{\omega }_{\mbox{d}}T_{0}}+$${\frac {bA{\overline {A}}}{{\omega }_{\mbox{a}}^{2}}}+cc$

(24)

Substituting Equations (22) and (24) into Equation (21) leads to:

${\begin{array}{c}D_{0}^{2}V_{13}+{\omega }_{\mbox{a}}^{2}V_{13}=\left(\chi i{\omega }_{\mbox{a}}A-2D_{2}i{\omega }_{\mbox{a}}A-{\frac {8b^{2}A^{2}{\overline {A}}}{3{\omega }_{\mbox{a}}^{2}}}-3cA^{2}{\overline {A}}-{\frac {bA}{{\omega }_{\mbox{a}}^{2}}}\right)e^{i{\omega }_{\mbox{a}}T_{0}}-\left({\frac {2b^{2}A^{3}}{3{\omega }_{\mbox{a}}^{2}}}+cA^{3}\right)e^{i3{\omega }_{\mbox{a}}T_{0}}\\{\begin{array}{ccccc}&&&&\end{array}}+\left({\frac {bif_{2}A}{{\omega }_{\mbox{d}}^{2}-{\omega }_{\mbox{a}}^{2}}}+{\frac {f_{1}}{2}}iA\right)U_{\mbox{d}}e^{i\left({\omega }_{\mbox{a}}+{\omega }_{\mbox{d}}\right)T_{0}}+\left({\frac {bif_{2}{\overline {A}}}{{\omega }_{\mbox{d}}^{2}-{\omega }_{\mbox{a}}^{2}}}+{\frac {f_{1}}{2}}i{\overline {A}}\right)U_{\mbox{d}}e^{i\left({\omega }_{\mbox{a}}+{\omega }_{\mbox{d}}\right)T_{0}}+cc\end{array}}$

(25)

Through practical observation [10], when excitation frequency ω_d approaches 2 times of the structure frequency ω_a, the parametric resonance may occur. A detuning parameter σ is introduced to quantify the deviation of ω_d from 2ω_a, and ω_d is described by:

${\omega }_{\mbox{d}}=2{\omega }_{\mbox{a}}+{\epsilon }^{2}\sigma$

(26)

Substituting Equation (26) into Equation (25) and eliminating secular terms yields：

$-2i{\omega }_{\mbox{a}}{\frac {\partial A}{\partial T_{2}}}+$$\chi i{\omega }_{\mbox{a}}A-{\frac {14b^{2}A^{2}{\overline {A}}}{3{\omega }_{\mbox{a}}^{2}}}-$$3cA^{2}{\overline {A}}+\left({\frac {bif_{2}{\overline {A}}}{{\omega }_{\mbox{d}}^{2}-{\omega }_{\mbox{a}}^{2}}}+\right.$$\left.{\frac {f_{1}}{2}}i{\overline {A}}\right)U_{\mbox{d}}e^{i\sigma T_{0}}=$$0$

(27)

Then, by expressing the function A in the polar form as:

$A={\frac {1}{2}}\alpha \left(t\right)e^{i\theta \left(t\right)}$

(28)

where α(t) and θ(t) are real functions. By substituting this transformation into Equation (27), and separating the resulting real and imaginary parts yield:

${\dot {\alpha }}={\frac {1}{2}}\chi \alpha +\left[{\frac {bf_{2}}{2\left({\omega }_{\mbox{d}}^{2}-{\omega }_{\mbox{a}}^{2}\right)}}+\right.$$\left.{\frac {f_{1}}{4}}\right]U_{\mbox{d}}\alpha sin\gamma$

(29)

$\alpha {\dot {\gamma }}=\sigma \alpha -{\frac {7b^{2}{\alpha }^{3}}{12{\omega }_{\mbox{a}}^{3}}}-$${\frac {3c{\alpha }^{3}}{8{\omega }_{\mbox{a}}}}+\left[{\frac {bf_{2}}{2\left({\omega }_{\mbox{d}}^{2}-{\omega }_{\mbox{a}}^{2}\right)}}+\right.$$\left.{\frac {f_{1}}{4}}\right]U_{\mbox{d}}\alpha sin\gamma$

(30)

where γ = ε²σT₀ - θ. The stationary solution equation α_s and γ_s can be obtained from Equations (29) and (30) by letting ${\dot {\alpha }}={\dot {\gamma }}=0$ . The vibration amplitude can be expressed as:

${\alpha }_{\mbox{s}}^{2}={\frac {2\sigma \pm {\sqrt {{\left({\frac {bf_{2}}{{\omega }_{\mbox{d}}^{2}-{\omega }_{\mbox{a}}^{2}}}+{\frac {f_{1}}{2}}\right)}^{2}U_{\mbox{d}}^{2}-{\chi }_{\mbox{a}}^{2}}}}{{\frac {7b^{2}}{6{\omega }_{\mbox{a}}^{3}}}+{\frac {3c}{4{\omega }_{\mbox{a}}}}}}$

(31)

For α_s is a real number, from Equation (31), the stability boundary is given by:

$\vert \sigma \vert \leq {\frac {1}{2}}{\sqrt {{\left({\frac {bf_{2}}{{\omega }_{\mbox{d}}^{2}-{\omega }_{\mbox{a}}^{2}}}+{\frac {f_{1}}{2}}\right)}^{2}U_{\mbox{d}}^{2}-{\chi }_{\mbox{a}}^{2}}}$ , $\left({\frac {bf_{2}}{{\omega }_{\mbox{d}}^{2}-{\omega }_{\mbox{a}}^{2}}}+\right.$$\left.{\frac {f_{1}}{2}}\right)U_{\mbox{d}}>{\chi }_{\mbox{a}}$

(32)

or

$\vert \sigma \vert >{\frac {1}{2}}{\sqrt {{\left({\frac {bf_{2}}{{\omega }_{\mbox{d}}^{2}-{\omega }_{\mbox{a}}^{2}}}+{\frac {f_{1}}{2}}\right)}^{2}U_{\mbox{d}}^{2}-{\chi }_{\mbox{a}}^{2}}}$ , $\left({\frac {bf_{2}}{{\omega }_{\mbox{d}}^{2}-{\omega }_{\mbox{a}}^{2}}}+\right.$$\left.{\frac {f_{1}}{2}}\right)U_{\mbox{d}}>{\chi }_{\mbox{a}}$

(33)

From Equations (32) and (33), emerging conditions of subharmonic parametric resonance can be obtained and many parameters like axial force, length of arms, etc. affect the structural motion stability region. It is a kind of way to get the dynamic instability region. But the stable solutions are functions of structural characteristics, which are not including effect of time factor. Therefore, the numerical method is further adopted to analyze.

4. Numerical analysis

In this section, the nonlinear responses of the radial gate arm are investigated. As a example, considering a radial gate arm with box-section, the cross-sectional dimensions are shown in Figure 4. The arm length L = 16 m, mass per unit length m = 386 kg/m, and modulus of elasticity E = 210 GPa.

Figure 4. Arm section (unit: mm).

Applying the fourth-order Runge–Kutta method to Equation (15), the radial gate arm subjected to the parametric excitation are investigated and the effects of the excitation and structural factors are evaluated.

4.1 Effect of water head

According to the Chinese and American codes, the layout of the main beam adopts equivalent load method, and the hydrostatic force H₀ is a function of water head. Table 1 summarizes the Chinese super high head radial gates (> 80 m). We first chose the initial disturbance V(0) = 10⁻⁴ m, excitation amplitude A_d = 0.01 m, time t = 500 s, damping coefficient μ = 0, and water head H = 0 as a basic case. The Figure 5 shows the influence of water head on the vibration amplitudes, λ_d denotes the ratio of the external excitation frequency to the natural frequency.

Table 1. Super high head radial gates in China.

	Engineering	Location	Size of face plane width-height (m)	Water head (m)
1	Xiaowan	Lantsang	5.0×7.0	163.00
2	Shuibuya	Qing River	6.0×7.0	154.00
3	Xiaolangdi	Yalong River	4.8×4.8	140.00
4	Jinping-I	Yalong River	5.0×6.0	133.00
5	Laxiwa	Yellow River	4.0×6.0	132.00
6	Pubugou	Dadu River	6.5×8.0	126.28
7	Nuozhadu	Lantsang	5.0×8.5	126.00
8	Tianshengqiao-I	Pearl River	6.4×7.5	120.00
9	Dongjiang	Dong Jiang Lake	6.4×7.5	120.00
10	Longyangxia	Yellow River	5.0×7.0	120.00
11	Xiluodu	Jinsha River	6.0×6.7	105.50
12	Jiangpinghe	Loushui River	6.0×6.0	105.14
13	Longtan	Hongshui River	5.0×8.0	90.00
14	Dagangshan	Dadu River	6.0×6.6	83.00
15	Goupitan	Wujiang River	10.0×9.0	81.00

When the principal parametric resonance occurs, the amplitude of the arm increases with the water head, and the maximum displacement reaches 0.57 m. In the case of subharmonic parametric resonance, the change is contrasting, and the maximum amplitude is 0.18 m. Overall, both the types of parametric resonance responses change slowly, so supposes the axial hydrostatic force generated by 80 m water head in further analysis.

Figure 5. Vibration amplitude–water head curves.

4.2 Effects of excitation frequency ratio and initial end bending moment

Figure 6 depicts the relationship curves plotted between the vibration amplitude and the excitation frequency ratio when subjected to the initial end bending moment of different water heads. As shown in the figure, at H = 0, the peak value of parametric resonance only appears at λ_d = 1 and λ_d = 2. This phenomenon conforms to the current design code “Design Code for Steel Gate in Water Resources and Hydropower Projects (SL74-2013)” (2013). Figure 7 shows that when the initial end bending moment is not considered, the beat vibration phenomenon can be observed, and the beat amplitude value increases over time at the subharmonic parametric resonance. However, considering the initial end bending moment (Figure 6), the amplitude peak shifts and the arm response amplitudes increase with the initial end bending moment level and the resonance region become larger, which facilitates parametric resonance excitation. According to the design retaining water height of the gate, 80 m is the critical value between high head gates and super high head gates, and maximum existing design head is 163 m (Table 1). The maximum displacements of parametric resonance are 1.44 m and 2.20 m at water heads of 80 m and 163 m, respectively. Through the above analysis, it is further elucidated that the parametric resonance excitation condition given in Chinese code is a special case in the parametric resonance region, which only considers the parametric resonance occurring at nearly twice the natural frequency. In order to obtain a higher accuracy response of the arm, the proposed calculation method of the instability regions in this study should be adopted.

(a) H < 100m	(b) H ≥ 100 m

Figure 6. Vibration amplitude–excitation frequency ratio curves.

(a) λd = 1	(b) λd = 2

Figure 7. Displacement time–history curves

4.3 Effect of excitation amplitude

Figure 8. Vibration amplitude–excitation amplitude curves.

Figure 8 shows the relationship curve between the displacement and the excitation amplitude. The graph shows that when parametric resonance occurs, the displacement response is closely related to the excitation amplitude. The response increases nonlinearly with an increase in the excitation amplitude. When the excitation amplitude is less than 0.08 m, the subharmonic parametric resonance response is greater than the principal parametric resonance response. On the other hand, when the excitation amplitude is greater than 0.08 m, both the parametric resonance responses are almost equal. When the excitation amplitude is less than 0.01 m, the arm swings in a smaller range, making it difficult to form a destructive vibration. Hence, the response of the arm should be completely considered with respect to the factor of excitation amplitudes.

4.4 Effect of inclined angle

The effect of inclined angle on the vibration amplitude is presented in Figure 9. The amplitude displacement increases significantly at subharmonic parametric resonance with the angle. However, at parametric resonance, when the angle is less than 45 degrees, the amplitude increases rapidly and then remains steady. Hence, the inclined angle should be considered along with the displacement constraint in the design of radial gates.

Figure 9. Vibration amplitude–inclined angle curves.

4.5 Effect of arm length and excitation time

According to the characteristics of conventional radial gate and special planar radial gate, as mentioned in the study by Chen, F.Z. and Yan, G.H. [18], the arm length from 5 to 25 m is examined in this study. The relationship curve is plotted between the amplitude and the length of arm, as shown in Figure 10.

Figure 10. Vibration amplitude–length curves.

It can be found that the calculated results are not in agreement with the description of the conventional partitioning method for unstable regions In the past review, with an increase in the length of the arm, its flexibility and response amplitude should increase simultaneously. However, the inflection point appears in Figure 10, and a nonlinear curve dependence is noted. Then, the time history analysis at a, b, c, and d on both sides of the inflection are given in Figure 11. When the arm length is greater than 45 m and 20 m at λ_d = 1 and λ_d = 2, a longer excitation time is needed to reach the peak response. This result proves that the gate may not show a strong vibration at the initial stage of the local opening operation, while the vibration amplitude will become considerably higher across time until the structure reaches integral instability failure. The results also prove that it is not advisable to ignore the time factor in the conventional design.

From the above analysis, it is noted that the time factor of the local opening should be considered in the initial design stage. Meanwhile, in the stages of completion and operation, the corresponding guidelines of local opening range and opening–closing time should be given in detail according to the vibration analysis results in a project. It is practical to choose the length of radial gate arm to be less then 20 m in 500 s according to the case presented in this exploration.

(a) λd= 1, L = 45 (Vmax= 0.2291)	(b) λd= 1, L = 60 (Vmax= 0.2427)
(c) λd = 2, L = 10 (Vmax = 0.8116)	(d) λd = 2, L = 20 (Vmax = 0.9604)

Figure 11. Displacement time-history curves

4.6 Effect of damping

(1) λd=1	(b) λd=2

Figure 12. Vibration amplitude–damping ratio curves.

Figure 12 illustrates the relationship between amplitude and damping ratio under different excitation amplitudes. As shown in the figure, the amplitude nonlinearly decreases with an increase in the damping ratio. Furthermore, the amplitude decreases rapidly at λ_d = 1, ξ < 0.01, and then remains steady. The damping decrement at λ_d = 2 is slightly higher than that at λ_d = 1. Overall, the damping effect to vibration control is not very obvious, once the parametric resonance excitation is achieved, the resonance amplitudes are significantly affected by the nonlinear term of system rather than the damping coefficient. Hence, when an energy transfer method, such as adding tuned mass dampers, is employed, the vibration suppressed effectiveness is better.

5. Conclusions

This paper presents the nonlinear parametric vibration of radial gates. The nonlinear differential equation of motion for the radial gate arm is derived by considering the geometrically nonlinear. Furthermore, the stability and vibration response of principal and subharmonic parametric resonance cases are investigated by using multi-scale and numerical methods. At the same time, the numerical analysis overcomes the limitations of the former methods, in which the dynamic instability regions are divided without considering the effects of response amplitude and time. The results are discussed as follows:

(1) The change in the hydrostatic force caused by the water head has a negligible effect on the parametric vibration, while the resonance region broadens with an increase in the initial bending moment.

(2) An increase in amplitude response leads to an increase in both the excitation amplitude and arm length and a decrease in the inclined angle.

(3) The vibration amplitude becomes significantly higher with time. Therefore, the time factor should be considered in the design and operation stages.

(4) Comparing with damping effect, the vibration response mainly depends on the nonlinear behavior. Thus, the energy transfer method should be adopted to suppress the vibration.

Acknowledgements

This research is supported by National Natural Science Foundation of China (Grant No. 51179164 and 51478354).

References

[1] CMWR (China Ministry of Water Resources). (2013) Design Code for Steel Gate in Water Resources and Hydropower Projects (SL74-2013). China Water & Power Press, Beijing.

[2] Niu, Z.G., Li, T.C. (2008) Research on Dynamic Stability of Steel Radial Gates. In: Earth and Space. Long Beach. pp. 1-8.

[3] Oh, L.S., Hoje, S., Won, K.J. (2018) Flow-induced vibration of a radial gate at various opening heights. Eng. Appl. Comp. Fluid, 12: 567-583.

[4] Brusewicz, K. Sterpejkowicz-Wersocki, W. Jankowski, R. (2017) Modal analysis of a steel radial gate exposed to different water levels. Arch. Hydro Eng. Environ. Mech., 64: 37-47.

[5] Wang, Z.Z., Zhang, X.C., Liu, J.L. (2017) Advances and developing trends in research of large hydraulic steel gates. J. Hydroelectr. Eng., 36: 1-18.

[6] Demsic, M., Uros, M., Lazarevic, A.J., Lazarevic, D. (2019) Resonance regions due to interaction of forced and parametric vibration of a parabolic cable. J. Sound. Vib., 447: 78-104.

[7] Zhang, X., Peng, J., Wang, L. (2014) Parametric resonances in the two-to-one resonant beams on elastic foundation. Nonlinear Dynam., 77: 339-352.

[8] Franzini, G.R., Pesce, C.P., Gonçalves, R.T., Fujarra, A.L.C., Mendes, P. (2018). An experimental investigation on concomitant vortex-induced vibration and axial top-motion excitation with a long flexible cylinder in vertical configuration. Ocean Engineering, 156: 596-612.

[9] Wei, M.H., Lin, K., Jin,L., Zou, D.J. (2016) Nonlinear dynamics of a cable-stayed beam driven by sub-harmonic and principal parametric resonance. Int. J. Mech. Sci., 110: 78-93.

[10] Zhang, J.G. (1985) Summary on the study of gate vibration in China. Water Power, 11: 38-44.

[11] Liu, J.L., Wang, Z.Z., Fang, X., Fang, H.M. (2011) Dynamic instability mechanism and vibration control of radial gate arms. Appl. Mech. Mater., 50: 309-313.

[12] Zhang, J., Xie, Z.X. (2011) Nonlinear vibration and chaos phenomena of arm structures in radial gate. Water Power, 50: 309-313.

[13] Cai, K. Zhang, C. (2011) An optimal construction of a hydropower arch gate. Adv. Mater. Res., 346: 109-115.

[14] Zhang, X.C., Wang, Z.Z., Sun, D.X. (2018) Research on rational layout of strut arms of tainter gate in vertical frame. Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 34: 1-6.

[15] Hu, M.S., Yang, Z.Z., Xu, J., Zhang, B. (2015) Modal parameter identification of radial gate based on stochastic subspace method. Water resour. power, 33: 164-167.

[16] Nayfeh, A. H., Mook, D. T. (1979) Nonlinear Oscillations. Wiley Publisher, New York.

[17] Chen, Y. S. (2002) Nonlinear Vibration, Higher Education Press, Beijing.

[18] Chen, F.Z., Yan, G.H. (2015) The key technology research on Chinese special type sluice. Hohai University Press, Nanjing.