[發(fā)明專利]一種基于動態(tài)響應(yīng)的電磁減振器輸出力計算方法

在審

• 申請?zhí)枺?/td>: 202110052036.4
• 申請日: 2021-01-15
• 公開（公告）號: CN112733299A
• 公開（公告）日: 2021-04-30
• 發(fā)明（設(shè)計）人: 李競;王子凡;王佩犇;關(guān)棟;張純;高悅文;沈輝;龔俊杰;楊征睿;朱仲文
• 申請（專利權(quán)）人: 揚州大學(xué)
• 主分類號: G06F30/17
• 分類號: G06F30/17;G06F30/20;G06F17/11;G06F119/14
• 代理公司: 揚州蘇中專利事務(wù)所(普通合伙) 32222
• 代理人: 沈志海

1.一種基于動態(tài)響應(yīng)的電磁減振器輸出力計算方法，其特征在于，所述計算方法包括如下步驟：

(A)建立電磁減振器的動力學(xué)等效模型

將電磁減振器中所有運動件的運動慣量等效為質(zhì)量M，將電磁減振器中所有的嚙合剛度和結(jié)構(gòu)剛度等效為剛度K，將電磁減振器中所有的干摩擦等效為內(nèi)部摩擦F_r，將電磁減振器中所有的粘性阻尼等效為阻尼C，因此電磁減振器可以等效為一個具有位移激勵的單自由度系統(tǒng)，該單自由度系統(tǒng)中質(zhì)量為M，位移激勵y通過剛度K與質(zhì)量相連，質(zhì)量M的下端為剛度K，側(cè)端為干摩擦F_r，上端為粘性阻尼C，粘性阻尼C的另一端為固定邊界；

(B)確定動力學(xué)等效模型中的等效參數(shù)

(B1)計算電磁減振器中所有運動件的等效質(zhì)量M，即：

其中，m_r表示齒條質(zhì)量，J_p、J_s、J_c、J_g、J_m分別表示齒輪、齒輪軸、聯(lián)軸器、齒輪箱內(nèi)齒輪和電機轉(zhuǎn)子的轉(zhuǎn)動慣量，η_r表示齒條齒輪之間的傳動效率，r表示齒條齒輪副中齒輪的半徑，i_g表示齒輪箱內(nèi)齒輪的總傳動比；

(B2)計算等效摩擦力F_r

通過對電磁減振器進行加載試驗，在超低頻加載條件下考察其輸出力與加載位移的關(guān)系，在每個作動沖程末期，電磁減振器的輸出力約等于內(nèi)部摩擦力，即F_r；

(B3)計算等效剛度K與等效阻尼C

測量單個沖程初期具有的輸出力波動中的波峰和波谷的輸出力F₁、F₂、F₃和F₄，測量各波峰波谷對應(yīng)的時刻t₁、t₂、t₃和t₄，則有

以及有阻尼自由振動的周期時長T_d，即

同時，相鄰兩個波峰或波谷的衰減系數(shù)為

式中ζ為阻尼比，ω_n為結(jié)構(gòu)的固有頻率，T_d為有阻尼自由振動的周期；表達式如下：

式(7)中，ω_d為有阻尼下結(jié)構(gòu)的固有頻率，表達式如下：

聯(lián)立式(2)至式(8)，最終求得等效剛度K和等效阻尼C；

(C)對動力學(xué)等效模型進行求解計算

由于系統(tǒng)中存在摩擦力，因此采用分段分析的方法對其動力學(xué)特性進行分析，按照質(zhì)量塊的運動狀態(tài)可將系統(tǒng)劃分為四個階段，即質(zhì)量塊靜止、質(zhì)量塊向上運動、質(zhì)量塊靜止、質(zhì)量塊向下運動，分別對這四個階段中的動力學(xué)方程與輸出力進行分析于計算；

(C1)質(zhì)量塊靜止時的輸出力計算

記第一階段的開始時間為t_0,n，t_0,n的角標(biāo)0表示第一階段開始，角標(biāo)n表示第n次加載，則第一階段的初始條件，即質(zhì)量塊所處位置與速度分別為：

式中x_0,n表示第n次加載循環(huán)中第一階段開始時刻質(zhì)量塊所處位置，t_4,n-1表示上一個加載循環(huán)第四階段結(jié)束時刻，對于第一個加載周期而言，其初始條件為：

假設(shè)此時激勵位移向上運動，持續(xù)壓縮彈簧，在此階段中，由于彈簧產(chǎn)生的壓縮彈性力小于結(jié)構(gòu)內(nèi)部的靜摩擦力，因此質(zhì)量塊位移保持不變，速度為0，外部輸出力等于彈簧壓縮力，即：

