李 榮,楊國(guó)來(lái),孫全兆,蕭 輝

(南京理工大學(xué) 機(jī)械工程學(xué)院,江蘇 南京 210094)

火炮后坐阻力的區(qū)間不確定性?xún)?yōu)化研究

李 榮,楊國(guó)來(lái),孫全兆,蕭 輝

(南京理工大學(xué) 機(jī)械工程學(xué)院,江蘇 南京 210094)

為減小火炮炮架受力,獲得更優(yōu)后坐阻力規(guī)律,進(jìn)行了后坐流液孔尺寸的區(qū)間不確定性?xún)?yōu)化研究。結(jié)合公差帶定義了三參數(shù)區(qū)間,構(gòu)建了基于三參數(shù)區(qū)間的不確定性?xún)?yōu)化模型及算法;以流液孔的幾何尺寸為不確定性變量,利用上述方法對(duì)節(jié)制桿形狀尺寸進(jìn)行了區(qū)間優(yōu)化設(shè)計(jì);優(yōu)化求解了不同公差等級(jí)的節(jié)制桿外徑。研究結(jié)果表明:優(yōu)化后的后坐阻力峰值區(qū)間更小,后坐阻力曲線更加平緩,充滿(mǎn)度更好。實(shí)現(xiàn)了節(jié)制桿外徑尺寸設(shè)計(jì)與公差設(shè)計(jì)的同步優(yōu)化,得到了兼顧工藝性和后坐阻力的節(jié)制桿外徑設(shè)計(jì)方案。

火炮;后坐阻力;結(jié)構(gòu)設(shè)計(jì);區(qū)間優(yōu)化;尺寸不確定性

通過(guò)設(shè)計(jì)火炮反后坐裝置的制退機(jī)液壓阻力可得到后坐阻力的規(guī)律,它極大地減小了火炮炮架受力。典型節(jié)制桿式制退機(jī)的液壓阻力是靠制退機(jī)活塞擠壓液體高速流過(guò)流液孔而產(chǎn)生,故后坐阻力規(guī)律實(shí)際是由流液孔面積的改變實(shí)現(xiàn)的[1]。宗士增[2]考慮了火炮后坐和全炮運(yùn)動(dòng)的相互影響,建立了火炮反后坐裝置動(dòng)力學(xué)耦合模型,以節(jié)制桿尺寸為設(shè)計(jì)變量進(jìn)行了結(jié)構(gòu)優(yōu)化。周成[3]以后坐阻力曲線與擬定的后坐力曲線的符合程度、后坐長(zhǎng)與設(shè)計(jì)值的接近程度為目標(biāo)函數(shù),建立以節(jié)制桿尺寸為設(shè)計(jì)變量的優(yōu)化模型。這些文獻(xiàn)所建立的模型都是確定性的模型,即認(rèn)為材料特性、載荷、幾何尺寸等參數(shù)在優(yōu)化過(guò)程中不存在攝動(dòng),忽略了參數(shù)的不確定性對(duì)目標(biāo)性能的影響[4]。然而,加工和裝配得到的流液孔尺寸是不確定的,它會(huì)影響火炮的后坐阻力和后坐長(zhǎng),因此需要研究流液孔尺寸的不確定性。

區(qū)間優(yōu)化方法僅需知道不確定性變量的上下限,因此,采用區(qū)間優(yōu)化理論及算法來(lái)研究不確定性問(wèn)題[5]。Jiang C[6]定義了公差指標(biāo),通過(guò)基于可靠度的區(qū)間可能度模型,提出了考慮公差的區(qū)間優(yōu)化方法,同時(shí)優(yōu)化得到了名義尺寸和公差。謝慧超[7]基于區(qū)間分析方法,建立了一種汽車(chē)懸架平順性的不確定性?xún)?yōu)化模型,在保證平順性的前提下使設(shè)計(jì)參數(shù)的對(duì)稱(chēng)公差最大化。上述區(qū)間優(yōu)化方法采用區(qū)間中點(diǎn)和半徑的形式來(lái)表征既是區(qū)間變量又是設(shè)計(jì)變量的量,但在機(jī)械工程實(shí)際中,尺寸常采用公差帶來(lái)表征,區(qū)間中點(diǎn)和半徑的形式使得優(yōu)化結(jié)果難以應(yīng)用。

