CN114890305A - Self-adaptive backstepping nonlinear swing suppression control method of double-pendulum rotary crane - Google Patents

Self-adaptive backstepping nonlinear swing suppression control method of double-pendulum rotary crane Download PDF

Info

Publication number
CN114890305A
CN114890305A CN202210589719.8A CN202210589719A CN114890305A CN 114890305 A CN114890305 A CN 114890305A CN 202210589719 A CN202210589719 A CN 202210589719A CN 114890305 A CN114890305 A CN 114890305A
Authority
CN
China
Prior art keywords
adaptive
double
rotary crane
self
pendulum
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210589719.8A
Other languages
Chinese (zh)
Other versions
CN114890305B (en
Inventor
朱鸿杰
郗焕
欧阳慧珉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Tech University
Original Assignee
Nanjing Tech University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Tech University filed Critical Nanjing Tech University
Priority to CN202210589719.8A priority Critical patent/CN114890305B/en
Publication of CN114890305A publication Critical patent/CN114890305A/en
Application granted granted Critical
Publication of CN114890305B publication Critical patent/CN114890305B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B66HOISTING; LIFTING; HAULING
    • B66CCRANES; LOAD-ENGAGING ELEMENTS OR DEVICES FOR CRANES, CAPSTANS, WINCHES, OR TACKLES
    • B66C13/00Other constructional features or details
    • B66C13/04Auxiliary devices for controlling movements of suspended loads, or preventing cable slack
    • B66C13/06Auxiliary devices for controlling movements of suspended loads, or preventing cable slack for minimising or preventing longitudinal or transverse swinging of loads
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B66HOISTING; LIFTING; HAULING
    • B66CCRANES; LOAD-ENGAGING ELEMENTS OR DEVICES FOR CRANES, CAPSTANS, WINCHES, OR TACKLES
    • B66C13/00Other constructional features or details
    • B66C13/04Auxiliary devices for controlling movements of suspended loads, or preventing cable slack
    • B66C13/08Auxiliary devices for controlling movements of suspended loads, or preventing cable slack for depositing loads in desired attitudes or positions
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B66HOISTING; LIFTING; HAULING
    • B66CCRANES; LOAD-ENGAGING ELEMENTS OR DEVICES FOR CRANES, CAPSTANS, WINCHES, OR TACKLES
    • B66C13/00Other constructional features or details
    • B66C13/18Control systems or devices
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B66HOISTING; LIFTING; HAULING
    • B66CCRANES; LOAD-ENGAGING ELEMENTS OR DEVICES FOR CRANES, CAPSTANS, WINCHES, OR TACKLES
    • B66C23/00Cranes comprising essentially a beam, boom, or triangular structure acting as a cantilever and mounted for translatory of swinging movements in vertical or horizontal planes or a combination of such movements, e.g. jib-cranes, derricks, tower cranes
    • B66C23/88Safety gear
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Jib Cranes (AREA)
  • Control And Safety Of Cranes (AREA)

Abstract

The invention discloses a self-adaptive backstepping nonlinear swing suppression control method of a double-pendulum rotary crane, which comprises the following steps: establishing a mathematical model of the double-pendulum rotary crane based on a Lagrange kinetic equation and analyzing characteristics; establishing an energy function of a mathematical model of the rotary crane according to the characteristics, establishing a Lyapunov function based on a backstepping method and a system control target, and establishing an adaptive controller based on the energy function and the Lyapunov function; and obtaining the self-adaptive control rate based on the Lyapunov function, and adding the self-adaptive control rate into the self-adaptive controller. The design of the self-adaptive control rate has good robustness on system parameters of different rotary cranes; the control method of the PD-like in the self-adaptive controller can realize accurate positioning, and external disturbance is restrained through the feedback of the swing angle information and the design of an error item; the two are combined to realize efficient track tracking and swing suppression of different crane systems, so that a control effect is achieved.

Description

Self-adaptive backstepping nonlinear swing suppression control method of double-pendulum rotary crane
Technical Field
The invention relates to the technical field of motion control of an under-actuated crane system, in particular to a self-adaptive backstepping nonlinear swing suppression control method of a double-pendulum rotary crane.
Background
Under-actuated systems, i.e. systems where the system inputs fewer degrees of freedom than the system. Among these, crane systems are typical under-actuated systems, where the angle of the load and hook is such that it cannot directly control the force. And the device has the advantages of simple structure, low power consumption, few actuating mechanisms, wide application occasions and the like. Generally, the conventional swing eliminating method of a crane is mostly operated by people, but this causes many problems because people's operation is different from machine operation, and accurate positioning and swing elimination cannot be realized. The imprecise operation of which can cause errors and can even be dangerous.
Disclosure of Invention
This section is for the purpose of summarizing some aspects of embodiments of the invention and to briefly introduce some preferred embodiments. In this section, as well as in the abstract and the title of the invention of this application, simplifications or omissions may be made to avoid obscuring the purpose of the section, the abstract and the title, and such simplifications or omissions are not intended to limit the scope of the invention.
The present invention has been made in view of the above-mentioned problems.
Therefore, the technical problem solved by the invention is as follows: under the conditions that fluctuation and rotation actions of a cantilever of the double-pendulum rotary crane influence the swinging of a load and the load cannot be directly controlled, how to simultaneously realize accurate positioning of the cantilever and fast restrain the swinging of a hook and the load.
In order to solve the technical problems, the invention provides the following technical scheme: a self-adaptive backstepping nonlinear swing suppression control method of a double-pendulum rotary crane comprises the following steps:
establishing a mathematical model of the double-pendulum rotary crane based on a Lagrange kinetic equation and analyzing characteristics;
establishing an energy function of a mathematical model of the rotary crane according to the characteristics, and establishing an adaptive controller based on the energy function so as to solve the problem of unknown model parameters;
based on a back stepping method and a system control target, establishing a Lyapunov function, and optimizing the self-adaptive controller based on the Lyapunov function; and controlling the swing suppression of the double-pendulum rotary crane through an adaptive controller.
As a preferable scheme of the adaptive backstepping nonlinear swing suppression control method of the double-pendulum rotary crane, the method comprises the following steps: the mathematical model of the rotary crane is expressed as:
order to
Figure BDA0003664617360000021
Figure BDA0003664617360000022
G(q)=[g 1 g z g3g4g5g6]T
N=[OOOOn5n6]T
F f =[f 1 f 2 f 3 f 4 f 5 f 6 ] T
Figure BDA0003664617360000023
q=[ζ 1 ζ 2 ζ 3 ζ 4 ζ 5 ζ 6 ] T
Wherein: m (q) is the inertia matrix of the double pendulum rotary crane system,
Figure BDA0003664617360000024
is a centripetal-Coriolis matrix, G (q) is a gravity vector, N is a control input vector, F f Is the mechanical friction force of the double-pendulum rotary crane system, q is the state variable of the double-pendulum rotary crane system,
Figure BDA0003664617360000025
in the form of the first derivative of the signal,
Figure BDA0003664617360000026
is its second derivative; m is h And m l Mass of hook and load, M 0 And M 1 Mass of the cantilever and ballast, respectively, /) 1 And l 2 Respectively the length from the suspension rope and the lifting hook to the center of mass of the load, L 1 L is the length of the ballast and the length of the cantilever, I x 、I y 、I z For the moment of inertia of the cantilever in the x, y, z axes, respectively, I b Is the moment of inertia of the ballast, g is the acceleration of gravity, and ζ is a generalized quantity of state describing a rotary crane system i (i 1.., 4) is a swing angle between the hook and the load, and ζ 5 is a cantilever heave direction driving torque, ζ 5 is a driving force/torque 6 For driving torque in the direction of rotation of the cantilever, f 5 、f 6 Mechanical friction in the direction of cantilever heave and rotation, d l And (1., 4) are air friction parameters.
As a preferable scheme of the adaptive backstepping nonlinear swing suppression control method of the double-pendulum rotary crane, the method comprises the following steps: the mathematical model of the rotary crane further comprises,
establishing a friction force feedforward compensation model to eliminate friction force generated by a driving mechanism of the double-pendulum spiral crane, wherein the friction force feedforward compensation model is expressed as follows:
Figure BDA0003664617360000027
Figure BDA0003664617360000028
wherein ,f51 、f 52 、f 61 、f 62 、ε 1 and ε2 Feedforward compensation model parameters for friction force, f 51 and f61 Value of (2) and maximum static friction forceCorresponds to f 52 and f62 Is the coefficient of viscous friction, ε 1 and ε2 Is the static coefficient of friction.
As a preferable scheme of the adaptive backstepping nonlinear swing suppression control method of the double-pendulum rotary crane, the method comprises the following steps: the analysis characteristics include analysis of the swing characteristics, load (lever translation characteristics), and the like.
As a preferable scheme of the adaptive backstepping nonlinear swing suppression control method of the double-pendulum rotary crane, the method comprises the following steps: the energy function is expressed as:
Figure BDA0003664617360000031
wherein :
Figure BDA0003664617360000032
respectively representing the cantilever load, hook, swing angle ζ i A speed signal of 1, 4,
Figure BDA0003664617360000033
is the kinetic energy part of a double-pendulum rotary crane system,
Figure BDA0003664617360000034
the potential energy portion of its load and hook.
As a preferable scheme of the adaptive backstepping nonlinear swing suppression control method of the double-pendulum rotary crane, the method comprises the following steps: the establishing of the Lyapunov function includes,
the lyapunov equation V is established based on a backstepping method and a system control target according to the dynamic rule of a rotary crane model 1
Figure BDA0003664617360000035
wherein :e1 =[α 1 α 2 α 3 α 4 α 5 α 6 ] T ,q d =[ζ 1d ζ 2d ζ 3d ζ 4d ζ 5d ζ 6d ] T ,e 1 =q-q d .e 1 The error of the swing angle and the error of the current position and the target position of the cantilever,
Figure BDA0003664617360000036
the error of the swing angular speed and the error of the rotation speed of the cantilever are obtained;
for the Lyapunov equation V 1 The derivation is carried out to obtain:
Figure BDA0003664617360000037
wherein :k1 Is a positive parameter, e 2 =Θ-Θ d =[β 1 β 2 β 3 β 4 β 5 β 6 ] T In order to make the derivative of the Lyapunov equation semi-negative, design
Figure BDA0003664617360000038
Then there are
Figure BDA0003664617360000039
Figure BDA00036646173600000310
Design Lyapunov equation V based on backstepping method 2 Expressed as:
Figure BDA00036646173600000311
wherein ,
Figure BDA0003664617360000041
to adapt the error values of the column vectors of the system matrix, there are
Figure BDA0003664617360000042
Figure BDA0003664617360000043
Figure BDA0003664617360000044
Is composed of
Figure BDA0003664617360000045
The error value of the adaptive system vector of (1) is
Figure BDA0003664617360000046
As a preferable scheme of the adaptive backstepping nonlinear swing suppression control method of the double-pendulum rotary crane, the method comprises the following steps: the adaptive controller is represented as:
Figure BDA0003664617360000047
lyapunov equation V by design 2 Positive definite, its derivative is semi-negative definite, and can obtain self-adaptive control rate, so that the system can form self-adaptive controller with stable convergence.
As a preferable scheme of the adaptive backstepping nonlinear swing suppression control method of the double-pendulum rotary crane, the method comprises the following steps: the adaptive control rate is expressed as:
Figure BDA0003664617360000048
as a preferable scheme of the adaptive backstepping nonlinear swing suppression control method of the double-pendulum rotary crane, the method comprises the following steps: the adaptive controller may further comprise a processor for controlling the adaptive controller,
by adjusting k 1 、k 2 Thereby obtaining a gain that quickly adapts to different systems, where k 1 and k2 Is a controller base parameter.
As a preferable scheme of the adaptive backstepping nonlinear swing suppression control method of the double-pendulum rotary crane, the method comprises the following steps: also comprises the following steps of (1) preparing,
the method comprises the following steps of performing tracking control on a double-pendulum rotary crane system by utilizing a reference track of a cantilever, wherein the reference track is a three-section type track:
Figure BDA0003664617360000049
wherein ,ζd5 =30[deg],ζ d6 =45[deg],ζ d1 =15[deg],ζ d2 =30[deg]And q is 5 and 6, the tracking control of 30 degrees is realized in the heave angle within 0 to 5 seconds, the tracking control of 45 degrees is realized in the rotation angle, the tracking control of 15 degrees is realized on the basis of 30 degrees in the heave angle within 5 to 10 seconds, the tracking control of 15 degrees is realized on the basis of 45 degrees in the rotation angle, and the tracking angle is reduced by 30 degrees simultaneously in the last 10 to 15 seconds.
The invention has the beneficial effects that: the design of the self-adaptive control rate has good robustness on system parameters of different rotary cranes; the control method of the PD-like in the self-adaptive controller can realize accurate positioning, and external disturbance is restrained through the feedback of the swing angle information and the design of an error item; the two are combined to realize efficient track tracking and swing suppression of different crane systems, so that a control effect is achieved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without inventive exercise. Wherein:
fig. 1 is an overall flowchart of an adaptive backstepping nonlinear swing suppression control method for a double-pendulum rotary crane according to a first embodiment of the present invention;
FIG. 2 is a schematic diagram of a mathematical model structure of an adaptive backstepping nonlinear swing suppression control method for a double-pendulum rotary crane according to a first embodiment of the present invention;
FIG. 3 is a schematic diagram of simulation results of an adaptive backstepping nonlinear swing suppression control method for a double-pendulum rotary crane according to a first embodiment of the present invention;
fig. 4 is a schematic diagram of a comparative controller btp (bottom projection planning) simulation result of an adaptive back-stepping nonlinear swing suppression control method for a double-pendulum rotary crane according to a second embodiment of the present invention.
Detailed Description
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, specific embodiments accompanied with figures are described in detail below, and it is apparent that the described embodiments are a part of the embodiments of the present invention, not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making creative efforts based on the embodiments of the present invention, shall fall within the protection scope of the present invention.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, but the present invention may be practiced in other ways than those specifically described and will be readily apparent to those of ordinary skill in the art without departing from the spirit of the present invention, and therefore the present invention is not limited to the specific embodiments disclosed below.
Furthermore, reference herein to "one embodiment" or "an embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one implementation of the invention. The appearances of the phrase "in one embodiment" in various places in the specification are not necessarily all referring to the same embodiment, nor are separate or alternative embodiments mutually exclusive of other embodiments.
The present invention will be described in detail with reference to the drawings, wherein the cross-sectional views illustrating the structure of the device are not enlarged partially in general scale for convenience of illustration, and the drawings are only exemplary and should not be construed as limiting the scope of the present invention. In addition, the three-dimensional dimensions of length, width and depth should be included in the actual fabrication.
Meanwhile, in the description of the present invention, it should be noted that the terms "upper, lower, inner and outer" and the like indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, and are only for convenience of describing the present invention and simplifying the description, but do not indicate or imply that the referred device or element must have a specific orientation, be constructed in a specific orientation and operate, and thus, cannot be construed as limiting the present invention. Furthermore, the terms first, second, or third are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
The terms "mounted, connected and connected" in the present invention are to be understood broadly, unless otherwise explicitly specified or limited, for example: can be fixedly connected, detachably connected or integrally connected; they may be mechanically, electrically, or directly connected, or indirectly connected through intervening media, or may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
Example 1
Referring to fig. 1 to 2, an embodiment of the present invention provides a method for controlling adaptive backstepping nonlinear swing suppression of a double-pendulum rotary crane, including:
s1: and establishing a mathematical model of the double-pendulum rotary crane based on a Lagrange kinetic equation and analyzing the characteristics.
Order to
Figure BDA0003664617360000061
Establishing a mathematical model of the double-pendulum rotary crane with distributed mass load:
Figure BDA0003664617360000062
G(q)=[g 1 g z g3g4g5g6lT
N=[OOOOn5n6]T
F f =[f 1 f 2 f 3 f 4 f 5 f 6 ] T
Figure BDA0003664617360000063
q=[ζ 1 ζ 2 ζ 3 ζ 4 ζ 5 ζ 6 ] T
wherein: m (q) is the inertia matrix of the double pendulum rotary crane system,
Figure BDA0003664617360000064
is a centripetal-Coriolis matrix, G (q) is a gravity vector, N is a control input vector, F f Is the mechanical friction force of the double-pendulum rotary crane system, q is the state variable of the double-pendulum rotary crane system,
Figure BDA0003664617360000065
in the form of the first derivative of the signal,
Figure BDA0003664617360000066
is its second derivative; m is h And m l Mass of hook and load, M 0 And M 1 Mass of the cantilever and ballast, respectively, /) 1 And l 2 Respectively the length from the suspension rope and the lifting hook to the center of mass of the load, L 1 L is the length of the ballast and the length of the cantilever, I x 、I y 、I z For the moment of inertia of the cantilever in the x, y, z axes, respectively, I b Is the moment of inertia of the ballast, g is the acceleration of gravity, and ζ is a generalized quantity of state describing a rotary crane system i (i 1.., 4.) is the swing angle of the hook and the load, and ζ is the driving force/torque 5 For cantilever heave direction drive torque, ζ 6 For driving torque in the direction of rotation of the cantilever, f 5 、f 6 Mechanical friction in the direction of cantilever heave and rotation, d i And (1., 4) are air friction parameters.
Further, a friction force feedforward compensation model is established to eliminate the friction force generated by the driving mechanism of the double-pendulum rotary crane, and the friction force feedforward compensation model is expressed as follows:
Figure BDA0003664617360000071
Figure BDA0003664617360000072
wherein ,f51 、f 52 、f 61 、f 62 、ε 1 and ε2 Feedforward compensation model parameters for friction force, f 51 and f61 The value of (A) corresponds to the maximum static friction force, f 52 and f62 Is the coefficient of viscous friction, ε 1 and ε2 Is the static coefficient of friction.
It should be noted that, the friction force feedforward compensation model is used to simply eliminate the friction force effect, so as to effectively avoid the adverse effect of the friction generated by the motion of the driving mechanism on the control effect, for example, the swing caused by the positioning lag/lead is generated to severely increase the difficulty of the anti-swing control.
It should be noted that the inertia matrix of the double pendulum rotary crane system is as follows:
Figure BDA0003664617360000073
Figure BDA0003664617360000076
Figure BDA0003664617360000077
m 13 =m l l 1 l 2 (1+ζ 1 ζ 3 ),
m 14 =m l l 1 l 2 ζ 1 ζ 4
m 15 =(m h +m l )l 1 L(C 51 S 5 ),
m 16 =-(m h +m l )l 1 2 ζ 2 +m l l 1 l 2 ζ 4
m 21 =(m h +m l )l 1 2 ζ 1 ζ 2
Figure BDA0003664617360000074
m 23 =m l l 1 l 2 ζ 2 ζ 3
m 24 =m 1 l 1 l 2 (1+ζ 2 ζ 4 ),
m 25 =-(m h +m l )l 12 S 5
m 26 =(m h +m l )l 1 2 ζ 1 +(m h +m 1 )l 1 LS 5 + m l 1 l 2 ζ 3
m 31 =m l l 1 l 2 (1+ζ 1 ζ 3 ),
m 32 =m l l 1 l 2 ζ 2 ζ 3
Figure BDA0003664617360000075
Figure BDA0003664617360000081
m 35 =m l l 2 L(C 53 S 5 ),
Figure BDA0003664617360000082
m 41 =m l l 1 l 2 ζ 1 ζ 4
m 42 =m l l 1 l 2 (1+ζ 2 ζ 4 ),
Figure BDA0003664617360000083
Figure BDA0003664617360000084
m 45 =-m l l 24 S 5
Figure BDA0003664617360000085
m 51 =(m h +m l )l 1 L(C 51 S 5 ),
m 52 =-(m h +m l )l 12 S 5
m 53 =m l l 2 L(C 53 S 5 ),
m 54 =-m l l 24 S 5
m 55 =((m h +m l )L 2 +I y ),
m 61 =-((m h +ml)l 1 2 ζ 2 +m l l 1 l 2 ζ 4 ),
m 62 =(m h +m l )l 1 LS 5 +m l l 1 l 2 ζ 3 +(m h +m l )l 1 2 ζ 1
Figure BDA0003664617360000086
Figure BDA0003664617360000087
m 65 =-((m h +m l )l 12 C 5 +m l l 24 S 5 ),
Figure BDA0003664617360000088
wherein ,mij The matrix coordinates are represented, i 1, 2 … 6, j 1, 2 … 6.
Centripetal coriolis matrix
Figure BDA0003664617360000089
The following were used:
Figure BDA00036646173600000810
Figure BDA0003664617360000091
Figure BDA0003664617360000092
Figure BDA0003664617360000093
Figure BDA0003664617360000094
Figure BDA0003664617360000095
Figure BDA0003664617360000096
Figure BDA0003664617360000097
Figure BDA0003664617360000098
Figure BDA0003664617360000099
Figure BDA00036646173600000910
Figure BDA00036646173600000911
Figure BDA00036646173600000912
Figure BDA00036646173600000913
Figure BDA00036646173600000914
Figure BDA00036646173600000915
Figure BDA00036646173600000916
Figure BDA00036646173600000917
Figure BDA00036646173600000918
Figure BDA00036646173600000919
Figure BDA00036646173600000920
Figure BDA00036646173600000921
Figure BDA00036646173600000922
Figure BDA00036646173600000923
Figure BDA00036646173600000924
Figure BDA00036646173600000925
Figure BDA00036646173600000926
Figure BDA00036646173600000927
Figure BDA00036646173600000928
c 55 =0,
Figure BDA00036646173600000929
c 61 =0,
c 62 =0,
c 63 =0,
c 64 =0,
Figure BDA00036646173600000930
Figure BDA00036646173600000931
wherein cij The matrix coordinates are represented, i is 1, 2 … 6, and j is 1, 2 … 6.
Furthermore, the characteristics are analyzed according to a mathematical model of the double-pendulum rotary crane.
It should be noted that the characteristics specifically include a swing characteristic (double pendulum system), a load (rod translation characteristic), and the like thereof, and the energy of the entire system can be analyzed by establishing a mathematical model based on the lagrangian equation and analyzing the characteristics thereof.
S2: and establishing an energy function of a mathematical model of the rotary crane according to the characteristics, establishing a Lyapunov function based on a backstepping method and a system control target, and establishing an adaptive controller based on the energy function and the Lyapunov function.
The energy function of the mathematical model of the double-pendulum tower crane is expressed as follows:
Figure BDA0003664617360000101
wherein :
Figure BDA0003664617360000102
respectively representing the cantilever load, hook, swing angle ζ i A speed signal of 1, 4,
Figure BDA0003664617360000103
is the kinetic energy part of a double-pendulum tower crane system,
Figure BDA0003664617360000104
the potential energy portion of its load and hook.
Further, a Lyapunov equation V is established based on a backstepping method and a system control target according to the dynamic rule of a rotary crane model 1
Figure BDA0003664617360000105
wherein :e1 =[α 1 α 2 α 3 α 4 α 5 α 6 ] T ,q d =[ζ 1d ζ 2d ζ 3d ζ 4d ζ 5d ζ 6d ] T ,e 1 =q-q d .e 1 The error of the swing angle and the error of the current position and the target position of the cantilever,
Figure BDA0003664617360000106
the error of the swing angular speed and the error of the rotation speed of the cantilever are obtained;
for the Lyapunov equation V 1 The derivation is carried out to obtain:
Figure BDA0003664617360000107
wherein :k1 Is a positive parameter, e 2 =Θ-Θ d =[β 1 β 2 β 3 β 4 β 5 β 6 ] T In order to make the derivative of the Lyapunov equation semi-negative, design
Figure BDA0003664617360000108
Then there are
Figure BDA0003664617360000109
Figure BDA00036646173600001010
Further, the Lyapunov equation V is designed based on the backstepping method 2 Expressed as:
Figure BDA0003664617360000111
wherein ,
Figure BDA0003664617360000112
to adapt the error values of the column vectors of the system matrix, there are
Figure BDA0003664617360000113
Figure BDA0003664617360000114
Figure BDA0003664617360000115
Is composed of
Figure BDA0003664617360000116
The error value of the adaptive system vector of (1) is
Figure BDA0003664617360000117
S3: and obtaining the self-adaptive control rate based on the Lyapunov function, and adding the self-adaptive control rate into the self-adaptive controller.
Based on Lyapunov equation V 2 And positive determination and semi-negative determination of a derivative of the positive determination are carried out, so that the self-adaptive control rate can be obtained, and the self-adaptive control rate is added into a self-adaptive controller, so that the system forms a controller with stable convergence, and the real-time regulation and control of the swing inhibition of the double-pendulum rotary crane through the self-adaptive controller are realized.
The adaptive controller is represented as:
Figure BDA0003664617360000118
wherein the adaptive control rate is expressed as:
Figure BDA0003664617360000119
further, by adjusting k 1 、k 2 Thereby obtaining a gain of the fast adaptation system, wherein k 1 and k2 Is a controller base parameter.
It should be noted that k in the adaptive controller 1 and k2 Under the condition of giving an initial value, the controller is used for real-time regulation and control, so that better control performance (accurate positioning and quick and effective oscillation elimination) can be obtained in the whole movement process.
In particular, the gain (k) of the PD-like part 1 and k2 ) All gains are positive and can be obtained by trial and error; finally, the correlation parameter f of the feedforward friction model 51 、f 52 、f 61 and f62 After offline recognition, without changing the value of the selection, ε 1 and ε2 Is the static coefficient of friction, which is selected to be 0.01.
Furthermore, the reference track of the fluctuation and rotation of the cantilever is utilized to carry out tracking control on the double-pendulum rotary crane system so as to verify the positioning and pendulum eliminating functions, wherein the reference track is a three-section track represented as:
Figure BDA00036646173600001110
wherein ,ζd5 =30[deg],ζ d6 =45[deg],ζ d1 =15[deg],ζ d2 =30[deg]And q is 5, 6. the tracking control of 30 degrees is realized in 0 to 5 seconds of the heave angle, the tracking control of 45 degrees is realized in the rotation angle, the tracking control of 15 degrees is realized on the basis of 30 degrees in 5 to 10 seconds of the heave angle, the tracking control of 15 degrees is realized on the basis of 45 degrees in the rotation angle, and finally the tracking angle is reduced by 30 degrees in 10 to 15 seconds.
It should be noted that tracking control of the double-pendulum rotary crane system is realized by simulating the pitch angle and the rotation angle of the jib and inputting different reference trajectories for the pitch angle and the rotation angle. The self-adaptive method is verified through simulation software, and the method has good positioning and pendulum-eliminating performance.
Example 2
Referring to fig. 3 to 4, an embodiment of the invention provides a self-adaptive backstepping nonlinear swing suppression control method for a double-pendulum rotary crane, and in order to verify the beneficial effects of the invention, scientific demonstration is performed through simulation experiments.
Compared simulation is carried out by utilizing the controller BTP and the controller using the control method, and the control formula of the controller BTP is as follows:
Figure BDA0003664617360000121
Figure BDA0003664617360000122
for a BTP controller, the state vector q ═ θ 1 θ 2 θ 3 θ 4 θ 5 θ 6 ] T Finally, the gain of the controller is k 5p =5,k 5d =5,k 6p =10,k 6d The simulation results obtained using the BTP controller constructed as described above are shown in table 1 below:
table 1: maximum amplitude experiments compare the results.
Figure BDA0003664617360000123
It can be seen by referring to fig. 2 and 3 that the proposed controller can completely track the target track and realize the positioning function, the BTP controller cannot realize positioning in the rotation and fluctuation directions with positioning errors, the tracking and positioning process of the method is smooth, the positioning task can be completed in 15.50 seconds, for the swing suppression aspect, the amplitude of the hook and the load caused by the controller of the method is not large and cannot exceed 1.95[ deg ], the amplitude of the hook and the load caused by the BTP controller is too large and can reach 2.93[ deg ] at most, and the swing of the method can be completely eliminated within 2-4 seconds after the positioning of the driving mechanism is completed, the suppression effect of the BTP is particularly poor, and the swing is not remained until 20 seconds after multiple violent oscillations, so the method has extremely high swing suppression efficiency, accurate positioning, and no overshoot and steady-state errors.
It should be noted that the above-mentioned embodiments are only for illustrating the technical solutions of the present invention and not for limiting, and although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions may be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, which should be covered by the claims of the present invention.

