陳 敏，李 盛，李新有，張 瑤

(1.中冶賽迪重慶信息技術(shù)有限公司，重慶 401122； 2.中冶賽迪技術(shù)研究中心有限公司，重慶 401100)

平整機的原理是通過軋輥對帶鋼施加一定的壓下量進行小壓下量的軋制變形，以改善產(chǎn)品結(jié)構(gòu)，提高帶鋼板形質(zhì)量，增強產(chǎn)品的市場競爭能力。為確保生產(chǎn)的穩(wěn)定性和產(chǎn)品質(zhì)量，對平整機的綜合性能，如軋制力控制精度以及壓下響應(yīng)時間等有明確的要求。系統(tǒng)的綜合性能是由機械系統(tǒng)、液壓系統(tǒng)、控制系統(tǒng)等交叉耦合作用的結(jié)果。傳統(tǒng)的單領(lǐng)域仿真往往從單學科的角度出發(fā)仿真，將涉及其他學科的子系統(tǒng)加以簡化甚至完全忽略。這種做法降低整個問題的復雜度，在實際應(yīng)用中也有一定效果，勢必將會造成系統(tǒng)整體性和相關(guān)性的丟失，從而導致仿真結(jié)果置信度的降低[1]，為此需要采用多領(lǐng)域的聯(lián)合仿真技術(shù)從系統(tǒng)角度進行研究與分析。

聯(lián)合仿真技術(shù)在國內(nèi)外已經(jīng)有較多的應(yīng)用，如楊勇等開發(fā)的多軸起重機聯(lián)合仿真平臺[2]；葉仁平等[3]、陳罡等[4]等開展的機器人手臂聯(lián)合仿真研究；胡驍檣等[5]進行的擺式車輛通過曲線研究；太原重工利用Amesim+Adams對挖掘機工作裝置和液壓系統(tǒng)進行了聯(lián)合仿真，模擬了大型礦用液壓挖掘機的作業(yè)過程，獲得了作業(yè)過程挖掘機的響應(yīng)特性，系統(tǒng)受力以及能量消耗；中煙機械集團公司利用MATLAB+Adams對超高速卷接機鼓輪系統(tǒng)進行了聯(lián)合仿真，研究了多軸鼓輪系統(tǒng)在控制系統(tǒng)控制下的同步性能。這些應(yīng)用都通過多體動力學軟件(如ADAMS、SIMPACK、RECDYN)構(gòu)建機械模型，忽略了機械系統(tǒng)內(nèi)部的非線性特性，存在較多的局限，難以滿足復雜裝備系統(tǒng)針對高精度、高響應(yīng)特性的仿真需求。

本文在MATLAB中建立了液壓與控制系統(tǒng)；在ABAQUS中建立了機械系統(tǒng)，考慮了接觸、局部大變形等非線性特性；進而實現(xiàn)了有限元模型和數(shù)學解析模型的聯(lián)合仿真。通過基于FMI標準的MATLAB與ABAQUS的聯(lián)合仿真技術(shù)，對軋鋼領(lǐng)域主機設(shè)備平整機進行分析模擬，獲得了軋制力控制精度、壓下響應(yīng)時間等系統(tǒng)綜合性能，并與現(xiàn)場測試數(shù)據(jù)進行了對比，驗證了方法的正確性。

1 液壓與控制系統(tǒng)

在本聯(lián)合仿真中，液壓壓下系統(tǒng)及控制算法在Matlab中建模。平整機AGC液壓伺服系統(tǒng)主要包括AGC壓下回路，AGC背壓回路、上支撐輥平衡回路、液壓泵站回路等，應(yīng)用Matlab/Simulink中的Simscape、SimMechanics、SimHydraulics等模塊，完成伺服控制回路中各環(huán)節(jié)和負載的物理模型描述，建立了液壓伺服系統(tǒng)仿真模型，并結(jié)合平整機AGC的兩種實際工作模式狀態(tài)：即位置控制模式(換輥或軋輥非接觸狀態(tài))及力控制模式(軋輥接觸狀態(tài))，進行了開環(huán)空載、位置閉環(huán)空載、力閉環(huán)帶載三種條件下的仿真，仿真結(jié)果均滿足設(shè)備的實際保證值，驗證了伺服系統(tǒng)的模型正確性。其中力閉環(huán)仿真模型及仿真結(jié)果如圖1、圖2所示。

圖1 力閉環(huán)控制仿真模型

圖2 力閉環(huán)控制仿真模型最大軋制力下穩(wěn)態(tài)精度仿真曲線

2 FMU模型

FMI(Functional Mock-up Interface)是一種工具或者獨立的標準，通過利用xml文件和編譯好的C代碼程序，提供動態(tài)模型的數(shù)據(jù)交換和聯(lián)合仿真。 FMI能夠?qū)⒛承┨厥夥抡孳浖⒌哪Ｐ娃D(zhuǎn)換為FMU(Functional Mock-up Unit)模型。當前FMI& FMU技術(shù)標準被78種工具軟件所支持，比如Abaqus、Dymola 、Matlab和Adams等[3]。

在轉(zhuǎn)換為FMU模型之前，在Matlab/Simulink 中的控制模型需要完成一些準備工作，因它經(jīng)常不符合FMI特征。首先，模型配置參數(shù)中的求解器類型必須是固定步類型，因為聯(lián)合仿真需要在某個確定的采樣時間中在不同的子系統(tǒng)之間完成信號交換；然后，所有的I/O接口需要設(shè)置為單向數(shù)值接口，以取代雙向物理連接接口；最后，一個合適的基礎(chǔ)采樣時間間隔需要設(shè)置，它一般取決于系統(tǒng)固有頻率、驅(qū)動頻率等等。在本液壓控制模型，固定步的求解器為ode3，xi/i/Fi(i=1，2，3，4)抽象為單向數(shù)值I/O接口，基礎(chǔ)采樣時間為1 ms。