針對(duì)火炮炮架受力大的問(wèn)題,從后坐流液孔幾何尺寸偏差的角度進(jìn)行了區(qū)間優(yōu)化研究。先結(jié)合公差帶定義了三參數(shù)區(qū)間,建立了三參數(shù)區(qū)間的不確定性?xún)?yōu)化模型及算法;然后以節(jié)制桿各段折點(diǎn)直徑、節(jié)制環(huán)內(nèi)徑、制退桿活塞與制退筒的間隙和節(jié)制環(huán)與節(jié)制桿起始位置的軸向裝配誤差作為不確定性變量,采用上述優(yōu)化方法進(jìn)行了后坐阻力的優(yōu)化,并定量評(píng)估了后坐阻力曲線的充滿(mǎn)度和平緩度。

1 三參數(shù)區(qū)間的不確定性?xún)?yōu)化

1.1 三參數(shù)區(qū)間定義

考慮結(jié)構(gòu)尺寸由加工和測(cè)量引起的不確定性,合格產(chǎn)品的尺寸會(huì)在設(shè)計(jì)的公差帶范圍內(nèi)攝動(dòng)。公差帶是上極限尺寸和下極限尺寸所限定的一個(gè)范圍,它的數(shù)學(xué)本質(zhì)是區(qū)間。公差帶一般有3種形式,區(qū)間一般有2種形式,如表1,表中Dmin為下極限尺寸,Dmax為上極限尺寸,D為公稱(chēng)尺寸,es為上極限偏差,ei為下極限偏差,ITn為公差等級(jí),xL為區(qū)間下限,xR為區(qū)間上限,xC為區(qū)間中點(diǎn),xC=(xL+xR)/2,xW為區(qū)間半徑,xW=(xR-xL)/2。

表1 公差帶和區(qū)間的形式及示例

公稱(chēng)尺寸是選型的標(biāo)準(zhǔn),而且確定性設(shè)計(jì)方法是以公稱(chēng)尺寸為基礎(chǔ);公差可表示加工的精度,還影響著零部件的成本。因此,機(jī)械工作者關(guān)心的尺寸參數(shù)是公稱(chēng)尺寸和公差。表1中公差帶的后2種形式均包含3個(gè)參數(shù),同時(shí)能有效地表示上述2個(gè)特征,然而現(xiàn)有兩參數(shù)的區(qū)間形式均不符合設(shè)計(jì)要求,其中區(qū)間中點(diǎn)和半徑的形式只是表1中公差帶第2種形式的特例。區(qū)間形式限制了區(qū)間優(yōu)化方法對(duì)不確定性結(jié)構(gòu)尺寸的研究。

由于公差大小不僅與公差等級(jí)相關(guān),還與公稱(chēng)尺寸的大小相關(guān),同時(shí)還是離散的,本文以公稱(chēng)尺寸和公差等級(jí)為基礎(chǔ)構(gòu)造三參數(shù)區(qū)間。

綜上,三參數(shù)區(qū)間記為

xI=〈xg,xIT n,xθ〉={x|xg-xθxt

xt,x、xg、xθ∈R}

(1)

1.2 三參數(shù)區(qū)間不確定性?xún)?yōu)化模型及算法

1.2.1 三參數(shù)區(qū)間設(shè)計(jì)變量轉(zhuǎn)化

區(qū)間優(yōu)化理論建立在兩參數(shù)的區(qū)間變量基礎(chǔ)上,故需將三參數(shù)區(qū)間轉(zhuǎn)化為兩參數(shù)區(qū)間。一般公差帶的實(shí)際約束如下:

①公差等級(jí)要求。以光軸為例,規(guī)定光軸的外圓公差等級(jí)為ITn,即xIT n為常數(shù)。

②行業(yè)習(xí)慣。如在鑄造生產(chǎn)中,行業(yè)習(xí)慣標(biāo)注對(duì)稱(chēng)公差,即xθ=0.5。

因此,三參數(shù)區(qū)間設(shè)計(jì)變量可通過(guò)公差帶實(shí)際等式約束h(x)=0轉(zhuǎn)化為兩參數(shù)區(qū)間。

1.2.2 目標(biāo)函數(shù)和約束函數(shù)的確定性轉(zhuǎn)化

