對稱間隙單自由度振動系統(tǒng)的等效電路仿真和實驗

汪 諍, 潘麗華, 劉文輝, 包發(fā)毅

(蘭州交通大學(xué) 機(jī)電工程學(xué)院，蘭州 730070)

對稱間隙單自由度振動系統(tǒng)的等效電路仿真和實驗

汪 諍, 潘麗華, 劉文輝, 包發(fā)毅

(蘭州交通大學(xué) 機(jī)電工程學(xué)院，蘭州 730070)

針對對稱間隙單自由度振動系統(tǒng)，為了解決仿真速度慢、實驗設(shè)計難度大的問題，設(shè)計了一種在數(shù)學(xué)模型上完全等效的電子電路，該電路僅使用了運放、電阻和電容等簡單的電子器件。通過MultiSim12.0軟件進(jìn)行仿真，結(jié)果和數(shù)值計算完全吻合。在仿真的基礎(chǔ)上，通過對非線性模塊的設(shè)計，搭建了硬件PCB電路，應(yīng)用信號發(fā)生器和示波器進(jìn)行了實驗，進(jìn)一步驗證了應(yīng)用等效電子電路進(jìn)行非線性振動實驗的合理性和正確性，并且具有數(shù)值仿真無可比擬的高速運算和動態(tài)參數(shù)調(diào)整的優(yōu)良特性，為非線性系統(tǒng)的仿真和實驗提供了一種參考。

振動系統(tǒng)；等效電子電路；仿真；實驗

含有分段線性、間隙、干摩擦和時滯等因素的非光滑機(jī)械系統(tǒng)是近年來機(jī)械工程領(lǐng)域研究的重點和熱點，該類系統(tǒng)在未接觸到邊界或未發(fā)生碰撞時，呈現(xiàn)出線性，當(dāng)接觸到邊界或出現(xiàn)干摩擦、碰撞等非線性因素時，將會引起噪聲和強(qiáng)烈的振動，伴隨著對接觸的部件造成嚴(yán)重的磨損和沖擊。這些都是由于非線性因素所引起，而線性系統(tǒng)是無法解釋的，如鞍結(jié)分岔、hopf分岔、同宿波分岔、倍周期分岔和其它局部非線性因素等，同時顯示非連續(xù)系統(tǒng)的非線性因素造成了很多非標(biāo)準(zhǔn)特性的發(fā)生，這就需要用復(fù)雜的非線性動力學(xué)來解釋[1-3]。

非線性機(jī)械振動的研究通常以理論研究為基礎(chǔ)，通過數(shù)值仿真來進(jìn)行驗證。數(shù)值仿真的缺點是受到計算機(jī)運算速度的限制，仿真時間漫長。在具體實驗方面，受到多個方面的約束，進(jìn)展很慢。論文設(shè)計了一種和含有對稱間隙的單自由度振動系統(tǒng)完全等效的模擬電子電路[4-5]，在MultiSim12.0軟件下進(jìn)行仿真，仿真結(jié)果和數(shù)值仿真結(jié)果完全相同，在此基礎(chǔ)上，用TL074CN運算放大器、電阻和電容等常用電子元件搭建了硬件PCB電路，應(yīng)用信號發(fā)生器和示波器完成了電路實驗，實驗結(jié)果與仿真結(jié)果完全等效。模擬電子電路不僅具有超高的運行速度，而且具有間隙、激勵頻率動態(tài)調(diào)整的特性，實驗表明，應(yīng)用等效模擬電子電路進(jìn)行非線性系統(tǒng)的仿真實驗具有正確性、高速性和動態(tài)調(diào)整的優(yōu)點。

1 物理模型和數(shù)學(xué)模型

典型的單自由度對稱間隙碰撞模型如圖1所示。質(zhì)量為M的振子通過剛度系數(shù)為K0的線性彈簧和阻尼系數(shù)為C的線性阻尼器相連接，振子在簡諧激振力的作用下沿著無摩擦的水平面運動。設(shè)線性彈簧在完全釋放狀態(tài)時，簡諧激振力為零點，以零點建立向右的一維坐標(biāo)系統(tǒng)，在距離零點位置為B的左右兩邊各有一個剛度系數(shù)為K1的線性彈簧。在激振力的作用下，振子M在運動位移小于間隙B時做簡諧振動，當(dāng)大于或等于B時，將會和K1作用，相互碰撞[6-7]。

圖1 單自由度間隙碰撞模型Fig.1 Single degree of freedom collision model with gap

1.1 數(shù)學(xué)模型

圖1所示對稱間隙單自由度振動系統(tǒng)為典型的非線性振動系統(tǒng)，主要含有分段線性彈性力F(X)，其數(shù)學(xué)模型為

(1)

非線性部分，即分段線性彈性力為

取無量綱化參數(shù)為

則方程的無量綱化形式為

(2)

式中

1.2 分岔圖

