楊志安，趙利沙

（1.唐山學(xué)院 唐山市結(jié)構(gòu)與振動(dòng)工程重點(diǎn)實(shí)驗(yàn)室，河北 唐山 063000；2.華北理工大學(xué) 機(jī)械工程學(xué)院，河北 唐山 063000）

有界窄帶激勵(lì)下的電磁開關(guān)系統(tǒng)主共振分析

楊志安1，趙利沙2

（1.唐山學(xué)院 唐山市結(jié)構(gòu)與振動(dòng)工程重點(diǎn)實(shí)驗(yàn)室，河北 唐山 063000；2.華北理工大學(xué) 機(jī)械工程學(xué)院，河北 唐山 063000）

以電磁開關(guān)系統(tǒng)為研究對象，研究電路與磁路耦合系統(tǒng)在有界窄帶激勵(lì)下的強(qiáng)非線性振動(dòng)問題。建立電磁開關(guān)強(qiáng)非線性系統(tǒng)在有界窄帶激勵(lì)下的隨機(jī)微分方程。應(yīng)用改進(jìn)的多尺度法得到系統(tǒng)的幅頻響應(yīng)方程，對窄帶激勵(lì)下系統(tǒng)主共振的穩(wěn)定性進(jìn)行計(jì)算，并分析系統(tǒng)各參數(shù)對系統(tǒng)主共振均方值的影響。結(jié)果表明，增大阻尼系數(shù)可以減小系統(tǒng)主共振的均方值；增大激勵(lì)電壓可以增大系統(tǒng)主共振的均方值；增大電阻可以減小系統(tǒng)主共振的均方值；增大隨機(jī)擾動(dòng)強(qiáng)度，極限環(huán)變化不大。

振動(dòng)與波；電磁開關(guān)；改進(jìn)多尺度法；窄帶激勵(lì)；均方值；強(qiáng)非線性振動(dòng)

近年來，隨著汽車市場的迅速發(fā)展，汽車已經(jīng)成為日常生活中不可分割的重要組成部分。汽車起動(dòng)機(jī)是汽車發(fā)動(dòng)機(jī)的關(guān)鍵零部件，現(xiàn)代汽車中幾乎每輛汽車都在使用，而電磁開關(guān)又是汽車起動(dòng)機(jī)中的一種重要的控制元件，其工作原理是線圈通電后產(chǎn)生電磁吸力，使活動(dòng)鐵芯移動(dòng)，接通勵(lì)磁繞組電路，產(chǎn)生勵(lì)磁，電樞工作，從而帶動(dòng)發(fā)動(dòng)機(jī)起動(dòng)。斷電后，活動(dòng)鐵芯在彈簧力下復(fù)位，等待下次啟動(dòng)[1]。電磁開關(guān)的設(shè)計(jì)質(zhì)量直接影響起動(dòng)機(jī)的可靠性，因此對電磁開關(guān)的研究具有十分重要的意義。

文獻(xiàn)[2]闡述了起動(dòng)機(jī)用電磁開關(guān)的工作原理，依據(jù)原理探討引起電磁開關(guān)觸片、觸點(diǎn)燒蝕的主要原因，并給出改善措施。文獻(xiàn)[3]闡述了電磁開關(guān)的組成及傳統(tǒng)結(jié)構(gòu)形式的利弊分析，并進(jìn)行改進(jìn)。文獻(xiàn)[4]研究非線性電阻電感型RLC串聯(lián)電路的非線性振動(dòng)，應(yīng)用多尺度法，得到系統(tǒng)滿足主共振條件的一次近似解以及對應(yīng)的定常解。文獻(xiàn)[5]研究環(huán)形極板機(jī)電耦合系統(tǒng)的強(qiáng)非線性問題。應(yīng)用多尺度法求得系統(tǒng)的主共振的幅頻響應(yīng)曲線，分析了不同系統(tǒng)參數(shù)對共振的影響。文獻(xiàn)[6]研究了Vander Pol-Duffing振子在窄帶隨機(jī)噪聲激勵(lì)下的響應(yīng)問題。用參數(shù)變換法使方程出現(xiàn)小參數(shù)，用多尺度法分離了系統(tǒng)的快變項(xiàng)，討論了系統(tǒng)的阻尼項(xiàng)、非線性項(xiàng)和隨機(jī)項(xiàng)等參數(shù)對系統(tǒng)的影響。文獻(xiàn)[7]從共振特征入手，對轉(zhuǎn)子的彎扭耦合振動(dòng)開展研究。

