某機(jī)載雷達(dá)框架結(jié)構(gòu)動(dòng)力學(xué)分析

盛文軍，毛 亮，王志海

(中國(guó)電子科技集團(tuán)公司 第三十八研究所，安徽 合肥 230088)

0 引言

機(jī)載雷達(dá)是整個(gè)飛機(jī)系統(tǒng)重要的組成部分，起著測(cè)繪、成像等重要作用[1]，在滑翔、加速和降落等飛行階段需要承受隨機(jī)振動(dòng)、沖擊等各種嚴(yán)酷載荷[2-3]。如何在滿足雷達(dá)重量、精度等系統(tǒng)指標(biāo)要求的同時(shí)，保證剛強(qiáng)度指標(biāo)的要求是雷達(dá)結(jié)構(gòu)設(shè)計(jì)中的重點(diǎn)。

機(jī)載雷達(dá)結(jié)構(gòu)的傳統(tǒng)動(dòng)力學(xué)特性驗(yàn)證方法基于振動(dòng)試驗(yàn)[4-5]，這種設(shè)計(jì)-試驗(yàn)-改進(jìn)的方法成本高、周期長(zhǎng)，機(jī)載設(shè)備的大試驗(yàn)周期、高振動(dòng)量級(jí)需求不能滿足當(dāng)前產(chǎn)品研發(fā)過程中降低成本、快速迭代的要求。隨著仿真技術(shù)的快速發(fā)展和計(jì)算機(jī)計(jì)算能力的快速提升，結(jié)構(gòu)動(dòng)力學(xué)仿真計(jì)算的準(zhǔn)確度和便捷性大大提高?；谟?jì)算機(jī)仿真的虛擬試驗(yàn)技術(shù)能夠?qū)崿F(xiàn)設(shè)計(jì)方案的快速驗(yàn)證，促進(jìn)產(chǎn)品研發(fā)的快速迭代，是機(jī)載設(shè)備動(dòng)力學(xué)特性驗(yàn)證方法的發(fā)展趨勢(shì)。本文基于有限元方法對(duì)某機(jī)載雷達(dá)框架結(jié)構(gòu)進(jìn)行了建模，并基于模態(tài)疊加法分析了結(jié)構(gòu)的耐振動(dòng)和耐沖擊動(dòng)力學(xué)特性，得到了結(jié)構(gòu)在載荷作用下的應(yīng)力、應(yīng)變?cè)茍D，并對(duì)結(jié)構(gòu)的剛強(qiáng)度進(jìn)行了校核。

1 機(jī)載環(huán)境條件

機(jī)載條件下雷達(dá)也承受沖擊載荷，沖擊采用半正弦波，峰值加速度為15g，沖擊脈寬為11 ms，沖擊沿坐標(biāo)軸的3個(gè)方向進(jìn)行。

圖1 隨機(jī)振動(dòng)功率譜密度

2 機(jī)載雷達(dá)結(jié)構(gòu)及有限元模型

機(jī)載雷達(dá)框架是雷達(dá)的主要承力結(jié)構(gòu)，其在振動(dòng)沖擊環(huán)境下的剛強(qiáng)度直接影響了雷達(dá)的輸出功率和波形等電訊指標(biāo)。機(jī)載雷達(dá)框架結(jié)構(gòu)模型如圖2所示，主要負(fù)重分布在中間橫隔板上，長(zhǎng)度方向分布有8個(gè)支耳，通過螺釘固定在平臺(tái)上。框架材料為5A06鋁，其密度為2 700 kg/m3，彈性模量為71 GPa，泊松比為0.33，屈服強(qiáng)度為155 MPa。

圖2 機(jī)載雷達(dá)框架結(jié)構(gòu)模型

考慮到結(jié)構(gòu)實(shí)際模型比較復(fù)雜，直接進(jìn)行有限元建模計(jì)算規(guī)模較大，為了提高分析效率，在進(jìn)行有限元分析之前，對(duì)模型進(jìn)行簡(jiǎn)化，去除對(duì)結(jié)構(gòu)剛強(qiáng)度影響不大的倒角、鉆孔等幾何特征，減小分析規(guī)模。采用實(shí)體單元六面體網(wǎng)格CTRIA6對(duì)結(jié)構(gòu)進(jìn)行網(wǎng)格劃分，框架與平臺(tái)的固定和連接用剛性連接模擬，建立的框架有限元模型如3圖所示，共64 884個(gè)單元，90 286個(gè)節(jié)點(diǎn)。

