First-order dynamic sliding mode variable structure-based bridge crane anti-swing method

Abstract

A bridge crane anti-swing method based on first-order dynamic sliding mode variable structure includes steps of: constructing a two-dimensional bridge crane system model and a crane system control model, respectively; performing differentiation on two sliding mode surfaces containing swing angle dynamic change and rope length dynamic change to obtain a crane position dynamic sliding surface s1 and a rope length dynamic sliding mode surface s2, respectively; combining a displacement x, a length l and a swing angle θ in the two-dimensional bridge crane system model with the crane position dynamic sliding surface s1 and rope length dynamic sliding mod surface s2 in the crane system control model to obtain a relationship among a horizontal traction force f1, an along-rope traction force f2, the displacement x, the length l and the swing angle θ.

Claims

1. A bridge crane anti-swing method based on first-order dynamic sliding mode variable structure, comprising steps of: monitoring parameters of a bridge crane system in real time, and establishing a two-dimensional bridge crane system model having a weight, and a bridge crane for hoisting the weight by a rope connecting between the weight and the bridge crane; determining, according to the two-dimensional bridge crane system model, position coordinates of the bridge crane and the weight, respectively by formulae of: $\{\begin{array}{c}{x}_M=x\\ y_M=0\\ x_m=x-l\sin \theta \\ y_m=l\cos \theta \end{array}$ where x.sub.M and y.sub.M are the X-axis coordinate and Y-axis coordinate of the bridge crane, respectively; x.sub.m and y.sub.m are the X-axis coordinate and Y-axis coordinate of the weight, respectively; θ is the swing angle of the weight during swing; l is the length of the rope; x is the displacement of the weight in a horizontal direction during swing; and, the displacement x, the length l and the swing angle θ are parameters of the bridge crane system to be monitored in real time; determining formulae of second-order derivatives of the displacement x, the length l and the swing angle θ, respectively, according to the two-dimensional bridge crane system model: $\stackrel{.Math.}{x}=\frac{1}{M}\left(f_1=D\stackrel{.}{x}+f_2\sin \theta \right)\stackrel{.Math.}{l}=g\cos \theta +l{\stackrel{.}{\theta }}^2+\frac{\left(f_1-D\stackrel{.}{x}\right)\sin \theta }{M}+\frac{f_2\left(M+m{\sin }^2\theta \right)}{\mathrm{Mm}}\stackrel{.Math.}{\theta }=\frac{\left(f_1-D\stackrel{.}{x}+f_2\sin \theta \right)}{\mathrm{Ml}}-2\frac{\stackrel{.}{l}\stackrel{.}{\theta }}{l}-\frac{g\sin \theta }{l}$ where M is the mass of the bridge crane, f.sub.1 is a horizontal traction force suffered by the bridge crane, D is a resistance coefficient, f.sub.2 is an along-rope traction force suffered by the weight, and m is the mass of the weight; constructing a crane system control model corresponding to the two-dimensional bridge crane system model, the crane system control model comprising a crane position dynamic sliding mode surface s.sub.1 containing a swing angle dynamic change and a rope length dynamic sliding mode surface s.sub.2 containing a rope length dynamic change, and the crane position dynamic sliding mode surface s.sub.1 and the rope length dynamic sliding mode surface s.sub.2 being constructed by the following formulae: $\{\begin{array}{c}{s}_1=\stackrel{.Math.}{x}+a\stackrel{.}{x}+b\left(x-x_d\right)+c\stackrel{.}{\theta }+d\theta \\ s_2=\stackrel{.Math.}{l}+a_1\stackrel{.}{l}+b_1\left(l-l_d\right)\end{array}$ where x.sub.d is a displacement reference value of the bridge crane and is a constant; l.sub.d is a rope length reference value of the rope and is a constant; and, a, b, c, d, a.sub.1 and b.sub.1 are all constants greater than 0; performing first-order derivation on the crane position dynamic sliding mode surface s.sub.1 and the rope length dynamic sliding mode surface s.sub.2 by using an exponential approach law control method to obtain a first-order crane position dynamic sliding mode surface {dot over (s)}.sub.1 and a first-order rope length dynamic sliding mode surface {dot over (s)}.sub.2, which are separately expressed by the following formulae: $\{\begin{array}{c}{\stackrel{.}{s}}_1=\stackrel{.Math.}{x}+a\stackrel{.Math.}{x}+b\stackrel{.}{x}+c\stackrel{.Math.}{\theta }+d\stackrel{.}{\theta }\\ {\stackrel{.}{s}}_2=\stackrel{.Math.}{l}+a_1\stackrel{.Math.}{l}+b_1\stackrel{.}{l}\end{array};$ combining the second-order derivatives of the displacement x, the length l and the swing angle θ with the first-order crane position dynamic sliding mode surface {dot over (s)}.sub.1 and a first-order rope length dynamic sliding mode surface {dot over (s)}.sub.2 to obtain formulae of first-order derivatives of the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2: $\{\begin{array}{c}{\stackrel{.}{f}}_1=-\sin \theta {\stackrel{.}{f}}_2+\frac{D}{M}\left(f_1-D\stackrel{.}{x}+f_2\sin \theta \right)-\stackrel{.}{\theta }\cos \theta f_2-\\ M\left(a\stackrel{.Math.}{x}+b\stackrel{.}{x}+c\stackrel{.Math.}{\theta }+d\stackrel{.}{\theta }+{\epsilon }_1\phi \left(\lambda ,s_1\right)+k_1{s}_1\right)\\ f_2=-\frac{m\sin \theta }{M+m{\sin }^2\theta }{f}_1-\frac{\mathrm{Mm}}{M+m{\sin }^2\theta }[\frac{\theta \sin 2\theta }{M}{f}_2+\\ \frac{M\left(f_1-D\stackrel{.}{x}\right)\stackrel{.}{\theta }\cos \theta -D\left(f_1-D\stackrel{.}{x}+f_2\sin \theta \right)\sin \theta }{M^2}-\\ g\stackrel{.}{\theta }\sin \theta +\stackrel{.}{l}{\stackrel{.}{\theta }}^2+2l\stackrel{.}{\theta }\stackrel{.Math.}{\theta }+a_1\stackrel{.Math.}{l}+b_1\stackrel{.}{l}+{\epsilon }_2\phi \left(\lambda ,s_2\right)+k_1{s}_2]\end{array};$ where λ is a positive constant, and the function $\phi \left(\lambda ,s_i\right)=\frac{1-e^{-{\lambda }_{s_i}}}{1+e^{-{\lambda }_{s_i}}}$ is a bipolar Sigmoid function which has the same definition as a sign function sgn(s.sub.i); integrating the first-order derivatives of the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2 by integrators, respectively, and cyclically adjusting and outputting the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2 in combination with the displacement x, the length l and the swing angle θ that change in real time; and applying, to the bridge crane and the weight, equivalent acting forces opposite to the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2, respectively, to prevent the bridge crane system from swinging.

2. The bridge crane anti-swing method based on first-order dynamic sliding mode variable structure according to claim 1, wherein, in the step of determining formulae of second-order derivatives of the displacement x, the length l and the swing angle θ, respectively, according to the two-dimensional bridge crane system model comprises: determining a Lagrange equation of the two-dimensional bridge crane system model according to the displacement x, the length l, the swing angle θ, the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2: $\{\begin{array}{c}\frac{d}{\mathrm{dt}}\left(\frac{\partial T}{\partial \stackrel{\text{'}}{x}}\right)-\frac{\partial T}{\partial x}=f_1-D\stackrel{.}{x}\\ \frac{d}{\mathrm{dt}}\left(\frac{\partial T}{\partial \stackrel{.}{l}}\right)-\frac{\partial T}{\partial l}=f_2+\mathrm{mg}\cos \theta \\ \frac{d}{\mathrm{dt}}\left(\frac{\partial T}{\partial \stackrel{.}{\theta }}\right)-\frac{\partial T}{\partial \theta }=-\mathrm{mg}\sin \theta \end{array};$ determining a nonlinear dynamic differential equation according to the Lagrange equation of the two-dimensional bridge crane system model: $\{\begin{array}{c}\left(M+m\right)\stackrel{.Math.}{x}-m\stackrel{.Math.}{l}\sin \theta -\mathrm{ml}\stackrel{.Math.}{\theta }\cos \theta \\ -2m\stackrel{.}{\theta }\stackrel{.}{l}\cos \theta +\mathrm{ml}{\stackrel{.}{\theta }}^2\sin \theta =f_1-D\stackrel{.}{x}\\ m\stackrel{.Math.}{l}-m\stackrel{.Math.}{x}\sin \theta -ml{\stackrel{.}{\theta }}^2-\mathrm{mg}\cos \theta =f_2\\ {\mathrm{ml}}^2\stackrel{.Math.}{\theta }-\mathrm{ml}\stackrel{.Math.}{x}\cos \theta +2\mathrm{ml}\stackrel{.}{l}\stackrel{.}{\theta }+\mathrm{mg}l\sin \theta =0\end{array};$ and performing formal transformation on the nonlinear dynamic differential equation of the two-dimensional bridge crane system model to obtain formulae of second-order derivatives of the displacement x, the length l and the swing angle θ, respectively.

3. The bridge crane anti-swing method based on first-order dynamic sliding mode variable structure according to claim 1, wherein: the crane system control model further comprises a first linear operation module, a first controller, a second linear operation module, a second controller, a crane system module, a first input module for inputting the displacement reference value x.sub.d and a second input module for inputting the rope length reference value l.sub.d; the first input module, the first linear operation module, the crane position dynamic sliding mode surface s.sub.1, the first controller and the crane system module are connected successively, and the crane system module feeds back and outputs the swing angle θ and the displacement x to the first linear operation module and the first input module, respectively; and, the second input module, the second linear operation module, the rope length dynamic sliding mode surface s.sub.2, the second controller and the crane system module are connected successively, and the crane system module feeds back and outputs the length l to the second input module.

4. The bridge crane anti-swing method based on first-order dynamic sliding mode variable structure according to claim 1, wherein, in the step of performing first-order derivation on the crane position dynamic sliding mode surface s.sub.1 and the rope length dynamic sliding mode surface s.sub.2, the relationship between the exponential approach law control method and the first-order crane position dynamic sliding mode surface {dot over (s)}.sub.1 and the first-order rope length dynamic sliding mode surface {dot over (s)}.sub.2 is expressed as: ${\stackrel{.}{s}}_i=-{\epsilon }_i\phi \left(\lambda ,s_i\right)-k_i{s}_{i}i=1,2\mathrm{where},\phi \left(\lambda ,s_i\right)=\{\begin{array}{cc}1& \mathrm{if}{s}_i>0\\ 0& \mathrm{if}{s}_i=0\\ -1& \mathrm{if}{s}_i<0\end{array}.$

5. The bridge crane anti-swing method based on first-order dynamic sliding mode variable structure according to claim 4, wherein, in the step of combining the second-order derivatives of the displacement x, the length l and the swing angle θ with the first-order crane position dynamic sliding mode surface {dot over (s)}.sub.1 and a first-order rope length dynamic sliding mode surface {dot over (s)}.sub.2 to obtain formulae of first-order derivatives of the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2 comprises: separately substituting the second-order derivatives of the displacement x, the length l and the swing angle θ into the formulae of the first-order crane position dynamic sliding mode surface {dot over (s)}.sub.1 and the first-order rope length dynamic sliding mode surface {dot over (s)}.sub.2 to obtain the following formulae: $\{\begin{array}{c}{\stackrel{.}{s}}_1=\frac{1}{M}\left({\stackrel{.}{f}}_1-D\stackrel{.Math.}{x}+{\stackrel{.}{f}}_2\sin \theta +f_2\stackrel{.}{\theta }\cos \theta \right)+a\stackrel{.Math.}{x}+b\stackrel{.}{x}+c\stackrel{.Math.}{\theta }+d\stackrel{.}{\theta }\\ {\stackrel{.}{s}}_2=-g\stackrel{.}{\theta }\sin \theta +\stackrel{.}{l}{\stackrel{.}{\theta }}^2+2l\stackrel{.}{\theta }\stackrel{.Math.}{\theta }+\frac{1}{M}\left[\left({\stackrel{.}{f}}_1-D\stackrel{.Math.}{x}\right)\sin \theta +\left(f_1-D\stackrel{.}{x}\right)\stackrel{.}{\theta }\cos \theta \right]+\\ \frac{1}{\mathrm{Mm}}\left[{\stackrel{.}{f}}_2\left(M+m{\sin }^2\theta \right)+2f_{2}m\stackrel{.}{\theta }\sin \theta \cos \theta \right]+a_1\stackrel{.Math.}{l}+b_1\stackrel{.}{l}\end{array};$ and separately substituting the second-order derivatives of the displacement x, the length l and the swing angle θ as well as the relationship between the exponential approach law control method and the first-order crane position dynamic sliding mode surface {dot over (s)}.sub.1 and the first-order rope length dynamic sliding mode surface {dot over (s)}.sub.2 into the formulae of the first-order crane position dynamic sliding mode surface {dot over (s)}.sub.1 and the first-order rope length dynamic sliding mode surface {dot over (s)}.sub.2 to obtain formulae of first-order derivatives of the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2.

6. The bridge crane anti-swing method based on first-order dynamic sliding mode variable structure according to claim 3, wherein: the integrators are arranged in the first controller and the second controller, respectively, and the first controller and the second controller output the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2 to the crane system module, respectively.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) The present disclosure will be further described below by embodiments with reference to the accompanying drawings.

(2) FIG. 1 is a schematic view of a two-dimensional bridge crane system model; and

(3) FIG. 2 is a schematic view of a crane system control model.

DETAILED DESCRIPTION

(4) With reference to FIGS. 1-2, the present disclosure provides a bridge crane anti-swing method based on first-order dynamic SMVS. In view of the control problem of a two-dimensional under-actuated bridge crane in the case of a long rope, two sliding mode surfaces containing a swing angle dynamic change and a rope length dynamic change are differentiated to obtain a crane position dynamic sliding surface s.sub.1 and a rope length dynamic sliding mode surface s.sub.2, respectively. Since the crane position dynamic sliding surface s.sub.1 and the rope length dynamic sliding mode surface s.sub.2 are related to the control inputs of the crane system control model, i.e., first-order derivatives of a horizontal traction force f.sub.1 and an along-rope traction force f.sub.2, a dynamic sliding mode control law that is continuous in time is obtained, so that the anti-swing positioning control of the bridge crane system is realized, and the chattering phenomenon of the SMVS control can be effectively weakened.

(5) In particular, the parameters of the bridge crane system are monitored in real time, and a two-dimensional bridge crane system model is established. With reference to FIG. 1, position coordinates of a bridge crane 2 and a weight 1 are as follows:

(6) $\{\begin{array}{c}{x}_M=x\\ y_M=0\\ x_m=x-l\sin \theta \\ y_m=l\cos \theta \end{array}$

(7) where x.sub.M and y.sub.M are the X-axis coordinate and Y-axis coordinate of the bridge crane 2, respectively; x.sub.m and y.sub.m are the X-axis coordinate and Y-axis coordinate of the weight 1, respectively; θ is the swing angle of the weight 1 during swing; l is the length of the rope; x is the displacement of the weight 1 in a horizontal direction during swing; and, the displacement x, the length l and the swing angle θ are parameters of the bridge crane system to be monitored in real time.

(8) Then, a Lagrange equation of the two-dimensional bridge crane system model can be obtained according to the two-dimensional bridge crane system model:

(9) $\{\begin{array}{c}\frac{d}{\mathrm{dt}}\left(\frac{\partial T}{\partial \stackrel{.}{x}}\right)-\frac{\partial T}{\partial x}=f_1-D\stackrel{.}{x}\\ \frac{d}{\mathrm{dt}}\left(\frac{\partial T}{\partial \stackrel{.}{l}}\right)-\frac{\partial T}{\partial l}=f_2+\mathrm{mg}\cos \theta \\ \frac{d}{\mathrm{dt}}\left(\frac{\partial T}{\partial \stackrel{.}{\theta }}\right)-\frac{\partial T}{\partial \theta }=-\mathrm{mg}\sin \theta \end{array};$

(10) Based on the above analysis, a nonlinear dynamic differential equation of the two-dimensional bridge crane system model can be obtained:

(11) $\{\begin{array}{c}\left(M+m\right)\stackrel{.Math.}{x}=m\stackrel{.Math.}{l}\sin \theta -ml\stackrel{.Math.}{\theta }\cos \theta \\ -2m\stackrel{.}{\theta }\stackrel{.}{l}\cos \theta +ml{\stackrel{.}{\theta }}^2\sin \theta =f_3-D\stackrel{.}{x}\\ m\stackrel{.Math.}{l}-m\stackrel{.Math.}{x}\sin \theta -ml{\stackrel{.}{\theta }}^2-\mathrm{mg}\cos \theta =f_2\\ {\mathrm{ml}}^2\stackrel{.Math.}{\theta }=\mathrm{ml}\stackrel{.Math.}{x}\cos \theta +2\mathrm{ml}\stackrel{.}{l}\stackrel{.}{\theta }+\mathrm{mg}l\sin \theta =0\end{array};$

(12) Formal transformation can be performed on the nonlinear dynamic differential equation to obtain formulae of second-order derivatives of the displacement x, the length l and the swing angle θ, respectively:

(13) $\stackrel{.Math.}{x}=\frac{1}{M}\left(f_1-D\stackrel{.}{x}+f_2\sin \theta \right)\stackrel{.Math.}{l}=g\cos \theta +l{\stackrel{.}{\theta }}^2+\frac{\left(f_i-D\stackrel{.}{x}\right)\sin \theta }{M}+\frac{f_2\left(M+m{\sin }^2\theta \right)}{\mathrm{Mm}}\stackrel{.Math.}{\theta }=\frac{\left(f_1-D\stackrel{.}{x}+f_2\sin \theta \right)\cos \theta }{Ml}-2\frac{\stackrel{.}{l}\stackrel{.}{\theta }}{l}-\frac{g\sin \theta }{l}$

(14) where M is the mass of the bridge crane 2, f.sub.1 is a horizontal traction force suffered by the bridge crane 2, D is a resistance coefficient, f.sub.2 is an along-rope traction force suffered by the weight 1, and m is the mass of the weight 1.

(15) Then, a crane system control model corresponding to the two-dimensional bridge crane system model is constructed. The crane system control model includes a crane position dynamic sliding mode surface s.sub.1 containing a swing angle dynamic change, a rope length dynamic sliding mode surface s.sub.2 containing a rope length dynamic change, a first linear operation module 3, a first controller 4, a second linear operation module 5, a second controller 6, a crane system module 7, a first input module 8 for inputting a displacement reference value x.sub.d and a second input module 9 for inputting a rope length reference value l.sub.d. The first input module 8, the first linear operation module 3, the crane position dynamic sliding mode surface s.sub.1, the first controller 4 and the crane system module 7 are connected successively, and the crane system module 7 feeds back and outputs the swing angle θ and the displacement x to the first linear operation module 3 and the first input module 8, respectively. The second input module 9, the second linear operation module 5, the rope length dynamic sliding mode surface s.sub.2, the second controller 6 and the crane system module 7 are connected successively, and the crane system module 7 feeds back and outputs the length l to the second input module 9. It can be seen from the structure of the crane system control model that there is mutual affection between the displacement x, the length l and the swing angle θ, and the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2. When equivalent acting forces opposite to the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2 are separately applied to the bridge crane 2 and the weight 1 to prevent the bridge crane system from swinging, the displacement x, the length l and the swing angle θ will correspondingly change. In this case, the feedback of the displacement x, the length l and the swing angle θ will affect the changes in the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2. Therefore, by cyclically adjusting and outputting the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2, the anti-swing positioning control of the bridge crane system is realized.

(16) The relationship between the displacement x, the length l, the swing angle θ and the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2 will be analyzed below.

(17) The crane position dynamic sliding mode surface s.sub.1 and the rope length dynamic sliding mode surface s.sub.2 are constructed by the following formulae:

(18) $\{\begin{array}{c}{s}_1=\stackrel{.Math.}{x}+a\stackrel{.}{x}+b\left(x-x_d\right)+c\stackrel{.}{\theta }+d\theta \\ s_2=\stackrel{.Math.}{l}+a_1\stackrel{.}{l}+b_1\left(l-l_d\right)\end{array}$

(19) where x.sub.d is the displacement reference value of the bridge crane 2 and is a constant; l.sub.d is the rope length reference value of the rope and is a constant; and, a, b, c, d, a.sub.1 and b.sub.1 are all constants greater than 0.

(20) First-order derivation is performed on the crane position dynamic sliding mode surface s.sub.1 and the rope length dynamic sliding mode surface s.sub.2 by using an exponential approach law control method to obtain a first-order crane position dynamic sliding mode surface {dot over (s)}.sub.1 and a first-order rope length dynamic sliding mode surface {dot over (s)}.sub.2, which are separately expressed by the following formulae:

(21) $\{\begin{array}{c}{\stackrel{.}{s}}_1=\stackrel{.Math.}{x}+a\stackrel{.Math.}{x}+b\stackrel{.}{x}+c\stackrel{.Math.}{\theta }+d\stackrel{.}{\theta }\\ {\stackrel{.}{s}}_2=\stackrel{.Math.}{l}+a_1\stackrel{.Math.}{l}+b_1\stackrel{.}{l}\end{array}$

(22) The relationship between the exponential approach law control method and the first-order crane position dynamic sliding mode surface {dot over (s)}.sub.1 and the first-order rope length dynamic sliding mode surface {dot over (s)}.sub.2 is expressed as:

(23) ${\stackrel{.}{s}}_i=-{\epsilon }_i\phi \left(\lambda ,s_i\right)-k_i{s}_{i}i=1,2\mathrm{where},\phi \left(\lambda ,s_i\right)=\frac{1-e^{-{\lambda }_{s_i}}}{1+e^{-{\lambda }_{s_i}}}=\{\begin{array}{cc}1& \mathrm{if}{s}_i>0\\ 0& \mathrm{if}{s}_i=0\\ -1& \mathrm{if}{s}_i<0\end{array};$
λ is a positive constant, and the function Φ(λ,s.sub.i) is a bipolar Sigmoid function which has the same definition as a sign function sgn(s.sub.i).

(24) Therefore, when the second-order derivatives of the displacement x, the length l and the swing angle θ are separately substituted into the formulae of the first-order crane position dynamic sliding mode surface {dot over (s)}.sub.1 and the first-order rope length dynamic sliding mode surface {dot over (s)}.sub.2, the following formulae can be obtained:

(25) $\{\begin{array}{c}{\stackrel{.}{s}}_1=\frac{1}{M}\left({\stackrel{.}{f}}_1-D\stackrel{.Math.}{x}+{\stackrel{.}{f}}_2\sin \theta +f_2\stackrel{.}{\theta }\cos \theta \right)+a\stackrel{.Math.}{x}+b\stackrel{.}{x}+c\stackrel{.Math.}{\theta }+d\stackrel{.}{\theta }\\ {\stackrel{.}{s}}_2=-g\stackrel{.}{\theta }\sin \theta +\stackrel{.}{l}{\stackrel{.}{\theta }}^2+2l\stackrel{.}{\theta }\stackrel{.Math.}{\theta }+\frac{1}{M}\left[\left({\stackrel{.}{f}}_1-D\stackrel{.Math.}{x}\right)\sin \theta +\left(f_1-D\stackrel{.}{x}\right)\stackrel{.}{\theta }\cos \theta \right]+\\ \frac{1}{\mathrm{Mm}}\left[{\stackrel{.}{f}}_2\left(M+m{\sin }^2\theta \right)+2f_{2}m\stackrel{.}{\theta }\sin \theta \cos \theta \right]+a_1\stackrel{.Math.}{l}+b_1\stackrel{.}{l}\end{array};$

(26) Then, the second-order derivatives of the displacement x, the length l and the swing angle θ as well as the relationship between the exponential approach law control method and the first-order crane position dynamic sliding mode surface {dot over (s)}.sub.1 and the first-order rope length dynamic sliding mode surface {dot over (s)}.sub.2 are separately substituted into the formulae of the first-order crane position dynamic sliding mode surface {dot over (s)}.sub.1 and the first-order rope length dynamic sliding mode surface {dot over (s)}.sub.2 to obtain formulae of first-order derivatives of the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2:

(27) $\{\begin{array}{c}{\stackrel{.}{f}}_1=-\sin \theta {\stackrel{.}{f}}_2+\frac{D}{M}\left(f_1-D\stackrel{.}{x}+f_2\sin \theta \right)-\stackrel{.}{\theta }\cos \theta f_2-\\ M\left(a\stackrel{.Math.}{x}+b\stackrel{.}{x}+c\stackrel{.Math.}{\theta }+d\stackrel{.}{\theta }+{\epsilon }_1\phi \left(\lambda ,s_1\right)+k_1{s}_1\right)\\ {\stackrel{.}{f}}_2=-\frac{m\sin \theta }{M+m{\sin }^2\theta }{\stackrel{.}{f}}_1-\frac{\mathrm{Mm}}{M+m{\sin }^2\theta }[\frac{\stackrel{.}{\theta }\sin 2\theta }{M}{f}_2+\\ \frac{M\left(f_1-D\stackrel{.}{x}\right)\stackrel{.}{\theta }\cos \theta -D\left(f_1-D\stackrel{.}{x}+f_2\sin \theta \right)\sin \theta }{M^2}-\\ g\stackrel{.}{\theta }\sin \theta +\stackrel{.}{l}{\stackrel{.}{\theta }}^2+2l\stackrel{.}{\theta }\stackrel{.Math.}{\theta }+a_3\stackrel{.Math.}{l}+b_1\stackrel{.}{l}+{\epsilon }_2\phi \left(\lambda ,s_2\right)+k_3{s}_2]\end{array};$

(28) At this time, the above formulae indicate relationship of the mutual affections among the displacement x, the length l, the swing angle θ, and the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2. Since the integrators are arranged in the first controller 4 and the second controller 6, and the first controller 4 and the second controller 6 output the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2 to the crane system module 7, respectively, the first-order derivatives of the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2 are integrated in the first controller 4 and the second controller 6 and then converted into the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2. Then, the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2 are input into the crane system module 7 to change the displacement x, the length l, the swing angle θ. At this time, the feedback of the displacement x, the length l, the swing angle θ will affect the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2, so that the output values of the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2 are cyclically adjusted continuously. At this time, equivalent acting forces opposite to the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2 are separately applied to the bridge crane 2 and the weight 1 according to the changes of the output values of the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2, so that the anti-swing positioning control of the bridge crane system is realized.

(29) In addition, since a conventional sign function is replaced with a bipolar Sigmoid function in the formulae of the first-order derivatives of the horizontal traction force f.sub.1 and the along-rope traction force f.sub.2, and the bipolar Sigmoid function is more derivable than the conventional sign function, smooth transition can be realized, and state points in a certain range are quickly attracted to a certain A neighborhood of a switching surface. This neighborhood is a boundary layer of the sliding mode switching surface. In this boundary layer, the quasi-sliding mode is not required to satisfy the conditions for the sliding mode, so the switching of the control structure on the switching surface is not required, and the chattering phenomenon can be fundamentally avoided or weakened.

(30) Although the preferred embodiments of the present disclosure have been described in detail above, the present disclosure is not limited the embodiments. Those skilled in the art may make various equivalent variations or substitutions without departing from the spirit of the present disclosure, and these equivalent variations or substitutions shall fall into the scope defined by the appended claims of the present application.