劉 濤，程相文，高軍霞

（1．河北聯(lián)合大學(xué) 機械工程學(xué)院，河北 唐山063009；2．唐山盾石電氣有限責(zé)任公司，河北 唐山063000；3．唐山學(xué)院 機電工程系，河北 唐山063000）

0 引言

板坯連鑄技術(shù)的核心部件是結(jié)晶器，盡管各種板坯連鑄機所用的結(jié)晶器構(gòu)造和形狀不同，但在振動形式上卻越來越趨于一致，即液壓振動，屬高頻率、小振幅的非正弦曲線振動方式。［1］采用計算機控制的電液伺服結(jié)晶器振動裝置可以很方便地產(chǎn)生各種振動，實現(xiàn)控制過程監(jiān)督、實時顯示，并根據(jù)拉坯速度實時修改振動參數(shù)［2］，提高連鑄坯質(zhì)量和金屬收得率，從而實現(xiàn)連鑄過程的自動化。

某鋼鐵公司熱軋帶鋼板坯連鑄機結(jié)晶器系統(tǒng)原理如圖1所示。液壓振動系統(tǒng)主要由伺服閥、液壓缸等組成，液壓缸的管路以并聯(lián)方式連接到液壓站上，液壓站用于提供液壓動力，管路上還配有控制伺服閥的截止閥，由壓力傳感器檢測油壓，伺服閥中有位移傳感器檢測閥芯的位置。另外，結(jié)晶器左、右兩邊分別裝設(shè)位移傳感器，檢測結(jié)晶器的實際位移。液壓缸的動作以理想的頻率、行程和波形造成結(jié)晶器的振動動作。

1 結(jié)晶器液壓伺服振動系統(tǒng)靜態(tài)設(shè)計

在設(shè)計中采用德馬克公司開發(fā)的非正弦振動曲線，描述為

式中：A為結(jié)晶器振動振幅；f為結(jié)晶器振動頻率；t為振動時間；h為非正弦振動曲線偏斜率。

參考實際生產(chǎn)工況，確定結(jié)晶器負載范圍為0～4 500kg。液壓伺服缸端非正弦振動參數(shù)設(shè)置為：A＝6mm，f＝4Hz，h＝0．4。結(jié)晶器預(yù)設(shè)非正弦振動曲線如圖2所示。

液壓缸加速度變化規(guī)律為

非正弦振動加速度曲線如圖3所示。由圖可看出各條曲線連續(xù)而光滑，加速度控制在5．2m／s2以內(nèi)，由此產(chǎn)生的沖擊在系統(tǒng)允許范圍內(nèi)。

在連鑄結(jié)晶器振動過程中，結(jié)晶器端總負載組成分為以下幾部分：結(jié)晶器裝置自身重量、負載慣性力以及振動過程中結(jié)晶器內(nèi)壁與鑄坯間的摩擦力，包括干摩擦和粘性摩擦，對結(jié)晶器振動端到液壓缸端的負載折算，總負載FL的表達式為

圖1 結(jié)晶器液壓振動系統(tǒng)原理圖

圖2 非正弦與正弦位置曲線對比

圖3 非正弦振動加速度曲線

式中Ff為結(jié)晶器內(nèi)壁與鑄坯表面的干摩擦力，這里取Ff＝10 000N；m為負載質(zhì)量；Bp為粘性阻尼系數(shù)；g為重力加速度。

根據(jù)對液壓缸一個行程周期內(nèi)的50個點采集負載力F（N）及負載速度v（m／s），利用 MATLAB軟件繪制出液壓缸v－F負載軌跡，如圖4所示。最大功率點功率Nmax＝15．1kW，此時負載力F＊L＝37 397．79N，負載振動速度V＊L＝0．404 6m／s。

以最大功率點為參照選擇系統(tǒng)壓力Ps＝16MPa。為了使液壓動力元件能與負載做到最佳匹配，取液壓缸輸出最大功率時的負載壓力為PL＊＝23Ps，則液壓缸面積為

圖4 負載軌跡

最大功率時所需的供油流量q＊L＝ApV＊L＝83．4L／min。

考慮泄漏等影響，將系統(tǒng)供油流量qL增大15%左右，取qL＝1．15q＊L＝96L／min。

根據(jù)閥口總壓降和所需供油流量，確定選擇MOOG公司生產(chǎn)的電液伺服閥，型號規(guī)格為D661G45H。該閥在閥壓降為1MPa時額定流量為45L／min，最大工作壓力28MPa，額定電流20mA。

2 結(jié)晶器液壓振動系統(tǒng)動態(tài)設(shè)計

伺服閥控液壓缸的動態(tài)工作過程可用下面3個方程描述［3］：

