JP2023517140A - A Low-Complexity Control Method for Asymmetric Servohydraulic Position Tracking Systems - Google Patents

A Low-Complexity Control Method for Asymmetric Servohydraulic Position Tracking Systems Download PDF

Info

Publication number
JP2023517140A
JP2023517140A JP2022579147A JP2022579147A JP2023517140A JP 2023517140 A JP2023517140 A JP 2023517140A JP 2022579147 A JP2022579147 A JP 2022579147A JP 2022579147 A JP2022579147 A JP 2022579147A JP 2023517140 A JP2023517140 A JP 2023517140A
Authority
JP
Japan
Prior art keywords
bounded
servo
max
hydraulic system
hydraulic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2022579147A
Other languages
Japanese (ja)
Other versions
JP7448692B2 (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.)
Yanshan University
Original Assignee
Yanshan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Yanshan University filed Critical Yanshan University
Publication of JP2023517140A publication Critical patent/JP2023517140A/en
Application granted granted Critical
Publication of JP7448692B2 publication Critical patent/JP7448692B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Abstract

本発明は、非対称サーボ油圧位置追跡システムのための低複雑度制御方法を開示し、サーボ油圧位置の制御分野に属し、該方法は、単一出力ロッドのサーボ油圧システムのモデルを確立するステップと、単一出力ロッドのサーボ油圧システムのモデルに基づいて、低複雑度制御ポリシーを用いて単一出力ロッドのサーボ油圧システムのコントローラを設計するステップと、単一出力ロッドのサーボ油圧システムのコントローラ及び単一出力ロッドのサーボ油圧システムのモデルに基づいて、単一出力ロッドのサーボ油圧システムの安定性を証明するステップと、を含み、油圧システムに存在する様々な不確実性問題及び未知の非線形問題を解決し、コントローラの設計は、正確な数学的モデルに依存せず、測定可能な状態信号のみが必要であり、制御率の計算は、既存のバックステッピング適応に基づいて開発されたアルゴリズムと比較して、計算プロセスが簡単であり、計算量が少なく、リアルタイム制御が容易であり、工学的実現が容易である。The present invention discloses a low-complexity control method for an asymmetric servo-hydraulic position tracking system and belongs to the field of servo-hydraulic position control, the method comprising the steps of establishing a model of a single output rod servo-hydraulic system; , designing a controller of the single output rod servohydraulic system with a low-complexity control policy based on a model of the single output rod servohydraulic system; proving the stability of the single-output rod servo-hydraulic system based on a model of the single-output rod servo-hydraulic system, including various uncertainties and unknown non-linear problems existing in the hydraulic system. , the controller design does not rely on an exact mathematical model, only measurable state signals are required, and the calculation of the control ratio is comparable to algorithms developed based on existing backstepping adaptations. Therefore, the calculation process is simple, the amount of calculation is small, the real-time control is easy, and the engineering implementation is easy.

Description

本発明は、サーボ油圧位置制御の分野に関し、特に非対称サーボ油圧位置追跡システムのための低複雑度制御方法に関する。 The present invention relates to the field of servohydraulic position control, and more particularly to a low complexity control method for asymmetric servohydraulic position tracking systems.

サーボ油圧システムは、負荷効率に優れ、電力比が小さく、応答速度が速いなどの利点を有するため、現代工業で広く応用されている(例えば、アクティブサスペンションシステム、油圧ショベル、マニピュレータ)。サーボ油圧システムは、二重出力ロッドの対称油圧アクチュエータ機構と単一出力ロッドの非対称油圧アクチュエータ機構との2種の油圧アクチュエータ機構を含む。アクティブサスペンション制御の研究では、2つのオイルキャビティの面積が等しいと見なされることが多く、対称油圧アクチュエータ機構が使用されている。非対称サーボ油圧システムについて、ピストンロッドが存在するため、油圧シリンダのロッドキャビティとロッドレスキャビティとの有効面積が等しくなくなり、同じ性能を実現できる対称油圧アクチュエータ機構は、常に非対称油圧アクチュエータ機構よりも体積が大きく、一般的に少数の特別な場合に応用されており、工業で広く応用されているのは非対称油圧アクチュエータ機構である。従って、空間が大きく制限されている車両サスペンションシステムに対称油圧アクチュエータ機構を配置することは不合理である。 Servo-hydraulic systems have advantages such as high load efficiency, small power ratio, and fast response speed, so they are widely applied in modern industry (eg, active suspension systems, hydraulic excavators, manipulators). Servohydraulic systems include two types of hydraulic actuator mechanisms: dual output rod symmetric hydraulic actuator mechanisms and single output rod asymmetric hydraulic actuator mechanisms. Active suspension control studies often assume that the areas of the two oil cavities are equal and use symmetrical hydraulic actuator mechanisms. For the asymmetric servo hydraulic system, due to the presence of the piston rod, the effective area of the rod cavity and the rodless cavity of the hydraulic cylinder is no longer equal, and the symmetric hydraulic actuator mechanism that can achieve the same performance always has more volume than the asymmetric hydraulic actuator mechanism. Large, generally applied in a few special cases, and widely applied in industry are asymmetric hydraulic actuator mechanisms. Therefore, it is unreasonable to place a symmetrical hydraulic actuator mechanism in a vehicle suspension system where space is highly restricted.

しかしながら、非対称サーボ油圧システムは、1つの複雑な非線形システムであり、同時に、負荷変化、パラメータの不確実性及び未知の非線形などの様々な不確実性に関し、モデルの確立及びコントローラの設計が困難になり、複雑になり、従って、非対称サーボ油圧システムについて、高精度の位置制御は大きな課題に直面している。パラメータの不確実性は、主に漏れ係数、油膜粘度及び粘性摩擦係数を含み、未知の非線形は、主にスプールの不感帯及び外乱を含み、これらの要因は、高性能コントローラの開発を妨げている。位置追跡システムの操作性能を満たすために、ニューラルネットワーク適応、ファジー論理制御、ロバスト適応制御、バックステッピング適応制御などの多くの制御技術が過去10年間に開発されている。これらの制御方式では、適応制御は、そのオンライン学習能力により、パラメータの不確実性の問題を処理でき、ロバスト適応制御は、未知の外乱の問題を処理できる。上記適応制御に基づくアルゴリズムは、良好なシミュレーション結果を取得できるが、1つの重要な問題としては、これらのアルゴリズムは計算負荷が高く、収束に達するまでに長い時間がかかるため、油圧システムでは、適応に基づく上記制御アルゴリズムを応用することは困難である。また、サーボ油圧システムの既存の制御方法のほとんどは、バックステッピング制御に基づいて開発されている。周知のように、高次システムのバックステッピング設計プロセスで、仮想制御関数の導出が繰り返されるため、計算爆発の問題が発生する。従って、設計のシステムへの依存性が低く、計算負荷が低く、関数近似機能がない新規の制御ポリシーを見つける必要がある。 However, an asymmetric servohydraulic system is a complex nonlinear system, and at the same time various uncertainties such as load changes, parameter uncertainties and unknown nonlinearities make model establishment and controller design difficult. , and complexity, and thus high-precision position control faces a major challenge for asymmetric servo-hydraulic systems. Parameter uncertainties mainly include leakage coefficient, oil film viscosity and viscous friction coefficient, and unknown non-linearities mainly include spool dead zone and disturbance, these factors hinder the development of high performance controllers. . Many control techniques, such as neural network adaptation, fuzzy logic control, robust adaptive control, and backstepping adaptive control, have been developed in the last decade to meet the operational performance of position tracking systems. In these control schemes, adaptive control can handle the problem of parameter uncertainty due to its online learning capability, and robust adaptive control can handle the problem of unknown disturbances. Algorithms based on the above adaptive control can obtain good simulation results, but one important problem is that these algorithms are computationally intensive and take a long time to reach convergence. It is difficult to apply the above control algorithm based on Also, most of the existing control methods for servo hydraulic systems are developed based on backstepping control. As is well known, in the backstepping design process of higher-order systems, the derivation of virtual control functions is repeated, causing the problem of computational explosion. Therefore, there is a need to find new control policies whose design is less system dependent, less computationally intensive, and free of function approximation capabilities.

また、サーボ油圧位置追跡システムでは、実際の応用観点から見ると、上記方法で潜在的に悪い過渡応答が存在する(オーバーシュートが大きすぎ、収束が遅い)場合、追跡性能が低下し、危険が発生し、ひいてはハードウェアを損傷する可能性がある。Bechlioulis CP、Rovithakis GAらは、最初に、追跡誤差の過渡性能を確保できる新規のコントローラ設計方法を提案している。該方法は、指定された性能関数を導入することにより、追跡誤差の過渡と定常状態との性能を指定し、誤差変換を導入することにより、元の制限のあるシステムを無制限のシステムと同等する。この思想は、さらに高次の非線形フォールトトレラント制御システム、油圧サーボシステム、サスペンションシステムに応用されている。しかしながら、指定された性能関数を使用する上記文献では、制御方案は、いずれも適応制御方式と組み合わせてシステムにおける未知の動力学問題を処理する。 In addition, in the servo hydraulic position tracking system, from a practical application point of view, if there is a potentially bad transient response (too large overshoot and slow convergence) in the above method, the tracking performance will be degraded and dangerous. can occur and eventually damage the hardware. Bechliolis CP, Rovithakis GA et al. first proposed a novel controller design method that can ensure the transient performance of the tracking error. The method specifies the transient and steady-state performance of the tracking error by introducing a specified performance function, and equates the original bounded system to the unbounded system by introducing an error transform. . This idea has been applied to higher-order nonlinear fault-tolerant control systems, hydraulic servo systems, and suspension systems. However, in the above references using a specified performance function, all of the control strategies are combined with adaptive control schemes to handle unknown dynamics problems in the system.

以上より、サーボ油圧システム及びその位置制御方法には以下の欠陥が存在し、一つ目は、単一出力ロッドの非線形サーボ油圧システムは、典型的な複雑な非線形システムであり、負荷変化、パラメータの不確実性及び未知の非線形などの問題に関するため、理論モデルの確立と実際のシステムとの間に大きなギャップがあり、コントローラの設計が困難になり、複雑になり、二つ目は、現在存在するサーボ油圧位置制御方式のほとんどは、バックステッピング制御適応に基づいて開発されているが、これらの制御方式のモデルへの依存性が高く、そして、高次システムのバックステッピング設計プロセスで、計算量が大きく、不確実性パラメータのオンライン学習が、実際の試験中のリアルタイム制御に不利であり、三つ目は、実際の応用観点から見ると、サーボ油圧位置追跡システムでは、潜在的に悪い過渡応答が存在する(オーバーシュートが大きすぎ、収束が遅い)場合、追跡性能が低下し、危険が発生し、ひいてはハードウェアを損傷する可能性がある。 From the above, the servo hydraulic system and its position control method have the following defects. There is a big gap between the establishment of the theoretical model and the actual system due to the uncertainties of and unknown nonlinear issues, which makes the design of the controller difficult and complicated. Most of the existing servo hydraulic position control schemes are developed based on backstepping control adaptation, but these control schemes are highly model dependent, and the backstepping design process of higher-order systems requires a large amount of computation. is large, the uncertainty of online learning of parameters is disadvantageous for real-time control during actual testing; If is present (too much overshoot and slow convergence), tracking performance is degraded and can be dangerous and can even damage the hardware.

従来技術に存在する課題に基づき、本発明は、単一出力ロッドのサーボ油圧システムのモデルを確立するステップS1と、単一出力ロッドのサーボ油圧システムのモデルに基づいて、低複雑度制御ポリシーを用いて単一出力ロッドのサーボ油圧システムのコントローラを設計するステップS2と、単一出力ロッドのサーボ油圧システムのコントローラ及び単一出力ロッドのサーボ油圧システムのモデルに基づいて、単一出力ロッドのサーボ油圧システムの安定性を証明するステップS3と、を含む非対称サーボ油圧位置追跡システムのための低複雑度制御方法を開示している。 Based on the problems existing in the prior art, the present invention provides a step S1 of establishing a model of a single output rod servo-hydraulic system, and based on the model of the single output rod servo-hydraulic system, a low-complexity control policy. step S2 of designing the controller of the single output rod servo-hydraulic system using the controller of the single output rod servo-hydraulic system and the model of the single output rod servo-hydraulic system based on the model of the single output rod servo A low-complexity control method for an asymmetric servo-hydraulic position tracking system is disclosed, including step S3 of verifying the stability of the hydraulic system.

さらに、単一出力ロッドのサーボ油圧システムのモデルを確立する前記プロセスは以下の通り、
ニュートンの第2法則に基づいて油圧シリンダの動力学モデルを確立し、

Figure 2023517140000002
(ただし、f(t)は様々な干渉を表し、xとmとはそれぞれ負荷の位置と質量とを表し、Bは粘性減衰係数であり、Kは負荷の等価ばね剛性であり、負荷が慣性負荷の場合、K=0であり、F=A*P-A*Pは油圧アクチュエータにより出力される主動力であり、ただし、P、Pは油圧シリンダの大きなキャビティと小さなキャビティとの圧力であり、A、Aは大きなキャビティと小さなキャビティとのピストンの有効面積である。)
3ポジション5ポートサーボバルブを用い、負荷圧力の動力学は下記式で表され、
Figure 2023517140000003
Figure 2023517140000004
(式中、V、Vはそれぞれロッドキャビティとロッドレスキャビティとの容積であり、V=V01+Axであり、V=V02-Axであり、V01とV02とはそれぞれピストンが初期位置にあるときのロッドレスキャビティとロッドキャビティとの容積であり、Cは油圧シリンダの内部漏れ係数であり、Cは油圧シリンダの外部漏れ係数であり、βはオイルの弾性率であり、Qはロッドキャビティの作動油の流量であり、Qはロッドレスキャビティの作動油の流量である。)
ただし、
Figure 2023517140000005
Figure 2023517140000006
(ただし、
Figure 2023517140000007
は流量ゲインであり、s(Γ(x))の式は
Figure 2023517140000008
であり、
は油圧システムの給油圧力であり、Pは油圧システムの油戻し圧力であり、Cはオリフィスの流量係数であり、wはスライドバルブの面積勾配であり、ρはオイルの密度であり、xはサーボバルブのスプールである。)
電圧又は電流入力uによりサーボバルブのスプールの変位xを制御して、さらに必要な対応する力を取得し、サーボバルブの動的特性は以下のように示され、
Figure 2023517140000009
(τはサーボバルブの動力学モデルの時定数であり、u(t)は電流入力である。)
未知の不感帯があるスプールの変位Γ(x)を考慮すると、その式は
Figure 2023517140000010
であり、
(パラメータmとmとは不感帯特性曲線の左右の傾きを表し、パラメータbとbとは入力非線形のブレークポイントを表す。)
スプールの不感帯はΓ(x)=m(t)x+d(t)で表され、m(t)、d(t)の式は以下のように示され、
Figure 2023517140000011
Figure 2023517140000012
(m、m、b、bが有界であるため、d(t)が有界であり、d(t)の上限を
Figure 2023517140000013
とする。)
状態変数を
Figure 2023517140000014
のように定義し、未知のスプールの不感帯を有する以下の位置追跡システムの状態方程式を取得する、ことを特徴とする。
Figure 2023517140000015
(ただし、
Figure 2023517140000016
。) Further, the process of establishing a model of a single output rod servo-hydraulic system is as follows:
Establishing a hydraulic cylinder dynamics model based on Newton's second law,
Figure 2023517140000002
(where f(t) represents the various interferences, x and m represent the position and mass of the load, respectively, B is the viscous damping coefficient, K is the equivalent spring stiffness of the load, and the load is inertial For the load, K=0 and F=A 1 *P 1 -A 2 *P 2 is the main force output by the hydraulic actuator, where P 1 , P 2 are the large cavity and small is the pressure with the cavity, and A 1 , A 2 are the effective areas of the piston with the large and small cavities.)
Using a 3-position 5-port servo valve, the dynamics of the load pressure is expressed by the following equation,
Figure 2023517140000003
Figure 2023517140000004
(where V 1 and V 2 are the volumes of the rod cavity and rodless cavity respectively, V 1 =V 01 +A 1 x, V 2 =V 02 -A 2 x, V 01 and V 02 is the volume of the rodless cavity and the rod cavity respectively when the piston is in the initial position, C t is the internal leakage coefficient of the hydraulic cylinder, C e is the external leakage coefficient of the hydraulic cylinder, β e is the elastic modulus of the oil, Q1 is the hydraulic fluid flow rate in the rod cavity, and Q2 is the hydraulic fluid flow rate in the rodless cavity.)
however,
Figure 2023517140000005
Figure 2023517140000006
(however,
Figure 2023517140000007
is the flow gain and the expression for s(Γ(x v )) is
Figure 2023517140000008
and
P s is the oil supply pressure of the hydraulic system, P r is the oil return pressure of the hydraulic system, C d is the flow coefficient of the orifice, w is the area gradient of the slide valve, and ρ is the density of the oil. , xv are the spools of the servovalves. )
Controlling the displacement xv of the spool of the servo valve by the voltage or current input u to obtain the corresponding force required, the dynamic characteristics of the servo valve are shown as:
Figure 2023517140000009
(τ is the time constant of the dynamic model of the servovalve and u(t) is the current input.)
Considering the spool displacement Γ(x v ) with an unknown dead band, the expression is
Figure 2023517140000010
and
(The parameters mr and ml represent the left and right slopes of the deadband characteristic curve, and the parameters br and bl represent the breakpoints of the input nonlinearity.)
The dead band of the spool is expressed by Γ( xv )=m(t) xv +d(t), and the formulas for m(t) and d(t) are shown below,
Figure 2023517140000011
Figure 2023517140000012
Since (m r , m l , b r , b l are bounded, d(t) is bounded and the upper bound of d(t) is given by
Figure 2023517140000013
and )
the state variable
Figure 2023517140000014
and obtain the following state equation of the position tracking system with unknown spool deadband:
Figure 2023517140000015
(however,
Figure 2023517140000016
. )

