本發(fā)明涉及雙激振電機驅(qū)動的雙質(zhì)體非線性振動系統(tǒng)技術(shù)領(lǐng)域,更具體地,涉及一種基于頻率俘獲的雙激振電機同步特性仿真分析方法。
背景技術(shù):
實際工程中大部分的機械振動都是非理想的非線性振動系統(tǒng),非線性振動通常是指恢復力與位移不成線性比例或阻尼力與速度不成線性比例的系統(tǒng)的振動。頻率俘獲是指非線性振動系統(tǒng)的激振頻率接近系統(tǒng)固有頻率時,激振頻率被固有頻率所俘獲的現(xiàn)象,即激振頻率與固有頻率進入同步狀態(tài),頻率俘獲情況下的系統(tǒng)能夠利用小的激振力獲得較大的振幅,提高其工作效率。
多激振電機驅(qū)動的振動壓實系統(tǒng)在壓實土壤過程中,土體對振動體的恢復力與其位移之間形成的滯回環(huán)是不對稱的,這種在應力應變關(guān)系曲線上的不對稱滯回特性對系統(tǒng)的振動響應有很大影響,須保證多激振電機的同步運轉(zhuǎn)才能實現(xiàn)系統(tǒng)同步運動穩(wěn)定性,進而保證振動機械的工作性能。
現(xiàn)有的振動實壓類工程機械,為保證多激振電機能實現(xiàn)同步運轉(zhuǎn)穩(wěn)定和系統(tǒng)振幅穩(wěn)定性,傳統(tǒng)的振動機械大部分都工作在遠超共振狀態(tài),現(xiàn)有振動系統(tǒng)的同步設計理論多是針對遠超共振單質(zhì)體振動系統(tǒng)的振動同步理論,而多激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的同步運轉(zhuǎn)特性分析缺少準確實用的分析理論。
目前針對多激振電機實現(xiàn)同步運轉(zhuǎn)和系統(tǒng)振幅穩(wěn)定性的分析理論推導過程繁瑣,人員工作量大,準確度不夠理想,而且不能及時有效地解決自同步系統(tǒng)在實際工作過程中的問題,使得工程實際中試驗的成本較高。
技術(shù)實現(xiàn)要素:
本發(fā)明提供一種克服上述問題或者至少部分地解決上述問題的基于頻率俘獲的雙激振電機同步特性仿真分析方法,該方法可解決現(xiàn)有多激振電機的同步運轉(zhuǎn)和系統(tǒng)振幅穩(wěn)定性分析問題。
根據(jù)本發(fā)明的技術(shù)方案,提供一種基于頻率俘獲的雙激振電機同步特性仿真分析方法,包括:
s1.根據(jù)雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng),建立雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的計算機仿真模型;
s2.基于計算機仿真模型,設置系統(tǒng)參數(shù)和初始輸入信號,對雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)進行模型仿真;
s3.根據(jù)仿真結(jié)果,對雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)在頻率俘獲情況下的雙激振電機同步運轉(zhuǎn)和系統(tǒng)穩(wěn)定性進行分析。
作為上述技術(shù)方案的進一步改進方案,所述步驟s1中,包括:
s11.根據(jù)雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng),建立雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的數(shù)學模型;
s12.根據(jù)雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的數(shù)學模型,建立雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的計算機仿真模型。
作為上述技術(shù)方案的進一步改進方案,所述步驟s11中,根據(jù)雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng),建立雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)模型,再根據(jù)系統(tǒng)模型建立雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的數(shù)學模型。
作為上述技術(shù)方案的進一步改進方案,所述雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)模型包括第一質(zhì)體和第二質(zhì)體,第二質(zhì)體包括振動體和設置在振動體內(nèi)的兩個偏心轉(zhuǎn)子,兩個偏心轉(zhuǎn)子分別繞各自的回轉(zhuǎn)軸心旋轉(zhuǎn);兩個偏心轉(zhuǎn)子的質(zhì)量、旋轉(zhuǎn)半徑、回轉(zhuǎn)阻尼及轉(zhuǎn)動慣量均相同。
具體地,所述雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)數(shù)學模型為:
式中,m1為第一質(zhì)體的質(zhì)量,m2為第二質(zhì)體的質(zhì)量,第二質(zhì)體的質(zhì)量m2包括振動體的質(zhì)量和兩個偏心轉(zhuǎn)子的偏心質(zhì)量,兩個偏心轉(zhuǎn)子的偏心質(zhì)量分別為m01和m02,兩個偏心轉(zhuǎn)子繞著各自的回轉(zhuǎn)軸心旋轉(zhuǎn)的半徑分別為r1和r2;第一質(zhì)體的豎直方向的位移為x1,速度為
作為上述技術(shù)方案的進一步改進方案,通過下述公式確定所述不對稱滯回力f(x2),
式中,x為非線性振動系統(tǒng)的振動位移,f(x)為不對稱滯回力,k′為滯回剛度,
作為上述技術(shù)方案的進一步改進方案,所述步驟s2中,系統(tǒng)參數(shù)包括質(zhì)體質(zhì)量、振動位移信號、轉(zhuǎn)動慣量、雙激振電機的轉(zhuǎn)速及轉(zhuǎn)速差,初始輸入信號包括初始轉(zhuǎn)速和初始相位。
作為上述技術(shù)方案的進一步改進方案,所述步驟s3中,根據(jù)仿真結(jié)果得到雙激振電機的轉(zhuǎn)速響應,判斷系統(tǒng)工作頻率是否被系統(tǒng)的固有頻率所俘獲,即系統(tǒng)發(fā)生頻率俘獲。
作為上述技術(shù)方案的進一步改進方案,計算機仿真模型采用變步長類型的ode45的求解方法,允相對誤差采用0.001。
作為上述技術(shù)方案的進一步改進方案,通過計算機仿真模型求解方法,得到系統(tǒng)的雙激振電機轉(zhuǎn)速響應圖以及初始轉(zhuǎn)速不同時雙激振電機的相位差和轉(zhuǎn)速差的相平面圖。
基于上述技術(shù)方案,本發(fā)明提出的基于頻率俘獲的雙激振電機同步特性仿真分析方法,以雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)模型為仿真研究對象,針對不對稱滯回雙質(zhì)體振動模型在頻率俘獲情況下的雙激振電機同步運轉(zhuǎn)特性進行仿真分析,利用計算機軟件可方便呈現(xiàn)多激振電機能實現(xiàn)同步運轉(zhuǎn)和系統(tǒng)振幅穩(wěn)定性,通過系統(tǒng)頻率俘獲情況下的同步特性仿真,使實現(xiàn)多激振電機能實現(xiàn)同步運轉(zhuǎn)穩(wěn)定和系統(tǒng)振幅穩(wěn)定性更加方便和簡單,也免去的傳統(tǒng)理論推導。同時還能夠呈現(xiàn)和解決自同步系統(tǒng)在實際工作過程中的問題,降低工程實際試驗的成本等,此外還可促進現(xiàn)有自同步振動系統(tǒng)的同步理論進一步豐富和發(fā)展。
附圖說明
圖1為根據(jù)本發(fā)明實施例的基于頻率俘獲的雙激振電機同步特性仿真分析方法流程圖;
圖2為根據(jù)本發(fā)明實施例的雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)模型的示意圖;
圖3為根據(jù)本發(fā)明實施例的不對稱滯回模型的示意圖;
圖4為根據(jù)本發(fā)明實施例的雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的計算機仿真模型的示意圖;
圖5為根據(jù)本發(fā)明實施例的雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的仿真模型中的子系統(tǒng)subsystem的示意圖;
圖6為根據(jù)本發(fā)明實施例的雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的仿真模型中的子系統(tǒng)subsystem4的示意圖;
圖7為根據(jù)本發(fā)明實施例的頻率俘獲情況下雙激振電機的轉(zhuǎn)速響應的示意圖;
圖8為根據(jù)本發(fā)明實施例的雙激振電機初始條件不同時相位差和轉(zhuǎn)速差的相平面圖。
具體實施方式
為使本發(fā)明實施例的目的、技術(shù)方案和優(yōu)點更加清楚,下面將結(jié)合本發(fā)明實施例中的附圖,對本發(fā)明實施例中的技術(shù)方案進行清楚、完整地描述,顯然,所描述的實施例是本發(fā)明一部分實施例,而不是全部的實施例?;诒景l(fā)明中的實施例及其之間任意組合,本領(lǐng)域普通技術(shù)人員在沒有作出創(chuàng)造性勞動前提下所獲得的所有其他實施例,都屬于本發(fā)明保護的范圍。
在根據(jù)本申請的一個實施例中,參考圖1,提供一種基于頻率俘獲的雙激振電機同步特性仿真分析方法,包括:
s1.根據(jù)雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng),建立雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的計算機仿真模型;
s2.基于計算機仿真模型,設置系統(tǒng)參數(shù)和初始輸入信號,對雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)進行模型仿真;
s3.根據(jù)仿真結(jié)果,對雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)在頻率俘獲情況下的雙激振電機同步運轉(zhuǎn)和系統(tǒng)穩(wěn)定性進行分析。
在根據(jù)本申請的一個實施例中,所述步驟s1中,進一步包括:
s11.根據(jù)雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng),建立雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的數(shù)學模型;
s12.根據(jù)雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的數(shù)學模型,建立雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的計算機仿真模型。
在根據(jù)本申請的一個實施例中,所述步驟s11中,根據(jù)雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng),建立雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)模型,如圖2所示,再根據(jù)系統(tǒng)模型建立雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的數(shù)學模型。
在根據(jù)本申請的一個實施例中,所述雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)模型包括第一質(zhì)體和第二質(zhì)體,第二質(zhì)體包括振動體和設置在振動體內(nèi)的兩個偏心轉(zhuǎn)子,兩個偏心轉(zhuǎn)子分別繞各自的回轉(zhuǎn)軸心旋轉(zhuǎn);兩個偏心轉(zhuǎn)子的質(zhì)量、旋轉(zhuǎn)半徑、回轉(zhuǎn)阻尼及轉(zhuǎn)動慣量均相同。
在根據(jù)本申請的一個實施例中,所述雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)數(shù)學模型的方程為:
式中,m1為第一質(zhì)體的質(zhì)量,m2為第二質(zhì)體的質(zhì)量,第二質(zhì)體的質(zhì)量m2包括振動體的質(zhì)量和兩個偏心轉(zhuǎn)子的偏心質(zhì)量,兩個偏心轉(zhuǎn)子的偏心質(zhì)量分別為m01和m02,兩個偏心轉(zhuǎn)子繞著各自的回轉(zhuǎn)軸心旋轉(zhuǎn)的半徑分別為r1和r2;第一質(zhì)體的豎直方向的位移為x1,速度為
雙質(zhì)體振動系統(tǒng)的邊界條件為第二質(zhì)體與固定地面之間存在非線性振動力fs,fs包括c2產(chǎn)生的阻尼力、k2產(chǎn)生的彈性力和不對稱滯回力f(x2)。系統(tǒng)的約束條件是第一質(zhì)體僅存在豎直方向的位移x1,第二質(zhì)體僅存在豎直方向的位移x2,而第一質(zhì)體和第二質(zhì)體的其他方向的位移都為0。
其中,不對稱滯回力f(x2),可通過下述的不對稱滯回模型確定:
式中,x為非線性振動系統(tǒng)的振動位移,f(x)為不對稱滯回力,k′為滯回剛度,
圖3為本發(fā)明實施例的不對稱滯回模型,其中a點的坐標為(xa,fa),b點的坐標為(xb,fb),c點的坐標為(xc,fc),d點的坐標為(0,0),且fa=fb,xa=xb-xc,d-a段與b-c段平行,a-b段與c-d段平行。
圖4為雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的計算機仿真模型。本發(fā)明實施例采用matlab/simlink計算機軟件,將系統(tǒng)的數(shù)學模型公式(1)和(2)轉(zhuǎn)換為計算機能處理的仿真模型。其中,輸出信號out1為第一質(zhì)體的x1的振動位移信號,out2為第二質(zhì)體的x2的振動位移信號,out3為振動系統(tǒng)上的雙激振電機的轉(zhuǎn)速差,out4為雙激振電機的相位差,out5為一激振電機的轉(zhuǎn)速,out6為另一激振電機的轉(zhuǎn)速。
圖4中,計算機仿真模型的函數(shù)fcn1的仿真函數(shù)為k1*(u(3)-u(1))+c1*(u(4)-u(2))。其中*表示乘號,k1為振動系統(tǒng)模型中的剛度k1,c1為振動系統(tǒng)模型中的阻尼c1,u(3)為振動系統(tǒng)模型中振動位移x2,u(1)為振動系統(tǒng)模型中振動位移x1,u(4)為系統(tǒng)模型中的
函數(shù)fcn4的仿真函數(shù)為-k1*(u(3)-u(1))-c1*(u(4)-u(2))-c2*u(4)-k2*u(3)-u(5)+m01*r1*(-u(11)*u(7)+u(10)*u(10)*u(6))-m02*r2*(u(13)*u(9)-u(12)*u(12)*u(8))。其中k1、k2、c1、c2、m01、m02、r1、r2分別表示振動系統(tǒng)模型中的k1、k2、c1、c2、m01、m02、r1、r2,u(1)代表系統(tǒng)模型中振動位移x1,u(2)代表系統(tǒng)模型中的
函數(shù)fcn2的仿真函數(shù)為m02*r2*u(2)*u(3)+m02*r2*g*u(3)+c02*u(1),其中,m02、r2、g、c02分別代表m02、r2、g、c20,u(1)代表振動系統(tǒng)模型中的
函數(shù)fcn3的仿真函數(shù)為m01*r1*u(2)*u(3)+m01*r1*g*u(3)+c01*u(1)。其中,m01、r1、c01分別代表m01、r1、c10,u(1)代表振動系統(tǒng)模型中的
仿真模型中的增益gain1為轉(zhuǎn)動慣量j20,gain4為轉(zhuǎn)動慣量j10,gain2為1/m1,gain3為1/m2。
圖5為雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的仿真模型中的子系統(tǒng)subsystem,對應式(2)所示的不對稱滯回力的仿真模型。
subsystem1和subsystem2為三相異步電機的仿真模型。
圖6為雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的仿真模型中的子系統(tǒng)subsystem4。
在根據(jù)本申請的一個實施例中,基于計算機仿真模型,設置系統(tǒng)參數(shù)和初始輸入信號,對雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)進行模型仿真,其中系統(tǒng)參數(shù)包括質(zhì)體質(zhì)量、振動位移信號、轉(zhuǎn)動慣量、雙激振電機的轉(zhuǎn)速及轉(zhuǎn)速差,初始輸入信號包括初始轉(zhuǎn)速和初始相位。
具體地,設定系統(tǒng)模型的相關(guān)參數(shù):m1=89kg,m2=56kg,k1=6500000n/m,k2=4000000n/m,c1=100nm·s/rad,c2=100nm·s/rad,m01=3.5kg,m02=3.5kg,r1=r2=0.08m,c01=0.01nm·s/rad,c02=0.01nm·s/rad,j01=0.01kg·m2,j02=0.01kg·m2,g=9.8m/s,另外滯回力參數(shù)選擇為:k′=400n/m,xa=0.01m,xb=0.02m,xc=0.01m。根據(jù)參數(shù)可知系統(tǒng)的固有頻率約為23hz。選擇激振電機為異步電機,其頻率為25hz(即157rad/s左右),也就是系統(tǒng)的固有頻率與激振電機的激振頻率比較接近。
根據(jù)仿真結(jié)果得到雙激振電機的轉(zhuǎn)速響應,判斷系統(tǒng)工作頻率是否被系統(tǒng)的固有頻率所俘獲,即系統(tǒng)發(fā)生頻率俘獲。
采用上述參數(shù),利用計算機工程軟件matlab/simlink對仿真模型進行仿真計算,計算機仿真模型采用變步長類型的ode45的求解方法,允相對誤差采用0.001。
通過計算機仿真模型求解方法,得到系統(tǒng)的雙激振電機轉(zhuǎn)速響應圖以及初始轉(zhuǎn)速不同時雙激振電機的相位差和轉(zhuǎn)速差的相平面圖。
在matlab/simlink計算機軟件中的仿真模型對系統(tǒng)進行各類初始條件不同情況下進行仿真模擬。首先,根據(jù)上述參數(shù),以及在雙激振電機的初始狀態(tài)相同的情況下,采用如圖4至圖6的仿真模型在計算機工程軟件matlab/simlink進行仿真,得到在正常工作狀態(tài)下系統(tǒng)的雙激振電機轉(zhuǎn)速響應,如圖7所示。由圖7可看出,雙激振電機的轉(zhuǎn)速最終穩(wěn)定在約142.6rad/s(22.7hz),與系統(tǒng)的固有頻率(23hz)接近,且可知系統(tǒng)的工作頻率已不是25hz(即157rad/s左右),意味著系統(tǒng)的工作頻率被系統(tǒng)的固有頻率所俘獲,即系統(tǒng)已發(fā)生頻率俘獲現(xiàn)象。
其次,取雙激振電機的初始轉(zhuǎn)速不同,即一激振電機的初始轉(zhuǎn)速滯后另一激振電機的初始轉(zhuǎn)速10rad/s,得到如圖8所示的雙激振電機的相位差和轉(zhuǎn)速差的相平面圖。由圖8可看出,當一激振電機的初始轉(zhuǎn)速滯后另一激振電機的初始轉(zhuǎn)速10rad/s時,雙激振電機的初始相位差為0rad和初始轉(zhuǎn)速差為-10rad/s,在經(jīng)歷了不規(guī)則的螺旋運動后,相位差和轉(zhuǎn)速差的相平面圖出現(xiàn)極限環(huán),即雙激振電機的轉(zhuǎn)速達到同步穩(wěn)定運轉(zhuǎn)和實現(xiàn)系統(tǒng)同步穩(wěn)定性。最終轉(zhuǎn)速差穩(wěn)定在0rad/s,相位差穩(wěn)定在4πrad(即0rad的2倍周期)左右。
本發(fā)明實施例在matlab/simlink計算機軟件中的仿真模型對系統(tǒng)進行各類初始條件不同情況下進行仿真實驗,由上述仿真方法分析得出,雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)選擇適當?shù)木€性阻尼時,系統(tǒng)的工作頻率(激振電機的頻率)能被系統(tǒng)的固有頻率所俘獲;系統(tǒng)在頻率俘獲的狀態(tài)下,即使雙激振電機存在一定范圍的初始差異,雙激振電機的轉(zhuǎn)速差和相位差都經(jīng)歷不規(guī)則運動后,最終轉(zhuǎn)速差趨于穩(wěn)定達到0rad/s,相位差最終穩(wěn)定在0rad(或2nπrad,其中n為整數(shù))位置左右。也就是說,雙激振電機存在差異時,雙質(zhì)體振動系統(tǒng)在頻率俘獲情況下仍能獲得雙激振電機的同步運轉(zhuǎn)和系統(tǒng)同步穩(wěn)定性。
本發(fā)明實施例提供雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)頻率俘獲情況下的同步特性仿真過程及分析方法,通過建立雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)模型,建立雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)數(shù)學模型,在將系統(tǒng)的數(shù)學模型轉(zhuǎn)換為計算機能處理的仿真模型,在
matlab/simlink計算機仿真模型的基礎上,給定系統(tǒng)的相關(guān)參數(shù)和初始條件的輸入信號,進行模型仿真和雙激振電機各類初始條件不同情況下進行仿真模擬。仿真分析結(jié)果能顯示出:雙激振電機的轉(zhuǎn)速被系統(tǒng)的固有頻率所俘獲,以及雙激振電機初始條件不同時,雙質(zhì)體振動系統(tǒng)在頻率俘獲情況下仍能獲得雙激振電機的同步運轉(zhuǎn)和系統(tǒng)同步穩(wěn)定性。
本發(fā)明實施例的仿真分析過程針對雙激振電機驅(qū)動的雙質(zhì)體振動系統(tǒng)的雙激振電機的同步運轉(zhuǎn)和系統(tǒng)的同步穩(wěn)定性問題,再現(xiàn)了該類自同步振動系統(tǒng)在頻率俘獲情況下的雙激振電機的同步運轉(zhuǎn)和系統(tǒng)的同步穩(wěn)定性性,為振動系統(tǒng)的同步設計提供了理論依據(jù)。
最后應說明的是:以上實施例僅用以說明本發(fā)明的技術(shù)方案,而非對其限制;盡管參照前述實施例對本發(fā)明進行了詳細的說明,本領(lǐng)域的普通技術(shù)人員應當理解:其依然可以對前述各實施例所記載的技術(shù)方案進行修改,或者對其中部分技術(shù)特征進行等同替換;而這些修改或者替換,并不使相應技術(shù)方案的本質(zhì)脫離本發(fā)明各實施例技術(shù)方案的精神和范圍。