摘" 要：靜力試驗(yàn)通過(guò)地面施加外載荷模擬飛機(jī)結(jié)構(gòu)真實(shí)受載情況，逐級(jí)加載并逐級(jí)測(cè)量結(jié)構(gòu)位移、應(yīng)變響應(yīng)。試驗(yàn)過(guò)程中數(shù)據(jù)的實(shí)時(shí)分析與預(yù)警，可以及時(shí)發(fā)現(xiàn)試驗(yàn)異常數(shù)據(jù)，監(jiān)控試驗(yàn)件狀態(tài)，提高靜力試驗(yàn)實(shí)施的安全性。該文研究相關(guān)性、斜率變化等試驗(yàn)數(shù)據(jù)分析方法，開(kāi)發(fā)靜力試驗(yàn)數(shù)據(jù)預(yù)警軟件，實(shí)現(xiàn)靜力試驗(yàn)中實(shí)時(shí)數(shù)據(jù)分析與預(yù)警，數(shù)據(jù)線(xiàn)性、重復(fù)性分析，正式試驗(yàn)快速結(jié)論等功能，提高靜力試驗(yàn)中響應(yīng)數(shù)據(jù)可視化、結(jié)果可視化能力。

關(guān)鍵詞：靜力試驗(yàn)；結(jié)構(gòu)響應(yīng)；數(shù)據(jù)分析；數(shù)據(jù)預(yù)警；軟件開(kāi)發(fā)；結(jié)果可視化

中圖分類(lèi)號(hào)：V216.1 " " 文獻(xiàn)標(biāo)志碼：A" " " " " 文章編號(hào)：2095-2945（2023）15-0061-04

Abstract： The static test simulates the real load condition of the aircraft structure by applying external load on the ground， and measures the structural displacement and strain response step by step. The real-time analysis and early warning of data in the test process can detect abnormal test data in time， monitor the state of test parts， and improve the safety of static test. The methods of test data analysis such as correlation analysis and slope change are studied， and the early warning software of static test data is developed， which realizes the functions of real-time data analysis and early warning， linear and repetitive analysis of data in static test， and rapid conclusion of formal test， thereby improves the ability of visualization of response data and results in static test.

Keywords： static test; structural response; data analysis; data early warning; software development; result visualization

結(jié)構(gòu)靜強(qiáng)度試驗(yàn)是考核飛機(jī)結(jié)構(gòu)強(qiáng)度的重要手段[1]。靜力試驗(yàn)中，在地面施加外載荷模擬飛機(jī)真實(shí)受載，并實(shí)時(shí)采集結(jié)構(gòu)的變形響應(yīng)[2-3]。通過(guò)試驗(yàn)中實(shí)測(cè)的位移、應(yīng)變數(shù)據(jù)獲取結(jié)構(gòu)真實(shí)變形情況，作為結(jié)構(gòu)響應(yīng)對(duì)比、有限元模型驗(yàn)證的基準(zhǔn)[4]。試驗(yàn)過(guò)程中對(duì)結(jié)構(gòu)響應(yīng)的實(shí)時(shí)分析與預(yù)警，可以及時(shí)發(fā)現(xiàn)高載時(shí)結(jié)構(gòu)變形異常部位，如結(jié)構(gòu)發(fā)生塑性變形或結(jié)構(gòu)出現(xiàn)局部失穩(wěn)?？煽康臄?shù)據(jù)采集與分析影響了試驗(yàn)的成敗，可以有效發(fā)現(xiàn)試驗(yàn)過(guò)程中的加載問(wèn)題或結(jié)構(gòu)問(wèn)題，提高試驗(yàn)實(shí)施的安全性和可靠性。

劉斯以等[5]提出了一種對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行最小二乘擬合，然后根據(jù)截距進(jìn)行零點(diǎn)偏移，消除試驗(yàn)件結(jié)構(gòu)、加載設(shè)備重量對(duì)測(cè)量結(jié)果影響的方法。王霞光等[6]提出了一種通過(guò)相關(guān)系數(shù)進(jìn)行數(shù)據(jù)篩選，然后選定線(xiàn)性區(qū)間利用最小二乘法進(jìn)行線(xiàn)性擬合的數(shù)據(jù)處理方法。柴東波等[7-8]提出了一種利用預(yù)試數(shù)據(jù)作為標(biāo)準(zhǔn)樣本構(gòu)建線(xiàn)性回歸方程和標(biāo)準(zhǔn)差，基于線(xiàn)性回歸方程和標(biāo)準(zhǔn)差構(gòu)建偏離度，快速識(shí)別非線(xiàn)性變化的試驗(yàn)數(shù)據(jù)的方法。