さらに、単一出力ロッドのサーボ油圧システムのモデルに基づいて、低複雑度制御ポリシーを用いて単一出力ロッドのサーボ油圧システムのコントローラを設計する前記プロセスは以下の通り、
誤差変数を定義し、
=x-x,z=x-α,z=F-α,z=x-α(17)
(ただし、α、α、αは仮想制御量であり、x1rは油圧位置追跡システムの指令信号であり、F(F=A-A)はアクチュエータにより出力される主動力であり、zは位置追跡誤差であり、z,i=2…4は制御誤差である。)
4つの正の滑らかな減少関数

Figure 2023517140000017
を指定された性能関数として選択し、正規化誤差ζを定義し、z,i=1…4に対して誤差変換を行い、変換された誤差
Figure 2023517140000018

Figure 2023517140000019
であり、
以下の仮想コントローラを選択し、
Figure 2023517140000020
Figure 2023517140000021
Figure 2023517140000022
サーボ油圧システムのコントローラは、
Figure 2023517140000023
である。 Further, the process of designing a controller for a single output rod servo-hydraulic system using a low-complexity control policy based on a model of a single output rod servo-hydraulic system is as follows:
Define an error variable,
z 1 =x 1 −x r ,z 2 =x 2 −α 1 ,z 3 =F −α 2 ,z 4 =x 5 −α 3 (17)
(where α 1 , α 2 , α 3 are virtual control variables, x 1r is the command signal of the hydraulic position tracking system, and F (F = A 1 x 3 - A 2 x 4 ) is output by the actuator. is the prime mover, z 1 is the position tracking error, and z i , i=2 . . . 4 are the control errors.)
4 positive smooth decreasing functions
Figure 2023517140000017
as the specified performance function, define the normalized errors ζ i , perform an error transformation on z i , i=1 . . .
Figure 2023517140000018
teeth
Figure 2023517140000019
and
Select the virtual controller below,
Figure 2023517140000020
Figure 2023517140000021
Figure 2023517140000022
The controller of the servo hydraulic system is
Figure 2023517140000023
is.

さらに、単一出力ロッドのサーボ油圧システムのコントローラ及び単一出力ロッドのサーボ油圧システムのモデルに基づいて、単一出力ロッドのサーボ油圧システムの安定性を証明する前記プロセスは以下の通り、S3-1:初期値定理を用い、正規化誤差ベクトルが期間t∈[0,τmax)にある場合、空でない開集合Ωζに最大解があることを証明し、S3-2:正規化誤差ベクトルが期間t∈[0,τmax)にあり、空でない開集合Ωζに最大解がある場合、仮想コントローラ及びサーボ油圧システムのコントローラは、t∈[0,τmax)に対して閉ループ信号が有界であることを確保でき、S3-3:初期値の提案を用いて、S3-2におけるτmax=+∞の場合、すべての閉ループ信号が有界であることが依然として正確であることを証明する。 Further, based on the controller of the single output rod servo-hydraulic system and the model of the single output rod servo-hydraulic system, the process of proving the stability of the single output rod servo-hydraulic system is as follows: S3- 1: Using the initial value theorem, prove that a non-empty open set Ω ζ has a maximum solution if the normalized error vector is in the period tε[0, τ max ), S3-2: Normalized error vector is in the period t ∈ [0, τ max ) and there is a maximum solution in the non-empty open set Ω ζ , the virtual controller and the controller of the servo-hydraulic system determine that for t ∈ [0, τ max ) the closed-loop signal is S3-3: Using the initial value proposal, it is still true that all closed-loop signals are bounded if τ max =+∞ in S3-2. Prove.

さらに、初期値定理を用い、正規化誤差ベクトルが期間t∈[0,τmax)にある場合、空でない開集合Ωζに最大解があることを証明する前記ステップにおいては、空でない開集合Ωζは、性能関数ρを選択する場合、

Figure 2023517140000024
を満たすことにより、
Figure 2023517140000025
を得ることができるため、正規化誤差ベクトルはζ(0)∈Ωζであり、所望の追跡奇跡x、性能関数ρ(t),i=1…4、中間制御信号α,i=1…3及び制御率uが連続的で滑らかに導出可能であり、スプールの不感帯の位置追跡システムの状態方程式における動力学変数が連続的に導出可能な関数であるため、式の正規化誤差ベクトルの関数L(t,ζ)は、時間tに関して区分的に連続し、開集合Ωζでは、ζはリプシッツ条件を満たし、初期値定理における条件がいずれも満たされるため、期間[0,τmax)で、正規化誤差ベクトルには唯一の最大解ζ(t)∈Ωζがあり、∀t∈[0,τmax)に対して、いずれもζ(t)∈Ωζ,i=1…4を確保できる。 Further, in the step of using the initial value theorem to prove that the non-empty open set Ω ζ has a maximum solution if the normalized error vector is in the period tε[0, τ max ), the non-empty open set Ω ζ , if we choose the performance function ρ i ,
Figure 2023517140000024
By satisfying
Figure 2023517140000025
, the normalized error vector is ζ(0)∈Ω ζ , the desired tracking miracle x r , the performance function ρ i (t), i=1 . . . 4, the intermediate control signals α i , i = 1 ... 3 and the control rate u are continuously and smoothly derivable, and the dynamics variables in the state equation of the position tracking system of the dead zone of the spool are continuously derivable functions, so the normalization error of the equation The vector function L(t, ζ) is piecewise continuous with respect to time t, and for an open set Ω ζ , ζ satisfies the Lipschitz condition and both conditions in the initial value theorem are satisfied, so the period [0, τ max ), the normalized error vector has a unique maximum solution ζ(t)∈Ω ζ , and for ∀t∈[0, τ max ), any ζ i (t)∈Ω ζ , i= 1...4 can be secured.

さらに、正規化誤差ベクトルが期間t∈[0,τmax)にあり、空でない開集合Ωζに最大解がある場合、仮想コントローラ及びサーボ油圧システムのコントローラは、t∈[0,τmax)に対して閉ループ信号が有界であることを確保できる前記プロセスは以下の通り、
S3-2-1:式(18)に基づいて

Figure 2023517140000026
の時間に対する導関数
Figure 2023517140000027
を求めて、
Figure 2023517140000028
(ただし、
Figure 2023517140000029
であり、1つの正の定数rM1があり、0>r>rM1となる。)
さらに、
Figure 2023517140000030
を得ることができ、
Figure 2023517140000031
とし、
Figure 2023517140000032
であり、
Figure 2023517140000033
が有界であることを考慮すると、極値定理から、aが有界であることが分かり、
(31)と(32)とから、
Figure 2023517140000034
を得て、
リアプノフ関数
Figure 2023517140000035
を選択し、
(33)と(34)とは、不等式補助定理の形態を満たし、不等式補助定理から分かるように、aが有界であり、すなわち∃A>0であり、|a|<A、0<Δ(t)<rM1となるため、
Figure 2023517140000036
であり、
Figure 2023517140000037
が有界であり、
Figure 2023517140000038
を得ることができるため、∀t∈[0,τmax)に対して、追跡誤差zが指定された境界内で減衰され、期間t∈[0,τmax)での追跡誤差zが指定された性能境界内で収束することを確保し、指令が有界正弦信号であるため、状態変数xが有界であり、
S3-2-2:式(18)に基づいて
Figure 2023517140000039
の時間に対する導関数
Figure 2023517140000040
を求めて、
Figure 2023517140000041
(ただし、
Figure 2023517140000042
であり、1つの正の定数rM2があり、0>r>rM2となる。)
さらに、
Figure 2023517140000043
を得ることができ、
Figure 2023517140000044
とし、
Figure 2023517140000045
であり、
Figure 2023517140000046
が有界であり、f(t)が未知の有界な外乱項であり、Bが有界な不確実性パラメータであることを考慮すると、極値定理から、aが有界であることが分かり、
(36)と(37)とから、
Figure 2023517140000047
を得て、
リアプノフ関数
Figure 2023517140000048
を選択し、
(38)と(39)とは、不等式補助定理の形態を満たし、不等式補助定理から分かるように、aが有界であり、すなわち∃A>0であり、|a|<A
Figure 2023517140000049
となるため、
Figure 2023517140000050
であり、
Figure 2023517140000051
が有界であり、
Figure 2023517140000052
を得ることができるため、∀t∈[0,τmax)に対して、制御誤差zが有界であり及び状態変数xが有界であることを確保し、
S3-2-3:式(18)に基づいて
Figure 2023517140000053
の時間に対する導関数
Figure 2023517140000054
を求めて、
Figure 2023517140000055
(ただし、
Figure 2023517140000056
であり、1つの正の定数rM3があり、0>r>rM3となる。)
式(10)、(24)、(11)及び(21)を(40)に代入して、
Figure 2023517140000057
を得て、
Figure 2023517140000058
とし、
Figure 2023517140000059
、ただし、V=V01+Ax、V=V02-Ax。)
サーボ油圧システムアクチュエータの安全な動作に鑑みて、安全マージンが残され、すなわち、物理的構造に基づいて、油圧シリンダが中央位置の付近に最大12センチだけ上下に変動でき、指令信号の振幅が10センチ以下であり、h、h、hが有界であることを確保し、すなわち、それぞれ3つの正数
Figure 2023517140000060
があり、
Figure 2023517140000061
となり、
Figure 2023517140000062
であり、
Figure 2023517140000063
が有界であり、m(t)がスプールの不感帯の未知の傾きであり、m(t)が有界である(すなわち∃M>0,0<m(t)<M)ことを考慮すると、極値定理から、aが有界であることが分かり、
(41)と(42)とから、
Figure 2023517140000064
を得て、
リアプノフ関数
Figure 2023517140000065
を選択し、
(43)と(44)とは、不等式補助定理の形態を満たし、不等式補助定理から分かるように、aが有界であり、すなわち∃A3>0であり、|a3|<A3、
Figure 2023517140000066
となるため、
Figure 2023517140000067
であり、
Figure 2023517140000068
が有界であり、
Figure 2023517140000069
を得ることができるため、∀t∈[0,τmax)に対して、制御誤差zが有界であり及び主動力uが有界であることを確保し、状態方程式の数式から見ると、状態変数xとxとが有界であり、いずれも油圧源Pの圧力値よりも小さく、
S3-2-4:式(18)に基づいて
Figure 2023517140000070
の時間に対する導関数
Figure 2023517140000071
を求めて、
Figure 2023517140000072
(ただし、
Figure 2023517140000073
であり、1つの正の定数rM4があり、0>r>rM4となる。)
さらに、
Figure 2023517140000074
を得て、
Figure 2023517140000075
とし、
Figure 2023517140000076
であり、
Figure 2023517140000077
が有界であり、kがコントローラを設計するときのオプションパラメータであることを考慮すると、aが有界であり、
(46)と(47)とから、
Figure 2023517140000078
を得て、
リアプノフ関数
Figure 2023517140000079
を選択し、
(48)と(49)とは、不等式補助定理の形態を満たし、不等式補助定理から分かるように、aが有界であり、すなわち∃A>0であり、|a|<A
Figure 2023517140000080
となるため、
Figure 2023517140000081
であり、
Figure 2023517140000082
が有界であり、
Figure 2023517140000083
を得ることができるため、∀t∈[0,τmax)に対して、制御誤差zが有界であり及び状態変数xが有界であることを確保する。 Furthermore, if the normalized error vector is in the period t ∈ [0, τ max ) and there is a maximum solution in the non-empty open set Ω ζ , then the virtual controller and the controller of the servo-hydraulic system will determine t ∈ [0, τ max ) The process that can ensure that the closed-loop signal is bounded for is as follows:
S3-2-1: based on formula (18)
Figure 2023517140000026
the derivative with respect to time of
Figure 2023517140000027
looking for
Figure 2023517140000028
(however,
Figure 2023517140000029
, with one positive constant r M1 , such that 0>r 1 >r M1 . )
moreover,
Figure 2023517140000030
can get
Figure 2023517140000031
year,
Figure 2023517140000032
and
Figure 2023517140000033
is bounded, we know from the extremum theorem that a 1 is bounded,
From (31) and (32),
Figure 2023517140000034
to get
Lyapunov function
Figure 2023517140000035
and select
(33) and (34) satisfy the form of the inequality lemma, and as can be seen from the inequality lemma, a 1 is bounded, i.e. ∃A 1 >0 and |a 1 |<A 1 , 0<Δ 1 (t)<r M1 k 1 , so
Figure 2023517140000036
and
Figure 2023517140000037
is bounded, and
Figure 2023517140000038
so that for ∀tε[0, τ max ), the tracking error z 1 decays within the specified bounds, and the tracking error z 1 over time period tε[0, τ max ) is ensuring convergence within the specified performance bounds, the state variable x1 is bounded because the command is a bounded sine signal,
S3-2-2: based on formula (18)
Figure 2023517140000039
the derivative with respect to time of
Figure 2023517140000040
looking for
Figure 2023517140000041
(however,
Figure 2023517140000042
and there is one positive constant r M2 , such that 0>r 2 >r M2 . )
moreover,
Figure 2023517140000043
can get
Figure 2023517140000044
year,
Figure 2023517140000045
and
Figure 2023517140000046
is bounded, f(t) is an unknown bounded disturbance term, and B is a bounded uncertainty parameter, then from the extremum theorem it follows that a 2 is bounded is understood,
From (36) and (37),
Figure 2023517140000047
to get
Lyapunov function
Figure 2023517140000048
and select
(38) and (39) satisfy the form of the inequality lemma, and as can be seen from the inequality lemma, a 2 is bounded, i.e. ∃A 1 >0 and |a 2 |<A 2 ,
Figure 2023517140000049
Because
Figure 2023517140000050
and
Figure 2023517140000051
is bounded, and
Figure 2023517140000052
ensuring that the control error z 2 is bounded and the state variable x 2 is bounded for ∀tε[0, τ max ),
S3-2-3: based on formula (18)
Figure 2023517140000053
the derivative with respect to time of
Figure 2023517140000054
looking for
Figure 2023517140000055
(however,
Figure 2023517140000056
and there is one positive constant r M3 such that 0>r 3 >r M3 . )
Substituting equations (10), (24), (11) and (21) into (40),
Figure 2023517140000057
to get
Figure 2023517140000058
year,
Figure 2023517140000059
, where V 1 =V 01 +A 1 x, V 2 =V 02 −A 2 x. )
In view of the safe operation of the servo-hydraulic system actuators, a safety margin is left, i.e., based on the physical construction, the hydraulic cylinder can move up and down about the center position by a maximum of 12 cm, and the command signal amplitude is 10 cm. centimeters or less and ensure that h 1 , h 2 , h 3 are bounded, i.e. each three positive numbers
Figure 2023517140000060
there is
Figure 2023517140000061
becomes,
Figure 2023517140000062
and
Figure 2023517140000063
is bounded, m(t) is the unknown slope of the dead zone of the spool, and m(t) is bounded (i.e. ∃M r >0,0<m(t)<M r ). Considering, from the extremum theorem, we know that a 3 is bounded,
From (41) and (42),
Figure 2023517140000064
to get
Lyapunov function
Figure 2023517140000065
and select
(43) and (44) satisfy the form of the inequality lemma, and as can be seen from the inequality lemma, a3 is bounded, i.e. ∃A3>0, |a3|<A3,
Figure 2023517140000066
Because
Figure 2023517140000067
and
Figure 2023517140000068
is bounded, and
Figure 2023517140000069
so that for ∀t∈[0, τ max ), we ensure that the control error z 3 is bounded and the main power u s is bounded, and from the equation of state formula we see , the state variables x3 and x4 are bounded, and both are smaller than the pressure value of the hydraulic source Ps ,
S3-2-4: based on formula (18)
Figure 2023517140000070
the derivative with respect to time of
Figure 2023517140000071
looking for
Figure 2023517140000072
(however,
Figure 2023517140000073
and there is one positive constant r M4 such that 0>r 4 >r M4 . )
moreover,
Figure 2023517140000074
to get
Figure 2023517140000075
year,
Figure 2023517140000076
and
Figure 2023517140000077
is bounded and k 4 is an optional parameter when designing the controller, a 4 is bounded and
From (46) and (47),
Figure 2023517140000078
to get
Lyapunov function
Figure 2023517140000079
and select
(48) and (49) satisfy the form of the inequality lemma, and as can be seen from the inequality lemma, a 4 is bounded, i.e. ∃A 4 >0 and |a 4 |<A 4 ,
Figure 2023517140000080
Because
Figure 2023517140000081
and
Figure 2023517140000082
is bounded, and
Figure 2023517140000083
to ensure that the control error z 4 is bounded and the state variable x 5 is bounded for ∀tε[0, τ max ).