準備工作完成后，液壓控制模型就可以轉(zhuǎn)化為FMU模型。在配置參數(shù)中，F(xiàn)MI版本設(shè)置為2.0，類型設(shè)置為聯(lián)合仿真。此時，F(xiàn)MU模型(舉例命名為hydcontrol.fmu)建立起來，這就是參與到聯(lián)合仿真中的實際控制模型。

3 機械系統(tǒng)

平整機主要由牌坊、支撐輥系、工作輥系、換輥裝置及壓下裝置等機構(gòu)構(gòu)成，為減少計算規(guī)模，建模時根據(jù)系統(tǒng)的質(zhì)量及剛度特性，僅保留結(jié)構(gòu)主要構(gòu)件，并根據(jù)力的轉(zhuǎn)遞路徑，對平整機進行了簡化，去掉倒角，小孔等小特征；對墊片、套筒、油缸等進行了簡化處理，簡化后的幾何模型如圖3所示。在ABAQUS中建立有限元模型如圖4所示，模型考慮了輥系之間，以及輥系與軸承座之間的接觸關(guān)系，并對接觸區(qū)域的網(wǎng)格進行了細化。

圖3 平整機幾何模型

圖4 平整機有限元模型

平整機工作時，上部兩個AGC液壓缸在液壓系統(tǒng)作用下，向上支撐輥和牌坊頂部施加作用力，使得支撐輥與工作輥，工作輥與帶鋼互相壓靠在一起。

通過壓下模擬獲得了系統(tǒng)的綜合剛度并與現(xiàn)場靜壓靠測試數(shù)據(jù)相比較，誤差在7%以內(nèi)，確認了有限元模型的正確性，如表1所示。

表1 平整機剛度對比結(jié)果 kN/mm

在平整機有限元模型中，由帶激勵特性(actuator)的幅值信號定義輸入，由帶傳感器特性(sensor)的歷史輸出定義輸出，實現(xiàn)和FMU模型的對接。

4 平整機聯(lián)合仿真

聯(lián)合仿真模擬平整機在恒軋制力閉環(huán)控制模式和給定工況條件下，模擬軋制力調(diào)整過程和時間。

系統(tǒng)框圖如圖5所示，平整機壓靠過程，液壓控制系統(tǒng)在平整機的操作端和傳動端輸出作用力(HC_ROLLWF和HC_ROLLDF)，驅(qū)動平整機(即對應(yīng)的AMP_ROLLWF和AMP_ROLLDF，屬于機械有限元模型)工作輥完成壓靠，機械系統(tǒng)向液壓控制系統(tǒng)輸出操作端和傳動端的位移和速度(AGC_ROLLWU、AGC_ROLLWV、AGC_ROLLDU和AGC_ROLLDV)，其中位移為計算過程值，速度數(shù)據(jù)與液壓控制系統(tǒng)輸出作用力和目標作用力的差值可獲得新的輸出力，以完成閉環(huán)控制，整個聯(lián)合仿真共有六條數(shù)據(jù)交換傳輸通道。

圖5 聯(lián)合仿真系統(tǒng)框圖

現(xiàn)場基本采用恒壓力控制模式生產(chǎn)，故本液壓控制模型也采用恒壓力閉環(huán)控制。液壓油經(jīng)過恒壓泵后形成高壓油，通過電液伺服閥控制后，驅(qū)動單向液壓缸(AGC缸)作用于上支承輥，完成平整工作；液壓缸回退采用3 MPa的恒定背壓方式。以液壓缸輸出力作為負反饋信號，通過標準PID控制恒定壓力。液壓控制模型如圖6所示。

圖6 液壓控制模型

采用技術(shù)保證值的設(shè)定工況，AGC缸先緩慢壓下支承輥，壓力值為300 t，然后階躍提升到360 t，要求軋制力控制精度為±1.0%，階躍響應(yīng)時間不大于75 ms(0到95%的上升時間)，最大超調(diào)量為不超過10%目標值。

聯(lián)合仿真模型的計算步長是涉及計算效率和收斂特性的關(guān)鍵參數(shù)，一般取參與聯(lián)合仿真的各個子模型最大固定步長的最大公約數(shù)，本文經(jīng)過測算，取0.5 ms較為合適。軋制力實際控制曲線和目標曲線如圖7所示：

圖7 軋制力實際控制曲線和目標曲線

仿真結(jié)果表明性能指標滿足保證值要求，與現(xiàn)場驗收情況一致。同時仿真獲得的控制指標略優(yōu)于現(xiàn)場結(jié)果，這主要是由于仿真中各項條件相對理想，忽略了實際現(xiàn)場中的各種誤差所致。

保證值仿真結(jié)果軋制力控制精度±1%-0.66%～0.58%壓下響應(yīng)時間階躍響應(yīng)時間75 ms階躍響應(yīng)時間43.7 ms

5 結(jié) 論

本文通過基于FMI的聯(lián)合仿真方法實現(xiàn)了考慮平整機系統(tǒng)非線性特性的機液控聯(lián)合仿真，實現(xiàn)了對裝備系統(tǒng)綜合性能的預測，該方法能夠充分考慮復雜裝備工作時的非線性特性如接觸、大變形、熱力耦合等，為復雜裝備系統(tǒng)的性能預測、子系統(tǒng)優(yōu)化設(shè)計提供了新的技術(shù)手段；如進一步將工藝過程進行模擬，可為工藝質(zhì)量的預測以及裝備綜合性能指標的確定提供技術(shù)依據(jù)。