CN111781836B - An adaptive asymptotic control method for hydraulic pressure preset performance - Google Patents
An adaptive asymptotic control method for hydraulic pressure preset performance Download PDFInfo
- Publication number
- CN111781836B CN111781836B CN202010710755.6A CN202010710755A CN111781836B CN 111781836 B CN111781836 B CN 111781836B CN 202010710755 A CN202010710755 A CN 202010710755A CN 111781836 B CN111781836 B CN 111781836B
- Authority
- CN
- China
- Prior art keywords
- formula
- adaptive
- control
- error
- bounded
- 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
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 14
- 238000013461 design Methods 0.000 claims abstract description 17
- 230000001052 transient effect Effects 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 239000003921 oil Substances 0.000 claims description 12
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 239000010727 cylinder oil Substances 0.000 claims description 3
- 239000010720 hydraulic oil Substances 0.000 claims description 3
- 230000003313 weakening effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000003723 Smelting Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004070 electrodeposition Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004154 testing of material Methods 0.000 description 1
Images
Classifications
-
- 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)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了液压控制系统技术领域的一种液压力预设性能自适应渐近控制方法,包括如下步骤:建立双出杆液压缸伺服系统模型;设计液压力预设性能自适应渐近控制器;调节上述步骤中的参数,使系统满足控制性能指标,通过鲁棒误差符号积分函数对系统干扰进行抑制,同时使用自适应控制器对系统参数进行逼近,大大削弱不确定性对系统的影响;此外,考虑到系统瞬态控制精度问题,通过预设性能函数的设计对系统的控制精度进行约束,理论上实现了系统的渐近跟踪,保证了双出杆液压缸伺服系统的输出力能准确地跟踪期望的力指令;本发明简化了控制器设计,更利于在工程实际中应用。
The invention discloses a hydraulic pressure preset performance adaptive asymptotic control method in the technical field of hydraulic control systems, comprising the following steps: establishing a double-rod hydraulic cylinder servo system model; designing a hydraulic pressure preset performance adaptive asymptotic controller ; Adjust the parameters in the above steps to make the system meet the control performance index, suppress the system disturbance through the robust error sign integral function, and at the same time use the adaptive controller to approximate the system parameters, greatly weakening the influence of uncertainty on the system; In addition, considering the transient control accuracy of the system, the control accuracy of the system is constrained by the design of the preset performance function, the asymptotic tracking of the system is theoretically realized, and the output force of the double-rod hydraulic cylinder servo system is guaranteed to be accurate. The desired force command can be tracked to the ground; the invention simplifies the design of the controller and is more conducive to the application in engineering practice.
Description
技术领域technical field
本发明涉及液压控制系统技术领域,具体为一种液压力预设性能自适应渐近控制方法。The invention relates to the technical field of hydraulic control systems, in particular to a hydraulic pressure preset performance adaptive asymptotic control method.
背景技术Background technique
电液伺服系统具有输出功率大、信号处理灵活等优点,易于实现各种参量的反馈,在功率和质量大的场合最为适合,其应用已遍及国防和工业的各个领域,比如飞机与船舶舵机的控制、雷达与火炮的控制、机床工作台的位置控制、板带轧机的板厚控制、电炉冶炼的电极位置控制、材料试验机及其他实验机的压力控制等等。然而,电液伺服系统中普遍存在的不确定性,增加了控制系统的设计难度,且可能会严重恶化能够取得的控制性能,从而导致低控制精度,极限环震荡,甚至不稳定性。为了提高电液系统的控制性能,设计人员对许多先进的非线性控制器进行了研究,如鲁棒自适应控制、自适应鲁棒控制(ARC)、滑模控制等等。虽然这些控制器能够保证良好的稳态性能,但并没有关注瞬态的控制性能,基于此,本发明设计了一种液压力预设性能自适应渐近控制方法,以解决上述问题。The electro-hydraulic servo system has the advantages of high output power and flexible signal processing, and is easy to realize the feedback of various parameters. It is most suitable for occasions with high power and quality. control, radar and artillery control, position control of machine tool table, plate thickness control of strip mill, electrode position control of electric furnace smelting, pressure control of material testing machine and other experimental machines, etc. However, the ubiquitous uncertainty in electro-hydraulic servo systems increases the design difficulty of the control system and may seriously deteriorate the achievable control performance, resulting in low control accuracy, limit cycle oscillations, and even instability. To improve the control performance of electro-hydraulic systems, designers have studied many advanced nonlinear controllers, such as robust adaptive control, adaptive robust control (ARC), sliding mode control, and so on. Although these controllers can ensure good steady-state performance, they do not pay attention to transient control performance. Based on this, the present invention designs an adaptive asymptotic control method for hydraulic pressure preset performance to solve the above problems.
发明内容SUMMARY OF THE INVENTION
本发明的目的在于提供一种液压力预设性能自适应渐近控制方法,以解决上述背景技术中提出的问题。The purpose of the present invention is to provide an adaptive asymptotic control method for hydraulic pressure preset performance, so as to solve the problems raised in the above background art.
为实现上述目的,本发明提供如下技术方案:一种液压力预设性能自适应渐近控制方法,包括如下步骤:In order to achieve the above purpose, the present invention provides the following technical solutions: a hydraulic pressure preset performance adaptive asymptotic control method, comprising the following steps:
S1:建立双出杆液压缸伺服系统模型;S1: Establish the servo system model of the double-rod hydraulic cylinder;
S2:设计液压力预设性能自适应渐近控制器;S2: Design hydraulic pressure preset performance adaptive asymptotic controller;
S3:调节步骤S2中的参数,使系统满足控制性能指标。S3: Adjust the parameters in step S2 so that the system meets the control performance index.
进一步的,所述步骤S1具体为:根据牛顿第二定律,双出杆液压缸惯性负载的动力学模型方程为:Further, the step S1 is specifically: according to Newton's second law, the dynamic model equation of the inertial load of the double-rod hydraulic cylinder is:
式中:y为负载位移,T表示驱动力,PL=P1-P2是负载驱动压力,P1和P2分别为液压缸两腔压力,A为活塞杆有效工作面积,b代表粘性摩擦系数,f为其他未建模干扰,包括非线性摩擦、外部干扰以及未建模动态,液压缸负载压力动态方程为:In the formula: y is the load displacement, T is the driving force, P L = P 1 -P 2 is the load driving pressure, P 1 and P 2 are the pressures of the two chambers of the hydraulic cylinder, A is the effective working area of the piston rod, and b is the viscosity Friction coefficient, f is other unmodeled disturbances, including nonlinear friction, external disturbances and unmodeled dynamics, the dynamic equation of hydraulic cylinder load pressure is:
式中:Vt为液压缸两腔总有效容积,Ct为液压缸泄露系数,QL=(Q1+Q2)/2是负载流量,Q1为液压缸进油腔供油流量,Q2为液压缸回油腔回油流量,q(t)为建模误差及未建模动态;In the formula: V t is the total effective volume of the two chambers of the hydraulic cylinder, C t is the leakage coefficient of the hydraulic cylinder, Q L = (Q 1 +Q 2 )/2 is the load flow, Q 1 is the oil supply flow of the hydraulic cylinder into the oil chamber, Q 2 is the return oil flow of the hydraulic cylinder oil return chamber, q(t) is the modeling error and unmodeled dynamics;
QL为伺服阀阀芯位移xv的函数:Q L is a function of servo valve spool displacement x v :
式中:为流量伺服阀的增益系数,Cd为伺服阀的流量系数,w为伺服阀的面积梯度;ρ为液压油的密度,Ps为供油压力,sign(xv)为:where: is the gain coefficient of the flow servo valve, C d is the flow coefficient of the servo valve, w is the area gradient of the servo valve; ρ is the density of the hydraulic oil, P s is the oil supply pressure, and sign(x v ) is:
假设伺服阀阀芯位移正比于控制输入u,即,xv=kiu,其中ki>0是比例系数,u是控制输入电压,因此,等式(3)可以转化为Assuming that the servo valve spool displacement is proportional to the control input u, that is, x v = ki u, where ki > 0 is the proportionality factor and u is the control input voltage, therefore, equation (3) can be transformed into
式中:kt=kqki表示总的流量增益。Where: k t = k q k i represents the total flow gain.
定义状态变量x1=T,那么整个系统可以写成如下状态空间形式:Define the state variable x 1 =T, then the whole system can be written in the following state space form:
定义未知参数集θ=[θ1,θ2,θ3,θ4]T,其中 θ4=B, Define the unknown parameter set θ = [θ 1 , θ 2 , θ 3 , θ 4 ] T , where θ 4 =B,
假设1:d(x,t)和d(x,t)的导数有界,即Assumption 1: The derivatives of d(x,t) and d(x,t) are bounded, i.e.
式中:ζ1、ζ2为已知常数。In the formula: ζ 1 and ζ 2 are known constants.
进一步的,所述步骤S2包括如下步骤:Further, the step S2 includes the following steps:
S2a:构建带速率限制的投影自适应律;S2a: Build a projection adaptive law with rate limit;
S2b:设计鲁棒控制器;S2b: Design a robust controller;
S2c:验证控制器自身系统的稳定性。S2c: Verify the stability of the controller's own system.
进一步的,所述步骤S2a具体为:令表示θ的估计,表示θ的估计误差,即建立一个投影函数如下:Further, the step S2a is specifically as follows: let represents the estimate of θ, represents the estimation error of θ, i.e. Create a projection function as follows:
式中:ζ∈R2,Γ(t)∈R2×2是一个随时间变化的正定对称矩阵,和分别表示Ωθ的内部和边界,表示时的外单位法向量;In the formula: ζ∈R 2 , Γ(t)∈R 2×2 is a time-varying positive definite symmetric matrix, and denote the interior and boundary of Ω θ , respectively, express The outer unit normal vector when ;
对于投影函数(8)式,在控制参数估计过程中,要用到预设的自适应限制速度,因而,建立一个饱和函数如下:For the projection function (8), in the control parameter estimation process, the preset adaptive limit speed is used, therefore, a saturation function is established as follows:
式中:是一个预先设置的限制速率,通过如下引理总结得出系统将要用到的参数估计算法的结构特性;where: is a preset limit rate, and the following lemma summarizes the structural characteristics of the parameter estimation algorithm that the system will use;
引理1:假设使用下面的投影型自适应律和预设的自适应限制速率更新估计参数 Lemma 1: Assume the following projective adaptive law and a preset adaptive limiting rate are used Update estimated parameters
式中:τ是自适应函数,Γ(t)>0是连续的可微正对称自适应率矩阵,由此自适应律,可得以下理想特性:In the formula: τ is the adaptive function, Γ(t)>0 is a continuous differentiable positive symmetric adaptive rate matrix, and from this adaptive law, the following ideal characteristics can be obtained:
P1)参数估计值总在已知有界的Ωθ集内,即对于任意t,总有因而由假设1可得 P1) The parameter estimates are always within the known bounded set of Ω θ , that is, for any t, there is always Therefore, from
P2) P2)
P3)参数变化律是一致有界的,即 P3) The law of parameter change is uniformly bounded, that is,
进一步的,所述步骤S2b具体为:定义液压缸的控制误差e=x1-x1d,假设其需满足以下性能指标:Further, the step S2b is specifically: defining the control error e=x 1 -x 1d of the hydraulic cylinder, assuming that it needs to meet the following performance indicators:
式中:δl,δu为待设计参数,用于辅助约束控制误差的上下限;ρ(t)为正的严格递增光滑函数,如下式所示:Where: δ l , δ u are parameters to be designed, which are used to assist the upper and lower limits of the control error; ρ(t) is a positive strictly increasing smooth function, as shown in the following formula:
式中:ρ0、ρ∞和k均为正的可设计参数;where: ρ 0 , ρ ∞ and k are all positive designable parameters;
式(11)中-δlρ0和δuρ0分别约束了输出力控制误差e(t)的最大下冲量和最大超调量,参数k约束了误差e(t)的收敛速度,ρ∞约束了误差的稳态界,从而,式(11)对输出力控制误差的性能给出了具体的规划,通过选择适当的参数ρ0、ρ∞、k、δl和δu等参数,便可对输出力控制误差的瞬态和稳定性能进行预先的规划,依据系统的实际需求完成瞬态性能的改进;In formula (11), -δ l ρ 0 and δ u ρ 0 constrain the maximum undershoot and maximum overshoot of the output force control error e(t), respectively, and the parameter k constrains the convergence speed of the error e(t), ρ ∞ constrains the steady-state boundary of the error, so Equation (11) gives a specific plan for the performance of the output force control error. By choosing appropriate parameters such as ρ 0 , ρ ∞ , k, δ l and δ u , The transient and stable performance of the output force control error can be planned in advance, and the transient performance can be improved according to the actual needs of the system;
建立如下递增函数:Create the following increment function:
式中:z1(t)为对应于控制误差e(t)的转换误差变量,经分析易知,式(13)等价于e(t)=ρ(t)S(z1(t)),且z1(t)有界时,预设性能特性式(8)始终满足;In the formula: z 1 (t) is the conversion error variable corresponding to the control error e(t). It is easy to know by analysis that formula (13) is equivalent to e(t)=ρ(t)S(z 1 (t) ), and when z 1 (t) is bounded, the preset performance characteristic formula (8) always satisfies;
满足特性公式(13)的递增函数S(z1)可选取如下:The increasing function S(z 1 ) satisfying the characteristic formula (13) can be selected as follows:
求取式(14)的反函数,可得:Taking the inverse function of equation (14), we can get:
再针对转换误差z1进行控制器设计;Then design the controller according to the conversion error z 1 ;
建立函数如下:The build function is as follows:
式中:k1为反馈增益;In the formula: k 1 is the feedback gain;
控制器设计如下:The controller is designed as follows:
式中:k2为反馈增益;In the formula: k 2 is the feedback gain;
将控制器(17)代入式(16)可得:Substitute controller (17) into equation (16) to get:
自适应律设计如下:The adaptive law is designed as follows:
式中: where:
然后可得:Then get:
设计鲁棒控制器us2如下:The robust controller u s2 is designed as follows:
式中:β1为待设计参数。In the formula: β 1 is the parameter to be designed.
进一步的,所述步骤S2c具体为:通过性能定理1:选取系统控制的初始条件满足-δlρ(0)<e(0)<δuρ(0),即-δl<λ(0)<δu;同时参数β1的选取满足如下不等式:Further, the step S2c is specifically: through the performance theorem 1: select the initial condition of the system control to satisfy -δ l ρ(0)<e(0)<δ u ρ(0), that is, -δ l <λ(0 )<δ u ; at the same time, the selection of parameter β 1 satisfies the following inequality:
同时设计足够大的参数k1和k2,使得如下矩阵Λ为正定矩阵At the same time, sufficiently large parameters k 1 and k 2 are designed, so that the following matrix Λ is a positive definite matrix
则可保证输出力的控制误差始终有界,输出力可实现较好的指令跟踪,且通过调节ρ0、ρ∞、k、δl和δu等参数,可保证控制误差满足式(11)所设计的预设性能需求;Then it can ensure that the control error of the output force is always bounded, the output force can achieve better command tracking, and by adjusting parameters such as ρ 0 , ρ ∞ , k, δ l and δ u , the control error can be guaranteed to satisfy the formula (11) The designed preset performance requirements;
证明如下:The proof is as follows:
定义如下Lyapunov函数:Define the following Lyapunov function:
进一步对V求导,并代入式(16)、(20),可得:Taking the derivative of V further and substituting into equations (16) and (20), we can get:
式中:Z=[z1,z2]T,矩阵Λ的定义如式(23),若通过合理的设计参数k1和k2使矩阵Λ为正定矩阵,可使下式满足:In the formula: Z=[z 1 , z 2 ] T , the definition of matrix Λ is as in formula (23). If the matrix Λ is a positive definite matrix through reasonable design parameters k 1 and k 2 , the following formula can be satisfied:
式中:λmin(Λ)表示矩阵Λ的最小特征值,分析式(26)可知Lyapunov函数有界,同时W积分有界,进而可知转换误差量z1和z2均有界,结合式(8)、(16)和(17)可知,系统中所有信号均有界,从而可知W的导数有界,由Barbalat引理可知,当时间趋于无穷大时,W趋近于零,也即转换误差量z1趋近于零,从而结合式(11)可知控制误差e(t)始终有界,进而证明控制器是收敛的,系统是稳定的。In the formula: λ min (Λ) represents the minimum eigenvalue of the matrix Λ. By analyzing formula (26), it can be seen that the Lyapunov function is bounded, and the W integral is bounded, and then it can be known that the conversion error quantities z 1 and z 2 are bounded, and the combined formula ( 8), (16) and (17), it can be known that all signals in the system are bounded, so the derivative of W is bounded. According to Barbalat's lemma, when time tends to infinity, W tends to zero, that is, the conversion The amount of error z 1 tends to zero, so combined with equation (11), it can be known that the control error e(t) is always bounded, which proves that the controller is convergent and the system is stable.
进一步的,所述步骤三具体为调节控制律u的参数k1、k2、ρ0、ρ∞、k、δl、δu、β1、Γ使系统满足控制性能指标。Further, the third step is to adjust the parameters k 1 , k 2 , ρ 0 , ρ ∞ , k, δ l , δ u , β 1 , and Γ of the control law u to make the system meet the control performance index.
与现有技术相比,本发明的有益效果是:本发明针对液压伺服系统的特点,建立了液压伺服系统模型和设计的液压系统输出力预设性能自适应渐近控制器,通过鲁棒误差符号积分函数对未建模干扰进行估计并进行前馈补偿,同时使用自适应控制器对系统参数进行估计,能有效解决电机伺服系统不确定非线性和参数不确定问题,并基于预设性能函数设计了预设性能控制器,最终通过李雅普诺夫证明了系统总体的稳定性,在现有干扰条件下,参数收敛好,系统的瞬态和稳态控制精度满足性能指标;同时本发明简化了控制器设计,最终仿真结果表明了其具有有效性,更利于在工程实际中应用。Compared with the prior art, the beneficial effects of the present invention are: aiming at the characteristics of the hydraulic servo system, the present invention establishes the hydraulic servo system model and the designed hydraulic system output force preset performance adaptive asymptotic controller, through the robust error The symbolic integral function estimates the unmodeled disturbance and performs feed-forward compensation. At the same time, the adaptive controller is used to estimate the system parameters, which can effectively solve the uncertain nonlinear and parameter uncertain problems of the motor servo system. Based on the preset performance function The preset performance controller is designed, and Lyapunov finally proves the overall stability of the system. Under the existing interference conditions, the parameters converge well, and the transient and steady-state control accuracy of the system meets the performance index; at the same time, the invention simplifies the The controller is designed, and the final simulation results show that it is effective and more conducive to the application in engineering practice.
附图说明Description of drawings
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to explain the technical solutions of the embodiments of the present invention more clearly, the following briefly introduces the accompanying drawings that are used in the description of the embodiments. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, other drawings can also be obtained from these drawings without any creative effort.
图1为本发明流程图;Fig. 1 is the flow chart of the present invention;
图2为本发明双出杆液压缸系统示意图;Fig. 2 is the schematic diagram of the double-rod hydraulic cylinder system of the present invention;
图3为本发明控制误差预设性能示意图;3 is a schematic diagram of the preset performance of the control error of the present invention;
图4为本发明函数S(z1)示意图;4 is a schematic diagram of the function S(z 1 ) of the present invention;
图5为本发明实施例一参数估计曲线图;5 is a parameter estimation curve diagram according to an embodiment of the present invention;
图6为本发明实施例一两种控制器控制误差对比曲线图。FIG. 6 is a comparison curve diagram of the control errors of two controllers according to an embodiment of the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.
请参阅图1-2,本发明提供一种技术方案:一种液压力预设性能自适应渐近控制方法,包括如下步骤:1-2, the present invention provides a technical solution: a hydraulic pressure preset performance adaptive asymptotic control method, including the following steps:
S1:建立双出杆液压缸伺服系统模型;S1: Establish the servo system model of the double-rod hydraulic cylinder;
S2:设计液压力预设性能自适应渐近控制器;S2: Design hydraulic pressure preset performance adaptive asymptotic controller;
S3:调节步骤S2中的参数,使系统满足控制性能指标。S3: Adjust the parameters in step S2 so that the system meets the control performance index.
其中,步骤S1具体为:根据牛顿第二定律,双出杆液压缸惯性负载的动力学模型方程为:Wherein, step S1 is specifically: according to Newton's second law, the dynamic model equation of the inertial load of the double-rod hydraulic cylinder is:
式中:y为负载位移,T表示驱动力,PL=P1-P2是负载驱动压力,P1和P2分别为液压缸两腔压力,A为活塞杆有效工作面积,b代表粘性摩擦系数,f为其他未建模干扰,包括非线性摩擦、外部干扰以及未建模动态,液压缸负载压力动态方程为:In the formula: y is the load displacement, T is the driving force, P L = P 1 -P 2 is the load driving pressure, P 1 and P 2 are the pressures of the two chambers of the hydraulic cylinder, A is the effective working area of the piston rod, and b is the viscosity Friction coefficient, f is other unmodeled disturbances, including nonlinear friction, external disturbances and unmodeled dynamics, the dynamic equation of hydraulic cylinder load pressure is:
式中:Vt为液压缸两腔总有效容积,Ct为液压缸泄露系数,QL=(Q1+Q2)/2是负载流量,Q1为液压缸进油腔供油流量,Q2为液压缸回油腔回油流量,q(t)为建模误差及未建模动态;In the formula: V t is the total effective volume of the two chambers of the hydraulic cylinder, C t is the leakage coefficient of the hydraulic cylinder, Q L = (Q 1 +Q 2 )/2 is the load flow, Q 1 is the oil supply flow of the hydraulic cylinder into the oil chamber, Q 2 is the return oil flow of the hydraulic cylinder oil return chamber, q(t) is the modeling error and unmodeled dynamics;
QL为伺服阀阀芯位移xv的函数:Q L is a function of servo valve spool displacement x v :
式中:为流量伺服阀的增益系数,Cd为伺服阀的流量系数,w为伺服阀的面积梯度;ρ为液压油的密度,Ps为供油压力,sign(xv)为:where: is the gain coefficient of the flow servo valve, C d is the flow coefficient of the servo valve, w is the area gradient of the servo valve; ρ is the density of the hydraulic oil, P s is the oil supply pressure, and sign(x v ) is:
假设伺服阀阀芯位移正比于控制输入u,即,xv=kiu,其中ki>0是比例系数,u是控制输入电压,因此,等式(3)可以转化为Assuming that the servo valve spool displacement is proportional to the control input u, that is, x v = ki u, where ki > 0 is the proportionality factor and u is the control input voltage, therefore, equation (3) can be transformed into
式中:kt=kqki表示总的流量增益。Where: k t = k q k i represents the total flow gain.
定义状态变量x1=T,那么整个系统可以写成如下状态空间形式:Define the state variable x 1 =T, then the whole system can be written in the following state space form:
定义未知参数集θ=[θ1,θ2,θ3,θ4]T,其中 θ4=B, Define the unknown parameter set θ = [θ 1 , θ 2 , θ 3 , θ 4 ] T , where θ 4 =B,
假设1:d(x,t)和d(x,t)的导数有界,即Assumption 1: The derivatives of d(x,t) and d(x,t) are bounded, i.e.
式中:ζ1、ζ2为已知常数。In the formula: ζ 1 and ζ 2 are known constants.
步骤S2包括如下步骤:Step S2 includes the following steps:
S2a:构建带速率限制的投影自适应律;S2a: Build a projection adaptive law with rate limit;
S2b:设计鲁棒控制器;S2b: Design a robust controller;
S2c:验证控制器自身系统的稳定性。S2c: Verify the stability of the controller's own system.
步骤S2a具体为:令表示θ的估计,表示θ的估计误差,即建立一个投影函数如下:Step S2a is specifically as follows: let represents the estimate of θ, represents the estimation error of θ, i.e. Create a projection function as follows:
式中:ζ∈R2,Γ(t)∈R2×2是一个随时间变化的正定对称矩阵,和分别表示Ωθ的内部和边界,表示时的外单位法向量;In the formula: ζ∈R 2 , Γ(t)∈R 2×2 is a time-varying positive definite symmetric matrix, and denote the interior and boundary of Ω θ , respectively, express The outer unit normal vector when ;
对于投影函数(8)式,在控制参数估计过程中,要用到预设的自适应限制速度,因而,建立一个饱和函数如下:For the projection function (8), in the control parameter estimation process, the preset adaptive limit speed is used, therefore, a saturation function is established as follows:
式中:是一个预先设置的限制速率,通过如下引理总结得出系统将要用到的参数估计算法的结构特性;where: is a preset limit rate, and the following lemma summarizes the structural characteristics of the parameter estimation algorithm that the system will use;
引理1:假设使用下面的投影型自适应律和预设的自适应限制速率更新估计参数 Lemma 1: Assume the following projective adaptive law and a preset adaptive limiting rate are used Update estimated parameters
式中:τ是自适应函数,Γ(t)>0是连续的可微正对称自适应率矩阵,由此自适应律,可得以下理想特性:In the formula: τ is the adaptive function, Γ(t)>0 is a continuous differentiable positive symmetric adaptive rate matrix, and from this adaptive law, the following ideal characteristics can be obtained:
P1)参数估计值总在已知有界的Ωθ集内,即对于任意t,总有因而由假设1可得 P1) The parameter estimates are always within the known bounded set of Ω θ , that is, for any t, there is always Therefore, from
P2) P2)
P3)参数变化律是一致有界的,即 P3) The law of parameter change is uniformly bounded, that is,
步骤S2b具体为:定义液压缸的控制误差e=x1-x1d,假设其需满足以下性能指标:Step S2b is specifically: defining the control error e=x 1 -x 1d of the hydraulic cylinder, assuming that it needs to meet the following performance indicators:
式中:δl,δu为待设计参数,用于辅助约束控制误差的上下限;ρ(t)为正的严格递增光滑函数,如下式所示:Where: δ l , δ u are parameters to be designed, which are used to assist the upper and lower limits of the control error; ρ(t) is a positive strictly increasing smooth function, as shown in the following formula:
式中:ρ0、ρ∞和k均为正的可设计参数;性能指标不等式(11)的大致曲线如图3所示;where: ρ 0 , ρ ∞ and k are all positive designable parameters; the approximate curve of the performance index inequality (11) is shown in Figure 3;
式(11)中-δlρ0和δuρ0分别约束了输出力控制误差e(t)的最大下冲量和最大超调量,参数k约束了误差e(t)的收敛速度,ρ∞约束了误差的稳态界,从而,式(11)对输出力控制误差的性能给出了具体的规划,通过选择适当的参数ρ0、ρ∞、k、δl和δu等参数,便可对输出力控制误差的瞬态和稳定性能进行预先的规划,依据系统的实际需求完成瞬态性能的改进;In formula (11), -δ l ρ 0 and δ u ρ 0 constrain the maximum undershoot and maximum overshoot of the output force control error e(t), respectively, and the parameter k constrains the convergence speed of the error e(t), ρ ∞ constrains the steady-state boundary of the error, so Equation (11) gives a specific plan for the performance of the output force control error. By choosing appropriate parameters such as ρ 0 , ρ ∞ , k, δ l and δ u , The transient and stable performance of the output force control error can be planned in advance, and the transient performance can be improved according to the actual needs of the system;
建立如下递增函数:Create the following increment function:
式中:z1(t)为对应于控制误差e(t)的转换误差变量,经分析易知,式(13)等价于e(t)=ρ(t)S(z1(t)),且z1(t)有界时,预设性能特性式(8)始终满足;In the formula: z 1 (t) is the conversion error variable corresponding to the control error e(t). It is easy to know by analysis that formula (13) is equivalent to e(t)=ρ(t)S(z 1 (t) ), and when z 1 (t) is bounded, the preset performance characteristic formula (8) always satisfies;
满足特性公式(13)的递增函数S(z1)可选取如下:The increasing function S(z 1 ) satisfying the characteristic formula (13) can be selected as follows:
递增函数S(z1)的曲线如图4所示;The curve of the increasing function S(z 1 ) is shown in Figure 4;
求取式(14)的反函数,可得:Taking the inverse function of equation (14), we can get:
再针对转换误差z1进行控制器设计;Then design the controller according to the conversion error z 1 ;
建立函数如下:The build function is as follows:
式中:k1为反馈增益;In the formula: k 1 is the feedback gain;
控制器设计如下:The controller is designed as follows:
式中:k2为反馈增益;In the formula: k 2 is the feedback gain;
将控制器(17)代入式(16)可得:Substitute controller (17) into equation (16) to get:
自适应律设计如下:The adaptive law is designed as follows:
式中: where:
然后可得:Then get:
设计鲁棒控制器us2如下:The robust controller u s2 is designed as follows:
式中:β1为待设计参数。In the formula: β 1 is the parameter to be designed.
步骤S2c具体为:通过性能定理1:选取系统控制的初始条件满足-δlρ(0)<e(0)<δuρ(0),即-δl<λ(0)<δu;同时参数β1的选取满足如下不等式:Step S2c is specifically: through the performance theorem 1: select the initial condition of the system control to satisfy -δ l ρ(0)<e(0)<δ u ρ(0), namely -δ l <λ(0)<δ u ; At the same time, the selection of parameter β 1 satisfies the following inequality:
同时设计足够大的参数k1和k2,使得如下矩阵Λ为正定矩阵At the same time, sufficiently large parameters k 1 and k 2 are designed, so that the following matrix Λ is a positive definite matrix
则可保证输出力的控制误差始终有界,输出力可实现较好的指令跟踪,且通过调节ρ0、ρ∞、k、δl和δu等参数,可保证控制误差满足式(11)所设计的预设性能需求;Then it can ensure that the control error of the output force is always bounded, the output force can achieve better command tracking, and by adjusting parameters such as ρ 0 , ρ ∞ , k, δ l and δ u , the control error can be guaranteed to satisfy the formula (11) The designed preset performance requirements;
证明如下:The proof is as follows:
定义如下Lyapunov函数:Define the following Lyapunov function:
进一步对V求导,并代入式(16)、(20),可得:Taking the derivative of V further and substituting into equations (16) and (20), we can get:
式中:Z=[z1,z2]T,矩阵Λ的定义如式(23),若通过合理的设计参数k1和k2使矩阵Λ为正定矩阵,可使下式满足:In the formula: Z=[z 1 , z 2 ] T , the definition of matrix Λ is as in formula (23). If the matrix Λ is a positive definite matrix through reasonable design parameters k 1 and k 2 , the following formula can be satisfied:
式中:λmin(Λ)表示矩阵Λ的最小特征值,分析式(26)可知Lyapunov函数有界,同时W积分有界,进而可知转换误差量z1和z2均有界,结合式(8)、(16)和(17)可知,系统中所有信号均有界,从而可知W的导数有界,由Barbalat引理可知,当时间趋于无穷大时,W趋近于零,也即转换误差量z1趋近于零,从而结合式(11)可知控制误差e(t)始终有界,进而证明控制器是收敛的,系统是稳定的。In the formula: λ min (Λ) represents the minimum eigenvalue of the matrix Λ. By analyzing formula (26), it can be seen that the Lyapunov function is bounded, and the W integral is bounded, and then it can be known that the conversion error quantities z 1 and z 2 are bounded, and the combined formula ( 8), (16) and (17), it can be known that all signals in the system are bounded, so the derivative of W is bounded. According to Barbalat's lemma, when time tends to infinity, W tends to zero, that is, the conversion The amount of error z 1 tends to zero, so combined with equation (11), it can be known that the control error e(t) is always bounded, which proves that the controller is convergent and the system is stable.
步骤三具体为调节控制律u的参数k1、k2、ρ0、ρ∞、k、δl、δu、β1、Γ使系统满足控制性能指标。The third step is to adjust the parameters k 1 , k 2 , ρ 0 , ρ ∞ , k, δ l , δ u , β 1 , and Γ of the control law u to make the system meet the control performance index.
实施例一:Example 1:
在仿真中取如下参数对系统进行建模:m=40kg,A=2×10-4m2,B=80N·s/m,βe=200Mpa,V01=1×10-3m3,V02=1×10-3m3,Ct=9×10-12m5/Ns,θmin=[0.01,0.1,1000,20]T,θmax=[10,30,500000,500]T,Г=diag{6,10,390,50},所选取的远离于参数的真值,以考核自适应控制律的效果,力输入信号单位KN;系统所加干扰为f=30sin(2πt)KN。In the simulation, the following parameters are taken to model the system: m=40kg, A=2×10 -4 m 2 , B=80N·s/m, β e =200Mpa, V 01 =1×10 -3 m 3 , V 02 =1×10 -3 m 3 , C t =9×10 -12 m 5 /Ns, θ min = [0.01, 0.1, 1000, 20] T , θ max = [10, 30, 500000, 500] T , Г=diag{6,10,390,50}, the selected away from the true value of the parameter to evaluate the effect of the adaptive control law, force the input signal The unit is KN; the interference added by the system is f=30sin(2πt)KN.
本发明提出的基于鲁棒误差符号积分的液压力自适应预设性能控制器(APFRISE),控制器设计参数控制器相关参数选取为:k1=200,k2=300,β=100,ρ0=800,ρ∞=500,k=0.0001,δl=δu=1;PID控制器参数为kp=500,ki=100,kd=0;The hydraulic pressure adaptive preset performance controller (APFRISE) based on the robust error sign integration proposed by the present invention, the controller design parameters The relevant parameters of the controller are selected as: k 1 =200, k 2 =300, β=100, ρ 0 = 800, ρ ∞ = 500, k = 0.0001, δ l = δ u = 1; PID controller parameters are k p = 500, ki = 100, k d = 0;
控制律作用效果如图5和图6所示,图5为参数估计曲线;图6为设计的控制器(APFRISE)和PID控制器的控制误差对比曲线。The effect of the control law is shown in Figure 5 and Figure 6, Figure 5 is the parameter estimation curve; Figure 6 is the control error comparison curve of the designed controller (APFRISE) and the PID controller.
由上图可知,本发明提出的算法在仿真环境下能够准确的估计系统参数,相比PID控制器,本发明设计的控制器能够取得良好的控制精度,能保证系统的预设瞬态和稳态控制精度要求。It can be seen from the above figure that the algorithm proposed by the present invention can accurately estimate the system parameters in the simulation environment. Compared with the PID controller, the controller designed by the present invention can achieve good control accuracy, and can ensure the preset transient and stable conditions of the system. state control accuracy requirements.
在本说明书的描述中,参考术语“一个实施例”、“示例”、“具体示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。In the description of this specification, description with reference to the terms "one embodiment," "example," "specific example," etc. means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one aspect of the present invention. in one embodiment or example. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
以上公开的本发明优选实施例只是用于帮助阐述本发明。优选实施例并没有详尽叙述所有的细节,也不限制该发明仅为所述的具体实施方式。显然,根据本说明书的内容,可作很多的修改和变化。本说明书选取并具体描述这些实施例,是为了更好地解释本发明的原理和实际应用,从而使所属技术领域技术人员能很好地理解和利用本发明。本发明仅受权利要求书及其全部范围和等效物的限制。The above-disclosed preferred embodiments of the present invention are provided only to help illustrate the present invention. The preferred embodiments do not exhaust all the details, nor do they limit the invention to only the described embodiments. Obviously, many modifications and variations are possible in light of the contents of this specification. These embodiments are selected and described in this specification in order to better explain the principles and practical applications of the present invention, so that those skilled in the art can well understand and utilize the present invention. The present invention is to be limited only by the claims and their full scope and equivalents.
Claims (3)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010710755.6A CN111781836B (en) | 2020-07-22 | 2020-07-22 | An adaptive asymptotic control method for hydraulic pressure preset performance |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010710755.6A CN111781836B (en) | 2020-07-22 | 2020-07-22 | An adaptive asymptotic control method for hydraulic pressure preset performance |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111781836A CN111781836A (en) | 2020-10-16 |
CN111781836B true CN111781836B (en) | 2022-05-27 |
Family
ID=72763813
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010710755.6A Active CN111781836B (en) | 2020-07-22 | 2020-07-22 | An adaptive asymptotic control method for hydraulic pressure preset performance |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111781836B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115629543B (en) * | 2022-10-21 | 2023-04-07 | 东北林业大学 | A method, system, and storage medium for intervening MDF continuous flattening three-branch decision-making collaborative control |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101846975A (en) * | 2010-05-28 | 2010-09-29 | 北京理工大学 | Servo system self-adaptive robust controller with dynamic frictional compensation |
CN104834220A (en) * | 2015-05-20 | 2015-08-12 | 南京理工大学 | Adaptive error symbol integration robust repetitive control method for electromechanical servo system |
CN105700347A (en) * | 2014-12-15 | 2016-06-22 | 南京理工大学 | Hydraulic motor preset performance tracking control method with hysteresis compensation |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106533299B (en) * | 2016-12-27 | 2018-10-02 | 山东大学 | Back-emf, the method for servo-controlling that load torque disturbs and its system are eliminated simultaneously |
-
2020
- 2020-07-22 CN CN202010710755.6A patent/CN111781836B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101846975A (en) * | 2010-05-28 | 2010-09-29 | 北京理工大学 | Servo system self-adaptive robust controller with dynamic frictional compensation |
CN105700347A (en) * | 2014-12-15 | 2016-06-22 | 南京理工大学 | Hydraulic motor preset performance tracking control method with hysteresis compensation |
CN104834220A (en) * | 2015-05-20 | 2015-08-12 | 南京理工大学 | Adaptive error symbol integration robust repetitive control method for electromechanical servo system |
Non-Patent Citations (1)
Title |
---|
王树波 等."基于预设性能转台伺服系统的".《北京理工大学学报》.2019,第39卷(第2期),第193-197页. * |
Also Published As
Publication number | Publication date |
---|---|
CN111781836A (en) | 2020-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112415891B (en) | Adaptive output feedback asymptotic control method for electro-hydraulic servo system | |
CN110794687B (en) | Electro-hydraulic servo system self-adaptive state constraint control method based on interference compensation | |
Feng et al. | Identification and compensation of non-linear friction for a electro-hydraulic system | |
CN112486021B (en) | A Low-Complexity Control Method for Asymmetric Servo-Hydraulic Position Tracking System | |
CN108303895B (en) | Multi-model robust self-adaptive control method for electro-hydraulic position servo system | |
CN105700347B (en) | A hydraulic motor preset performance tracking control method with hysteresis compensation | |
CN104317198B (en) | Nonlinear Robust Position Control Method for Electro-hydraulic Servo System with Time-Varying Output Constraints | |
CN108873702A (en) | A kind of linear active disturbance rejection control method and device of electro-hydraulic position servo control system | |
CN104898428B (en) | A kind of electrohydraulic servo system adaptive robust control method based on Interference Estimation | |
Shen et al. | Fuzzy sliding mode control with state estimation for velocity control system of hydraulic cylinder using a new hydraulic transformer | |
CN106154833B (en) | A kind of electro-hydraulic load simulator output feedback ontrol method | |
CN113359485B (en) | Method for controlling output feedback preset performance of electro-hydraulic servo system | |
CN110456646B (en) | Design method of anchor rod drilling machine drill arm swing angle self-adaptive sliding mode control system based on EKF | |
CN114967466B (en) | A preset performance output feedback control method for an electro-hydraulic servo system based on disturbance compensation and full-state constraints | |
CN110308651B (en) | Electro-hydraulic servo system all-state constraint control method based on extended state observer | |
CN107165892B (en) | A kind of sliding-mode control of electrohydraulic servo system | |
CN108869420A (en) | A kind of adaptive backstepping control method of electrohydraulic servo system filtered based on instruction | |
CN108958023A (en) | A kind of electro-hydraulic position servo control system, terminal, storage medium | |
CN106483844A (en) | The implementation method of the electrohydraulic servo system adaptive location controller based on non linear robust | |
CN109426150B (en) | Load simulator backstepping control method based on extended state observer | |
CN111308889B (en) | Adaptive integral robust control method of spray rod system | |
CN105093936B (en) | The continuous sliding-mode control of finite time of electro-hydraulic position servo system | |
CN115143165A (en) | Electro-hydraulic proportional servo valve position shaft control method considering unknown dynamic compensation of system | |
CN105068426B (en) | The continuous sliding-mode control of electro-hydraulic position servo system based on interference compensation | |
CN111781836B (en) | An adaptive asymptotic control method for hydraulic pressure preset performance |
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 |