本文研究總結(jié)了靜力試驗(yàn)中數(shù)據(jù)分析及實(shí)時(shí)預(yù)警方法，在此基礎(chǔ)上，利用Python編程語(yǔ)言開(kāi)發(fā)了靜力試驗(yàn)數(shù)據(jù)預(yù)警軟件，實(shí)現(xiàn)了靜力試驗(yàn)實(shí)時(shí)數(shù)據(jù)分析與預(yù)警，數(shù)據(jù)線(xiàn)性、重復(fù)性分析，正式試驗(yàn)快速結(jié)論等功能，提高了靜力試驗(yàn)中響應(yīng)數(shù)據(jù)可視化、結(jié)果可視化能力。

1" 靜力試驗(yàn)數(shù)據(jù)分析方法

1.1" 線(xiàn)性分析

靜力試驗(yàn)中數(shù)據(jù)測(cè)量通常包含位移測(cè)量和應(yīng)變測(cè)量，按照試驗(yàn)大綱中規(guī)定的加載程序逐級(jí)加載，并同步采集應(yīng)變和位移數(shù)據(jù)。低載狀態(tài)下，結(jié)構(gòu)處于線(xiàn)彈性范圍，結(jié)構(gòu)的應(yīng)變、位移響應(yīng)與載荷呈線(xiàn)性關(guān)系，當(dāng)進(jìn)入到高載后，可能出現(xiàn)塑性變形或局部區(qū)域傳力路徑改變，此時(shí)部分區(qū)域響應(yīng)與載荷不再是完全線(xiàn)性。因此通過(guò)數(shù)據(jù)線(xiàn)性度變化識(shí)別結(jié)構(gòu)危險(xiǎn)部位，可以有效提高試驗(yàn)實(shí)施的安全性和可靠性。對(duì)響應(yīng)數(shù)據(jù)線(xiàn)性度判斷有以下幾種方法。

1.1.1" 線(xiàn)性相關(guān)系數(shù)計(jì)算

線(xiàn)性相關(guān)系數(shù)是研究變量之間線(xiàn)性相關(guān)程度的量[9]，通過(guò)計(jì)算載荷量級(jí)和各結(jié)構(gòu)響應(yīng)（位移、應(yīng)變）的線(xiàn)性相關(guān)系數(shù)，判斷結(jié)構(gòu)響應(yīng)是否與載荷線(xiàn)性相關(guān)。

2" 數(shù)據(jù)預(yù)警軟件開(kāi)發(fā)

靜力試驗(yàn)中，數(shù)據(jù)采集系統(tǒng)實(shí)時(shí)采集并轉(zhuǎn)發(fā)結(jié)構(gòu)響應(yīng)數(shù)據(jù)包。在試驗(yàn)過(guò)程中實(shí)時(shí)讀取數(shù)據(jù)包，并開(kāi)展數(shù)據(jù)分析，進(jìn)行可視化顯示，可以有效發(fā)現(xiàn)試驗(yàn)中存在的安全風(fēng)險(xiǎn)。試驗(yàn)完成后，快速對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行回收，并進(jìn)行數(shù)據(jù)分析，可快速得到試驗(yàn)結(jié)論，發(fā)現(xiàn)存在異常的數(shù)據(jù)?；赑ython編程語(yǔ)言，開(kāi)發(fā)了靜力試驗(yàn)數(shù)據(jù)預(yù)警軟件，主要包含數(shù)據(jù)監(jiān)控、數(shù)據(jù)預(yù)警、數(shù)據(jù)分析、可視化顯示等功能。軟件主界面如圖1所示。

2.1" 實(shí)時(shí)監(jiān)控與預(yù)警