本文應(yīng)用改進(jìn)多尺度法研究電磁開關(guān)強(qiáng)非線性系統(tǒng)主共振在窄帶隨機(jī)激勵(lì)下的響應(yīng)[8]，得到主共振的幅頻響應(yīng)曲線和系統(tǒng)隨機(jī)均方響應(yīng)的近似表達(dá)式，考慮了隨機(jī)項(xiàng)對電磁開關(guān)系統(tǒng)響應(yīng)的影響。應(yīng)用4階龍格-庫塔法計(jì)算出系統(tǒng)的時(shí)間響應(yīng)曲線和相圖，對電磁開關(guān)系統(tǒng)的設(shè)計(jì)與優(yōu)化具有重要的參考意義。

1 有界窄帶激勵(lì)電磁開關(guān)系統(tǒng)強(qiáng)非線性動(dòng)力學(xué)方程

窄帶隨機(jī)過程是指頻譜密度集中在中心頻率附近相對較窄的頻率分布過程，它是工程中經(jīng)常遇到的隨機(jī)過程[9–10]，特點(diǎn)是頻譜主要集中在某段區(qū)間上。

電磁開關(guān)的數(shù)學(xué)模型如圖1所示，電源的電動(dòng)勢為U(t)，電阻為r，電路電流為i，線圈的自感為L(x)，電源、電阻及電感組成回路。電磁鐵的質(zhì)量為m，彈簧剛度為k，電磁鐵線圈的自感為Lm(x)，x為自彈簧平衡位置算起的水平位移，電磁鐵上線圈激勵(lì)電壓為Um，電磁鐵上線圈電阻為R，電磁鐵上線圈電流為im，M(x)為電路與電磁鐵的互感，rv為線圈與磁鐵之間的摩擦阻力系數(shù)。

圖1 電磁開關(guān)模型

選位移x，電路電流為i，電磁鐵上線圈電流為im為廣義坐標(biāo)。

將拉格朗日方程和損耗函數(shù)代入拉格朗日麥克斯韋方程[11]，建立系統(tǒng)的非線性動(dòng)力學(xué)方程為

在式（2）中，隨機(jī)電壓U(t)控制電磁開關(guān)在靜平衡位置附近的振動(dòng)位移，進(jìn)而分析電磁開關(guān)在有界窄帶激勵(lì)下的力學(xué)響應(yīng)。U(t)為有界隨機(jī)噪聲，它是由電路中信號(hào)發(fā)動(dòng)器發(fā)出信號(hào)，可表示為

式中H＞0為隨機(jī)激勵(lì)幅值，Ω為隨機(jī)激勵(lì)的中心頻率，W(t)為標(biāo)準(zhǔn)的Winner過程，γ為隨機(jī)的擾動(dòng)強(qiáng)度，γW(t)代表了電壓的隨機(jī)擾動(dòng)項(xiàng)，即U(t)可以看作為確定性諧和激勵(lì)的頻率Ω受到了擾動(dòng)強(qiáng)度為γW(t)的隨機(jī)擾動(dòng)。

