CN113110023B - Structure correction control method based on Hamiltonian model of diesel engine - Google Patents

Structure correction control method based on Hamiltonian model of diesel engine Download PDF

Info

Publication number
CN113110023B
CN113110023B CN202110360225.8A CN202110360225A CN113110023B CN 113110023 B CN113110023 B CN 113110023B CN 202110360225 A CN202110360225 A CN 202110360225A CN 113110023 B CN113110023 B CN 113110023B
Authority
CN
China
Prior art keywords
diesel engine
matrix
hamiltonian
actuator
control
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110360225.8A
Other languages
Chinese (zh)
Other versions
CN113110023A (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.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
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 Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN202110360225.8A priority Critical patent/CN113110023B/en
Publication of CN113110023A publication Critical patent/CN113110023A/en
Application granted granted Critical
Publication of CN113110023B publication Critical patent/CN113110023B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/36Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)

Abstract

The invention discloses a structural correction control method based on a Hamiltonian model of a diesel engine, which belongs to the technical field of control, wherein the output control of the diesel engine is a key problem in the operation of the diesel engine.

Description

Structure correction control method based on Hamiltonian model of diesel engine
Technical Field
The invention belongs to the technical field of diesel engine control, and particularly relates to a structural correction control method based on a Hamiltonian model of a diesel engine.
Background
The diesel engine belongs to mechanical equipment for providing power moment, and the control research on the output characteristics of the diesel engine is an important point in the research of the diesel engine. At present, control constructed based on a diesel engine model can be roughly divided into two types, namely, PID parameter optimization based on traditional PID control and combined with a new control theory, and new controller replacement PID control is constructed, and the design ideas of the two types of control are different, but the aim is to realize the inhibition of overshoot in the transient process and the reduction of transient oscillation time. The generalized Hamiltonian system has outstanding characteristics in solving the control problem of a nonlinear system, and the structure and the damping matrix of the generalized Hamiltonian system can display the dynamic information related between the internal parameters of the model, so that the research on the control characteristics of the diesel engine based on the Hamiltonian model is a new thought.
Disclosure of Invention
The invention adopts the hamilton structure correction method to design the additional control of the diesel engine, designs the control on the basis of the hamilton model of the diesel engine through the hamilton structure correction method,the diesel engine control law is obtained through structural matrix factor design under the condition of ensuring the stability of a diesel engine system, feedback stabilization control is realized, and the diesel engine is taken as a research object by the control method through simulation verification, so that the purpose of finding a corrected control variable is achieved; firstly, directly giving a Hamiltonian model of a diesel engine, and selecting a modified Hamiltonian function Hd and a modification matrix J on the basis of the Hamiltonian model α And R is α Calculating a vector function K (x), and deducing a control variable u according to the Hamiltonian correction theory α And M
The invention is based on the Hamiltonian of the diesel engine:
wherein: h (x) is Hamiltonian, m 1 (kg) is the mass of the moving part of the diesel engine actuator, x 1 (m/s) is the displacement of the actuator output; x is x 2 (m) is the actuator movement speed; k (k) 1 Is the elastic coefficient of the spring in the actuator; t (T) j Is an inertial time constant; omega 1 Is the electrical angular velocity increment (per-sign value) of the rotation axis of the diesel engine; omega B =314 (rad/s) is a reference value of the electrical angular velocity;
a non-linear system port controlled hamiltonian system representation can be written as:
x is the state variable matrix, i.e. x= [ x ] 1 x 2 ω 1 ] T W (x) is a control variable, J (x) is an antisymmetric matrix, i.e., J (x) = -J T (x) R (x) is a damping symmetry matrix, i.e. R (x) =r T (x) More than or equal to 0, g (x) is an input matrix, and y is an output matrix;
wherein:
wherein: k (k) x1 、k u 、a 1 、d 1 、c 1 、k ω Is a diesel engine actuator parameter; d (D) d The equivalent damping coefficient of the diesel engine; w (x) is a control variable; m is M L Controlling the moment for the load; t (T) j Is an inertial time constant(s); u is the control amount of the diesel engine.
The design of the invention is based on the following:
theorem: given a balance point x 0 And desired structural modification J a (x) Damping correction R a (x) It is assumed that a control α (x) and a vector function K (x) can be found that satisfy:
and causing:
(i)J d (x)=J(x)+J α (x)=-[J(x)+J α (x)] T ;R d (x)=R(x)+R α (x)=[R(x)+R α (x)] T
(ii)
(iii) At the equilibrium point x 0 The method comprises the following steps:
(iv) At the equilibrium point x 0 The method comprises the following steps:
(iv)
the closed loop system at u (x) =α (x) can be:
wherein H is d (x)=H(x)+H α (x),x 0 Is the (local) stable equilibrium point of the closed loop system, H a (x) Is an energy correction function; h d (x) Is a quasi-Lyapunov function of the equilibrium point.
Based on the conditions, the invention provides a structure correction control method based on a Hamiltonian model of a diesel engine, wherein the diesel engine comprises an electromagnetic actuator and a body; on the basis of the Hamiltonian model of the diesel engine, a structural matrix J which is expected to be corrected is given α Finding the control α (x) and the vector function K (x) allows the system to obtain the additional control variable u under stable conditions α And M The method comprises the following specific steps:
the first step: selecting Hamiltonian and correction matrix of corrected system
Selecting Hamiltonian H of the corrected system d The following are provided:
wherein x is 0 =[x 10 x 20 ω 1(0) ] T Is the initial value of the state variable;
and a second step of: design structure correction matrix J α (x)
Since there are three variables in the Hamiltonian model, the structural matrix J is sought α (x) For a 3×3 matrix, only the correction of the principal element is considered in the correction design, due to the actuator displacement x 1 And actuator speed x 2 Rotational speed omega of shaft 1 Is the influence of two main factors, so in the antisymmetric matrix J a (x) Setting correction element J at corresponding position of (2) 13 、J 23 Wherein J 13 Representing the effect of displacement on output shaft speed, J 23 The other positions are set to zero, representing the effect of actuator speed on shaft speed. Due to not matching R α The matrix is modified and may be given as a constant in design. J (J) α And R is α The following are provided:
and a third step of: calculating vector matrix function K (x)
Vector function K (x):
according to the relation that K should satisfy in theorem (iii), the following is calculated:
let the feedback control α be:
substituting (1) (2-1) (2-2) (4) (5) (6) (7) into (3), two expressions can be obtained as follows:
from (8) and (9) the control variable u can be derived α And a load moment M 2 The expression of (2) is as follows:
wherein:representing control before correction;
representing additional control resulting from the structural modification;
representing the load control torque before correction;
representing additional load control torque generated by the structural modification;
u α and M Is a modified control variable, and after such design, formulas (10) and (11) contain J 13 And J 23 That is, the structure correction factor J is realized 13 And J 23 Is involved in the control of the diesel engine.
The beneficial effects of the invention are as follows:
the invention discloses a control method for structural modification on the basis of a Hamiltonian model of a diesel engine, which does not need to change the physical structure of the diesel engine, and achieves the purpose of modifying output by modifying structural parameters of a model of the diesel engine. The method is realized by means of software.
Drawings
FIG. 1 is J 23 Influence on diesel engine output;
FIG. 2 is J 13 Impact on diesel engine output.
Detailed Description
The invention will be further illustrated with reference to specific examples.
Example 1
The diesel engine hamilton function according to this embodiment:
wherein: h (x) is Hamiltonian, m 1 (kg) is the mass of the moving part of the diesel engine actuator, x 1 (m) is the displacement of the actuator output; x is x 2 (m/s) is the actuator movement speed; k (k) 1 Is the elastic coefficient of the spring in the actuator; t (T) j Is an inertial time constant; omega 1 Is the electrical angular velocity increment (per-sign value) of the rotation axis of the diesel engine; omega B =314 (rad/s) is a reference value of the electrical angular velocity;
a non-linear system port controlled hamiltonian system representation can be written as:
x is the state variable matrix, i.e. x= [ x ] 1 x 2 ω 1 ] T W (x) is a control variable, J (x) is an antisymmetric matrix, i.e., J (x) = -J T (x) R (x) is a damping symmetry matrix, i.e. R (x) =r T (x) More than or equal to 0, g (x) is an input matrix, and y is an output matrix;
wherein:
wherein: k (k) x1 、k u 、a 1 、d 1 、c 1 、k ω Is a diesel engine actuator parameter; d (D) d The equivalent damping coefficient of the diesel engine; w (x) is a control variable; m is M L Controlling the moment for the load; t (T) j Is an inertial time constant(s); u is the control amount of the diesel engine.
In the above prior art, the design of this embodiment is based on the following:
design theorem: given a balance point x 0 And desired structural modification J a (x) Damping correction R a (x) Assuming that a control α (x) and a vector function K (x) can be found, the following needs to be satisfied:
and causing:
(i)J d (x)=J(x)+J α (x)=-[J(x)+J α (x)] T ;R d (x)=R(x)+R α (x)=[R(x)+R α (x)] T
(ii)
(iii) At the equilibrium point x 0 The method comprises the following steps:
(iv) At the equilibrium point x 0 The method comprises the following steps:
(iv)
the closed loop system at u (x) =α (x) can be:
wherein H is d (x)=H(x)+H α (x),
x 0 Is the (local) stable equilibrium point of the closed loop system, H a (x) Is an energy correction function; h d (x) Is a quasi-Lyapunov function of the equilibrium point.
Based on the above conditions, the embodiment provides a structure correction control method based on a Hamiltonian model of a diesel engine, wherein the diesel engine comprises an electromagnetic actuator and a body; on the basis of the Hamiltonian model of the diesel engine, a structural matrix J alpha which is expected to be corrected is given, and the control alpha (x) and the vector function K (x) are searched, so that the system obtains an additional control variable u under the stable condition α And M The method comprises the following specific steps:
the first step: selecting Hamiltonian and correction matrix of closed-loop system
Selecting Hamiltonian H of the corrected system d The following are provided:
wherein x is 0 =[x 10 x 20 ω 1(0) ] T Is the initial value of the state variable;
and a second step of: design structure correction matrix J α (x)
Since there are three variables in the Hamiltonian model, the structural matrix J is sought α (x) For a 3×3 matrix, only the correction of the principal element is considered in the correction design, due to the actuator displacement x 1 And actuator speed x 2 Rotational speed omega of shaft 1 Is the influence of two main factors, so in the antisymmetric matrix J a (x) Setting correction element J at corresponding position of (2) 13 、J 23 Wherein J 13 Representing the effect of displacement on output shaft speed, J 23 Representing the effect of actuator speed on shaft speed, other positions are set to zero, since R is not α The matrix is corrected, and can be given as a constant in design, J α And R is α The following are provided:
and a third step of: calculating vector matrix function K (x)
Vector function K (x):
according to the relation that K should satisfy in theorem (iii), the following is calculated:
let the feedback control α be:
substituting (1) (2-1) (2-2) (4) (5) (6) (7) into (3), two expressions can be obtained as follows:
from (8) and (9) the control variable u can be derived α And a load moment M The expression of (2) is as follows:
wherein:representing control before correction;
representing additional control resulting from the structural modification;
representing the load control torque before correction;
representing additional load control torque generated by the structural modification; u (u) α And M Is a modified control variable, and after such design, formulas (10) and (11) contain J 13 And J 23 That is, the structure correction factor J is realized 13 And J 23 Is involved in the control of the diesel engine.
Example 2
In this embodiment, the conclusion of embodiment 1 is simulated by using the data of a certain diesel engine, and the specific process is as follows:
the load carried by the diesel engine is selected as a generator under the simulation working condition, and the diesel generator is startedThe initial operating conditions are as follows: the active power output by the generator is p e =0.8 (normalized value), the generator power factor is 0.8, a step disturbance of 0.1 (normalized value) occurs at 0.1s generator load, and 0.2s disturbance disappears.
Simulation parameters: the rated power of the diesel engine is 1250kW omega B Rated rotation speed n=1500r/min, mass m 1 =0.8 kg, mechanical damping coefficient a 1 =-1245,c 1 =10.0N·s/mm,d 1 =15042,k c =0.2,k d =5.54,k ω =20,k x1 =200,k u =40, spring rate k 1 =3610; output shaft maximum stroke 10mm, moment of inertia j=71.822 kg.m 2 The method comprises the steps of carrying out a first treatment on the surface of the The number of pole pairs of the generator p=2; equivalent damping coefficient D of diesel engine d =2.1753。
The Hamiltonian models of the diesel engine in example 1, namely (2), (2-1), (2-2), (2-3), (2-4) and the deduced control expressions (10) and (11), are combined to obtain a diesel engine simulation model with correction control, and the correction factor J 13 And J 23 The results are shown in fig. 1 and 2 for the diesel engine.
As can be seen from the simulation results, in FIG. 1, change J 23 Different oscillation curves can be obtained, which shows that J23 plays different roles in inhibiting oscillation; in FIG. 2, change J 13 The oscillation of the load mutation can be effectively restrained, and the restraining effect on the active oscillation can be obtained by selecting proper structural factors.

