[go: up one dir, main page]

CN114735140A - 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
CN114735140A
CN114735140A CN202210380956.3A CN202210380956A CN114735140A CN 114735140 A CN114735140 A CN 114735140A CN 202210380956 A CN202210380956 A CN 202210380956A CN 114735140 A CN114735140 A CN 114735140A
Authority
CN
China
Prior art keywords
increment
value
formula
interference
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.)
Granted
Application number
CN202210380956.3A
Other languages
Chinese (zh)
Other versions
CN114735140B (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

本发明提出一种风电桩登乘栈桥的干扰速度补偿方法、设备和介质。本发明通过引入干扰速度(船舶运动速度)补偿的方式与干扰位置补偿相结合,提高主动式运动补偿栈桥的快速性和补偿效果,并改善栈桥搭接过程的使用体验,提高恶劣天气下风电桩的可达性,增加登乘人员的安全性。

Figure 202210380956

The present invention provides a method, equipment and medium for compensating for interference speed of a wind power pile boarding trestle. The invention combines the interference speed (vessel motion speed) compensation method with the interference position compensation, improves the rapidity and compensation effect of the active motion compensation trestle, improves the use experience of the trestle overlapping process, and improves the wind power pile in bad weather. improve the accessibility and increase the safety of boarding personnel.

Figure 202210380956

Description

一种风电桩登乘栈桥的干扰速度补偿方法、设备和介质Interference speed compensation method, equipment and medium for wind power pile boarding trestle

技术领域technical field

本发明属于具有波浪补偿功能的海上风电登乘栈桥技术领域,特别是涉及一种风电桩登乘栈桥的干扰速度补偿方法、设备和介质。具体说是一种对船舶运动速度进行补偿,以提高运动补偿精度,实现高海况条件下栈桥顺利搭接海上风电桩的控制方法。The invention belongs to the technical field of offshore wind power boarding trestle with wave compensation function, and particularly relates to a disturbance speed compensation method, equipment and medium for a wind power pile boarding trestle. Specifically, it is a control method for compensating the motion speed of the ship to improve the motion compensation accuracy and realize the smooth lap of the trestle to the offshore wind power pile under high sea conditions.

背景技术Background technique

风能是清洁的可再生能源。相比陆地风力发电,海上风电具有不占用陆地面积和对鸟类影响小等优点,此外海上风电还具有静风期少,风力大等优点。我国在近海和中远海均具有丰富的风能,发展海上风电对优化我国能源结构具有重要的意义。Wind energy is clean and renewable energy. Compared with onshore wind power, offshore wind power has the advantages of not occupying land area and less impact on birds. In addition, offshore wind power has the advantages of less quiet wind period and large wind power. my country has abundant wind energy in both offshore and COSCO seas. The development of offshore wind power is of great significance to optimizing my country's energy structure.

然而,海上风电的安装和维护成本要比陆地风电高。在恶劣的天气条件下,运维船剧烈的摇晃使得运维人员登乘风电桩变得困难和危险。主动式运动补偿(Active MotionCompensation,AMC)栈桥通过补偿某几个自由度的船舶运动,使得栈桥末端相对大地坐标系(固定于大地的风电桩)保持静止或者减少海浪对船舶运动的影响,用以提高恶略天气下人员对风电桩的可达性。However, installation and maintenance costs for offshore wind are higher than for onshore wind. In severe weather conditions, the violent shaking of the operation and maintenance vessel makes it difficult and dangerous for the operation and maintenance personnel to board the wind power piles. Active Motion Compensation (AMC) trestle bridge can keep the end of the trestle stationary relative to the geodetic coordinate system (wind power pile fixed to the ground) or reduce the influence of waves on the motion of the ship by compensating the ship motion of certain degrees of freedom. Improve accessibility to wind farms in bad weather.

三自由度AMC栈桥如图1所示,主要由(B1)底座、(B2)支座、(B3)俯仰桥、(B4) 伸缩桥和位于伸缩桥末端的(B5)插入板组成。当插入板插入风电机组平台的基础爬梯后维护人员通过支座、俯仰桥和伸缩桥到达爬梯,登上海上风电机组平台。根据安全规范要求,当运维船与海上风电机组爬梯的距离不超过成人一半的步长时,才能登乘风电平台的爬梯。然而海上潮汐和海浪的影响,会导致船舶和海上风电平台的爬梯在垂直方向和水平方向上存在一定的距离,该距离会随着船舶的摇晃而变化。The three-degree-of-freedom AMC trestle is shown in Figure 1, which is mainly composed of (B1) base, (B2) support, (B3) pitch bridge, (B4) telescopic bridge and (B5) insert plate at the end of the telescopic bridge. When the insert board is inserted into the foundation ladder of the wind turbine platform, the maintenance personnel reach the ladder through the support, pitch bridge and telescopic bridge, and board the offshore wind turbine platform. According to the requirements of safety regulations, only when the distance between the operation and maintenance vessel and the ladder of the offshore wind turbine is not more than half the step length of an adult can you board the ladder of the wind power platform. However, the influence of sea tides and waves will cause a certain distance between the ship and the ladder of the offshore wind power platform in the vertical and horizontal directions, and the distance will change with the shaking of the ship.

安装在船舶上的运动测量单元,会测量船舶六个自由度的运动,通过(Y1)回转液压缸、(Y2)俯仰液压缸和(Y3)伸缩液压缸补偿船舶的首摇、纵摇和垂荡运动导致(B5) 插入板相对爬梯的位移变化。同时整个栈桥也可以保持相对稳定的状态。The motion measurement unit installed on the ship will measure the motion of the ship's six degrees of freedom, and compensate the ship's yaw, pitch and sag through (Y1) slewing hydraulic cylinder, (Y2) pitching hydraulic cylinder and (Y3) telescopic hydraulic cylinder. The oscillating motion results in a change in the displacement of the (B5) insert plate relative to the ladder. At the same time, the entire trestle can also maintain a relatively stable state.

当海上风浪很小时,操作人员可以关闭主动式运动补偿功能,通过操作手柄,控制比例阀组实现对回转、俯仰和伸缩液压缸活塞杆运动速度的控制,能够顺利地使伸缩桥末端的插入板插入爬梯,完成人员通道的建立。而大风浪状态下,位于伸缩桥末端的插入板会放大船舶的摇晃,使得插入爬梯的过程变的困难。此时,需要开启主动补偿功能,消除或减小船舶摇晃对插入板影响,当操作手柄无动作时,插入板末端的坐标保持相对大地坐标系的静止或减小船舶运动的影响。为了完成主动运动补偿,需要利用运动参考单元(MotionReference Unit,MRU)测量出船舶在六个自由度的运动,进而通过转换矩阵计算出栈桥在回转自由度、俯仰自由度和伸缩自由度的补偿量,来补偿船舶在首摇、纵摇和垂荡方向的运动。为了保持插入板相对大地坐标系的静止,主动式运动补偿栈桥需要补偿船舶运动导致的插入板位置的变化,因此栈桥在回转、俯仰和伸缩三个自由度上的控制都是位置伺服系统,根据位置误差进行液压缸活塞杆伸出量的控制。然而这种补偿方式,在主动运动补偿的过程中主要存在以下问题:When the offshore wind and waves are very small, the operator can turn off the active motion compensation function, and control the proportional valve group to control the movement speed of the slewing, pitching and telescopic hydraulic cylinder piston rod through the operating handle, which can smoothly make the insertion plate at the end of the telescopic bridge. Insert the ladder to complete the establishment of the personnel passage. In the case of strong wind and waves, the insertion plate at the end of the telescopic bridge will magnify the shaking of the ship, making the process of inserting the ladder difficult. At this time, it is necessary to turn on the active compensation function to eliminate or reduce the influence of ship shaking on the insertion plate. When the operating handle does not move, the coordinates of the end of the insertion plate remain relatively static in the geodetic coordinate system or reduce the influence of ship motion. In order to complete the active motion compensation, it is necessary to use the Motion Reference Unit (MRU) to measure the motion of the ship in six degrees of freedom, and then calculate the compensation amount of the trestle in the rotational degrees of freedom, pitching degrees of freedom and telescopic degrees of freedom through the transformation matrix. , to compensate for the motion of the ship in the yaw, pitch and heave directions. In order to keep the insertion plate still relative to the geodetic coordinate system, the active motion compensation trestle needs to compensate for the change of the position of the insertion plate caused by the movement of the ship. Therefore, the control of the trestle in the three degrees of freedom of rotation, pitch and telescopic is a position servo system. According to The position error controls the extension of the hydraulic cylinder piston rod. However, this compensation method mainly has the following problems in the process of active motion compensation:

(1)由于船舶运动随机性,MRU的采样延迟,液压缸运动速度的限制,栈桥安装的机械间隙等影响,在主动运动补偿栈桥工作过程中,对于相对高频的船舶运动,仅仅依靠测量的船舶位置进行补偿,会导致跟踪滞后,补偿效果欠佳。(1) Due to the randomness of ship motion, the sampling delay of MRU, the limitation of hydraulic cylinder movement speed, the mechanical clearance of trestle installation, etc., during the active motion compensation trestle work process, for relatively high-frequency ship motion, only the measured Compensation for ship position will lead to tracking lag and poor compensation effect.

(2)高海况下操作手柄驱动插入板插入风电机组平台的基础爬梯过程中,仅仅将操作手柄的速度控制与主动运动补偿的位置控制进行简单叠加,不能获得与平静海面下关闭补偿功能时操作手柄过程一致的体验,操作不平顺。(2) In the process of inserting the operation handle to drive the insert plate into the basic ladder of the wind turbine platform under high sea conditions, only the speed control of the operation handle and the position control of active motion compensation are simply superimposed, and the operation when the compensation function is turned off under calm sea cannot be obtained. The experience of the handle process is consistent, and the operation is not smooth.

发明内容SUMMARY OF THE INVENTION

本发明目的是为了解决现有技术中的问题,提出了一种风电桩登乘栈桥的干扰速度补偿方法、设备和介质。本发明通过引入干扰速度(船舶运动速度)补偿的方式与干扰位置补偿相结合,提高主动式运动补偿栈桥的快速性和补偿效果,并改善栈桥搭接过程的使用体验,提高恶劣天气下风电桩的可达性,增加登乘人员的安全性。The purpose of the present invention is to solve the problems in the prior art, and proposes a method, equipment and medium for compensating the disturbance speed of a wind power pile boarding trestle. The invention combines the interference speed (vessel motion speed) compensation method with the interference position compensation, improves the rapidity and compensation effect of the active motion compensation trestle, improves the use experience of the trestle overlapping process, and improves the wind power pile in bad weather. improve the accessibility and increase the safety of boarding personnel.

本发明是通过以下技术方案实现的,本发明提出一种风电桩登乘栈桥的干扰速度补偿方法,所述方法具体包括以下步骤:The present invention is achieved through the following technical solutions, and the present invention proposes a method for compensating for the disturbance speed of a wind power pile boarding trestle, the method specifically includes the following steps:

步骤1、设插入板的位置为x=[x1,x2,x3]T,x1,x2,x3分别为插入板在大地坐标系O-XYZ三个坐标轴的坐标值,则该坐标值表示为:Step 1. Set the position of the insertion board as x=[x 1 , x 2 , x 3 ] T , where x 1 , x 2 , and x 3 are the coordinate values of the insertion board in the three coordinate axes of the geodetic coordinate system O-XYZ, respectively, Then the coordinate value is expressed as:

Figure BDA0003592964980000021
Figure BDA0003592964980000021

式中,A、B、C和D为系统矩阵,u(k)=[r(k),p(k),s(k)]T分别为栈桥的回转值、俯仰值和伸缩值;d(k)=[rs(k),ps(k),hs(k)]T分别为船舶的首摇值、纵摇值和垂荡值,输出值 y=[x1,x2,x3]T,k表示时间序列;In the formula, A, B, C and D are system matrices, u(k)=[r(k), p(k), s(k)] T is the slewing value, pitch value and telescopic value of the trestle, respectively; d (k)=[r s (k), p s (k), h s (k)] T are the yaw value, pitch value and heave value of the ship respectively, and the output value y=[x 1 , x 2 ,x 3 ] T , k represents time series;

步骤2、将方程(1)变为差分形式,设置包换位置和速度跟踪误差的目标函数J(k),通过目标函数确定控制量的增量;Step 2. Change equation (1) into a differential form, set the objective function J(k) of the replacement position and velocity tracking error, and determine the increment of the control amount through the objective function;

步骤3、通过扩张状态观测器ESO获得船舶运动干扰的增量;Step 3. Obtain the increment of ship motion disturbance through the expanded state observer ESO;

步骤4、综合控制量增量和干扰量增量获得最终的控制量;Step 4. The final control amount is obtained by synthesizing the control amount increment and the disturbance amount increment;

步骤5、将开启AMC功能时刻的系统状态记为x(0),当操作手柄无动作时,插入板的目标位置xr(k)=x(0),插入板的目标速度

Figure BDA0003592964980000031
此时的控制量用来补偿船舶运动导致的插入板坐标的变化;当操作手柄动作时,插入板的目标速度
Figure BDA0003592964980000032
与操作手柄的输出电压成比例,xr(k)为
Figure BDA0003592964980000033
的积分,以跟踪手柄的操作。Step 5. Denote the system state when the AMC function is turned on as x(0). When the operating handle does not move, the target position of the insertion board is xr (k)=x(0), and the target speed of the insertion board is
Figure BDA0003592964980000031
The control amount at this time is used to compensate for the change of the coordinates of the insertion plate caused by the movement of the ship; when the operation handle moves, the target speed of the insertion plate
Figure BDA0003592964980000032
Proportional to the output voltage of the joystick, x r (k) is
Figure BDA0003592964980000033
Integral to track the operation of the handle.

进一步地,将方程(1)变成差分形式:Further, transform equation (1) into differential form:

Figure BDA0003592964980000034
Figure BDA0003592964980000034

式中,

Figure BDA0003592964980000035
分别为状态量、控制量和干扰量的增量;系统矩阵中
Figure BDA0003592964980000036
Figure BDA0003592964980000037
O为零矩阵,I为单位矩阵。In the formula,
Figure BDA0003592964980000035
are the increments of state quantity, control quantity and disturbance quantity respectively; in the system matrix
Figure BDA0003592964980000036
Figure BDA0003592964980000037
O is a zero matrix and I is an identity matrix.

进一步地,考虑控制变量u(k)的幅值约束和速度约束,通过解下述目标函数来获得控制变量的增量:Further, considering the amplitude constraint and velocity constraint of the control variable u(k), the increment of the control variable is obtained by solving the following objective function:

Figure BDA0003592964980000038
Figure BDA0003592964980000038

Figure BDA0003592964980000039
Figure BDA0003592964980000039

式中,Np为预测时域,矩阵Q和S为半正定矩阵,矩阵R和P为正定矩阵,

Figure BDA00035929649800000310
Figure BDA00035929649800000311
代表状态误差,xr(k)为插入板的目标位置,xd(k)为x(k)的微分即插入板的速度,因此误差中既包含了位置误差也包含了速度误差;
Figure BDA00035929649800000312
为控制量增量,增加的自由度
Figure BDA00035929649800000313
用来模拟干扰量的增量。In the formula, N p is the prediction time domain, matrices Q and S are semi-positive definite matrices, matrices R and P are positive definite matrices,
Figure BDA00035929649800000310
Figure BDA00035929649800000311
represents the state error, x r (k) is the target position of the inserted board, x d (k) is the differential of x (k), that is, the speed of the inserted board, so the error includes both the position error and the speed error;
Figure BDA00035929649800000312
To control the increment, the added degrees of freedom
Figure BDA00035929649800000313
Increment used to simulate the amount of disturbance.

进一步地,最小化目标函数J(k)获得如下的控制序列:Further, minimize the objective function J(k) to obtain the following control sequence:

Figure BDA00035929649800000314
Figure BDA00035929649800000314

因此,在每一个采样周期解目标函数问题可以转化为如下的二次规划:Therefore, solving the objective function problem at each sampling period can be transformed into the following quadratic programming:

Figure BDA00035929649800000315
Figure BDA00035929649800000315

Subject to Lη≤bSubject to Lη≤b

式中,L和b用来限制控制量的幅值和速度;矩阵F和H如下定义,并通过实时计算获得:In the formula, L and b are used to limit the amplitude and speed of the control variable; the matrices F and H are defined as follows and obtained by real-time calculation:

Figure BDA0003592964980000041
Figure BDA0003592964980000041

Figure BDA0003592964980000042
Figure BDA0003592964980000042

Figure BDA0003592964980000043
Figure BDA0003592964980000043

Figure BDA0003592964980000044
Figure BDA0003592964980000044

状态

Figure BDA0003592964980000045
可表示为:state
Figure BDA0003592964980000045
can be expressed as:

Figure BDA0003592964980000046
Figure BDA0003592964980000046

进一步地,采用式(5)实时优化η*,可以获得控制量增量的优化值

Figure BDA0003592964980000047
同时,也可以获得优化的干扰值增量
Figure BDA0003592964980000048
因此控制量增量可以通过下式获得:Further, using formula (5) to optimize η * in real time, the optimized value of the control amount increment can be obtained
Figure BDA0003592964980000047
At the same time, an optimized disturbance value increment can also be obtained
Figure BDA0003592964980000048
Therefore, the control amount increment can be obtained by the following formula:

Figure BDA0003592964980000049
Figure BDA0003592964980000049

式中,μ(k)用来补偿干扰

Figure BDA00035929649800000410
其通过
Figure BDA00035929649800000411
获得;最小化误差E(k),μ(k)满足where μ(k) is used to compensate for interference
Figure BDA00035929649800000410
its passed
Figure BDA00035929649800000411
Obtain; minimize the error E(k), μ(k) satisfies

Figure BDA00035929649800000412
Figure BDA00035929649800000412

式中,

Figure BDA00035929649800000413
是干扰的估计值,同时μ(k)通过下式计算:In the formula,
Figure BDA00035929649800000413
is an estimate of the interference, while μ(k) is calculated by:

Figure BDA00035929649800000414
Figure BDA00035929649800000414

进一步地,为了获得干扰量增量

Figure BDA00035929649800000415
的估计值
Figure BDA00035929649800000416
采用如下的扩张状态观测器ESO 来观测干扰:Further, in order to obtain the disturbance amount increment
Figure BDA00035929649800000415
estimated value of
Figure BDA00035929649800000416
The disturbance is observed using the extended state observer ESO as follows:

Figure BDA0003592964980000051
Figure BDA0003592964980000051

式中,

Figure BDA0003592964980000052
代表*的估计值,
Figure BDA0003592964980000053
代表x微分的估计值,
Figure BDA0003592964980000054
代表x二次微分的估计值;Ts为采样周期,θ是观测器增益,ko1,ko2,ko3和ko4为观测器系数;函数gi(i=1,2,3,4) 表示为:In the formula,
Figure BDA0003592964980000052
represents the estimated value of *,
Figure BDA0003592964980000053
represents the estimated value of the derivative of x,
Figure BDA0003592964980000054
represents the estimated value of the second derivative of x; T s is the sampling period, θ is the observer gain, k o1 , k o2 , k o3 and k o4 are the observer coefficients; the function gi (i=1,2,3,4 ) Expressed as:

Figure BDA0003592964980000055
Figure BDA0003592964980000055

式中,α1=γ,α2=2γ-1,α3=3γ-2,α4=4γ-3,γ∈(3/4,1);K为正定矩阵;基于式(10)对d(k)的估计,获得

Figure BDA0003592964980000056
并通过式(9)计算μ(k);进而优化的控制变量为In the formula, α 1 =γ,α 2 =2γ-1,α 3 =3γ-2,α 4 =4γ-3,γ∈(3/4,1); K is a positive definite matrix; Estimate of d(k), get
Figure BDA0003592964980000056
And calculate μ(k) by formula (9); then the optimized control variable is

Figure BDA0003592964980000057
Figure BDA0003592964980000057

进一步地,通过液压缸行程和机械结构,确定控制量的最大最小值分别为usU、usL;通过比例阀流量确定控制量增量的最大最小值分别为

Figure BDA0003592964980000058
即Further, according to the hydraulic cylinder stroke and mechanical structure, the maximum and minimum values of the control variables are determined to be u sU and u sL respectively; the maximum and minimum values of the control variable increments determined by the proportional valve flow are respectively:
Figure BDA0003592964980000058
which is

Figure BDA0003592964980000059
Figure BDA0003592964980000059

式(9)中的μ(k)可被重新描述为:μ(k) in equation (9) can be re-described as:

Figure BDA00035929649800000510
Figure BDA00035929649800000510

式中,

Figure BDA00035929649800000511
因此控制量的幅值约束可采用如下不等式表示:In the formula,
Figure BDA00035929649800000511
Therefore, the amplitude constraint of the control variable can be expressed by the following inequality:

Figure BDA00035929649800000512
Figure BDA00035929649800000512

式中,In the formula,

Figure BDA0003592964980000061
Figure BDA0003592964980000061

同样控制量增量的约束可采用如下不等式表示:Similarly, the constraint on the increment of the control quantity can be expressed by the following inequality:

Figure BDA0003592964980000062
Figure BDA0003592964980000062

式中,

Figure BDA0003592964980000063
Figure BDA0003592964980000064
分别为
Figure BDA0003592964980000065
Figure BDA0003592964980000066
的增量;In the formula,
Figure BDA0003592964980000063
and
Figure BDA0003592964980000064
respectively
Figure BDA0003592964980000065
and
Figure BDA0003592964980000066
increment;

Figure BDA0003592964980000067
Figure BDA0003592964980000067

本发明提出一种电子设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现所述一种风电桩登乘栈桥的干扰速度补偿方法的步骤。The present invention provides an electronic device comprising a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the processor implements the steps of the method for compensating for the disturbance speed of a wind power pile boarding trestle.

本发明提出一种计算机可读存储介质,用于存储计算机指令,所述计算机指令被处理器执行时实现所述一种风电桩登乘栈桥的干扰速度补偿方法的步骤。The present invention provides a computer-readable storage medium for storing computer instructions, and when the computer instructions are executed by a processor, the steps of the method for compensating for the disturbance speed of a wind power pile boarding trestle are realized.

本发明的有益效果为:The beneficial effects of the present invention are:

(1)采用栈桥模型的差分形式,并在误差中加入跟踪速度误差,提高了主动式运动补偿的快速性和动态补偿精度;(1) The differential form of the trestle model is adopted, and the tracking speed error is added to the error, which improves the rapidity and dynamic compensation accuracy of active motion compensation;

(2)将与操作手柄输出电压成比例的值作为跟踪速度,将跟踪速度的积分作目标位置,改善手柄的操纵体验;(2) The value proportional to the output voltage of the operating handle is used as the tracking speed, and the integral of the tracking speed is used as the target position to improve the operating experience of the handle;

(3)控制量求解控制中考虑液压缸运动的速度限制和栈桥的位置约束,使获得的控制量是物理上的可行控制,避免了依靠限位开关限制运动幅值导致的紧急制动,保证了在极限位置的运动平稳性。(3) The speed limit of the hydraulic cylinder movement and the position constraint of the trestle bridge are considered in the control quantity solution control, so that the obtained control quantity is a physically feasible control, avoiding the emergency braking caused by the limit switch to limit the movement amplitude, ensuring that the The smoothness of motion at the extreme position.

附图说明Description of drawings

图1为AMC栈桥组成示意图;其中Y1回转液压缸;Y2俯仰液压缸;Y3伸缩液压缸;B1底座;B2支座;B3俯仰桥;B4伸缩桥;B5插入板;Figure 1 is a schematic diagram of the composition of the AMC trestle; Y1 rotary hydraulic cylinder; Y2 pitching hydraulic cylinder; Y3 telescopic hydraulic cylinder; B1 base; B2 support; B3 pitching bridge; B4 telescopic bridge; B5 insert plate;

图2为栈桥搭接风电桩示意图;其中B4伸缩桥;B5插入板;C1风电机组平台;C2 爬梯;Figure 2 is a schematic diagram of the trestle lapped wind power pile; wherein B4 telescopic bridge; B5 insert plate; C1 wind turbine platform; C2 climbing ladder;

图3为干扰速度补偿功能原理框图。Figure 3 is a functional block diagram of the interference speed compensation.

具体实施方式Detailed ways

下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

结合图1-图3,本发明提出一种风电桩登乘栈桥的干扰速度补偿方法,栈桥的搭接过程如图2所示,插入板插入风电机组平台的基础爬梯,完成人员通道的建立;所述方法具体包括以下步骤:In conjunction with Fig. 1-Fig. 3, the present invention proposes a method for compensating for the interference speed of a wind power pile boarding trestle. The lap joint process of the trestle is shown in Fig. 2, and the insert plate is inserted into the foundation ladder of the wind turbine platform to complete the establishment of the personnel passage; The method specifically includes the following steps:

步骤1、设插入板的位置为x=[x1,x2,x3]T,x1,x2,x3分别为插入板在大地坐标系O-XYZ三个坐标轴的坐标值,则该坐标值表示为:Step 1. Set the position of the insertion board as x=[x 1 , x 2 , x 3 ] T , where x 1 , x 2 , and x 3 are the coordinate values of the insertion board in the three coordinate axes of the geodetic coordinate system O-XYZ, respectively, Then the coordinate value is expressed as:

Figure BDA0003592964980000071
Figure BDA0003592964980000071

式中,A、B、C和D为系统矩阵,u(k)=[r(k),p(k),s(k)]T分别为栈桥的回转值、俯仰值和伸缩值;d(k)=[rs(k),ps(k),hs(k)]T分别为船舶的首摇值、纵摇值和垂荡值,输出值 y=[x1,x2,x3]T,k表示时间序列;In the formula, A, B, C and D are system matrices, u(k)=[r(k), p(k), s(k)] T is the slewing value, pitch value and telescopic value of the trestle, respectively; d (k)=[r s (k), p s (k), h s (k)] T are the yaw value, pitch value and heave value of the ship respectively, and the output value y=[x 1 , x 2 ,x 3 ] T , k represents time series;

步骤2、为了处理快速变化的船舶运动,将方程(1)变为差分形式,设置包换位置和速度跟踪误差的目标函数J(k),通过目标函数确定控制量的增量;Step 2. In order to deal with the rapidly changing ship motion, the equation (1) is changed into a differential form, the objective function J(k) of the replacement position and speed tracking error is set, and the increment of the control amount is determined by the objective function;

步骤3、通过扩张状态观测器ESO获得船舶运动干扰的增量;Step 3. Obtain the increment of ship motion disturbance through the expanded state observer ESO;

步骤4、综合控制量增量和干扰量增量获得最终的控制量;Step 4. The final control amount is obtained by synthesizing the control amount increment and the disturbance amount increment;

步骤5、将开启AMC功能时刻的系统状态记为x(0),当操作手柄无动作时,插入板的目标位置xr(k)=x(0),插入板的目标速度

Figure BDA0003592964980000072
此时的控制量用来补偿船舶运动导致的插入板坐标的变化,误差中包含了速度误差;当操作手柄动作时,插入板的目标速度
Figure BDA0003592964980000073
与操作手柄的输出电压成比例,xr(k)为
Figure BDA0003592964980000074
的积分,以跟踪手柄的操作。当操作手柄时,式(3)中插入板的目标位置
Figure BDA0003592964980000075
(数字控制系统中以累加和代替积分值),插入板末端的目标速度
Figure BDA0003592964980000076
(vr为操作手柄三个自由度的输出电压,kv为比例系数)。此时的控制量除了用来补偿船舶运动导致的插入板坐标的变化,还同时响应手柄的操作。通过调整kv获得与通过vr直接控制比例阀开口一致的速度,使手柄操作更加平顺,改善开启AMC功能后手柄的操作体验。Step 5. Denote the system state when the AMC function is turned on as x(0). When the operating handle does not move, the target position of the insertion board is xr (k)=x(0), and the target speed of the insertion board is
Figure BDA0003592964980000072
The control value at this time is used to compensate for the change of the coordinates of the insertion plate caused by the movement of the ship, and the error includes the speed error; when the operation handle moves, the target speed of the insertion plate
Figure BDA0003592964980000073
Proportional to the output voltage of the joystick, x r (k) is
Figure BDA0003592964980000074
Integral to track the operation of the handle. When the handle is operated, the target position of the insertion plate in formula (3)
Figure BDA0003592964980000075
(In the digital control system, the accumulated sum replaces the integral value), insert the target speed of the end of the board
Figure BDA0003592964980000076
(v r is the output voltage of the three degrees of freedom of the operating handle, and k v is the proportional coefficient). The control amount at this time is not only used to compensate the change of the coordinates of the insertion plate caused by the movement of the ship, but also responds to the operation of the handle. By adjusting k v to obtain the same speed as directly controlling the opening of the proportional valve through v r , the operation of the handle is smoother and the operation experience of the handle after the AMC function is turned on is improved.

将方程(1)变成差分形式:Transform equation (1) into differential form:

Figure BDA0003592964980000081
Figure BDA0003592964980000081

式中,

Figure BDA0003592964980000082
分别为状态量、控制量和干扰量的增量;系统矩阵中
Figure BDA0003592964980000083
Figure BDA0003592964980000084
O为零矩阵,I为单位矩阵。In the formula,
Figure BDA0003592964980000082
are the increments of state quantity, control quantity and disturbance quantity respectively; in the system matrix
Figure BDA0003592964980000083
Figure BDA0003592964980000084
O is a zero matrix and I is an identity matrix.

考虑控制变量u(k)的幅值约束和速度约束,通过解下述目标函数来获得控制变量的增量:Considering the amplitude constraint and velocity constraint of the control variable u(k), the increment of the control variable is obtained by solving the following objective function:

Figure BDA0003592964980000085
Figure BDA0003592964980000085

Figure BDA0003592964980000086
Figure BDA0003592964980000086

式中,Np为预测时域,矩阵Q和S为半正定矩阵,矩阵R和P为正定矩阵,

Figure BDA0003592964980000087
Figure BDA0003592964980000088
代表状态误差,xr(k)为插入板的目标位置,xd(k)为x(k)的微分即插入板的速度,因此误差中既包含了位置误差也包含了速度误差;
Figure BDA0003592964980000089
为控制量增量,增加的自由度
Figure BDA00035929649800000810
用来模拟干扰量的增量。In the formula, N p is the prediction time domain, matrices Q and S are semi-positive definite matrices, matrices R and P are positive definite matrices,
Figure BDA0003592964980000087
Figure BDA0003592964980000088
represents the state error, x r (k) is the target position of the inserted board, x d (k) is the differential of x (k), that is, the speed of the inserted board, so the error includes both the position error and the speed error;
Figure BDA0003592964980000089
To control the increment, the added degrees of freedom
Figure BDA00035929649800000810
Increment used to simulate the amount of disturbance.

目标函数中包含了速度误差和干扰量的增量(等效于干扰速度),因此能够反映快速变化的船舶运动。The velocity error and the increment of disturbance amount (equivalent to disturbance velocity) are included in the objective function, so it can reflect the rapidly changing ship motion.

最小化目标函数J(k)获得如下的控制序列:Minimizing the objective function J(k) obtains the following control sequence:

Figure BDA00035929649800000811
Figure BDA00035929649800000811

因此,在每一个采样周期解目标函数问题可以转化为如下的二次规划:Therefore, solving the objective function problem at each sampling period can be transformed into the following quadratic programming:

Figure BDA00035929649800000812
Figure BDA00035929649800000812

Subject to Lη≤bSubject to Lη≤b

式中,L和b用来限制控制量的幅值和速度;矩阵F和H如下定义,并通过实时计算获得:In the formula, L and b are used to limit the amplitude and speed of the control variable; the matrices F and H are defined as follows and obtained by real-time calculation:

Figure BDA0003592964980000091
Figure BDA0003592964980000091

Figure BDA0003592964980000092
Figure BDA0003592964980000092

Figure BDA0003592964980000093
Figure BDA0003592964980000093

Figure BDA0003592964980000094
Figure BDA0003592964980000094

状态

Figure BDA0003592964980000095
可表示为:state
Figure BDA0003592964980000095
can be expressed as:

Figure BDA0003592964980000096
Figure BDA0003592964980000096

采用式(5)实时优化η*,可以获得控制量增量的优化值

Figure BDA0003592964980000097
同时,也可以获得优化的干扰值增量
Figure BDA0003592964980000098
因此控制量增量可以通过下式获得:Using formula (5) to optimize η * in real time, the optimized value of the control amount increment can be obtained
Figure BDA0003592964980000097
At the same time, an optimized disturbance value increment can also be obtained
Figure BDA0003592964980000098
Therefore, the control amount increment can be obtained by the following formula:

Figure BDA0003592964980000099
Figure BDA0003592964980000099

式中,μ(k)用来补偿干扰

Figure BDA00035929649800000910
其通过
Figure BDA00035929649800000911
获得;最小化误差E(k),μ(k)满足where μ(k) is used to compensate for interference
Figure BDA00035929649800000910
its passed
Figure BDA00035929649800000911
Obtain; minimize the error E(k), μ(k) satisfies

Figure BDA00035929649800000912
Figure BDA00035929649800000912

式中,

Figure BDA00035929649800000913
是干扰的估计值,同时μ(k)通过下式计算:In the formula,
Figure BDA00035929649800000913
is an estimate of the interference, while μ(k) is calculated by:

Figure BDA00035929649800000914
Figure BDA00035929649800000914

为了获得干扰量增量

Figure BDA00035929649800000915
的估计值
Figure BDA00035929649800000916
采用如下的扩张状态观测器ESO来观测干扰:In order to obtain the amount of disturbance increment
Figure BDA00035929649800000915
estimated value of
Figure BDA00035929649800000916
The disturbance is observed using the extended state observer ESO as follows:

Figure BDA0003592964980000101
Figure BDA0003592964980000101

式中,

Figure BDA0003592964980000102
代表*的估计值,
Figure BDA0003592964980000103
代表x微分的估计值,
Figure BDA0003592964980000104
代表x二次微分的估计值;Ts为采样周期,θ是观测器增益,ko1,ko2,ko3和ko4为观测器系数;函数gi(i=1,2,3,4) 表示为:In the formula,
Figure BDA0003592964980000102
represents the estimated value of *,
Figure BDA0003592964980000103
represents the estimated value of the derivative of x,
Figure BDA0003592964980000104
represents the estimated value of the second derivative of x; T s is the sampling period, θ is the observer gain, k o1 , k o2 , k o3 and k o4 are the observer coefficients; the function gi (i=1,2,3,4 ) Expressed as:

Figure BDA0003592964980000105
Figure BDA0003592964980000105

式中,α1=γ,α2=2γ-1,α3=3γ-2,α4=4γ-3,γ∈(3/4,1);K为正定矩阵;基于式(10)对d(k)的估计,获得

Figure BDA0003592964980000106
并通过式(9)计算μ(k);进而优化的控制变量为In the formula, α 1 =γ,α 2 =2γ-1,α 3 =3γ-2,α 4 =4γ-3,γ∈(3/4,1); K is a positive definite matrix; Estimate of d(k), get
Figure BDA0003592964980000106
And calculate μ(k) by formula (9); then the optimized control variable is

Figure BDA0003592964980000107
Figure BDA0003592964980000107

通过液压缸行程和机械结构,确定控制量的最大最小值分别为usU、usL;通过比例阀流量确定控制量增量的最大最小值分别为

Figure BDA0003592964980000108
即According to the hydraulic cylinder stroke and mechanical structure, the maximum and minimum values of the control variables are determined as u sU and u sL respectively; the maximum and minimum values of the control variable increments determined by the proportional valve flow are respectively
Figure BDA0003592964980000108
which is

Figure BDA0003592964980000109
Figure BDA0003592964980000109

式(9)中的μ(k)可被重新描述为:μ(k) in equation (9) can be re-described as:

Figure BDA00035929649800001010
Figure BDA00035929649800001010

式中,

Figure BDA00035929649800001011
因此控制量的幅值约束可采用如下不等式表示:In the formula,
Figure BDA00035929649800001011
Therefore, the amplitude constraint of the control variable can be expressed by the following inequality:

Figure BDA00035929649800001012
Figure BDA00035929649800001012

式中,In the formula,

Figure BDA0003592964980000111
Figure BDA0003592964980000111

同样控制量增量的约束可采用如下不等式表示:Similarly, the constraint on the increment of the control quantity can be expressed by the following inequality:

Figure BDA0003592964980000112
Figure BDA0003592964980000112

式中,

Figure BDA0003592964980000113
Figure BDA0003592964980000114
分别为
Figure BDA0003592964980000115
Figure BDA0003592964980000116
的增量;In the formula,
Figure BDA0003592964980000113
and
Figure BDA0003592964980000114
respectively
Figure BDA0003592964980000115
and
Figure BDA0003592964980000116
increment;

Figure BDA0003592964980000117
Figure BDA0003592964980000117

如图3所示,本发明基于干扰速度补偿的主动式运动补偿方法的系统包括操作手柄积分单元,二次规划单元,控制量计算单元,扩张状态观测器单元,栈桥被控对象和MRU。As shown in FIG. 3 , the system of the active motion compensation method based on interference speed compensation of the present invention includes an operation handle integration unit, a quadratic planning unit, a control quantity calculation unit, an expansion state observer unit, a trestle controlled object and an MRU.

本发明提出一种电子设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现所述一种风电桩登乘栈桥的干扰速度补偿方法的步骤。The present invention provides an electronic device comprising a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the processor implements the steps of the method for compensating for the disturbance speed of a wind power pile boarding trestle.

本发明提出一种计算机可读存储介质,用于存储计算机指令,所述计算机指令被处理器执行时实现所述一种风电桩登乘栈桥的干扰速度补偿方法的步骤。The present invention provides a computer-readable storage medium for storing computer instructions, and when the computer instructions are executed by a processor, the steps of the method for compensating for the disturbance speed of a wind power pile boarding trestle are realized.

以上对本发明所提出的一种风电桩登乘栈桥的干扰速度补偿方法、设备和介质进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。The interference speed compensation method, equipment and medium for a wind power pile boarding trestle proposed by the present invention have been described above in detail. In this paper, specific examples are used to illustrate the principles and implementations of the present invention. The description of the above embodiments It is only used to help understand the method of the present invention and its core idea; at the same time, for those of ordinary skill in the art, according to the idea of the present invention, there will be changes in the specific embodiments and application scope. In summary, The contents of this specification should not be construed as limiting the present invention.

Claims (9)

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 an insertion plate as x ═ x1,x2,x3]T,x1,x2,x3The coordinate values of the insert plate in three coordinate axes of the geodetic coordinate system O-XYZ are respectively expressed as follows:
Figure FDA0003592964970000011
wherein A, B, C and D are system matrices, u (k) ═ r (k), p (k), s (k)]TRespectively a rotation value, a pitching value and a stretching value of the trestle; d (k) ═ rs(k),ps(k),hs(k)]TRespectively the ship's yaw value, pitch value and heave value, and the output value y ═ x1,x2,x3]TK represents a time series;
step 2, changing the equation (1) into a differential form, setting a target function J (k) for converting the tracking errors of the position and the speed, and determining the increment of the control quantity through the target function;
step 3, obtaining the increment of ship motion interference 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 operatedr(k) X (0), target speed of insert plate
Figure FDA0003592964970000012
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 FDA0003592964970000013
Proportional to the output voltage of the operating handle, xr(k) Is composed of
Figure FDA0003592964970000014
To track the operation of the handle.
2. The method of claim 1, wherein equation (1) is transformed into a differential form:
Figure FDA0003592964970000015
in the formula (I), the compound is shown in the specification,
Figure FDA0003592964970000016
respectively increment of the state quantity, the control quantity and the interference quantity; in the system matrix
Figure FDA0003592964970000017
Figure FDA0003592964970000018
O is a zero matrix and I is an identity matrix.
3. A method according to claim 2, characterized in that 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 FDA0003592964970000021
Figure FDA0003592964970000022
in the formula, NpFor predicting the time domain, matrices Q and S are semi-positive definite matrices, matrices R and P are positive definite matrices,
Figure FDA0003592964970000023
Figure FDA0003592964970000024
representing a state error, xr(k) For inserting into the target position of the plate, xd(k) The differential of x (k), i.e., the velocity of the insert plate, thus the error includes both position and velocity errors;
Figure FDA0003592964970000025
increased degree of freedom for controlled quantity increments
Figure FDA0003592964970000026
To simulate the increase in the amount of interference.
4. A method according to claim 3, characterized in that minimizing the objective function j (k) results in a control sequence as follows:
Figure FDA0003592964970000027
therefore, solving the objective function problem at each sampling period can be converted to a quadratic programming as follows:
Figure FDA0003592964970000028
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 FDA0003592964970000029
Figure FDA00035929649700000210
Figure FDA00035929649700000211
Figure FDA00035929649700000212
status of state
Figure FDA00035929649700000213
Can be expressed as:
Figure FDA0003592964970000031
5. the method of claim 4, wherein η is optimized in real time using equation (5)*An optimized value of the control quantity increment can be obtained
Figure FDA0003592964970000032
At the same time, optimized interference value increment can be obtained
Figure FDA0003592964970000033
The control amount increment can thus be obtained by the following equation:
Figure FDA0003592964970000034
in which mu (k) is used to compensate for interference
Figure FDA0003592964970000035
Which pass through
Figure FDA0003592964970000036
Obtaining; minimize error E (k), μ (k) satisfies
Figure FDA0003592964970000037
In the formula (I), the compound is shown in the specification,
Figure FDA0003592964970000038
is an estimate of interference, while μ (k) is calculated by:
Figure FDA0003592964970000039
6. method according to claim 5, characterized in that, in order to obtain the interference increment
Figure FDA00035929649700000310
Is estimated value of
Figure FDA00035929649700000311
The disturbance is observed using the extended state observer ESO as follows:
Figure FDA00035929649700000312
in the formula (I), the compound is shown in the specification,
Figure FDA00035929649700000313
an estimate of the value of x is represented,
Figure FDA00035929649700000314
representing the estimated value of the differential of x,
Figure FDA00035929649700000315
represents an estimate of x second derivative; t is a unit ofsFor the sampling period, θ is the observer gain, ko1,ko2,ko3And ko4Is the observer coefficient; function gi(i ═ 1,2,3,4) is represented by:
Figure FDA00035929649700000316
in the formula, alpha1=γ,α2=2γ-1,α3=3γ-2,α44 γ -3, γ ∈ (3/4, 1); k is a positive definite matrix; obtaining based on the estimation of d (k) by equation (10)
Figure FDA00035929649700000317
And calculating μ (k) by equation (9); and then the optimized control variable is
Figure FDA00035929649700000318
7. Method according to claim 6, 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 structuresU、usL(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 FDA0003592964970000041
Namely, it is
Figure FDA0003592964970000042
μ (k) in equation (9) may be re-described as:
Figure FDA0003592964970000043
in the formula (I), the compound is shown in the specification,
Figure FDA0003592964970000044
the magnitude constraint of the control quantity can therefore be expressed by the following inequality:
Figure FDA0003592964970000045
in the formula (I), the compound is shown in the specification,
Figure FDA0003592964970000046
also, the constraint on the control quantity increment can be expressed by the following inequality:
Figure FDA0003592964970000047
in the formula (I), the compound is shown in the specification,
Figure FDA0003592964970000048
and
Figure FDA0003592964970000049
are respectively as
Figure FDA00035929649700000410
And
Figure FDA00035929649700000411
an increment of (d);
Figure FDA00035929649700000412
8. an electronic device comprising a memory and a processor, the memory storing a computer program, wherein the processor implements the steps of the method according to any one of claims 1-7 when executing the computer program.
9. 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 to 7.
CN202210380956.3A 2022-04-12 2022-04-12 Disturbance speed compensation method, equipment and medium for wind power pile boarding trestle bridge Active CN114735140B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210380956.3A CN114735140B (en) 2022-04-12 2022-04-12 Disturbance speed compensation method, equipment and medium for wind power pile boarding trestle bridge

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210380956.3A CN114735140B (en) 2022-04-12 2022-04-12 Disturbance speed compensation method, equipment and medium for wind power pile boarding trestle bridge

Publications (2)

Publication Number Publication Date
CN114735140A true CN114735140A (en) 2022-07-12
CN114735140B 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 Disturbance speed compensation method, equipment and medium for wind power pile boarding trestle bridge

Country Status (1)

Country Link
CN (1) CN114735140B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116184841A (en) * 2023-04-28 2023-05-30 山东大学 A Model Predictive Control Method of Offshore Trestle Based on Extremum Search Algorithm
CN117163219A (en) * 2023-09-15 2023-12-05 哈尔滨理工大学 Shipborne trestle feedforward fuzzy control method considering constraint between long rods

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100230370A1 (en) * 2008-05-21 2010-09-16 Klaus Schneider Crane control with active heave compensation
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
US20180244505A1 (en) * 2017-02-28 2018-08-30 J. Ray Mcdermott S.A. Offshore ship-to-ship lifting with target tracking assistance
CN111045332A (en) * 2019-12-27 2020-04-21 哈尔滨工程大学 An Unmanned Vehicle 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 大连海事大学 An offshore transfer trestle bridge wave compensation control system and its working method
CN113431814A (en) * 2021-06-17 2021-09-24 江苏科技大学 Synchronous control method based on parallel motion of multiple hydraulic cylinders of heave compensation platform

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100230370A1 (en) * 2008-05-21 2010-09-16 Klaus Schneider Crane control with active heave compensation
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
US20160144932A1 (en) * 2013-07-19 2016-05-26 Icd Software As 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
US20180244505A1 (en) * 2017-02-28 2018-08-30 J. Ray Mcdermott S.A. Offshore ship-to-ship lifting with target tracking assistance
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 哈尔滨工程大学 An Unmanned Vehicle 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 大连海事大学 An offshore transfer trestle bridge wave compensation control system and its working method
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
苏长青等: "一种用于登乘栈桥的主动波浪补偿方法", 《船舶与海洋工程》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116184841A (en) * 2023-04-28 2023-05-30 山东大学 A Model Predictive Control Method of Offshore Trestle Based on Extremum Search Algorithm
CN117163219A (en) * 2023-09-15 2023-12-05 哈尔滨理工大学 Shipborne trestle feedforward fuzzy control method considering constraint between long rods
CN117163219B (en) * 2023-09-15 2024-05-14 哈尔滨理工大学 A feedforward fuzzy control method for shipborne trestle considering the constraints between rod length and rod interval

Also Published As

Publication number Publication date
CN114735140B (en) 2023-03-10

Similar Documents

Publication Publication Date Title
CN114735140A (en) Method, equipment and medium for compensating disturbance speed of wind power pile boarding trestle
CN107434010B (en) An electric wave active compensation boarding system and its control method
CN112405497B (en) A hybrid mechanism system based on passive compensation and its motion decomposition method
WO2012110173A2 (en) Blade load reduction for wind turbines
CN115826606B (en) Dynamic positioning control method for self-elevating ship platform
CN108459504A (en) The cooperative self-adapted iterative learning control method of multipoint mooring
CN111665722A (en) Sea wave active compensation control method based on predictive control
CN113387276A (en) Ship crane control method for improving LQR
CN117163219B (en) A feedforward fuzzy control method for shipborne trestle considering the constraints between rod length and rod interval
Borg et al. Elastic deformations of floaters for offshore wind turbines: dynamic modelling and sectional load calculations
Liu et al. Dynamic response of a SWATH vessel for installing pre-assembled floating wind turbines
CN110422284A (en) A kind of active compensation method and system based on Ship Motion forecast
Meng et al. Experimental study of floating wind turbine control on a TetraSub floater with tower velocity feedback gain
CN108303870B (en) A kind of control method, the control system of multiple degrees of freedom change in displacement compensating platform
Zhang et al. Wave heave compensation based on an optimized backstepping control method
Xu et al. Review of heave compensation systems: Design and control strategies
CN114940237A (en) Control method for offshore platform heave compensation and tensioner device thereof
CN119165763A (en) Active control method of tension tendon
Qiao et al. Motion tracking control of actuator in the active-truncated model test of deep-water mooring system
CN115639751B (en) Active compensation method for hysteresis response of hybrid platform
Jaouen et al. How does barge-master compensate for the barge motions: experimental and numerical study
CN111994201A (en) Compensation method based on sea wave compensation device
CN114879504A (en) An adaptive nonlinear control method for a four-degree-of-freedom marine rotary crane
CN111170176B (en) Active stabilization control method suitable for offshore or onshore load hoisting
Zhang et al. Control Performance Analysis of Mining Ship Heave Compensation System Based on Fuzzy Logic Algorithm.

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