CN109334892B - Simplified robust self-adaptive pitching reduction control method for multi-hull vessel - Google Patents
Simplified robust self-adaptive pitching reduction control method for multi-hull vessel Download PDFInfo
- Publication number
- CN109334892B CN109334892B CN201811100020.0A CN201811100020A CN109334892B CN 109334892 B CN109334892 B CN 109334892B CN 201811100020 A CN201811100020 A CN 201811100020A CN 109334892 B CN109334892 B CN 109334892B
- Authority
- CN
- China
- Prior art keywords
- heave
- pitch
- motion model
- channel
- vessel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000009467 reduction Effects 0.000 title claims abstract description 18
- 238000010168 coupling process Methods 0.000 claims abstract description 14
- 238000005859 coupling reaction Methods 0.000 claims abstract description 14
- 230000008878 coupling Effects 0.000 claims abstract description 9
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 4
- 239000011888 foil Substances 0.000 claims description 18
- 230000003044 adaptive effect Effects 0.000 claims description 15
- 238000013016 damping Methods 0.000 claims description 15
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 230000001133 acceleration Effects 0.000 claims description 7
- 230000008859 change Effects 0.000 abstract description 3
- 230000009471 action Effects 0.000 description 7
- 230000006641 stabilisation Effects 0.000 description 7
- 238000011105 stabilization Methods 0.000 description 7
- 238000013461 design Methods 0.000 description 4
- 238000005096 rolling process Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000011217 control strategy Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 239000013535 sea water Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000000476 body water Anatomy 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
- B63B—SHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING
- B63B39/00—Equipment to decrease pitch, roll, or like unwanted vessel movements; Apparatus for indicating vessel attitude
- B63B39/06—Equipment to decrease pitch, roll, or like unwanted vessel movements; Apparatus for indicating vessel attitude to decrease vessel movements by using foils acting on ambient water
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive 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/042—Adaptive 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
Landscapes
- Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- Ocean & Marine Engineering (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses a simplified robust self-adaptive pitching reduction control method for a multi-hull vessel, which comprises the steps of establishing a coupled motion model of heaving and pitching of the multi-hull vessel, converting the model into a state space form, and decomposing the coupled motion model of heaving and pitching of the multi-hull vessel into a decoupled heaving motion model and a decoupled pitching motion model; respectively designing an extended state observer aiming at the decoupled heave motion model and pitch motion model, wherein the extended state observers are used for estimating the motion state and coupling terms of the multi-hull vessel; and combining the steps, synthesizing the basic control quantity of the multi-hull ship and the estimated coupling term to obtain a virtual control quantity, and distributing the virtual control quantity to the attack angle of the T-shaped wing and the wave pressing plate. The change of the heave height and the pitch angle of the multi-hull vessel is controlled by force and moment generated by an attack angle, the control method provided by the invention is simple and reliable, is easy to realize, and can effectively reduce 20-35% of the heave of the multi-hull vessel and 40-50% of the pitch.
Description
Technical Field
The invention belongs to the field of ship navigation stability control, and particularly relates to a simplified robust self-adaptive pitching reduction control method for a multi-hull ship.
Background
The multi-hull ship is an important development direction of modern high-performance ships and has good transverse stability, wave resistance and maneuverability. However, when the ship sails at a high speed, the unique line type and structure of the multi-hull ship result in that hydrodynamic force at different sailing speeds can generate longitudinal overturning moment with different degrees on the underwater hull, the moment is increased along with the increase of the pitch angle and is dynamically changed along with the sailing speed; on the other hand, the longitudinal restoring moment of the multi-hull vessel is small, so that the pitching/heaving motion amplitude of the multi-hull vessel is too large, the longitudinal motion is easy to destabilize, the navigation performance is affected by non-negligible influence, and the seaworthiness of the multi-hull vessel is seriously affected. Therefore, the pitching reduction control of the multi-hulled vessel is required.
At present, a multihull ship at home and abroad is mainly provided with two types of anti-rolling equipment, namely T-shaped wings and wave pressing plates, and a reasonable control strategy is designed to reduce the heave, the pitch and the heave of the multihull ship. The multi-hull ship motion model has the characteristics of fast time-varying property, strong nonlinearity, multiple uncertainties and the like, and the design of the controller has certain challenges. The traditional longitudinal stabilization control adopts the separation of a pitching channel and a heaving channel, and the gain scheduling proportional-differential control is designed for stabilization, but the control method needs a large amount of time for off-line parameter debugging, the robustness is weak, and the stabilization effect is general. On the other hand, although the fuzzy control and robust control improve the rolling reduction robust performance of the multi-hull vessel, the fuzzy control and robust control have the defects of over conservative controller design, poor adaptive capacity, limited application range, complex design and the like, and are difficult to apply to engineering. At present, no research is available in China on the integration of feedback control law and adaptive online compensation to improve the adaptivity, stability and the like of the anti-rolling control.
Disclosure of Invention
According to the defects and shortcomings of the prior art, the invention provides a simplified robust adaptive pitching reduction control method for a multi-hull vessel, which can improve the adaptivity and stability of pitching reduction control by combining a feedback control law and adaptive online compensation, thereby improving the stability of the multi-hull vessel.
The technical scheme adopted by the invention is as follows:
a simplified robust adaptive pitching reduction control method for a multi-hull vessel comprises the following steps:
step 4, respectively solving basic control quantities of the heave motion model and the pitch motion model by adopting a proportional-differential control method;
and 5, combining the steps, synthesizing the basic control quantity of the multi-hull ship and the estimated coupling term to obtain a virtual control quantity, and distributing the virtual control quantity to the attack angle of the T-shaped wing and the wave pressing plate.
Further, the coupled motion model of the heaving and pitching of the multi-hull vessel in the step 1 is as follows:
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;respectively representing the heave velocity and the pitch angular velocity;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 disturbance force and the moment of the sea waves;
the coupled kinematics model of heave and pitch is transformed by an equivalent mathematical transformation into the following state space form:
Further, the heave motion model is represented as:
wherein x is1=x3,x1Which is indicative of the displacement of the heave,x2representing heave velocity, the amount of pitch coupled to the heave channel as an uncertainty of the heave channel, i.e.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,is the gain value.
The pitch motion model is represented as:
wherein x is11Representing pitch angle, x22Denotes pitch angular velocity, and x11=x5,The amount of heave coupled to the pitch channel motion is used as uncertainty of the pitch channel, i.e. uncertaintyInput moment M ═ MT-foil+Mflap,MT-foilRespectively representing T-shaped hydrofoil lifting moment, MflapRespectively representing the moment provided by the press plates, MwaveRespectively represent the disturbance moment of the sea waves,is the gain value.
Further, the following extended state observer is designed for the heave motion model:
parameter βiBy adopting the configuration method based on the bandwidth, the following conditions are met:
[β1,β2,β3]=[ω0α1,ω0 2α2,ω0 3α3];
wherein, βiFor adjusting the parameters, i is 1, 2, 3, ω0Selecting a gain coefficient α corresponding to the bandwidth of the heave channel extended state observeri=3!/i!×(3-i)!,i=1,2,3,e1Is the system state x1With the estimated state z of the observer1Error of (2), z1、z2Is the system state x1、x2Corresponding observer estimated state, z3Is an estimate of the lumped interference of the system, i.e. z1→x1,z2→x2,z3→x3=f2Error g1i(e1)=e1,i=1、2、3;Respectively being estimated states z1、z2、z3First derivative of b1For the gain value, U, corresponding to the heave channel1Designing the virtual control quantity for the heave channel;
designing the following extended state observer aiming at the model of the pitching motion:
parameter βiiBy adopting the configuration method based on the bandwidth, the following conditions are met:
[β11,β22,β33]=[ω1α1,ω1 2α2,ω1 3α3]
wherein, βiiFor adjusting the parameters, i is 1, 2, 3, ω1Is the bandwidth corresponding to the pitch channel extended state observer, e2Is the system state x11With the estimated state z of the observer11Error of (2), z11、z22Is the system state x11,x22Estimate of z33Is an estimate of the lumped interference of the system, i.e. z11→x11,z22→x22,z33→x33=f22Error g2i(e2)=e2,i=1、2、3; Respectively being estimated states z11、z22、z33B is the gain value corresponding to the pitch channel, U2Designing the obtained virtual control quantity for the pitching channel;
further, the ratio term kp of the heave motion model is obtained1And a derivative term kd1The method comprises the following steps:
wherein, wnLet epsilon be damping ratio, take epsilon 0.85, m be the mass of the multihull vessel, a33As an additional mass of the multi-hulled vessel, c33Is the coefficient of restitution of the system, K1The magnitude of the force generated by the press wave plate;
the method for solving the proportional term and the differential term of the pitching motion model comprises the following steps:
further, the pitching channel is solved to obtain a proportional term kp2And a derivative term kd2The method comprises the following steps:
wherein, wnTaking epsilon as a damping ratio, I as the natural frequency of the system, and taking epsilon as 0.855Is the moment of inertia of the multihull vessel about the y-axis, a55To add moment of inertia, c55Is the coefficient of restitution of the system, K2The magnitude of the moment generated by the T-shaped wing;
further, the method for obtaining the comprehensive control quantity of the virtual control law and the coupling term of the multi-hull vessel comprises the following steps:
step 5.1, combining the proportional term and the differential term of the heave motion model and the pitch motion model to respectively obtain the basic control input β of the heave motion model and the pitch motion model1、β2;
Step 5.2, two basic control input comprehensive disturbance lumped interference estimated values z obtained based on the heave and pitch motion models3And z33To determine the final virtual control quantity U of the heave motion model and the pitch motion model1、U2;
Step 5.3, according to the virtual control quantity U of the heave motion model and the pitch motion model1、U2And obtaining the comprehensive control quantity of the multi-hull ship, namely the attack angle input of the T-shaped wing and the wave pressing plate.
Further, the basic control inputs β of the heave motion model and the pitch motion model in said step 5.11、β2Expressed as: the basic control quantities of the heave channel are:
the basic control quantities for the pitch channel are:
wherein x is1For heave displacement, x2For the heave velocity, a33、a55For additional mass and additional moment of inertia of the multihull vessel, b33、b55Is the damping coefficient of the system, c33、c55The coefficient of the restoring force of the system, m is the mass of the multi-hull ship, and the heave motion model proportion term kp1And a derivative term kd1,I55Is the moment of inertia, x, of the multihull vessel about the y-axis11Pitch angle, x22Is the pitch angular velocity; scaling term kp of multi-hull ship pitching channel2And a derivative term kd2。
Further, the total virtual control amount is:
U1=β1-z3/b1;
U2=β2-z33/b2;
wherein, β1As a basic control input for the heave channel, β2As a basic control input for the pitch channel, z3(t)、z33(t) is the lumped disturbance estimate for the heave and pitch channels,
further, the method for obtaining the attack angle input of the T-shaped wing and the wave pressing plate comprises the following steps:
wherein, α1Angle of attack for pressed plates, α2Angle 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.
The invention has the beneficial effects that: the method realizes the robust adaptive stabilization control of heave and pitch stabilization of the multi-hull vessel and improves the comfort of the multi-hull vessel by the feedback stabilization of a decoupled proportional-differential control law and the online feedforward compensation of the adaptive extended observer.
In the invention, the designed observable states of the feedback control law and the feedforward compensation are a pitch angle, an angular velocity, a heave height and a heave speed, the synthesized total control quantity is calculated into an attack angle of the T-shaped wing and the wave pressing plate, and the change of the heave height and the pitch angle of the multi-hull ship is controlled through the force and the moment generated by the attack angle. By adopting the control algorithm, the control law is simple and reliable, the programming is easy, the engineering is easy to realize, and the heaving and swaying of the multi-hull ship are reduced by 20-35%, and the pitching and swaying are reduced by 40-50%.
Drawings
FIG. 1 is a control schematic of a multi-hulled vessel;
FIG. 2 is a heave displacement and observer estimation value of a multi-hull vessel under the action of a controller under the interference of sea waves and a state quantity curve without the action of the controller;
FIG. 3 is a heave velocity and observer estimation value of a multi-hull vessel under the action of sea wave after a controller is added and a state quantity curve without the action of the controller;
FIG. 4 is a curve of a pitch angle and an observer estimation value of a multi-hull vessel under the action of sea waves after a controller is added and a state quantity without the action of the controller;
FIG. 5 is a curve of the pitch angular velocity and observer estimation value of a multihull vessel under the action of sea wave disturbance with the addition of a controller and state quantities without the addition of the controller.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The method adopts a simplified feedforward and feedback composite control method to robustly pitch and roll the multi-hull vessel, the control principle is shown in figure 1, and a pitch and roll control system is designed by taking the multi-hull vessel as an object. The control system comprises a DSP controller, a navigation device, an execution structure T-shaped wing and a wave pressing plate. The navigation device comprises a gyroscope and a vertical accelerometer which are arranged near the center of a hull of the multi-hull ship and used for measuring the posture and the position of the multi-hull ship, wherein the pitch angle is obtained through the measurement of the gyroscope, the pitch acceleration is obtained through accurate digital differentiation, the heave speed is obtained through the integration of the accelerometer, and the heave height is obtained through the twice integration of the accelerometer. The DSP controller obtains the heave angle and the pitch angle of the multi-hull vessel through a navigation device, compares the obtained values with expected values, calculates control input according to a pre-designed control strategy and on-line compensation, and distributes the control input to the T-shaped wing and the wave pressing plate actuating mechanism; the execution structure executes the control command to swing the attack angle, control the motion attitude of the multi-hull vessel and perform stabilization. The controller adopts a control structure based on the combination of proportional-differential control and an extended state observer, the controlled observation state is a pitch angle, an angular velocity, a heave height and a heave rate, the synthesized virtual control quantity is calculated into the attack angle of a T-shaped wing and a wave pressing plate, and the change of the heave height and the pitch angle of the multi-hull vessel is controlled through the force and the moment generated by the attack angle.
The specific process of the invention is as follows:
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;respectively representing the heave velocity and the pitch angular velocity;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.
Through equivalent mathematical transformation, the coupled kinematic model of heave and pitch is converted into the following state space form through equivalent mathematical transformation:
the heave motion model is represented as:
wherein x is1Representing heave displacement, x2Represents the heave velocity, and x1=x3,The amount of pitch coupled to the heave channel being an uncertainty of the heave channel, i.e.Input F ═ FT-foil+FflapT-foil stands for T-shaped wing, flap stands for press wave plate,is the gain value.
The pitch motion model is represented as:
wherein x is11Representing pitch angle, x22Denotes pitch angular velocity, and x11=x5,The amount of heave coupled to the pitch channel motion is used as uncertainty of the pitch channel, i.e. uncertaintyInputting M ═ MT-foil+MflapT-foil stands for T-shaped wing, flap stands for press wave plate,is the gain value.
the following extended state observer is designed for the heave motion model:
parameter βiBy adopting the configuration method based on the bandwidth, the following conditions are met:
[β1,β2,β3]=[ω0α1,ω0 2α2,ω0 3α3]; (7)
wherein, βiFor adjusting the parameters, i is 1, 2, 3, ω0Selecting a gain coefficient α corresponding to the bandwidth of the heave channel extended state observeri=3!/i!×(3-i)!,i=1,2,3,e1Is the system state x1With the estimated state z of the observer1Error of (2), z1、z2Is the system state x1、x2Corresponding observer estimated state, z3Is an estimate of the lumped interference of the system, i.e. z1→x1,z2→x2,z3→x3=f2,g1i(e1)=e1,i=1、2、3;Respectively being estimated states z1、z2、z3First derivative of b1For the gain value, U, corresponding to the heave channel1The obtained virtual control quantity is designed for the heave channel.
Designing the following extended state observer aiming at the model of the pitching motion:
parameter βiiBy adopting the configuration method based on the bandwidth, the following conditions are met:
[β11,β22,β33]=[ω1α1,ω1 2α2,ω1 3α3](9)
wherein, βiiFor adjusting the parameters, i is 1, 2, 3, ω1Is the bandwidth corresponding to the pitch channel extended state observer, e2Is the system state x11With the estimated state z of the observer11Error of (2), z11、z22Is the system state x11,x22Estimate of z33Is an estimate of the lumped interference of the system, i.e. z11→x11,z22→x22,z33→x33=f22,g2i(e2)=e2,i=1、2、3; Respectively being estimated states z11、z22、z33B is the gain value corresponding to the pitch channel, U2And designing the obtained virtual control quantity for the pitching channel.
Step 4, respectively solving basic control quantities of the heave motion model and the pitch motion model by adopting a proportional-differential control method; the method comprises the following specific steps:
in a multi-hull ship heave motion model, a proportional term kp in a control law1Can be determined by the system natural frequency wnFound and the natural frequency of the multihull vessel is known, so the proportional term kp1The values of (A) are:
derivative term kd1The value is determined by the system natural frequency wnAnd damping ratio ε, so the derivative term kd1The values of (A) are:
wherein, wnLet epsilon be damping ratio, take epsilon 0.85, m be the mass of the multihull vessel, a33As an additional mass of the multi-hulled vessel, c33To be the coefficient of restitution of the system, K1The magnitude of the force generated for pressing the corrugated board.
Obtaining the proportional term kp in the pitching motion model by referring to the solution of the proportional term and the differential term in the heave motion model2And a derivative term kd2:
wherein, wnTaking epsilon as a damping ratio, I as the natural frequency of the system, and taking epsilon as 0.855Is the moment of inertia of the multihull vessel about the y-axis, a55To add moment of inertia, c55Is the coefficient of restitution of the system, K2The magnitude of the moment generated by the T-shaped wing;
and 5, combining the steps, synthesizing the basic control quantity of the multi-hull ship and the estimated coupling term to obtain a virtual control quantity, and distributing the virtual control quantity to the attack angle of the T-shaped wing and the wave pressing plate.
Step 5.1, combining the proportional term and the differential term of the heave motion model and the pitch motion model to respectively obtain the basic control input β of the heave motion model and the pitch motion model1、β2;
The basic control quantities of the heave channel are:
the basic control quantities for the pitch channel are:
wherein x is1For heave displacement, x2For the heave velocity, a33、a55For additional mass and additional moment of inertia of the multihull vessel, b33、b55Is the damping coefficient of the system, c33、c55The coefficient of the restoring force of the system, m is the mass of the multi-hull ship, and the heave motion model proportion term kp1And a derivative term kd1,I55Is the moment of inertia, x, of the multihull vessel about the y-axis11Pitch angle, x22Is the pitch angular velocity; scaling term kp of multi-hull ship pitching channel2And a derivative term kd2。
Step 5.2, two basic control input comprehensive disturbance lumped interference estimated values z obtained based on the heave and pitch motion models3And z33To determine the final virtual control quantity U of the heave motion model and the pitch motion model1、U2;
U1=β1-z3/b1; (16)
U2=β2-z33/b2; (17)
Wherein, β1As a basic control input for the heave channel, β2As a basic control input for the pitch channel, z3(t)、z33(t) is the lumped disturbance estimate for the heave and pitch channels,
step 5.3, according to the virtual control quantity U of the heave motion model and the pitch motion model1、U2And obtaining the comprehensive control quantity of the multi-hull ship, namely the attack angle input of the T-shaped wing and the wave pressing plate.
Wherein, α1Angle of attack for pressed plates, α2Angle 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.
The calculation simulation and the water pool experiment show that the sea condition grade of the multihull ship with the navigational speed of 14kn, the course angle of 180 degrees and the ship body water line length of about 48.4m is 4 grades. According to the images of the relationship between heave displacement and time and the images of heave speed and time shown in the attached figures 2 and 3, it is obvious that the heave of the multi-hulled vessel is reduced by 20-35% compared with the multi-hulled vessel without a controller; with reference to fig. 4 and 5, the pitch angle and pitch angular velocity versus time images, it is evident that with the controller, the pitch of the multihull vessel is reduced by 40% -50% compared to without the controller. The stability of the system can be obviously improved, and the rolling reduction performance of the multi-hull ship is improved. Besides, the observer can well estimate the state of the system. The invention is equally applicable to other multihull vessel controls.
The above embodiments are only used 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 accordingly, and the protection scope of the present invention is not limited to the above embodiments. Therefore, all equivalent changes and modifications made in accordance with the principles and concepts disclosed herein are intended to be included within the scope of the present invention.
Claims (8)
1. A simplified robust adaptive pitch reduction control method for a multihull vessel is characterized by comprising the following steps:
step 1, establishing a coupled motion model of heave and pitch of a multi-hull ship, and converting the model into a state space form;
step 2, decomposing the coupled motion model of heave and pitch of the multi-hull vessel into a decoupled heave motion model and a decoupled pitch motion model;
step 3, aiming at the decoupled heave motion model and pitch motion model, respectively designing an extended state observer for estimating the motion state and the coupling term of the multi-hull vessel; in the step 3, the following extended state observer is designed for the heave motion model:
parameter βiBy adopting the configuration method based on the bandwidth, the following conditions are met:
[β1,β2,β3]=[ω0α1,ω0 2α2,ω0 3α3];
wherein, βiFor adjusting the parameters, i is 1, 2, 3, ω0Selecting a gain coefficient α corresponding to the bandwidth of the heave channel extended state observeri=3!/i!×(3-i)!,i=1,2,3,x1Representing heave displacement, x2Represents the heave velocity, e1Is x1With the estimated state z of the observer1Error of (2), z1、z2Is x1、x2Corresponding observer estimated state, z3Is an estimate of the lumped interference of the system, i.e. z1→x1,z2→x2,z3→x3=f2,f2The amount of motion coupled to the heave channel for pitch and roll is used as an uncertainty for the heave channel, m being the mass of the multihull vessel; a is33Is an additional mass of the multi-hulled vessel; a is35、b35、c35Coupling term coefficient of force and moment; x is the number of11Represents the pitch angle; x is the number of22Representing pitch angular velocity, FwaveRepresenting the disturbance force of sea waves; error g1i(e1)=e1,i=1、2、3;Respectively being estimated states z1、z2、z3First derivative of b1For the gain value, U, corresponding to the heave channel1Designing the virtual control quantity for the heave channel;
designing the following extended state observer aiming at the model of the pitching motion:
parameter βiiBy adopting the configuration method based on the bandwidth, the following conditions are met:
[β11,β22,β33]=[ω1α1,ω1 2α2,ω1 3α3];
wherein, βiiFor adjusting the parameters, i is 1, 2, 3, ω1Is the bandwidth, x, corresponding to the pitch channel extended state observer11Representing pitch angle, x22Representing pitch angular velocity, e2Is x11With the estimated state z of the observer11Error of (2), z11、z22Is x11、x22Estimate of z33Is an estimate of the lumped interference of the system, i.e. z11→x11,z22→x22,z33→x33=f22,f22An amount of motion coupled to the pitch channel for heave as an uncertainty for the pitch channel; i is55Is the moment of inertia of the multihull vessel about the y-axis, a55Additional mass and additional moment of inertia for the multihull vessel; a is53、b53、c53Coupling term coefficient of force and moment; mwaveRepresenting the wave moment; error g2i(e2)=e2,i=1、2、3; Respectively being estimated states z11、z22、z33B is the gain value corresponding to the pitch channel, U2Designing the obtained virtual control quantity for the pitching channel;
step 4, respectively solving basic control quantities of the heave motion model and the pitch motion model by adopting a proportional-differential control method;
and 5, combining the steps, synthesizing the basic control quantity of the multi-hull ship and the estimated coupling term to obtain a virtual control quantity, and distributing the virtual control quantity to the attack angle of the T-shaped wing and the wave pressing plate.
2. The simplified robust adaptive pitch reduction control method for a multihull vessel according to claim 1, wherein the coupled motion model of the multihull vessel heave and pitch in step 1 is:
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;respectively representing the heave velocity and the pitch angular velocity;respectively representing heave acceleration and pitch angular acceleration;respectively 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.
4. the simplified robust adaptive pitch reduction control method for a multihull vessel according to claim 2, wherein in step 2, the heave motion model is expressed as:
wherein x is1=x3,x1Which is indicative of the displacement of the heave,x2representing heave velocity, the amount of pitch coupled to the heave channel as an uncertainty, i.e. uncertainty, of the heave channelInput force F ═ FT-foil+Fflap,Is the gain value;
the pitch motion model is represented as:
5. The simplified robust adaptive pitch reduction control method for a multihull vessel according to claim 1, wherein the basic control quantities in step 4 are:
heave motion model scale term kp1And a derivative term kd1Comprises the following steps:
wherein, wnLet epsilon be damping ratio, take epsilon 0.85, m be the mass of the multihull vessel, a33As an additional mass of the multi-hulled vessel, b33Is the damping coefficient of the system, c33To be the coefficient of restitution of the system, K1The magnitude of the force generated by the press wave plate;
scaling term kp for pitching channel2And a derivative term kd2The method comprises the following steps:
wherein, wnTaking epsilon as a damping ratio, I as the natural frequency of the system, and taking epsilon as 0.855Is the moment of inertia of the multihull vessel about the y-axis, a55To add moment of inertia, b55Is the damping coefficient of the system, c55Is the coefficient of restitution of the system, K2The magnitude of the moment generated for the T-shaped wing.
6. The simplified robust adaptive pitch reduction control method for a multihull vessel according to claim 1, wherein the procedure of step 5 is:
step 5.1, combining the heave motion model and the pitch motionThe proportional term and the differential term of the motion model respectively obtain basic control quantities β of the heave motion model and the pitch motion model1、β2;
Step 5.2, two basic control quantity comprehensive lumped interference estimated values z obtained based on the heave and pitch motion models3And z33To determine the final virtual control quantity U of the heave motion model and the pitch motion model1、U2;
Step 5.3, according to the virtual control quantity U of the heave motion model and the pitch motion model1、U2And obtaining the comprehensive control quantity of the multi-hull ship, namely the attack angle input of the T-shaped wing and the wave pressing plate.
7. The simplified robust adaptive pitch reduction control method for a multihull vessel according to claim 6, wherein the basic control quantities β of the heave motion model and the pitch motion model in step 5.1 are1、β2Expressed as: the basic control quantities of the heave channel are:
the basic control quantities for the pitch channel are:
wherein x is1For heave displacement, x2For the heave velocity, a33、a55For additional mass and additional moment of inertia of the multihull vessel, b33、b55Is the damping coefficient of the system, c33、c55The coefficient of the restoring force of the system, m is the mass of the multi-hull ship, and the heave motion model proportion term kp1And a derivative term kd1,I55Is the moment of inertia, x, of the multihull vessel about the y-axis11Pitch angle, x22Is the pitch angular velocity; scaling term kp of multi-hull ship pitching channel2And a derivative term kd2。
8. The simplified robust adaptive pitch reduction control method for a multihull vessel according to claim 6, wherein the virtual control quantity is:
U1=β1-z3/b1;
U2=β2-z33/b2;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811100020.0A CN109334892B (en) | 2018-09-20 | 2018-09-20 | Simplified robust self-adaptive pitching reduction control method for multi-hull vessel |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811100020.0A CN109334892B (en) | 2018-09-20 | 2018-09-20 | Simplified robust self-adaptive pitching reduction control method for multi-hull vessel |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109334892A CN109334892A (en) | 2019-02-15 |
CN109334892B true CN109334892B (en) | 2020-06-26 |
Family
ID=65305836
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811100020.0A Active CN109334892B (en) | 2018-09-20 | 2018-09-20 | Simplified robust self-adaptive pitching reduction control method for multi-hull vessel |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109334892B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111332424B (en) * | 2020-03-24 | 2021-02-02 | 上海大学 | Water surface robot full-navigational-speed stabilization and stabilization reduction method in lift-drag combined mode |
CN111498037B (en) * | 2020-04-05 | 2022-05-20 | 哈尔滨工程大学 | High-speed catamaran longitudinal stabilization method based on variable structure active disturbance rejection control |
CN111552186B (en) * | 2020-06-22 | 2022-06-07 | 大连海事大学 | Rudder and fin combined anti-rolling control method and system based on CAN bus |
CN112068437B (en) * | 2020-09-17 | 2022-06-21 | 哈尔滨工程大学 | Single-step prediction control anti-rolling method for high-speed multi-hull ship |
CN114802594B (en) * | 2022-06-02 | 2022-10-18 | 哈尔滨理工大学 | Three-degree-of-freedom shipborne stable platform with heave compensation |
SE2350689A1 (en) * | 2023-06-07 | 2024-06-18 | Volvo Penta Corp | Flight control systems for hydrofoiling marine vessels |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5511504A (en) * | 1995-08-09 | 1996-04-30 | Martin; John R. | Computer controlled fins for improving seakeeping in marine vessels |
CN101825871A (en) * | 2010-04-09 | 2010-09-08 | 哈尔滨工程大学 | Intelligent adaptive control method for heave and pitch device for oblique rudder ship |
CN103935482A (en) * | 2014-05-04 | 2014-07-23 | 中国舰船研究设计中心 | Automatic navigation control system improving seakeeping performance of ship |
CN103955226A (en) * | 2014-04-30 | 2014-07-30 | 中国舰船研究设计中心 | Pitching and heaving proportion-differential control method for wave-piercing catamaran |
CN103963938A (en) * | 2014-05-14 | 2014-08-06 | 哈尔滨工程大学 | Intelligent self-adaptive prediction control system based on slant-rudder ship anti-pitching device |
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 |
CN107255919A (en) * | 2017-06-13 | 2017-10-17 | 哈尔滨工程大学 | A kind of trimaran pitching stabilization decoupling controller and design method |
-
2018
- 2018-09-20 CN CN201811100020.0A patent/CN109334892B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5511504A (en) * | 1995-08-09 | 1996-04-30 | Martin; John R. | Computer controlled fins for improving seakeeping in marine vessels |
CN101825871A (en) * | 2010-04-09 | 2010-09-08 | 哈尔滨工程大学 | Intelligent adaptive control method for heave and pitch device for oblique rudder ship |
CN103955226A (en) * | 2014-04-30 | 2014-07-30 | 中国舰船研究设计中心 | Pitching and heaving proportion-differential control method for wave-piercing catamaran |
CN103935482A (en) * | 2014-05-04 | 2014-07-23 | 中国舰船研究设计中心 | Automatic navigation control system improving seakeeping performance of ship |
CN103963938A (en) * | 2014-05-14 | 2014-08-06 | 哈尔滨工程大学 | Intelligent self-adaptive prediction control system based on slant-rudder ship anti-pitching device |
CN107255919A (en) * | 2017-06-13 | 2017-10-17 | 哈尔滨工程大学 | A kind of trimaran pitching stabilization decoupling controller and design method |
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 |
Also Published As
Publication number | Publication date |
---|---|
CN109334892A (en) | 2019-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109334892B (en) | Simplified robust self-adaptive pitching reduction control method for multi-hull vessel | |
CN109814392B (en) | Fault-tolerant control method for fault robustness of under-actuated underwater robot actuator | |
CN106292287A (en) | A kind of UUV path following method based on adaptive sliding-mode observer | |
CN111580523A (en) | Unmanned ship path tracking active disturbance rejection control method based on sideslip angle compensation | |
CN112036099B (en) | Method for determining minimum turning radius based on attack angle of sail of sailing boat | |
CN113359785B (en) | Microminiature AUV underwater motion and hovering control method | |
CN109739248B (en) | Ship-borne three-degree-of-freedom parallel stable platform stability control method based on ADRC | |
CN103935482A (en) | Automatic navigation control system improving seakeeping performance of ship | |
CN109240289B (en) | Wave glider bow information adaptive filtering method | |
Wille et al. | Modeling and course control of sailboats | |
CN114611206A (en) | Environment-optimal heading positioning method of dynamic positioning ship | |
Matusiak | Dynamics of a Rigid Ship-with applications | |
CN112434428B (en) | Ship transverse-swinging nonlinear dynamics analysis method in regular wave | |
CN102910264B (en) | A kind of antirolling apparatus of ship mooring state | |
Tang et al. | Simulation of optimal integral sliding mode controller for the depth control of AUV | |
CN117250971A (en) | Control method of microminiature AUV | |
CN103955226A (en) | Pitching and heaving proportion-differential control method for wave-piercing catamaran | |
CN115616920B (en) | Active coupling interference suppression method in multi-unmanned-ship system | |
CN115783159A (en) | Ship power anti-rolling device and anti-rolling control method | |
Li et al. | Dynamics modeling and experiments of wave driven robot | |
CN116048090A (en) | Sail navigation aid ship path tracking control method with energy consumption optimization effect | |
CN113848915A (en) | Polar ship course control method based on linear active disturbance rejection technology | |
Liu | Adaptive sliding mode control for ship autopilot with speed keeping | |
CN111498037B (en) | High-speed catamaran longitudinal stabilization method based on variable structure active disturbance rejection control | |
CN115686034A (en) | Unmanned underwater vehicle track tracking control method considering failure of speed sensor |
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 |