CN114735140B - Method, equipment and medium for compensating disturbance speed of wind power pile boarding trestle - Google Patents

Method, equipment and medium for compensating disturbance speed of wind power pile boarding trestle Download PDF

Info

Publication number
CN114735140B
CN114735140B CN202210380956.3A CN202210380956A CN114735140B CN 114735140 B CN114735140 B CN 114735140B CN 202210380956 A CN202210380956 A CN 202210380956A CN 114735140 B CN114735140 B CN 114735140B
Authority
CN
China
Prior art keywords
increment
value
interference
formula
control
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.)
Active
Application number
CN202210380956.3A
Other languages
Chinese (zh)
Other versions
CN114735140A (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.)
Harbin Harbin Ship Rolling Reduction Automation Equipment Co ltd
Harbin Engineering University
Original Assignee
Harbin Harbin Ship Rolling Reduction Automation Equipment Co ltd
Harbin Engineering 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 Harbin Harbin Ship Rolling Reduction Automation Equipment Co ltd, Harbin Engineering University filed Critical Harbin Harbin Ship Rolling Reduction Automation Equipment Co ltd
Priority to CN202210380956.3A priority Critical patent/CN114735140B/en
Publication of CN114735140A publication Critical patent/CN114735140A/en
Application granted granted Critical
Publication of CN114735140B publication Critical patent/CN114735140B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B63SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63BSHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B27/00Arrangement of ship-based loading or unloading equipment for cargo or passengers
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B63SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63BSHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B17/00Vessels parts, details, or accessories, not otherwise provided for
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B63SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
    • B63BSHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING 
    • B63B17/00Vessels parts, details, or accessories, not otherwise provided for
    • B63B2017/0072Seaway compensators
    • 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
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind turbines with rotation axis in wind direction
    • 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
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/727Offshore wind turbines

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • Ocean & Marine Engineering (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Feedback Control In General (AREA)

Abstract

The invention provides a method, equipment and medium for compensating the interference speed of a wind power pile boarding trestle. According to the method, the interference speed (ship motion speed) compensation mode is introduced and is combined with the interference position compensation, so that the rapidity and the compensation effect of the active motion compensation trestle are improved, the use experience of the trestle lapping process is improved, the accessibility of the wind power pile in severe weather is improved, and the safety of passengers is improved.

Description

Method, equipment and medium for compensating interference speed of wind power pile boarding trestle
Technical Field
The invention belongs to the technical field of offshore wind power boarding trestles with a wave compensation function, and particularly relates to a method, equipment and medium for compensating the interference speed of a wind power pile boarding trestle. In particular to a control method for compensating the movement speed of a ship to improve the movement compensation precision and realize the smooth lap joint of an offshore wind power pile on a trestle under the condition of high sea.
Background
Wind energy is a clean renewable energy source. Compared with land wind power generation, offshore wind power has the advantages of no land occupation, small influence on birds and the like, and in addition, offshore wind power also has the advantages of less calm period, large wind power and the like. China has abundant wind energy in offshore and middle and far seas, and the development of offshore wind power has important significance for optimizing energy structures in China.
However, offshore wind power is more expensive to install and maintain than terrestrial wind power. Under severe weather conditions, the severe shaking of the operation and maintenance ship makes it difficult and dangerous for operation and maintenance personnel to board the wind power pile. An Active Motion Compensation (AMC) trestle enables the tail end of the trestle to keep static relative to a geodetic coordinate system (a wind power pile fixed on the ground) or reduces the influence of sea waves on the movement of the trestle by compensating the movement of a ship with a plurality of degrees of freedom, so that the accessibility of personnel to the wind power pile in severe weather is improved.
The three-degree-of-freedom AMC trestle is mainly composed of a base (B1), a support (B2), a pitching bridge (B3), a telescopic bridge (B4) and an inserting plate (B5) positioned at the tail end of the telescopic bridge, as shown in figure 1. When the inserting plate is inserted into the basic ladder stand of the wind turbine generator platform, maintenance personnel can reach the ladder stand through the support, the pitching bridge and the telescopic bridge to climb to the offshore wind turbine generator platform. According to the requirement of safety regulations, when the distance between the operation and maintenance ship and the ladder stand of the offshore wind turbine generator set does not exceed half of the step length of an adult, the ladder stand of the wind turbine generator set can be used for climbing. However, the influence of sea tides and waves can cause the ladder of the ship and the offshore wind power platform to have a certain distance in the vertical direction and the horizontal direction, and the distance can change along with the shaking of the ship.
The motion measuring unit arranged on the ship can measure the motion of the ship with six degrees of freedom, and the yawing, pitching and heaving motions of the ship are compensated through the (Y1) rotary hydraulic cylinder, the (Y2) pitching hydraulic cylinder and the (Y3) telescopic hydraulic cylinder, so that the displacement change of the (B5) inserting plate relative to the ladder stand is caused. Meanwhile, the whole trestle can also keep a relatively stable state.
When the sea stormy waves are small, an operator can close the active motion compensation function, the proportional valve group is controlled to control the rotation speed, the pitching speed and the motion speed of the piston rod of the telescopic hydraulic cylinder through the operating handle, and the inserting plate at the tail end of the telescopic bridge can be smoothly inserted into the ladder stand to complete the establishment of a personnel channel. In a heavy storm state, the swing of the ship can be amplified by the inserting plate positioned at the tail end of the telescopic bridge, so that the process of inserting the ladder becomes difficult. At this time, it is necessary to turn on the active compensation function, to eliminate or reduce the influence of the roll of the vessel on the insertion plate, and when the operation handle is not operated, the coordinates of the end of the insertion plate are kept stationary with respect to the geodetic coordinate system or the influence of the motion of the vessel is reduced. In order to complete active Motion compensation, motion of a ship in six degrees of freedom needs to be measured by a Motion Reference Unit (MRU), and compensation quantities of a trestle in a rotation degree of freedom, a pitching degree of freedom and a stretching degree of freedom are calculated through a conversion matrix so as to compensate motions of the ship in the directions of yawing, pitching and heaving. In order to keep the insertion plate static relative to a geodetic coordinate system, the active motion compensation trestle needs to compensate the position change of the insertion plate caused by the movement of a ship, so the trestle is controlled by a position servo system in three degrees of freedom of rotation, pitching and stretching, and the extension of a piston rod of a hydraulic cylinder is controlled according to position errors. However, the following problems mainly exist in the active motion compensation process in this compensation method:
(1) Due to the influences of ship motion randomness, MRU sampling delay, hydraulic cylinder motion speed limitation, mechanical clearance of trestle installation and the like, in the working process of the active motion compensation trestle, for relatively high-frequency ship motion, only the measured ship position is used for compensation, tracking delay can be caused, and the compensation effect is poor.
(2) In the process that the operating handle drives the insertion plate to be inserted into the basic ladder of the wind turbine platform under the high sea condition, only the speed control of the operating handle and the position control of active motion compensation are simply superposed, and the experience consistent with the process of the operating handle when the compensation function is closed under the calm sea cannot be obtained, so that the operation is not smooth.
Disclosure of Invention
The invention aims to solve the problems in the prior art and provides a method, equipment and medium for compensating the disturbance speed of a wind power pile boarding trestle. According to the method, the interference speed (ship motion speed) compensation mode is introduced and is combined with the interference position compensation, so that the rapidity and the compensation effect of the active motion compensation trestle are improved, the use experience of the trestle lapping process is improved, the accessibility of the wind power pile in severe weather is improved, and the safety of passengers is improved.
The invention is realized by the following technical scheme, and provides an interference speed compensation method for a wind power pile boarding trestle, which specifically comprises the following steps:
step 1, setting the position of an insertion plate as x = [ x ] 1 ,x 2 ,x 3 ] T ,x 1 ,x 2 ,x 3 The coordinate values of the insert plate in three coordinate axes of the geodetic coordinate system O-XYZ are respectively expressed as follows:
Figure BDA0003592964980000021
wherein A, B, C and D are system matrices, u (k) = [ r (k), p (k), s (k)] T Respectively a rotation value, a pitching value and a stretching value of the trestle; d (k) = [ r s (k),p s (k),h s (k)] T Respectively the ship's yaw value, pitch value and heave value, and the output value y = [ x = [ [ x ] 1 ,x 2 ,x 3 ] T K represents a time series;
step 2, changing the equation (1) into a differential form, setting a target function J (k) for wrapping and changing the position and speed tracking error, and determining the increment of the controlled variable through the target function;
step 3, obtaining the increment of the ship motion disturbance through an Extended State Observer (ESO);
step 4, synthesizing the control quantity increment and the interference quantity increment to obtain a final control quantity;
step 5, recording the system state at the moment when the AMC function is started as x (0), and inserting the target position x of the board when the operating handle is not operated r (k) = x (0), target speed of insertion plate
Figure BDA0003592964980000031
The control quantity at this time is used for compensating the change of the coordinates of the inserting plate caused by the movement of the ship; when in useTarget speed of insert plate when operating handle is actuated
Figure BDA0003592964980000032
Proportional to the output voltage of the operating handle, x r (k) Is composed of
Figure BDA0003592964980000033
To track the operation of the handle.
Further, equation (1) is changed to differential form:
Figure BDA0003592964980000034
in the formula (I), the compound is shown in the specification,
Figure BDA0003592964980000035
respectively increment of the state quantity, the control quantity and the interference quantity; in the system matrix
Figure BDA0003592964980000036
Figure BDA0003592964980000037
O is a zero matrix and I is an identity matrix.
Further, the increments of the control variables are obtained by solving the following objective function, taking into account the magnitude constraint and the velocity constraint of the control variables u (k):
Figure BDA0003592964980000038
Figure BDA0003592964980000039
in the formula, N p For predicting the time domain, matrices Q and S are semi-positive definite matrices, matrices R and P are positive definite matrices,
Figure BDA00035929649800000310
Figure BDA00035929649800000311
representing a state error, x r (k) For inserting into the target position of the plate, x d (k) The differential of x (k), i.e., the velocity of the insert plate, thus the error includes both position and velocity errors;
Figure BDA00035929649800000312
increased degrees of freedom for controlled quantity increments
Figure BDA00035929649800000313
To simulate the increase in the amount of interference.
Further, minimizing the objective function J (k) yields the following control sequence:
Figure BDA00035929649800000314
therefore, solving the objective function problem at each sampling period can be converted to a quadratic programming as follows:
Figure BDA00035929649800000315
Subject to Lη≤b
wherein L and b are used to limit the magnitude and speed of the control quantity; the matrices F and H are defined as follows and are obtained by real-time calculation:
Figure BDA0003592964980000041
Figure BDA0003592964980000042
Figure BDA0003592964980000043
Figure BDA0003592964980000044
status of state
Figure BDA0003592964980000045
Can be expressed as:
Figure BDA0003592964980000046
further, formula (5) is adopted to optimize eta in real time * An optimized value of the control quantity increment can be obtained
Figure BDA0003592964980000047
At the same time, optimized interference value increment can be obtained
Figure BDA0003592964980000048
The control amount increment can thus be obtained by the following equation:
Figure BDA0003592964980000049
in which mu (k) is used to compensate for interference
Figure BDA00035929649800000410
Which pass through
Figure BDA00035929649800000411
Obtaining; minimize error E (k), μ (k) satisfies
Figure BDA00035929649800000412
In the formula (I), the compound is shown in the specification,
Figure BDA00035929649800000413
is an estimate of the interference, while mu (k) is onCalculated by the following formula:
Figure BDA00035929649800000414
further, in order to obtain the interference increment
Figure BDA00035929649800000415
Is estimated by
Figure BDA00035929649800000416
The disturbance is observed using the extended state observer ESO as follows:
Figure BDA0003592964980000051
in the formula (I), the compound is shown in the specification,
Figure BDA0003592964980000052
an estimate of the value of the symbol x,
Figure BDA0003592964980000053
representing the estimated value of the differential of x,
Figure BDA0003592964980000054
represents an estimate of x second derivative; t is s For the sampling period, θ is the observer gain, k o1 ,k o2 ,k o3 And k o4 Is the observer coefficient; function g i (i =1,2,3,4) is expressed as:
Figure BDA0003592964980000055
in the formula, alpha 1 =γ,α 2 =2γ-1,α 3 =3γ-2,α 4 4 γ -3, γ ∈ (3/4,1); k is a positive definite matrix; based on the estimation of d (k) by equation (10), the
Figure BDA0003592964980000056
And calculating μ (k) by equation (9); and then the optimized control variable is
Figure BDA0003592964980000057
Further, the maximum and minimum values of the control quantity are respectively determined to be u through the stroke of the hydraulic cylinder and the mechanical structure sU 、u sL (ii) a The maximum value and the minimum value of the increment of the control quantity are respectively determined by the flow of the proportional valve
Figure BDA0003592964980000058
Namely, it is
Figure BDA0003592964980000059
μ (k) in equation (9) may be re-described as:
Figure BDA00035929649800000510
in the formula (I), the compound is shown in the specification,
Figure BDA00035929649800000511
the magnitude constraint of the control quantity can therefore be expressed by the following inequality:
Figure BDA00035929649800000512
in the formula (I), the compound is shown in the specification,
Figure BDA0003592964980000061
also, the constraint on the control quantity increment can be expressed by the following inequality:
Figure BDA0003592964980000062
in the formula (I), the compound is shown in the specification,
Figure BDA0003592964980000063
and
Figure BDA0003592964980000064
are respectively as
Figure BDA0003592964980000065
And
Figure BDA0003592964980000066
an increment of (d);
Figure BDA0003592964980000067
the invention provides electronic equipment, which comprises a memory and a processor, wherein the memory stores a computer program, and the processor realizes the steps of the interference speed compensation method for the wind power pile boarding trestle when executing the computer program.
The invention provides a computer-readable storage medium for storing computer instructions, which when executed by a processor implement the steps of the disturbance speed compensation method for a wind power pile boarding trestle.
The invention has the beneficial effects that:
(1) The difference form of the trestle model is adopted, and the tracking speed error is added into the error, so that the rapidity and the dynamic compensation precision of the active motion compensation are improved;
(2) A value proportional to the output voltage of the operating handle is used as a tracking speed, and the integral of the tracking speed is used as a target position, so that the operating experience of the handle is improved;
(3) The speed limit of the movement of the hydraulic cylinder and the position constraint of the trestle are considered in the control quantity solving control, so that the obtained control quantity is physically feasible control, the emergency brake caused by limiting the movement amplitude by a limit switch is avoided, and the movement stability at the limit position is ensured.
Drawings
FIG. 1 is a schematic diagram of an AMC trestle; wherein Y1 is a rotary hydraulic cylinder; a Y2 pitching hydraulic cylinder; a Y3 telescopic hydraulic cylinder; b1, a base; b2, a support; b3, a pitching bridge; b4, a telescopic bridge; b5 insert plate;
FIG. 2 is a schematic diagram of a trestle lapped wind power pile; wherein B4 telescopic bridge; b5 insert plate; c1, a wind turbine platform; c2 Climbing a ladder;
fig. 3 is a functional block diagram of disturbance speed compensation.
Detailed Description
The technical solutions in the embodiments of the present invention will be described below clearly and completely with reference to the accompanying drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
With reference to fig. 1 to 3, the invention provides a disturbance speed compensation method for a wind power pile boarding trestle, wherein the lapping process of the trestle is as shown in fig. 2, and an insertion plate is inserted into a basic ladder of a wind turbine platform to complete the establishment of a personnel channel; the method specifically comprises the following steps:
step 1, setting the position of an insertion plate as x = [ x ] 1 ,x 2 ,x 3 ] T ,x 1 ,x 2 ,x 3 The coordinate values of the insert plate in three coordinate axes of the geodetic coordinate system O-XYZ are respectively expressed as follows:
Figure BDA0003592964980000071
wherein A, B, C and D are system matrices, u (k) = [ r (k), p (k), s (k)] T Respectively a rotation value, a pitching value and a stretching value of the trestle; d (k) = [ r = s (k),p s (k),h s (k)] T Respectively yaw value, pitch value and heave value of the ship, and an output value y = [ x = 1 ,x 2 ,x 3 ] T K represents a time series;
step 2, in order to process the rapidly changing ship motion, changing the equation (1) into a difference form, setting a target function J (k) for wrapping and changing position and speed tracking errors, and determining the increment of the control quantity through the target function;
step 3, obtaining the increment of the ship motion disturbance through an Extended State Observer (ESO);
step 4, synthesizing the control quantity increment and the interference quantity increment to obtain a final control quantity;
step 5, recording the system state at the moment when the AMC function is started as x (0), and inserting the target position x of the board when the operating handle is not operated r (k) = x (0), target speed of insertion plate
Figure BDA0003592964980000072
The control quantity at the moment is used for compensating the change of the coordinates of the inserting plate caused by the movement of the ship, and the error comprises a speed error; inserting the target speed of the plate when the operating handle is actuated
Figure BDA0003592964980000073
Proportional to the output voltage of the operating handle, x r (k) Is composed of
Figure BDA0003592964980000074
To track the operation of the handle. Target position of the insert plate in equation (3) when the handle is operated
Figure BDA0003592964980000075
(integrated values are replaced by accumulated sums in a digitally controlled system), target velocity of the end of the plate is interpolated
Figure BDA0003592964980000076
(v r Output voltage k for three degrees of freedom of the operating handle v Is a scale factor). The control amount at this time is used to compensate for the change in the coordinates of the insert plate due to the movement of the vessel, and is also responsive to the operation of the handle. By adjusting k v Is obtained by v r Directly control the speed of the proportional valve opening to make the handle operate more smoothly and improve the function of opening the AMCThe operation experience of the rear handle.
Equation (1) is changed to differential form:
Figure BDA0003592964980000081
in the formula (I), the compound is shown in the specification,
Figure BDA0003592964980000082
respectively increment of the state quantity, the control quantity and the interference quantity; in the system matrix
Figure BDA0003592964980000083
Figure BDA0003592964980000084
O is a zero matrix and I is an identity matrix.
The increments of the control variables are obtained by solving the following objective function, taking into account the magnitude constraint and the velocity constraint of the control variables u (k):
Figure BDA0003592964980000085
Figure BDA0003592964980000086
in the formula, N p For predicting the time domain, matrices Q and S are semi-positive definite matrices, matrices R and P are positive definite matrices,
Figure BDA0003592964980000087
Figure BDA0003592964980000088
representing a state error, x r (k) For inserting into the target position of the plate, x d (k) The differential of x (k) is the speed of the insert plate, so the error includes both position error and speed error;
Figure BDA0003592964980000089
increased degree of freedom for controlled quantity increments
Figure BDA00035929649800000810
To simulate the increase in the amount of interference.
The objective function contains the speed error and the increment of the disturbance quantity (equivalent to the disturbance speed), so that the rapidly changing ship motion can be reflected.
Minimizing the objective function J (k) yields the following control sequence:
Figure BDA00035929649800000811
therefore, solving the objective function problem at each sampling period can be converted to a quadratic programming as follows:
Figure BDA00035929649800000812
Subject to Lη≤b
wherein L and b are used to limit the magnitude and speed of the control quantity; the matrices F and H are defined as follows and are obtained by real-time calculation:
Figure BDA0003592964980000091
Figure BDA0003592964980000092
Figure BDA0003592964980000093
Figure BDA0003592964980000094
state of state
Figure BDA0003592964980000095
Can be expressed as:
Figure BDA0003592964980000096
real-time optimization of eta using equation (5) * An optimum value of the control amount increment can be obtained
Figure BDA0003592964980000097
At the same time, optimized interference value increment can be obtained
Figure BDA0003592964980000098
The control amount increment can thus be obtained by the following equation:
Figure BDA0003592964980000099
in which mu (k) is used to compensate for interference
Figure BDA00035929649800000910
Which pass through
Figure BDA00035929649800000911
Obtaining; minimizing the error E (k), μ (k) satisfies
Figure BDA00035929649800000912
In the formula (I), the compound is shown in the specification,
Figure BDA00035929649800000913
is an estimate of interference, while μ (k) is calculated by:
Figure BDA00035929649800000914
to obtain an increase in the amount of interference
Figure BDA00035929649800000915
Is estimated value of
Figure BDA00035929649800000916
The disturbance is observed using the extended state observer ESO as follows:
Figure BDA0003592964980000101
in the formula (I), the compound is shown in the specification,
Figure BDA0003592964980000102
an estimate of the value of x is represented,
Figure BDA0003592964980000103
represents an estimate of the differential of x,
Figure BDA0003592964980000104
an estimate representing the x second derivative; t is a unit of s For the sampling period, θ is the observer gain, k o1 ,k o2 ,k o3 And k o4 Is the observer coefficient; function g i (i =1,2,3,4) is expressed as:
Figure BDA0003592964980000105
in the formula, alpha 1 =γ,α 2 =2γ-1,α 3 =3γ-2,α 4 =4 γ -3, γ ∈ (3/4,1); k is a positive definite matrix; based on the estimation of d (k) by equation (10), the
Figure BDA0003592964980000106
And calculating μ (k) by equation (9); and then the optimized control variable is
Figure BDA0003592964980000107
Passing liquidStroke and mechanical structure of pressure cylinder, and maximum and minimum values of control quantity are respectively determined as u sU 、u sL (ii) a The maximum value and the minimum value of the control quantity increment are respectively determined by the flow of the proportional valve
Figure BDA0003592964980000108
Namely that
Figure BDA0003592964980000109
μ (k) in equation (9) can be re-described as:
Figure BDA00035929649800001010
in the formula (I), the compound is shown in the specification,
Figure BDA00035929649800001011
the magnitude constraint of the control quantity can therefore be expressed by the following inequality:
Figure BDA00035929649800001012
in the formula (I), the compound is shown in the specification,
Figure BDA0003592964980000111
also, the constraint on the control quantity increment can be expressed by the following inequality:
Figure BDA0003592964980000112
in the formula (I), the compound is shown in the specification,
Figure BDA0003592964980000113
and
Figure BDA0003592964980000114
are respectively as
Figure BDA0003592964980000115
And
Figure BDA0003592964980000116
the increment of (d);
Figure BDA0003592964980000117
as shown in fig. 3, the system of the active motion compensation method based on disturbance velocity compensation of the present invention includes an operation handle integration unit, a quadratic programming unit, a control quantity calculation unit, an extended state observer unit, a controlled object of a trestle, and an MRU.
The invention provides electronic equipment which comprises a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to realize the steps of the interference speed compensation method for the wind power pile boarding trestle.
The invention provides a computer-readable storage medium for storing computer instructions, which when executed by a processor, implement the steps of the disturbance speed compensation method for a wind power pile boarding trestle.
The method, the device and the medium for compensating the disturbance speed of the wind power pile boarding trestle are introduced in detail, a specific example is applied in the method for explaining the principle and the implementation mode of the method, and the description of the embodiment is only used for helping to understand the method and the core idea of the method; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.