該階段將持續(xù)到壓縮彈性力等于靜摩擦力為止，即：

F(t_1,n)＝K(y(t_1,n)-x_0,n)＝μF_r (12)

式中μ表示靜動摩擦力系數(shù)，一般取1.05；t_1,n即表示第n次加載循環(huán)中第一階段結(jié)束時刻；隨后，質(zhì)量塊將開始運動，系統(tǒng)將進入第二階段；此時，有第二階段的初始條件：

(C2)質(zhì)量塊向上運動時的輸出力計算

由于質(zhì)量塊開始運動，因此摩擦力從靜摩擦力變?yōu)閯幽Σ亮?，且摩擦力的方向與質(zhì)量塊運動方向相反，而質(zhì)量塊在彈簧的彈性力，粘性阻尼的阻尼力和摩擦力的共同作用下產(chǎn)生運動，因此有動力學(xué)方程：

該方程的解為：

x(t)＝x_P(t)+x_H(t) (15)

其中，x_P為方程的通解，x_H為方程的特解，特解x_H可以表示為：

其中

通解x_P的具體形式與阻尼大小有關(guān)，即：

式(18)中的A₁和A₂則與第二階段的初始條件有關(guān)，將初始條件(13)與式(18)和式(16)共同帶入方程(15)中即求解出A₁和A₂，再將求解出的A₁和A₂代入式(15)即可得到質(zhì)量塊的位移隨時間的變化規(guī)律，此時，系統(tǒng)對外的輸出力有：

F＝K(y(t)-x(t)) (19)

隨著向上位移激勵的加載速度慢慢降低并最終反向，同時加上摩擦阻尼與粘性阻尼的影響，質(zhì)量塊的速度也將下降至0，此時第二階段結(jié)束，設(shè)此時刻為t_2,n，即：

求解式(20)可得t_2,n，將t_2,n代入式(15)可得第二階段結(jié)束時質(zhì)量塊所處位置，由此可得第三階段的初始條件，即：

隨后系統(tǒng)進入第三階段；

(C3)質(zhì)量塊靜止時的輸出力計算

位移激勵應(yīng)開始向下運動，彈簧的狀態(tài)逐漸由壓縮狀態(tài)改變?yōu)槔鞝顟B(tài)，進而產(chǎn)生拉伸彈性力；當(dāng)拉伸彈性力小于靜摩擦力時，則質(zhì)量塊繼續(xù)保持靜止?fàn)顟B(tài)；此時質(zhì)量塊位移保持不變，速度為0，外部輸出力等于彈簧壓縮力，即：

該階段將持續(xù)到拉伸彈性力等于靜摩擦力為止，即：

F(t_3,n)＝K(y(t_3,n)-x_2,n)＝μF_r (23)

該t_3,n即表示第n次加載循環(huán)中第三階段結(jié)束時刻；此時，可得第四階段的初始條件，即：

隨后質(zhì)量塊將開始運動，系統(tǒng)將進入第四階段；

(C4)質(zhì)量塊向下運動時的輸出力計算

第四階段與第二階段類似，不同點在于質(zhì)量塊的運動方向與摩擦力的方向均與第二階段相反；因此，此時的動力學(xué)方程為：

該方程的解的形式如式(15)，即：

x(t)＝x_P(t)+x_H(t) (26)

其中，x_P為方程的通解，x_H為方程的特解，特解x_H可以表示為：

其通解x_P的形式與式(18)相同，通解中的系數(shù)A₁和A₂由第四階段的初始條件有關(guān)，將初始條件(24)與式(18)和式(27)代入方程(26)中即求解出A₁和A₂，再將求解出的A₁和A₂代入式(26)即可得到質(zhì)量塊的位移隨時間的變化規(guī)律，此時，系統(tǒng)對外的輸出力有：

F＝K(y(t)-x(t)) (28)

隨著向下位移激勵的加載速度慢慢降低并最終反向，同時加上摩擦阻尼與粘性阻尼的影響，質(zhì)量塊的速度也將下降至0，此時第四階段結(jié)束，設(shè)此時刻為t_4,n，即：

求解式(29)可得t_4,n，將t_4,n代入式(26)可得第四階段結(jié)束時質(zhì)量塊所處位置，即：

x(t_4,n)＝x_P(t_4,n)+x_H(t_4,n) (30)

隨后，系統(tǒng)進入下一個加載周期的第一階段，則下一個加載周期的第一階段的初始條件為：

至此，一個加載周期中四個階段的輸出力均計算完成，通過設(shè)定計算時長，可以進一步算出n個加載周期下輸出力的變化規(guī)律。