文中研究窄帶隨機(jī)過程中的Winner隨機(jī)過程，Winner過程在應(yīng)用中為標(biāo)準(zhǔn)信號(hào)，包含了γ=0單頻激勵(lì)。應(yīng)用改進(jìn)多尺度法研究γ≥0遠(yuǎn)小于1時(shí)的情形。多尺度法分析隨機(jī)窄帶噪聲激勵(lì)下的響應(yīng)問題，已有一些成果[12–13]。所以式（3）變?yōu)?/p>

2 窄帶激勵(lì)主共振分析

引入非小量參數(shù)ε

將式（8）、式（9）代入式（5）后，同時(shí)對比α同次冪的系數(shù)，并對α0方程設(shè)解帶入到α1的方程中，根據(jù)消除長期項(xiàng)的條件，分離實(shí)虛部可得

假設(shè)γ=0，則式（10）變?yōu)?/p>

對于穩(wěn)態(tài)振動(dòng)a=a0，ν=ν0，有D1a=D1ν=0。得系統(tǒng)（11）的幅頻響應(yīng)方程

當(dāng)γ≠0時(shí)，由于γ很小，可設(shè)方程（10）的解為

對于穩(wěn)態(tài)振動(dòng)a=a0，ν=ν0，有D1a=D1ν=0。代入式（11）可以得到a0和φ0，其中a1和φ1是攝動(dòng)項(xiàng)，把式（13）代入式（10）中，不計(jì)關(guān)于a1和φ1的高階項(xiàng)，便得方程（11）在點(diǎn)（a0，φ0）的線性方程為

依據(jù)Hurwitz[14]定理可知，該線性系統(tǒng)穩(wěn)定的必要條件為

方程（14）可寫成對應(yīng)的Ito微分方程

可用矩法[14]求出式（19）的1階矩為

式（16）可寫成如下形式

以下各式中Q1、Q2和Q3同上。

由Ito微分法則可得

則方程（19）可轉(zhuǎn)化為

根據(jù)式（21）的系數(shù)矩陣，則

式（22）的特征方程為

式（25）即為2階矩存在的充分必要條件，并與前面確定性系統(tǒng)周期解得穩(wěn)定性條件是一致的。a0的值與Ea2的值是相互對應(yīng)的。

結(jié)合式（13）、式（17）、式（26）可得

3 數(shù)值模擬

圖2為電磁開關(guān)系統(tǒng)的幅頻響應(yīng)曲線，方程（5）中βx3項(xiàng)與項(xiàng)系數(shù)進(jìn)行比較β相對較大，屬于強(qiáng)非線性問題。由圖2可知，當(dāng)系統(tǒng)產(chǎn)生共振時(shí)振幅達(dá)到8 mm左右，電磁開關(guān)系統(tǒng)的間隙較大，符合工況要求。當(dāng)電磁開關(guān)系統(tǒng)的間隙小于8 mm時(shí)，不符合工況要求，可以通過修改系統(tǒng)參數(shù)減小振幅，也可以在對電磁開關(guān)系統(tǒng)進(jìn)行參數(shù)設(shè)計(jì)時(shí)避開固有頻率與外激頻率相等以減小振幅。

圖3至圖7為電磁開關(guān)系統(tǒng)主共振的均方響應(yīng)曲線。其中實(shí)線部分為穩(wěn)定解，虛線部分為不穩(wěn)定解。

對比圖2和圖3，可得窄帶隨機(jī)激勵(lì)下主共振系統(tǒng)均方響應(yīng)曲線與γ=0時(shí)主共振系統(tǒng)幅頻響應(yīng)曲線的拓?fù)浣Y(jié)構(gòu)相同。

圖3為電磁開關(guān)系統(tǒng)主共振隨阻尼系數(shù)rv變化的均方響應(yīng)曲線，由圖3可知，增大線圈與磁鐵之間的摩擦阻力系數(shù)，均方值減小。

圖4、圖5為電磁開關(guān)系統(tǒng)主共振隨激勵(lì)電壓U0和Um變化的均方響應(yīng)曲線，由圖4、圖5可知，隨著激勵(lì)電壓U0、Um的增大，均方值增大。