さらに、前記不等式補助定理の式は、

Figure 2023517140000084
であり、
aが有界であり(すなわち∃A>0,|a|<A)、
Figure 2023517140000085
である場合、xが必ず有界であり、
Figure 2023517140000086
である。 Furthermore, the expression of the above inequality lemma is
Figure 2023517140000084
and
a is bounded (i.e. ∃A>0, |a|<A),
Figure 2023517140000085
, then x must be bounded, and
Figure 2023517140000086
is.

さらに、初期値の提案を用いて、S3-2におけるτmax=+∞の場合、すべての閉ループ信号が有界であることが依然として正確であることを証明する前記プロセスは以下の通り、

Figure 2023517140000087
であるため、
Figure 2023517140000088
は、空でない集合であり、
Figure 2023517140000089
であり、τmax<+∞を仮定すると、初期値の提案では1つの時定数t∈[0,τmax)があることが強調され、
Figure 2023517140000090
となり、矛盾を引き起こし、得られた仮定がτmax=+∞であり、いずれかのt∈[0,+∞)に対して、すべての閉ループシステム信号が有界であり、追跡誤差が指定された境界内で収束することを確保できる。 Furthermore, the above process to prove that it is still correct that all closed-loop signals are bounded for τ max =+∞ in S3-2 using the initial value proposal is as follows:
Figure 2023517140000087
Because
Figure 2023517140000088
is a non-empty set and
Figure 2023517140000089
, and assuming τ max <+∞, the initial value proposal emphasizes that there is one time constant t 1 ε[0, τ max ),
Figure 2023517140000090
, which causes a contradiction, and the assumption obtained is that τ max =+∞, that for any tε[0,+∞) all closed-loop system signals are bounded, and the tracking error is specified can be guaranteed to converge within the bounds

上記技術的解決手段が用いられるため、本発明に係る非対称サーボ油圧位置追跡システムのための低複雑度制御方法は、油圧システムに存在する様々な不確実性問題(例えば、未知の摩擦効果、パラメータの不確実性及び負荷変化)及び未知の非線形問題(例えば、スプールの不感帯、外乱)を解決でき、コントローラの設計は、正確な数学的モデルに依存せず、測定可能な状態信号のみが必要であり、制御率の計算は、バックステッピング適応に基づいて開発された既存のアルゴリズムと比較して、計算プロセスが簡単であり、計算量が少なく、リアルタイム制御が容易であり、工学的実現が容易であり、本発明は、追跡誤差の収束速度と定常状態の精度とを確保でき、最後に、実験結果によると、従来のpid制御方法と比較して、本発明の位置追跡効果は、定常状態の精度が高くなり、追跡変位の位相遅れ程度が小さくなることが証明されている。所望の軌跡信号の周波数が向上する場合、pidアルゴリズムの制御により得られた実際の追跡変位の位相遅れ程度が大きいほど、振幅が大幅に減衰され、変位追跡誤差が大きくなるが、本発明の制御アルゴリズムは、所望の軌跡信号の周波数が向上する場合、得られた実際の制御変位の遅れ程度が基本的に変化せず、追跡誤差が常に指定された境界内で収束し、振幅が基本的に減衰されない。所望の信号の周波数が変化せず、振幅が大きくなる場合、pidアルゴリズムの変位追跡誤差が大きくなるが、本発明の制御アルゴリズムは、依然として追跡誤差が指定された境界内で収束することを確保でき、位相遅れ程度が小さく、振幅が基本的に減衰されない。 Because the above technical solution is used, the low-complexity control method for the asymmetric servo-hydraulic position tracking system according to the present invention can be applied to various uncertainties existing in the hydraulic system (e.g. unknown friction effect, parameter uncertainty and load changes) and unknown non-linear problems (e.g. spool deadband, disturbances), the controller design does not rely on exact mathematical models and only requires measurable state signals. , and the control rate calculation has a simple calculation process, a small amount of calculation, an easy real-time control, and an easy engineering realization, compared with existing algorithms developed based on backstepping adaptation. There is, the present invention can ensure the tracking error convergence speed and steady state accuracy. Finally, according to the experimental results, compared with the conventional pid control method, the position tracking effect of the present invention is It has been proven that the accuracy is higher and the phase lag of the tracking displacement is smaller. When the frequency of the desired trajectory signal is improved, the greater the degree of phase lag of the actual tracking displacement obtained by the control of the pid algorithm, the greater the attenuation of the amplitude and the greater the displacement tracking error. The algorithm ensures that when the frequency of the desired trajectory signal increases, the lag degree of the obtained actual control displacement remains essentially unchanged, the tracking error always converges within the specified bounds, and the amplitude is essentially not attenuated. If the desired signal's frequency does not change and its amplitude increases, the displacement tracking error of the pid algorithm will increase, but the control algorithm of the present invention can still ensure that the tracking error converges within the specified bounds. , the phase lag is small and the amplitude is basically not attenuated.

本願の実施例又は従来技術における技術的解決手段をより明確に説明するために、以下、実施例又は従来技術の説明で使用される必要がある図面を簡単に説明し、明らかに、以下に説明される図面は、本願に記載のいくつかの実施例に過ぎず、当業者であれば、創造的な労働を必要とせずに、これらの図面に基づいて他の図面を取得することができる。 In order to describe the embodiments of the present application or the technical solutions in the prior art more clearly, the following briefly describes the drawings that need to be used in the description of the embodiments or the prior art. The drawings shown are only some examples described in this application, and those skilled in the art can obtain other drawings based on these drawings without creative efforts.

本発明に係る非対称サーボ油圧位置追跡システムのための低複雑度制御方法のフローチャートである。1 is a flowchart of a low-complexity control method for an asymmetric servohydraulic position tracking system according to the present invention; 単一出力ロッドのサーボ油圧システムのモデルのブロック図である。1 is a block diagram of a model of a single output rod servo-hydraulic system; FIG. 実験プラットフォームの構造構成図Iである。1 is a structural diagram I of an experimental platform; FIG. 実験プラットフォームの構造構成図IIである。Fig. 2 is a structural diagram II of the experimental platform; 実験プラットフォームの構造構成図IIIである。Fig. III is the structural configuration diagram of the experimental platform; 実験プラットフォームの構造構成図IVである。FIG. IV is a structural configuration diagram of the experimental platform; 本発明のコントローラの作用下での誤差収束のシミュレーション曲線図である。FIG. 4 is a simulation curve diagram of error convergence under the action of the controller of the present invention; 外乱作用下での追跡誤差収束のシミュレーション曲線図である。FIG. 4 is a simulation curve diagram of tracking error convergence under disturbance action; 未知のスプールの不感帯の作用下での誤差収束のシミュレーション曲線図である。FIG. 10 is a simulation curve diagram of error convergence under the action of unknown spool deadband; 未知のスプールの不感帯の作用下でのスプール変位のシミュレーション曲線図である。FIG. 11 is a simulation curve diagram of spool displacement under the action of an unknown spool dead band; 本発明の制御方法と、指定された性能を満たす、スプールの未知の不感帯を有する適応バックステッピング制御方式(SPPFBSA)と、の誤差収束の比較シミュレーション曲線図である。FIG. 10 is a comparative simulation curve diagram of the error convergence of the control method of the present invention and an adaptive backstepping control scheme with unknown spool deadband (SPPFBSA) meeting specified performance; 本発明の制御方法とSPPFBSAとのシミュレーション計算時間及び誤差調整時間の統計図である。FIG. 4 is a statistical diagram of simulation calculation time and error adjustment time between the control method of the present invention and SPPFBSA; 指令信号がTS=3s:x=100sin((2*pi*t)/3)mmであるときの、本発明の制御方法とpidとの変位追跡の比較実験図である。FIG. 10 is a comparative experiment diagram of displacement tracking between the control method of the present invention and pid when the command signal is TS=3s: xr =100sin((2*pi*t)/3)mm; 指令信号がTS=3s:x=100sin((2*pi*t)/3)mmであるときの、本発明の制御方法の追跡誤差の実験図である。FIG. 10 is an experimental diagram of the tracking error of the control method of the present invention when the command signal is TS=3s: xr =100sin((2*pi*t)/3)mm; 指令信号がTS=2s:x=100sin(pi*t)mmであるときの、本発明の制御方法とpidとの変位追跡の比較実験図である。FIG. 10 is a comparison experimental diagram of displacement tracking between the control method of the present invention and pid when the command signal is TS=2s: xr =100sin(pi*t)mm; 指令信号がTS=2s:x=100sin(pi*t)mmであるときの、本発明の制御方法の追跡誤差の実験図である。FIG. 10 is an experimental diagram of the tracking error of the control method of the present invention when the command signal is TS=2s: xr =100sin(pi*t)mm; 指令信号がTS=1s:x=50sin(2*pi*t)mmであるときの、本発明の制御方法とpidとの変位追跡の比較実験図である。FIG. 10 is a comparison experimental diagram of displacement tracking between the control method of the present invention and pid when the command signal is TS=1 s:x r =50 sin(2*pi*t) mm; 指令信号がTS=1s:x=50sin(2*pi*t)mmであるときの、本発明の制御方法の追跡誤差の実験図である。FIG. 5 is an experimental diagram of the tracking error of the control method of the present invention when the command signal is TS=1 s:x r =50 sin(2*pi*t) mm; 指令信号がTS=2s:x=100sin(pi*t)mmであるときの、本発明の制御方法とpidとの変位追跡の比較実験図である。FIG. 10 is a comparison experimental diagram of displacement tracking between the control method of the present invention and pid when the command signal is TS=2s: xr =100sin(pi*t)mm; 指令信号がTS=2s:x=100sin(pi*t)mmであるときの、本発明の制御方法の追跡誤差の実験図である。FIG. 10 is an experimental diagram of the tracking error of the control method of the present invention when the command signal is TS=2s: xr =100sin(pi*t)mm; 指令信号がTS=2s:x=80sin(pi*t)mmであるときの、本発明の制御方法とpidとの変位追跡の比較実験図である。FIG. 10 is a comparison experimental diagram of displacement tracking between the control method of the present invention and pid when the command signal is TS=2s: xr =80sin(pi*t)mm; 指令信号がTS=2s:x=80sin(pi*t)mmであるときの、本発明の制御方法の追跡誤差の実験図である。FIG. 10 is an experimental diagram of the tracking error of the control method of the present invention when the command signal is TS=2s: xr =80sin(pi*t)mm; 指令信号がTS=2s:x=60sin(pi*t)mmであるときの、本発明の制御方法とpidとの変位追跡の比較実験図である。FIG. 10 is a comparison experiment diagram of displacement tracking between the control method of the present invention and pid when the command signal is TS=2s: xr =60sin(pi*t)mm; 指令信号がTS=2s:x=60sin(pi*t)mmであるときの、本発明の制御方法の追跡誤差の実験図である。FIG. 10 is an experimental diagram of the tracking error of the control method of the present invention when the command signal is TS=2s: xr =60sin(pi*t)mm;

本発明の技術的解決手段及び利点をより明確にするために、以下、本発明の実施例における図面を参照しながら、本発明の実施例における技術的解決手段を明確、完全に説明する。 In order to make the technical solutions and advantages of the present invention clearer, the following clearly and completely describes the technical solutions in the embodiments of the present invention with reference to the drawings in the embodiments of the present invention.

以下、図面を参照しながら、本発明の具体的な実施形態をさらに説明する。 Specific embodiments of the present invention will be further described below with reference to the drawings.

図1は本発明に係る非対称サーボ油圧位置追跡システムのための低複雑度制御方法のフローチャートであり、非対称サーボ油圧位置追跡システムのための低複雑度制御方法は、単一出力ロッドのサーボ油圧システムのモデルを確立するステップS1と、単一出力ロッドのサーボ油圧システムのモデルに基づいて、低複雑度制御ポリシーを用いて単一出力ロッドのサーボ油圧システムのコントローラを設計するステップS2と、単一出力ロッドのサーボ油圧システムのコントローラ及び単一出力ロッドのサーボ油圧システムのモデルに基づいて、単一出力ロッドのサーボ油圧システムの安定性を証明するステップS3と、を含み、図2は単一出力ロッドのサーボ油圧システムのモデルのブロック図である。 FIG. 1 is a flowchart of a low-complexity control method for an asymmetric servohydraulic position tracking system according to the present invention. designing a controller for the single output rod servo-hydraulic system using a low-complexity control policy based on the model for the single output rod servo-hydraulic system; and step S3 of proving the stability of the single output rod servo-hydraulic system based on the controller of the output rod servo-hydraulic system and the model of the single output rod servo-hydraulic system, wherein FIG. 1 is a block diagram of a model of a rod's servo-hydraulic system; FIG.

単一出力ロッドのサーボ油圧システムのモデルを確立する前記プロセスは以下の通り、
先ず、ニュートンの第2法則に基づいて油圧シリンダの動力学モデルを確立し、