Claims (10)

1. A self-adaptive backstepping nonlinear swing suppression control method of a double-pendulum rotary crane is characterized by comprising the following steps: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
establishing a mathematical model of the double-pendulum rotary crane based on a Lagrange kinetic equation and analyzing characteristics;
establishing an energy function of a mathematical model of the rotary crane according to the characteristics, establishing a Lyapunov function based on a backstepping method and a system control target, and establishing an adaptive controller based on the energy function and the Lyapunov function;
and obtaining the self-adaptive control rate based on the Lyapunov function, and adding the self-adaptive control rate into the self-adaptive controller.
2. The adaptive backstepping nonlinear oscillation suppression control method of the double-pendulum rotary crane according to claim 1, characterized in that: the mathematical model of the double-pendulum rotary crane is expressed as follows:
order to
Figure FDA0003664617350000011
Figure FDA0003664617350000012
G(q)=[g 1 g 2 g 3 g 4 g 5 g 6 ] T
N=[0 0 0 0 n 5 n 6 ] T
F f =[f 1 f 2 f 3 f 4 f 5 f 6 ] T
Figure FDA0003664617350000013
q=[ζ 1 ζ 2 ζ 3 ζ 4 ζ 5 ζ 6 ] T
Wherein: m (q) is the inertia matrix of the double pendulum rotary crane system,
Figure FDA0003664617350000014
is a centripetal-Coriolis matrix, G (q) is a gravity vector, N is a control input vector, F f Is the mechanical friction force of the double-pendulum rotary crane system, and q is the state variable of the double-pendulum rotary crane system,
Figure FDA0003664617350000015
In the form of the first derivative of the signal,
Figure FDA0003664617350000016
is its second derivative; m is h And m l Mass of hook and load, M 0 And M 1 Mass of the cantilever and ballast, respectively, k 1 And l 2 Respectively the length from the suspension rope and the lifting hook to the center of mass of the load, L 1 L is the length of the ballast and the length of the cantilever, I x 、I y 、I z For the moment of inertia of the cantilever in the x, y, z axes, respectively, I b Is the moment of inertia of the ballast, g is the acceleration of gravity, and ζ is a generalized quantity of state describing a rotary crane system i (i 1.., 4.) is the swing angle of the hook and the load, and ζ is the driving force/torque 5 For cantilever heave direction drive torque, ζ 6 For driving torque in the direction of rotation of the cantilever, f 5 、f 6 Mechanical friction in the direction of cantilever heave and rotation, d l And (1., 4) are air friction parameters.
3. The adaptive backstepping nonlinear oscillation suppression control method of the double-pendulum rotary crane according to claim 2, characterized in that: the mathematical model of the rotary crane further comprises,
establishing a friction feedforward compensation model, wherein the friction feedforward compensation model is expressed as follows:
Figure FDA0003664617350000021
Figure FDA0003664617350000022
wherein ,f51 、f 52 、f 61 、f 62 、ε 1 and ε2 Feedforward compensation model parameters for friction force, f 51 and f61 The value of (A) corresponds to the maximum static friction force, f 52 and f62 Is the coefficient of viscous friction, ε 1 and ε2 Is the static coefficient of friction.
4. The adaptive backstepping nonlinear oscillation suppression control method of the double-pendulum rotary crane according to claim 3, characterized in that: the analysis of the characteristics comprises the analysis of the swing characteristics and the load.
5. The adaptive back-stepping nonlinear swing suppression control method of a double-pendulum rotary crane according to claim 1 or 2, characterized in that: the energy function is expressed as:
Figure FDA0003664617350000023
wherein :
Figure FDA0003664617350000024
respectively representing the cantilever load, hook, swing angle ζ i A speed signal of 1, 4,
Figure FDA0003664617350000025
is the kinetic energy part of a double-pendulum rotary crane system,
Figure FDA0003664617350000026
the potential energy portion of its load and hook.
6. The adaptive backstepping nonlinear oscillation suppression control method of the double-pendulum rotary crane according to claim 5, characterized in that: the establishing of the Lyapunov function includes,
the lyapunov equation V is established based on a backstepping method and a system control target according to the dynamic rule of a rotary crane model 1
Figure FDA0003664617350000027
wherein :e1 =[α 1 α 2 α 3 α 4 α 5 α 6 ] T ,q d =[ζ 1d ζ 2d ζ 3d ζ 4d ζ 5d ζ 6d ] T ,e 1 =q-q d .e 1 The error of the swing angle and the error of the current position and the target position of the cantilever,
Figure FDA0003664617350000028
the error of the swing angular speed and the error of the rotation speed of the cantilever are obtained;
for the Lyapunov equation V 1 The derivation is carried out to obtain:
Figure FDA0003664617350000029
wherein :k1 Is a positive parameter, e 2 =Θ-Θ d =[β 1 β 2 β 3 β 4 β 5 β 6 ] T In order to make the derivative of the Lyapunov equation semi-negative, design
Figure FDA0003664617350000031
Then there are
Figure FDA0003664617350000032
Figure FDA0003664617350000033
Design Lyapunov equation V based on backstepping method 2 Expressed as:
Figure FDA0003664617350000034
wherein ,
Figure FDA0003664617350000035
to adapt the error values of the column vectors of the system matrix, there are
Figure FDA0003664617350000036
Figure FDA0003664617350000037
Figure FDA0003664617350000038
Is composed of
Figure FDA0003664617350000039
The error value of the adaptive system vector of (1) is
Figure FDA00036646173500000310
7. The adaptive backstepping nonlinear oscillation suppression control method of the double-pendulum rotary crane according to claim 6, characterized in that: the adaptive controller is represented as:
Figure FDA00036646173500000311
lyapunov equation V by design 2 And positive determination, wherein the derivative of the positive determination is semi-negative determination, so that a self-adaptive control rate is obtained, and the self-adaptive control rate is added into a self-adaptive controller, so that the system can form a self-adaptive controller with stable convergence.
8. The adaptive backstepping nonlinear oscillation suppression control method of the double-pendulum rotary crane according to claim 7, characterized in that: the adaptive control rate is expressed as:
Figure FDA00036646173500000312
9. the adaptive backstepping nonlinear oscillation suppression control method of the double-pendulum rotary crane according to claim 7, characterized in that: the adaptive controller may further comprise a processor for controlling the adaptive controller,
by adjusting k 1 、k 2 Thereby obtaining a gain that quickly adapts to different systems, where k 1 and k2 Is a controller base parameter.
10. The adaptive backstepping nonlinear oscillation suppression control method of the double-pendulum rotary crane according to claim 7, characterized in that: also comprises the following steps of (1) preparing,
the method comprises the following steps of performing tracking control on a double-pendulum rotary crane system by utilizing a reference track of a cantilever, wherein the reference track is a three-section type track:
Figure FDA00036646173500000313
wherein ,ζd5 =30[deg],ζ d6 =45[deg],ζ d1 =15[deg],ζ d2 =30[deg]And q is 5, the heave angle realizes 30-degree tracking control in 6.0 to 5 seconds, the rotation angle realizes 45-degree tracking control, the heave angle realizes 15-degree tracking control on the basis of 30 degrees in 5 to 10 seconds, the rotation angle realizes 15-degree tracking control on the basis of 45 degrees, and finally the tracking angle is reduced by 30 degrees in 10 to 15 seconds.
CN202210589719.8A 2022-05-26 2022-05-26 Self-adaptive backstepping nonlinear swing inhibition control method for double-swing rotary crane Active CN114890305B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210589719.8A CN114890305B (en) 2022-05-26 2022-05-26 Self-adaptive backstepping nonlinear swing inhibition control method for double-swing rotary crane

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210589719.8A CN114890305B (en) 2022-05-26 2022-05-26 Self-adaptive backstepping nonlinear swing inhibition control method for double-swing rotary crane