圖3 框架有限元模型

3 模態(tài)分析

結(jié)構(gòu)動(dòng)力學(xué)分析基于線性模態(tài)疊加，模態(tài)分析是動(dòng)力學(xué)分析的基礎(chǔ)[6]。模態(tài)分析將振動(dòng)動(dòng)力學(xué)微分方程進(jìn)行解耦，成為一組以模態(tài)坐標(biāo)及模態(tài)參數(shù)為描述的獨(dú)立方程。

對(duì)于多自由度系統(tǒng)，如果考慮黏性阻尼，則其受迫振動(dòng)的微分方程為：

(1)

由于阻尼對(duì)系統(tǒng)固有頻率的影響較小，在計(jì)算系統(tǒng)的固有頻率時(shí)，暫不考慮阻尼帶來的影響?？紤]自由振動(dòng)的狀態(tài)，系統(tǒng)的受迫振動(dòng)微分方程變?yōu)椋?/p>

(2)

求解式(2)，得到系統(tǒng)的特征值方程為：

其中：ωi為系統(tǒng)的第i階固有圓頻率；φi為系統(tǒng)的第i階振型。

至此，系統(tǒng)動(dòng)力學(xué)微分方程完成了解耦，系統(tǒng)對(duì)正弦激勵(lì)的響應(yīng)可以等效于各階振型的疊加，考慮到任意激勵(lì)都可以等效為不同頻率正弦激勵(lì)的疊加，從而實(shí)現(xiàn)系統(tǒng)對(duì)任意激勵(lì)動(dòng)力學(xué)響應(yīng)的計(jì)算。本文計(jì)算中為了減少計(jì)算量，采用振型截?cái)喾?，只?jì)算有效質(zhì)量總占比90%的前r階振型的響應(yīng)。

本文中有限元模型的建立基于很多等效和簡(jiǎn)化，為了得到更準(zhǔn)確的分析結(jié)果，如圖4(a)所示，對(duì)鋁框架進(jìn)行了500 Hz～1 000 Hz的正弦掃頻試驗(yàn)，設(shè)置控制點(diǎn)和監(jiān)測(cè)點(diǎn)兩個(gè)加速度傳感器，試驗(yàn)結(jié)果如圖4(b)所示。監(jiān)測(cè)點(diǎn)在873.2 Hz出現(xiàn)峰值，則結(jié)構(gòu)的基頻為873.2 Hz，基于此，對(duì)有限元模型進(jìn)行了進(jìn)一步的修正和確認(rèn)，從而保證模型的準(zhǔn)確性[7]。

圖4正弦掃頻試驗(yàn)及結(jié)果

對(duì)框架結(jié)構(gòu)的有限元模型進(jìn)行模態(tài)分析，得到的前4階模態(tài)如圖5所示。由圖5可以看出:1階模態(tài)為橫支撐板上下振動(dòng)，2階模態(tài)為橫支撐板左右部分分別上下振動(dòng)，3、4階模態(tài)為側(cè)面擋板前后振動(dòng)；基頻為837.1 Hz，與試驗(yàn)值873.2 Hz相當(dāng)，確認(rèn)了有限元模型的準(zhǔn)確性。

圖5 框架結(jié)構(gòu)的前4階模態(tài)

圖6為框架結(jié)構(gòu)在X、Y、Z三個(gè)方向的前30階模態(tài)的有效質(zhì)量。動(dòng)力學(xué)分析時(shí)截取有效質(zhì)量占比前90%的模態(tài)，以保證計(jì)算的準(zhǔn)確性。

圖6 X、Y、Z三向模態(tài)有效質(zhì)量分布

4 隨機(jī)振動(dòng)

