張培忠, 鄭建福, 米中賀, 陳國利, 寧金貴, 黃彥昌
(中國白城兵器試驗中心, 吉林 白城 137001)
以非滿管代替滿管進行轉(zhuǎn)管炮身管壽命試驗研究
張培忠, 鄭建福, 米中賀, 陳國利, 寧金貴, 黃彥昌
(中國白城兵器試驗中心, 吉林 白城 137001)
為減少轉(zhuǎn)管炮身管壽命試驗的用彈量,研究用非滿管代替滿管進行轉(zhuǎn)管炮身管壽命試驗的技術。建立炮口振動位移計算方程,利用高階傅里葉級數(shù)擬合炮膛合力,獲得炮膛合力的頻域函數(shù);利用模態(tài)實驗獲得了轉(zhuǎn)管炮的模態(tài)參數(shù),構(gòu)建轉(zhuǎn)管炮起落部分的頻響函數(shù);利用頻響函數(shù)計算了轉(zhuǎn)管炮在滿管射擊、非滿管射擊方式下的炮口振動位移。以滿管射擊、非滿管射擊方式下炮口振動位移的一致性為原則,確定身管壽命試驗方案。通過試驗驗證了滿管射擊、滿非管射擊方式下的身管壽命試驗結(jié)果具有一致性,從而獲得用彈量較少的身管壽命試驗方法,節(jié)省了試驗的彈藥消耗。
兵器科學與技術; 火炮; 身管壽命; 試驗技術
轉(zhuǎn)管炮是一種高射速防空/反導武器,其身管數(shù)量較多,有6~11根管,射擊時各個身管在旋轉(zhuǎn)中交替、循環(huán)射擊,在火炮無故障、無瞎火彈的情況下,各個身管累計射擊的彈數(shù)基本相等。在靶場定型試驗時必須進行身管壽命考核,按照傳統(tǒng)的身管壽命試驗方法,身管累積射擊彈數(shù)必須達到身管壽命指標規(guī)定的彈數(shù),再依照軍標規(guī)定的判據(jù)判斷身管壽命是否達到指標,因此試驗總用彈量=單管壽命指標彈數(shù)×身管數(shù)[1]。以某A型轉(zhuǎn)管炮為例,按此方法進行試驗,用彈量為2萬多發(fā),可見彈藥的消耗巨大。
為了減少彈藥的消耗,身管壽命試驗一般結(jié)合其他試驗項目一并進行,例如強度試驗、隨動精度測試、有效射程試驗、淋雨試驗、炮口沖擊波測試等,所用彈數(shù)全部計入身管壽命中。當全部計入的彈數(shù)仍然不能達到身管壽命指標數(shù)時,則補充一部分彈數(shù),進行磨損射擊試驗,直至達到身管壽命指標數(shù)為止。此方法的用彈量也很大,因此需要研究新方法以減少彈藥消耗。
在“試驗總用彈量=單管壽命指標彈數(shù)×身管數(shù)”公式中,單管壽命指標彈數(shù)是上級規(guī)定的需要考核的指標,在定型試驗中不可更改,因此只能采用更改試驗用身管數(shù)的方法,通過使用少部分身管(定義為非滿管)來代替全部身管(定義為滿管),以完成身管壽命試驗,從而減少試驗用彈量。
轉(zhuǎn)管炮在靶場進行定型試驗時,分為前期、中期和后期以及在必要時補充進行的磨損射擊試驗階段,共4個階段。
在定型試驗的前期和中期,主要進行初速檢查、彈帶性能測試、引信早炸與瞎火率測試、立靶密集度試驗、強度試驗、射速檢查、隨動精度測試、有效射程測試等項目,其中:初速檢查、彈帶性能測試、引信早炸與瞎火率測試的射擊方式采取每根身管獨立進行、單發(fā)射擊的方式,立靶密集度試驗、隨動精度測試等項目均采取全部身管(滿管)連發(fā)射擊的方式。射擊的彈數(shù)全部計入身管壽命中,到此階段結(jié)束時,每根身管累計射擊的彈數(shù)大約為800~1 000發(fā),處于壽命中期。
在定型試驗后期,首先要進行初速、彈帶性能、引信早炸與瞎火率、立靶密集度等檢查項目,依照國家軍用標準規(guī)定的判據(jù)分析身管壽命是否終止,若終止則停止試驗,否則繼續(xù)進行定型試驗剩余項目,采用滿管連發(fā)射擊方式進行隨動熱態(tài)射擊、炮口沖擊波測試、振動加速度測試、機構(gòu)動可靠性檢查、淋雨試驗等。射擊彈數(shù)計入身管壽命中,到此階段結(jié)束時,每根身管累計射擊的彈數(shù)大約為1 500~1 700發(fā),處于壽命后期,距離身管壽命指標數(shù)相差300~500發(fā),需要補充一部分彈藥,進行磨損射擊試驗。
在磨損射擊試驗階段,僅進行初速、彈帶性能、引信早炸與瞎火率、立靶密集度檢查,依照國家軍用標準規(guī)定的判據(jù)分析身管壽命是否終止,判據(jù)如下[1]:1)彈帶削光;2)引信連續(xù)瞎火2發(fā)以上,或者彈丸早炸;3)立靶密集度超過戰(zhàn)技指標的8倍,或者橫彈數(shù)量超過50%;4)彈丸距離散布增大1.5%(轉(zhuǎn)管炮無該指標);5)彈丸初速下降超過戰(zhàn)技指標規(guī)定的要求,其中任何一項超標均可判定為身管壽命終止,從而停止試驗。否則重復進行4個階段的試驗,直至每根身管累計射擊彈數(shù)達到身管壽命指標數(shù)為止。以某A型轉(zhuǎn)管炮為例,在磨損射擊試驗階段需要補充射擊彈數(shù)約3 000~5 000發(fā)。
轉(zhuǎn)管炮定型試驗的項目眾多,研究全部項目均采用以非滿管代替滿管進行射擊是很不現(xiàn)實的,因此本文僅研究在磨損射擊試驗階段用非滿管代替滿管射擊以節(jié)省彈藥的方法。
在磨損射擊試驗之前,各個身管的射擊彈數(shù)、彈藥種類、發(fā)射速度、維護保養(yǎng)等情況基本一致,各個身管的內(nèi)膛磨損程度也基本相等。因此,允許任意抽取其中部分身管作為樣本,進行獨立、單發(fā)射擊,各個身管之間互不影響,分別獲得初速、彈帶性能、引信早炸與瞎火率數(shù)據(jù),用以表征全部身管的初速、彈帶性能、引信早炸與瞎火率。但是,由于立靶密集度試驗采取連發(fā)射擊方式進行,各個身管之間相互影響,抽取的身管數(shù)量、各身管所處位置都會對轉(zhuǎn)管炮產(chǎn)生不同頻率的激勵,該頻率與轉(zhuǎn)管炮(尤其彈性較大的起落部分)固有頻率的接近(或者遠離)程度不同,引起的炮口振動位移也不同,最終造成立靶密集度試驗結(jié)果有差異。因此,用非滿管代替滿管進行立靶密集度試驗時,抽取的身管數(shù)量及其所處位置必須保證使非滿管射擊方式與滿管射擊方式獲得的立靶密集度結(jié)果具有一致性,同時保證初速、彈帶性能、引信早炸與瞎火率結(jié)果的一致性,由此獲得的身管壽命值才能與傳統(tǒng)方法獲得的身管壽命值等效。由于轉(zhuǎn)管炮使用內(nèi)能源或者內(nèi)+外能源驅(qū)動,轉(zhuǎn)速基本不變,在射擊過程中,每根身管旋轉(zhuǎn)一周才射擊一發(fā)彈,用非滿管代替滿管方式,參與射擊的身管自身的射速基本不變,而且每組每管射擊≤10發(fā)彈,組與組的間隔時間≥4 h[1],身管溫升也不因此而顯著改變,故其磨損程度不會因此而改變。
目前常用動力學仿真法求解火炮的發(fā)射過程[2-4],但是該方法需要經(jīng)過仿真建模、校模和驗模等一系列復雜過程,專業(yè)性很強,靶場難以掌握,且費用高、周期長,在武器定型試驗中重復應用較困難。為此,下面介紹頻響函數(shù)法,該方法利用模態(tài)實驗測試火炮的模態(tài)參數(shù)、構(gòu)建頻響函數(shù),建立頻域炮膛合力函數(shù),基于頻響函數(shù)與頻域炮膛合力函數(shù)計算火炮振動響應,費用少,周期短,便于靶場掌握和重復運用。
2.1 轉(zhuǎn)管炮的振動力學方程
設轉(zhuǎn)管炮為N個自由度的系統(tǒng),在物理坐標系中的方程[5]為
(1)
該系統(tǒng)的固有頻率矩陣和固有振型矩陣分別為
(2)
在模態(tài)坐標系中用Φ作為坐標空間的基向量矩陣,令X=Φq,q為模態(tài)坐標向量,將 (1) 式變換為
(3)
令:f=Fejωt,q=Qejωt,x=Xejωt,由(3)式[6]得
X=Φ(kr-ω2mr+jωcr)-1ΦTF=
ΦYrΦTF=HF,
(4)
將Yr展開為
(5)
系統(tǒng)的激勵輸入點I、響應輸出點O之間的頻響函數(shù)為
(6)
(7)
2.2 炮口橫向振動位移函數(shù)
轉(zhuǎn)管炮的炮口橫向振動位移(垂直于身管軸線的方位向、高低向振動位移)的大小直接關系到立靶密集度的大小,為了研究用非滿管代替滿管進行立靶密集度試驗的方法,并使兩者的結(jié)果一致,必須保證炮口橫向振動位移一致。
轉(zhuǎn)管炮的起落部分由兩個耳軸、搖架、后坐部分、兩個制退緩沖器組成,如圖1所示。后坐部分由回轉(zhuǎn)部分、馬達、供彈口和排殼口組成?;剞D(zhuǎn)部分由自動機、身管組構(gòu)成。起落部分的質(zhì)心C、后坐部分的質(zhì)心C1都在回轉(zhuǎn)部分的軸線L上,起落部分的質(zhì)心C至最下方身管的軸線L1的距離為e. 在設計上軸線L與軸線L1在1 km處相交,近似認為兩者平行,且垂直于轉(zhuǎn)管炮的方位回轉(zhuǎn)軸和高低回轉(zhuǎn)軸。
圖1 轉(zhuǎn)管炮起落部分受力分析Fig.1 Forces of revolving barrels cannon elevating part
轉(zhuǎn)管炮射擊是由處于最下方位置的身管完成的,起落部分受到的炮膛合力的作用線在最下方位置的身管軸線L1上,炮膛合力對起落部分的質(zhì)心C形成力矩。后坐部分因后坐運動而產(chǎn)生慣性力,其作用線經(jīng)過質(zhì)心C1、在后坐部分的軸線L上。起落部分依靠耳軸支撐、隨動系統(tǒng)穩(wěn)定,由兩個耳軸提供支撐力,抵制炮膛合力、后坐部分的慣性力和重力,由隨動系統(tǒng)提供方位方向、高低方向控制力矩。轉(zhuǎn)管炮起落部分受力分析如圖1所示。
轉(zhuǎn)管炮射擊時,回轉(zhuǎn)部分在馬達驅(qū)動下進行旋轉(zhuǎn),自動機從供彈口輸入炮彈,完成裝填、推彈、閉鎖、擊發(fā)、開鎖和抽殼等一系列動作,從排殼口排殼。當某一身管處于最下方位置時發(fā)射炮彈、產(chǎn)生后坐力,推動后坐部分沿著搖架的前套箍和后滑動導軌運動,搖架上的兩個制退緩沖器產(chǎn)生阻力,吸收后坐動能。起落部分動力學方程為
(8)
式中:Fz(s)為耳軸頻域支撐力;Fp(s)為頻域炮膛合力;Fg(s)為后坐部分的頻域慣性力;G為起落部分的重力;Ma(s)為隨動系統(tǒng)方位向頻域控制力矩;Me(s)為隨動系統(tǒng)高低向頻域控制力矩;M(s)為炮膛合力對起落部分產(chǎn)生的頻域力矩;MG為起落部分的重力矩。
視起落部分為柔性懸臂梁,由于主動力矩M(s)的作用而產(chǎn)生炮口振動位移y(s),即炮口橫向振動位移,M(s)的力學作用等效于在炮口處橫向作用一個力F(s),造成起落部分在炮口處的橫向振動位移y(s). 假設耳軸中心O至炮口的距離為l,則:
lF(s)=M(s)=eFp(s).
(9)
距離e、長度l是定值。由于Fp、l近似平行,且都垂直于轉(zhuǎn)管炮的方位回轉(zhuǎn)軸、高低回轉(zhuǎn)軸,由(9)式有:
(10)
式中:ea和ee分別為距離e在方位回轉(zhuǎn)軸、高低回轉(zhuǎn)軸上的分量;Fa和Fe分別為F(s)在方位回轉(zhuǎn)軸、高低回轉(zhuǎn)軸上的分量。
令jω=s[6-7],則炮口橫向振動頻域位移為
(11)
式中:X(s)、Y(s)分別為炮口方位方向、高低方向頻域振動位移;Ha(s)、He(s)分別為方位方向、高低方向頻響函數(shù)。
(11)式經(jīng)過拉普拉斯反變換得到時域炮口橫向振動位移:
(12)
3.1 建立炮膛合力頻域函數(shù)
為了計算(12)式的轉(zhuǎn)管炮炮口振動位移,根據(jù)(11)式,需要給出炮膛合力頻域函數(shù)。近似認為火炮發(fā)射過程中的炮膛合力是一系列等幅度周期性激勵,滿管射擊方式、每間隔1管射擊1發(fā)的非滿管射擊方式、每間隔2管射擊1發(fā)的非滿管射擊方式的炮膛合力激勵如圖2~圖4所示,3種射擊方式力的幅值相等,但炮膛合力作用的頻率不同。鑒于本文研究的問題是滿管射擊方式、非滿管射擊方式下炮口振動位移的差別問題,不關心計算結(jié)果與實際值的符合程度,為了便于計算,在此將炮膛合力按最大幅值Fpm做歸一化處理,然后利用高階傅里葉級數(shù)擬合炮膛合力,獲得級數(shù)形式的時域函數(shù)。經(jīng)過多次嘗試,確定用8階傅里葉級數(shù)來擬合轉(zhuǎn)管炮的炮膛合力較適宜,兩者之間的符合程度較高,曲線如圖5~圖7所示。
圖2 滿管射擊時的炮膛合力Fig.2 Resultant force in bore for full loading firing
圖3 每間隔1管射擊1發(fā)的炮膛合力Fig.3 Resultant force in bore for firing one round every one barrel
圖4 每間隔2管射擊1發(fā)的炮膛合力Fig.4 Resultant force in bore for firing one round every 2 barrels
圖5 8階傅里葉級數(shù)擬合滿管射擊時的炮膛合力曲線Fig.5 Resultant force in bore fit by eighth-order Fourier series for full loading firing
圖6 8階傅里葉級數(shù)擬合每間隔1管射擊1發(fā)的 炮膛合力曲線Fig.6 Resultant force in bore fit by eighth-order Fourier series for firing one round every one barrel
圖7 8階傅里葉級數(shù)擬合每間隔2管射擊1發(fā)的炮膛 合力曲線Fig.7 Resultant force in bore fit by eighth-order Fourier series for firing one round every 2 barrels
利用8階傅里葉級數(shù)擬合3種射擊方式的炮膛合力,形成的時域函數(shù)表達式分別為
Fp(t)=0.272 3-
0.000 153 31cos(464.7t)+2.324×10-5sin(464.7t)+
0.082 5cos(2×464.7t)+0.327 4sin(2×464.7t)-
0.000 167 8cos(3×464.7t)+7.527×10-5sin(3×464.7t)-
0.078 71cos(4×464.7t)+0.170 2sin(4×464.7t)-
0.000 204 5cos(5×464.7t)+0.000 146 9sin(5×464.7t)-
0.104 9cos(6×464.7t)+0.066 65sin(6×464.7t)-
0.000 292 8cos(7×464.7t)+0.000 298 1sin(7×464.7t)-
0.071 97cos(8×464.7t)-0.002 224sin(8×464.7t),
(13)
Fp(t)=0.170 1+
0.158 3cos(511.5t)+0.194 2sin(511.5t)+
0.029 1cos(2×511.5t)+0.170 6sin(2×511.5t)-
0.002 68cos(3×511.5t)+0.120 2sin(3×511.5t)-
0.045 93cos(4×511.5t)+0.080 43sin(4×511.5t)-
0.054 26cos(5×511.5t)+0.050 44sin(5×511.5t)-
0.050 33cos(6×511.5t)+0.022 74sin(6×511.5t)-
0.041 43cos(7×511.5t)+0.003 614sin(7×511.5t)-
0.029 14cos(8×511.5t)-0.004 327sin(8×511.5t),
(14)
Fp(t)=0.117 8+0.152 1cos(344t)+
0.119sin(344t)-0.069 78cos(2×344t)+
0.129 8sin(2×344t)-0.019 69cos(3×344t)+
0.114sin(3×344t)-0.008 87cos(4×344t)+
0.091 22sin(4×344t)-0.022 98cos(5×344t)+
0.070 2sin(5×344t)-0.029 9cos(6×344t)+
0.053 79sin(6×344t)-0.034 17cos(7×344t)+
0.040 27sin(7×344t)-0.036 17cos(8×344t)+
0.027 63sin(8×344t).
(15)
將以上炮膛合力時域函數(shù)經(jīng)過拉普拉斯變換Fp(s)=L[Fp(t)],得到3種射擊方式的炮膛合力頻域函數(shù),表達式分別為
Fp(s)=
0.272/s-0.153×10-3s/(s2+0.216×106)+
0.108×10-1/(s2+0.216×106)+0.825×
10-1s/(s2+0.864×106)+304/(s2+0.864×106)-
0.168×10-3s/(s2+0.194×107)+0.105/(s2+
0.194×107)-0.787×10-1s/(s2+0.346×107)+
316/(s2+0.346×107)-0.204×10-3s/(s2+0.540×
107)+0.348/(s2+0.540×107)-0.105s/(s2+
0.777×107)+186/(s2+0.777×107)-
0.293×10-3s/(s2+0.106×108)+
0.970/(s2+0.106×108)-0.720×10-1s/(s2+
0.138×108)-8.27/(s2+0.138×108),
(16)
Fp(s)=0.170/s+0.158s/(s2+0.262×106)+
99.3/(s2+0.262×106)+0.291×10-1s/(s2+
0.105×107)+175/(s2+0.105×107)-0.268×
10-1s/(s2+0.235×107)+184/(s2+
0.235×107)-0.459×10-1s/(s2+0.419×107)+
165/(s2+0.419×107)-0.543×10-1s/(s2+
0.654×107)+129/(s2+0.654×107)-0.530×
10-1s/(s2+0.942×107)+69.8/(s2+
0.942×107)-0.414×10-1s/(s2+0.128×108)+
12.9/(s2+0.128×108)-0.291×10-1s/(s2+
0.167×108)-17.7/(s2+0.167×108),
(17)
Fp(s)=0.12/s+0.15s/(s2+0.12×106)+
41/(s2+0.12×106)+0.07s/(s2+0.47×106)+
89/(s2+0.47×106)+0.02s/(s2+0.11×107)+
120/(s2+0.11×107)-0.0089s/(s2+
0.19×107)+130/(s2+0.19×107)-
0.023s/(s2+0.30×107)+120/(s2+
0.30×107-0.03/(s2+0.43×107)+
110/(s2+0.43×107)-0.034s/(s2+0.58×
107)+97/(s2+0.58×107)-0.036s/(s2+
0.76×107)+75/(s2+0.76×107).
(18)
3.2 構(gòu)建轉(zhuǎn)管炮的頻響函數(shù)
為了計算(12)式的轉(zhuǎn)管炮炮口振動位移,根據(jù)(11)式,需要給出頻響函數(shù)。采用力錘法進行轉(zhuǎn)管炮的模態(tài)實驗(傳感器布置見圖8),獲得模態(tài)參數(shù)[8],利用模態(tài)參數(shù)構(gòu)建因炮膛合力引起炮口振動的頻響函數(shù)??紤]到M(s)的作用等效于在炮口處橫向作用一個力F(s),炮口橫向振動的力激勵與位移響應處在同一點、同一方向,在測試高低方向的振動模態(tài)時,力錘敲擊選在身管組最上方的身管、靠近炮口套箍的位置垂直向下敲,如圖9所示。在測試方位方向的振動模態(tài)時,力錘敲擊選在身管組最左方位置的身管、靠近炮口套箍的位置水平向右敲。部分模態(tài)如圖10所示,測試結(jié)果如表1~表3所示。
圖8 傳感器布置圖Fig.8 Sensor layout
圖9 測試轉(zhuǎn)管炮的模態(tài)Fig.9 Measuring modal of revolving barrels cannon
忽略高階模態(tài)的影響,轉(zhuǎn)管炮的方位方向力激勵、位移響應之間的頻響函數(shù)近似取前5階:
(19)
轉(zhuǎn)管炮的高低方向力激勵、位移響應之間的頻響函數(shù)為近似取前4階:
(20)
圖10 轉(zhuǎn)管炮的高低向模態(tài)Fig.10 Elevation modals of revolving barrels cannon
表1 模態(tài)實驗測試結(jié)果
3.3 計算炮口振動位移
本文僅關注滿管射擊方式、非滿管射擊方式下炮口振動位移的差別,而不關注其實際振動位移值,因此距離e、長度l無需計算。經(jīng)過計算可得某A型轉(zhuǎn)管炮分別按滿管射擊方式、每間隔1管射擊1發(fā)的非滿管射擊方式、每間隔2管射擊1發(fā)的非滿管射擊方式下的炮口橫向振動位移,如圖11和圖12所示。
表2 高低方向模態(tài)留數(shù)
表3 方位方向模態(tài)留數(shù)
圖11 滿管射擊、每間隔1管射擊1發(fā)、每間隔2管 射擊1發(fā)的方式下炮口方位方向振動位移Fig.11 Azimuth vibration displacement of muzzle for full loading firing, firing one round every one barrel, and firing one round every 2 barrels
圖12 滿管射擊、每間隔1管射擊1發(fā)、每間隔2管 射擊1發(fā)的方式下炮口高低方向振動位移Fig.12 Elevation vibration displacements of muzzle for full loading firing, firing one round every one barrel, and firing one round every 2 barrels
從圖11和圖12看出:滿管射擊方式、每間隔1管射擊1發(fā)的方式引起炮口的方位方向、高低方向振動幅值沒有顯著差別,因此理論上可以利用每間隔1管射擊1發(fā)的方式代替滿管射擊方式,進行某A型轉(zhuǎn)管炮立靶密集度試驗。而滿管射擊方式、每間隔2管射擊1發(fā)的方式引起炮口的方位方向、高低方向振動幅值有著明顯差別,理論上不能利用每間隔2管射擊1發(fā)的方式代替滿管射擊方式進行立靶密集度試驗,這是因為該方式的射擊頻率(3 000發(fā)/min)接近了轉(zhuǎn)管炮2階模態(tài)的固有頻率(47.9 Hz)。
某B型轉(zhuǎn)管炮分別按滿管射擊方式、每間隔1管射擊1發(fā)的非滿管射擊方式下的炮口橫向振動位移的計算結(jié)果如圖13和圖14所示。
圖13 滿管射擊、每間隔1管射擊1發(fā)的方式下炮口 方位方向振動位移Fig.13 Azimuth vibration displacements of muzzle for full loading firing, firing one round every one barrel, and firing one round every 2 barrels
圖14 滿管射擊、每間隔1管射擊1發(fā)的方式下 炮口高低方向振動位移Fig.14 Elevation vibration displacements of muzzle for full loading firing, firing one round every one barrel, and firing one round every 2 barrels
從圖13和圖14看出:滿管射擊方式、每間隔1管射擊1發(fā)的方式引起炮口的方位方向、高低方向振動幅值沒有明顯的差別,因此可以利用每間隔1管射擊1發(fā)的方式代替滿管射擊方式,進行某B型轉(zhuǎn)管炮的立靶密集度試驗。
首先在某A型轉(zhuǎn)管炮定型試驗的最后階段進行,此時該炮每根身管已經(jīng)累計射擊1 248發(fā),距離身管壽命指標還有752發(fā),轉(zhuǎn)入磨損射擊階段。在磨損射擊階段,選擇身管1、身管3、身管5、身管7等(也可以選擇身管2、身管4、身管6、身管8等)進行脫殼穿甲彈初速檢查、榴彈彈帶性能檢查、榴彈引信早炸與瞎火率測試以及脫殼穿甲彈立靶密集度試驗,裝填方法如下:用搖柄手動反向旋轉(zhuǎn)自動機,將彈鼓與自動機聯(lián)動,在被選定的身管對應的撥彈齒中裝入實彈,其他位置裝教練彈,將彈旋轉(zhuǎn)進入彈鼓中。當火炮射擊時,彈鼓與自動機聯(lián)動正向旋轉(zhuǎn),將彈送入自動機,被選定的身管射擊。其中前3項為每管獨立、單發(fā)射擊(工況1),立靶密集度為連發(fā)射擊,按2種方式實施:滿管射擊(工況2)的立靶密集度1組、每間隔1根身管射擊1發(fā)的非滿管連發(fā)射擊(工況3)1組,結(jié)果如表4和表5所示。
表4 某A型轉(zhuǎn)管炮單發(fā)射擊方式的試驗結(jié)果
從表4可以看出,在工況1中,滿管射擊方式、每間隔1根身管射擊1發(fā)的非滿管射擊方式下,脫殼穿甲彈初速、榴彈彈帶性能、榴彈引信早炸與瞎火率測試結(jié)果具有一致性。
表5 某A型轉(zhuǎn)管炮兩種射擊方式的立靶密集度試驗結(jié)果
利用F檢驗法[4,9],給定顯著性水平α=0.10,自由度ν1=54,ν2=54,查表有Fα/2(ν1,ν2)=1.433、F1-α/2(ν2,ν1)=0.697 8. 方位方向F(ν1,ν2)=0.610/0.494=1.234,高低方向F(ν1,ν2)=0.941/0.884=1.064.
由于F1-α/2(ν2,ν1)≤F(ν1,ν2)≤Fα/2(ν1,ν2),可見兩種立靶密集度試驗結(jié)果具有一致性,且立靶上均無橫彈,從而驗證了用非滿管代替滿管進行某A型轉(zhuǎn)管炮的身管壽命試驗是可行的。
在某B型轉(zhuǎn)管炮定型試驗的最后階段進行驗證試驗,此時該炮每根身管已經(jīng)累計射擊1 667發(fā),距離身管壽命指標還有333發(fā)。在磨損射擊階段,選擇身管1、身管3、身管5以工況1進行榴彈初速檢查、榴彈彈帶性能測試、榴彈引信早炸與瞎火率測試,再以工況2進行6組立靶密集度試驗,以工況3進行2組立靶密集度試驗,結(jié)果如表6~表8所示。
表6 某B型轉(zhuǎn)管炮單發(fā)射擊方式的試驗結(jié)果
從表6可以看出,在工況1中,滿管射擊方式、每間隔1根身管射擊1發(fā)的非滿管射擊方式下,榴彈初速、榴彈彈帶性能、榴彈引信早炸與瞎火率測試結(jié)果具有一致性。
利用F檢驗,給定顯著性水平α=0.10,自由度ν1=14,ν2=29,查表有:Fα/2(ν1,ν2)=1.75、F1-α/2(ν2,ν1)=0.57. 方位方向F(ν1,ν2)=
表7 某B型轉(zhuǎn)管炮滿管射擊方式的立靶密集度試驗結(jié)果
表8 某B型轉(zhuǎn)管炮非滿管射擊方式的立靶密集度試驗結(jié)果
由于F1-α/2(ν2,ν1)≤F(ν1,ν2)≤Fα/2(ν1,ν2),可見兩種立靶密集度的試驗結(jié)果具有一致性,且立靶上均無橫彈,從而驗證了用非滿管代替滿管進行某B型轉(zhuǎn)管炮的身管壽命試驗是可行的。
實際運用本方法在A、B兩型轉(zhuǎn)管炮的試驗中,共計節(jié)省彈藥3 759發(fā)。
本文研究了用非滿管代替滿管進行轉(zhuǎn)管炮身管壽命試驗、減少試驗用彈量的原理,分析了炮膛合力引起炮口振動的機理;利用高階傅里葉級數(shù)擬合炮膛合力,利用模態(tài)實驗測試模態(tài)參數(shù),構(gòu)建了頻響函數(shù),運用頻響函數(shù)、炮膛合力頻域函數(shù)計算了轉(zhuǎn)管炮射擊過程中的炮口振動位移;從理論上分析了轉(zhuǎn)管炮以非滿管代替滿管進行立靶密集度試驗的可行性,并在兩種轉(zhuǎn)管炮上進行了驗證試驗,獲得了初速、彈帶性能、引信早炸與瞎火率、立靶密集度試驗結(jié)果。研究所得結(jié)論如下:
1)利用高階傅里葉級數(shù)擬合轉(zhuǎn)管炮的炮膛合力,經(jīng)拉普拉斯變換后,得到炮膛合力頻域函數(shù)。
2)利用模態(tài)實驗與分析技術,獲得轉(zhuǎn)管炮的模態(tài)頻率、阻尼和留數(shù),構(gòu)建因炮膛合力引起炮口振動的頻響函數(shù),運用頻響函數(shù)與炮膛合力頻域函數(shù)的積得到炮口振動位移的頻域函數(shù)。
3)比較轉(zhuǎn)管炮在不同射擊方式下的炮口振動位移差別,以炮口振動位移一致性為原則,可以獲得具有一致性的立靶密集度試驗方案。
4)在轉(zhuǎn)管炮試驗的磨損射擊階段,利用非滿管代替滿管進行初速檢查、彈帶性能檢查、引信早炸與瞎火率測試、立靶密集度試驗,完成身管壽命試驗,可以顯著節(jié)省彈藥。
References)
[1] GJB349.23—1990常規(guī)兵器定型試驗方法:高射炮和海軍炮[S].北京:國防科學技術工業(yè)委員會,1990:10,21,29.
GJB349.23—1990 Conventional weapon approval test method : antiaircraft gun and navy gun[S]. Beijing: the Committee of Science, Technology and Industry for National Defense,1990: 10,21,29.(in Chinese)
[2] 米中賀,王毅,李戰(zhàn)旗,等. 轉(zhuǎn)管武器炮膛合力在ADAMS中的實現(xiàn)方法[J].火炮發(fā)射與控制學報,2011,32(4):46-49.
MI Zhong-he,WANG Yi, LI Zhan-qi, et al. Realization method of revolving barrels cannon resultant force based on ADAMS software [J]. Journal of Gun Launch and Control, 2011,32(4):46-49.(in Chinese)
[3] 楊明華,薛靜,韓小平,等.某高射速多管炮剛體動力學建模仿真[J].四川兵工學報,2011,32(12): 24-28.
YANG Ming-hua,XUE Jing, HAN Xiao-ping, et al. Dynamic simulation analyzing for high speed revolving barrels cannon[J]. Sichuan Ordnance Journal,2011,32(12) : 24-28.(in Chinese)
[4] 芮筱亭,陸毓琪,王國平,等.多管火箭發(fā)射動力學仿真與試驗測試方法[M].北京:國防工業(yè)出版社,2001: 188-189,217-235.
RUI Xiao-ting, LU Yu-qi, WANG Guo-ping,et al. Simulation and measurement method for multi-barrels rocket gun[M]. Beijing: National Defense Institute Press,2001:188-189,217-235.(in Chinese)
[5] 廖日東.I-DEAS實例教程[M].北京:北京理工大學出版社,2003:97-100.
LIAO Ri-dong.Course of I-DEAS examples [M].Beijing: Beijing Institute of Technology Press,2003:97-100.(in Chinese)
[6] 李德葆,陸秋海. 實驗模態(tài)分析及其應用[M].北京:科學出版社,2001:62-66,87-91.
LI De-bao, LU Qiu-hai. Experiment model analysis and utility[M].Beijing: Science Press,2001:62-66,87-91.(in Chinese)
[7] 沃德·海倫,斯蒂芬·拉門茲,波爾·薩斯.模態(tài)分析理論與試驗[M].白化同,李繼忠, 譯. 北京:北京理工大學出版社,2003: 5-6.
Heylen W, Lammens S, Szasz B. Modal analysis theory and experiment[M]. BAI Hua-tong, LI Ji-zhong, translated. Beijing:Beijing Institute of Technology Press,2003:5-6.(in Chinese)
[8] 張海航,狄長安.某型多管火炮的振動模態(tài)分析[J].兵工學報,2008,29(12): 1514-1517.
ZHANG Hai-hang, DI Chang-an. Analysis of model of a multi barrel cannon[J]. Acta Armamentarii,2008,29(12):1514-1517.(in Chinese)
[9] 張學鋒,馬大為,樂貴高,等.對某型多管火箭武器減少試驗用彈量的仿真[J].兵工學報,2006,27(5): 792-796.
ZHANG Xue-feng, MA Da-wei, LE Gui-gao, et al. Simulation on decreasing ammunition consumption for test of some multile launch rocket system[J]. Acta Armamentarii,2006,27(5):792-796.(in Chinese)
ResearchonBarrelLifeTestofRevolvingBarrelCannonbySubstitutingFullChargeLoadingwithNon-fullChargeLoading
ZHANG Pei-zhong, ZHENG Jian-fu, MI Zhong-he, CHEN Guo-li, NING Jin-gui, HUANG Yan-chang
(China Baicheng Ordnance Test Center, Baicheng 137001, Jilin, China)
The technology of substituting full charge loading with non-full charge loading is researched to decrease the ammunition consumption in the barrel life test of revolving barrel cannon. The computing equations of muzzle vibration displacement are established. High-order Fourier series is used to fit the resultant force in bore, and then the frequency equations of resultant force are established. The modal parameters of cannon are obtained by model measuring and analysis, and a frequency response equation of the cannon elevating part is established. The frequency response equation is used to calculate the muzzle vibration displacements in the way of non-full charge loading and full charge loading. The differences among these muzzle vibration displacements are compared. Based on the law that the muzzle vibration displacements of non-full charge loading and full charge loading are consistent, the barrel life test scheme is confirmed. In verification experiment, the barrel life’s consistency is proved. The barrel life test method of less amount ammunition is designed to decrease the ammunition consumption in test.
ordnance science and technology; cannon; barrel life; test technology
2017-02-14
總裝備部試驗技術研究項目(2010SY4205001)
張培忠(1966—),男,研究員,博士。E-mail:zhangpeizhong9966@163.com
寧金貴(1981—),男,工程師,碩士。E-mail:njgnjgnjg19810118@163.com
TJ306+.1
A
1000-1093(2017)10-1899-10
10.3969/j.issn.1000-1093.2017.10.004