采用罰函數(shù)將目標(biāo)和約束區(qū)間處理為增廣目標(biāo)區(qū)間,那么目標(biāo)和約束函數(shù)的確定性轉(zhuǎn)化的實(shí)質(zhì)是增廣目標(biāo)區(qū)間的排序。先采用姜潮[4]提出的區(qū)間可能度模型進(jìn)行區(qū)間的兩兩比較,然后建立可能度矩陣,并用徐澤水[8]提出的排序公式進(jìn)行排序。區(qū)間可能度定量描述了一區(qū)間優(yōu)于另一區(qū)間的具體程度,模型如下:

(2)

式中:AI和BI為區(qū)間,AR和AL分別為區(qū)間AI的上、下限,BR和BL分別為區(qū)間BI的上、下限。

(3)

1.2.3 三參數(shù)區(qū)間不確定性?xún)?yōu)化模型

通過(guò)上述確定性轉(zhuǎn)化,一般的優(yōu)化模型可轉(zhuǎn)化為如下三參數(shù)區(qū)間不確定性?xún)?yōu)化模型:

(4)

式中:f為目標(biāo)函數(shù),X為設(shè)計(jì)向量,m為設(shè)計(jì)變量的個(gè)數(shù),g(XI)為不等式約束函數(shù)。

1.2.4 算法流程

上述不確定性?xún)?yōu)化問(wèn)題的求解一般是兩層嵌套優(yōu)化問(wèn)題,其中外層用于設(shè)計(jì)向量的尋優(yōu),內(nèi)層用于計(jì)算不確定目標(biāo)和約束區(qū)間。本文外層采用NSGA-Ⅱ遺傳算法,為兼顧計(jì)算效率和精度,在不確定程度較小的情況下,內(nèi)層采用對(duì)稱(chēng)拉丁超立方試驗(yàn)設(shè)計(jì)方法(SLHS)[9]計(jì)算目標(biāo)和約束區(qū)間,內(nèi)層不再需不斷尋優(yōu),因此兩層嵌套優(yōu)化問(wèn)題被解耦為單層優(yōu)化問(wèn)題,算法流程如圖2所示。

1.3 懸臂梁算例

如圖3所示矩形截面懸臂梁結(jié)構(gòu)[6],長(zhǎng)度L=1 000mm,截面尺寸為b和h,頂端受水平載荷Fx=50kN和垂直向下載荷Fy=25kN的作用。

懸臂梁最大應(yīng)力σ在固定端,計(jì)算公式為

(5)

以截面尺寸b和h為設(shè)計(jì)變量,σ滿(mǎn)足許用應(yīng)力[σ]=250 MPa的約束,以體積V最小為優(yōu)化目標(biāo),建立確定性?xún)?yōu)化模型。

考慮鑄造懸臂梁的制造誤差,有3個(gè)不確定變量b、h和L,三者公差等級(jí)相同,均采用鑄造公差標(biāo)準(zhǔn)(GB/T6414-1999),且根據(jù)鑄造行業(yè)習(xí)慣,無(wú)特殊要求的尺寸均用對(duì)稱(chēng)公差形式。采用三參數(shù)區(qū)間不確定性?xún)?yōu)化方法建立優(yōu)化模型:

(6)

根據(jù)砂型鑄造大批量生產(chǎn)達(dá)到的公差等級(jí)范圍,懸臂梁鑄造公差等級(jí)選定CT8、CT9、CT10、CT11、CT12、CT13和CT14分別進(jìn)行優(yōu)化求解,并采用遺傳算法進(jìn)行確定性?xún)?yōu)化,尺寸圓整后的結(jié)果均列入表2。表中,ΔV為懸臂梁體積攝動(dòng)百分率,σI為懸臂梁固定端應(yīng)力,Δσ為懸臂梁固定端應(yīng)力攝動(dòng)百分率。

表2 確定性?xún)?yōu)化結(jié)果和區(qū)間不確定性?xún)?yōu)化結(jié)果對(duì)比

