CN112068437A - Single-step prediction control anti-rolling method for high-speed multi-hull ship - Google Patents

Single-step prediction control anti-rolling method for high-speed multi-hull ship Download PDF

Info

Publication number
CN112068437A
CN112068437A CN202010978084.1A CN202010978084A CN112068437A CN 112068437 A CN112068437 A CN 112068437A CN 202010978084 A CN202010978084 A CN 202010978084A CN 112068437 A CN112068437 A CN 112068437A
Authority
CN
China
Prior art keywords
heave
control
pitch
motion
channel
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
CN202010978084.1A
Other languages
Chinese (zh)
Other versions
CN112068437B (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 Engineering University
Original Assignee
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 Engineering University filed Critical Harbin Engineering University
Priority to CN202010978084.1A priority Critical patent/CN112068437B/en
Publication of CN112068437A publication Critical patent/CN112068437A/en
Application granted granted Critical
Publication of CN112068437B publication Critical patent/CN112068437B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/048Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators using a predictor
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0875Control of attitude, i.e. control of roll, pitch, or yaw specially adapted to water vehicles

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Evolutionary Computation (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Feedback Control In General (AREA)

Abstract

The invention provides a single-step prediction control stabilizing method for a high-speed multi-hull ship, which comprises the following steps: establishing a vertical motion model of multi-hull ship coupling; performing mathematical transformation on the coupled vertical motion model to obtain a single-input single-output decoupling heave and pitch motion model; designing a finite time observer, and estimating the coupling term of the heave motion and the pitch motion of the multi-hull vessel on line to obtain the anti-rolling compensation quantity; designing a single-step prediction control law, and solving the stabilization feedback control quantity of heave motion and pitch motion; and adding the compensation quantity and the feedback control quantity, and performing inverse operation on a control distribution matrix of the roll reducing attachment to obtain the control attack angle of the T-shaped wing and the wave pressing plate. The invention obviously reduces the heave and pitch of the multi-hull vessel in motion. The method realizes the application of the predictive control method with the finite time extended observer to the longitudinal rolling of the high-speed multi-hull ship. The multi-hull ship is suitable for multi-interference and high-speed marine navigation environments, the influence of sea waves on the stability of the multi-hull ship can be effectively reduced, and the performance of the multi-hull ship is improved.

Description

Single-step prediction control anti-rolling method for high-speed multi-hull ship
Technical Field
The invention relates to a single-step prediction control stabilization method, in particular to a single-step prediction control stabilization method of a high-speed multi-hull vessel, and belongs to the field of high-speed multi-hull vessel longitudinal motion stabilization control.
Background
The high-speed multi-hull ship has good transverse stability, wave resistance and maneuverability and wide application range. However, due to the unique linearity and structure, the pitching and heaving amplitudes are too large during navigation, and the performance of the vehicle is seriously influenced, so that the suppression of the heaving and pitching is particularly important.
At present, a multi-hull ship mostly takes T-shaped wings and wave pressing plates as anti-rolling attachments, and a reasonable control strategy is designed for anti-rolling control.
The multi-hull ship is influenced by different sea conditions and ship motion states encountered in actual navigation, the motion of the multi-hull ship is complex, and the motion model has the characteristics of strong coupling, time variation, uncertainty, nonlinearity and the like.
In the traditional stabilization control, the pitching channel and the heaving channel are separated to further control, and a PD stabilization control method is designed, but a large amount of time is needed for off-line parameter debugging, the robustness is weak, and the stabilization effect is general.
On the other hand, the optimal control robustness of the modern control theory is poor, buffeting phenomenon exists in sliding mode control, and therefore the anti-rolling performance is limited, the self-adaptive capacity is poor, and the application range is limited.
The predictive control breaks through the limitation of other robust control ideas, comprehensively utilizes the historical information and the model information of the controlled object, continuously performs rolling optimization on the multi-target function, and outputs and corrects the predictive model according to the actually measured object, so that the control effect is improved, and the conservatism is reduced.
Disclosure of Invention
The invention aims to provide a single-step predictive control anti-rolling method of a high-speed multi-hull vessel for improving the stability of longitudinal motion of the high-speed multi-hull vessel.
The purpose of the invention is realized as follows:
a single-step predictive control anti-rolling method for a high-speed multi-hull ship comprises the following steps:
step one, establishing a vertical motion model of multi-hull ship coupling;
step two, carrying out mathematical transformation on the coupled vertical motion model to obtain a single-input single-output decoupling heave and pitch motion model;
designing a finite time observer, and estimating the coupling terms of the heave motion and the pitch motion of the multi-hull vessel on line to obtain a roll reduction compensation quantity;
designing a single-step prediction control law, and solving the stabilization feedback control quantity of the heave motion and the pitch motion;
and step five, adding the compensation quantity and the feedback control quantity, and obtaining the control attack angle of the T-shaped wing and the wave suppression plate through inverse operation of a control distribution matrix of the roll reduction attachment.
The invention also includes such features:
the longitudinal motion model of the multi-hull vessel in the step one is as follows:
obtaining a multihull ship pitching and heaving motion equation set by the Dalnbell theorem:
Figure BDA0002686553090000021
Figure BDA0002686553090000022
wherein: m is the mass of the multihull vessel; i is55Is the moment of inertia of the multihull vessel about the y-axis; a is33,a55Additional mass and additional moment of inertia for the multihull vessel; b33,b55The damping coefficient of the system; c. C33,c55Is the coefficient of restitution force of the system; a is35,a53,b35,b53,c35,c53Coupling term coefficient of force and moment; x is the number of3,x5Respectively representing heave displacement and pitch angle;
Figure BDA0002686553090000023
individual watchShowing heave speed and pitch angular speed;
Figure BDA0002686553090000024
respectively representing heave acceleration and pitch angular acceleration; fT-foil,MT-foilRespectively representing the lift force and the lifting moment of the T-shaped hydrofoil; fflap,MflapRespectively representing the force and moment provided by the press corrugated plate; fwave,MwaveRespectively representing the wave disturbance force and moment.
The second step is specifically as follows:
decoupling model of heave movement: order to
Figure BDA0002686553090000025
Representing heave displacement and heave velocity, respectively, for a decoupled heave channel model expressed as the following form with single input and single output:
Figure BDA0002686553090000026
wherein
Figure BDA0002686553090000027
The amount of motion of the pitch channel coupled to the heave channel is used as a comprehensive uncertainty item of the decoupling heave channel; input force F ═ FT-foil+Fflap,FT-foilRespectively representing the lift of a T-shaped hydrofoil, FflapRespectively representing the force provided by the press plates, FwaveRespectively represent the disturbance force of the sea waves,
Figure BDA0002686553090000028
is the gain value;
a pitching motion decoupling model: order to
Figure BDA0002686553090000029
Respectively representing pitch angle and pitch angular velocity
Figure BDA00026865530900000210
Wherein
Figure BDA0002686553090000031
The amount of motion of a heave channel coupled to a pitch channel is used as a comprehensive uncertainty item of a decoupling pitch motion model; input torque M-MT-foil+Mflap
Figure BDA0002686553090000032
Is the gain value.
The third step is specifically as follows: step three, regarding the coupling quantity of the pitching and the heaving as an uncertain quantity, and establishing a finite time observer of the decoupled pitching and heaving channel:
introduction of expansion state x in decoupled heave channel3=Δf1Memory for recording
Figure BDA0002686553090000033
Expand into a new system
Figure BDA0002686553090000034
Designing an extended state observer with non-uniform convergence in a finite time, i.e.
Figure BDA0002686553090000035
In the formula: e.g. of the type1=z1-x1;0<a<1;knMore than 0(n is 1,2,3) is an adjustable parameter; z is a radical of1,z2Is the system state x1,x2Estimate of z3Is a systematic lumped interference estimate, sig (-) is a sign function, letm=zm-xmThe coupling control compensation of (m ═ 1,2,3) heave channel is configured as follows:
uc-heave=-z3·(m+a33)
introducing an expansion state to the pitch motion model
Figure BDA0002686553090000036
Introducing an expansion state x to the pitch motion model33=Δf2To and from
Figure DA00026865530957996
Can be expanded into a new system
Figure BDA0002686553090000037
Designed corresponding finite time extended state observer
Figure BDA0002686553090000041
In the formula: e.g. of the type11=z11-x11;0<<1;knn> 0 (n-1, 2,3) is an adjustable parameter, z11,z22Is the system state x11,x22Estimated value of z33Estimating system lumped interference; sig (·) is a sign function;
the coupling control compensation of the pitching channel is configured as follows:
uc-pitch=-z33.(I55+a55)。
the fourth step is specifically as follows: adopting a single-step prediction controller as a roll-reducing attached body control system of the multi-hull vessel, taking the heave and pitch of the multi-hull vessel as control objects, and controlling the control objects by using the prediction controller;
predicted control quantity of heave channel
Figure BDA0002686553090000042
Figure BDA0002686553090000043
Where 0.05, 0.255 and 1 × 10-6Then the control quantity of heave channel feedback
Figure BDA0002686553090000044
The feedback quantity u of the pitching channel can be obtained in the same way2:
Figure BDA0002686553090000045
Figure BDA0002686553090000046
Where 0.05, 0.255 and 1 × 10-6
The fifth step is specifically as follows: adding the compensation quantity and the feedback control quantity, and obtaining a control attack angle of the T-shaped wing and the wave pressing plate through inverse operation of a control distribution matrix of the anti-rolling appendage:
step 1, combining the calculated estimated coupling compensation amount uc-heaveVirtual comprehensive control quantity can be obtained;
U1=u1+uc-heave
step 2, obtaining the virtual comprehensive control quantity U of the pitching channel by the same method2=u2+uc-pitch
Step 3, resolving the virtual control quantity through a stabilizing appendage control distribution matrix to obtain the actual control attack angle of the T-shaped wing and the wave pressing plate
Figure BDA0002686553090000051
Figure BDA0002686553090000052
Wherein alpha is1Angle of attack, α, for press-formed sheets2Of T-shaped wingAngle of attack, fflapForce generated for pressing the corrugated board, fT-foilForce generated by the T-shaped wing, rho is sea water density, A is T-shaped hydrofoil area, CLIs the coefficient of lift of the hydrofoil, V is the velocity of the fluid relative to the hydrofoil, CL1Is the lift coefficient of the press wave plate, S is the effective area of the press wave plate, lflap、lT-foilRespectively is the arm of force of pressing unrestrained board and T type wing.
Compared with the prior art, the invention has the beneficial effects that:
the control method designed by the invention aims at longitudinal rolling reduction in sailing of the most common high-speed multi-hull ship. Obviously reducing the heave and pitch of the multi-hull vessel in motion. The method realizes the application of the predictive control method with the finite time extended observer to the longitudinal rolling of the high-speed multi-hull ship. The multi-hull ship navigation device is suitable for multi-interference and high-ship-speed marine navigation environments, the influence of sea waves on the stability of the multi-hull ship can be effectively reduced, and the performance of the multi-hull ship is improved.
The invention designs a new single-step prediction control method to improve the anti-rolling performance, designs a finite time observer to improve the robustness and comprehensively improves the heave and pitch inhibition capability of the high-speed multi-hull vessel.
Drawings
FIG. 1 is a wave force disturbance diagram;
FIG. 2 is a wave moment disturbance diagram;
FIG. 3 is a graph of the pitch motion of a multihull vessel;
FIG. 4 is a heave motion profile for a multihull vessel;
FIG. 5 is a diagram of the actual and estimated states of heave motion;
fig. 6 is a diagram of the actual state and the estimated state of the pitching motion.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
The invention aims to provide a single-step prediction control anti-rolling method of a high-speed multi-hull vessel, which can improve the stability of longitudinal motion of the high-speed multi-hull vessel.
A single-step prediction control anti-rolling method for a high-speed multi-hull ship specifically comprises the following steps:
s1, establishing a multi-hull ship coupled vertical motion model;
s2, performing mathematical transformation on the coupled vertical motion model to obtain a single-input single-output decoupling heave and pitch motion model;
s3, designing a finite time observer, and estimating the coupling term of the heave motion and the pitch motion of the multi-hull vessel on line to obtain the anti-rolling compensation amount;
s4, designing a single-step prediction control law, and solving the stabilization feedback control quantity of heave motion and pitch motion;
and S5, adding the compensation quantity and the feedback control quantity, and obtaining the control attack angle of the T-shaped wing and the wave suppression plate through the inverse operation of the control distribution matrix of the roll reduction attachment.
The single-step prediction control anti-rolling method for the high-speed multi-hull ship comprises the steps of performing mathematical transformation on two coupled vertical motion models to obtain a single-input single-output decoupling heave and pitch motion model;
designing a finite time observer, and estimating a coupling term of heave motion and pitch motion of the multi-hull vessel on line to obtain a roll reduction compensation quantity;
the single-step predictive control anti-rolling method for the high-speed multi-hull vessel comprises the following steps of designing a single-step predictive control law in the fourth step, and solving anti-rolling feedback control quantity of heave motion and pitch motion;
the single-step prediction control stabilization method for the high-speed multi-hull vessel comprises the fifth step of adding the compensation quantity and the feedback control quantity, and obtaining the control attack angle of the T-shaped wing and the wave suppression plate through inverse operation of the control distribution matrix of the stabilization appendage.
As shown in fig. 3, a longitudinal roll-reducing decoupling control scheme diagram of a multi-hull vessel is shown, and the specific flow is as follows:
step one, establishing a vertical motion model of multi-hull ship coupling;
due to the special hull construction of a multi-hulled vessel, the longitudinal motion of the multi-hulled vessel is not coupled with the transverse motion. Only heave and pitch motion need to be considered, and the T-shaped wings and the wave pressing plates generate restoring force and restoring moment to offset disturbing force and disturbing moment generated by sea waves. Assuming that the multihull vessel is sailing in an infinitely deep water area with a stable course and a constant speed, the sheets of the underwater part are sufficiently slender, and the motion of the hull caused by wave disturbance is slightly radial, without considering the influence of wind and flow on the motion. Under the action of sea wave disturbance, the system of the heave and pitch coupled motion equation of the multi-hull ship can be obtained as follows:
Figure BDA0002686553090000061
Figure BDA0002686553090000062
wherein: m is the mass of the multihull vessel; i is55Is the moment of inertia of the multihull vessel about the y-axis; a isiiAdditional mass and additional moment of inertia for the multihull vessel; biiThe damping coefficient of the system; c. CiiIs the coefficient of restitution force of the system; a isij,bij,cijCoupling term coefficient of force and moment; x is the number of3,x5Respectively representing heave displacement and pitch angle;
Figure BDA0002686553090000063
respectively representing heave velocity and pitch angular velocity;
Figure BDA0002686553090000064
respectively representing heave acceleration and pitch angular acceleration; fT-foil,MT-foilRespectively representing the restoring force and the restoring moment of the T-shaped wing; fflap,MflapRespectively representing the restoring force and the restoring moment provided by the press wave plate; fwave,MwaveRespectively representing the sea wave disturbance force and moment.
Step two, decoupling the longitudinal motion model to obtain decoupled heave motion and pitch motion models;
decoupling model of heave movement: order to
Figure BDA0002686553090000071
Representing heave displacement and heave velocity separately, for decoupled heave channel models as follows
Decoupling model of heave movement: order to
Figure BDA0002686553090000072
Representing heave displacement and heave velocity separately, expressed as the following form of input single output for decoupled heave channel model
Figure BDA0002686553090000073
Wherein
Figure BDA0002686553090000074
The amount of motion of the pitch channel coupled to the heave channel is used as a comprehensive uncertainty item of the decoupling heave channel; input force F ═ FT-foil+Fflap,FT-foilRespectively representing the lift of a T-shaped hydrofoil, FflapRespectively representing the force provided by the press plates, FwaveRespectively represent the disturbance force of the sea waves,
Figure BDA0002686553090000075
is the gain value.
Order to
Figure BDA0002686553090000076
Respectively representing pitch angle and pitch angular velocity
Figure BDA0002686553090000077
Wherein
Figure BDA0002686553090000078
Is the amount of motion that the heave channel couples to the pitch channel as a comprehensive uncertainty for the decoupled pitch motion model. Input torque M-MT-foil+Mflap
Figure BDA0002686553090000079
Is the gain value.
Step three, regarding the coupling quantity of the pitching and the heaving as an uncertain quantity, and establishing a finite time observer of the decoupled pitching and heaving channel;
regarding the decoupled heave motion model, considering the amount of motion of pitch coupled to the heave channel and the amount of external random disturbance as uncertain quantities, introducing a dilated state in the decoupled heave channel
Figure BDA00026865530900000710
Note the book
Figure BDA00026865530900000711
Expand into a new system
Figure BDA0002686553090000081
An extended state observer with non-uniform convergence in a finite time is designed, i.e.
Figure BDA0002686553090000082
In the formula: e.g. of the type1=z1-x1;0<a<1;knMore than 0(n is 1,2,3) is an adjustable parameter; z is a radical of1,z2Is the system state x1,x2Estimate of z3Is a systematic lumped interference estimate, sig (·) is a sign function. Note em=zm-xm(m=1,2,3),
The coupling control compensation of the heave channel is configured as follows:
uc-heave=-z3(m+a33) (7)
for the decoupling pitch motion model, the motion quantity of heave coupled to the pitch channel and the external random disturbance quantity are regarded as uncertain quantities, and an expansion state is introduced into the decoupling pitch channel
Figure BDA0002686553090000083
Note the book
Figure BDA0002686553090000084
Expand into a new system
Figure BDA0002686553090000085
Designed corresponding finite time extended state observer
Figure BDA0002686553090000086
In the formula: e.g. of the type11=z11-x11;0<<1;knn> 0 (n-1, 2,3) is an adjustable parameter, z11,z22Is the system state x11,x22Estimated value of z33Is a systematic lumped interference estimate, sig (·) is a sign function.
The coupling control compensation of the pitching channel is configured as follows:
uc-pitch=-z33.(I55+a55) (10)
designing a single-step prediction control law, and solving the stabilization feedback control quantity of the heave motion and the pitch motion; predicted control quantity of heave channel
Figure BDA0002686553090000091
Figure BDA0002686553090000092
Where 0.05, 0.255 and 1 × 10-6Then the control quantity of heave channel feedback
Figure BDA0002686553090000093
The feedback quantity u of the pitching channel can be obtained in the same way2:
Figure BDA0002686553090000094
Figure BDA0002686553090000095
Where 0.05, 0.255 and 1 × 10-6
And step five, adding the compensation quantity and the feedback control quantity, and obtaining the control attack angle of the T-shaped wing and the wave suppression plate through inverse operation of a control distribution matrix of the roll reduction attachment.
Combining the calculated estimated coupling compensation uc-heave,uc-pitchVirtual integrated control quantity can be obtained
U1=u1+uc-heave (14)
U2=u2+uc-pitch (15)
By utilizing the control distribution matrix of the multi-hull ship anti-rolling appendage, the actual control attack angle of the T-shaped wing and the wave pressing plate can be obtained as follows:
Figure BDA0002686553090000096
Figure BDA0002686553090000097
wherein alpha is1Angle of attack, α, for press-formed sheets2Angle of attack, f, for T-shaped wingflapForce generated for pressing the corrugated board, fT-foilForce generated by the T-shaped wing, rho is sea water density, A is T-shaped hydrofoil area, CLIs the coefficient of lift of the hydrofoil, V is the velocity of the fluid relative to the hydrofoil, CL1Is the lift coefficient of the press wave plate, S is the effective area of the press wave plate, lflap、lT-foilRespectively is the arm of force of pressing unrestrained board and T type wing.
Simulation and experiments show that when the multi-hull ship sails at 14 knots and is subjected to wave-facing sailing, the encounter frequency is 1.3rad/s, the course angle is 180 degrees, the hull waterline is usually 48m, and the sea wave grade is 4. Fig. 1 and 2 show the disturbance of sea wave disturbance force and disturbance moment, and fig. 4 shows that the pitching amount of the multi-hull vessel is reduced by 50-60% when a controller is arranged, compared with the controller without the controller; FIG. 5 can see a 40% -50% reduction in heave by the controller compared to no controller; the stability of the system can be obviously improved, and the performance of the multi-hull ship is improved. In addition, as can be seen from fig. 6, the observer can well estimate the state of the system, and can be popularized and applied to control of other multi-hull vessels.
The above embodiments are only for illustrating the design idea and features of the present invention, and the purpose of the present invention is to enable those skilled in the art to understand the content of the present invention and implement the present invention, and the protection scope of the present invention is not limited to the above embodiments. Therefore, all the equivalent changes or modifications made according to the principles and design ideas disclosed by the present invention are within the protection scope of the present invention.
In summary, the following steps: the invention discloses a single-step prediction control anti-rolling method for a high-speed multi-hull ship. And establishing a single-input single-output decoupling model of the heaving and pitching motions of the multi-hull vessel, and converting the heaving and pitching coupled motions into a comprehensive uncertain item. Aiming at decoupled heave motion and pitch motion models, designing a finite time observer to estimate a time-varying coupling term on line and giving a real-time compensation amount; and designing a new single-step prediction control law to give feedback control quantity. And (4) integrating the feedback control quantity and the compensation quantity, obtaining the attack angle of the T-shaped wing and the wave pressing plate through matrix inversion, and automatically adjusting the heave height and pitch angle change of the multi-hull ship. The control method provided by the invention is simple, easy to realize, strong in self-adaptation, and capable of effectively reducing the heaving amount of the multi-hull ship by 40-50% and the pitching amount by 50-60%.

Claims (6)

1. A single-step prediction control anti-rolling method of a high-speed multi-hull ship is characterized by comprising the following steps:
step one, establishing a vertical motion model of multi-hull ship coupling;
step two, carrying out mathematical transformation on the coupled vertical motion model to obtain a single-input single-output decoupling heave and pitch motion model;
designing a finite time observer, and estimating the coupling term of the heave motion and the pitch motion of the multi-hull vessel on line to obtain the anti-rolling compensation quantity;
designing a single-step prediction control law, and solving the stabilization feedback control quantity of the heave motion and the pitch motion;
and step five, adding the compensation quantity and the feedback control quantity, and obtaining the control attack angle of the T-shaped wing and the wave suppression plate through inverse operation of a control distribution matrix of the roll reduction attachment.
2. The single-step predictive control roll-reducing method for a high-speed multi-hull vessel according to claim 1, wherein the longitudinal motion model of the multi-hull vessel in the first step is:
obtaining a multihull ship pitching and heaving motion equation set by the Dalnbell theorem:
Figure FDA0002686553080000011
Figure FDA0002686553080000012
wherein: m is the mass of the multihull vessel; i is55Is the moment of inertia of the multihull vessel about the y-axis; a is33,a55Additional mass and additional moment of inertia for the multihull vessel; b33,b55The damping coefficient of the system; c. C33,c55Is the coefficient of restitution force of the system; a is35,a53,b35,b53,c35,c53Coupling term coefficient of force and moment; x is the number of3,x5Respectively representing heave displacement and pitch angle;
Figure FDA0002686553080000015
respectively representing the heave speedAnd pitch angular velocity;
Figure FDA0002686553080000014
respectively representing heave acceleration and pitch angular acceleration; fT-foil,MT-foilRespectively representing the lift force and the lifting moment of the T-shaped hydrofoil; fflap,MflapRespectively representing the force and moment provided by the press corrugated plate; fwave,MwaveRepresenting wave disturbance forces and moments, respectively.
3. The single-step predictive control roll-reducing method for a high-speed multi-hull vessel according to claim 1, wherein the second step is specifically:
decoupling model of heave movement: let x1=x3,
Figure FDA0002686553080000016
Heave displacement and heave velocity are represented separately, and for the decoupled heave channel model, the following form is represented for input single output:
Figure FDA0002686553080000013
wherein
Figure FDA0002686553080000021
The amount of motion of the pitch channel coupled to the heave channel is used as a comprehensive uncertainty item of the decoupling heave channel; input force F ═ FT-foil+Fflap,FT-foilRespectively representing the lift of a T-shaped hydrofoil, FflapRespectively representing the force provided by the press plates, FwaveRespectively represent the disturbance force of the sea waves,
Figure FDA0002686553080000022
is the gain value;
a pitching motion decoupling model: order to
Figure FDA0002686553080000029
Respectively representing pitch angle and pitch angular velocity
Figure FDA0002686553080000023
Wherein
Figure FDA0002686553080000024
The amount of motion of a heave channel coupled to a pitch channel is used as a comprehensive uncertainty item of a decoupling pitch motion model; input torque M-MT-foil+Mflap
Figure FDA0002686553080000025
Is the gain value.
4. The single-step predictive control roll-reducing method for a high-speed multi-hull vessel according to claim 1, wherein the third step is specifically: step three, regarding the coupling quantity of the pitching and the heaving as an uncertain quantity, and establishing a finite time observer of the decoupled pitching and heaving channel:
introduction of expansion state x in decoupled heave channel3=Δf1Memory for recording
Figure FDA0002686553080000028
Expand into a new system
Figure FDA0002686553080000026
Designing an extended state observer with non-uniform convergence in a finite time, i.e.
Figure FDA0002686553080000027
In the formula: e.g. of the type1=z1-x1;0<a<1;knMore than 0(n is 1,2,3) is adjustable ginsengCounting; z is a radical of1,z2Is the system state x1,x2Estimate of z3Is a systematic lumped interference estimate, sig (-) is a sign function, letm=zm-xm(m=1,2,3)
The coupling control compensation of the heave channel is configured as follows:
uc-heave=-z3·(m+a33)
introducing an expansion state to the pitch motion model
Figure FDA0002686553080000031
Introducing an expansion state x to the pitch motion model33=Δf2To and from
Figure FDA0002686553080000038
Can be expanded into a new system
Figure FDA0002686553080000032
Designed corresponding finite time extended state observer
Figure FDA0002686553080000033
In the formula: e.g. of the type11=z11-x11;0<<1;knn> 0 (n-1, 2,3) is an adjustable parameter, z11,z22Is the system state x11,x22Estimate of z33Estimating system lumped interference; sig (·) is a sign function;
the coupling control compensation of the pitching channel is configured as follows:
uc-pitch=-z33.(I55+a55)。
5. the single-step predictive control roll-reducing method for a high-speed multi-hull vessel according to claim 1, wherein said step four is specifically: adopting a single-step prediction controller as a roll-reducing attached body control system of the multi-hull vessel, taking the heave and pitch of the multi-hull vessel as control objects, and controlling the control objects by using the prediction controller;
predicted control quantity of heave channel
Figure FDA0002686553080000034
Figure FDA0002686553080000035
Where 0.05, 0.255 and 1 × 10-6Then the control quantity of heave channel feedback
Figure FDA0002686553080000036
The feedback quantity u of the pitching channel can be obtained in the same way2:
Figure FDA0002686553080000037
Figure FDA0002686553080000041
Where 0.05, 0.255 and 1 × 10-6
6. The single-step predictive control roll-reducing method for a high-speed multi-hull vessel according to claim 1, wherein said step five is specifically: adding the compensation quantity and the feedback control quantity, and obtaining a control attack angle of the T-shaped wing and the wave pressing plate through inverse operation of a control distribution matrix of the anti-rolling appendage:
step 1, combining the calculated estimated coupling compensation amount uc-heaveVirtual comprehensive control quantity can be obtained;
U1=u1+uc-heave
step 2, obtaining the virtual comprehensive control quantity U of the pitching channel by the same method2=u2+uc-pitch
Step 3, resolving the virtual control quantity through a stabilizing appendage control distribution matrix to obtain the actual control attack angle of the T-shaped wing and the wave pressing plate
Figure FDA0002686553080000042
Figure FDA0002686553080000043
Wherein alpha is1Angle of attack, α, for press-formed sheets2Angle of attack, f, for T-shaped wingflapForce generated for pressing the corrugated board, fT-foilForce generated by the T-shaped wing, rho is sea water density, A is T-shaped hydrofoil area, CLIs the coefficient of lift of the hydrofoil, V is the velocity of the fluid relative to the hydrofoil, CL1Is the lift coefficient of the press wave plate, S is the effective area of the press wave plate, lflap、lT-foilRespectively is the arm of force of pressing unrestrained board and T type wing.
CN202010978084.1A 2020-09-17 2020-09-17 Single-step prediction control anti-rolling method for high-speed multi-hull ship Active CN112068437B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010978084.1A CN112068437B (en) 2020-09-17 2020-09-17 Single-step prediction control anti-rolling method for high-speed multi-hull ship

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010978084.1A CN112068437B (en) 2020-09-17 2020-09-17 Single-step prediction control anti-rolling method for high-speed multi-hull ship

Publications (2)

Publication Number Publication Date
CN112068437A true CN112068437A (en) 2020-12-11
CN112068437B CN112068437B (en) 2022-06-21

Family

ID=73680635

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010978084.1A Active CN112068437B (en) 2020-09-17 2020-09-17 Single-step prediction control anti-rolling method for high-speed multi-hull ship

Country Status (1)

Country Link
CN (1) CN112068437B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050247251A1 (en) * 2004-05-07 2005-11-10 Scism Randy M Integrated marine performance system
CN107145074A (en) * 2017-06-27 2017-09-08 哈尔滨工程大学 A kind of high-speed trimaran pitching stabilization control method based on sliding moding structure convergence law
CN107247413A (en) * 2017-06-27 2017-10-13 哈尔滨工程大学 A kind of high-speed trimaran pitching stabilization control method based on fuzzy control
CN109334892A (en) * 2018-09-20 2019-02-15 江苏大学 A kind of simplification robust adaptive pitching stabilization control method of multi-hull ship
CN109878633A (en) * 2019-03-14 2019-06-14 哈尔滨工程大学 A kind of high-speed catamaran longitudinal direction stabilization method based on Active Disturbance Rejection Control
CN111498037A (en) * 2020-04-05 2020-08-07 哈尔滨工程大学 High-speed catamaran longitudinal stabilization method based on variable structure active disturbance rejection control

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050247251A1 (en) * 2004-05-07 2005-11-10 Scism Randy M Integrated marine performance system
CN107145074A (en) * 2017-06-27 2017-09-08 哈尔滨工程大学 A kind of high-speed trimaran pitching stabilization control method based on sliding moding structure convergence law
CN107247413A (en) * 2017-06-27 2017-10-13 哈尔滨工程大学 A kind of high-speed trimaran pitching stabilization control method based on fuzzy control
CN109334892A (en) * 2018-09-20 2019-02-15 江苏大学 A kind of simplification robust adaptive pitching stabilization control method of multi-hull ship
CN109878633A (en) * 2019-03-14 2019-06-14 哈尔滨工程大学 A kind of high-speed catamaran longitudinal direction stabilization method based on Active Disturbance Rejection Control
CN111498037A (en) * 2020-04-05 2020-08-07 哈尔滨工程大学 High-speed catamaran longitudinal stabilization method based on variable structure active disturbance rejection control

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DALEI QIAO.ETC: "Marine Vessel Re-Identification: A Large-Scale Dataset and Global-and-Local Fusion-Based Discriminative Feature Learning", 《IEEE ACCESS》 *
刘娜: "基于减摇附体的高速多体船减纵摇控制规律的研究", 《中国优秀硕士学位论文全文数据库(电子期刊)》 *

Also Published As

Publication number Publication date
CN112068437B (en) 2022-06-21

Similar Documents

Publication Publication Date Title
Matusiak Dynamics of a rigid ship
CN109334892B (en) Simplified robust self-adaptive pitching reduction control method for multi-hull vessel
CN104527943B (en) Energy optimization based integrated stabilization device dual neural network self-tuning PID (Proportion Integration Differentiation) control method
CN110333658B (en) Under-actuated ship course control method with input constraint and drift angle compensation
CN109739248B (en) Ship-borne three-degree-of-freedom parallel stable platform stability control method based on ADRC
CN116911135B (en) Nonlinear water elasticity time domain calculation and evaluation method considering slamming load under waves
CN111198502A (en) Unmanned ship track tracking control method based on interference observer and fuzzy system
CN110937076A (en) Ship comfort control system based on model prediction design of two-channel parameters and control method thereof
Mao et al. Force model of flapping foil stabilizers based on CFD parameterization
Matusiak Dynamics of a Rigid Ship-with applications
Huang et al. Numerical and experimental studies on a predictive control approach for pitch stabilization in heading waves
CN112068437B (en) Single-step prediction control anti-rolling method for high-speed multi-hull ship
CN108681245B (en) Fin stabilizer controller for solving problem of limited actual lift force
CN109398594B (en) Wave climbing control method for hydrofoil ship
CN110457827B (en) General calculation method for resistance of yacht movement mathematical model in yacht simulator
CN110008505B (en) Numerical model method for solving problem of response divergence of horizontal motion of ship in oblique waves
CN111498037B (en) High-speed catamaran longitudinal stabilization method based on variable structure active disturbance rejection control
Guan et al. Concise robust fin roll stabilizer design based on integrator backstepping and CGSA
Li et al. SMC for rudder/fin joint roll stabilization based on input-output feedback linearization
CN110837254B (en) Longitudinal motion control method for hydrofoil craft
CN113848915A (en) Polar ship course control method based on linear active disturbance rejection technology
CN113901578A (en) Distortion rudder optimization design method considering uncertainty factors
CN110096745B (en) Wave abundance depth calculation method based on wave encounter period
JP2002284087A (en) Anti-rolling device
Li et al. Ship roll stabilization using supervision control based on inverse model wavelet neural network

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