基于MATLAB/Simulink的單頭螺桿馬達(dá)運動仿真

于虹 崔勇 陳引生 陶德清 劉晶晶

摘 要：介紹了螺桿馬達(dá)的工作原理，研究分析了單頭螺桿馬達(dá)轉(zhuǎn)子的自轉(zhuǎn)及公轉(zhuǎn)規(guī)律。利用MATLAB仿真工具中的Simulink軟件包仿真了馬達(dá)的運動規(guī)律并導(dǎo)出運動參數(shù)，再將數(shù)據(jù)導(dǎo)入到ADAMS中建立AKISPL樣條函數(shù)，驅(qū)動馬達(dá)轉(zhuǎn)子運轉(zhuǎn)，模擬出了螺桿馬達(dá)在工況狀態(tài)下的運轉(zhuǎn)狀態(tài)。

關(guān)鍵詞：MATLAB/Simulink；單頭螺桿馬達(dá)；ADAMS；運動仿真

中圖分類號：TH137 文獻(xiàn)標(biāo)識碼：A

Single head Screw motor motion simulation based on MATLAB/Simulink

Yu Hong CuiYong ChenYinsheng TaoDeqing LiuJingjing

Jiangsu Agrianimal Husbandry Vocational College JiangsuTaizhou 225300