從表2可知,隨著公差等級(jí)的增大,懸臂梁體積和固定端應(yīng)力的攝動(dòng)也在增大,最大分別達(dá)到22.1%和32.2%,表明尺寸的不確定性在設(shè)計(jì)時(shí)必須加以考慮。與確定性?xún)?yōu)化結(jié)果對(duì)比,不確定性?xún)?yōu)化值偏大,因?yàn)樗紤]了結(jié)構(gòu)尺寸的攝動(dòng)對(duì)懸臂梁體積和應(yīng)力的影響,具有更好的穩(wěn)健性。體積最小為輕量化要求,公差等級(jí)與制造工藝性相關(guān)。如果設(shè)計(jì)者對(duì)懸臂梁質(zhì)量較關(guān)注,可以選擇第1組和第2組數(shù)據(jù);如果對(duì)工藝性較關(guān)注,可以選擇第6組和第7組這兩組;如果需綜合考慮質(zhì)量和工藝性,則可選擇第3組、第4組和第5組這三組方案。

懸臂梁算例表明,三參數(shù)區(qū)間的不確定性?xún)?yōu)化聯(lián)合了結(jié)構(gòu)尺寸設(shè)計(jì)與公差設(shè)計(jì),得到了截面尺寸及其公差帶,并得到穩(wěn)健性更好的設(shè)計(jì)方案。

2 基于三參數(shù)區(qū)間的后坐阻力優(yōu)化

2.1 火炮后坐受力分析

火炮射擊時(shí),其受力十分復(fù)雜,取后坐部分為研究對(duì)象,在后坐過(guò)程中,主動(dòng)力Fpt是火炮發(fā)射時(shí)產(chǎn)生后坐運(yùn)動(dòng)的原因,而后坐阻力FR是阻止后坐運(yùn)動(dòng)的原因。后坐阻力包括制退機(jī)的液壓阻力Fh、復(fù)進(jìn)機(jī)力Ff和它們緊塞裝置的摩擦力F、搖架導(dǎo)軌摩擦力FT、后坐部分重力沿炮膛軸線分力mhgsinφ。則后坐運(yùn)動(dòng)微分方程為

(7)

2.2 后坐反面問(wèn)題計(jì)算模型驗(yàn)證與分析

加工和裝配是流液孔幾何尺寸產(chǎn)生不確定性的重要原因。以節(jié)制桿各段折點(diǎn)直徑dk、節(jié)制環(huán)內(nèi)徑dp、制退桿活塞與制退筒間隙Δδ1、節(jié)制環(huán)與節(jié)制桿起始位置的軸向裝配誤差Δl作為不確定性變量,其中前兩者主要由加工引起,后兩者主要由裝配導(dǎo)致。典型節(jié)制桿式制退機(jī)結(jié)構(gòu)如圖4所示。

以某火炮為例,在保持火炮系統(tǒng)基本參數(shù)和反后坐裝置其他結(jié)構(gòu)尺寸不變的基礎(chǔ)上,流液孔尺寸的不確定性變量范圍采用該火炮的設(shè)計(jì)參數(shù),如表3所示。

表3 流液孔尺寸的不確定參數(shù)范圍

以不確定性變量的公稱(chēng)尺寸為初始值,Δδ1=0.315mm為初始值,Δl=0為初始值,運(yùn)用對(duì)稱(chēng)拉丁超立方試驗(yàn)設(shè)計(jì)方法在總體中抽500個(gè)樣本。對(duì)每個(gè)樣本,用式(7)建立的后坐運(yùn)動(dòng)微分方程在正常射擊條件下進(jìn)行后坐反面問(wèn)題仿真。得到后坐阻力峰值FRmax的初始值、最小值和最大值分別為190 982N,182 391N和197 072N,后坐阻力峰值區(qū)間相對(duì)其初始值的變化率為7.69%,得到后坐長(zhǎng)λ的初始值、最小值和最大值分別為865.8mm、853.2mm和889.3mm,后坐長(zhǎng)λ的最小值和最大值相差36.1mm。

該火炮設(shè)計(jì)后坐長(zhǎng)為910mm,后坐長(zhǎng)λ的仿真值中最大為889.3mm,其與設(shè)計(jì)值的相對(duì)誤差僅為2.27%,說(shuō)明了仿真模型的可信性。

計(jì)算結(jié)果表明,計(jì)及流液孔尺寸不確定性,后坐阻力峰值波動(dòng)為7.69%,后坐長(zhǎng)波動(dòng)為36.1mm。在傳統(tǒng)的確定性后坐阻力優(yōu)化中沒(méi)有考慮流液孔尺寸的不確定性帶來(lái)的目標(biāo)和約束的攝動(dòng),這會(huì)導(dǎo)致優(yōu)化結(jié)果魯棒性差。同時(shí)傳統(tǒng)優(yōu)化不能將結(jié)構(gòu)尺寸設(shè)計(jì)與公差設(shè)計(jì)聯(lián)合求解。因此,需進(jìn)行計(jì)及流液孔尺寸不確定性的區(qū)間優(yōu)化研究。