靜力試驗(yàn)中隨著載荷逐級(jí)施加，實(shí)時(shí)同步采集結(jié)構(gòu)響應(yīng)數(shù)據(jù)，如從大量數(shù)據(jù)中篩選出可能存在異常的數(shù)據(jù)，可以有效降低試驗(yàn)風(fēng)險(xiǎn)。試驗(yàn)數(shù)據(jù)量巨大，連續(xù)加載過(guò)程中人工很難發(fā)現(xiàn)異常數(shù)據(jù)，因此需要通過(guò)實(shí)時(shí)分析，自動(dòng)識(shí)別異常數(shù)據(jù)并進(jìn)行預(yù)警。

軟件中采用基于斜率變化的方法判斷可能存在異常的數(shù)據(jù)。根據(jù)采集系統(tǒng)輸出數(shù)據(jù)編碼格式，將十六進(jìn)制數(shù)據(jù)轉(zhuǎn)化為易操作的十進(jìn)制數(shù)據(jù)。數(shù)據(jù)顯示包含2個(gè)窗口，一個(gè)是數(shù)據(jù)監(jiān)控窗口，如圖2所示，根據(jù)理論分析得到的危險(xiǎn)部位和重點(diǎn)監(jiān)測(cè)部位，選擇相應(yīng)的位移、應(yīng)變測(cè)點(diǎn)進(jìn)行實(shí)時(shí)跟蹤顯示。在數(shù)據(jù)監(jiān)控窗口中，還可以導(dǎo)入預(yù)試數(shù)據(jù)或理論數(shù)據(jù)，方便進(jìn)行數(shù)據(jù)對(duì)比。另一個(gè)是數(shù)據(jù)預(yù)警窗口，如圖3所示，在試驗(yàn)加載過(guò)程中實(shí)時(shí)計(jì)算所有響應(yīng)曲線(xiàn)的當(dāng)前斜率與前序斜率，如果斜率變化超過(guò)設(shè)定值，則在主界面中顯示出異常的應(yīng)變片號(hào)及斜率變化值，同時(shí)在斜率預(yù)警窗口中自動(dòng)顯示該響應(yīng)曲線(xiàn)，便于及時(shí)查找原因，排查試驗(yàn)安全風(fēng)險(xiǎn)。在2個(gè)觀(guān)察窗口中，可以捕捉測(cè)得點(diǎn)顯示測(cè)點(diǎn)信息，實(shí)時(shí)測(cè)量曲線(xiàn)上2點(diǎn)間斜率。

2.2" 預(yù)試數(shù)據(jù)分析

正式靜力試驗(yàn)前通常需要進(jìn)行至少2遍低載狀態(tài)的預(yù)試，通過(guò)預(yù)試一方面檢查試驗(yàn)加載設(shè)備、測(cè)量設(shè)備的可靠性，加載系統(tǒng)安裝的準(zhǔn)確性；另一方面與理論數(shù)據(jù)對(duì)比，預(yù)測(cè)正式試驗(yàn)結(jié)果。全面有效的預(yù)試數(shù)據(jù)分析關(guān)乎最終試驗(yàn)的成敗與安全。

2.2.1" 線(xiàn)性分析

預(yù)試完成后，對(duì)所有測(cè)量數(shù)據(jù)進(jìn)行線(xiàn)性分析。計(jì)算所有位移、應(yīng)變響應(yīng)與載荷之間的相關(guān)系數(shù)。然后對(duì)相關(guān)系數(shù)從小到大進(jìn)行分區(qū)，統(tǒng)計(jì)每個(gè)區(qū)間內(nèi)元素出現(xiàn)的個(gè)數(shù)或概率密度，由此得到所有響應(yīng)數(shù)據(jù)的線(xiàn)性相關(guān)系數(shù)分布（根據(jù)設(shè)置的截止值將小應(yīng)變和小位移剔除）。某試驗(yàn)60%限制載荷預(yù)試中（2遍）位移、應(yīng)變線(xiàn)性度及相關(guān)系數(shù)統(tǒng)計(jì)如圖4所示，線(xiàn)性相關(guān)系數(shù)設(shè)置了0.05的置信區(qū)間，在此區(qū)間內(nèi)的響應(yīng)數(shù)據(jù)認(rèn)為線(xiàn)性度良好。

2.2.2" 重復(fù)性分析

預(yù)試完成后，快速判斷多遍預(yù)試數(shù)據(jù)的重復(fù)性結(jié)果，重復(fù)性越高，說(shuō)明試驗(yàn)加載、測(cè)量越穩(wěn)定可靠。