式中：Kq為伺服閥流量增益；Kc為伺服閥流量系數(shù)；mt為負載質(zhì)量；Bp為粘性阻尼系數(shù)；K為負載彈簧剛度；Ctp為液壓缸總泄漏系數(shù)；Fg為外加干擾負載力。

由上式3個基本方程得到閥芯位移xv到液壓缸輸出位移xp的傳遞函數(shù)為

閥控缸輸出位移xp對Fg的傳遞函數(shù)為

式中Kce＝Kc＋Ctp，指的是考慮泄漏因素后的總壓力流量系數(shù)。

系統(tǒng)有效體積彈性模量βe＝700MPa；液壓缸有效面積34．34×10－4m2；液壓缸活塞行程L＝25mm，則液壓缸最大容積Vt＝LAp＝8．59×10－5m3；動力元件的阻尼比一般根據(jù)實測確定，這里取ζh＝0．2；忽略液壓缸泄漏，總流量－壓力系數(shù)Kce＝Kc，由閥的零位泄漏決定，取Kce＝3．73×10－12m3／（s·Pa）。

當(dāng)電液伺服閥的頻寬與液壓動力元件的固有頻率相近時，電液伺服閥的傳遞函數(shù)可以看成是二階環(huán)節(jié)。設(shè)計的液壓動力元件的固有頻率

則液壓動力元件的固有頻率ω0＝292rad／s。

查MOOGD661G45H伺服閥樣本得出，伺服閥的固有頻率ωsv＝571．2rad／s，阻尼比ζsv＝0．7。由于該閥在閥壓降1MPa時額定流量為45L／min，因此其最大空載流量為

閥的額定電流為in＝20mA＝0．02A，則伺服閥的增益為。

電液伺服閥的傳遞函數(shù)為

伺服放大器、位移傳感器的響應(yīng)很快，因此一般忽略它們對系統(tǒng)的動態(tài)影響，看成比例環(huán)節(jié)，位移傳感器的傳遞函數(shù)可表示為；伺服放大器的傳遞函數(shù)可以表示為。

結(jié)合以上分析，得出系統(tǒng)傳遞函數(shù)方塊示意圖如圖5所示。

圖5 系統(tǒng)傳遞函數(shù)方塊圖

3 仿真分析

為了對系統(tǒng)輸出性能進行預(yù)估，利用如圖6所示的系統(tǒng)方框圖，采用MATLAB仿真軟件對系統(tǒng)進行仿真［4］。仿真時代入有關(guān)數(shù)據(jù)，可直接得出系統(tǒng)在時域的輸入、輸出曲線，如圖7所示。系統(tǒng)的輸出基本能跟隨輸入信號的變化，但是存在幅值超調(diào)，振幅最大誤差值為1×10－5m，位置控制精度達不到設(shè)計要求。

圖6 仿真模型圖

圖7 輸入、輸出曲線

利用Simulink模塊強大的功能，構(gòu)建加入PID控制器后的仿真模型如圖8所示。

不同PID控制參數(shù)下的伺服控制系統(tǒng)輸入輸出性能仿真曲線如圖9－11所示。通過研究仿真曲線，以求確定最佳PID控制參數(shù)。

根據(jù)圖9－11，綜合考慮系統(tǒng)的響應(yīng)速度與位置控制精度，確定最佳PID控制參數(shù)為Kp＝1，Ki＝1，Kd＝0。

圖8 加入PID控制器后的仿真模型（非正弦輸入）

圖9 仿真曲線（Kp＝0．6，Ki＝0．2，Kd＝0）

圖10 仿真曲線（Kp＝0．6，Ki＝1，Kd＝0）

圖11 仿真曲線（Kp＝1，Ki＝1，Kd＝0）

4 結(jié)論

本文對結(jié)晶器的液壓伺服非正弦振動系統(tǒng)的設(shè)計、仿真優(yōu)化、控制等進行了詳細地闡述說明，并利用 MATLAB／Simulink軟件確定了最佳PID控制參數(shù)。仿真結(jié)果表明，在非正弦振動控制系統(tǒng)中增加PID控制器可大大改善系統(tǒng)的控制精度，本文可為類似帶彈性負載的液壓伺服非正弦振動設(shè)計和PID控制性能仿真分析提供借鑒和參考。

［1］ 孔德才．結(jié)晶器振動系統(tǒng)在連鑄中的應(yīng)用與研究［J］．連鑄，2007（3）：23－25．

［2］ 吳曉明．板坯連鑄結(jié)晶器液壓伺服振動裝置控制策略的研究［D］．秦皇島：燕山大學(xué)，2007．

［3］ 高軍霞，楊國權(quán)．基于 MATLAB／SIMULINK的篦式冷卻機液壓系統(tǒng)改造［J］．液壓與氣動，2011（9）：101－103．

［4］ 高軍霞，楊國權(quán)．房屋液壓減震裝置的設(shè)計與研究［J］．機床與液壓，2011，30（12）：54－56．