2.3 計(jì)及流液孔尺寸不確定性的后坐阻力優(yōu)化

本文通過(guò)對(duì)節(jié)制桿形狀進(jìn)行優(yōu)化來(lái)調(diào)整流液孔面積變化規(guī)律,從而實(shí)現(xiàn)后坐阻力的優(yōu)化。

本文在后坐長(zhǎng)基本不變的情況下減小后坐阻力峰值,根據(jù)2.2節(jié)后坐長(zhǎng)計(jì)算結(jié)果,故將后坐長(zhǎng)λ約束在λ0=895mm內(nèi)。同時(shí)為避免節(jié)制桿根部與節(jié)制環(huán)卡滯,節(jié)制桿與節(jié)制環(huán)配合間隙不小于Δδ2。

在第1節(jié)的基礎(chǔ)上,建立基于三參數(shù)區(qū)間的流液孔尺寸不確定的后坐阻力優(yōu)化模型:

(8)

式中:dk,min和dk,max分別為節(jié)制桿各段折點(diǎn)初始公稱(chēng)直徑的最小值和最大值,λI為考慮流液孔尺寸不確定后得到的后坐長(zhǎng)區(qū)間。

2.4 后坐阻力曲線評(píng)估

后坐阻力不僅要求其峰值小,還要求后坐阻力隨后坐行程變化的曲線(FR-s)充滿(mǎn)度好和更加平緩。目前在后坐阻力曲線的評(píng)價(jià)上,后坐阻力峰值采用定量評(píng)估,但是只有定性分析充滿(mǎn)度和平緩度,沒(méi)有建立定量的評(píng)價(jià)指標(biāo)。

1)充滿(mǎn)度。

后坐阻力FR在后坐長(zhǎng)度λ內(nèi)所做的總負(fù)功等于炮膛合力對(duì)后坐部分的總功,因此理想上的后坐阻力為常數(shù)FRC。如圖6所示,充滿(mǎn)度指陰影面積充滿(mǎn)FRC-s直線作功面積的程度。

充滿(mǎn)度α計(jì)算公式為

α=Sy/S

(9)

式中:Sy為設(shè)計(jì)的FR-s曲線和s軸所圍成的圖形與

FRC-s直線和s軸所圍圖形的交集所對(duì)應(yīng)的面積,即陰影面積;S為FRC-s直線形成的矩形面積。

充滿(mǎn)度α∈(0,1),α越接近1代表FR-s曲線的充滿(mǎn)度越好,也越接近理想后坐阻力。

2)平緩度。

后坐阻力的變化越平緩越好,它有利于全炮受力和減小振動(dòng),因此平緩度是評(píng)價(jià)后坐阻力的另一重要指標(biāo)。在FR-s曲線中,平緩度指后坐阻力隨后坐行程變化的快慢程度,它是對(duì)整條曲線全局的描述;而導(dǎo)數(shù)正是反映函數(shù)因變量隨自變量變化的快慢程度,但導(dǎo)數(shù)是函數(shù)的局部性質(zhì)。因此本文采用導(dǎo)數(shù)的平均值來(lái)描述FR-s曲線的平緩度β:

(10)

式中:q為在整個(gè)后坐行程上等距且連續(xù)的取點(diǎn)個(gè)數(shù);kw為上述第w個(gè)取點(diǎn)在FR-s曲線的導(dǎo)數(shù)。β越小,FR-s曲線的平緩度越好。

2.5 結(jié)果及分析

節(jié)制桿外徑公差等級(jí)保持IT10不變,采用1.2.4節(jié)算法求解計(jì)及流液孔尺寸不確定的后坐阻力優(yōu)化模型,優(yōu)化前后的不確定性變量取公稱(chēng)值對(duì)應(yīng)的FR-s曲線,如圖7所示。繼續(xù)取節(jié)制桿外徑公差等級(jí)為IT8、IT9、IT11和IT12進(jìn)行數(shù)值計(jì)算,結(jié)果圓整后計(jì)算得到后坐阻力與后坐長(zhǎng),如表4所示。