選擇同一載荷量級(jí)（選擇預(yù)試最大載荷）進(jìn)行重復(fù)性比較，根據(jù)設(shè)置的位移和應(yīng)變截止值，分別計(jì)算相對(duì)誤差和絕對(duì)誤差，并分類(lèi)進(jìn)行統(tǒng)計(jì)。

2.2.3" 異常數(shù)據(jù)輸出

對(duì)線(xiàn)性相關(guān)系數(shù)低于0.95的數(shù)據(jù)輸出，對(duì)重復(fù)性大于5或5%的數(shù)據(jù)輸出，對(duì)響應(yīng)數(shù)據(jù)一直為零的數(shù)據(jù)判斷為壞片并輸出，輸出所有異常數(shù)據(jù)，方便進(jìn)一步排查原因。

2.3" 正式試驗(yàn)快速結(jié)論

正式試驗(yàn)結(jié)束后，快速得到試驗(yàn)結(jié)論，包含了線(xiàn)性分析與統(tǒng)計(jì)、殘余變形統(tǒng)計(jì)，如圖5所示。

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

研究了靜力試驗(yàn)中測(cè)量得到的響應(yīng)數(shù)據(jù)分析方法及試驗(yàn)過(guò)程中數(shù)據(jù)實(shí)時(shí)預(yù)警方法，基于Python編程語(yǔ)言開(kāi)發(fā)了靜力試驗(yàn)數(shù)據(jù)預(yù)警軟件，實(shí)現(xiàn)了靜力試驗(yàn)中實(shí)時(shí)數(shù)據(jù)分析與預(yù)警，數(shù)據(jù)線(xiàn)性、重復(fù)性分析，正式試驗(yàn)快速結(jié)論等功能，對(duì)數(shù)據(jù)進(jìn)行高效分析并分類(lèi)統(tǒng)計(jì)，結(jié)果可視化輸出，使試驗(yàn)實(shí)施更加安全高效，并在工程中得到應(yīng)用推廣。

參考文獻(xiàn)：

[1] 國(guó)防科學(xué)技術(shù)工業(yè)委員會(huì).飛機(jī)結(jié)構(gòu)靜強(qiáng)度試驗(yàn)通用要求：HB 7713—2002[S].

[2] 鄭建軍，唐吉運(yùn)，王彬文.C919飛機(jī)全機(jī)靜力試驗(yàn)技術(shù)[J].航空學(xué)報(bào)，2019，40（1）：210-221.

[3] 強(qiáng)寶平.全尺寸飛機(jī)結(jié)構(gòu)試驗(yàn)技術(shù)[J].航空科學(xué)技術(shù)，2012（6）：10-13.

[4] 趙峻峰，李三平，李強(qiáng).民用飛機(jī)機(jī)體結(jié)構(gòu)靜強(qiáng)度驗(yàn)證[J].民用飛機(jī)設(shè)計(jì)與研究，2020（2）：1-5.

[5] 劉斯以，王崇，張俊愿.懸吊式靜力試驗(yàn)應(yīng)變測(cè)量數(shù)據(jù)快速處理方法[J].直升機(jī)技術(shù)，2018（2）：55-59.

[6] 王霞光，左佳.飛機(jī)結(jié)構(gòu)靜力試驗(yàn)中應(yīng)變測(cè)量數(shù)據(jù)線(xiàn)性化處理方法研究[C]//第六屆中國(guó)航空學(xué)會(huì)青年科技論壇文集（上冊(cè)），2014：461-464.

[7] 柴東波，郭曉冬，楊東濤.全機(jī)靜力試驗(yàn)預(yù)警方法的研究[J].工程與試驗(yàn)，2020，60（4）：13-15.

[8] 柴東波，張建鋒，郭曉冬，等.飛機(jī)結(jié)構(gòu)靜強(qiáng)度試驗(yàn)實(shí)時(shí)預(yù)警方法研究[J].工程與試驗(yàn)，2022，62（1）：36-37，53.

[9] 孟慶寬，劉剛，張漫，等.基于線(xiàn)性相關(guān)系數(shù)約束的作物行中心線(xiàn)檢測(cè)方法[J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2013，44（S1）：216-223.

[10] 田垅，劉宗田.最小二乘法分段直線(xiàn)擬合[J].計(jì)算機(jī)科學(xué)，2012，39（S1）：482-484.