機(jī)載環(huán)境下工作時(shí)，雷達(dá)會(huì)受到自身或外部的隨機(jī)振動(dòng)載荷的影響，導(dǎo)致其產(chǎn)生內(nèi)部應(yīng)力和變形，應(yīng)力應(yīng)變的幅度應(yīng)當(dāng)滿足設(shè)計(jì)指標(biāo)，從而保證雷達(dá)的剛強(qiáng)度要求。通常情況下，認(rèn)為隨機(jī)振動(dòng)是穩(wěn)態(tài)的隨機(jī)過程，符合正態(tài)分布的統(tǒng)計(jì)規(guī)律，振動(dòng)過程的激勵(lì)和響應(yīng)的統(tǒng)計(jì)特性如振幅的均值和方差都是穩(wěn)定的，功率譜密度給出的振動(dòng)譜值為1σ值(σ為正態(tài)分布的均方根值，幅值小于1σ的概率為68.3%，幅值小于3σ的概率為99.7%)，相應(yīng)的剛強(qiáng)度計(jì)算結(jié)果也是1σ值，為了得到較高的結(jié)果置信度，剛強(qiáng)度評(píng)價(jià)時(shí)取3σ值。本文基于模態(tài)疊加法進(jìn)行計(jì)算，取模態(tài)阻尼為0.03，截取有效質(zhì)量占比前90%的模態(tài)計(jì)算框架結(jié)構(gòu)對(duì)隨機(jī)振動(dòng)三個(gè)方向激勵(lì)的響應(yīng)進(jìn)行計(jì)算。框架結(jié)構(gòu)隨機(jī)振動(dòng)時(shí)1σ的應(yīng)力、變形云圖如圖7所示。由圖7可知：Z方向振動(dòng)時(shí)，應(yīng)力值最大，3σ值為106.2 MPa，位于橫隔板支撐點(diǎn)，取安全系數(shù)為1.3，5A06的屈服強(qiáng)度為155 MPa，許用應(yīng)力為119 MPa，大于材料最大3σ應(yīng)力值106.2 MPa，結(jié)構(gòu)滿足強(qiáng)度設(shè)計(jì)要求；Z方向振動(dòng)時(shí)變形最大，3σ變形值為0.25 mm，位于橫隔板中間，小于設(shè)計(jì)要求的結(jié)構(gòu)最大變形0.5 mm，滿足剛強(qiáng)度設(shè)計(jì)要求。

圖7 框架結(jié)構(gòu)隨機(jī)振動(dòng)時(shí)1σ的應(yīng)力、變形云圖

5 沖擊振動(dòng)

沖擊是機(jī)載環(huán)境下常見的載荷之一，考察系統(tǒng)受到短暫的脈沖載荷激勵(lì)下的響應(yīng)。本文對(duì)框架結(jié)構(gòu)進(jìn)行沖擊動(dòng)力學(xué)分析，結(jié)果如圖8所示。對(duì)比三向振動(dòng)結(jié)構(gòu)應(yīng)力值，Z方向沖擊時(shí)應(yīng)力最大為53.0 MPa，位于橫隔板支撐點(diǎn)，考慮結(jié)構(gòu)要求最大許用應(yīng)力為119 MPa。對(duì)比三向振動(dòng)結(jié)構(gòu)變形值，Z方向沖擊時(shí)最大變形為0.12 mm，位于橫隔板中間，考慮設(shè)計(jì)要求最大變形為0.5 mm，滿足剛強(qiáng)度要求。

圖8 框架結(jié)構(gòu)沖擊振動(dòng)應(yīng)力、變形云圖

6 結(jié)論

本文針對(duì)機(jī)載雷達(dá)框架結(jié)構(gòu)進(jìn)行了動(dòng)力學(xué)分析，分析結(jié)果表明隨機(jī)振動(dòng)時(shí)，最大3σ應(yīng)力106.2 MPa，最大3σ變形0.25 mm；沖擊振動(dòng)時(shí)，最大應(yīng)力53.0 MPa，最大變形0.12 mm，均滿足結(jié)構(gòu)設(shè)計(jì)剛強(qiáng)度要求。本文的分析結(jié)果能夠有力地輔助結(jié)構(gòu)設(shè)計(jì)，減少環(huán)境試驗(yàn)次數(shù)，降低開發(fā)成本，是機(jī)載雷達(dá)研發(fā)過程中必不可少的技術(shù)手段。