Publications (2)

Publication Number Publication Date
CN114890305A true CN114890305A (en) 2022-08-12
CN114890305B CN114890305B (en) 2023-05-12

Family

ID=82726959

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210589719.8A Active CN114890305B (en) 2022-05-26 2022-05-26 Self-adaptive backstepping nonlinear swing inhibition control method for double-swing rotary crane

Country Status (1)

Country Link
CN (1) CN114890305B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116605784A (en) * 2023-05-25 2023-08-18 南京工业大学 Self-adaptive fault-tolerant control method and system for tower crane
CN117105096A (en) * 2023-09-25 2023-11-24 哈尔滨理工大学 Sliding mode control method suitable for rope-length-variable double-swing type ship crane
CN118239386A (en) * 2024-05-29 2024-06-25 山东鲁能特种设备检验检测有限公司 Crane control method, system, medium, equipment and product

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102929150A (en) * 2012-11-13 2013-02-13 湖南航天机电设备与特种材料研究所 Spoiler self-adaptive control method based on discrete control model
CN107678277A (en) * 2017-03-16 2018-02-09 南京工业大学 Sliding mode control method for nonlinear sliding mode surface of double-pendulum bridge crane
CN110980521A (en) * 2019-12-25 2020-04-10 武汉科技大学 System passivity-based under-actuated bridge crane anti-interference anti-swing method
CN111142384A (en) * 2019-12-31 2020-05-12 济南大学 Adaptive neural network tracking control method and system for two-stage pendulum tower crane
CN112875509A (en) * 2021-01-13 2021-06-01 南京工业大学 Positioning and swing eliminating method for under-actuated tower crane with load lifting motion
CN112897338A (en) * 2021-01-13 2021-06-04 南京工业大学 Trajectory tracking and swing suppression control method for under-actuated double-pendulum tower crane
CN113353799A (en) * 2021-06-08 2021-09-07 南京工业大学 Swing suppression control method for double-pendulum tower crane with distributed mass loads
CN113942934A (en) * 2021-11-08 2022-01-18 南开大学 Container bridge crane accurate positioning and anti-swing control method based on speed control

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102929150A (en) * 2012-11-13 2013-02-13 湖南航天机电设备与特种材料研究所 Spoiler self-adaptive control method based on discrete control model
CN107678277A (en) * 2017-03-16 2018-02-09 南京工业大学 Sliding mode control method for nonlinear sliding mode surface of double-pendulum bridge crane
CN110980521A (en) * 2019-12-25 2020-04-10 武汉科技大学 System passivity-based under-actuated bridge crane anti-interference anti-swing method
CN111142384A (en) * 2019-12-31 2020-05-12 济南大学 Adaptive neural network tracking control method and system for two-stage pendulum tower crane
CN112875509A (en) * 2021-01-13 2021-06-01 南京工业大学 Positioning and swing eliminating method for under-actuated tower crane with load lifting motion
CN112897338A (en) * 2021-01-13 2021-06-04 南京工业大学 Trajectory tracking and swing suppression control method for under-actuated double-pendulum tower crane
CN113353799A (en) * 2021-06-08 2021-09-07 南京工业大学 Swing suppression control method for double-pendulum tower crane with distributed mass loads
CN113942934A (en) * 2021-11-08 2022-01-18 南开大学 Container bridge crane accurate positioning and anti-swing control method based on speed control

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116605784A (en) * 2023-05-25 2023-08-18 南京工业大学 Self-adaptive fault-tolerant control method and system for tower crane
CN116605784B (en) * 2023-05-25 2023-10-24 南京工业大学 Self-adaptive fault-tolerant control method and system for tower crane
CN117105096A (en) * 2023-09-25 2023-11-24 哈尔滨理工大学 Sliding mode control method suitable for rope-length-variable double-swing type ship crane
CN117105096B (en) * 2023-09-25 2024-02-23 哈尔滨理工大学 Sliding mode control method suitable for rope-length-variable double-swing type ship crane
CN118239386A (en) * 2024-05-29 2024-06-25 山东鲁能特种设备检验检测有限公司 Crane control method, system, medium, equipment and product