圖6、圖7為電磁開關(guān)系統(tǒng)主共振隨電阻r和R變化的均方響應(yīng)曲線，由圖6、圖7可知，隨著電阻r和R的增大，均方值減小。

圖2 幅頻響應(yīng)曲線

圖3 均方響應(yīng)曲線

圖4 均方響應(yīng)曲線

圖5 均方響應(yīng)曲線

圖6 均方響應(yīng)曲線

圖7 均方響應(yīng)曲線

4 主共振時(shí)間歷程分析

用4階龍格-庫塔法計(jì)算出系統(tǒng)（8）的時(shí)間響應(yīng)曲線和相圖。

圖8－圖11為γ=0.01、γ=0.05時(shí)電磁開關(guān)系統(tǒng)主共振時(shí)間響應(yīng)曲線。圖9，11分別為γ=0.01、γ=0.05時(shí)電磁開關(guān)系統(tǒng)主共振相圖。對比圖8和10，圖9和11可知，當(dāng)γ變化時(shí)電磁開關(guān)系統(tǒng)的時(shí)間響應(yīng)曲線和極限環(huán)變化不大。

5 結(jié)語

根據(jù)拉格朗日-麥克斯韋方程建立電磁開關(guān)系統(tǒng)動(dòng)力學(xué)模型。應(yīng)用改進(jìn)多尺度法研究電磁開關(guān)強(qiáng)非線性系統(tǒng)主共振在窄帶隨機(jī)激勵(lì)下的響應(yīng)，得到系統(tǒng)主共振的幅頻響應(yīng)方程和隨機(jī)均方響應(yīng)的近似表達(dá)式，考慮了隨機(jī)項(xiàng)對電磁開關(guān)系統(tǒng)響應(yīng)的影響。結(jié)果表明，當(dāng)系統(tǒng)產(chǎn)生共振時(shí)振幅達(dá)到8 mm左右，當(dāng)電磁開關(guān)系統(tǒng)的間隙較大時(shí)符合工況要求。當(dāng)電磁開關(guān)系統(tǒng)的間隙小于8 mm時(shí)，不符合工況要求，可以通過修改系統(tǒng)參數(shù)減小振幅，也可以在對電磁開關(guān)系統(tǒng)進(jìn)行參數(shù)設(shè)計(jì)時(shí)使固有頻率與外激頻率避開以減小振幅。隨著阻尼系數(shù)的增大，系統(tǒng)的均方值減??；隨著激勵(lì)電壓幅值增大，系統(tǒng)的均方值增大；隨著電阻增大，系統(tǒng)的均方值減小。通過改變電磁開關(guān)系統(tǒng)參數(shù)，得到系統(tǒng)均方響應(yīng)曲線的變化規(guī)律，找到了振動(dòng)的控制策略。

圖8 γ=0.01時(shí)間響應(yīng)曲線

圖9 γ=0.01的相圖

圖10 γ=0.05時(shí)間響應(yīng)曲線

圖11 γ=0.05的相圖圖9

[1]于成偉，馬秀艷.起動(dòng)機(jī)用電磁開關(guān)[J].汽車電器，2009 (2)：29-33.

[2]張勝川，邊建華.汽車起動(dòng)機(jī)電磁開關(guān)觸片觸點(diǎn)燒蝕原因探析[J].汽車電器，2010(4)：28-30.

[3]黃文鋒.起動(dòng)機(jī)電磁開關(guān)的結(jié)構(gòu)改進(jìn)[J].機(jī)電工程技術(shù)，2014(8)：130-131.

[4]楊志安，催一輝.非線性電阻電感型RLC串聯(lián)電路主共振分析[J].天津大學(xué)學(xué)報(bào)，2007，40(5)：579-583.

