本發(fā)明涉及一種發(fā)射藥壓伸參數(shù)的動態(tài)計算方法，適用于發(fā)射藥產(chǎn)品制造過程中壓伸工序工藝參數(shù)的實時修正計算。
背景技術(shù)：
：發(fā)射藥產(chǎn)品制備的工序包含吸收、驅(qū)水、造粒、膠化、壓伸、烘干等，其中壓伸工序是發(fā)射藥的成型工序，發(fā)射藥的尺寸精度、密實度、表面光澤性等均由壓伸過程的質(zhì)量決定，而壓伸工序的核心就是壓伸工藝參數(shù)的設計，在模具一定的條件下，通過合適的工藝參數(shù)才能獲得理想的發(fā)射藥產(chǎn)品。現(xiàn)有的發(fā)射藥壓伸工藝參數(shù)主要依靠經(jīng)驗數(shù)據(jù)，人工推斷獲得，如制備三基發(fā)射藥，則需要根據(jù)以往的試驗數(shù)據(jù)進行調(diào)整，若沒有試驗數(shù)據(jù)，則需要從最小的壓伸壓力開始，逐步上調(diào)壓力，進行壓伸，直到獲得認為合適的壓伸壓力為止?，F(xiàn)有的發(fā)射藥壓伸工藝參數(shù)的計算方法存在如下問題：(一)計算過程需要依靠已有的試驗數(shù)據(jù)，通過人工經(jīng)驗進行修正，經(jīng)驗豐富的人員與沒有經(jīng)驗的人員所得結(jié)果差異較大；(二)壓伸工藝參數(shù)的計算值不能進行實時修正，若壓力參數(shù)不合適，無法進行快速修正計算；(三)在對壓伸參數(shù)進行修正時，需要人工手動修正，修正的數(shù)值大小，全憑經(jīng)驗進行，無法考慮發(fā)射藥的配方特性及當前工藝條件的綜合影響。技術(shù)實現(xiàn)要素：為了解決現(xiàn)有技術(shù)存在的問題，本發(fā)明提供了一種能夠?qū)崿F(xiàn)不依靠經(jīng)驗數(shù)據(jù)及人工經(jīng)驗，在發(fā)射藥壓伸過程中實時利用當前獲得的壓伸數(shù)據(jù)進行動態(tài)計算，給出較合適的下一時間點的壓伸工藝參數(shù)的動態(tài)計算方法。計算步驟如下：(1)建立發(fā)射藥壓伸工藝參數(shù)與出藥速度、外徑等的數(shù)據(jù)庫；(2)輸入發(fā)射藥配方及目標出藥速度，查詢數(shù)據(jù)庫并采用線性插值法進行壓伸工藝參數(shù)與目標出藥速度的匹配計算；(3)將目標出藥速度vm與對應的發(fā)射藥壓伸壓力pc輸入到壓伸機；(4)實時測量出藥速度與對應的壓伸壓力值，獲得第一個時間點的p1、v1值；(5)對比第一個時間點p1、v1與vm、pc的數(shù)據(jù)，有dp＝p1×(vm-v1)/v1×y，其中dp為壓力調(diào)整幅度，p1為第一時間點的壓伸壓力，v1為第一時間點的出藥速度，y為非線性因子，由不同發(fā)射藥的配方及藥料性質(zhì)決定，一般取0.5-0.8；(6)輸入目標出藥速度vm和壓伸壓力pc1到壓伸機，pc1＝pc+dp，若dp大于零則壓伸壓力增大，若dp小于零則壓伸壓力減小；(7)實時測量出藥速度v與對應的壓伸壓力值p，獲得第二個時間點的p2、v2值；(8)對比第二個時間點p2、v2與vm、pc的數(shù)據(jù)，有dp＝p2×(vm-v2)/v2×y；(9)輸入目標出藥速度vm和壓伸壓力pc2到壓伸機，pc2＝pc1+dp；(10)對于以后的每一點執(zhí)行7、8、9步驟；(11)在經(jīng)歷n步以后(n≥1)，獲得滿足目標出藥速度vm的壓伸工藝參數(shù)。因此，采用數(shù)據(jù)庫及插值算法，前后時間點數(shù)據(jù)的修正算法，不同藥料特性的非線性因子，實時數(shù)據(jù)點的循環(huán)迭代，能夠?qū)崿F(xiàn)不依靠經(jīng)驗數(shù)據(jù)及人工經(jīng)驗，在發(fā)射藥壓伸過程中實時利用當前獲得的壓伸數(shù)據(jù)進行動態(tài)計算，給出較合適的下一時間點的壓伸工藝參數(shù)的動態(tài)計算方法。本發(fā)明的有益效果主要體現(xiàn)在以下幾個方面。(一)采用數(shù)據(jù)庫及插值算法，能夠根據(jù)已有相關(guān)數(shù)據(jù)獲得所需數(shù)據(jù)的初始值。與現(xiàn)有技術(shù)相比，對于任意發(fā)射藥的配方，均能給出初始的壓伸工藝參數(shù)及出藥速度的匹配數(shù)值。(二)采用前后時間點數(shù)據(jù)的修正算法，能夠利用實時測量的真實試驗數(shù)據(jù)進行修正，逐步逼近合適的壓伸工藝參數(shù)值。與現(xiàn)有技術(shù)相比，能夠充分利用本次試驗的真實數(shù)據(jù)進行修正，精度更高。(三)采用不同藥料特性的非線性因子，能夠針對不同發(fā)射藥進行修正計算，保證能夠?qū)Σ煌愋桶l(fā)射藥進行較好的匹配計算。與現(xiàn)有技術(shù)相比，計算的適用面更寬，計算的步數(shù)較少。(四)采用實時數(shù)據(jù)點的循環(huán)迭代方法，保證初始輸入值在一定的范圍內(nèi)都能夠獲得較合適的工藝參數(shù)。與現(xiàn)有技術(shù)相比，即使初始輸入值與目標值的偏差較大，也能夠保證獲得較合適壓伸工藝參數(shù)。附圖說明圖1是發(fā)射藥壓伸工藝參數(shù)動態(tài)計算步驟示意圖。具體實施方式下面對本發(fā)明的具體實施方式作進一步的詳細說明。根據(jù)圖1，以三基發(fā)射藥為例，該發(fā)射藥的壓伸工藝參數(shù)動態(tài)計算步驟如下。(1)建立發(fā)射藥壓伸工藝參數(shù)與出藥速度、外徑等的數(shù)據(jù)庫；(2)輸入發(fā)射藥配方及目標出藥速度25mm/s，發(fā)射藥配方如下nc(％)ng(％)rdx(％)ngu(％)c2(％)tio2(％)32.3±2.026.2±1.227.6±1.512.0±1.01.4±0.30.5±0.2查詢數(shù)據(jù)庫后采用線性插值法計算結(jié)果如下：pc＝10mpa，vm＝25mm/s(3)將目標出藥速度vm與對應的發(fā)射藥壓伸壓力pc輸入到壓伸機；(4)實時測量出藥速度v與對應的壓伸壓力值p值如下：p＝10mpa，v＝15mm/s；(5)對比第一個時間點p、v與vm、pc的數(shù)據(jù)，有dp＝p1×(vm-v1)/v1×ydp為壓力調(diào)整幅度，p1為第一時間點的壓伸壓力，v1為第一時間點的出藥速度，y為非線性因子，由不同發(fā)射藥的配方及藥料性質(zhì)決定，一般取0.5-0.8。該配方發(fā)射藥取y＝0.5，計算得dp＝3.33mpa。(6)采用如下公式計算：pc1＝pc+dp，pc1＝13.33mpa輸入目標出藥速度vm和壓伸壓力pc1到壓伸機(7)實時測量出藥速度v與對應的壓伸壓力值p，有p＝13.33mpa，v＝18mm/s；(8)對比第二個時間點p、v與vm、pc的數(shù)據(jù)，有dp＝p2×(vm-v2)/v2×y經(jīng)計算得dp＝2.59mpa(9)采用如下公式計算：pc2＝pc1+dp，pc2＝15.92mpa輸入目標出藥速度vm和壓伸壓力pc2到壓伸機(10)對于以后的每一點執(zhí)行7、8、9步驟；(11)在經(jīng)歷n步以后(n≥1)，pcn＝18mpa，vm＝25mm/s此時即為滿足目標出藥速度vm的壓伸工藝參數(shù)。對于該配方的發(fā)射藥，采用以上計算步驟，經(jīng)過n步以后(n≥1)，計算得到滿足目標出藥速度vm的壓伸工藝參數(shù)，pcn＝18mpa，vm＝25mm/s。因此，本發(fā)明的一種壓伸工藝參數(shù)的動態(tài)計算方法，能夠?qū)崿F(xiàn)不依靠經(jīng)驗數(shù)據(jù)及人工經(jīng)驗，在發(fā)射藥壓伸過程中實時利用當前獲得的壓伸數(shù)據(jù)進行動態(tài)計算，并最后獲得滿足目標出藥速度的壓伸工藝參數(shù)。當前第1頁12