Claims (1)

1. The structure correction control method based on the Hamiltonian model of the diesel engine is characterized in that the diesel engine comprises an electromagnetic actuator and a body; on the basis of the Hamiltonian model of the diesel engine, a structural matrix J which is expected to be corrected is given α Finding the control α (x) and the vector function K (x) allows the system to obtain the additional control variable u under stable conditions α And M
The method comprises the following specific steps:
the first step: selecting Hamiltonian H of closed loop system d
Selecting and correctingHamiltonian function H of (1) d The following are provided:
wherein: x is a state variable; x is x 1 The displacement of the diesel engine actuator is m; x is x 2 For the movement speed of the actuator, m/s; omega 1 The value of the angular velocity increment of the main shaft of the diesel engine is shown as a per-mark value; x= [ x ] 1 x 2 ω 1 ] T ;x 0 Is the initial value of the state variable; x is x 0 =[x 10 x 20 ω 1(0) ] T k 1 N/m is the elastic coefficient of the actuator spring; m is m 1 Kg, the mass of the moving part of the actuator; t (T) j Is an inertial time constant, s; omega B =314 rad/s, which is a reference value of angular velocity;
and a second step of: design structure correction matrix J α (x)
Since the number of variables in the diesel engine model is three, the structural matrix J is sought α (x) For a 3 x 3 matrix, only the correction of the principal element is considered in the correction design, due to the actuator displacement x 1 And actuator speed x 2 Rotational speed omega of shaft 1 Is the influence of two main factors, so in the antisymmetric matrix J a (x) Setting correction element J at corresponding position of (2) 13 、J 23 Wherein J 13 Representing the effect of displacement on output shaft speed, J 23 Representing the effect of actuator speed on shaft speed, other positions are set to zero for R α The matrix can be given as constant, J, since it is not modified α And R is α The form is as follows:
and a third step of: calculating vector matrix function K (x)
According to H d And H (x) calculation:wherein H is a (x) Is an energy correction function, H α =H d -H (x); h (x) is a Hamiltonian of the diesel engine and is obtained from a Hamiltonian model of the diesel engine;
the invention is based on the Hamiltonian of the diesel engine:
wherein: h (x) is the diesel hamiltonian;
the hamilton model is:
w (x) is a control variable, J (x) is an antisymmetric matrix, i.e., J (x) = -J T (x) R (x) is a damping symmetry matrix, i.e. R (x) =r T (x) More than or equal to 0, g (x) is an input matrix, and y is an output matrix;
wherein:
wherein: k (k) x1 、k u 、a 1 、d 1 、c 1 、k ω Is a diesel engine actuator parameter; d (D) d The equivalent damping coefficient of the diesel engine; w (x) is a control variable; u is the control quantity of the diesel engine; m is M L Controlling the moment for the load;
let the feedback control be alpha-let,hamiltonian H (x), J (x), R (x) and H in Hamiltonian model of diesel engine d 、K(x)、J α 、R α Alpha is substituted into the following formula:
obtaining the control variable u a And a load moment M The expression of (2) is as follows:
CN202110360225.8A 2021-04-02 2021-04-02 Structure correction control method based on Hamiltonian model of diesel engine Active CN113110023B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110360225.8A CN113110023B (en) 2021-04-02 2021-04-02 Structure correction control method based on Hamiltonian model of diesel engine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110360225.8A CN113110023B (en) 2021-04-02 2021-04-02 Structure correction control method based on Hamiltonian model of diesel engine