表4 節(jié)制桿外徑取不同公差等級(jí)的優(yōu)化結(jié)果

結(jié)果表明,優(yōu)化后的后坐阻力峰值區(qū)間較優(yōu)化前有較大幅度降低,后坐阻力曲線更加平緩,充滿(mǎn)度更好,后坐長(zhǎng)基本不變,而后坐阻力得到了優(yōu)化,還實(shí)現(xiàn)了節(jié)制桿外徑尺寸設(shè)計(jì)與公差設(shè)計(jì)的同步優(yōu)化,使考慮設(shè)計(jì)方案的工藝性變得可行。

節(jié)制桿屬于多錐度細(xì)長(zhǎng)桿,加工過(guò)程中由于它剛性差和多錐度變截面的問(wèn)題,導(dǎo)致加工高質(zhì)量的節(jié)制桿不易,因此節(jié)制桿的工藝性需重點(diǎn)考慮。5個(gè)優(yōu)化方案中,從工藝性角度考慮,以公差等級(jí)越大工藝性越好為準(zhǔn)則;從后坐阻力的角度考慮,以后坐阻力峰值區(qū)間越小后坐阻力越好為準(zhǔn)則。當(dāng)設(shè)計(jì)者較為關(guān)注后坐阻力時(shí),可選擇節(jié)制桿公差等級(jí)為IT8的優(yōu)化設(shè)計(jì)方案;當(dāng)設(shè)計(jì)者較為關(guān)注加工工藝性時(shí),可選擇公差等級(jí)為IT11或IT12的優(yōu)化設(shè)計(jì)方案;設(shè)計(jì)者需綜合考慮工藝性和后坐阻力時(shí),可選擇節(jié)制桿公差等級(jí)為IT9或IT10的優(yōu)化設(shè)計(jì)方案。

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

本文提出了一種三參數(shù)區(qū)間的不確定性?xún)?yōu)化模型及算法,定義了三參數(shù)區(qū)間來(lái)描述結(jié)構(gòu)尺寸的不確定性,研究了三參數(shù)區(qū)間的轉(zhuǎn)化,采用NSGA-Ⅱ和試驗(yàn)設(shè)計(jì)方法相結(jié)合的算法,將兩層嵌套優(yōu)化解耦為單層優(yōu)化來(lái)求解。懸臂梁算例結(jié)果表明,該方法將結(jié)構(gòu)尺寸優(yōu)化設(shè)計(jì)與公差設(shè)計(jì)相結(jié)合,得到了關(guān)鍵尺寸及公差帶,并能得到穩(wěn)健性更好的設(shè)計(jì)方案。

利用上述提出的方法進(jìn)行了計(jì)及流液孔尺寸不確定的后坐阻力優(yōu)化,并對(duì)后坐阻力曲線的2個(gè)定性評(píng)價(jià)指標(biāo)充滿(mǎn)度和平緩度進(jìn)行了定量化表征。結(jié)果表明,優(yōu)化后的后坐阻力峰值區(qū)間較優(yōu)化前有較大幅度降低,后坐阻力曲線更加平緩,充滿(mǎn)度更好。實(shí)現(xiàn)了節(jié)制桿外徑尺寸設(shè)計(jì)與公差設(shè)計(jì)的同步優(yōu)化,得到了兼顧制造和性能的設(shè)計(jì)方案。

[1] 高樹(shù)滋,陳運(yùn)生,張?jiān)铝?等.火炮反后坐裝置設(shè)計(jì)[M].北京:兵器工業(yè)出版社,1995. GAO Shu-zi,CHEN Yun-sheng,ZHANG Yue-lin,et al.Design of recoil system for gun[M].Beijing:Ordnance Industry Press,1995.(in Chinese)

[2] 宗士增,錢(qián)林方,徐亞棟.火炮反后坐裝置動(dòng)力學(xué)耦合分析與優(yōu)化[J].兵工學(xué)報(bào),2007,28(3):272-275. ZONG Shi-zeng,QIAN Lin-fang,XU Ya-dong.Dynamic coupling analysis and optimization of gun recoil mechamism[J].Acta Armamentarii,2007,28(3):272-275.(in Chinese)

