蔣國慶， 李家文, 唐國金

(國防科學技術大學 航天科學與工程學院，長沙 410073)

基于矢量圖的力狀態(tài)映射法辨識參數(shù)誤差分析

蔣國慶， 李家文, 唐國金

(國防科學技術大學 航天科學與工程學院，長沙 410073)

力狀態(tài)映射法被常用于辨識非線性連接結(jié)構的等效動力學參數(shù)，該方法需要用到實驗儀器測量得到或者經(jīng)過微分或積分計算得到的激勵力、速度、位移、加速度等數(shù)據(jù)。而往往這些數(shù)據(jù)的相位會出現(xiàn)失真，這會直接導致通過力狀態(tài)映射辨識得到的系統(tǒng)參數(shù)出現(xiàn)誤差。用矢量圖法計算得到了單自由度系統(tǒng)在強迫振動和自由振動兩種情況下出現(xiàn)測量數(shù)據(jù)或計算數(shù)據(jù)小幅度相位失真時的誤差計算公式，結(jié)果表明這種小幅度相位失真對系統(tǒng)阻尼系數(shù)的影響明顯大于對系統(tǒng)剛度系數(shù)的影響，且當激勵力頻率等于固有頻率時辨識參數(shù)的誤差最小。

矢量圖；力狀態(tài)映射法；參數(shù)辨識；誤差分析

現(xiàn)代航天器大量使用連接結(jié)構，這些連接結(jié)構是系統(tǒng)非線性和無源被動阻尼的主要來源之一，連接結(jié)構線性模型常常不能反映其非線性特性，必須采用非線性模型，而非線性模型的參數(shù)辨識問題是必須解決的關鍵問題之一，一直以來都是學者們的研究熱點[1]。

力狀態(tài)映射法是一種用于辨識非線性結(jié)構的常用方法，該方法將連接結(jié)構的回復力表示為速度和位移的函數(shù)，能夠直接顯示系統(tǒng)的剛度特性和阻尼特性。該方法最初由MARSI等[2-3]提出，并經(jīng)CRAWLEY等[4-5]加以發(fā)展和利用。AL-HADID等[6]選用一般多項式對多自由度集總系統(tǒng)進行了辨識，發(fā)現(xiàn)一般多項式與正交多項式相比具有物理意義明確的優(yōu)點。蔡力鋼等[7]采用Tikhonov正則化方法和迭代算法對非線性結(jié)合部等效力學模型參數(shù)進行了辨識，辨識精度高。吳爽等[8]采用力狀態(tài)映射法對太陽翼板間鉸鏈結(jié)構非線性動力學模型進行了參數(shù)辨識，實驗測試結(jié)果和數(shù)值仿真結(jié)果取得了較好的一致性。

力狀態(tài)映射法是一種時域內(nèi)的方法，需要通過實驗測量得到必要的數(shù)據(jù)，然而，由于實驗儀器沒有校準或者對實驗數(shù)據(jù)進行積分或微分處理不恰當會導致各檢測量的相位出現(xiàn)失真，從而導致通過力狀態(tài)映射法辨識得到的系統(tǒng)參數(shù)出現(xiàn)誤差。WRIGHT等[9]獲得了用力狀態(tài)映射法辨識得到的正弦激勵載荷作用下線性系統(tǒng)參數(shù)的誤差計算表達式。MESKELL等[10]采用矢量圖法分別建立了單自由度系統(tǒng)在強迫振動和自由振動兩種情況下的模型，分析了激勵力、位移或者速度中某一項發(fā)生相位偏差時系統(tǒng)剛度參數(shù)和阻尼參數(shù)的誤差。但在實際情況下，可能是某兩項的相位發(fā)生偏差。本文采用矢量圖法推導得到了單自由度系統(tǒng)在強迫振動和自由振動兩種情況下某兩項發(fā)生相位偏差時辨識得到的系統(tǒng)參數(shù)誤差公式，誤差公式可以退化至Meskell和Fitzpatrick得到的公式。

1 強迫振動系統(tǒng)參數(shù)誤差分析

圖1為單自由度系統(tǒng)強迫振動示意圖，其中F為強迫激勵，其表達式為F=fsin(ωt)，m為質(zhì)量塊的質(zhì)量，k為彈簧剛度，c為阻尼系數(shù)。系統(tǒng)振動微分方程為

(1)

圖1 單自由度系統(tǒng)強迫振動示意圖Fig.1 Schematic of a linear single degree of freedom system

設y0=Asin(ωt-θ)為式(1)的特解，經(jīng)過換算可以將式(1)變?yōu)?/p>

(-mω2+jcω+k)ejωt=fejθejωt

(2)

k=fcosθ+mω2

(3)

c=fsinθ

(4)

mω2=kr2

(5)

由式(3)和(5)可以得到

f=(1-r2)k/cosθ

(6)

如果由于試驗時測量設備之間相位沒有校準等原因造成測量或者計算得到的各力的相位出現(xiàn)一定的失真，那么通過力狀態(tài)映射法辨識得到的剛度系數(shù)和阻尼系數(shù)便會出現(xiàn)誤差。文獻[10]中已經(jīng)利用矢量圖法計算得到了某一項的相位出現(xiàn)偏差時系統(tǒng)剛度系數(shù)和阻尼系數(shù)的誤差，但實際情況下不能確定發(fā)生相位偏差的力的個數(shù)，而可能是測量得到的激勵力、速度和位移兩兩出現(xiàn)相位偏差，分別如圖2(b)～(d)所示。當測量得到的三個力均出現(xiàn)相位偏差時，可以通過數(shù)據(jù)處理(如把所有的相位偏差均減去最小的相位偏差)使得其轉(zhuǎn)換為兩力出現(xiàn)相位偏差時的情況。

圖2 單自由度強迫振動系統(tǒng)矢量圖Fig.2 Vector diagrams of a linear signal degree of freedom system for forced response

圖2(b)為振動系統(tǒng)測量得到的激勵力和位移的相位均出現(xiàn)偏差時的矢量圖，其中激勵力相位偏差為τ，位移相位偏差為γ，kt和ct為中間剛度系數(shù)和阻尼系數(shù)，k′和c′為辨識得到的剛度系數(shù)和阻尼系數(shù)。

由ΔOAB和ΔOCD，可以得到

kt=k+fcos(θ-τ)-fcos(θ)

(7)

ct=fsin(θ-τ)

(8)

由ΔOEF，可以得到

k′=kt/cosγ

(9)

c′=ct+kttanγ

(10)

在此，定義參數(shù)誤差比

εk=(k′-k)/k

(11)

εc=(c′-c)/c

(12)

由式(6)、(7)、(9)和(11)可得

(13)

若γ=0，則

εk=(1-r2)(cosτ-1)+2ζrsinτ

(14)

若τ=0，則

(15)

將式(4)、(6)、(8)、(10)和(12)可以得到

εc=(cosτ-1)-((1-r2)/2ζr)sinτ+

(r/2ζ+sinτ+((1-r2)/2ζr)cosτ)tanγ

(16)

若γ=0，則

εc=(cosτ-1)-sinτ(1-r2)/2ζr

(17)

若τ=0，則

εc=tanγ/2ζr

(18)

式(14)、(15)、(17)(18)與文獻[10]中的結(jié)果一致。

同理可以得到其余兩種情況下系統(tǒng)剛度系數(shù)和阻尼系數(shù)的誤差公式，如表1所示。

表1 強迫振動系統(tǒng)剛度系數(shù)和阻尼系數(shù)誤差Tab.1 Stiffness errors and damper errors of signal degree of freedom for forced response

2 自由振動系統(tǒng)參數(shù)誤差分析

圖1中激勵力F=0時該系統(tǒng)即為自由振動系統(tǒng)，其振動微分方程為

(19)

經(jīng)過推導可以得到如下式子

(20)

從而有

(21)

根據(jù)式(21)可以得到圖3(a)所示矢量圖。

圖3 單自由度自由振動系統(tǒng)矢量圖Fig.3 Vector diagrams of a linear signal degree of freedom system for free response

3 算例分析

對于圖1中的強迫振動系統(tǒng)，取ζ=0.2，τ,γ的變化范圍均為[-4°,4°]，r分別取0.5、1、2，調(diào)用MATLAB對式(13)和(16)分別進行計算，計算結(jié)果如圖4～圖6所示。

表2 自由振動系統(tǒng)剛度系數(shù)和阻尼系數(shù)誤差Tab.2 Stiffness errors and damper errors of signal degree of freedom for free response

對圖4～圖6中的計算結(jié)果(r=2)統(tǒng)計見表3。

表3 圖4～圖6計算結(jié)果統(tǒng)計表Tab.3 Data from Fig.4 to Fig.6(r=2)

由圖4～圖6及表3可知：

(1)激勵力、位移和速度中任意兩項同時發(fā)生相位偏差會導致通過力狀態(tài)映射法辨識得到的系統(tǒng)剛度系數(shù)和阻尼系數(shù)出現(xiàn)一定誤差，誤差與激勵力、位移或者速度均為近似線性關系，且這種偏差對系統(tǒng)剛度系數(shù)的影響明顯小于對系統(tǒng)阻尼系數(shù)的影響；

(2)激勵力和速度相位偏差對剛度系數(shù)誤差的影響相近，均大于位移相位偏差對剛度系數(shù)誤差的影響；激勵力相位偏差對阻尼系數(shù)誤差的影響明顯大于其余兩者的影響；

圖4 激勵力和位移相位出現(xiàn)偏差時估計參數(shù)誤差Fig.4Parametererrorwhenphasedistortionhappenedonexcitationforceanddisplacement圖5 激勵力和速度相位出現(xiàn)偏差時估計參數(shù)誤差Fig.5Parametererrorwhenphasedistortionhappenedonexcitationforceandvelocity圖6 位移和速度相位出現(xiàn)偏差時估計參數(shù)誤差Fig.6Parametererrorwhenphasedistortionhappenedondisplacementandvelocity

(3)當激勵力頻率大于固有頻率時，剛度系數(shù)和阻尼系數(shù)誤差均隨激勵力頻率與固有頻率的比值的增加而增大；當激勵力頻率小于固有頻率時，剛度系數(shù)和阻尼系數(shù)誤差均隨激勵力頻率與固有頻率比值的較小而增大；當激勵力頻率等于固有頻率時，剛度系數(shù)和阻尼系數(shù)誤差最小。

4 結(jié) 論

本文構造了單自由度系統(tǒng)在強迫振動和自由振動兩種情況下的矢量圖模型，并推導得到了這兩種情況下系統(tǒng)某兩個檢測量出現(xiàn)相位偏差時系統(tǒng)辨識參數(shù)的誤差，當退化至某一檢測量發(fā)生相位偏差時誤差公式與文獻中對應結(jié)果一致。

研究結(jié)果表明：

(1)若干檢測量的相位偏差會導致系統(tǒng)辨識參數(shù)出現(xiàn)誤差，其中系統(tǒng)剛度系數(shù)受到的影響要明顯小于系統(tǒng)阻尼系數(shù)受到的影響；

(2)當激勵力頻率等于固有頻率時，辨識參數(shù)的誤差最小。

目前矢量圖法只適用于單自由度系統(tǒng)，對于多自由度系統(tǒng)，需要開展進一步研究。

[1] 吳志剛,王本利,馬興瑞.在軌航天器連接結(jié)構動力學及其參數(shù)辨識[J].宇航學報,1998,19(3):103-109. WU Zhigang, WANG Benli, MA Xingrui. Dynamic and parameters identification of joints in spacecraft [J]. Journal of Astronautics,1998,19(3):103-109.

[2] MARSI S F, CAUGHEY T K. Anon-parametric identification technique for nonlinear dynamic problems [J]. Journal of Applied Mechanics, 1979, 46: 433-447.

[3] MARSI S F, SASSI H, CAUGHEY T K. Nonparametric identification of nearly arbitrary nonlinear systems[J]. Journal of AppliedMechanics, 1982, 49:619-627.

[4] CRAWLEY E F, O’DONNELL K J. Force-state mapping identification of nonlinear joint[J]. AIAA Journal, 1987,25:1003-1010.

[5] CRAWLEY E F, AUBERT A C. Identification of nonlinear structural elements by force-statemapping [J]. AIAA Journal, 1986, 24:155-162.

[6] Al-HADID M A, WRIGHT J R. Developments in the force-state mapping technique for nonlinear systems and the extension to the location of nonlinear elements in a lumpe-parameter system [J]. Mechanical Systems and Signal Processing, 1989 ,3(3):269-290.

[7] 蔡力鋼,李玲,郭鐵能,等.基于力狀態(tài)映射法辨識非線性結(jié)合部動態(tài)參數(shù)[J].機械工程學報,2011,47(7):65-72. CAI Ligang, LI Ling, GUO Tieneng, et al. Identification of nonlinear joint parameters with force-state mapping method [J]. Journal of Mechanical Engineering, 2011, 47(7):65-72.

[8] 吳爽,趙壽根,吳大方,等.太陽翼鉸鏈結(jié)構的動力學試驗與非線性動力學建模[J].宇航學報,2013,34(12):1550-1556. WU Shuang, ZHAO Shougen, WU Dafang, et al. Dynamic experiments and nonlinear dynamics modeling of joints in solar array [J]. Journal of Astronautics, 2013,34(12):1550-1556.

[9] WRIGHT J R, Al-HADID M A. Sensitivity of the force-state mapping approach to measurement errors [J]. International Journal of Analytical and Experimental Modal Analysis, 1991, 6: 89-103.

[10] MESKELL C, FITZPATRICK J A. Errors in parameter estimates from the force-state mapping technique for free response due to phase distortion [J]. Journal of Sound and Vibration, 2002,252(5):967-974.

Error analysis in parameters from the force-state mapping method based on the vector diagram

JIANG Guoqing，LI Jiawen, TANG Guojin

(College of Aerospace Science and Engineering, National University of Defense Technology, Changsha 410073, China)

The force-state mapping is usually used to identify the equivalent dynamic parameters of a nonlinear joint, which needs lots of data from instruments or from differential or integral, such as excitation force, velocity, displacement, and acceleration. If those signals happen to suffer phase distortion, parameters identified from the force-state mapping will carry errors. Error formulations of identified parameters were derived by the vector diagram when the single degree of freedom system suffered a narrow range of phase distortion. The results indicate that the system stiffness is influenced heavily by the narrow range of phase distortion compared with the system damping.

vector diagram; force-state mapping; parameter identification; error analysis

國家自然科學基金(11502290)；工業(yè)裝備結(jié)構分析國家重點實驗室項目(GZ15104)

2015-06-29 修改稿收到日期：2015-11-20

蔣國慶 男，博士生，1987年8月生

李家文 男，博士，講師，1982年8月生

O327

A

10.13465/j.cnki.jvs.2016.24.021