Publications (2)

Publication Number Publication Date
CN113110023A CN113110023A (en) 2021-07-13
CN113110023B true CN113110023B (en) 2024-02-27

Family

ID=76713760

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110360225.8A Active CN113110023B (en) 2021-04-02 2021-04-02 Structure correction control method based on Hamiltonian model of diesel engine

Country Status (1)

Country Link
CN (1) CN113110023B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115929534B (en) * 2022-12-12 2023-09-05 昆明理工大学 Method for inhibiting pressure pulsation induced active oscillation based on Hamiltonian model

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101813053A (en) * 2010-03-10 2010-08-25 昆明理工大学 Correcting and controlling method of interconnecting structure for improving power-angle oscillation of hydroelectric generating set
CN103066902A (en) * 2013-01-24 2013-04-24 河海大学 Method for realizing passive control law of direct-current motor based on load observation
JP2015095126A (en) * 2013-11-12 2015-05-18 富士通株式会社 Information processing method, device, and program
CN108365788A (en) * 2018-04-17 2018-08-03 哈尔滨理工大学 A kind of Matrix Converter-Permanent Magnetic Synchronous Machine governing system and method based on passive coherent locating
CN111668867A (en) * 2019-03-05 2020-09-15 南京理工大学 Passive sliding mode control method for wind power plant through VSC-HVDC system grid connection
CN112464391A (en) * 2020-10-15 2021-03-09 昆明理工大学 Method for establishing Hamilton model of diesel generator set

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10197040B2 (en) * 2016-04-14 2019-02-05 National Technology & Engineering Solutions Of Sandia, Llc Optimal control of wave energy converters
US10538174B2 (en) * 2017-09-29 2020-01-21 Nio Usa, Inc. Real-time nonlinear receding horizon control of batteries for power systems

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101813053A (en) * 2010-03-10 2010-08-25 昆明理工大学 Correcting and controlling method of interconnecting structure for improving power-angle oscillation of hydroelectric generating set
CN103066902A (en) * 2013-01-24 2013-04-24 河海大学 Method for realizing passive control law of direct-current motor based on load observation
JP2015095126A (en) * 2013-11-12 2015-05-18 富士通株式会社 Information processing method, device, and program
CN108365788A (en) * 2018-04-17 2018-08-03 哈尔滨理工大学 A kind of Matrix Converter-Permanent Magnetic Synchronous Machine governing system and method based on passive coherent locating
CN111668867A (en) * 2019-03-05 2020-09-15 南京理工大学 Passive sliding mode control method for wind power plant through VSC-HVDC system grid connection
CN112464391A (en) * 2020-10-15 2021-03-09 昆明理工大学 Method for establishing Hamilton model of diesel generator set

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
哈密顿结构修正的控制设计方法及其应用;曾云等;电机与控制学报;第18卷(第3期);正文的第93-100页 *

Also Published As

Publication number Publication date
CN113110023A (en) 2021-07-13

Similar Documents

Publication Publication Date Title
Mao et al. Design and implementation of continuous finite-time sliding mode control for 2-DOF inertially stabilized platform subject to multiple disturbances
CN103701368B (en) The energy-conservation anti-backlash control method of bi-motor
CN107121932B (en) Motor servo system error symbol integral robust self-adaptive control method
Ma et al. Polynomial-method-based design of low-order controllers for two-mass systems
CN109581868B (en) Reconfigurable robot decentralized nerve optimal control method based on judgment and identification structure
CN110707981A (en) Permanent magnet synchronous motor speed controller based on novel extended state observer
CN108599649A (en) PMSM positional servosystem High order Plant controller designs and parameter determination method
CN114839880B (en) Self-adaptive control method based on flexible joint mechanical arm
Cho et al. An indirect model reference adaptive fuzzy control for SISO Takagi–Sugeno model
CN113110023B (en) Structure correction control method based on Hamiltonian model of diesel engine
CN111679576B (en) Variable cycle engine controller design method based on improved deterministic strategy gradient algorithm
CN102629293A (en) Simulation calculation method for additional mechanical damping of hydropower set
CN116317794A (en) High-precision control method for electric actuator of aero-engine
CN110842911B (en) Flexible mechanical arm combined modeling and sliding mode control method considering joint motor characteristics
CN113283090B (en) Friction compensation method for electric servo system of airplane steering engine based on bee colony algorithm
CN103117692B (en) Control method of mechanical elastic energy storing permanent magnet motor group under various external disturbances
CN107092189A (en) Multivariable based on Model Predictive Control inputs the control method of EHA systems
Kumar et al. An adaptive particle swarm optimization algorithm for robust trajectory tracking of a class of under actuated system
CN117283562A (en) Flexible mechanical arm control method based on preset time stabilization and high-order sliding mode algorithm
CN109884881A (en) A kind of design for surely taking aim at servo controller based on nonlinear PID controller technology
Huang et al. State feedback control for stabilization of the ball and plate system
Yonezawa et al. Model-free active vibration control approach using proof-mass actuator with uncertainty
Dulau et al. Design of robust control for single machine infinite bus system
Lai et al. Modeling and control for stability and rotation velocity of a rotary inverted pendulum
Li et al. Adaptive backstepping sliding mode control for the oscillation displacement system of continuous casting mold with mismatched disturbances

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