Claims (4)

1. The method for compensating the disturbance speed of the wind power pile boarding trestle is characterized by comprising the following steps of:
step 1, setting the position of the insertion plateIs x = [) 1 ,x 2 ,x 3 ] T ,x 1 ,x 2 ,x 3 The coordinate values of the insert plate on three coordinate axes of the geodetic coordinate system O-XYZ are respectively shown as follows:
Figure FDA0003985546890000011
wherein A, B, C and D are system matrices, u (k) = [ r (k), p (k), s (k)] T Respectively a rotation value, a pitching value and a stretching value of the trestle; d (k) = [ r = s (k),p s (k),h s (k)] T Respectively a yawing value, a pitching value and a heaving value of the ship, wherein the output value is y (k), and k represents a time sequence; u (k) is a control variable, and d (k) is an interference quantity;
step 2, changing the equation (1) into a differential form, setting a target function J (k) containing position and speed tracking errors, and determining the increment of the controlled variable through the target function;
step 3, obtaining the increment of the ship motion disturbance through an Extended State Observer (ESO);
step 4, synthesizing the control quantity increment and the interference quantity increment to obtain a final control quantity;
step 5, recording the system state of the moment when the active motion compensation AMC function is started as x (0), and when the operating handle is not operated, inserting the target position x of the board r (k) = x (0), target speed of insertion plate
Figure FDA0003985546890000012
The control quantity at this time is used for compensating the change of the coordinates of the inserting plate caused by the movement of the ship; inserting the target speed of the plate when the operating handle is actuated
Figure FDA0003985546890000013
Proportional to the output voltage of the operating handle, x r (k) Is composed of
Figure FDA0003985546890000014
Is integrated withOperation of the tracking handle;
changing equation (1) to differential form:
Figure FDA0003985546890000015
in the formula (I), the compound is shown in the specification,
Figure FDA0003985546890000016
respectively increment of the state quantity, the control quantity and the interference quantity; in the system matrix
Figure FDA0003985546890000017
Figure FDA0003985546890000018
O is a zero matrix, and I is an identity matrix;
the increments of the control variables are obtained by solving the following objective function, taking into account the magnitude constraint and the velocity constraint of the control variables u (k):
Figure FDA0003985546890000021
in the formula, N p For predicting the time domain, matrices Q and S are semi-positive definite matrices, matrices R and P are positive definite matrices,
Figure FDA0003985546890000022
Figure FDA0003985546890000023
representing a state error, x r (k) For inserting into the target position of the plate, x d (k) The differential of x (k), i.e., the velocity of the insert plate, thus the error includes both position and velocity tracking errors;
Figure FDA0003985546890000024
increased degrees of freedom for controlled quantity increments
Figure FDA0003985546890000025
To model the increment of the interference amount;
minimizing the objective function J (k) yields the following control sequence:
Figure FDA0003985546890000026
therefore, solving the objective function problem at each sampling period translates into quadratic programming as follows:
Figure FDA0003985546890000027
wherein L and b are used to limit the magnitude and speed of the control quantity; the matrices F and H are defined as follows and obtained by real-time calculation:
Figure FDA0003985546890000028
Figure FDA0003985546890000029
Figure FDA00039855468900000210
Figure FDA00039855468900000211
state of state
Figure FDA00039855468900000212
Expressed as:
Figure FDA00039855468900000213
real-time optimization of eta using equation (5) * Obtaining an optimized value of the control quantity increment
Figure FDA0003985546890000031
At the same time, an optimized interference value increment is obtained
Figure FDA00039855468900000317
The control amount increment is thus obtained by the following formula:
Figure FDA0003985546890000032
in which mu (k) is used to compensate for interference
Figure FDA0003985546890000033
Which is passed through
Figure FDA0003985546890000034
Obtaining; minimize error E (k), μ (k) satisfies
Figure FDA0003985546890000035
In the formula (I), the compound is shown in the specification,
Figure FDA0003985546890000036
is an estimate of interference, while μ (k) is calculated by:
Figure FDA0003985546890000037
to obtain an increase in the amount of interference
Figure FDA0003985546890000038
Is estimated value of
Figure FDA0003985546890000039
The disturbance is observed using the extended state observer ESO as follows:
Figure FDA00039855468900000310
in the formula (I), the compound is shown in the specification,
Figure FDA00039855468900000311
an estimate of the value of the symbol x,
Figure FDA00039855468900000312
representing the estimated value of the differential of x,
Figure FDA00039855468900000313
an estimate representing the x second derivative; t is s For the sampling period, θ is the observer gain, k o1 ,k o2 ,k o3 And k o4 Is the observer coefficient; function g i (i =1,2,3,4) is expressed as:
Figure FDA00039855468900000314
in the formula, alpha 1 =γ,α 2 =2γ-1,α 3 =3γ-2,α 4 4 γ -3, γ ∈ (3/4,1); k is a positive definite matrix; based on the estimation of d (k) by equation (10), the
Figure FDA00039855468900000315
And calculating μ (k) by equation (9); the control variable to be optimized is
Figure FDA00039855468900000316
2. Method according to claim 1, characterized in that the maximum and minimum values of the controlled variable are determined as u, respectively, by the stroke of the hydraulic cylinder and the mechanical structure sU 、u sL (ii) a The maximum value and the minimum value of the control quantity increment are respectively determined by the flow of the proportional valve
Figure FDA0003985546890000041
Namely that
Figure FDA0003985546890000042
μ (k) in equation (9) is re-described as:
Figure FDA0003985546890000043
in the formula (I), the compound is shown in the specification,
Figure FDA0003985546890000044
the amplitude constraint of the control quantity is therefore expressed by the following inequality:
Figure FDA0003985546890000045
in the formula (I), the compound is shown in the specification,
Figure FDA0003985546890000046
also, the constraint of the control quantity increment is expressed by the following inequality:
Figure FDA0003985546890000047
in the formula (I), the compound is shown in the specification,
Figure FDA0003985546890000048
and
Figure FDA0003985546890000049
are respectively as
Figure FDA00039855468900000410
And
Figure FDA00039855468900000411
an increment of (d);
Figure FDA00039855468900000412
3. an electronic device comprising a memory and a processor, the memory storing a computer program, wherein the processor, when executing the computer program, performs the steps of the method according to any of claims 1-2.
4. A computer-readable storage medium storing computer instructions, which when executed by a processor, perform the steps of the method of any one of claims 1-2.
CN202210380956.3A 2022-04-12 2022-04-12 Method, equipment and medium for compensating disturbance speed of wind power pile boarding trestle Active CN114735140B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210380956.3A CN114735140B (en) 2022-04-12 2022-04-12 Method, equipment and medium for compensating disturbance speed of wind power pile boarding trestle

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210380956.3A CN114735140B (en) 2022-04-12 2022-04-12 Method, equipment and medium for compensating disturbance speed of wind power pile boarding trestle