式(2)中有ζ、μ和δ三個參數(shù)，選擇不同的參數(shù)會出現(xiàn)不同的運動特性。在此取μ=30，ζ=0. 2，δ=0.04時對系統(tǒng)進(jìn)行數(shù)值計算，得到系統(tǒng)的分岔圖如圖2所示，其中橫坐標(biāo)為激振頻率ω，縱坐標(biāo)為瞬時速度，激振頻率選擇在1.82～2.02之間，從分岔圖中可以看出，隨著激振頻率的變化，系統(tǒng)出現(xiàn)周期運動、倍化分岔，周期三分岔和混沌等特性。

圖2 分岔圖Fig.2 Bifurcation diagram

圖3 系統(tǒng)在不同激勵頻率下的相圖Fig.3The phase diagram of system on different excitation frequency

2 等效電路模型及仿真

振動系統(tǒng)的數(shù)值分析和計算均是用計算機(jī)在離散情況下進(jìn)行積分計算，通常使用龍哥庫塔算法，算法需要若干次的迭代運算，因此在速度和精度上存在很大的矛盾，特別是當(dāng)改變參數(shù)(或者初值)，運算又得重新開始，在非線性系統(tǒng)的分析和研究中，嚴(yán)重的受到計算速度的制約，并且計算為靜態(tài)計算。

為了實現(xiàn)高速仿真和實驗，設(shè)計了一種和振動系統(tǒng)完全等效的模擬電路，該電路的微分方程和振動系統(tǒng)的微分方程完全相同，電路由運算放大器、電阻和電容組成，具有數(shù)值分析無法比擬的高速性、參數(shù)可實時變化的優(yōu)點。

2.1 等效電路設(shè)計

(3)

實現(xiàn)式(3)的電路有很多種，在此選用由TL074CN運算放大器、電容和電阻組成的積分電路、加法電路和反相器來實現(xiàn)式(3)，設(shè)計如圖4所示。電路采用了3個模塊部分，圖4(a)中的第一個(U1A和U1C)和第二個(U1B和U1D)模塊用來實現(xiàn)式(3)中的兩個方程，圖4(b)為第三個模塊(非線性模塊)，用來實現(xiàn)f(x1)，圖4中R1～R9分別為10 kΩ、100 kΩ、100 kΩ、43.3 kΩ、30 kΩ、25 kΩ、100 kΩ、100 kΩ和100 kΩ，C1和C2均為22 nF，供電電壓為±15 V。該電路在Multisim12.0中設(shè)計并完成了仿真。

圖4 等效電路及非線性模塊Fig.4 The equivalent circuit and nonlinear module

為了確定圖4電路中電容的容值和電阻的阻值，寫出該電路的數(shù)學(xué)模型

將式(4)轉(zhuǎn)換為二階微分方程，則有

(5)

式(5)在形式上和式(1)完全相同，為了完全等效，必須確定式(5)中各個參數(shù)的數(shù)值。根據(jù)μ=30，ζ=0. 2，δ=0.04進(jìn)行計算，首先取R1=10 kΩ，R7=100 kΩ，則得到R4=43.3 kΩ、R5=30 kΩ、R6=25 kΩ。運放U1C和U1D單元僅用于反相，因此選取R2=R3=R8=R9=100 kΩ。

正弦波激勵信號的電壓幅值和間隙δ相關(guān)，根據(jù)δ=0.04可得到

激勵頻率f和ω之間的關(guān)系為

得到f=417.673 4ω

2.2 非線性模塊設(shè)計[7]

非線性模塊在工作過程中的輸入和輸出必須符合式(2)中f(x1)的要求，即

(6)

根據(jù)式(6)的要求，可以設(shè)計2個飽和模塊進(jìn)行實現(xiàn)，第一個模塊在±δ處出現(xiàn)飽和，第二個模塊的飽和參數(shù)遠(yuǎn)大于±δ。由于運放在放大部分是線性放大，因此利用2個放大倍數(shù)完全相同且相位相反、飽和參數(shù)不同的模塊，通過加法電路后可以實現(xiàn)式(6)。

考慮到TL074CN的電源電壓為±18 V，飽和電壓約為±16 V，x11通過分壓電阻降低U4A與U4B組成的同相放大電路的放大倍數(shù)，并且與U4D組成的反相放大電路的放大倍數(shù)相同。也就是說x11和x12放大倍數(shù)相同，但相互反相。因此在U4A與U4B組成的同相放大電路飽和之前，線路輸出f(x)為0，而在其飽和后直至U1D組成的反相放大電路飽和之前，輸出f(x)為輸入x1的25倍。整體非線性模塊電路設(shè)計如圖5所示，圖中R30、R33、R34、R35和R36分別為20 kΩ、500 kΩ、4 kΩ、9 kΩ和1 kΩ，其余電阻均為100 kΩ。

圖5 非線性模塊電路Fig.5 Nonlinear module circuit

2.2.1 U4A和U4B

放大倍數(shù)為

輸出x11為

(7)

即

2.2.2 U4D

U4D組成的反相放大電路的放大倍數(shù)A2為

考慮到U4D的飽和輸入電壓遠(yuǎn)遠(yuǎn)高于x1的實際輸入最大電壓，因此只選取其放大部分，即當(dāng)輸入x1滿足-0.64 V≤x1≤0.64 V時，輸出x12為

(8)

2.2.3 U4C

U4C的作用是將x11和x12相加并且反向，即將式(7)和式(8)相加，當(dāng)R38=R39=R40時，輸出f(x1)為

(9)

由于飽和電壓約為±16 V，所以飽和輸入電壓約為

式(9)可以簡化為

該式完全滿足式(9)的要求，其中間隙δ為0.064 V。在MultiSim12.0軟件下仿真的輸入輸出特性如圖6所示(橫軸為CH1，50 mV/Div，縱軸為CH2，2 V/Div)。

圖6 非線性電路輸入輸出特性仿真Fig.6 The coordinate figure of input and output characteristics in nonlinear simulation module circuit

2.3 Multisim12.0仿真

3 硬件電子電路實驗

圖7 MultiSim12.0軟件中仿真的相圖Fig.7 The simulation phase diagram use software of MultiSim12.0

圖8 硬件PCB電路和測試圖Fig.8 The photograph both PCB circuit and test

圖9 不同頻率下示波器輸出的相圖Fig.9 The phase diagram photograph on different frequency by oscilloscope testing

4 結(jié) 論

論文以含有對稱間隙的單自由度振動系統(tǒng)模型為對照模型，應(yīng)用TL074CN運算放大器、電容和電阻，設(shè)計了完全等效的電子電路模型。電路中的非線性模塊應(yīng)用了運算放大器的飽和特性，通過壓縮、反相和加法電路實現(xiàn)了分段線性特性。對電路模型在MultiSim12.0軟件中進(jìn)行仿真，并根據(jù)電路原理搭建了硬件PCB電路，實現(xiàn)了等效硬件電路的實驗過程。

通過仿真和實驗，表明使用模擬電路可以完全等效非線性動力學(xué)模型，模擬電路不僅具有運算速度快的特點，最主要的是可以方便的調(diào)整間隙、激勵頻率等參數(shù)實現(xiàn)動態(tài)的仿真和實驗。仿真和實驗結(jié)果與數(shù)值仿真結(jié)果完全符合，進(jìn)一步驗證了等效電路對于動力學(xué)實驗的可靠性和等同性。類似的實驗可以應(yīng)用與生物、化學(xué)和經(jīng)濟(jì)等不同學(xué)科的非線性模型中，為非線性動力學(xué)的仿真和實驗提供了一種高效快速、參數(shù)實時調(diào)整的仿真和實驗方法。

[1] SHAW S W, HOLMES P J. A periodically forced piecewise linear oscillator[J]. Journal of Sound and Vibration, 1983,90(1): 129-155.

[2] PETERKA F,VACIK J. Transition to chaotic motion in mechanical systems with impacts[J]. Journal of Sound and Vibration,1992,154(1): 95-115.

[4] PEHLIVAN I, MOROZ I M, VAIDYANATHAN S. Analysis, synchronization and circuit design of a novel butterfly attractor[J]. Journal of Sound and Vibration,2014,333:5077-5096.

[5] RUIZ P, GUTIERREZ J M, GUEMEZ J. Experimental mastering of nonlinear dynamics in circuits by sporadic pulses[J]. Chaos, Solitons and Fractals,2008,36: 635-645.

[6] 李萬祥,牛衛(wèi)中.一類含間隙系統(tǒng)的分岔與混沌的形成過程[J]. 振動與沖擊,2005,24(3):47-49． LI Wangxiang, NIU Weizhong. Bifurcation and chaos of a system with a pair of symmetric set-up elastic stops[J]. Journal of Vibration and Shock, 2005,24(3):47-49.

[7] 張有強(qiáng),丁旺才,孫闖.單自由度含間隙和干摩擦碰撞振動系統(tǒng)的分岔與混沌[J].振動與沖擊, 2008,27(7):102-130． ZHANG Youqiang, DING Wangcai, SUN Chuang. Bifurcation and chaos of one degree of freedom impact vibration system with clearance and dry friction[J]. Journal of Vibration and Shock, 2008,27(7):102-130.

Equivalent circuits simulation and tests for a single-DOF vibration system with symmetrical gap

WANG Zheng, PAN Lihua, LIU Wenhui, BAO Fayi

(School of Mechatronic Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China)

Based on a single-DOF vibration system with symmetric gap, in order to solve problems of slow simulation speed, and difficult to design tests, an electronic circuit with a full equivalent mathematical model was designed. The circuit consisted of simple electronic components, such as, operational amplifier’s, resistances and capacitances. With the software MultiSim12.0, the simulation was done and the simulation results agreed well with those of numerical calculation. On the basis of simulation and designing nonlinear modules, the circuit board of PCB was made. Tests were conducted with devices of signal generator and oscillograph. The rationality and correctness of nonlinear vibration tests were further verified with the equivalent electronic circuit. It was shown that the circuit has excellent characteristics in high speed computing and dynamic parameter adjustment, they are incomparably beyond the characteristics of numerical simulation. The results provided a reference for simulation and tests of nonlinear systems.

vibration system; equivalent electronic circuit; simulation; test

2015-09-07 修改稿收到日期：2015-12-31

汪諍 男，碩士，副教授，1970年生

N33

A

10.13465/j.cnki.jvs.2017.01.021