JP3896081B2 - 油圧サーボ系のパラメータ同定方法 - Google Patents

油圧サーボ系のパラメータ同定方法 Download PDF

Info

Publication number
JP3896081B2
JP3896081B2 JP2003000081A JP2003000081A JP3896081B2 JP 3896081 B2 JP3896081 B2 JP 3896081B2 JP 2003000081 A JP2003000081 A JP 2003000081A JP 2003000081 A JP2003000081 A JP 2003000081A JP 3896081 B2 JP3896081 B2 JP 3896081B2
Authority
JP
Japan
Prior art keywords
equation
servo system
output
hydraulic servo
identified
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.)
Expired - Fee Related
Application number
JP2003000081A
Other languages
English (en)
Other versions
JP2004209523A (ja
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.)
Nippon Steel Corp
Original Assignee
Nippon Steel Corp
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 Nippon Steel Corp filed Critical Nippon Steel Corp
Priority to JP2003000081A priority Critical patent/JP3896081B2/ja
Publication of JP2004209523A publication Critical patent/JP2004209523A/ja
Application granted granted Critical
Publication of JP3896081B2 publication Critical patent/JP3896081B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Description

【0001】
【発明の属する技術分野】
本発明は、例えば圧延設備の油圧圧下装置や油圧自動板位置調整機構において、シリンダの位置をサーボ制御するために用いられている油圧サーボ系のパラメータ同定方法に関するものである。
【0002】
【従来の技術】
【特許文献1】
特開平3-174910号公報
【特許文献2】
特公平6−2290号公報
【0003】
油圧サーボ系は、切替弁のスプールをサーボ機構制御部からの弁開度指令値(基準値)により移動させて弁開度を制御し、シリンダに流入する油量を制御している。サーボバルブの特性を時定数Tp、ゲインKpで表すと、系全体のブロック線図は、制御のための比例ゲインをKcとして、図1のように表すことができる。ここで入力u(例えばシリンダ位置指令)から出力y(同様にシリンダ位置実績)への伝達関数を一次遅れで表現した場合、数1の(1)で表される。同様に二次遅れで表現した場合、数1の(2)で表される。ただし数1中のTは時定数、K1,K2は定常ゲイン、ωは共振周波数、ζは減衰係数である。これらのパラメータT,K1,K2,ω、ζを同定できれば、油圧サーボ系の挙動を正確に把握することができ、診断や管理を的確に行うことができる。
【0004】
【数1】
Figure 0003896081
【0005】
油圧サーボ系の診断方法としては、ステップ状の入力信号を与え、それに対する図2のようなステップ応答を取る方法が基本的な方法である。しかし稼動中の圧延設備などに対してステップ入力を与えることは困難であり、設備を停止してオフラインでしか診断を行うことができない。そこで上記の特許文献1では、圧延設備のワークロールに被圧延材が噛み込まれる際、もしくはワークロールから被圧延材が抜け出す際にサーボ系に入力されるインパルス状の信号と、それに対するサーボ系の出力信号をサンプリングし、これに基づいて減衰係数ζを同定して診断する方法を取っている。この方法は、インパルス信号が理論上すべての周波数成分を持つことから、その疑似信号を利用し、オンラインでパラメータの同定をする方法である。
【0006】
しかし減衰係数ζが許容範囲にあっても、必ずしも油圧サーボ系が正常とはいえない場合がある。すなわち、図3は油圧サーボ系の減衰係数ζとステップ応答との関係を示したグラフであり、▲3▼と▲4▼はともに減衰係数ζの値自体は正常とされる0.701に近い許容範囲内にあるにもかかわらず、▲4▼は応答が遅く、設備によっては要求される仕様を満足しないため、正常とはいえない。従って、特許文献1に開示された減衰係数ζだけを同定する方法では、不十分な面がある。
【0007】
また特許文献2には、稼動中の入出力データにフーリエ変換を施してボード線図を作成し、ゲイン曲線の所定周波数区間を近似することにより、油圧サーボ系の減衰係数ζを同定する方法が開示されている。この方法も設備を停止させる必要はないが、特許文献1の方法と同様にパラメータとして減衰係数ζのみに着目しているため、油圧サーボ系の診断を正しく行う上では不十分である。
【0008】
【発明が解決しようとする課題】
本発明は上記した従来の問題点を解決して、油圧サーボ系の減衰係数ζのみならず、ω、Tなどの他のパラメータをも、設備を停止することなく正確に同定することができる油圧サーボ系のパラメータ同定方法を提供するためになされたものである。
【0009】
【課題を解決するための手段】
上記の課題を解決するためになされた本発明は、同定対象となる油圧サーボ系の挙動を、未知パラメータを含まない状態遷移方程式と、同定すべき未知パラメータを含む出力方程式からなる状態方程式で表現し、この状態方程式と同定対象の入力及び出力とを用いて推定器を構成し、この推定器の推定出力と同定対象の出力との二乗積分誤差を最小化するように未知パラメータを決定することを特徴とするものである。また、本発明は、圧延設備などを止めることなく稼動中のまま、入力信号に特別な特徴を要求することなく、油圧サーボ系の入出力データを用い、ζ、ω、Tなどの未知パラメータを正確に同定することができる。
【0010】
このため異常診断、傾向管理の簡便化、迅速化、精度向上などを実現でき、設備稼働率向上の効果も得ることができる。
【0011】
【発明の実施の形態】
以下に本発明を具体的に説明する。
先ず本発明の概要を図4に示す。本発明においては予め同定モデルを定め、これを漸近安定な状態遷移方程式と出力方程式とからなる状態方程式で表現する。後記する数学的手法により状態遷移方程式は同定すべき未知パラメータを含まない形とし、出力方程式に同定すべき未知パラメータをすべて含ませる。そして同定対象となる油圧サーボ系の入出力信号u,yを用いて、図3に示す推定器を構成する。この推定器は同定対象となる油圧サーボ系の出力を推定するものである。
【0012】
そして同一の入力信号uに対する同定対象となる油圧サーボ系の出力yと、上記した推定器の出力ypとの誤差eを演算し、図5にハッチングで示す面積、すなわちある区間の誤差eの二乗が最小となるように最小二乗法により出力方程式に含ませたζ、ω、Tなどの未知パラメータを同定する。この方法によれば、入力信号uとして特別な特徴は要求されず、単一周波数の入力信号uを用いることもできる。これに対して特許文献1の方法ではインパルス状の入力信号を要求し、また特許文献2の方法ではボード線図から同定を行うため、幅広い周波数を持つ入力信号が要求されるため、この点においても本発明法は従来法と明確に相違している。
【0013】
以下に、本発明において利用した数学的手法をより詳細に説明する。
前記したように、本発明においては予め同定モデルを定め、それを状態方程式で表現する。図1に示したようにサーボ系の特性は時定数TpでゲインK1、Kcで表される。また入力uから出力yへの伝達関数の近似は、一次遅れの場合に数1(1)、二次遅れの場合に数1(2)によって表される。まず同定モデルを一次遅れ系とした場合、xを状態とし、uを入力、yを出力とすると、状態方程式は数2の(3),(4)式となる。ここで同定すべき未知パラメータはT、K1である。後に同定モデルを二次遅れとした場合も説明するが、同定モデルはこの二つに限定されるものではなく、変更することも可能である。
【0014】
【数2】
Figure 0003896081
【0015】
この時点で、同定すべき未知パラメータがシステム行列、入力行列に含まれている。それに対し、出力行列には未知パラメータが含まれていない。ただしこのモデルの出力yは未知パラメータT、K1に依存するため推定値となる。この値をサーボ系の出力と区別するためypとし、数2の(4)を数3の(5)に書き改めておく。この状態方程式が推定器の内容である。
【0016】
【数3】
Figure 0003896081
【0017】
ここで数2の(3)式の変形を行い、数4の(6)式のようにシステム行列(ここではスカラ)を未知パラメータを含まない適当な負の実数係数部分と、未知パラメータを含む部分とに分ける。式中のaは適当な負の実数である。
【0018】
【数4】
Figure 0003896081
【0019】
そして状態変数の一部を、油圧サーボ系の出力yに置き換えられることに注目して、更に変形する。置き換えは、未知パラメータを含んだ部分のみ行う。ここではx=yと置き換え、数5の(7)式を得る。
【0020】
【数5】
Figure 0003896081
【0021】
この時点で、このモデルに対する入力は油圧サーボ系の入力uと出力yになる。またこのモデルの出力方程式は数3の(5)と同じものである。また未知パラメータは入力行列(ここではスカラ)にのみ含まれる形になる。
【0022】
さらに、状態方程式の未知パラメータをp=1/T、q=K1/Tと置き換え、状態方程式を入力行列に関して、未知パラメータpに対応する部分およびqに対応する部分、未知パラメータに対応しない部分に細分化し、新たな状態変数をm1,m2,m3を導入することにより数6の状態方程式(8),(9)を得る(d1、d2、bは式の変形に伴う係数である)。この状態遷移方程式(8)には未知パラメータが含まれず、出力方程式に未知パラメータが含まれる形となっている。また、入力は装置の入出力u,yであり、出力はxとなる。状態遷移方程式に未知パラメータが含まれていないので、装置の入出力u,yのみから状態m1,m2,m3の応答を求めることができる。
【0023】
【数6】
Figure 0003896081
【0024】
推定器の出力ypは式(5),(9)により得られる。更に未知パラメータを含む変数p,qは状態m1,m2,m3に線形に係っているのみである。従って、以下に述べるように最小二乗法を応用することで未知パラメータを含む変数p,qを決定することができる。まず適当な時間m1,m2,m3,y,ypをサンプリングし、それらをm1 i,m2 i,m3 i,yi,yp i(i=1、…、N)とする。そして図5に示すように、yi,とyp iの二重積分誤差Σ(yi−yp i2を最小化するように未知パラメータを決定することができる。具体的には、下記のとおりである。
【0025】
式(5),(9)よりypを書き下すと、数7の(10)となる。そしてy−ypは(11)で表されるから、(12)のようにY、M1,M2,M3を定め、(13)のようにW,Z,θを置けば、最小二乗法によりθの推定値は(ZTZ)-1ZWで求められる。そしてT=1/p、K1=q/pとして未知パラメータをすべて決定することができる。
【0026】
【数7】
Figure 0003896081
【0027】
同定モデルを二次遅れとした場合、状態方程式は数8の(14),(15)と表現できる。そして、一次遅れで使用した方法と同じ方法を適用することによって状態方程式は(16),(17)のようになる。式中でAは適当な2×2の安定行列、D1,D2,D3,Bはそれぞれ式変形に伴う2×1の係数行列、m1,m2,m3,m4はそれぞれ2×1の新たな状態変数ベクトル、Iは2×2の単位行列、p,q,rはそれぞれp=ω2、q=ζω、r=K2ω2と置き換えたものである。(16)から明らかなように、状態遷移方程式に未知パラメータは含まれていない。よって、入出力u,yのみからm1,m2,m3,m4の応答を求めることができるので、一次遅れの場合と同様に最小二乗法により未知パラメータω、ζ、K2を決定することができる。
【0028】
【数8】
Figure 0003896081
【0029】
このようにして求められたT、ω、ζ、K1、K2に基づき、油圧サーボ系の診断を行うことができる。また本発明の方法によれば、入出力u,yにある程度のノイズが乗っていても、同定精度がそれほど悪くならないことを確認している。
【0030】
【発明の効果】
以上に説明したように、本発明の油圧サーボ系のパラメータ同定方法によれば、圧延設備などを止めることなく、稼動中の油圧サーボ系の入出力データを用いてζ、ω、K2などの未知パラメータを正確に同定することができる。このため従来のように特別な入力波形を必要としない。未知パラメータを正確に同定することができるので、油圧サーボ系の異常診断、傾向管理の簡便化、迅速化、精度向上などを実現することができる。さらに設備を停止する必要がないため、稼働率を低下させることもない。
【図面の簡単な説明】
【図1】油圧サーボ系のブロック線図である。
【図2】油圧サーボ系のステップ応答図である。
【図3】油圧サーボ系の減衰係数ζとステップ応答との関係を示したグラフである。
【図4】本発明の概要を示すブロック図である。
【図5】油圧サーボ系の出力yと、推定器の出力ypとの誤差eを示すグラフである。

Claims (1)

  1. 同定対象となる油圧サーボ系の挙動を、未知パラメータを含まない状態遷移方程式と、同定すべき未知パラメータを含む出力方程式からなる状態方程式で表現し、この状態方程式と同定対象の入力及び出力とを用いて推定器を構成し、この推定器の推定出力と同定対象の出力との二乗積分誤差を最小化するように未知パラメータを決定することを特徴とする油圧サーボ系のパラメータ同定方法。
JP2003000081A 2003-01-06 2003-01-06 油圧サーボ系のパラメータ同定方法 Expired - Fee Related JP3896081B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2003000081A JP3896081B2 (ja) 2003-01-06 2003-01-06 油圧サーボ系のパラメータ同定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2003000081A JP3896081B2 (ja) 2003-01-06 2003-01-06 油圧サーボ系のパラメータ同定方法

Publications (2)

Publication Number Publication Date
JP2004209523A JP2004209523A (ja) 2004-07-29
JP3896081B2 true JP3896081B2 (ja) 2007-03-22

Family

ID=32818517

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2003000081A Expired - Fee Related JP3896081B2 (ja) 2003-01-06 2003-01-06 油圧サーボ系のパラメータ同定方法

Country Status (1)

Country Link
JP (1) JP3896081B2 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102174967A (zh) * 2011-03-01 2011-09-07 江苏科技大学 电液角位移伺服系统中执行机构参数识别装置及识别方法
CN102563181A (zh) * 2012-01-10 2012-07-11 江苏科技大学 电液线位移伺服系统中伺服控制器控制参数的确定方法
CN102563182A (zh) * 2012-01-10 2012-07-11 江苏科技大学 电液线位移伺服系统中伺服控制器控制参数的调整方法

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008102714A (ja) * 2006-10-18 2008-05-01 Tokyo Univ Of Agriculture & Technology 多軸工作機械の送り駆動系用最適設計支援装置およびその装置用プログラム
CN102252126B (zh) * 2011-04-07 2013-04-24 江苏科技大学 电液角位移伺服系统中伺服对象参数识别方法
CN102252127B (zh) * 2011-04-07 2012-12-05 江苏科技大学 电液线速度伺服系统中伺服对象参数识别方法
CN102331797B (zh) * 2011-04-07 2012-08-29 江苏科技大学 电液角速度伺服系统中伺服对象参数识别方法
CN102128304B (zh) * 2011-04-07 2012-07-18 江苏科技大学 电液线位移伺服系统中伺服对象参数识别方法
CN112486021B (zh) * 2020-12-07 2021-10-08 燕山大学 一种针对非对称伺服液压位置跟踪系统的低复杂控制方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102174967A (zh) * 2011-03-01 2011-09-07 江苏科技大学 电液角位移伺服系统中执行机构参数识别装置及识别方法
CN102174967B (zh) * 2011-03-01 2012-10-24 江苏科技大学 电液角位移伺服系统中执行机构参数识别装置及识别方法
CN102563181A (zh) * 2012-01-10 2012-07-11 江苏科技大学 电液线位移伺服系统中伺服控制器控制参数的确定方法
CN102563182A (zh) * 2012-01-10 2012-07-11 江苏科技大学 电液线位移伺服系统中伺服控制器控制参数的调整方法

Also Published As

Publication number Publication date
JP2004209523A (ja) 2004-07-29

Similar Documents

Publication Publication Date Title
DE102004058238B4 (de) Adaptive, multivariable Prozesssteuerung, die Modellschaltung und Attribut-Interpolation nutzt
US6459939B1 (en) Performance assessment of model predictive controllers
EP2108139B1 (de) Verfahren zur rechnergestützten regelung und/oder steuerung eines technischen systems, insbesondere einer gasturbine
EP1715352B1 (de) Verfahren und vorrichtung zur fehlerdiagnose mechatronischer systeme
JP3896081B2 (ja) 油圧サーボ系のパラメータ同定方法
EP2987039B1 (de) Verfahren und vorrichtung zur co-simulation von zwei teilsystemen
DE4213795A1 (de) Motor-servosystem-regelung
JP4396541B2 (ja) 電動機制御装置の制御パラメータ感度解析装置および電動機制御装置の制御パラメータ設定方法
KR102114068B1 (ko) 계산 토크 방식 제어기, 이의 파라미터 결정 및 성능 평가 방법
EP0401847A2 (en) Method and apparatus for controlling engine in engine characteristics test
DE10020448A1 (de) Verfahren und Vorrichtung zur Optimierung des Betriebs eines Verbrennungsmotors
DE112013002374T5 (de) Kontroller und verfahren zum steuerneines betriebs eines systems
JP2006227793A5 (ja)
DE102020205031A1 (de) Werkzeugmaschinensteuerung und Verfahren zur kennfeldbasierten Fehlerkompensation an einer Werkzeugmaschine
DE102016007404A1 (de) Servosteuervorrichtung mit Funktion zur Online-Optimierung der Regelungsverstärkung unter Verwendung einer Bewertungsfunktion
JP4908433B2 (ja) 制御パラメータ調整方法および制御パラメータ調整プログラム
DE102010025916A1 (de) Verfahren und Vorrichtung zur Ermittlung von Modellparametern zur Regelung eines Dampfkraftwerksblocks, Regeleinrichtung für einen Dampferzeuger und Computerprogrammprodukt
JP3935079B2 (ja) 油圧サーボ系を含む装置の異常診断・傾向管理方法
DE102018003244A1 (de) Numerische Steuerung
DE602006000731T2 (de) Autoadaptives Einstellungsgerät zur Positionssteuerung von Aktuatoren in einem Antriebssystem mittels Gleitmodusverfahren und entsprechendes Betriebsverfahren
JP2005275588A (ja) 電動機制御装置制御パラメータ感度解析装置
DE102016222732A1 (de) Verfahren und Vorrichtung zum Durchführen einer Positionsregelung für eine Stellgebereinheit
CN113641181A (zh) 基于在线性能评估飞行器增益自适应姿态控制方法和系统
JPH0651805A (ja) プラントの適応制御方法およびそれを実現する装置
CN110538881B (zh) 一种基于改进型内模控制器的热连轧厚度控制方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20050913

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20061117

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: 20061208

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20061215

R151 Written notification of patent or utility model registration

Ref document number: 3896081

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101222

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101222

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111222

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111222

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121222

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121222

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131222

Year of fee payment: 7

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131222

Year of fee payment: 7

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131222

Year of fee payment: 7

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees