徐承軍 羅 東 常亞瓊

（武漢理工大學(xué)物流工程學(xué)院 武漢 430063）

0 引 言

自振弦式傳感器發(fā)明以來，一直受到工程界的重視，由于其自身特殊的工作機(jī)理和優(yōu)越性，目前已經(jīng)在多個(gè)領(lǐng)域得到廣泛應(yīng)用．振弦傳感器是基于鋼弦頻率隨張力的變化而變化的原理而工作的，與工程、科研中普遍應(yīng)用的電阻應(yīng)變計(jì)等其他傳感器相比，其優(yōu)越性主要表現(xiàn)在［1］：（l）振弦傳感器以振弦頻率的變化量來表征受力的大小，因而具有良好的抗干擾能力和遠(yuǎn)距離傳輸能力．相對(duì)于電阻應(yīng)變計(jì)等傳感器的模擬量輸出來說，能更為簡(jiǎn)單方便地進(jìn)行數(shù)據(jù)采集、傳輸、處理和存儲(chǔ)，實(shí)現(xiàn)高精度的自動(dòng)測(cè)試．（2）從目前工程現(xiàn)場(chǎng)的測(cè)試要求和惡劣環(huán)境來看，由于振弦傳感器自身獨(dú)特的結(jié)構(gòu)形式、工作方式以及安裝方式等原因，相對(duì)于其他傳感器來說，具有優(yōu)越的適應(yīng)環(huán)境的能力和長(zhǎng)期穩(wěn)定性．有資料表明，目前振弦式傳感器的最長(zhǎng)連續(xù)工作時(shí)間已經(jīng)達(dá)到了27 a．

此外，振弦式傳感器還具有受電參數(shù)影響小、零點(diǎn)飄移小、受溫度影響小、性能穩(wěn)定可靠、耐震動(dòng)等優(yōu)點(diǎn)．目前振弦式傳感器主要應(yīng)用于巖土工程進(jìn)行靜應(yīng)力監(jiān)測(cè)，并且取得了很好的結(jié)果．大量的國(guó)內(nèi)外文獻(xiàn)查閱表明，尚未發(fā)現(xiàn)其在結(jié)構(gòu)應(yīng)力動(dòng)態(tài)測(cè)試中的應(yīng)用．

本文對(duì)振弦式傳感器的工作原理進(jìn)行了研究，設(shè)計(jì)了基于振弦式傳感器的結(jié)構(gòu)應(yīng)力動(dòng)態(tài)信號(hào)采集系統(tǒng)，并利用試制樣機(jī)與經(jīng)計(jì)量的信號(hào)測(cè)試分析系統(tǒng)進(jìn)行了對(duì)比分析驗(yàn)證．實(shí)驗(yàn)分析結(jié)果初步證實(shí)了振弦式傳感器應(yīng)用于結(jié)構(gòu)應(yīng)力動(dòng)態(tài)測(cè)試中的可行性，為其在結(jié)構(gòu)應(yīng)力動(dòng)態(tài)監(jiān)測(cè)的工程應(yīng)用打下了基礎(chǔ)．

1 振弦式傳感器動(dòng)態(tài)數(shù)據(jù)采集系統(tǒng)

1．1 硬件模塊

根據(jù)振弦式傳感器的工作原理，采用目前較為成熟的單片機(jī)89C51芯片［2］，利用該芯片的I／O口實(shí)現(xiàn)對(duì)附加電路的控制和信號(hào)的采集等功能．以該芯片內(nèi)的計(jì)數(shù)器作為計(jì)數(shù)單元，對(duì)鋼弦傳感器頻率信號(hào)計(jì)數(shù)，并進(jìn)行分析處理．利用單片機(jī)的I／O口對(duì)傳感器觸發(fā)電路進(jìn)行控制，發(fā)高壓脈沖觸發(fā)傳感器，采集傳感器信號(hào)．由于試驗(yàn)樣機(jī)利用的單線圈振弦式傳感器，在單片機(jī)的外圍電路中采用光耦電路來實(shí)現(xiàn)對(duì)激勵(lì)環(huán)節(jié)和采集環(huán)節(jié)的切換．

對(duì)于頻率測(cè)量模塊來說，傳感器每2 ms輸出一次測(cè)量結(jié)果，因而動(dòng)態(tài)數(shù)據(jù)采集數(shù)據(jù)量很大，不適合使用ASCII碼通信，因而在此改用了二進(jìn)制碼以減小數(shù)據(jù)量．系統(tǒng)設(shè)計(jì)的通信數(shù)據(jù)格式見表1．

表1 通信數(shù)據(jù)格式

系統(tǒng)的通信方式則使用目前比較成熟的串口通信，設(shè)置波特率為115 200 bit／s，這樣也便于和其他控制通信設(shè)備進(jìn)行對(duì)接．

1．2 軟件模塊

上位機(jī)軟件和系統(tǒng)界面選用目前做常用的軟件編程工具VC6．0來完成．其架構(gòu)見圖1．

圖1 軟件架構(gòu)模塊圖

軟件模塊由通信模塊、數(shù)據(jù)顯示模塊、數(shù)據(jù)處理模塊和數(shù)據(jù)存儲(chǔ)模塊組成．具體工作過程為：上位機(jī)發(fā)送接收從機(jī)的地址，等待接收從機(jī)的應(yīng)答信號(hào)，然后當(dāng)下位機(jī)接到開始信號(hào)后，開始進(jìn)行數(shù)據(jù)采集并通過串口向上位機(jī)傳輸數(shù)據(jù)．上位機(jī)收到數(shù)據(jù)后，立即處理數(shù)據(jù)并記錄目前系統(tǒng)時(shí)間，同時(shí)做出應(yīng)力信號(hào)的時(shí)間歷程曲線圖，直到數(shù)據(jù)采集結(jié)束后，按停止按鈕結(jié)束采集．

2 對(duì)比驗(yàn)證

為了驗(yàn)證振弦式傳感器動(dòng)態(tài)測(cè)試系統(tǒng)數(shù)據(jù)的準(zhǔn)確性，本文使用經(jīng)專業(yè)校驗(yàn)單位進(jìn)行過計(jì)量標(biāo)準(zhǔn)認(rèn)定的“東華DHDA（5923＿1394）信號(hào)測(cè)試分析系統(tǒng)”（以下簡(jiǎn)稱“比對(duì)系統(tǒng)”）來對(duì)其測(cè)試數(shù)據(jù)的準(zhǔn)確性進(jìn)行對(duì)比．東華DHDA信號(hào)測(cè)試分析系統(tǒng)可使用金屬應(yīng)變片來進(jìn)行金屬結(jié)構(gòu)的應(yīng)力測(cè)試，具有較長(zhǎng)的工業(yè)應(yīng)用背景．因此，使用該測(cè)試分析系統(tǒng)對(duì)振弦式傳感器動(dòng)態(tài)測(cè)試數(shù)據(jù)進(jìn)行驗(yàn)證分析，具有代表性和科學(xué)性．

2．1 試驗(yàn)方法與步驟

試驗(yàn)采用了懸臂梁來模擬機(jī)械的金屬結(jié)構(gòu)，傳感器的布置見圖2．

圖2 傳感器安裝示意圖

在懸臂梁根部同一截面的上、下表面分別安裝了振弦式傳感器與金屬應(yīng)變片，并連接振弦式傳感器動(dòng)態(tài)采集系統(tǒng)試驗(yàn)樣機(jī)和比對(duì)系統(tǒng)．待兩系統(tǒng)均調(diào)試準(zhǔn)備好后，同時(shí)開始進(jìn)行數(shù)據(jù)采集．由于兩個(gè)測(cè)試系統(tǒng)分別位于同一截面的上下兩個(gè)表面，從理想狀況上來說，其采集的數(shù)據(jù)應(yīng)該是大小相同，方向相反的，在應(yīng)力數(shù)據(jù)曲線上應(yīng)顯示為衰減幅度和速度相同，而且幅值上下對(duì)稱的波形．

2．2 實(shí)驗(yàn)數(shù)據(jù)

如圖3是比對(duì)系統(tǒng)的數(shù)據(jù)采集曲線，圖4為振弦式傳感器數(shù)據(jù)采集信號(hào)系統(tǒng)的數(shù)據(jù)采集曲線，橫軸均為時(shí)間，縱軸為結(jié)構(gòu)應(yīng)變．以圖3為例，試驗(yàn)開始較平穩(wěn)的階段為無載荷狀態(tài)，曲線中A點(diǎn)時(shí)刻開始施加載荷，A點(diǎn)到B點(diǎn)之間通過手工隨意施加緩慢交變載荷，B點(diǎn)時(shí)刻在懸臂端給予初始位移后突然釋放，使得懸臂梁進(jìn)行自由阻尼振蕩，由于阻尼使懸臂梁振幅逐漸減小，C點(diǎn)時(shí)刻重復(fù)施加緩慢交變載荷，然后再次自由振蕩，依次循環(huán)若干次．

圖3 比對(duì)系統(tǒng)時(shí)間應(yīng)變圖

圖4 振弦式傳感器數(shù)據(jù)采集信號(hào)系統(tǒng)時(shí)間應(yīng)變圖

2．3 振弦式傳感器采集系統(tǒng)準(zhǔn)確性分析

從肉眼直觀看來，上述2套系統(tǒng)采集的數(shù)據(jù)曲線具有很高的相似性，將其中一組數(shù)據(jù)進(jìn)行反相，并將兩曲線的橫軸變成相同的時(shí)間刻度，放到一個(gè)坐標(biāo)系下繪制出來，兩者的吻合性也明顯．為了從科學(xué)的角度驗(yàn)證兩者的相似性，本文采用Matlab的互相關(guān)分析方法來進(jìn)行驗(yàn)證［3］，從而確定振弦式傳感器采集系統(tǒng)的準(zhǔn)確性．

由于2個(gè)系統(tǒng)的采樣方式不同，因此兩組數(shù)據(jù)的采集點(diǎn)不是一一對(duì)應(yīng)的，所以在進(jìn)行相關(guān)性分析之前，首先要進(jìn)行數(shù)據(jù)插值，使兩數(shù)據(jù)具有相同的步長(zhǎng)和取值范圍，然后才能夠進(jìn)行相關(guān)性分析．

2．3．1 數(shù)據(jù)插值處理

數(shù)據(jù)插值是在離散數(shù)據(jù)的基礎(chǔ)上補(bǔ)插連續(xù)函數(shù)，構(gòu)成連續(xù)曲線，而這條連續(xù)曲線就可以通過全部要求的離散數(shù)據(jù)點(diǎn)［4］．它的原理是根據(jù)在有限離散點(diǎn)處的取值情況，估算出函數(shù)在其他點(diǎn)處的近似值．Matlab自帶數(shù)據(jù)插值處理的函數(shù)，進(jìn)行數(shù)據(jù)插值計(jì)算的Matlab命令如下．

該命令用多項(xiàng)式技術(shù)計(jì)算實(shí)現(xiàn)一維插值．

式中：Y為函數(shù)值矢量；X為自變量取值范圍；xi為插值點(diǎn)的自變量矢量；method為插值方法選項(xiàng)．

Matlab實(shí)現(xiàn)一維插值共有4種方法：

臨近點(diǎn)插值：method＝＇nearest＇

線性插值：method＝＇linear＇

三次樣條插值：method＝＇spline＇

立方插值：method＝＇pchip＇or＇cubic＇

對(duì)于［min｛xi｝，max｛xi｝］外的值，Matlab使用外推的方法計(jì)算數(shù)值．

本文試驗(yàn)中2種系統(tǒng)采集到的數(shù)據(jù)都是離散值，用Matlab完成兩者數(shù)據(jù)插值的代碼如下：

A＝load（＇D：＼AA．txt＇）：%A 是振弦傳感器系統(tǒng)測(cè)得的值B＝load（＇D：＼BB．txt＇）：%B 是東華測(cè)試儀測(cè)得的值

ax＝ A（：，2）：

ay＝A（：，3）：

bx＝（0．000：0．005：162．79）＋0．15：

by＝－B（：，2）＊0．21：

plot（ax，ay，bx，by）：

xlabel（＇時(shí)間（s）＇），

ylabel（＇應(yīng)力（MP）＇）

： title（＇振弦式傳感器數(shù)據(jù)跟東華金屬應(yīng)變片數(shù)據(jù)原始數(shù)據(jù)＇）：

legend（＇振弦式應(yīng)變計(jì)數(shù)據(jù)＇，＇金屬應(yīng)變片數(shù)據(jù)＇）：

ha＝1：0．001：19：

hb＝ha：

ta＝interp1（ax，ay，ha，＇linear＇）：

figure（2）：

subplot（3，1，1）：

plot（ha，ta）：

xlabel（＇時(shí)間（s）＇），ylabel（＇）：

title（＇插值后的振弦式傳感器數(shù)據(jù)＇）：

tb＝interp1（bx，by，hb，＇linear＇）：subplot（3，1，2）：

plot（hb，tb，＇g＇）：

xlabel（＇時(shí)間（s）＇），

ylabel（＇應(yīng)力（MP）＇）

：title（＇插值后的金屬應(yīng)變片數(shù)據(jù)＇）：

subplot（3，1，3）：

plot（ha，ta，hb，tb）：

xlabel（＇時(shí)間（s）＇），ylabel（＇應(yīng)力（MP）＇）：

title（＇插值后振弦式傳感器數(shù)據(jù)跟東華金屬應(yīng)變片數(shù)據(jù)＇）：

legend（＇振弦式應(yīng)變計(jì)數(shù)據(jù)＇，＇金屬應(yīng)變片數(shù)據(jù)＇）：

圖5a）是應(yīng)用Matlab繪制的振弦式傳感器采集系統(tǒng)與比對(duì)系統(tǒng)獲得的原始數(shù)據(jù)曲線．為了便于進(jìn)行相關(guān)性分析，對(duì)利用振弦式傳感器獲得的數(shù)據(jù)進(jìn)行了幅值反向．圖5b）從上至下分別是用Matlab插值后繪制的振弦式傳感器數(shù)據(jù)曲線、比對(duì)系統(tǒng)數(shù)據(jù)曲線以及兩者在同一坐標(biāo)下的曲線合并．

圖5 兩采集系統(tǒng)的數(shù)據(jù)曲線

2．3．2 數(shù)據(jù)相關(guān)性分析

在信號(hào)分析里，互相關(guān)分析主要分析的是兩個(gè)時(shí)間序列之間的相關(guān)程度．在進(jìn)行互相關(guān)分析時(shí)，兩種數(shù)據(jù)的相關(guān)程度可以用他們的相關(guān)系數(shù)來表示．相關(guān)函數(shù)的相關(guān)系數(shù)只是一個(gè)比率，無單位量綱，一般取小數(shù)點(diǎn)后兩位來表示．相關(guān)系數(shù)的絕對(duì)值表示相關(guān)的程度，正負(fù)號(hào)表示相關(guān)的方向［5］．相關(guān)系數(shù)跟相關(guān)程度的關(guān)系見表2．

表2 相關(guān)系數(shù)與相關(guān)程度對(duì)照關(guān)系表

下面是進(jìn)行互相關(guān)性分析的Matlab代碼：

［C，lags］＝ xcorr（ta，tb，＇coeff＇）；figure（3）；

plot（lags，C）；

xlabel（＇偏移個(gè)數(shù)＇），ylabel（＇互相關(guān)系數(shù)＇）；%x軸y 軸標(biāo)簽

title（＇A 和B 的互相關(guān)圖形＇，＇Color＇，＇r＇）；

2．4 結(jié)論

根據(jù)圖6所示結(jié)果，振弦式傳感器采集系統(tǒng)獲得的曲線與比對(duì)系統(tǒng)采集的數(shù)據(jù)曲線的相關(guān)系數(shù)達(dá)到了0．952 7，超過了0．95，屬于高度相關(guān)，即振弦式傳感器所采集到的數(shù)據(jù)與比對(duì)系統(tǒng)所采集的數(shù)據(jù)非常相似．

圖6 2組數(shù)據(jù)的互相關(guān)曲線

3 結(jié)束語(yǔ)

本文通過Matlab對(duì)測(cè)試數(shù)據(jù)的相關(guān)性分析，表明所研制的振弦式傳感器數(shù)據(jù)采集系統(tǒng)可行，如實(shí)反映了結(jié)構(gòu)的真實(shí)應(yīng)力．此外據(jù)資料統(tǒng)計(jì)，港口機(jī)械鋼結(jié)構(gòu)的振動(dòng)頻率一般在10 Hz左右［6］，而本文設(shè)計(jì)的振弦式傳感器數(shù)據(jù)采集系統(tǒng)的采樣頻率可達(dá)到200 Hz以上，因此，本振弦式傳感器動(dòng)態(tài)數(shù)據(jù)采集系統(tǒng)可以滿足港口機(jī)械結(jié)構(gòu)的動(dòng)態(tài)應(yīng)力檢測(cè)的要求．

［1］鄧鐵六，趙勝利，趙振遠(yuǎn)，等．高準(zhǔn)確度振弦式傳感器研究［J］．山東科技大學(xué)學(xué)報(bào)：自然科學(xué)版，2010，29（1）：52－57

［2］遲欽河，趙仲生，喬桂芳，等．89C51單片機(jī)在多通道數(shù)據(jù)采集系統(tǒng)中的應(yīng)用［J］．自動(dòng)化儀表，2000，21（6）：33－39．

［3］劉文白，徐海俠．砂土宏觀力學(xué)特性與細(xì)觀結(jié)構(gòu)的相關(guān)性試驗(yàn)研究［J］．武漢理工大學(xué)學(xué)報(bào)：交通科學(xué)與工程版，2011，35（4）：683－686．

［4］陳海龍，李 宏．基于 MATLAB的偽隨機(jī)序列的產(chǎn)生和分析［J］．計(jì)算機(jī)仿真，2005（5）：98－100．

［5］莊楚強(qiáng) 吳亞森，應(yīng)用數(shù)理統(tǒng)計(jì)基礎(chǔ)［M］．3版．廣州：華南理工大學(xué)出版社，2006．

［6］徐承軍．基于無線局域網(wǎng)的集裝箱碼頭機(jī)械調(diào)度系統(tǒng)的仿真、優(yōu)化與監(jiān)控［D］．武漢：武漢理工大學(xué)物流工程學(xué)院，2007．