[5]楊志安，李熙，孟佳佳.改進(jìn)多尺度法求解環(huán)形極板機(jī)電耦合強(qiáng)非線性系統(tǒng)主共振的研究[J].振動(dòng)與沖擊，2015，34(19)：208-212.

[6]戎海武，孟光，王向東，等.窄帶隨機(jī)噪聲作用下強(qiáng)非線性系統(tǒng)的響應(yīng)[J].振動(dòng)工程學(xué)報(bào)，2003，16(1)：81-84.

[7]關(guān)卓懷，蔡敢為.轉(zhuǎn)子彎扭耦合振動(dòng)共振特征分析[J].噪聲與振動(dòng)控制，2016，36(1)：5-10.

[8]韓修品，蔣豐.具有非線性減振器的彈性雙層隔振系統(tǒng)建模[J].噪聲與振動(dòng)控制，2016，36(2)：31-34.

[9]DI LIUA,JING LIB,YONG XUC.Principal resonance responsesofSDOF systemswith smallfractional derivative damping under narrow-band random parametric excitation[J].Commun Nonlinear,2014,19(10):3642-3652.

[10]FENG Z H,LAN X J,ZHU X D.Explanation on the importance of narrow-band random excitation characters in the response of a cantilever beam[J].Journal of Sound and Vibration.2009,325(4):923-937.

[11]YANG ZHIAN,LI WENLAN,QIU JIAJUN.Lagrangemaxwellequationsand parameterresonanceofthe generator under magnetic saturation[J].Applied Mathematics and Mechanics,2007,28(11):1379-1386.

[12]DAVIS H G,LIU Q.The response probability density function of a Duffing oscillator with random narrow band excitation[J].Journal of Sound and Vibration,1990, 139:1-8.

[13]HUANG Z L,ZHU W Q,NI Y Q,et al.Stochastic averaging of strongly non-linear oscillators under bound-ed noise excitation[J].Journal of Sound and Vibration,2002,25（4）:245-267.

[14]朱位秋.隨機(jī)振動(dòng)[M].北京：科學(xué)出版社，1992.

Analysis of Primary Resonance of Nonlinear Vibration of Solenoid Switch Systems under Narrow-band Excitation

YANG Zhi-an1,ZHAO Li-sha2
(1.Key Laboratory of Structure and Vibration Engineering of Tangshan,Tangshan College, Tangshan 063000,Hebei China; 2.School of Mechanical Engineering,North China University of Technology, Tangshan 063000,Hebei China)

The strong nonlinear vibration of the electric circuit and magnetic circuit coupled solenoid switch system under narrow-band excitation is investigated.The stochastic differential equation of the solenoid switch system subjected to narrow-band excitation is established.The frequency response equation of the system is obtained based on the modified multi-scale method.The stability of the primary resonance of the system is calculated.And the influence of the parameters of the system on the mean square value of the primary resonance of the system is analyzed.The results show that the mean square value of the primary resonance of the system reduces with the increasing of the damping coefficient,increases with the increasing of the excitation voltage amplitude,and reduces with the increasing of the resistance.The change of limit cycle of the system is small with the increasing of the intensity of the random disturbance.

vibration and wave;solenoid switch;the modified method of multiple scales;narrow-band random excitation;mean square value;strongly nonlinear vibration

0324

：ADOI編碼：10.3969/j.issn.1006-1335.2017.01.007

1006-1355(2017)01-0030-05

2016-05-05

河北省自然科學(xué)基金項(xiàng)目資助（A200900097）

楊志安（1963－），男，河北省秦皇島市人，博士，教授，碩士生導(dǎo)師，研究方向?yàn)闄C(jī)構(gòu)學(xué)與機(jī)械動(dòng)力學(xué)。E-mail:yangzhi_an@163.com

趙利沙（1988－），女，河北省邯鄲市人，碩士研究生，研究方向?yàn)闄C(jī)構(gòu)學(xué)與機(jī)械動(dòng)力學(xué)。E-mail:948281620@qq.com