Also Published As

Publication number Publication date
CN114890305B (en) 2023-05-12

Similar Documents

Publication Publication Date Title
CN114890305A (en) Self-adaptive backstepping nonlinear swing suppression control method of double-pendulum rotary crane
US8839967B2 (en) Crane for handling a load hanging on a load cable
CN113353799A (en) Swing suppression control method for double-pendulum tower crane with distributed mass loads
Yang et al. Motion trajectory-based transportation control for 3-D boom cranes: Analysis, design, and experiments
CN107942670A (en) A kind of double-flexibility space manipulator Fuzzy Robust Controller sliding formwork, which is cut, trembles motion control method
CN112875509A (en) Positioning and swing eliminating method for under-actuated tower crane with load lifting motion
CN112897338B (en) Under-actuated double-pendulum tower crane track tracking and swing inhibition control method
US11029650B2 (en) Machine learning device, control system, and machine learning method
CN110376902A (en) A kind of design method of Underactuated Mechanical Systems Servo Restriction tracking control unit
CN109018237A (en) A kind of ship fin/water tank joint rollstabilization intelligence control system and its control method
CN114890314A (en) Fault-tolerant control method for double-pendulum tower crane with online track correction
Cunha et al. Dynamic modeling and stability analysis of model-scale helicopters with bell-hiller stabilizing bar
Kumar et al. An adaptive particle swarm optimization algorithm for robust trajectory tracking of a class of under actuated system
Kim et al. Oscillation damping control of pendulum-like manipulation platform using moving masses
US11243501B2 (en) Machine learning device, control system, and machine learning
CN110134135B (en) Four-rotor aircraft control method based on improved MPC-PID
CN114879504B (en) Self-adaptive nonlinear control method of four-degree-of-freedom marine rotary crane
Roy Robust backstepping control for small helicopter
CN113110023B (en) Structure correction control method based on Hamiltonian model of diesel engine
Wai et al. Development of adaptive sliding-mode control for nonlinear dual-axis inverted-pendulum system
WO2010079801A1 (en) Jib crane control method and device
Tum et al. Swing-up control of a single inverted pendulum on a cart with input and output constraints
CN107728470B (en) Mechanical arm attitude control method without speed measurement
Monir Analyzing and designing control system for an inverted pendulum on a cart
CN110254741A (en) A kind of design method of flight control system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant