胡宇達(dá)， 王　彤

(燕山大學(xué) 河北省重型裝備與大型結(jié)構(gòu)力學(xué)可靠性重點(diǎn)實(shí)驗(yàn)室，河北 秦皇島　066004)

磁場中導(dǎo)電旋轉(zhuǎn)圓板的磁彈性非線性共振

胡宇達(dá)， 王彤

(燕山大學(xué) 河北省重型裝備與大型結(jié)構(gòu)力學(xué)可靠性重點(diǎn)實(shí)驗(yàn)室，河北 秦皇島066004)

研究旋轉(zhuǎn)運(yùn)動(dòng)圓形薄板在磁場中受到機(jī)械載荷作用時(shí)的磁彈性非線性共振問題。根據(jù)哈密頓原理推導(dǎo)出旋轉(zhuǎn)運(yùn)動(dòng)圓板在磁場中的磁彈性非線性振動(dòng)方程，基于電磁理論給出了旋轉(zhuǎn)板所受電磁力的表達(dá)式;通過位移函數(shù)的設(shè)定并應(yīng)用伽遼金積分法，得到橫向磁場中旋轉(zhuǎn)導(dǎo)電圓板的磁彈性軸對稱振動(dòng)微分方程。應(yīng)用平均法對系統(tǒng)非線性主共振問題進(jìn)行求解，得到穩(wěn)態(tài)運(yùn)動(dòng)下的幅頻響應(yīng)方程。通過數(shù)值計(jì)算，得到固支邊界條件下圓板的幅頻特性曲線以及振幅隨磁感應(yīng)強(qiáng)度、轉(zhuǎn)速、激勵(lì)力等參數(shù)的變化規(guī)律曲線圖，分析了不同參數(shù)對旋轉(zhuǎn)板共振幅值及非線性特性的影響。

導(dǎo)電圓板；旋轉(zhuǎn)運(yùn)動(dòng)；主共振；磁場；平均法

在現(xiàn)代工程實(shí)際中，旋轉(zhuǎn)運(yùn)動(dòng)和電磁結(jié)構(gòu)廣泛存在于航空航天、機(jī)械、土木等領(lǐng)域中的高端設(shè)備及運(yùn)動(dòng)結(jié)構(gòu)中。處于電磁場環(huán)境中的磁彈性構(gòu)件，因受到電磁場和機(jī)械場等因素的作用，往往會(huì)產(chǎn)生較為明顯的振動(dòng)，并影響著系統(tǒng)的正常運(yùn)行。因此，對電磁場中旋轉(zhuǎn)運(yùn)動(dòng)結(jié)構(gòu)的磁彈性振動(dòng)問題進(jìn)行研究具有理論和實(shí)際意義。近年來，國內(nèi)外學(xué)者針對圓板及磁彈性板的振動(dòng)做了許多有意義的理論研究工作。文獻(xiàn)[1-3]分別用微分變換法、切比雪夫里茨法和漸進(jìn)法對圓板的非線性振動(dòng)進(jìn)行了研究與計(jì)算。文獻(xiàn)[4]分析了導(dǎo)電板的電磁熱機(jī)械行為，采取用于渦流分析的T方法得到了導(dǎo)電圓板的解析解。文獻(xiàn)[5-7]研究了導(dǎo)電薄板在磁場中的組合共振和諧波共振等問題，建立了軸向運(yùn)動(dòng)導(dǎo)電板的磁熱彈性耦合動(dòng)力學(xué)理論模型。文獻(xiàn)[8]研究了電磁場環(huán)境下導(dǎo)電圓板的磁彈性強(qiáng)迫振動(dòng)問題。文獻(xiàn)[9-11]研究了旋轉(zhuǎn)運(yùn)動(dòng)板的非線性振動(dòng)，得到了不同條件下固有頻率的計(jì)算式。文獻(xiàn)[12]用半解析法分析了功能梯度圓薄板的非線性自由振動(dòng)與強(qiáng)迫振動(dòng)。文獻(xiàn)[13]運(yùn)用有限元方法分析了旋轉(zhuǎn)夾層板的非軸對稱振動(dòng)與穩(wěn)定問題，討論了不同參數(shù)對振動(dòng)的影響。文獻(xiàn)[14]針對含有黏彈性夾芯層旋轉(zhuǎn)圓板在氣動(dòng)載荷作用下的行波動(dòng)力學(xué)及穩(wěn)定性問題進(jìn)行了研究。文獻(xiàn)[15]建立了高速旋轉(zhuǎn)柔性矩形薄板的耦合動(dòng)力學(xué)模型，并應(yīng)用模態(tài)截?cái)喾ㄟM(jìn)行了求解。

本文研究旋轉(zhuǎn)運(yùn)動(dòng)圓板在磁場中的非線性共振問題，得到磁彈性非線性振動(dòng)方程和幅頻特性方程，給出共振幅值隨不同參數(shù)變化規(guī)律曲線圖，分析磁場、轉(zhuǎn)速、激勵(lì)力等參數(shù)的影響。

1旋轉(zhuǎn)運(yùn)動(dòng)圓板磁彈性振動(dòng)方程

1.1動(dòng)能與變形能

對于處于磁場中做旋轉(zhuǎn)運(yùn)動(dòng)的圓形薄板，設(shè)薄板內(nèi)任意一點(diǎn)的位移矢量在柱坐標(biāo)系三個(gè)方向的分量為：

ur1(r,θ,z,t)=ur(r,θ,t)+zu1(r,θ,t)

(1)

uθ1(r,θ,z,t)=uθ(r,θ,t)+zv1(r,θ,t)

(2)

uz(r,θ,z,t)=w(r,θ,t)

(3)

在式(1),式(2),式(3)基礎(chǔ)上，通過求導(dǎo)運(yùn)算，得到旋轉(zhuǎn)運(yùn)動(dòng)薄板的速度分量表達(dá)式：

(4)

式中：Ω為旋轉(zhuǎn)速度。進(jìn)而得到旋轉(zhuǎn)運(yùn)動(dòng)圓板動(dòng)能的表達(dá)式為：

(5)

旋轉(zhuǎn)運(yùn)動(dòng)薄板由彎曲變形引起的形變勢能Uε1表達(dá)式為：

(6)

同時(shí)，旋轉(zhuǎn)運(yùn)動(dòng)薄板的中面應(yīng)變勢能為；

(7)

式中：Mr、Mθ為彎矩，Mrθ為扭矩，κr、κθ為曲率，κrθ為扭率，Nr、Nθ、Nrθ為中面內(nèi)力，εr、εθ、εrθ為中面應(yīng)變。

1.2外力所作功

假設(shè)旋轉(zhuǎn)運(yùn)動(dòng)薄板的橫向位移w發(fā)生了位移邊界條件所容許的微小變化量，即虛位移δw，則外加橫向強(qiáng)迫力Pz在虛位移上所作的虛功可以表示為：

(8)

同樣，電磁力Fr,Fθ,Fz與電磁力矩mr，mθ在虛位移上所作虛功為如下形式：

δUF=

(9)

1.3哈密頓原理建立振動(dòng)方程

根據(jù)哈密頓原理有：

(10)

式中：t0和t1為兩個(gè)固定時(shí)刻。

將式(5)、式(6)、式(7)、式(8)及式(9)代入式(10)中，并考慮幾何非線性情況及軸對稱振動(dòng)問題，經(jīng)過變分運(yùn)算，推得旋轉(zhuǎn)運(yùn)動(dòng)圓形薄板的磁彈性軸對稱振動(dòng)方程：

(11)

(12)

2電磁力表達(dá)式

在磁場中做旋轉(zhuǎn)運(yùn)動(dòng)的薄板將受到Lorenz力的作用，薄板單位體積內(nèi)電磁力表達(dá)式為：

f(fr,fθ,fz)=J×B0=fri+fθj+fzk

(13)

式中：J=σ0(e+V×B0)，e為電場強(qiáng)度矢量，V為板內(nèi)各點(diǎn)速度矢量，σ0為電導(dǎo)率。

對式(13)沿板厚方向進(jìn)行積分，可得到式(11),式(12)中薄板軸對稱振動(dòng)時(shí)所受電磁力和力矩的表達(dá)式：

(14)

(15)

(16)

3橫向磁場中旋轉(zhuǎn)板的軸對稱主共振

針對在恒定橫向磁場(即B0r=0、B0θ=0)環(huán)境中做勻速旋轉(zhuǎn)運(yùn)動(dòng)的圓形薄板，考慮其軸對稱振動(dòng)情況。將相應(yīng)彎矩內(nèi)力、中面內(nèi)力和電磁力表達(dá)式代入振動(dòng)方程中，并忽略縱向位移對薄板橫向振動(dòng)的影響(即ur=0)。由于在橫向磁場中Fz=0，得出旋轉(zhuǎn)運(yùn)動(dòng)薄板在橫向磁場中關(guān)于撓度w的磁彈性振動(dòng)方程：

(17)

周邊夾支圓板的邊界條件為：

周邊簡支圓板的邊界條件為：

r=R時(shí)，w=0，Mr=0

設(shè)滿足邊界條件的關(guān)于式(17)的動(dòng)位移解為：

w(r,t)=W(r)T0(t)=

(18)

將式(18)代入式(17)，應(yīng)用伽遼金積分法推得系統(tǒng)的非線性強(qiáng)迫振動(dòng)微分方程：

(19)

式中：“撇號(hào)”表示對時(shí)間t的一階和二階導(dǎo)數(shù)(下同)；

(20)

應(yīng)用平均法對方程式(20)進(jìn)行解析求解，并將系統(tǒng)的解及其導(dǎo)數(shù)分別取為：

q=acos(ωt-θ)

(21)

q′=-aωsin(ωt-θ)

(22)

這里，a和θ為關(guān)于時(shí)間t的慢變函數(shù)。

由式(21)，式(22)可得：

a′cosψ+aθ′sinψ=0

(23)

將式(21)，式(22)代入式(20)中，經(jīng)過運(yùn)算得到：

(24)

聯(lián)立式(23)，式(24)，考慮a和θ的慢變規(guī)律，應(yīng)用平均法得到系統(tǒng)的平均化方程：

(25)

(26)

由式(25)，式(26)可得穩(wěn)態(tài)運(yùn)動(dòng)下主共振的幅頻響應(yīng)方程：

(27)

4算例分析

對于橫向磁場中法向周期動(dòng)載作用的鋁制薄板，給定物理參數(shù)：密度ρ=2 670 kg/m3，電導(dǎo)率σ0=3.63×107(Ω·m)-1，泊松系數(shù)μ=0.34，彈性模量E=71 GPa。

圖1、圖2分別給出了不同旋轉(zhuǎn)速度和不同板厚條件下旋轉(zhuǎn)板的幅頻響應(yīng)特性曲線圖。圖中曲線表明，在共振區(qū)域(εσ≈0附近)振幅明顯增大且出現(xiàn)解的多值性；同時(shí)由兩圖可看出，不同轉(zhuǎn)速或不同厚度所對應(yīng)的曲線之間存在多個(gè)相交點(diǎn)，表明隨著轉(zhuǎn)速或板厚的增大，共振幅值既有隨之增大的區(qū)域，也有隨之減小的區(qū)域，非線性特征非常明顯。圖3、圖4分別給出了不同磁感應(yīng)強(qiáng)度和激勵(lì)力幅值條件下的幅頻響應(yīng)特性曲線圖。圖中曲線也呈現(xiàn)了共振區(qū)域幅值明顯增大和解的多值性，同時(shí)表明，隨磁感應(yīng)強(qiáng)度的增大和激勵(lì)力的減小，共振曲線呈現(xiàn)明顯內(nèi)縮趨勢。

圖1　速度變化時(shí)幅頻特性曲線圖Fig.1 The charactertic curve of amplitude-frequency with different velocity(B0z=3 T,R=0.8 m,h=0.015 m,q0=60 N/m2)

圖2　板厚變化時(shí)幅頻特性曲線圖Fig.2 The charactertic curve of amplitude-frequency with different thickness(q0=60 N/m2,B0z=1 T,R=0.3 m,Ω=1 200 r/min)

圖5、圖6分別給出了不同調(diào)諧參數(shù)和不同激勵(lì)力條件下振幅隨磁感應(yīng)強(qiáng)度的變化曲線圖。圖中曲線均呈現(xiàn)相對磁感應(yīng)強(qiáng)度B0z=0處的左右對稱形式，且在B0z=0附近共振被激發(fā)，幅值較大并有多值性。當(dāng)磁感應(yīng)強(qiáng)度增大到一定值時(shí)，振幅明顯減小，多值性也會(huì)消失。同時(shí)，由圖5可知，隨著調(diào)諧參數(shù)值εσ的不斷增大，共振區(qū)域幅值呈現(xiàn)從單值過渡到三個(gè)值，并且共振曲線逐漸對稱內(nèi)縮并最終分離出上部的封閉曲線的

變化規(guī)律，其中εσ=0.03和εσ=0.037為對應(yīng)兩個(gè)臨界狀態(tài)下的臨界值。圖6中曲線也具有隨激勵(lì)力減小，曲線呈現(xiàn)對稱內(nèi)縮且共振幅值減小的變化規(guī)律。

圖3　磁場變化時(shí)幅頻特性曲線圖Fig.3 The charactertic curve of amplitude-frequency with different magnetic induction intensity(q0=60 N/m2,Ω=1 200 r/min,R=0.3 m,h=0.004 m)

圖4　激勵(lì)力變化時(shí)幅頻特性曲線圖Fig.4 The charactertic curve of amplitude-frequency with different exitation force(B0z=1 T,R=0.3 m,Ω=1 200 r/min,h=0.004 m)

圖7、圖8分別給出了不同磁感應(yīng)強(qiáng)度和不同轉(zhuǎn)速條件下振幅隨激勵(lì)力幅值的變化曲線圖。圖中曲線表明，較小的激勵(lì)力幅值即可激發(fā)系統(tǒng)的主共振現(xiàn)象，并且首先會(huì)出現(xiàn)多值解，之后激勵(lì)力增大到一定值后，解退化為較大的單值。同時(shí)，由圖7可見，在單值區(qū)域內(nèi)，磁感應(yīng)強(qiáng)度越大，共振振幅越小。由圖8可知，隨轉(zhuǎn)速的增大，在力幅小的單值區(qū)域內(nèi)，共振振幅呈增大變化趨勢，而在力幅大的單值區(qū)域內(nèi)，共振振幅呈減小變化趨勢。

圖5　調(diào)諧參數(shù)變化時(shí)振幅-磁感應(yīng)強(qiáng)度特性曲線圖Fig.5 The charactertic curve of amplitude-magnetic induction intensity with different tuning parameters(q0=60 N/m2,R=0.3 m,Ω=1 200 r/min,h=0.004 m)

圖6 激勵(lì)力變化時(shí)振幅-磁感應(yīng)強(qiáng)度特性曲線圖Fig.6Thecharacterticcurveofamplitude-magneticinductionintensitywithdifferentexitationforce(εσ=0.03,R=0.3m,Ω=1200r/min,h=0.004m)圖7 磁感應(yīng)強(qiáng)度變化時(shí)振幅-激勵(lì)力特性曲線圖Fig.7Thecharacterticcurveofamplitude-exitationforcewithdifferentmagneticinductionintensity(εσ=0.02,R=0.3m,Ω=1200r/min,h=0.004m)圖8 速度變化時(shí)振幅-激勵(lì)力特性曲線圖Fig.8Thecharacterticcurveofamplitude-exitationforcewithdifferentvelocity(εσ=0.4,R=0.8m,B0z=3T,h=0.015m)

5結(jié)論

本文針對磁場中旋轉(zhuǎn)運(yùn)動(dòng)圓形板，推導(dǎo)出了其磁彈性振動(dòng)方程，得到了圓板的幅頻響應(yīng)方程。通過數(shù)值算例，對夾支邊界條件圓板的共振問題進(jìn)行了分析，結(jié)果表明：

(1)共振區(qū)域振幅明顯增大并有多值性，板的旋轉(zhuǎn)速度、厚度大小等參數(shù)對共振曲線均有顯著影響。

(2)共振振幅隨激勵(lì)力力幅的增大呈增大趨勢，而激勵(lì)力頻率的改變直接影響著解的形態(tài)。

(3)磁感應(yīng)強(qiáng)度增大到一定數(shù)值時(shí)共振振幅明顯下降，即可通過調(diào)整外加磁感應(yīng)強(qiáng)度的大小控制系統(tǒng)的振動(dòng)。

[1] Yalcin H S, Arikoglu A, Ozkol I. Free vibration analysis of circular plates by differential transformation method [J]. Applied Mathematics and Computation, 2009, 212(2): 377-386.

[2] Zhou D, Au F T K, Cheung Y K, et al. Three-dimensional vibration analysis of circular and annular plates via the chebyshev-ritz method [J]. International Journal of Solids and Structures, 2003, 40(12): 3089-3105.

[3] Jeong K H. Free vibration of two identical circular plates coupled with bounded fluid [J]. Journal of Sound and Vibration, 2003, 260(4): 653-670.

[4] Gao Yuan-wen, Xu Bang, Huh H. Electromagneto-thermo-mechanical behaviors of conductive circular plate subject to time-dependent magnetic fields [J]. Acta Mechanics, 2010, 210(1/2): 99-116．

[5] Hu Yu-da, Li Jing. Magneto-elastic combination resonances analysis of current-conducting thin plate [J]. Applied Mathematic and Mechanics, 2008, 29(8): 1053-1066．

[6] Hu Yu-da, Li Jing. The magneto-elastic subharmonic resonances of current-conducting thin plate in magnetic field [J]. Journal of Sound and Vibration, 2009, 319(3/4/5): 1107-1120．

[7] 胡宇達(dá), 張金志. 軸向運(yùn)動(dòng)載流導(dǎo)電板磁熱彈性耦合動(dòng)力學(xué)方程[J]. 力學(xué)學(xué)報(bào), 2013, 45(5): 792-796.

HU Yu-da, ZHANG Jin-zhi. Magneto-thermo-elastic coupled dynamic equations of axially moving carry current plate in magnetic field [J]. Acta Mechanica Sinica, 2013, 45(5): 792-796.

[8] 胡宇達(dá), 杜國君. 磁場環(huán)境下導(dǎo)電圓形薄板的磁彈性強(qiáng)迫振動(dòng)[J]. 工程力學(xué), 2007, 24(7): 184-187.

HU Yu-da, DU Guo-jun. Forced vibration of a thin round conductive plate in magnetic field [J]. Engineering Mechanics, 2007, 24(7): 184-187.

[9] Hashemi S H, Farhadi S, Carra S. Free vibration analysis of rotating thick plates [J]. Journal of Sound and Vibration, 2009, 323(1/2): 366-384.

[10] Maretic R. Transverse vibration and stability of an eccentric rotating circular plate [J]. Journal of Sound and Vibration, 2005, 280(3/4/5): 467-478.

[11] Bauer H F, Eidel W. Transverse vibration and stability of spinning circular plates of constant thickness and different boundary conditions [J]. Journal of Sound and Vibration, 2007, 300(3/4/5): 877-895.

[12] Allahverdizadeh A, Naei M H, Nikkhah B M. Nonlinear free and forced vibration analysis of thin circular functionally graded plates [J]. Journal of Sound and Vibration, 2008, 310(4/5): 966-984.

[13] Chen Y R, Chen L W. Vibration and stability of rotating polar orthotropic sandwich annular plates with a viscoelastic core layer [J]. Composite Structures, 2007, 78(1): 45-57.[14] 李龍飛, 王省哲, 周又和. 黏彈性夾芯層合旋轉(zhuǎn)圓板的氣動(dòng)彈性動(dòng)力穩(wěn)定性分析[J]. 固體力學(xué)學(xué)報(bào), 2013, 34(4):325-332.

LI Long-fei, WANG Xing-zhe, ZHOU You-he. Analysis on aeroelastic dynamic stability of a rotating sandwich annular plate with viscoelastic core layer [J]. Acta Mechanic Solida Sinica, 2013, 34(4): 325-332.

[15] 趙飛云, 洪嘉振, 劉錦陽, 等. 高速旋轉(zhuǎn)柔性矩形薄板的動(dòng)力學(xué)建模和近似算法[J]. 振動(dòng)工程學(xué)報(bào), 2006, 19(3): 416-421.

ZHAO Fei-yun, HONG Jia-zhen, LIU Jin-yang, et al. Dynamic modeling and modal truncation approach for a high-speed rotating thin elastic rectangular plate[J]. Journal of Vibration Engineering, 2006, 19(3): 416-421.

[16] 徐芝綸. 彈性力學(xué)[M]. 北京: 高等教育出版社, 2006.

Nonlinear resonance of a conductive rotating circular plate in the magnetic field

HU Yu-da， WANG Tong

(Hebei Provincial Key Laboratory of Mechanical Reliability for Heavy Equipments and Large Structures，Yanshan University，Qinhuangdao 066004，China)

The nonlinear resonance of a conductive rotating thin circular plate subjected to mechanical loads in a magnetic field was investigated. The nonlinear vibration equation, related to the spinning round plate, was derived according to the Hamilton principle. Based on electromagnetic theory, the expressions of electromagnetic forces were derived. According to the set of a displacement function, the magnetoelastic forced vibration differential equation of the round plate was obtained through the application of the Galerkin integral method. By means of an averaging method, the amplitude-frequency response equation in a steady state was established. The amplitude frequency characteristic curves and the relationship curves of an amplitude that changes with the magnetic induction intensity, the speed of rotation and the excitation force of the plate with a fixed boundary condition were obtained according to the numerical calculation. The influence of different parameters on the amplitude and nonlinear characteristics of the spinning plate was analyzed finally.

conductive round plate; rotary motion; primary resonance; magnetic field; averaging method

10.13465/j.cnki.jvs.2016.12.028

國家自然科學(xué)基金(11472239)；河北省自然科學(xué)基金(A2015203023)；河北省高等學(xué)校科學(xué)技術(shù)研究重點(diǎn)項(xiàng)目(ZD20131055)

2014-05-20修改稿收到日期：2015-01-08

胡宇達(dá) 男，博士，教授，1968年11月生

O322;O442

A