Publications (2)

Publication Number Publication Date
CN114735140A CN114735140A (en) 2022-07-12
CN114735140B true CN114735140B (en) 2023-03-10

Family

ID=82282201

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210380956.3A Active CN114735140B (en) 2022-04-12 2022-04-12 Method, equipment and medium for compensating disturbance speed of wind power pile boarding trestle

Country Status (1)

Country Link
CN (1) CN114735140B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116184841B (en) * 2023-04-28 2023-07-21 山东大学 Marine trestle model predictive control method based on extremum searching algorithm
CN117163219B (en) * 2023-09-15 2024-05-14 哈尔滨理工大学 Shipborne trestle feedforward fuzzy control method considering constraint between long rods

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014131884A (en) * 2013-01-04 2014-07-17 Shimonoseki Ryoju Engineering Kk Embarking-disembarking device
CN105408199A (en) * 2013-07-19 2016-03-16 Icd软件股份公司 Apparatus and method for providing active motion compensation control of an articulated gangway
WO2016193713A1 (en) * 2015-06-02 2016-12-08 Marine Electrical Consulting Limited Method and apparatus for adaptive motion compensation
CN107434010A (en) * 2017-09-26 2017-12-05 哈尔滨工程大学 A kind of electronic wave Active Compensation, which is stepped on, multiplies system and its control method
CN107675607A (en) * 2016-10-28 2018-02-09 福建省新能海上风电研发中心有限公司 A kind of six degree of freedom Active Compensation formula offshore platform steps on the application method for multiplying trestle
CN111045332A (en) * 2019-12-27 2020-04-21 哈尔滨工程大学 Unmanned ship path tracking guidance strategy and disturbance compensation method
CN113031429A (en) * 2021-02-26 2021-06-25 北京星光凯明智能科技有限公司 Shipborne carrier stabilizing platform and control method
CN113104153A (en) * 2021-04-25 2021-07-13 大连海事大学 Marine transfer trestle wave compensation control system and working method thereof
CN113431814A (en) * 2021-06-17 2021-09-24 江苏科技大学 Synchronous control method based on parallel motion of multiple hydraulic cylinders of heave compensation platform

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102008024513B4 (en) * 2008-05-21 2017-08-24 Liebherr-Werk Nenzing Gmbh Crane control with active coast sequence
AU2018227805B2 (en) * 2017-02-28 2021-05-27 J. Ray Mcdermott, S.A. Offshore ship-to-ship lifting with target tracking assistance

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014131884A (en) * 2013-01-04 2014-07-17 Shimonoseki Ryoju Engineering Kk Embarking-disembarking device
CN105408199A (en) * 2013-07-19 2016-03-16 Icd软件股份公司 Apparatus and method for providing active motion compensation control of an articulated gangway
WO2016193713A1 (en) * 2015-06-02 2016-12-08 Marine Electrical Consulting Limited Method and apparatus for adaptive motion compensation
CN107675607A (en) * 2016-10-28 2018-02-09 福建省新能海上风电研发中心有限公司 A kind of six degree of freedom Active Compensation formula offshore platform steps on the application method for multiplying trestle
CN107434010A (en) * 2017-09-26 2017-12-05 哈尔滨工程大学 A kind of electronic wave Active Compensation, which is stepped on, multiplies system and its control method
CN111045332A (en) * 2019-12-27 2020-04-21 哈尔滨工程大学 Unmanned ship path tracking guidance strategy and disturbance compensation method
CN113031429A (en) * 2021-02-26 2021-06-25 北京星光凯明智能科技有限公司 Shipborne carrier stabilizing platform and control method
CN113104153A (en) * 2021-04-25 2021-07-13 大连海事大学 Marine transfer trestle wave compensation control system and working method thereof
CN113431814A (en) * 2021-06-17 2021-09-24 江苏科技大学 Synchronous control method based on parallel motion of multiple hydraulic cylinders of heave compensation platform

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
一种用于登乘栈桥的主动波浪补偿方法;苏长青等;《船舶与海洋工程》;20170825;第33卷(第04期);第22-25页 *