Figure 2023517140000091
(力学方程式(1)の追加項f(t)は様々な干渉(例えば、モデル化されていない摩擦効果、モデル化されていない動力学、外乱など)を表し、xとmとはそれぞれ負荷の位置と質量とを表し、Bは粘性減衰係数であり、Kは負荷の等価ばね剛性であり、負荷が主に慣性負荷であるため、K=0であり、F=A*P-A*Pは油圧アクチュエータにより出力される主動力であり、ただし、P、Pは油圧シリンダの大きなキャビティと小さなキャビティとの圧力であり、圧力センサにより測定でき、A、Aは大きなキャビティと小さなキャビティとのピストンの有効面積である。)
仮定1:干渉項f(t)が有界であり、すなわち∃f>0であり、|f(t)<f|となり、粘性減衰係数Bが不確実な有界な正のパラメータであり、すなわち
Figure 2023517140000092
であり、
Figure 2023517140000093
となり、本発明では、3ポジション5ポートバルブのピストンアクチュエータを用い、負荷圧力の動力学は以下のように示され、
Figure 2023517140000094
Figure 2023517140000095
(式中、V、Vはそれぞれロッドキャビティとロッドレスキャビティとの容積(V=V01+Ax、V=V02-Ax)であり、V01とV02とはそれぞれピストンが初期位置にあるときのロッドレスキャビティとロッドキャビティとの容積であり、Cは油圧シリンダの内部漏れ係数でありであり、Cは油圧シリンダの外部漏れ係数であり、βはオイルの弾性率であり、Qはロッドキャビティの給油(油戻し)流量であり、Qはロッドレスキャビティの油戻し(給油)流量である。)、QとQとの計算式は、
Figure 2023517140000096
Figure 2023517140000097
であり、
(ただし、
Figure 2023517140000098
は流量ゲインであり、s(Γ(x))の式は
Figure 2023517140000099
であり、
は油圧システムの給油圧力であり、Pは油圧システムの油戻し圧力であり、Cはオリフィスの流量係数であり、wはスライドバルブの面積勾配であり、ρはオイルの密度であり、xはサーボバルブのスプールの変位である。)
仮定2:不確実性パラメータの内部漏れ係数Cが有界であり、すなわち
Figure 2023517140000100
であり、
Figure 2023517140000101
となり、不確実性パラメータの内部漏れ係数Cが有界であり、すなわち
Figure 2023517140000102
であり、
Figure 2023517140000103
となり、不確実性パラメータのオイルの弾性率βが有界であり、すなわち
Figure 2023517140000104
であり、
Figure 2023517140000105
となり、
油圧アクチュエータの操作では、実際には、電圧又は電流入力uによりサーボバルブのスプールの変位xを制御して、必要な対応する力をさらに取得し、サーボバルブの動的特性は以下のように示され、
Figure 2023517140000106
(τはサーボバルブの動力学モデルの時定数であり、u(t)は電流入力である。)
本明細書では、未知の不感帯があるスプールの変位Γ(x)を考慮すると、その式は
Figure 2023517140000107
であり、
(パラメータmとmとは不感帯特性曲線の左右の傾きを表し、パラメータbとbとは入力非線形のブレークポイントを表す。)
仮定3:パラメータm、m、b、bが未知の正数であり、m=mであり、本明細書では、傾きの上限はmmaxであり、bの上限は
Figure 2023517140000108
であり、bの上限は
Figure 2023517140000109
であり、スプールの不感帯はΓ(x)=m(t)x+d(t)で表されてもよく、m(t)、d(t)の式は以下のように示され、
Figure 2023517140000110
Figure 2023517140000111
(m、m、b、bが有界であるため、d(t)が有界であり、d(t)の上限を
Figure 2023517140000112
とする。)
周知のように、サーボバルブは、電気油圧アクチュエータにおける重要な機械部材であり、サーボバルブのスプールの変位を電流又は電圧制御し、さらにオイルキャビティへの作動油の抽入又は抽出を制御し、最後に、アクチュエータは対応する移動を実行し、明らかに、サーボ油圧位置追跡システムでは、スプールの不感帯の非線形問題が必然的に存在し、従って、この問題の悪影響を考慮してより良好なシステム性能を取得する必要があり、また、実際の応用では、不感帯モデルの正確な傾き及び間隔ポイントを取得することは困難であることを考慮するため、この問題を解決するためのロバスト性の高い新規の制御ポリシーが提案されており、
油圧非対称シリンダの位置追跡システムを状態空間式の形態に導出するために、状態変数を
Figure 2023517140000113
のように定義し、(1)~(6)の動力学方程式によりスプールの不感帯を有する以下の状態方程式を取得できる。
Figure 2023517140000114
(ただし、
Figure 2023517140000115
。) The process of establishing a model of a single output rod servo-hydraulic system is as follows:
First, establish a dynamics model of the hydraulic cylinder based on Newton's second law,
Figure 2023517140000091
(The additional term f(t) in the dynamics equation (1) represents various interferences (e.g., unmodeled frictional effects, unmodeled dynamics, disturbances, etc.), and x and m are the loads, respectively. Denoting position and mass, B is the viscous damping coefficient, K is the equivalent spring stiffness of the load, and since the load is primarily an inertial load, K=0 and F=A 1 *P 1 −A 2 * P2 is the main power output by the hydraulic actuator, where P1 , P2 are the pressures of the large and small cavities of the hydraulic cylinder, which can be measured by pressure sensors, and A1 , A2 are is the effective area of the piston between the large and small cavities.)
Assumption 1: The interference term f(t) is bounded, i.e. ∃f m >0, such that |f(t)<f m |, and the viscous damping coefficient B is an uncertain bounded positive parameter Yes, i.e.
Figure 2023517140000092
and
Figure 2023517140000093
Therefore, in the present invention, a 3-position 5-port valve piston actuator is used, and the dynamics of the load pressure is shown as follows,
Figure 2023517140000094
Figure 2023517140000095
(In the formula, V 1 and V 2 are the volumes of the rod cavity and the rodless cavity (V 1 =V 01 +A 1 x, V 2 =V 02 −A 2 x), respectively, and V 01 and V 02 are are the volumes of the rodless cavity and rod cavity respectively when the piston is in the initial position, Ct is the internal leakage coefficient of the hydraulic cylinder, Ce is the external leakage coefficient of the hydraulic cylinder, and βe is is the elastic modulus of oil, Q1 is the rod cavity oil supply (oil return) flow rate, and Q2 is the rodless cavity oil return (oil return) flow rate.), the formula for Q1 and Q2 is ,
Figure 2023517140000096
Figure 2023517140000097
and
(however,
Figure 2023517140000098
is the flow gain and the expression for s(Γ(x v )) is
Figure 2023517140000099
and
P s is the oil supply pressure of the hydraulic system, P r is the oil return pressure of the hydraulic system, C d is the flow coefficient of the orifice, w is the area gradient of the slide valve, and ρ is the density of the oil. , xv is the displacement of the spool of the servovalve. )
Assumption 2: The internal leakage coefficient Ct of the uncertainty parameter is bounded, i.e.
Figure 2023517140000100
and
Figure 2023517140000101
and the internal leakage coefficient C e of the uncertainty parameter is bounded, i.e.
Figure 2023517140000102
and
Figure 2023517140000103
and the uncertainty parameter β e of the elastic modulus of oil is bounded, i.e.
Figure 2023517140000104
and
Figure 2023517140000105
becomes,
In the operation of the hydraulic actuator, in practice, the displacement x v of the spool of the servo valve is controlled by the voltage or current input u to further obtain the required corresponding force, and the dynamic characteristics of the servo valve are as follows: shown,
Figure 2023517140000106
(τ is the time constant of the dynamic model of the servovalve and u(t) is the current input.)
Here, considering the spool displacement Γ(x v ) with an unknown dead band, the expression is
Figure 2023517140000107
and
(The parameters mr and ml represent the left and right slopes of the deadband characteristic curve, and the parameters br and bl represent the breakpoints of the input nonlinearity.)
Assumption 3: The parameters m r , m l , b r , b l are unknown positive numbers, m r = m l , and here the upper bound of the slope is m max and the upper bound of b r is
Figure 2023517140000108
and the upper limit of b l is
Figure 2023517140000109
, and the dead band of the spool may be expressed as Γ(x v )=m(t)x v +d(t), and the formulas for m(t), d(t) are shown as follows,
Figure 2023517140000110
Figure 2023517140000111
Since (m r , m l , b r , b l are bounded, d(t) is bounded and the upper bound of d(t) is given by
Figure 2023517140000112
and )
As is well known, a servovalve is an important mechanical member in an electrohydraulic actuator, which controls current or voltage control of the displacement of the spool of the servovalve, further controls the inflow or extraction of hydraulic oil into the oil cavity, and finally On the other hand, the actuator performs a corresponding movement, and obviously, in a servo hydraulic position tracking system, the non-linear problem of the spool dead band inevitably exists, so considering the adverse effects of this problem leads to better system performance. Considering that it is necessary to obtain and that in practical applications it is difficult to obtain the exact slope and interval points of the dead zone model, a novel robust control is proposed to solve this problem. A policy is proposed
In order to derive the position tracking system for a hydraulic asymmetric cylinder in the form of a state-space formula, the state variables are
Figure 2023517140000113
, and the following equation of state with the dead zone of the spool can be obtained from the dynamic equations (1) to (6).
Figure 2023517140000114
(however,
Figure 2023517140000115
. )

単一出力ロッドのサーボ油圧システムのモデルに基づいて、低複雑度制御ポリシーを用いて単一出力ロッドのサーボ油圧システムのコントローラを設計する前記プロセスは以下の通り、
誤差変数を定義し、
=x-x,z=x-α,z=F-α,z=x-α(17)
(ただし、α、α、αは後続の証明プロセスで取得された仮想制御量であり、X1rは油圧位置追跡システムの指令信号であり、F(F=A-A)はアクチュエータにより出力される主動力であり、zは位置追跡誤差であり、z,i=2…4は制御誤差である。)
追跡誤差の過渡と定常状態との問題を解決し、仮想制御量をスムーズに選択するために、4つの正の滑らかな減少関数

Figure 2023517140000116
を指定された性能関数として選択し、正規化誤差ζを定義し、z,i=2…4に対して誤差変換を行い、変換された誤差
Figure 2023517140000117

Figure 2023517140000118
であり、
以下の仮想制御関数を選択し、
Figure 2023517140000119
Figure 2023517140000120
Figure 2023517140000121
サーボ油圧システムのコントローラは、
Figure 2023517140000122
である。 The process of designing a controller for a single output rod servo-hydraulic system using a low-complexity control policy based on a model of a single output rod servo-hydraulic system is as follows:
Define an error variable,
z 1 =x 1 −x r ,z 2 =x 2 −α 1 ,z 3 =F −α 2 ,z 4 =x 5 −α 3 (17)
(where α 1 , α 2 , α 3 are virtual control variables obtained in the subsequent proof process, X 1r is the command signal of the hydraulic position tracking system, F(F=A 1 x 3 −A 2 x 4 ) is the prime force output by the actuator, z 1 is the position tracking error, and z i , i=2...4 is the control error.)
To solve the tracking error transient and steady-state problem and select the virtual controlled variable smoothly, four positive smooth decreasing functions
Figure 2023517140000116
as the specified performance function, define the normalized errors ζ i , perform an error transformation on z i , i=2 .
Figure 2023517140000117
teeth
Figure 2023517140000118
and
Select the virtual control function below,
Figure 2023517140000119
Figure 2023517140000120
Figure 2023517140000121
The controller of the servo hydraulic system is
Figure 2023517140000122
is.

さらに、単一出力ロッドのサーボ油圧システムのコントローラ及び単一出力ロッドのサーボ油圧システムのモデルに基づいて、単一出力ロッドのサーボ油圧システムの安定性を証明する前記プロセスは以下の通り、S3-1:初期値定理を用い、正規化誤差ベクトルが期間t∈[0,τmax)にある場合、空でない開集合Ωζに最大解があることを証明し、S3-2:正規化誤差ベクトルが期間t∈[0,τmax)にあり、空でない開集合Ωζに最大解がある場合、仮想コントローラ及びサーボ油圧システムのコントローラは、t∈[0,τmax)に対して閉ループ信号が有界であることを確保でき、S3-3:初期値の提案を用いて、S3-2におけるτmax=+∞の場合、すべての閉ループ信号が有界であることが依然として正確であることを証明する。コントローラの安定性の証明をスムーズに完了するために、補助定理が提案されており、該補助定理の正確さが証明される。 Further, based on the controller of the single output rod servo-hydraulic system and the model of the single output rod servo-hydraulic system, the process of proving the stability of the single output rod servo-hydraulic system is as follows: S3- 1: Using the initial value theorem, prove that a non-empty open set Ω ζ has a maximum solution if the normalized error vector is in the period tε[0, τ max ), S3-2: Normalized error vector is in the period t ∈ [0, τ max ) and there is a maximum solution in the non-empty open set Ω ζ , the virtual controller and the controller of the servo-hydraulic system determine that for t ∈ [0, τ max ) the closed-loop signal is S3-3: Using the initial value proposal, it is still true that all closed-loop signals are bounded if τ max =+∞ in S3-2. Prove. In order to smoothly complete the proof of controller stability, a lemma is proposed and the correctness of the lemma is proved.

前記不等式補助定理の式は、