Abstract：`Introduces the working principle of screw motor，conducting the research and the revolution of the single screw motor rotor's autobiography. Application of MATLAB dynamic simulation software package Simulink tools for motor motion simulation and derived motion data， import it into ADAMS to establish AKISPL spline function ， drive motor rotor operation， simulated the screw motor running state in the working condition.

Key words：MATLAB/Simulink；single screw motor；ADAMS；Motion simulation

螺桿馬達(dá)主要由定子、轉(zhuǎn)子組成，它屬于容積式動力機。定子和轉(zhuǎn)子互相嚙合，液體通過擠壓兩者之間導(dǎo)程差形成的螺旋密封腔，推動轉(zhuǎn)子在定子內(nèi)部不斷旋轉(zhuǎn)。它也可以作為一個可以將壓力能轉(zhuǎn)化為動能能量的轉(zhuǎn)換機構(gòu)。

1 單頭螺桿馬達(dá)的運動規(guī)律分析

（1） 單頭螺桿馬達(dá)轉(zhuǎn)子的公轉(zhuǎn)。

從單頭螺桿馬達(dá)X=0的平面（圖1）來看，轉(zhuǎn)子的公轉(zhuǎn)是按照正弦的運動規(guī)律沿著Z軸上下往復(fù)移動，在兩端速度為0，在中間位置時速度達(dá)到最大。

整體來看螺桿馬達(dá)運動情況，轉(zhuǎn)子的公轉(zhuǎn)是轉(zhuǎn)子質(zhì)心繞定子軸心做勻速圓周運動的一個過程。

（2） 單頭螺桿馬達(dá)轉(zhuǎn)子的自轉(zhuǎn)。

螺桿馬達(dá)轉(zhuǎn)子的自轉(zhuǎn)是轉(zhuǎn)子以自身質(zhì)心為軸線進(jìn)行勻速轉(zhuǎn)動。為研究轉(zhuǎn)子的自轉(zhuǎn)和公轉(zhuǎn)之間的關(guān)系，將馬達(dá)轉(zhuǎn)子簡化為圖2所示模型。該模型顯示的是半個螺距轉(zhuǎn)子，轉(zhuǎn)子的兩端分別標(biāo)號1和2。這半個螺距的轉(zhuǎn)子在自轉(zhuǎn)的同時公轉(zhuǎn)，最后表現(xiàn)為端點1沿著縱軸往復(fù)運動，端點2沿著橫軸往復(fù)運動。

當(dāng)簡化轉(zhuǎn)子初始位置處在垂直處時，即端點1位于縱軸頂點，端點2位于中心點時，若質(zhì)心公轉(zhuǎn)運動速度為π12rad/s，則端點1的運動規(guī)律為Asinπ12t，第6s時端點1將移動至中心點，其位移為L。

L=∫60Asinπ12tdtA=π12L（1）

則：

S=π12L∫t0sinπ12tdt=L-Lcosπ12t（2）

由幾何關(guān)系可知：

S=L-Lcosαcosα=L-SLα=arccos（L-SL）

ω=α′=[arccos（L-SL）]'=π12rad/s（3）

由上述推導(dǎo)可知，馬達(dá)轉(zhuǎn)子自轉(zhuǎn)速度為π12rad/s，與公轉(zhuǎn)速度相同。

2 數(shù)學(xué)模型的建立

（1） 線性化流量方程。

線性化流量方程可表示為：

對式（4）進(jìn)行拉普拉斯變換得：

式中油液液阻Rv=1Zv。

線性化扭矩方程可表示為：

液壓馬達(dá)輸出扭矩M用來抑制負(fù)載扭矩Mz，即：

式中ΔPu為實際工作壓降。

扭矩方程可改寫為：

式中kp=k2πVt=（Vt2π+k′）2πVt=1+2πk′Vt，液阻Rp = 2πbn V2t ，液感L=（2πVt）2J。

對式（10）進(jìn)行拉普拉斯變換得：

由于扭矩?fù)p失系數(shù)k′很小，可略去，故將式（6）（11）整理為：

式中Th為液壓時間常數(shù)（Th=RvC），Tm為機械時間常數(shù)（Tm=LRp）。

將式（12）（13）繪制成如圖3所示的系統(tǒng)動態(tài)圖。

由圖3 可得，液壓馬達(dá)的傳遞函數(shù)Fqq（s）為：

式中放大系數(shù)Kqq=11+ZVRp；時間常數(shù)TM=ThTmKm；阻尼系數(shù)aM=（Th+Tm）Km2Tm；系數(shù)Km=ZVRp1+ZVRp。

3 系統(tǒng)仿真

（1） 仿真模型的建立。

根據(jù)圖3建立Simulink仿真模型，如圖4所示。

將螺桿馬達(dá)的仿真參數(shù)（如下表）輸入到Simulink仿真框圖中，運行Simulink后得到動態(tài)特性仿真曲線，如圖5所示。

由圖5可知，液壓馬達(dá)的速度波動發(fā)生在啟動初期0～0.01s過程中，在0.002s時達(dá)到最大速度8.6r/s。0.01s之后，馬達(dá)將達(dá)到穩(wěn)定轉(zhuǎn)速6.8 r/s。

（2） 仿真數(shù)據(jù)的導(dǎo)出。

應(yīng)用圖4中的To Workplace模塊導(dǎo)出仿真數(shù)據(jù)，并保存至文本文檔，同時檢查數(shù)據(jù)的完整性及連續(xù)性。

4 虛擬樣機仿真

（1） 虛擬樣機的建立。

運用三維建模軟件CATIA對螺桿馬達(dá)零件進(jìn)行實體建模，并在裝配模塊中將所建零件體進(jìn)行裝配。將裝配好的模型文件轉(zhuǎn)換為ADAMS可識別的Parasolid格式。

（2） 運動仿真。

將模型導(dǎo)入到ADAMS，按照螺桿馬達(dá)的實際運動狀態(tài)對各部分進(jìn)行約束，并定義相應(yīng)的運動副，保證轉(zhuǎn)子在定子內(nèi)按照設(shè)定的軌跡運動。

將文本文檔中的數(shù)據(jù)文件導(dǎo)入到ADAMS中建立樣條線SPLINE_1。將該樣條線編輯成轉(zhuǎn)動驅(qū)動函數(shù)：AKISPL（time，0，SPLINE_1， 0）加載至相應(yīng)的運動副上，最終建立如圖6所示的樣機。然后在ADAMS求解器中求解，模擬出螺桿馬達(dá)從啟動到穩(wěn)工作狀態(tài)的動態(tài)變化圖。

運用ADAMS后處理模塊（PostPocess）繪制出虛擬樣機仿真結(jié)果如圖7所示，結(jié)果圖中顯示了轉(zhuǎn)子角速度和角加速度變化情況。藍(lán)色虛線顯示了角加速度的變化。在啟動初期0～001s過程中角加速度波動較大，在0.001s時，角加速度達(dá)到峰值7000r/s2左右。0.01s之后，馬達(dá)進(jìn)入穩(wěn)定工作狀態(tài)。

5 結(jié)語

分析了單頭螺桿馬達(dá)轉(zhuǎn)子的運動規(guī)律，證明了轉(zhuǎn)子公轉(zhuǎn)和自轉(zhuǎn)具有等速的特性，為運動仿真過程中各運動副驅(qū)動函數(shù)的設(shè)置提供了依據(jù)。

應(yīng)用MATLAB中的Simulink模塊對馬達(dá)工況條件下的運動規(guī)律進(jìn)行了仿真，將仿真的結(jié)果作為驅(qū)動函數(shù)驅(qū)動在ADAMS中建立的虛擬樣機運動。通過ADAMS虛擬樣機仿真更直觀的觀察到單頭螺桿馬達(dá)的運動狀況，得出的仿真結(jié)果圖也能更準(zhǔn)確地展示螺桿馬達(dá)從啟動到穩(wěn)定狀態(tài)的角速度和角加速度的變化。

參考文獻(xiàn)：

[1]蔣文斌，劉壽康，鄒柏華.基于MATLAB/Simulink的液壓馬達(dá)低速穩(wěn)定性仿真研究[J].礦冶工程，2008.1：9496.

[2]張德虎，朱建公，肖躍軍.基于MATLAB的液壓系統(tǒng)數(shù)字仿真研究[J].機床與液壓，2006，8：202203.

[3]劉春芳，吳盛林，賈錦虹.液壓仿真轉(zhuǎn)臺中液壓爬行現(xiàn)象分析及消除措施[J].液壓與氣動， 2002，5：1819.

[4]蔡廷文.液壓系統(tǒng)現(xiàn)代建模方法[M].中國標(biāo)準(zhǔn)出版社，2002，12.