[3] 周成,顧克秋,邵躍林.火炮制退機(jī)節(jié)制桿優(yōu)化設(shè)計(jì)[J].火炮發(fā)射與控制學(xué)報(bào),2012(4):68-70. ZHOU Cheng,GU Ke-qiu,SHAO Yue-lin.Optimal design of throttling bar in gun recoil mechanism[J].Journal of Gun Launch & Control,2012(4):68-70.(in Chinese)

[4] 姜潮.基于區(qū)間的不確定性?xún)?yōu)化理論與算法[D].長(zhǎng)沙:湖南大學(xué),2008:15-18. JIANG Chao.Theories and algorithms of uncertain optimization based on interval[D].Changsha:Hunan University,2008:15-18.(in Chinese)

[5] 李榮,楊國(guó)來(lái),劉寧,等.火炮制退機(jī)關(guān)鍵結(jié)構(gòu)尺寸及公差優(yōu)化研究[J].彈道學(xué)報(bào),2016,28(4):80-85. LI Rong,YANG Guo-lai,LIU Ning,et al.Research on optimization for key structure dimension of gun recoil mechanism considering tolerance[J].Journal of Ballistics,2016,28(4):80-85.(in Chinese)

[6] JIANG C,XIE H C,ZHANG Z G,et al.A new interval optimization method considering tolerance design[J].Engineering Optimization,2015,47(12):1 637-1 650.

[7] 謝慧超,姜潮,張智罡,等.基于區(qū)間分析的汽車(chē)平順性?xún)?yōu)化[J].汽車(chē)工程,2014,36(9):1 127-1 131. XIE Hui-chao,JIANG Chao,ZHANG Zhi-gang,et al.Vehicle ride comfort optimization based on interval analysis[J].Automotive Engineering,2014,36(9):1 127-1 131.(in Chinese)

[8] 徐澤水,達(dá)慶利.區(qū)間數(shù)排序的可能度法及其應(yīng)用[J].系統(tǒng)工程學(xué)報(bào),2003,18(1):67-70. XU Ze-shui,DA Qing-li.Possibility degree method for ranking interval numbers and its application[J].Journal of Systems Engineering,2003,18(1):67-70.(in Chinese)

[9] 王群.基于對(duì)稱(chēng)拉丁超立方設(shè)計(jì)的多目標(biāo)進(jìn)化算法[D].西安:西安電子科技大學(xué),2011:29-43. WANG Qun.Multiobjective evolutionary algorithm based on symmetric Latin hypercube designs[D].Xi’an:Xidian University,2011:29-43.(in Chinese)

Optimization on Interval Uncertainty of Artillery Recoil Resistance

LI Rong,YANG Guo-lai,SUN Quan-zhao,XIAO Hui

(School of Mechanical Engineering,Nanjing University of Science and Technology,Nanjing 210094,China)

In order to reduce the force of artillery carriage and obtain the better recoil resistance,the optimization on interval uncertainty of recoil fluid-hole were conducted.The three-parameter interval was defined combining with tolerance,and then the uncertainty optimization model and algorithm based on three-parameter interval were established.Taking the dimension of fluid hole as the uncertainty variable,the interval optimization design of the shape size of control rod was carried out by using the method mentioned above.The outer diameter of the control rod with different tolerance-grades was optimized.Research results show that the peak interval of recoil resistance is reduced after the optimization,and the recoil resistance curve is more gentler with a better fullness.The synchronous optimization of the tolerance design and the design of the outer diameter of control rod was realized.The optimum design-scheme of the outer diameter of the control rod considering the workmanship and the recoil resistance was obtained.

artillery;recoil resistance;structure design;interval optimization;uncertainty dimension

2016-12-16

國(guó)防973計(jì)劃項(xiàng)目(1503613249);國(guó)家自然科學(xué)基金項(xiàng)目(11572158)

李榮(1991- ),男,碩士研究生,研究方向?yàn)榛鹋诮Y(jié)構(gòu)優(yōu)化與優(yōu)選。E-mail:lirongwork@qq.com。

楊國(guó)來(lái)(1968- ),男,教授,博士生導(dǎo)師,研究方向?yàn)榛鹋诂F(xiàn)代設(shè)計(jì)理論與方法。E-mail:yyanggl@mail.njust.edu.cn。

TJ303

A

1004-499X(2017)02-0078-07