Figure 2023517140000123
であり、
aが有界であり(すなわち∃A>0,|a|<A)、
Figure 2023517140000124
である場合、xが必ず有界であり、
Figure 2023517140000125
であり、
Figure 2023517140000126
であるため、
Figure 2023517140000127
であり、
|x|>Aの場合、
Figure 2023517140000128
であるため、
Figure 2023517140000129
であることが証明され、
位置誤差、及び誤差変数(17)と正規化誤差ζ(18)との関係に基づいて、x、x、F、xを、x=ζρ+x,x=ζρ+α,F=ζρ+α,x=ζρ+α(24)のように書き直すことができ、
(18)におけるζ(i=1…4)を直接導出して、
Figure 2023517140000130
Figure 2023517140000131
Figure 2023517140000132
(ただし、
Figure 2023517140000133

Figure 2023517140000134
を得ることができ、
正規化誤差ベクトルζ=[ζ,ζ,ζ,ζを定義することにより、(25)~(28)を以下の形態に書くことができ、
Figure 2023517140000135
開集合
Figure 2023517140000136
を定義する。 The formula for the inequality lemma is
Figure 2023517140000123
and
a is bounded (i.e. ∃A>0, |a|<A),
Figure 2023517140000124
, then x must be bounded, and
Figure 2023517140000125
and
Figure 2023517140000126
Because
Figure 2023517140000127
and
If |x|>A, then
Figure 2023517140000128
Because
Figure 2023517140000129
is proved to be
Based on the position error and the relationship between the error variables (17) and the normalized error ζ i (18), let x 1 , x 2 , F, x 5 be x 1 = ζ 1 ρ 1 + x r , x 2 = ζ 2 ρ 2 + α 1 , F = ζ 3 ρ 3 + α 2 , x 5 = ζ 4 ρ 4 + α 3 (24) can be rewritten as
Directly deriving ζ i (i=1 . . . 4) in (18),
Figure 2023517140000130
Figure 2023517140000131
Figure 2023517140000132
(however,
Figure 2023517140000133
)
Figure 2023517140000134
can get
By defining the normalized error vector ζ=[ζ 1 , ζ 2 , ζ 3 , ζ 4 ] T , (25)-(28) can be written in the form
Figure 2023517140000135
open set
Figure 2023517140000136
Define

さらに、初期値定理を用い、正規化誤差ベクトルが期間t∈にある場合、空でない開集合Ωζに最大解があることを証明するプロセスは以下の通り、Ωζは空でない開集合であり、性能関数ρを選択する場合、

Figure 2023517140000137
を満たすことにより、
Figure 2023517140000138
を導出することができるため、式(29)に対してζ(0)∈Ωζであり、また、所望の追跡奇跡x、性能関数ρ(t),i=1…4、中間制御信号α,i=1…3及び制御率uが連続的で滑らかに導出可能であり、スプールの不感帯の位置追跡システムの状態方程式における動力学変数が連続的に導出可能な関数であるため、式の正規化誤差ベクトルの関数L(t,ζ)は、時間tに関して区分的に連続し、開集合Ωζでは、ζはリプシッツ条件を満たし、初期値定理における条件がいずれも満たされるため、期間[0,τmax)で、正規化誤差ベクトルには唯一の最大解ζ(t)∈Ωζがあり、∀t∈[0,τmax)に対して、いずれもζ(t)∈Ωζ,i=1…4を確保できる。 Furthermore, using the initial value theorem, the process of proving that a non-empty open set Ω ζ has a maximum solution if the normalized error vector is in the period tε is as follows, Ω ζ is a non-empty open set and , if we choose the performance function ρ i ,
Figure 2023517140000137
By satisfying
Figure 2023517140000138
(29) so that ζ(0)∈Ω ζ for Eq. (29), and the desired tracking miracle x r , the performance function ρ i (t), i=1 . . . 4, the intermediate control Since the signals α i , i=1 . Since the function L(t, ζ) of the normalized error vector in Eq. is piecewise continuous with respect to time t, and for an open set Ω ζ , ζ satisfies the Lipschitz condition and both conditions in the initial value theorem are satisfied, In the period [0, τ max ), the normalized error vector has a unique maximum solution ζ(t)∈Ω ζ , and for ∀t∈[0, τ max ), any ζ i (t)∈ Ω ζ , i=1 . . . 4 can be secured.

さらに、正規化誤差ベクトルが期間t∈[0,τmax)にあり、空でない開集合Ωζに最大解がある場合、仮想コントローラ及びサーボ油圧システムのコントローラは、t∈[0,τmax)に対して閉ループ信号が有界であることを確保できる前記プロセスは以下の通り、
S3-2-1:式(18)に基づいて

Figure 2023517140000139
の時間に対する導関数
Figure 2023517140000140
を求めて、
Figure 2023517140000141
(ただし、
Figure 2023517140000142
であり、1つの正の定数rM1があり、0>r>rM1となる。)
式(10)、(24)及び(19)を(30)に代入して、
Figure 2023517140000143
を得て、
Figure 2023517140000144
とし、
Figure 2023517140000145
であり、
Figure 2023517140000146
が有界であることを考慮すると、極値定理から、aが有界であることが分かり、
(31)と(32)とから、
Figure 2023517140000147
を得て、
リアプノフ関数
Figure 2023517140000148
を選択し、
(33)と(34)とは、不等式補助定理の形態を満たし、不等式補助定理から分かるように、aが有界∃A>0であり、|a|<A、0<Δ(t)<rM1となるため、
Figure 2023517140000149
であり、
Figure 2023517140000150
が有界であり、(13)と(14)とから
Figure 2023517140000151
を得ることができるため、∀t∈[0,τmax)に対して、追跡誤差zが指定された境界内で減衰され、従って、期間t∈[0,τmax)での追跡誤差zが指定された性能境界内で収束することを確保でき、指令が有界正弦信号であるため、状態変数xが有界であり、
S3-2-2:式(18)に基づいて
Figure 2023517140000152
の時間に対する導関数
Figure 2023517140000153
を求めて、
Figure 2023517140000154
(ただし、
Figure 2023517140000155
であり、1つの正の定数rM2があり、0>r>rM2となる。)
式(10)、(24)及び(20)を(35)に代入して、
Figure 2023517140000156
を得て、
Figure 2023517140000157
とし、
Figure 2023517140000158
であり、
Figure 2023517140000159
が有界であり、f(t)が未知の有界な外乱項であり、Bが有界な不確実性パラメータであることを考慮すると、極値定理から、aが有界であることが分かり、
(36)と(37)とから、
Figure 2023517140000160
を得て、
リアプノフ関数
Figure 2023517140000161
を選択し、
(38)と(39)とは、不等式補助定理の形態を満たし、不等式補助定理から分かるように、aが有界∃A>0であり、|a|<A
Figure 2023517140000162
となるため、
Figure 2023517140000163
であり、
Figure 2023517140000164
が有界であり、(13)と(14)とから
Figure 2023517140000165
を得ることができるため、∀t∈[0,τmax)に対して、制御誤差zが有界であり及び状態変数xが有界であることを確保でき、
S3-2-3:式(18)に基づいて
Figure 2023517140000166
の時間に対する導関数
Figure 2023517140000167
を求めて、
Figure 2023517140000168
(ただし、
Figure 2023517140000169
であり、1つの正の定数rM3があり、0>r>rM3となる。)
式(10)、(24)、(11)及び(21)を(40)に代入して、
Figure 2023517140000170
を得て、
Figure 2023517140000171
とし、
(ただし、
Figure 2023517140000172
、ただし、V=V01+Ax、V=V02-Ax。)、実際の問題では、アクチュエータの安全な動作を考慮し、安全マージンが残され(すなわち、物理的構造に基づいて、油圧シリンダが中央位置の付近に最大12センチだけ上下に変動でき、与えられた指令信号が、一般的に10センチを超えない。)、このように、h,h,hが有界であることを確保でき、すなわち、それぞれ3つの正数
Figure 2023517140000173
があり、
Figure 2023517140000174
となり、
Figure 2023517140000175
であり、
Figure 2023517140000176
が有界であり、m(t)がスプールの不感帯の未知の傾きであり、m(t)が有界である(すなわち∃M>0,0<m(t)<M)ことを考慮すると、極値定理から、aが有界であることが分かり、
(41)と(42)とから、
Figure 2023517140000177
を得て、
リアプノフ関数
Figure 2023517140000178
を選択し、
(43)と(44)とは、不等式補助定理の形態を満たし、不等式補助定理から分かるように、aが有界∃A>0であり、|a|<A
Figure 2023517140000179
となるため、
Figure 2023517140000180
であり、
Figure 2023517140000181
が有界であり、(13)と(14)とから
Figure 2023517140000182
を得ることができるため、に対して、制御誤差zが有界であり及び主動力uが有界であることを確保でき、数学的モデルのHとHとの数式から見ると、状態変数xとxとが有界であり、いずれも油圧源Pの圧力値よりも小さく、
S3-2-4:式(18)に基づいて
Figure 2023517140000183
の時間に対する導関数
Figure 2023517140000184
を求めて、
Figure 2023517140000185
(ただし、
Figure 2023517140000186
であり、1つの正の定数rM4があり、0>r>rM4となる。)
式(10)、(24)及び(22)を(45)に代入して、
Figure 2023517140000187
を得て、
Figure 2023517140000188
とし、
Figure 2023517140000189
であり、
Figure 2023517140000190
が有界であり、kがコントローラを設計するときのオプションパラメータであることを考慮すると、aが有界であり、
(46)と(47)とから、
Figure 2023517140000191
を得て、
リアプノフ関数
Figure 2023517140000192
を選択し、
(48)と(49)とは、不等式補助定理の形態を満たし、不等式補助定理から分かるように、aが有界∃A>0であり、|a|<A
Figure 2023517140000193
となるため、
Figure 2023517140000194
であり、
Figure 2023517140000195
が有界であり、(13)と(14)とから
Figure 2023517140000196
を得ることができるため、∀t∈[0,τmax)に対して、制御誤差zが有界であり及び状態変数xが有界であることを確保できる。 Furthermore, if the normalized error vector is in the period t ∈ [0, τ max ) and there is a maximum solution in the non-empty open set Ω ζ , then the virtual controller and the controller of the servo-hydraulic system will determine t ∈ [0, τ max ) The process that can ensure that the closed-loop signal is bounded for is as follows:
S3-2-1: based on formula (18)
Figure 2023517140000139
the derivative with respect to time of
Figure 2023517140000140
looking for
Figure 2023517140000141
(however,
Figure 2023517140000142
, with one positive constant r M1 , such that 0>r 1 >r M1 . )
Substituting equations (10), (24) and (19) into (30),
Figure 2023517140000143
to get
Figure 2023517140000144
year,
Figure 2023517140000145
and
Figure 2023517140000146
is bounded, we know from the extremum theorem that a 1 is bounded,
From (31) and (32),
Figure 2023517140000147
to get
Lyapunov function
Figure 2023517140000148
and select
(33) and (34) satisfy the form of the inequality lemma, and as can be seen from the inequality lemma, a 1 is bounded ∃ A 1 > 0 and |a 1 | Since 1 (t)<r M1 k 1 ,
Figure 2023517140000149
and
Figure 2023517140000150
is bounded, and from (13) and (14),
Figure 2023517140000151
so that for ∀tε [0, τ max ), the tracking error z 1 is damped within the specified bounds, thus the tracking error z 1 can be guaranteed to converge within the specified performance bounds, and the command is a bounded sine signal, so the state variable x 1 is bounded,
S3-2-2: based on formula (18)
Figure 2023517140000152
the derivative with respect to time of
Figure 2023517140000153
looking for
Figure 2023517140000154
(however,
Figure 2023517140000155
and there is one positive constant r M2 , such that 0>r 2 >r M2 . )
Substituting equations (10), (24) and (20) into (35),
Figure 2023517140000156
to get
Figure 2023517140000157
year,
Figure 2023517140000158
and
Figure 2023517140000159
is bounded, f(t) is an unknown bounded disturbance term, and B is a bounded uncertainty parameter, then from the extremum theorem it follows that a 2 is bounded is understood,
From (36) and (37),
Figure 2023517140000160
to get
Lyapunov function
Figure 2023517140000161
and select
(38) and ( 39) satisfy the form of the inequality lemma, and as can be seen from the inequality lemma, a 2 is bounded ∃A 2 > 0 ,
Figure 2023517140000162
Because
Figure 2023517140000163
and
Figure 2023517140000164
is bounded, and from (13) and (14),
Figure 2023517140000165
, we can ensure that the control error z 2 is bounded and the state variable x 2 is bounded for ∀tε[0, τ max ),
S3-2-3: based on formula (18)
Figure 2023517140000166
the derivative with respect to time of
Figure 2023517140000167
looking for
Figure 2023517140000168
(however,
Figure 2023517140000169
and there is one positive constant r M3 such that 0>r 3 >r M3 . )
Substituting equations (10), (24), (11) and (21) into (40),
Figure 2023517140000170
to get
Figure 2023517140000171
year,
(however,
Figure 2023517140000172
, where V 1 =V 01 +A 1 x, V 2 =V 02 −A 2 x. ), in a practical matter, considering the safe operation of the actuator, a safety margin is left (i.e., based on the physical construction, the hydraulic cylinder can fluctuate up and down by a maximum of 12 cm around the center position, given command signals generally do not exceed 10 cm.), thus ensuring that h 1 , h 2 , and h 3 are bounded, i.e. each three positive numbers
Figure 2023517140000173
there is
Figure 2023517140000174
becomes,
Figure 2023517140000175
and
Figure 2023517140000176
is bounded, m(t) is the unknown slope of the dead zone of the spool, and m(t) is bounded (i.e. ∃M r >0,0<m(t)<M r ). Considering, from the extremum theorem, we know that a 3 is bounded,
From (41) and (42),
Figure 2023517140000177
to get
Lyapunov function
Figure 2023517140000178
and select
(43) and ( 44) satisfy the form of the inequality lemma, and as can be seen from the inequality lemma, a 3 is bounded ∃A 3 > 0 ,
Figure 2023517140000179
Because
Figure 2023517140000180
and
Figure 2023517140000181
is bounded, and from (13) and (14),
Figure 2023517140000182
, it can be ensured that the control error z 3 is bounded and the main power u s is bounded, and from the expressions of H 1 and H 2 in the mathematical model, , the state variables x3 and x4 are bounded, both smaller than the pressure value of the hydraulic source Ps ,
S3-2-4: based on formula (18)
Figure 2023517140000183
the derivative with respect to time of
Figure 2023517140000184
looking for
Figure 2023517140000185
(however,
Figure 2023517140000186
and there is one positive constant r M4 such that 0>r 4 >r M4 . )
Substituting equations (10), (24) and (22) into (45), we get
Figure 2023517140000187
to get
Figure 2023517140000188
year,
Figure 2023517140000189
and
Figure 2023517140000190
is bounded and k 4 is an optional parameter when designing the controller, a 4 is bounded and
From (46) and (47),
Figure 2023517140000191
to get
Lyapunov function
Figure 2023517140000192
and select
(48) and ( 49) satisfy the form of the inequality lemma, and as can be seen from the inequality lemma, a 4 is bounded ∃A 4 > 0 ,
Figure 2023517140000193
Because
Figure 2023517140000194
and
Figure 2023517140000195
is bounded, and from (13) and (14),
Figure 2023517140000196
, we can ensure that the control error z 4 is bounded and the state variable x 5 is bounded for ∀tε[0, τ max ).

さらに、初期値の提案を用いて、S3-2におけるrmax=+∞の場合、すべての閉ループ信号が有界であることが依然として正確であることを証明する前記プロセスは以下の通り、上証明から、

Figure 2023517140000197
が有界であり、すなわち
Figure 2023517140000198
であることが分かり、式(23)から、
Figure 2023517140000199
が空でない集合であり、
Figure 2023517140000200
が明らかになることが分かり、従って、τmax<+∞を仮定すると、初期値の提案では1つの時定数t∈[0,τmax)があることが強調され、
Figure 2023517140000201
となり、矛盾を顕著に引き起こす。従って、得られた正確な仮定がτmax=+∞であるため、いずれかのt∈[0,+∞)に対して、すべての閉ループシステム信号が有界であり、追跡誤差が指定された境界内で収束することを確保できる。 Further, using the initial value proposal, the above process of proving that it is still correct that all closed-loop signals are bounded for r max =+∞ in S3-2 follows the above proof from,
Figure 2023517140000197
is bounded, i.e.
Figure 2023517140000198
and from equation (23),
Figure 2023517140000199
is a non-empty set and
Figure 2023517140000200
is found to become apparent, thus emphasizing that the initial value proposal has one time constant t 1 ε[0, τ max ), assuming τ max <+∞,
Figure 2023517140000201
, which leads to significant contradiction. Therefore, for any tε[0,+∞), all closed-loop system signals are bounded, and the tracking error was specified because the exact assumption obtained is τ max =+∞ It can ensure convergence within bounds.

初期値問題を考慮すると、
ζ(t)=ν(t、ζ)、ζ(0)=ζ∈Ωζ(50)
(ただし、ν:R+×Ωζ→Rは連続関数ベクトルであり、Ωζ∈Rは空でない開集合である。)
定義1:初期値問題(50)の1つの解ζ(t)が、適切に大きくならない場合、該解は最大であり、
初期値定理としては、初期値問題(12)に対して、(1)t>0の場合、ζ、ν(t、ζ)に関して、ν(t、ζ)が部分Lipschitz条件を満たし、(2)ζ(t)∈Ωζに関して、ν(t、ζ)が区分的連続を満たし、(3)ζ(t)∈Ωζに関して、ν(t、ζ)がtに関して局所的に積分可能である場合、期間t∈[0、τmax)で、初期値問題(50)には1つの解ζ(t)∈Ωζがあり、ただしτmax>0である。初期値の提案としては、初期値定理が成立すると仮定し、期間[0、τmax)における最大解ζ(t)かつ集合

Figure 2023517140000202
に対して、τmax<∞の場合、時定数t∈[0、τmax)があり、
Figure 2023517140000203
となる。 Considering the initial value problem,
ζ(t)=ν(t, ζ), ζ(0)=ζ 0 ∈Ω ζ (50)
(where ν:R++Ω ζ →R n is a continuous function vector and Ω ζ εR n is a non-empty open set.)
Definition 1: A solution ζ(t) of the initial value problem (50) is maximal if it does not grow adequately, and
As the initial value theorem, for the initial value problem (12), (1) if t>0, ν(t, ζ) satisfies the partial Lipschitz condition for ζ, ν(t, ζ), and (2 ) for ζ(t)∈Ω ζ , ν(t, ζ) satisfies piecewise continuity, and (3) for ζ(t) ∈ Ω ζ , ν(t, ζ) is locally integrable with respect to t. If there is a period tε[0, τ max ), the initial value problem (50) has one solution ζ(t)εΩ ζ with τ max >0. As an initial value proposal, assuming that the initial value theorem holds, the maximum solution ζ(t) in the period [0, τ max ) and the set
Figure 2023517140000202
, there is a time constant t 1 ε[0, τ max ) for τ max <∞, and
Figure 2023517140000203
becomes.

さらに、前記不等式補助定理の式は、

Figure 2023517140000204
であり、
aが有界であり(すなわち∃A4>0,|a|<A)、
Figure 2023517140000205
である場合、xが必ず有界であり、
Figure 2023517140000206
である。 Furthermore, the expression of the above inequality lemma is
Figure 2023517140000204
and
a is bounded (i.e. ∃A>0, |a|<A),
Figure 2023517140000205
, then x must be bounded, and
Figure 2023517140000206
is.

本発明のコントローラの設計の正確さを検証するために、実験プラットフォームが構築され、実験ソフトウェアプログラムをデバッグし、実験プラットフォームの構築は、主にソフトウェア設計とハードウェア回路接続との2つの部分に分けられ、図3(a)は実験プラットフォームの構造構成図Iであり、図3(b)は実験プラットフォームの構造構成図IIであり、図3(c)は実験プラットフォームの構造構成図IIIであり、図3(d)は実験プラットフォームの構造構成図IVであり、実験プラットフォームのハードウェアは、主に3つの部分に分けられ、第1部分はアクチュエータ機構であり、主にアキュムレータ、油圧ポンプなどの油圧源システム、サーボバルブ及び単一出力ロッドの油圧シリンダアクチュエータを含み、本発明の実験プラットフォームでは、型番FD234-01K004VSX2Aの3ポジション5ポートサーボバルブが用いられる。第2部分は信号収集機構であり、信号収集機構のハードウェアは、主に信号変換ボード及びA/Dボードカードであり、信号変換ボードの機能は、電圧電流信号を互いに変換することであり、具体的には、センサ信号である4-20maの電流信号を1-5vの電圧信号に変換し、±5vの電圧信号を±10maの電流信号に変換し、本発明の実験プラットフォームで用いられるA/DボードカードはADT882であり、機能は、アナログ量の連続信号とデジタル量の散布度信号とを互いに変換することを実現することである。第3部分は制御機構であり、制御機構の中核は産業用制御コンピュータであり、本発明の実験プラットフォームでは産業用コンピュータpc104が用いられ、機能は、ソフトウェアプラットフォームの構築を実現し及び制御アルゴリズムを実現することである。ハードウェア回路接続とは、主に、各センサ信号線及びサーボバルブの制御電流信号線と信号変換ボードとの接続、信号変換ボードとADT882ボードカードとの接続を指す。ソフトウェア設計環境は、VC++6.0、windowオペレーティングシステムであり、ソフトウェアプログラムは、主に、インタフェース関数の編集、A/Dボードカードの配置、信号収集を完了するための割り込みプログラムの設定、制御量の計算及び出力を含む。 In order to verify the correctness of the design of the controller of the present invention, an experimental platform is built and the experimental software program is debugged. The construction of the experimental platform is mainly divided into two parts: software design and hardware circuit connection. 3(a) is a structural diagram I of the experimental platform, FIG. 3(b) is a structural diagram II of the experimental platform, FIG. 3(c) is a structural diagram III of the experimental platform, Figure 3(d) is the structural configuration diagram IV of the experimental platform, the hardware of the experimental platform is mainly divided into three parts, the first part is the actuator mechanism, mainly the hydraulic A 3-position, 5-port servovalve model number FD234-01K004VSX2A is used in the experimental platform of the present invention, including a source system, a servovalve and a single output rod hydraulic cylinder actuator. The second part is the signal acquisition mechanism, the hardware of the signal acquisition mechanism is mainly the signal conversion board and the A/D board card, the function of the signal conversion board is to convert voltage current signals to each other, Specifically, the current signal of 4-20 ma, which is a sensor signal, is converted into a voltage signal of 1-5 v, the voltage signal of ±5 v is converted into a current signal of ±10 ma, and the A The /D board card is ADT882, the function is to realize the conversion between analog quantity continuous signal and digital quantity dispersion signal. The third part is the control mechanism, the core of the control mechanism is the industrial control computer, the experimental platform of the present invention uses the industrial computer pc104, the function is to realize the construction of the software platform and the control algorithm It is to be. The hardware circuit connection mainly refers to the connection between each sensor signal line and the control current signal line of the servo valve and the signal conversion board, and the connection between the signal conversion board and the ADT882 board card. The software design environment is VC++ 6.0, window operating system. Includes calculations and outputs.

さらに、実験結果が予期される制御効果に達するまで、実験をデバッグし、パラメータを最適化し、ステップ1、本発明のアルゴリズムの指定された境界の境界初期値ρi0、誤差収束レートの上限h及び収束誤差の定常状態残差セットの上限値ρi∞を調整し、ただしi=1…4であり、境界初期値ρi0を可能な限り大きくし、誤差収束レートの上限hを可能な限り小さくし、収束誤差の定常状態残差セットの上限値ρi∞を可能な限り大きくし、ステップ2、制御率ゲインkが適切な値に調整されるまで、仮想制御率ゲインk、k、k、kを調整し、ステップ3、予期される制御効果に達するまで、境界初期値ρi0を緩慢かつ適切に小さくし、誤差収束レートの上限hを大きくし、収束誤差の定常状態残差セットの上限値ρi∞を小さくする。 Further, until the experimental results reach the expected control effect, debug the experiment, optimize the parameters, step 1, bound initial values ρ i0 for the specified bounds of our algorithm, upper bound h i for the error convergence rate and the upper bound ρ i of the steady -state residual set of convergence errors, where i=1 . as small as possible, the upper limit ρ i∞ of the steady-state residual set of convergence errors is made as large as possible, step 2, until the control ratio gains k i are adjusted to appropriate values, the virtual control ratio gains k 1 , Adjust k 2 , k 3 , k 4 , step 3, slowly and appropriately decrease the boundary initial value ρ i0 , increase the upper bound h i of the error convergence rate, and increase the convergence error Decrease the upper bound ρ i∞ of the steady-state residual set of .

図4は本発明のコントローラの作用下での誤差収束のシミュレーション曲線図であり、図4から、本発明のコントローラは、変位追跡誤差の収束速度と制御精度とを確保できることが分かり、図5は外乱作用下での追跡誤差収束のシミュレーション曲線図であり、図5では、システムモデルの強い外乱f(t)=9000sin(10t)の作用下での本発明のコントローラとバックステッピングコントローラとの変位追跡誤差の比較曲線が示されており、図4と図5の比較から、強い外乱の存在は、バックステッピングコントローラの制御効果に大きく影響され、変位追跡誤差を大きくし、変動周波数を高くすることが分かる。しかしながら、図4及び図5から、本発明のコントローラは、強い外乱に対する抑制作用が強く、依然として収束誤差の過渡と定常状態との性能を確保できることも分かる。 Figure 4 is a simulation curve diagram of error convergence under the action of the controller of the present invention. Fig. 5 is a simulation curve diagram of the tracking error convergence under the action of disturbance, Fig. 5 shows the displacement tracking of the controller of the present invention and the backstepping controller under the action of a strong disturbance f(t) = 9000 sin(10t) of the system model; A comparison curve of errors is shown, and from a comparison of FIG. 4 and FIG. 5, the presence of a strong disturbance is greatly affected by the control effect of the backstepping controller, and can increase the displacement tracking error and increase the fluctuation frequency. I understand. However, from FIGS. 4 and 5, it can also be seen that the controller of the present invention has a strong suppression effect against strong disturbances, and can still ensure the transient and steady state performance of the convergence error.

図6は未知のスプールの不感帯の作用下での誤差収束のシミュレーション曲線図であり、図7は未知のスプールの不感帯の作用下でのスプール変位のシミュレーション曲線図であり、図6及び図7では、本明細書で提案されている低複雑度制御方案が未知のスプールの不感帯を処理するときにロバスト性が高いことが検証されている。システムモデルにスプールの未知の不感帯の非線形が追加され、スプールの不感帯の未知の非線形特性をより良好に具現化するために、本明細書では、時変不感帯モデルの傾きm(t)=1+0.3sin(2t)、時変不感帯モデルの不連続ポイントbl=br=0.3|sin(2t)|が用いられ、図7から、本発明のコントローラは、スプールの未知の不感帯の悪影響を効果的に補償できることが分かり、図6から、システムには未知のスプールの不感帯問題が存在するが、本発明のコントローラは、依然として追跡誤差の収束速度と定常状態の精度とを良好に確保できることが分かる。 FIG. 6 is a simulation curve diagram of error convergence under the action of unknown spool dead band, and FIG. 7 is a simulation curve diagram of spool displacement under the action of unknown spool dead band. , it has been verified that the low-complexity control scheme proposed herein is highly robust when dealing with unknown spool deadbands. The nonlinearity of the unknown spool deadband is added to the system model, and in order to better embody the unknown nonlinearity of the spool deadband, the slope of the time-varying deadband model m(t)=1+0 . 3sin(2t), the discontinuity point bl=br=0.3|sin(2t)| of the time-varying deadband model is used and from FIG. From FIG. 6, it can be seen that although the system has an unknown spool dead-band problem, the controller of the present invention still provides good tracking error convergence speed and steady-state accuracy.

図8は本発明の制御方法と、指定された性能を満たす、スプールの未知の不感帯を有する適応バックステッピング制御方式(SPPFBSA)と、の誤差収束の比較シミュレーション曲線図であり、図9は本発明の制御方法とSPPFBSAとのシミュレーション計算時間及び誤差調整時間の統計図であり、図8から、2種のコントローラはいずれも変位追跡誤差の収束速度と定常状態の精度とを確保できることが分かるが、図9から、本発明のコントローラは、SPPFBSAコントローラと比較して、シミュレーション実行時間を91.7%短縮し、誤差収束の調整時間を95.5%短縮することが分かる。本発明のコントローラの設計は、モデルへの依存性が低いため、バックステッピング適応に基づいて開発されたアルゴリズムと比較して、計算量が少なく、オンライン学習が不要であるため、リアルタイム制御が容易であり、工学的実現が容易である。 FIG. 8 is a comparison simulation curve diagram of error convergence between the control method of the present invention and an adaptive backstepping control scheme with unknown spool deadband (SPPFBSA) that satisfies the specified performance, and FIG. It is a statistical chart of the simulation calculation time and error adjustment time of the control method and SPPFBSA, and from FIG. From FIG. 9, it can be seen that the controller of the present invention reduces the simulation run time by 91.7% and the error convergence adjustment time by 95.5% compared to the SPPFBSA controller. The design of the controller of the present invention is less model-dependent, thus less computationally intensive compared to algorithms developed based on backstepping adaptation, and does not require online learning, thus facilitating real-time control. and is easy to implement engineeringally.

図10は指令信号がTS=3s:x=100sin((2*pi*t)/3)mmであるときの、本発明の制御方法とpidとの変位追跡の比較実験図であり、図11は指令信号がTS=3s:x=100sin((2*pi*t)/3)mmであるときの、本発明の制御方法の追跡誤差の実験図であり、図12は指令信号がTS=2s:x=100sin(pi*t)mmであるときの、本発明の制御方法とpidとの変位追跡の比較実験図であり、図13は指令信号がTS=2s:x=100sin(pi*t)mmであるときの、本発明の制御方法の追跡誤差の実験図であり、図14は指令信号がTS=1s:x=50sin(2*pi*t)mmであるときの、本発明の制御方法とpidとの変位追跡の比較実験図であり、図15は指令信号がTS=1s:x=50sin(2*pi*t)mmであるときの、本発明の制御方法の追跡誤差の実験図であり、図16は指令信号がTS=2s:x=100sin(pi*t)mmであるときの、本発明の制御方法とpidとの変位追跡の比較実験図であり、図17は指令信号がTS=2s:x=100sin(pi*t)mmであるときの、本発明の制御方法の追跡誤差の実験図であり、図18は指令信号がTS=2s:x=80sin(pi*t)mmであるときの、本発明の制御方法とpidとの変位追跡の比較実験図であり、図19は指令信号がTS=2s:x=80sin(pi*t)mmであるときの、本発明の制御方法の追跡誤差の実験図であり、図20は指令信号がTS=2s:x=60sin(pi*t)mmであるときの、本発明の制御方法とpidとの変位追跡の比較実験図であり、図21は指令信号がTS=2s:x=60sin(pi*t)mmであるときの、本発明の制御方法の追跡誤差の実験図であり、図11、図13、図15、図17、図19及び図21から分かるように、指令信号の周波数及び振幅がどのように変化するかにかかわらず、本発明のコントローラは、追跡誤差が指定された性能境界内で収束することを確保でき、従って、本発明のコントローラは、追跡誤差の収束速度と定常状態の精度とを確保できる。図10、図12、図14、図16、図18及び図20のいずれかから分かるように、本発明のコントローラにより得られた変位追跡曲線は、シリンダの上昇プロセスでほとんど小さい位相遅れしかなく、シリンダの上昇プロセスと比較して、シリンダが下降する場合、変位追跡曲線の位相遅れ程度が大きいが、本発明のコントローラの制御効果について、変位追跡のプロセス全体で、変位追跡曲線は、いずれもpidの変位追跡曲線よりも位相遅れ程度が小さくなる。 FIG. 10 is a comparison experimental diagram of displacement tracking between the control method of the present invention and pid when the command signal is TS=3s: xr =100sin ((2*pi*t)/3)mm. 11 is an experimental diagram of the tracking error of the control method of the present invention when the command signal is TS=3s: xr =100sin ((2*pi*t)/3) mm, and FIG. FIG. 13 is a comparative experiment diagram of displacement tracking between the control method of the present invention and pid when TS=2s: xr =100sin(pi*t) mm, and FIG . FIG. 14 is an experimental diagram of the tracking error of the control method of the present invention when 100 sin(pi*t) mm, FIG. 14 is the command signal TS=1 s:x r =50 sin(2*pi*t) mm FIG . 15 is a comparison experiment diagram of displacement tracking between the control method of the present invention and pid when the FIG. 16 is a comparison of the displacement tracking between the control method of the present invention and pid when the command signal is TS=2s: xr =100sin(pi*t)mm FIG. 17 is an experimental diagram of the tracking error of the control method of the present invention when the command signal is TS=2s: xr =100sin(pi*t) mm, and FIG. FIG. 19 is a comparative experiment diagram of displacement tracking between the control method of the present invention and pid when TS=2s: xr =80sin(pi*t) mm, and FIG . FIG. 20 is an experimental diagram of the tracking error of the control method of the present invention when the command signal is TS=2s:x r =60 sin(pi*t) mm; 21 is a comparative experiment diagram of displacement tracking between the control method of the present invention and pid , and FIG. 11, 13, 15, 17, 19 and 21 are experimental plots of tracking error, regardless of how the frequency and amplitude of the command signal varies. The controller can ensure that the tracking error converges within the specified performance bounds, and thus the controller of the present invention can ensure tracking error convergence speed and steady-state accuracy. As can be seen from any of Figures 10, 12, 14, 16, 18 and 20, the displacement tracking curves obtained by the controller of the present invention have very little phase lag in the cylinder lift process, Compared with the cylinder rising process, the phase lag of the displacement tracking curve is large when the cylinder is descending, but for the control effect of the controller of the present invention, the displacement tracking curve is pid The phase lag is smaller than the displacement tracking curve of .

図10~図13を組み合わせると、指令信号の周波数が高いほど、pid追跡曲線の位相遅れ程度が高くなり、振幅が大幅に減衰される。定常状態の場合、本発明のコントローラの制御効果について、シリンダの上昇プロセスで、位相遅れ程度が基本的に変化せず、シリンダの下降プロセスで、位相遅れ程度が小さくなり、振幅が基本的に減衰されない。指令信号の周波数が高いほど、pid追跡誤差が大きくなるが、本発明のコントローラの追跡誤差が基本的に変化せず、指定された性能境界内で収束する。 Combining FIGS. 10-13, the higher the frequency of the command signal, the higher the degree of phase lag in the pid tracking curve and the greater the amplitude attenuation. In the steady state, regarding the control effect of the controller of the present invention, the degree of phase lag basically does not change during the upward process of the cylinder, and the degree of phase lag becomes smaller and the amplitude is basically attenuated during the downward process of the cylinder. not. The higher the frequency of the command signal, the larger the pid tracking error, but the tracking error of the controller of the present invention remains essentially unchanged and converges within the specified performance bounds.

図17、図19及び図21を組み合わせると、指令信号の周波数が高いほど、pid追跡誤差が大きくなるが、本発明のコントローラの追跡誤差が基本的に変化せず、指定された性能境界内で収束する。 Combining Figures 17, 19 and 21, the higher the frequency of the command signal, the greater the pid tracking error, but the tracking error of the controller of the present invention remains essentially unchanged and within the specified performance bounds. converge.

以上は、本発明の好ましい具体的な実施形態に過ぎないが、本発明の保護範囲はこれに限定されず、当業者が、本発明に開示されている技術範囲内で、本発明の技術的解決手段及びその発明構想に基づいて行った等価置換や変更は、いずれも本発明の保護範囲内に含まれるべきである。 The above are only preferred and specific embodiments of the present invention, but the protection scope of the present invention is not limited thereto. Any equivalent replacement or modification made on the basis of the solution and its inventive concept shall fall within the protection scope of the present invention.

Claims (8)

非対称サーボ油圧位置追跡システムのための低複雑度制御方法であって、
単一出力ロッドのサーボ油圧システムのモデルを確立するステップS1と、
単一出力ロッドのサーボ油圧システムのモデルに基づいて、低複雑度制御ポリシーを用いて単一出力ロッドのサーボ油圧システムのコントローラを設計するステップS2と、
単一出力ロッドのサーボ油圧システムのコントローラ及び単一出力ロッドのサーボ油圧システムのモデルに基づいて、単一出力ロッドのサーボ油圧システムの安定性を証明するステップS3と、を含む、ことを特徴とする非対称サーボ油圧位置追跡システムのための低複雑度制御方法。
A low-complexity control method for an asymmetric servohydraulic position tracking system, comprising:
Step S1 of establishing a model of a single output rod servo-hydraulic system;
designing a controller of the single output rod servo-hydraulic system using a low-complexity control policy based on the model of the single output rod servo-hydraulic system;
and verifying the stability of the single output rod servo-hydraulic system based on the controller of the single output rod servo-hydraulic system and the model of the single output rod servo-hydraulic system. low-complexity control method for an asymmetric servohydraulic position tracking system.
単一出力ロッドのサーボ油圧システムのモデルを確立する前記プロセスは以下の通り、
ニュートンの第2法則に基づいて油圧シリンダの動力学モデルを確立し、
Figure 2023517140000207
(ただし、f(t)は様々な干渉を表し、xとmとはそれぞれ負荷の位置と質量とを表し、Bは粘性減衰係数であり、Kは負荷の等価ばね剛性であり、負荷が慣性負荷の場合、K=0であり、F=A*P-A*Pは油圧アクチュエータにより出力される主動力であり、ただし、P、Pは油圧シリンダの大きなキャビティと小さなキャビティとの圧力であり、A、Aは大きなキャビティと小さなキャビティとのピストンの有効面積である。)
3ポジション5ポートサーボバルブを用い、負荷圧力の動力学は下記式で表され、
Figure 2023517140000208
Figure 2023517140000209
(式中、V、Vはそれぞれロッドキャビティとロッドレスキャビティとの容積であり、V=V01+Axであり、V=V02-Axであり、V01とV02とはそれぞれピストンが初期位置にあるときのロッドレスキャビティとロッドキャビティとの容積であり、Cは油圧シリンダの内部漏れ係数であり、Cは油圧シリンダの外部漏れ係数であり、βはオイルの弾性率であり、Qはロッドキャビティの作動油の流量であり、Qはロッドレスキャビティの作動油の流量である。)
ただし、
Figure 2023517140000210
Figure 2023517140000211
(ただし、
Figure 2023517140000212
は流量ゲインであり、s(Γ(x))の式は
Figure 2023517140000213
であり、
は油圧システムの給油圧力であり、Pは油圧システムの油戻し圧力であり、Cはオリフィスの流量係数であり、wはスライドバルブの面積勾配であり、ρはオイルの密度であり、xはサーボバルブのスプールである。)
電圧又は電流入力uによりサーボバルブのスプールの変位xを制御して、必要な対応する力をさらに取得し、サーボバルブの動的特性は以下のように示され、
Figure 2023517140000214
(τはサーボバルブの動力学モデルの時定数であり、u(t)は電流入力である。)
未知の不感帯があるスプールの変位Γ(x)を考慮すると、その式は
Figure 2023517140000215
であり、
(パラメータmとmとは不感帯特性曲線の左右の傾きを表し、パラメータbとbとは入力非線形のブレークポイントを表す。)
スプールの不感帯はΓ(x)=m(t)x+d(t)で表され、m(t)、d(t)の式は以下のように示され、
Figure 2023517140000216
Figure 2023517140000217
(m、m、b、bが有界であるため、d(t)が有界であり、d(t)の上限を
Figure 2023517140000218
とする。)
状態変数を
Figure 2023517140000219
のように定義し、未知のスプールの不感帯を有する以下の位置追跡システムの状態方程式を取得する
Figure 2023517140000220
(ただし、
Figure 2023517140000221
。)、ことを特徴とする請求項1に記載の非対称サーボ油圧位置追跡システムのための低複雑度制御方法。
The process of establishing a model of a single output rod servo-hydraulic system is as follows:
Establishing a hydraulic cylinder dynamics model based on Newton's second law,
Figure 2023517140000207
(where f(t) represents the various interferences, x and m represent the position and mass of the load, respectively, B is the viscous damping coefficient, K is the equivalent spring stiffness of the load, and the load is inertial For the load, K=0 and F=A 1 *P 1 -A 2 *P 2 is the main force output by the hydraulic actuator, where P 1 , P 2 are the large cavity and small is the pressure with the cavity, and A 1 , A 2 are the effective areas of the piston with the large and small cavities.)
Using a 3-position 5-port servo valve, the dynamics of the load pressure is expressed by the following equation,
Figure 2023517140000208
Figure 2023517140000209
(where V 1 and V 2 are the volumes of the rod cavity and rodless cavity respectively, V 1 =V 01 +A 1 x, V 2 =V 02 -A 2 x, V 01 and V 02 is the volume of the rodless cavity and the rod cavity respectively when the piston is in the initial position, C t is the internal leakage coefficient of the hydraulic cylinder, C e is the external leakage coefficient of the hydraulic cylinder, β e is the elastic modulus of the oil, Q1 is the hydraulic fluid flow rate in the rod cavity, and Q2 is the hydraulic fluid flow rate in the rodless cavity.)
however,
Figure 2023517140000210
Figure 2023517140000211
(however,
Figure 2023517140000212
is the flow gain and the expression for s(Γ(x v )) is
Figure 2023517140000213
and
P s is the oil supply pressure of the hydraulic system, P r is the oil return pressure of the hydraulic system, C d is the flow coefficient of the orifice, w is the area gradient of the slide valve, and ρ is the density of the oil. , xv are the spools of the servovalves. )
Controlling the displacement x v of the spool of the servo valve by voltage or current input u to further obtain the required corresponding force, the dynamic characteristics of the servo valve are shown as:
Figure 2023517140000214
(τ is the time constant of the dynamic model of the servovalve and u(t) is the current input.)
Considering the spool displacement Γ(x v ) with an unknown dead band, the expression is
Figure 2023517140000215
and
(The parameters mr and ml represent the left and right slopes of the deadband characteristic curve, and the parameters br and bl represent the breakpoints of the input nonlinearity.)
The dead band of the spool is expressed by Γ( xv )=m(t) xv +d(t), and the formulas for m(t) and d(t) are shown below,
Figure 2023517140000216
Figure 2023517140000217
Since (m r , m l , b r , b l are bounded, d(t) is bounded and the upper bound of d(t) is given by
Figure 2023517140000218
and )
the state variable
Figure 2023517140000219
and obtain the state equation of the following position tracking system with unknown spool deadband
Figure 2023517140000220
(however,
Figure 2023517140000221
. ).
単一出力ロッドのサーボ油圧システムのモデルに基づいて、低複雑度制御ポリシーを用いて単一出力ロッドのサーボ油圧システムのコントローラを設計する前記プロセスは以下の通り、
誤差変数を定義し、
=x-x,z=x-α,z=F-α,z=x-α(17)
(ただし、α、α、αは仮想制御量であり、X1rは油圧位置追跡システムの指令信号であり、F(F=A-A)はアクチュエータにより出力される主動力であり、zは位置追跡誤差であり、z,i=2…4は制御誤差である。)
4つの正の滑らかな減少関数
Figure 2023517140000222
を指定された性能関数として選択し、正規化誤差ζを定義し、z,i=1…4に対して誤差変換を行い、変換された誤差
Figure 2023517140000223

Figure 2023517140000224
であり、
以下の仮想コントローラを選択し、
Figure 2023517140000225
Figure 2023517140000226
Figure 2023517140000227
サーボ油圧システムのコントローラは、
Figure 2023517140000228
である、ことを特徴とする請求項1に記載の非対称サーボ油圧位置追跡システムのための低複雑度制御方法。
The process of designing a controller for a single output rod servo-hydraulic system using a low-complexity control policy based on a model of a single output rod servo-hydraulic system is as follows:
Define an error variable,
z 1 =x 1 −x r ,z 2 =x 2 −α 1 ,z 3 =F −α 2 ,z 4 =x 5 −α 3 (17)
(where α 1 , α 2 , α 3 are virtual control variables, X 1r is the command signal of the hydraulic position tracking system, and F (F = A 1 x 3 - A 2 x 4 ) is output by the actuator. is the prime mover, z 1 is the position tracking error, and z i , i=2 . . . 4 are the control errors.)
4 positive smooth decreasing functions
Figure 2023517140000222
as the specified performance function, define the normalized errors ζ i , perform an error transformation on z i , i=1 . . .
Figure 2023517140000223
teeth
Figure 2023517140000224
and
Select the virtual controller below,
Figure 2023517140000225
Figure 2023517140000226
Figure 2023517140000227
The controller of the servo hydraulic system is
Figure 2023517140000228
A low-complexity control method for an asymmetric servohydraulic position tracking system according to claim 1, characterized in that:
単一出力ロッドのサーボ油圧システムのコントローラ及び単一出力ロッドのサーボ油圧システムのモデルに基づいて、単一出力ロッドのサーボ油圧システムの安定性を証明する前記プロセスは以下の通り、
S3-1:初期値定理を用い、正規化誤差ベクトルが期間t∈[0,τmax)にある場合、空でない開集合Ωζに最大解があることを証明し、
S3-2:正規化誤差ベクトルが期間t∈[0,τmax)にあり、空でない開集合Ωζに最大解がある場合、仮想コントローラ及びサーボ油圧システムのコントローラは、t∈[0,τmax)に対して閉ループ信号が有界であることを確保でき、
S3-3:初期値の提案を用いて、S3-2におけるτmax=+∞の場合、すべての閉ループ信号が有界であることが依然として正確であることを証明する、ことを特徴とする請求項1に記載の非対称サーボ油圧位置追跡システムのための低複雑度制御方法。
Based on the controller of the single output rod servo-hydraulic system and the model of the single output rod servo-hydraulic system, the above process of proving the stability of the single output rod servo-hydraulic system is as follows:
S3-1: Using the initial value theorem, prove that a non-empty open set Ω ζ has a maximum solution if the normalized error vector is in the period tε[0, τ max ),
S3-2: If the normalized error vector is in the period t ∈ [0, τ max ) and there is a maximum solution in the non-empty open set Ω ζ , the virtual controller and the controller of the servo-hydraulic system, t ∈ [0, τ max ) to ensure that the closed-loop signal is bounded,
S3-3: Using the initial value proposition, prove that it is still correct that all closed-loop signals are bounded if τ max =+∞ in S3-2. A low-complexity control method for an asymmetric servohydraulic position tracking system of claim 1.
初期値定理を用い、正規化誤差ベクトルが期間t∈[0,τmax)にある場合、空でない開集合Ωζに最大解があることを証明する前記ステップにおいては、
空でない開集合Ωζは、性能関数ρを選択する場合、
Figure 2023517140000229
を満たすことにより、
Figure 2023517140000230
を得ることができるため、正規化誤差ベクトルはζ(0)∈Ωζであり、所望の追跡奇跡x、性能関数ρ(t),i=1…4、中間制御信号α,i=1…3及び制御率uが連続的で滑らかに導出可能であり、スプールの不感帯の位置追跡システムの状態方程式における動力学変数が連続的に導出可能な関数であるため、式の正規化誤差ベクトルの関数L(t,ζ)は、時間tに関して区分的に連続し、開集合Ωζでは、ζはリプシッツ条件を満たし、初期値定理における条件がいずれも満たされるため、期間[0,τmax)で、正規化誤差ベクトルには唯一の最大解ζ(t)∈Ωζがあり、∀t∈[0,τmax)に対して、いずれもζ(t)∈Ωζ,i=1…4を確保できる、ことを特徴とする請求項4に記載の非対称サーボ油圧位置追跡システムのための低複雑度制御方法。
In the step of using the initial value theorem to prove that a non-empty open set Ω ζ has a maximum solution if the normalized error vector is in the period tε[0, τ max ):
A non-empty open set Ω ζ , if we choose a performance function ρ i ,
Figure 2023517140000229
By satisfying
Figure 2023517140000230
, the normalized error vector is ζ(0)∈Ω ζ , the desired tracking miracle x r , the performance function ρ i (t), i=1 . . . 4, the intermediate control signals α i , i = 1 ... 3 and the control rate u are continuously and smoothly derivable, and the dynamics variables in the state equation of the position tracking system of the dead zone of the spool are continuously derivable functions, so the normalization error of the equation The vector function L(t, ζ) is piecewise continuous with respect to time t, and for an open set Ω ζ , ζ satisfies the Lipschitz condition and both conditions in the initial value theorem are satisfied, so the period [0, τ max ), the normalized error vector has a unique maximum solution ζ(t)∈Ω ζ , and for ∀t∈[0, τ max ), any ζ i (t)∈Ω ζ , i= A low-complexity control method for an asymmetric servohydraulic position tracking system according to claim 4, characterized in that 1...4 can be ensured.
正規化誤差ベクトルが期間t∈[0,τmax)にあり、空でない開集合Ωζに最大解がある場合、仮想コントローラ及びサーボ油圧システムのコントローラは、t∈[0,τmax)に対して閉ループ信号が有界であることを確保できる前記プロセスは以下の通り、
S3-2-1:式(18)に基づいて
Figure 2023517140000231
の時間に対する導関数
Figure 2023517140000232
を求めて、
Figure 2023517140000233
(ただし、
Figure 2023517140000234
であり、1つの正の定数rM1があり、0>r>rM1となる。)
さらに、
Figure 2023517140000235
を得ることができ、
Figure 2023517140000236
とし、
Figure 2023517140000237
であり、
Figure 2023517140000238
が有界であることを考慮すると、極値定理から、aが有界であることが分かり、
(31)と(32)とから、
Figure 2023517140000239
を得て、
リアプノフ関数
Figure 2023517140000240
を選択し、
(33)と(34)とは、不等式補助定理の形態を満たし、不等式補助定理から分かるように、aが有界であり、すなわち∃A>0であり、|a|<A、0<Δ(t)<rM1となるため、
Figure 2023517140000241
であり、
Figure 2023517140000242
が有界であり、
Figure 2023517140000243
を得ることができるため、∀t∈[0,τmax)に対して、追跡誤差zが指定された境界内で減衰され、期間t∈[0,τmax)での追跡誤差zが指定された性能境界内で収束することを確保し、指令が有界正弦信号であるため、状態変数xが有界であり、
S3-2-2:式(18)に基づいて
Figure 2023517140000244
の時間に対する導関数
Figure 2023517140000245
を求めて、
Figure 2023517140000246
(ただし、
Figure 2023517140000247
であり、1つの正の定数rM2があり、0>r>rM2となる。)
さらに、
Figure 2023517140000248
を得ることができ、
Figure 2023517140000249
とし、
Figure 2023517140000250
であり、
Figure 2023517140000251
が有界であり、f(t)が未知の有界な外乱項であり、Bが有界な不確実性パラメータであることを考慮すると、極値定理から、aが有界であることが分かり、
(36)と(37)とから、
Figure 2023517140000252
を得て、
リアプノフ関数
Figure 2023517140000253
を選択し、
(38)と(39)とは、不等式補助定理の形態を満たし、不等式補助定理から分かるように、aが有界であり、すなわち∃A>0であり、|a|<A
Figure 2023517140000254
となるため、
Figure 2023517140000255
であり、
Figure 2023517140000256
が有界であり、
Figure 2023517140000257
を得ることができるため、∀t∈[0,τmax)に対して、制御誤差zが有界であり及び状態変数xが有界であることを確保し、
S3-2-3:式(18)に基づいて
Figure 2023517140000258
の時間に対する導関数
Figure 2023517140000259
を求めて、
Figure 2023517140000260

(ただし、
Figure 2023517140000261
であり、1つの正の定数rM3があり、0>r>rM3となる。)
式(10)、(24)、(11)及び(21)を(40)に代入して、
Figure 2023517140000262
を得て、
Figure 2023517140000263
とし、

Figure 2023517140000264
、ただし、V=V01+Ax、V=V02-Ax。)
サーボ油圧システムアクチュエータの安全な動作に鑑みて、安全マージンが残され、すなわち、物理的構造に基づいて、油圧シリンダが中央位置の付近に最大12センチだけ上下に変動でき、指令信号の振幅が10センチ以下であり、h,h,hが有界であることを確保し、すなわち、それぞれ3つの正数
Figure 2023517140000265
があり、
Figure 2023517140000266
となり、
Figure 2023517140000267
であり、
Figure 2023517140000268
が有界であり、m(t)がスプールの不感帯の未知の傾きであり、m(t)が有界である(すなわち∃M>0,0<m(t)<M)ことを考慮すると、極値定理から、aが有界であることが分かり、
(41)と(42)とから、
Figure 2023517140000269
を得て、
リアプノフ関数
Figure 2023517140000270
を選択し、
(43)と(44)とは、不等式補助定理の形態を満たし、不等式補助定理から分かるように、aが有界であり、すなわち∃A>0であり、|a|<A
Figure 2023517140000271
となるため、
Figure 2023517140000272
であり、
Figure 2023517140000273
が有界であり、
Figure 2023517140000274
を得ることができるため、∀t∈[0,τmax)に対して、制御誤差zが有界であり及び主動力uが有界であることを確保し、状態方程式の数式から見ると、状態変数xとxとが有界であり、いずれも油圧源Pの圧力値よりも小さく、
S3-2-4:式(18)に基づいて
Figure 2023517140000275
の時間に対する導関数
Figure 2023517140000276
を求めて、
Figure 2023517140000277

(ただし、
Figure 2023517140000278
であり、1つの正の定数rM4があり、0>r>rM4となる。)
さらに、
Figure 2023517140000279
を得て、
Figure 2023517140000280
とし、
Figure 2023517140000281
であり、
Figure 2023517140000282
が有界であり、kがコントローラを設計するときのオプションパラメータであることを考慮すると、aが有界であり、
(46)と(47)とから、
Figure 2023517140000283
を得て、
リアプノフ関数
Figure 2023517140000284
を選択し、
(48)と(49)とは、不等式補助定理の形態を満たし、不等式補助定理から分かるように、aが有界であり、すなわち∃A>0であり、|a|<A
Figure 2023517140000285
となるため、
Figure 2023517140000286
であり、
Figure 2023517140000287
が有界であり、
Figure 2023517140000288
を得ることができるため、に対して、制御誤差zが有界であり及び状態変数xが有界であることを確保する、ことを特徴とする請求項4に記載の非対称サーボ油圧位置追跡システムのための低複雑度制御方法。
If the normalized error vector is in the period t ∈ [0, τ max ) and the non-empty open set Ω ζ has a maximum solution, then the virtual controller and the controller of the servo-hydraulic system, for t ∈ [0, τ max ) The process that can ensure that the closed-loop signal is bounded is as follows:
S3-2-1: based on formula (18)
Figure 2023517140000231
the derivative with respect to time of
Figure 2023517140000232
looking for
Figure 2023517140000233
(however,
Figure 2023517140000234
, with one positive constant r M1 , such that 0>r 1 >r M1 . )
moreover,
Figure 2023517140000235
can get
Figure 2023517140000236
year,
Figure 2023517140000237
and
Figure 2023517140000238
is bounded, we know from the extremum theorem that a 1 is bounded,
From (31) and (32),
Figure 2023517140000239
to get
Lyapunov function
Figure 2023517140000240
and select
(33) and (34) satisfy the form of the inequality lemma, and as can be seen from the inequality lemma, a 1 is bounded, i.e. ∃A 1 >0 and |a 1 |<A 1 , 0<Δ 1 (t)<r M1 k 1 , so
Figure 2023517140000241
and
Figure 2023517140000242
is bounded, and
Figure 2023517140000243
so that for ∀tε[0, τ max ), the tracking error z 1 decays within the specified bounds, and the tracking error z 1 over time period tε[0, τ max ) is ensuring convergence within the specified performance bounds, the state variable x1 is bounded because the command is a bounded sine signal,
S3-2-2: based on formula (18)
Figure 2023517140000244
the derivative with respect to time of
Figure 2023517140000245
looking for
Figure 2023517140000246
(however,
Figure 2023517140000247
and there is one positive constant r M2 , such that 0>r 2 >r M2 . )
moreover,
Figure 2023517140000248
can get
Figure 2023517140000249
year,
Figure 2023517140000250
and
Figure 2023517140000251
is bounded, f(t) is an unknown bounded disturbance term, and B is a bounded uncertainty parameter, then from the extremum theorem it follows that a 2 is bounded is understood,
From (36) and (37),
Figure 2023517140000252
to get
Lyapunov function
Figure 2023517140000253
and select
(38) and (39) satisfy the form of the inequality lemma, and as can be seen from the inequality lemma, a 2 is bounded, i.e. ∃A 2 >0 and |a 2 |<A 2 ,
Figure 2023517140000254
Because
Figure 2023517140000255
and
Figure 2023517140000256
is bounded, and
Figure 2023517140000257
ensuring that the control error z 2 is bounded and the state variable x 2 is bounded for ∀tε[0, τ max ),
S3-2-3: based on formula (18)
Figure 2023517140000258
the derivative with respect to time of
Figure 2023517140000259
looking for
Figure 2023517140000260

(however,
Figure 2023517140000261
and there is one positive constant r M3 such that 0>r 3 >r M3 . )
Substituting equations (10), (24), (11) and (21) into (40),
Figure 2023517140000262
to get
Figure 2023517140000263
year,
(
Figure 2023517140000264
, where V 1 =V 01 +A 1 x, V 2 =V 02 −A 2 x. )
In view of the safe operation of the servo-hydraulic system actuators, a safety margin is left, i.e., based on the physical construction, the hydraulic cylinder can move up and down about the center position by a maximum of 12 cm, and the command signal amplitude is 10 cm. Centimeters or less and ensure that h 1 , h 2 , and h 3 are bounded, i.e. each three positive numbers
Figure 2023517140000265
there is
Figure 2023517140000266
becomes,
Figure 2023517140000267
and
Figure 2023517140000268
is bounded, m(t) is the unknown slope of the dead zone of the spool, and m(t) is bounded (i.e. ∃M r >0,0<m(t)<M r ). Considering, from the extremum theorem, we know that a 3 is bounded,
From (41) and (42),
Figure 2023517140000269
to get
Lyapunov function
Figure 2023517140000270
and select
(43) and (44) satisfy the form of the inequality lemma, and as can be seen from the inequality lemma, a 3 is bounded, i.e. ∃A 3 >0 and |a 3 |<A 3 ,
Figure 2023517140000271
Because
Figure 2023517140000272
and
Figure 2023517140000273
is bounded, and
Figure 2023517140000274
so that for ∀t∈[0, τ max ), we ensure that the control error z 3 is bounded and the main power u s is bounded, and from the equation of state formula we see , the state variables x3 and x4 are bounded, and both are smaller than the pressure value of the hydraulic source Ps ,
S3-2-4: based on formula (18)
Figure 2023517140000275
the derivative with respect to time of
Figure 2023517140000276
looking for
Figure 2023517140000277

(however,
Figure 2023517140000278
and there is one positive constant r M4 such that 0>r 4 >r M4 . )
moreover,
Figure 2023517140000279
to get
Figure 2023517140000280
year,
Figure 2023517140000281
and
Figure 2023517140000282
is bounded and k 4 is an optional parameter when designing the controller, a 4 is bounded and
From (46) and (47),
Figure 2023517140000283
to get
Lyapunov function
Figure 2023517140000284
and select
(48) and (49) satisfy the form of the inequality lemma, and as can be seen from the inequality lemma, a 4 is bounded, i.e. ∃A 4 >0 and |a 4 |<A 4 ,
Figure 2023517140000285
Because
Figure 2023517140000286
and
Figure 2023517140000287
is bounded, and
Figure 2023517140000288
5. The asymmetric servohydraulic position of claim 4, ensuring that the control error z4 is bounded and the state variable x5 is bounded, for A low-complexity control method for tracking systems.
前記不等式補助定理の式は、
Figure 2023517140000289
であり、
aが有界であり(すなわち∃A>0,|a|<A)、
Figure 2023517140000290
である場合、xが必ず有界であり、
Figure 2023517140000291
である、ことを特徴とする請求項6に記載の非対称サーボ油圧位置追跡システムのための低複雑度制御方法。
The formula for the inequality lemma is
Figure 2023517140000289
and
a is bounded (i.e. ∃A>0, |a|<A),
Figure 2023517140000290
, then x must be bounded, and
Figure 2023517140000291
A low-complexity control method for an asymmetric servohydraulic position tracking system according to claim 6, characterized in that:
初期値の提案を用いて、S3-2におけるτmax=+∞の場合、すべての閉ループ信号が有界であることが依然として正確であることを証明する前記プロセスは以下の通り、
Figure 2023517140000292
であるため、
Figure 2023517140000293
は、空でない集合であり、
Figure 2023517140000294
であり、
τmax<+∞を仮定すると、初期値の提案では1つの時定数t∈[0,τmax)があることが強調され、
Figure 2023517140000295
となり、矛盾を引き起こし、得られた仮定がτmax=+∞であり、いずれかのt∈[0,+∞)に対して、すべての閉ループシステム信号が有界であり、追跡誤差が指定された境界内で収束することを確保できる、ことを特徴とする請求項4に記載の非対称サーボ油圧位置追跡システムのための低複雑度制御方法。
Using the initial value proposal, the above process to prove that it is still correct that all closed-loop signals are bounded for τ max =+∞ in S3-2 is as follows:
Figure 2023517140000292
Because
Figure 2023517140000293
is a non-empty set and
Figure 2023517140000294
and
Assuming τ max <+∞, the initial value proposal emphasizes that there is one time constant t 1 ε[0, τ max ),
Figure 2023517140000295
, which causes a contradiction, and the assumption obtained is that τ max =+∞, that for any tε[0,+∞) all closed-loop system signals are bounded, and the tracking error is specified 5. A low complexity control method for an asymmetric servohydraulic position tracking system according to claim 4, characterized in that it can ensure convergence within a bounded boundary.
JP2022579147A 2020-12-07 2021-10-19 Low complexity control method for asymmetric servohydraulic position tracking system Active JP7448692B2 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN202011438457.2 2020-12-07
CN202011438457.2A CN112486021B (en) 2020-12-07 2020-12-07 Low-complexity control method for asymmetric servo hydraulic position tracking system
PCT/CN2021/124567 WO2022121507A1 (en) 2020-12-07 2021-10-19 Low-complexity control method for asymmetric servo hydraulic position tracking system

Publications (2)

Publication Number Publication Date
JP2023517140A true JP2023517140A (en) 2023-04-21
JP7448692B2 JP7448692B2 (en) 2024-03-12

Family

ID=74941242

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2022579147A Active JP7448692B2 (en) 2020-12-07 2021-10-19 Low complexity control method for asymmetric servohydraulic position tracking system

Country Status (3)

Country Link
JP (1) JP7448692B2 (en)
CN (1) CN112486021B (en)
WO (1) WO2022121507A1 (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112486021B (en) * 2020-12-07 2021-10-08 燕山大学 Low-complexity control method for asymmetric servo hydraulic position tracking system
CN113942354B (en) * 2021-10-09 2023-09-22 燕山大学 State feedback control method for active suspension of whole vehicle with electro-hydraulic actuator
CN115629543B (en) * 2022-10-21 2023-04-07 东北林业大学 Interventional MDF continuous flat pressing three-branch decision cooperative control method, system and storage medium
CN116300420B (en) * 2022-12-13 2024-03-15 南方海洋科学与工程广东省实验室(湛江) Servo control method, system, device, terminal equipment and storage medium
CN115903748A (en) * 2022-12-14 2023-04-04 燕山大学 Electro-hydraulic servo displacement tracking test device based on reference model and control method
CN116300565A (en) * 2022-12-27 2023-06-23 南京理工大学 Preset performance control method with boundary crossing preventing function
CN117389156B (en) * 2023-12-09 2024-03-22 中信重工机械股份有限公司 Hydraulic mechanical arm self-adaptive integral robust control method based on friction compensation

Family Cites Families (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000000415A (en) * 1998-06-16 2000-01-07 Mitsubishi Heavy Ind Ltd Seal structure of ceramic filter
JP4247699B2 (en) 2001-02-28 2009-04-02 富士電機システムズ株式会社 Feedback control method for nonlinear control object
JP2004192594A (en) * 2002-12-09 2004-07-08 Takara Sangyo:Kk Servo cylinder feeder
JP3896081B2 (en) * 2003-01-06 2007-03-22 新日本製鐵株式会社 Parameter identification method for hydraulic servo system
AT502348B1 (en) 2005-08-17 2008-09-15 Voest Alpine Ind Anlagen CONTROL METHOD AND REGULATOR FOR A MECHANICAL-HYDRAULIC SYSTEM WITH A MECHANICAL FREEDOM DEGREE PER HYDRAULIC ACTUATOR
DK1924669T3 (en) * 2005-09-16 2012-10-22 Halliburton Energy Serv Inc Polymer blends for crosslinked fluids
CN101216695B (en) * 2007-12-26 2010-06-02 燕山大学 Pressure setting, prediction and self learning method for temper rolling
NL2003800C2 (en) * 2009-11-13 2011-05-16 Baggerwerken Decloedt & Zn N V Device for dredging soil material under water.
RU2559998C2 (en) * 2010-04-12 2015-08-20 Марк Энтони КУЧЕЛ Method for soil treatment
CN104635490B (en) 2014-12-15 2017-06-27 南京理工大学 A kind of output feedback ontrol method of asymmetric servo cylinder positional servosystem
CN104698844B (en) * 2015-02-09 2017-04-19 南京理工大学 Uncertainty compensatory sliding-mode control method of hydraulic position servo system
CN105068426B (en) * 2015-08-24 2018-03-09 南京理工大学 The continuous sliding-mode control of electro-hydraulic position servo system based on interference compensation
CN105697463B (en) * 2016-03-22 2017-08-25 西安理工大学 A kind of Hydraulic Position Servo exports feedback adaptive control method
JP6487872B2 (en) * 2016-03-30 2019-03-20 日立建機株式会社 Drive control device for work machine
CN106406088A (en) * 2016-08-30 2017-02-15 上海交通大学 Switching MMSLA-based hydraulic servo system control method and control system
CN106640809B (en) * 2016-11-07 2018-04-17 同济大学 Variable displacement and variable power adjusts load matched electro-hydraulic position tracking and controlling method
JP6544408B2 (en) * 2017-03-03 2019-07-17 日亜化学工業株式会社 Optical component and method of manufacturing optical component
CN109901387B (en) * 2017-12-11 2023-01-10 上海航空电器有限公司 Self-adaptive flight trajectory prediction method for aircraft automatic near-ground collision avoidance system
CN108995495B (en) * 2018-08-09 2020-04-03 燕山大学 Anti-saturation self-adaptive control method and system for nonlinear active suspension
CN110081046B (en) * 2019-05-27 2020-07-31 电子科技大学 Multi-hydraulic servo actuator tracking synchronous control method based on backstepping control
CN110262232B (en) * 2019-05-30 2021-12-17 中国矿业大学 Control method for high tracking precision of underwater missile launching manhole cover opening device
CN110397634B (en) * 2019-07-17 2020-07-07 太原理工大学 Low-energy-consumption high-dynamic pump valve combined position servo system and control method thereof
CN110928182B (en) * 2019-11-05 2022-12-13 南京理工大学 Robust self-adaptive repetitive control method of hydraulic servo system based on state estimation
CN110703608B (en) * 2019-11-18 2021-08-06 南京工业大学 Intelligent motion control method for hydraulic servo actuator
CN110794687B (en) * 2019-12-02 2022-05-31 安徽工业大学 Electro-hydraulic servo system self-adaptive state constraint control method based on interference compensation
CN111577711B (en) * 2020-04-30 2022-08-16 南京理工大学 Active-disturbance-rejection robust control method for double-rod hydraulic cylinder position servo system
CN112486021B (en) * 2020-12-07 2021-10-08 燕山大学 Low-complexity control method for asymmetric servo hydraulic position tracking system

Also Published As

Publication number Publication date
CN112486021B (en) 2021-10-08
CN112486021A (en) 2021-03-12
JP7448692B2 (en) 2024-03-12
WO2022121507A1 (en) 2022-06-16

Similar Documents

Publication Publication Date Title
JP2023517140A (en) A Low-Complexity Control Method for Asymmetric Servohydraulic Position Tracking Systems
Deng et al. Robust adaptive precision motion control of hydraulic actuators with valve dead-zone compensation
Feng et al. Identification and compensation of non-linear friction for a electro-hydraulic system
CN105700347B (en) A kind of hydraulic motor default capabilities tracking and controlling method containing Hysteresis compensation
CN110928182B (en) Robust self-adaptive repetitive control method of hydraulic servo system based on state estimation
CN112415891B (en) Adaptive output feedback asymptotic control method for electro-hydraulic servo system
CN104345638B (en) A kind of active disturbance rejection self-adaptation control method of hydraulic motor positional servosystem
CN108869420A (en) A kind of adaptive backstepping control method of electrohydraulic servo system filtered based on instruction
CN110703608B (en) Intelligent motion control method for hydraulic servo actuator
CN108415249B (en) Electro-hydraulic servo system adaptive robust control method based on low-frequency learning
Luo et al. Extended-state-observer-based output feedback adaptive control of hydraulic system with continuous friction compensation
CN110308651B (en) Electro-hydraulic servo system all-state constraint control method based on extended state observer
CN110647111B (en) Output-discreteness-considered non-linear active disturbance rejection control method for electro-hydraulic servo system
Niu et al. Neural network-based finite-time command-filtered adaptive backstepping control of electro-hydraulic servo system with a three-stage valve
CN114967466A (en) Electro-hydraulic servo system preset performance output feedback control method based on interference compensation and all-state constraint
CN114185276B (en) Nonlinear strict system output feedback control method based on multidimensional Taylor network
CN115236974A (en) Composite anti-interference controller and control parameter optimization method thereof
CN114077192B (en) Nonlinear system output feedback intelligent control method
Wang et al. Active fault-tolerant control for the dual-valve hydraulic system with unknown dead-zone
CN111781836A (en) Self-adaptive asymptotic control method for hydraulic pressure preset performance
Carneiro et al. VSC approach angle based boundary layer thickness: a new variation law and its stability proof
CN114995127A (en) Self-adaptive integral control method of uncertain hydraulic position servo system
Nguyen et al. Synthesizing an Adaptive Controller to Enhance Movement Quality of a Flight Control Actuator under External Load
CN109426140A (en) Load simulator parameter Impact analysis method based on SIMULINK
CN117823496A (en) Robust output feedback composite control method for pilot type proportional servo valve

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220906

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230831

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230904

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231204

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20240205

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240229

R150 Certificate of patent or registration of utility model

Ref document number: 7448692

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150