謝 歡, 王紅巖, 郝丙飛, 李善樂, 芮 強
(陸軍裝甲兵學院 車輛工程系,北京 100072)
坦克炮在射擊時,受到外界激勵的作用,火炮身管將發(fā)生彎曲變形,而炮口的振動變形將對坦克炮的射擊精度產(chǎn)生很大影響[1].模態(tài)分析是一種研究物體固有屬性的方法,通過模態(tài)分析可以得到物體的模態(tài)參數(shù),為進一步研究火炮身管的振動特性提供基礎.
針對坦克火炮身管的振動特性,李猛等[2]研究了火炮身管振動性能與射擊頻率的匹配關系,通過對身管模型進行模態(tài)分析,得到其固有頻率和固有振型,并對模型進行諧響應分析.徐亞棟等[3]研究了復合材料火炮身管的動力學性能,在模態(tài)分析的基礎上,對火炮身管進行強迫響應分析,并與金屬材料的火炮身管動力學性能進行了比較.郝雷[4]對某火炮身管進行了有限元模態(tài)分析和試驗模態(tài)分析,提出通過改變身管相應部位的剛度來減小其振動情況.劉雷等[5]為考慮身管的彈性變形,采用有限質(zhì)量段的方法,將火炮身管離散成若干段剛體,研究了彈丸與身管的振動耦合問題.
文中針對火炮身管的振動變形特點,對坦克火炮身管進行有限元約束模態(tài)分析,并基于模態(tài)試驗測試系統(tǒng)對其進行試驗驗證;根據(jù)分析結(jié)果,為提高動力學分析的計算效率,建立了火炮身管的有限段振動簡化模型.
坦克在越野路面上行駛時,受到復雜的外部激擾作用,底盤將產(chǎn)生劇烈振動,并通過懸掛裝置、車體、炮塔等部件對火炮身管的振動產(chǎn)生影響[6].
火炮身管、炮尾和搖架連接成一體,通過炮塔上的耳軸旋轉(zhuǎn)付和液壓驅(qū)動機構(gòu)安裝在炮塔上,當坦克高速行駛時,底盤產(chǎn)生大幅的俯仰運動,當火炮相對坦克底盤的擺動超過最大俯仰角范圍(-6°~14°),此時炮控系統(tǒng)液壓缸閉鎖,耳軸處等效為一固定端,火炮身管與底盤聯(lián)結(jié)為一體.因此可將火炮身管近似簡化為一懸臂梁結(jié)構(gòu),簡化示意圖如圖1所示.
圖1 火炮身管簡化示意圖
圖1中,y為絕對坐標下火炮的振動位移;x為測量點到炮尾末端的距離;t為響應時間;xt為耳軸點到炮尾末端的距離;yt為耳軸處垂向振動位移;θ為火炮繞耳軸點的俯仰角位移.
考慮身管的彈性變形,火炮炮口的振動情況表達如下:
(1)
火炮身管是一連續(xù)彈性,理論上具有無數(shù)個自由度,具有無數(shù)階模態(tài),在進行模態(tài)分析時,需要考慮提取多少階模態(tài)是足夠的,針對這一問題,引入模態(tài)有效質(zhì)量,用以判斷系統(tǒng)模態(tài)的重要程度.火炮身管的振動是微幅振動,近似屬于時不變定常系統(tǒng),因此,根據(jù)經(jīng)典模態(tài)理論,可設系統(tǒng)的微分方程[6]為
(2)
式中:[M1]、[C1]、[K1]分別為系統(tǒng)質(zhì)量矩陣、阻尼矩陣、剛度矩陣;{θ1}為火炮身管的位移;{F}為激振力矩陣.
對上式進行特征值求解即可求得系統(tǒng)的固有頻率和固有振型φ.系統(tǒng)廣義質(zhì)量矩陣定義如下:
(3)
(4)
對于第i階模態(tài),其模態(tài)有效質(zhì)量meff,i定義如下:
(5)
所有模態(tài)有效質(zhì)量之和等于結(jié)構(gòu)的總質(zhì)量,通過模態(tài)有效質(zhì)量可以判斷各階模態(tài)的重要性.
針對圖1中的懸臂梁模型,設其長度為L,單位長度上的質(zhì)量為ρ,彈性模量為E,橫截面積為A.則可求出前四階模態(tài)頻率和模態(tài)有效質(zhì)量如表1所示.
表1 懸臂梁模態(tài)參數(shù)
從表1中可以看出,研究懸臂梁模型的振動問題,提取前五階模態(tài)即可,并且模態(tài)階數(shù)越低,模態(tài)有效質(zhì)量越大,前兩階模態(tài)的模態(tài)有效質(zhì)量占80.14%,說明這兩階模態(tài)對懸臂梁的彎曲振動貢獻量最大.
查閱相關資料,火炮身管模型的幾何參數(shù)、質(zhì)量屬性、材料屬性主要參數(shù)如表2所示.
表2 身管模型結(jié)構(gòu)尺寸和材料參數(shù)
根據(jù)火炮身管的主要參數(shù),并對身管連接處的倒角、螺紋等細節(jié)特征進行簡化處理,建立火炮身管的有限元模型如圖2所示.根據(jù)火炮身管的約束情況,在炮尾末端和耳軸處添加固定副約束,建立約束狀態(tài)下的火炮身管計算模型.
圖2 身管有限元網(wǎng)絡劃分模型
火炮身管,屬于基本對稱結(jié)構(gòu),振動主要沿著火炮軸向的縱向振動和垂直縱向軸線的橫向振動.由于火炮在射擊時,炮口的振動主要集中在高低平面內(nèi),因此文中考慮的主要是火炮橫向振動中的高低向振動,而不考慮水平向的振動[7].
求解得到各階模態(tài)的頻率如下表3所示,前5階振型如圖3所示.
圖3 身管前五階模態(tài)振型圖
為進一步分析各階模態(tài)的重要性,對火炮身管進行頻響分析.在模型上選取典型節(jié)點位置,沿高低向施加200 N的掃頻集中力,選取頻帶0~1 000 Hz,求解炮口振動的頻率響應[7],結(jié)果如圖4所示.
表3 火炮身管前5階固有頻率和振動周期
圖4 炮口振動頻率響應分析
由圖4分析可得,各階模態(tài)的模態(tài)參與因子沿固有頻率方向依次減小.經(jīng)相關計算可知,1階和2階模態(tài)對炮口振動響應的貢獻達85%,因此,在誤差允許的范圍內(nèi),3階以上的模態(tài)振動貢獻量可以忽略不計[7].
為了驗證有限元模態(tài)分析的結(jié)果,對火炮身管進行模態(tài)試驗測試.試驗方法為錘擊法,試驗裝置如圖5所示,包括LMS測試系統(tǒng)、加速度傳感器、力錘、筆記本電腦、電源等.
圖5 試驗裝置及傳感器布置
試驗工況:身管熱護套和排氣裝置拆除,炮口系統(tǒng)閉鎖,身管軸線仰角設置在很小角度,此時身管與坦克底盤固結(jié)在一起,處于完全約束狀態(tài).
試驗假設:僅考慮火炮身管的橫向振動特征[8],身管簡化為懸臂梁,假設梁各截面的中心主慣性矩在同一XOY平面內(nèi),外載荷也作用在同一平面.在低頻振動時,忽略剪切變形以及截面繞中心軸轉(zhuǎn)動慣量的影響.
試驗步驟:
(1)選取合適的測量點,對火炮身管進行取點編號;
(2)合理設置測試系統(tǒng)量程,調(diào)試測試儀器;
(3)由于傳感器數(shù)量和儀器通道數(shù)目的限制,不能一次測得所有點的響應信號,故采用分組多次測量的方式.
(4)試驗過程中,為保證測得的頻響函數(shù)真實可信,應實時監(jiān)測響應信號和激勵信號之間的相干函數(shù),對相干函數(shù)值低于0.9的測點需要剔除,重新進行試驗.
測得典型測量點54號點的振動頻響函數(shù)如圖6所示.
圖6 54號測量點高低向振動頻響函數(shù)
圖7為試驗模態(tài)穩(wěn)定圖,試驗采用PolyMAX法對模態(tài)參數(shù)進行識別.PolyMAX是基于系統(tǒng)模態(tài)穩(wěn)定圖從計算得到的一組模態(tài)結(jié)果中提取所要關注的模態(tài).
計算得到身管約束狀態(tài)的前5階固有頻率和阻尼比如表4所示,前四階振型如圖8所示.
圖7 試驗模態(tài)穩(wěn)定圖
表4 身管前5階模態(tài)頻率和阻尼比
階次模態(tài)頻率(Hz)阻尼比(%)備注110.9142.73橫向振動252.2304.08橫向振動3144.7352.33橫向振動4210.2622.97徑向振動5223.5560.3橫向振動
圖8 前四階模態(tài)振型圖
為驗證系統(tǒng)模態(tài)參數(shù)識別的效果,利用模態(tài)置信判據(jù)MAC矩陣進行各階模態(tài)的相關性分析.前5階模態(tài)的MAC矩陣如圖9所示,從圖中可知對角線值都接近100,非對角線值接近0,說明各階模態(tài)振型是相互獨立的.
圖9 前5階模態(tài)MAC矩陣分析圖
通過模態(tài)試驗得到火炮身管約束狀態(tài)下的固有頻率和主振型,取前3階有限元模態(tài)頻率與試驗結(jié)果對比如下表5所示.
表5 模態(tài)頻率計算結(jié)果對比
通過表5分析可得,有限元模態(tài)分析結(jié)果和試驗模態(tài)分析結(jié)果基本吻合,誤差均在10%以內(nèi),驗證了所建火炮身管有限元模型的準確性.
為了更加精確地描述火炮炮口的振動響應,剛性多體系統(tǒng)動力學建模理論過于簡單,僅將火炮整體視為剛性體不能很好的表達炮口的振動特征,而將火炮身管考慮成柔性體,進行動力學分析時又會帶來計算效率低,求解難度大的問題.因此,為解決以上問題,采用有限質(zhì)量段的方法對火炮身管模型進行簡化.
有限段方法是基于多剛體動力學理論發(fā)展起來的柔性體動力學研究方法,其主要思想是將連續(xù)柔性體結(jié)構(gòu)件離散成有限質(zhì)量段的剛體,各質(zhì)量段之間通過彈簧阻尼相連,連續(xù)體的柔性特征通過等效彈簧阻尼來表征[9-10].
根據(jù)火炮身管模態(tài)分析結(jié)果,在誤差允許的范圍內(nèi),忽略高階模態(tài)影響,用分段函數(shù)擬合身管歸一化的2階振型,建立火炮身管約束狀態(tài)下的2階線性振動系統(tǒng)模型[11].在實際工作過程中,坦克火炮身管一端懸臂,另一端和炮控系統(tǒng)的搖架和炮尾相連,近似為一個支撐懸臂梁結(jié)構(gòu).建立分段函數(shù):
(6)
利用MATLAB工具箱中的Curve Fitting函數(shù)擬合歸一化的階振型如圖10所示.
圖10 分段函數(shù)擬合歸一化的2階振型
分段函數(shù)表達式為
(7)
通過截點位置計算火炮三段有限段模型參數(shù)包括質(zhì)量、慣量、質(zhì)心位置、等效彈簧阻尼特性參數(shù).圖11為火炮有限段分段結(jié)果示意圖,模型分段計算結(jié)果如表6所示.
圖11 火炮有限段分段結(jié)果示意圖
圖11中,L1為炮尾段火炮身管的長度;L2為中間段火炮身管的長度;L3為炮口段火炮身管的長度;η1為炮尾段質(zhì)心到耳軸點的距離;η2為中間段質(zhì)心到界面1的距離;η3為炮口段質(zhì)心到界面2的距離;k12、c12分別為炮尾段與中間段之間的連接剛度和阻尼;k23、c12分別為炮口段與中間段之間的連接剛度和阻尼;m1,m2,m3分別為3個質(zhì)量段的質(zhì)量;J1,J2,J3分別為3個質(zhì)量段繞火炮身管軸線的轉(zhuǎn)動慣量.
結(jié)合所擬合的分段函數(shù)和模態(tài)振型,假設炮尾端質(zhì)量段固定不動,根據(jù)所建立的三質(zhì)量段火炮身管模型,對模型進行動力學分析,建立微分方程[6]:
(8)
連接處的阻尼矩陣根據(jù)經(jīng)驗公式計算:
[C]=β[K],β=0.000 1. (9)
由式8計算可得,前兩階固有頻率分別為10.88 Hz和52.77 Hz,與有限元前兩階固有頻率比較如表7.
表7 有限元模態(tài)頻率與有限段模態(tài)頻率對比
從表中可知,利用有限段方法對火炮身管進行簡化,建立的三質(zhì)量段模型的固有特性計算結(jié)果與有限元模型差別不大,說明有限段方法應用于火炮身管是有效的,所建立的三質(zhì)量段火炮身管模型為其動力學分析提供了基礎.
文中從火炮身管的固有特性出發(fā),將其簡化為一懸臂梁模型,建立了其有限元模型,分別進行了有限元模態(tài)分析和試驗模態(tài)分析,并基于模態(tài)分析結(jié)果,建立了火炮身管有限段2階質(zhì)量模型.具體結(jié)論如下:
1)根據(jù)火炮身管的實際約束狀態(tài),可將其簡化為一懸臂梁模型,模態(tài)分析中前五階模態(tài)的貢獻量最大;
2)火炮身管有限元模態(tài)分析結(jié)果和試驗模態(tài)分析結(jié)果相近,說明了所建立的火炮身管有限元模型的準確性;
3)火炮身管的低階模態(tài)振動主要是橫向振動,以前兩階的模態(tài)貢獻率最大,根據(jù)有限段方法建立的火炮身管2階質(zhì)量模型能夠較好地表征火炮身管的振動特性,為下一步仿真計算提供簡化模型基礎.