Also Published As

Publication number Publication date
CN114735140A (en) 2022-07-12

Similar Documents

Publication Publication Date Title
CN114735140B (en) Method, equipment and medium for compensating disturbance speed of wind power pile boarding trestle
CN107434010B (en) A kind of electronic wave Active Compensation, which is stepped on, multiplies system and its control method
EP2489872B1 (en) Blade load reduction for wind turbine
CN106828789B (en) A kind of high-rise pile cap steps on the wave compensating device by operation
Ren et al. Underactuated control and analysis of single blade installation using a jackup installation vessel and active tugger line force control
Schipf et al. Collective pitch feedforward control of floating wind turbines using lidar
CN104817019B (en) A kind of ship's crane heave compensation method based on the forecast of hanging object heave movement
CN109739248B (en) Ship-borne three-degree-of-freedom parallel stable platform stability control method based on ADRC
CN116627043B (en) Regional power positioning control method of combined anchoring system
CN115826606B (en) Dynamic positioning control method for self-elevating ship platform
CN105253264A (en) Ocean wave compensation device of deepwater semisubmersible drilling platform and control method thereof
Zheng et al. Sea trial test on offshore integration of an oscillating buoy wave energy device and floating breakwater
CN117163219B (en) Shipborne trestle feedforward fuzzy control method considering constraint between long rods
CN115495935B (en) Modeling method and device of floating wind turbine generator
CN115107935A (en) Rapid and accurate wave compensation control method for offshore corridor bridge
CN116126003A (en) Wave compensation system modeling and pose control method based on Stewart platform
CN115560796A (en) Floating type fan real-time monitoring and intelligent control system based on digital twin and environmental test
CN112849359B (en) Method for determining main scale of ultra-large piling ship based on wave resistance calculation
CN114879504A (en) Self-adaptive nonlinear control method of four-degree-of-freedom marine rotary crane
Zhang et al. Wave heave compensation based on an optimized backstepping control method
CN111170176B (en) Active stabilization control method suitable for offshore or onshore load hoisting
Sreenivasan et al. Numerical simulations of offshore wind turbine component lifts to determine installation operability
Yu et al. Model Test and Validation of the Crown Floating Offshore Wind Turbine
Xu et al. Interactive visual reality of the offshore hoisting operation and numerical modeling
CN109726519A (en) A kind of movement load calculation method of cutter suction dredger under longitudinal buffer effect

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