CN106290233B - 双基推进剂装药安定期快速预估方法 - Google Patents
双基推进剂装药安定期快速预估方法 Download PDFInfo
- Publication number
- CN106290233B CN106290233B CN201610614757.9A CN201610614757A CN106290233B CN 106290233 B CN106290233 B CN 106290233B CN 201610614757 A CN201610614757 A CN 201610614757A CN 106290233 B CN106290233 B CN 106290233B
- Authority
- CN
- China
- Prior art keywords
- effective content
- double base
- sample
- temperature
- base propellant
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000003380 propellant Substances 0.000 title claims abstract description 96
- 239000000843 powder Substances 0.000 title claims abstract description 87
- 238000000034 method Methods 0.000 title claims abstract description 66
- PZIMIYVOZBTARW-UHFFFAOYSA-N centralite Chemical compound C=1C=CC=CC=1N(CC)C(=O)N(CC)C1=CC=CC=C1 PZIMIYVOZBTARW-UHFFFAOYSA-N 0.000 claims abstract description 104
- 230000032683 aging Effects 0.000 claims abstract description 60
- 230000001133 acceleration Effects 0.000 claims abstract description 42
- 238000003860 storage Methods 0.000 claims abstract description 24
- 238000005086 pumping Methods 0.000 claims abstract description 16
- 238000001228 spectrum Methods 0.000 claims abstract description 12
- 238000001514 detection method Methods 0.000 claims abstract description 11
- 238000003333 near-infrared imaging Methods 0.000 claims abstract description 8
- 239000000126 substance Substances 0.000 claims description 39
- 239000010410 layer Substances 0.000 claims description 27
- 239000003795 chemical substances by application Substances 0.000 claims description 26
- 238000004448 titration Methods 0.000 claims description 18
- 238000002329 infrared spectrum Methods 0.000 claims description 16
- 230000006641 stabilisation Effects 0.000 claims description 11
- 238000011105 stabilization Methods 0.000 claims description 11
- 238000012360 testing method Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000009795 derivation Methods 0.000 claims description 7
- 238000004128 high performance liquid chromatography Methods 0.000 claims description 6
- 239000011229 interlayer Substances 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 5
- 238000013213 extrapolation Methods 0.000 claims description 4
- 238000007689 inspection Methods 0.000 claims description 4
- 238000012417 linear regression Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 238000010561 standard procedure Methods 0.000 claims description 4
- 238000004811 liquid chromatography Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 239000003814 drug Substances 0.000 claims description 2
- 239000003381 stabilizer Substances 0.000 claims description 2
- 238000003556 assay Methods 0.000 claims 1
- 239000000523 sample Substances 0.000 description 69
- 230000015556 catabolic process Effects 0.000 description 20
- 238000006731 degradation reaction Methods 0.000 description 20
- 238000003483 aging Methods 0.000 description 19
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 12
- RTZKZFJDLAIYFH-UHFFFAOYSA-N Diethyl ether Chemical compound CCOCC RTZKZFJDLAIYFH-UHFFFAOYSA-N 0.000 description 6
- 238000005070 sampling Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 239000002360 explosive Substances 0.000 description 5
- 230000003595 spectral effect Effects 0.000 description 5
- 238000004880 explosion Methods 0.000 description 4
- 239000002351 wastewater Substances 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 238000002512 chemotherapy Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 239000004615 ingredient Substances 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000002360 preparation method Methods 0.000 description 3
- 238000004611 spectroscopical analysis Methods 0.000 description 3
- CPELXLSAUQHCOX-UHFFFAOYSA-M Bromide Chemical compound [Br-] CPELXLSAUQHCOX-UHFFFAOYSA-M 0.000 description 2
- SNIOPGDIGTZGOP-UHFFFAOYSA-N Nitroglycerin Chemical compound [O-][N+](=O)OCC(O[N+]([O-])=O)CO[N+]([O-])=O SNIOPGDIGTZGOP-UHFFFAOYSA-N 0.000 description 2
- 238000005844 autocatalytic reaction Methods 0.000 description 2
- 238000009835 boiling Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 229960003711 glyceryl trinitrate Drugs 0.000 description 2
- 239000007791 liquid phase Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 229920001220 nitrocellulos Polymers 0.000 description 2
- 238000009659 non-destructive testing Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- WKBOTKDWSSQWDR-UHFFFAOYSA-N Bromine atom Chemical compound [Br] WKBOTKDWSSQWDR-UHFFFAOYSA-N 0.000 description 1
- 229910004679 ONO2 Inorganic materials 0.000 description 1
- 230000031709 bromination Effects 0.000 description 1
- 238000005893 bromination reaction Methods 0.000 description 1
- GDTBXPJZTBHREO-UHFFFAOYSA-N bromine Substances BrBr GDTBXPJZTBHREO-UHFFFAOYSA-N 0.000 description 1
- 229910052794 bromium Inorganic materials 0.000 description 1
- 239000013068 control sample Substances 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000011068 loading method Methods 0.000 description 1
- 150000002823 nitrates Chemical class 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 238000010238 partial least squares regression Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000002203 pretreatment Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000036632 reaction speed Effects 0.000 description 1
- 238000010992 reflux Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/359—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N31/00—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods
- G01N31/16—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods using titration
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Molecular Biology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
一种双基推进剂装药安定期快速预估方法,利用热加速老化法制备Ⅱ号中定剂有效含量近红外定量模型建模样品,采用近红外漫反射光谱仪采集光谱并建立Ⅱ号中定剂有效含量近红外定量模型;以可靠温度系数r0建立安定期预估方程;近红外方法跟踪检测单一温度装药热加速老化过程中Ⅱ号中定剂有效含量,获得有效含量到达临界判据消耗的临界时间,由安定期预估方程外推工况及贮存环境温度下的双基推进剂装药安定期。本方法只需要近红外方法无损跟踪一个温度下的一个装药样品不同老化时间有效含量便可获得安定期,样品量减少96%,安全、低成本。
Description
技术领域
本发明属于火炸药安全性评估技术领域,涉及一种火炸药工况、贮存过程安定期快速预估方法。特别是一种利用近红外漫反射光谱法检测双基推进剂(ST)装药单一温度热加速老化样品Ⅱ号中定剂有效含量,获得Ⅱ号中定剂有效含量到达临界判据消耗的时间,根据安定期预估方程快速预估工况及贮存环境温度下双基推进剂装药安定期的方法。
背景技术
火炸药化学安定性是火炸药在工况、贮存条件下保持其化学性质变化不超过允许范围以致引发安全事故的能力。
硝化棉、硝化甘油是双基推进剂中的能量成分,硝化棉和硝化甘油分子结构中都含有-C-ONO2结构,在受热或贮存过程中双基推进剂发生缓慢分解,O-NO2键断裂释放能量,同时产生具有自催化作用的氮氧化物,氮氧化物进一步催化能量成分分解,持续缓慢放热并产生热积累,降低化学安定性,导致燃烧、爆炸,影响生产及贮存安全性。
提高双基推进剂化学安定性有效的方法是在双基推进剂中加入Ⅱ号中定剂作为安定剂,它能吸收双基推进剂分解放出的氮氧化物,从而抑制氮氧化物对双基推进剂分解的自催化作用,同时生成仍然具有安定作用Ⅱ号中定剂衍生物,延缓能量成分分解,提高双基推进剂化学安定性。Ⅱ号中定剂及其衍生物总含量称为Ⅱ号中定剂有效含量,工况及贮存条件下当Ⅱ号中定剂有效含量下降至临界判据时经历的时间为双基推进剂安定期,在安定期内双基推进剂化学安定性良好。
双基推进剂采用热加速老化试验跟踪Ⅱ号中定剂有效含量的方法获得安定期,常用化学滴定法(溴与安定剂反应生成溴化物的化学方法)或高效液相色谱法测定Ⅱ号中定剂有效含量,采用乙醚回流24小时甚至更长的时间提取粉碎试样Ⅱ号中定剂有效含量,低沸点溶剂乙醚易燃易爆,及通过机械粉碎双基推进剂老化试样制备乙醚提取用试样时安全性差,操作繁琐、耗时且化学滴定法产生的溴化物带来污染。通过热加速老化试验获得的不同温度下老化试样Ⅱ号中定剂有效含量变化规律,以阿累尼乌斯方程获得安定期时,由于热加速老化试验至少采用4个温度点(65℃、75℃、85℃、95℃),每个温度点取样次数至少为6次,需要的老化样品量大,试验持续时间6个月以上,极为耗时,且长时间高温热加速老化危险性高,难以满足新配方研发及工艺过程中化学安定性设计和评价需求。
对于工况及贮存过程具有一定尺寸的双基推进剂装药,安定期预估时,通过热加速老化试验受热过程会产生热积累,一方面导致温场不均匀分布,引起安定剂含量分布不均匀,另一方面双基推进剂装药多温度点热加速老化试验极易因热积累发生燃烧及爆炸事故,因此,难以通过双基推进剂装药多个温度点老化试验并采用化学滴定法或高效液相色谱法获得安定剂含量判定安定期,难以满足双基推进剂装药贮存过程中化学安定性评价的需求。
因此,双基推进剂装药需要采用新的方法获得安定期。
近红外光谱技术是基于有机物分子中X-H(X-C、N、O)的倍频、合频信息对物质进行定性、定量分析的光谱技术,具有便捷、无损、绿色的优点。双基推进剂装药安定期快速预估方法,能安全、省时、低成本,环保的获得Ⅱ号中定剂安定的双基推进剂装药安定期。
发明内容
针对上述现有双基推进剂装药安定期技术存在的缺陷或不足,本发明的目的在于,提供一种双基推进剂装药安定期快速预估方法,该方法能安全、省时、低成本,环保的获得双基推进剂装药安定期。
为了实现上述任务,本发明采取如下的技术方案:
一种双基推进剂装药安定期快速预估方法,其特征在于,该方法利用热加速老化法制备双基推进剂装药Ⅱ号中定剂有效含量近红外定量模型建模样品,采用化学滴定法或高效液相色谱法获取Ⅱ号中定剂有效含量化学值,利用近红外漫反射光谱仪采集光谱,通过化学计量学建立Ⅱ号中定剂有效含量近红外定量模型;进行双基推进剂装药热加速老化试验,通过Ⅱ号中定剂有效含量变化规律获得可靠温度系数r0,建立装药安定期预估方程;然后开展单一温度双基推进剂装药热加速老化试验,近红外方法无损跟踪检测老化过程中双基推进剂装药Ⅱ号中定剂有效含量,获得Ⅱ号中定剂有效含量到达临界判据消耗的临界时间,根据安定期预估方程及可靠温度系数外推工况及贮存环境温度下双基推进剂装药的安定期。具体按下列步骤进行:
(一)热加速老化法制备Ⅱ号中定剂有效含量近红外建模样品
双基推进剂药饼三明治式夹心结构试样,由上层药片、夹心层、下层药片依次叠加组成药饼,药饼装入带有毛细排气管的磨口减量瓶中,于单一温度下Tn下热加速老化试验,间隔时间取夹心层,获得不同Ⅱ号中定剂有效含量的双基推进剂装药样品,分为校正集样品和验证集样品。
(二)Ⅱ号中定剂有效含量近红外定量模型
1、采用化学滴定法或高效液相色谱法获取双基推进剂装药Ⅱ号中定剂有效含量化学值,确定含量范围;
2、采用近红外光谱仪对校正集样品和验证集样品采集近红外光谱。采样方式为漫反射,采样波段为700nm-2500nm,优化光谱扫描参数,确定最佳分辨率、扫描次数及样品重复测量次数。采集到的近红外光谱传输至计算机。
3、选择特征谱带,采用归一化法、基线平滑和一阶求导二阶求导、多元散射校正或其组合形式对校正集样品光谱进行预处理,继而采用偏最小二乘法,即PLS法进行回归拟合,建立Ⅱ号中定剂有效含量近红外定量模型。用验证集样品对所建模型进行外部验证,依据模型的内部检验和外部验证相结合的评价体系,逐步优化模型。
根据建立的Ⅱ号中定剂有效含量近红外定量模型,采用近红外光谱仪对选定的样品进行Ⅱ号中定剂有效含量测定,获得近红外预测值。化学值与近红外预测值的最大误差应小于标准方法重复性误差。
(三)安定期预估方程
硝酸酯类含能材料分解引起双基推进剂Ⅱ号中定剂有效含量变化的化学反应,服从阿累尼乌斯方程,表示k-T关系的较准确的经验式为:
定义反应速度的温度系数r为温度每上升10℃反应速度常数的变化倍率:
r=kT+10/kT 式(2)
式中:r表示温差为10℃的反应速率温度系数;
kT表示T℃时的速率常数;
kT+10表示(T+10)℃时的速率常数。
由于在T1和T2两个温度下,分别进行同一反应,aA+bB+…→lL+…,都从同一初始浓度进行到相同的转化率,所需时间分别为t1和t2,k1与k2分别为T1和T2时的速率常数,则有
t2/t1=k1/k2 式(3)
设:Tn-Tn-1=10℃
由式(1)(2)可得:
式中:r表示温差为10℃的反应速率温度系数;
t0表示工况或贮存环境温度安定期,d;
tn表示高温热加速老化时间,d;
Tn表示高温热加速老化试验温度,℃;
T0表示工况或贮存环境温度,℃。
(四)可靠温度系数r0
对公式(4)两边取对数,整理后得到:
Tn=A+Blgtn 式(5)
其中,A=T0+10×lgt0/lgr、B=-10/lgr
温度系数r=10-10/B,采用近红外方法跟踪测试样品不同温度Tni下热加速老化试验不同老化时间的Ⅱ号中定剂有效含量,获得不同温度Ⅱ号中定剂有效含量随时间变化关系曲线,以Ⅱ号中定剂有效含量消耗50%作为临界判据,获得不同温度Tni对应的时间tni,经代入(5)式进行一元线性回归,采用最小二乘法获得系数B,进而计算得到温度系数r。
采用上述方法获得多种双基推进剂的温度系数。
双基推进剂安定期预估方程中,可靠温度系数保证了安定期内化学安定性良好。可靠温度系数取值偏小,获得安定期短,可靠性高,但取值过小,安定期过于保守,提前销毁造成浪费;反之,可靠温度系数的取值过大,获得安定期长,超过安全贮存寿命时,引起热爆炸。因此,以多种典型双基推进剂温度系数的均值与标准偏差差值的0.9倍作为可靠温度系数r0。
(五)热加速老化试验与安定期
将双基推进剂药柱沿径向等分为两段,将两段药柱沿轴向上下叠加制作成两段式药柱(附图3),获得装药两段式老化试样,于单一温度下Tn下热加速老化试验,选取下段药柱叠合面的中心区域及药柱表面进行近红外光谱扫描,获得药柱在不同老化时间、不同部位的Ⅱ号中定剂有效含量,以Ⅱ号中定剂有效含量消耗最快的区域作为关重部位,以关重部位Ⅱ号中定剂有效含量消耗50%作为临界判据,获得温度Tn下安定剂消耗的时间tn,依公式(4)快速预估工况及贮存环境温度下双基推进剂装药的安定期。
本发明的双基推进剂装药安定期快速预估方法,带来的有益技术效果体现在以下几个方面:
1、双基推进剂Ⅱ号中定剂发挥安定作用的过程中,不仅Ⅱ号中定剂自身具有安定效果,Ⅱ号中定剂衍生物仍然具有安定效果,是Ⅱ号中定剂有效含量的重要组成部分,因此,双基推进剂装药老化过程中的安定剂含量为Ⅱ号中定剂有效含量。由于Ⅱ号中定剂衍生物种类多,结构复杂,难以获得,通过热加速老化试验获得近红外安定剂含量建模样品的方法,克服了多种安定剂衍生物制备及配比问题,回避了复杂危险性高的双基推进剂制造过程,便捷、安全、低成本、环保。
2、热加速老化法制备的近红外Ⅱ号中定剂有效含量建模样品通过控制取样时间获得不同Ⅱ号中定剂有效含量的样品,达到在变化范围内分布均匀,解决了工艺制造带来的低含量样品分布不均匀的缺陷。
3、化学滴定法和高效液相色谱法为安定剂含量的有损检测方法,只能对双基推进剂材料进行安定剂含量检测,无法对装药安定剂含量进行检测,不能评定装药瞬态化学安定性。近红外光线具有很强的穿透能力,在检测样品时,直接对样品做近红外光谱的扫描,这种无损检测技术,无需破坏装药便可以直接测定装药各部位安定剂含量,评价双基推进剂装药瞬态化学安定性。此外,近红外安定剂含量建模后进行Ⅱ号中定剂有效含量检测时,可在2~5min内各给出检测结果,回避了化学滴定法和高效液相色谱法中Ⅱ号中定剂分离时样品需要在低沸点乙醚溶剂中回流提取24h的易燃工序,解决了耗时长、效率低及污染等问题,操作简便、低耗、环保、安全、快速。
4、以往获得双基推进剂的温度系数时,采用化学滴定法进行破坏性检测,跟踪不同温度不同老化时间Ⅱ号中定剂有效含量,以获得至少4个温度下6个老化时间与Ⅱ号中定剂有效含量的对应关系,进而拟合得到温度系数,可见,采用化学滴定法获得温度系数,至少需要25个以上样品测定Ⅱ号中定剂有效含量。近红外有效安定剂含量建模后,跟踪4个温度热加速老化试样安定剂含量变化时,只需4个样品进行非破坏性检测,回避了化学滴定法需要多次取出老化试样进行破坏性检测,样品量减少84%,同时回避了多次(每个温度点至少6次)取样进行机械粉碎的危险操作过程,更为快速、便捷,实现本质安全的同时大幅降低了成本。
5、通过热加速老化试验获得多种典型双基推进剂的温度系数,取温度系数均值与标准偏差差值的0.9作为可靠温度系数r0。近红外跟踪双基推进剂装药不同老化时间Ⅱ号中定剂有效含量,通过对比多种典型的温度系数,以Ⅱ号中定剂有效含量消耗50%作为临界判据,采用预估方程获得装药安定期,建立了双基推进剂装药安定期快速预估方法。这种采用单一高温热加速老化并对单一样品进行Ⅱ号中定剂有效含量近红外无损检测的方法,与GJB770B-2005方法506.1预估贮存寿命热加速老化法相比,无需进行装药样品多个温度热加速老化并且多次取样进行Ⅱ号中定剂有效含量滴定的破坏性检测,只需要近红外方法跟踪一个温度下的一个装药样品不同老化时间Ⅱ号中定剂有效含量便可获得安定期,更为快速、便捷,样品量减少96%,实现双基推进剂装药安定期快速预估。
附图说明
图1双基推进剂药饼三明治式夹心结构试样.1-上层药片,2-夹心层药片,3-下层药片。
图2毛细排气管磨口减量瓶.4-毛细排气管,5-磨口减量瓶盖,6-磨口减量瓶体。
图3双基推进剂两段式药柱结构试样.7-上层药柱8-中心区域9-下层药柱叠合面10-下层药柱。
图4双基推进剂ST-1校正集样品光谱图。
图5ST-1不同温度Ⅱ号中定剂有效含量随时间变化关系曲线。
具体实施方式
本发明的双基推进剂装药安定期快速预估方法,其特征在于,该方法利用热加速老化法制备Ⅱ号中定剂有效含量近红外定量模型建模样品,依据化学滴定法或高效液相色谱法获取Ⅱ号中定剂有效含量化学值,利用近红外漫反射光谱仪采集光谱,通过化学计量学建立Ⅱ号中定剂有效含量近红外定量模型;以Ⅱ号中定剂有效含量消耗50%作为安定期的临界判据,开展双基推进剂装药热加速老化试验,通过Ⅱ号中定剂有效含量变化规律获得温度系数,以多种典型双基推进剂温度系数的均值与标准偏差差值的0.9倍作为可靠温度系数r0,获得安定期预估方程;然后开展单一温度热加速老化试验,近红外方法跟踪检测老化过程中样品Ⅱ号中定剂有效含量,获得Ⅱ号中定剂有效含量到达临界判据消耗的时间,根据安定期预估方程及可靠温度系数外推工况及贮存环境温度下双基推进剂装药的安定期。具体按下列步骤进行:
(一)热加速老化法制备Ⅱ号中定剂有效含量近红外建模样品
双基推进剂药饼三明治式夹心结构试样,由上层药片1、夹心层药片2、下层药片3依次叠加组成药饼(附图1),药饼装入带有毛细排气管的磨口减量瓶中(半密闭条件,附图2),于单一温度下热加速老化,间隔时间取夹心层药片2,获得不同Ⅱ号中定剂有效含量的双基推进剂样品,分为校正集样品和验证集样品。
(二)Ⅱ号中定剂有效含量近红外定量模型
1、使用的装置包括近红外光谱仪、计算机、高效液相色谱法、化学滴定法测定Ⅱ号中定剂有效含量的装置和化学计量学软件等。
依据化学滴定法或高效液相色谱法测定所有双基推进剂样品的Ⅱ号中定剂有效含量,确定含量范围。
2、采用近红外光谱仪对校正集样品和验证集样品采集近红外光谱。采样方式为漫反射,采样波段为700nm-2500nm,优化光谱扫描参数,确定最佳分辨率、扫描次数及样品重复测量次数。采集到的近红外光谱传输至计算机。
3、采用化学计量学方法对校正集样品近红外光谱进行预处理和线性拟合。选择特征谱带,在归一化、基线平滑、一阶求导、二阶求导、多元散射校正光谱预处理方法中选择单独或组合方式对采集光谱进行预处理,继而采用偏最小二乘法即PLS法进行回归拟合,建立校正模型,用验证集样品对所建模型进行外部验证,依据模型的内部检验和外部验证相结合的评价体系,逐步优化模型。
根据建立的Ⅱ号中定剂有效含量近红外定量模型,采用手近红外光谱仪对选定的样品进行Ⅱ号中定剂有效含量测定,获得近红外预测值。化学值与近红外预测值的最大误差小于标准方法重复性误差。
(三)安定期预估方程
双基推进剂装药安定性预估方程
式中:r0表示温差为10℃的反应速率可靠温度系数;
t0表示贮存环境温度安定期,d;
tn表示高温热加速老化时间,d;
Tn表示高温热加速老化试验温度,℃;
T0表示工况或贮存环境温度,℃。
(四)可靠温度系数r0
1、调试安全型油浴烘箱,使其在4个不同温度恒定,温度间隔为10℃。
2、将双基推进剂药饼(附图1)放在毛细排气管磨口减量瓶(半密闭条件,图2),分别置于4个不同温度安全型恒温油浴烘箱中进行热加速老化,根据建立的Ⅱ号中定剂有效含量近红外定量模型,采用近红外光谱仪对Ⅱ号中定剂有效含量进行跟踪测定,以Ⅱ号中定剂有效含量消耗50%作为安定期临界判据,获得4个温度与高温热加速老化临界时间对应关系,应用公式(5),拟合得到温度系数:
r=10-10/B
采用上述方法获得多种典型双基推进剂的温度系数。可靠温度系数的取值偏小,获得安定期短,可靠性高,但取值过小,安定期过于保守,提前销毁造成浪费;可靠温度系数的取值过大,获得安定期长,超过安全贮存寿命时,引起热爆炸。因此,以多种典型双基推进剂温度系数的均值与标准偏差差值的0.9倍作为可靠温度系数r0。
(五)热加速老化试验与安定期
将双基推进剂装药药柱沿径向等分为两段,将两段药柱沿轴向上下叠加制作成两段式药柱(附图3),于单一温度下Tn下热加速老化试验,获得装药两段式老化试样,选取下段药柱叠合面9的中心区域8及药柱表面进行近红外光谱扫描,获得药柱在不同老化时间、不同部位的Ⅱ号中定剂有效含量,以Ⅱ号中定剂有效含量消耗最快的区域作为关重部位,以关重部位Ⅱ号中定剂有效含量消耗50%作为临界判据,获得温度Tn下安定剂消耗的时间tn,依公式(6)外推工况或贮存环境温度下双基推进剂装药安定期。
实施例1
以典型双基推进剂ST-1为代表进行说明。
(1)热加速老化法制备Ⅱ号中定剂有效含量近红外建模样品
典型双基推进剂ST-1药饼三明治式夹心结构试样,由上层药片1、夹心层药片2、下层药片3依次叠加组成药饼(Φ=30mm,上层药片1H1=5mm,夹心层2H2=10mm,下层药片3H3=5mm,附图1),药饼装入带有毛细排气管的磨口减量瓶中(半密闭条件,附图2),于95℃下热加速老化,每5小时取样,获得不同Ⅱ号中定剂有效含量的双基推进剂装药样品,收集60个ST-1样品,其中29个样品(编号为1-29)作为校正集,进行线性拟合和内部交叉验证,31个样品作为外部验证集(编号为30-60)用于外部验证。
(二)Ⅱ号中定剂有效含量近红外定量模型
1、使用的装置包括赛默飞世尔科技近红外光谱仪AntarisⅡ、计算机、化学滴定法测定Ⅱ号中定剂有效含量的装置,以及采用赛默飞世尔科技化学计量学软件TQ Analyst等。
2、采用国家军用标准(GJB770B-2005方法210.1)《中定剂溴化法》测定所有ST-1样品的Ⅱ号中定剂有效含量,其范围为1.53%~0.20%。
3、采用近红外光谱仪对上述选定的双基推进剂装药校正集样品采集近红外光谱。采样方式为漫反射,采样波段为700nm-2500nm,分辨率为8cm-1,扫描64次,每个样品重复测量5次。采集到的近红外光谱通过USB数据线传输至计算机(附图4)。
采用化学计量学软件TQ Analyst对校正集样品近红外光谱进行预处理和线性拟合。选择特征谱带,采用归一化法、基线平滑、一阶求导和多元散射校正等方法进行光谱预处理方法,采用偏最小二乘法(PLS)进行线性回归,建立校正模型,用验证集样品对所建模型进行外部验证,依据模型的内部检验和外部验证相结合的评价体系,逐步优化模型。表1给出了几种不同预处理方法、不同维数时较好光谱模型的参数。
表1较好光谱模型的参数
由表1可以看出,1600nm-2400nm范围内采用归一化、基线平滑和一阶求导组合预处理的模型最优。最终采用该预处理方法,维数选用5时建立的偏最小二乘回归模型为最佳,模型的评价指标R=0.9997,RMSEP=0.055,该回归作为双基推进剂ST-1装药Ⅱ号中定剂有效含量近红外定量模型。
根据建立的双基推进剂ST-1装药Ⅱ号中定剂有效含量近红外定量预测模型,用近红外光谱仪对已用化学滴定法精确测定Ⅱ号中定剂有效含量的10个样品进行预测,详细结果见表2,化学值与近红外预测值的最大误差不大于0.03%。化学值与近红外预测值的最大误差应小于标准方法重复性误差。
表2双基推进剂ST-1装药10个待测样品Ⅱ号中定剂有效含量预测结果
样品编号 | 化学值/% | 近红外预测值/% | 误差/% |
1 | 1.53 | 1.52 | -0.01 |
2 | 1.45 | 1.47 | 0.02 |
3 | 1.37 | 1.37 | 0.00 |
4 | 1.22 | 1.25 | 0.03 |
5 | 1.11 | 1.09 | -0.02 |
6 | 0.98 | 0.97 | 0.01 |
7 | 0.85 | 0.88 | 0.03 |
8 | 0.76 | 0.74 | -0.02 |
9 | 0.68 | 0.70 | 0.02 |
10 | 0.55 | 0.52 | -0.02 |
(三)可靠温度系数r0获得
1、调试安全型油浴烘箱,使其温度恒定在95±1℃、85±1℃、75±1℃和65±1℃
2、将双基推进剂ST-1药饼三明治式夹心结构试样(Φ=30mm,上层药片1H1=5mm,夹心层2H2=10mm,下层药片3H3=5mm,附图1),装入带有毛细排气管的磨口减量瓶中(半密闭条件,附图2),分别置于95±1℃、85±1℃、75±1℃和65±1℃安全型油浴烘箱中进行热加速老化,对于不同老化时间的样品,采用双基推进剂ST-1装药Ⅱ号中定剂有效含量近红外定量模型跟踪测定Ⅱ号中定剂有效含量变化(见表3)。获得4个温度下老化时间与Ⅱ号中定剂有效含量的对应关系。
表3双基推进剂ST-1 4个温度下老化时间与Ⅱ号中定剂有效含量
图5给出了以Ⅱ号中定剂有效含量消耗50%作为安定期的临界判据,获得4组温度与热加速老化临界时间对应关系。
表4不同温度与热加速老化临界时间对应关系
温度/℃ | 95 | 85 | 75 | 65 |
安定期/d | 3.4974 | 15.3673 | 58.5000 | 233.6667 |
应用贝瑟洛特方程拟合,得到:
T=104.21-16.5080lgt(R2=0.9996)
温度系数:
r10=4.03
按照上述方法获得14种典型双基推进剂的温度系数,如下表5所示。14种典型双基推进剂温度系数的均值为4.12,标准偏差0.37,取均值与标准偏差差值的0.9作为可靠温度系数,r0=3.38。
表5多种双基推进剂的温度系数
样品代码 | ST-1 | ST-2 | ST-3 | ST-4 | ST-6 | ST-6 | ST-7 |
温度系数 | 4.03 | 4.36 | 3.53 | 3.55 | 4.23 | 4.55 | 4.53 |
样品代码 | ST-8 | ST-9 | ST-10 | ST-11 | ST-12 | ST-13 | ST-14 |
温度系数 | 4.26 | 4.13 | 4.11 | 4.09 | 4.27 | 4.53 | 3.44 |
(四)热加速老化试验与安定期
双基推进剂装药安定期预估方程
将双基推进剂药柱沿径向等分为两段,将两段药柱沿轴向上下叠加制作成两段式药柱(见附图3),于71℃下热加速老化试验,获得装药两段式老化试样,选取下段药柱叠合面9的中心区域8及药柱表面进行近红外光谱扫描,获得药柱在不同老化时间、不同部位的Ⅱ号中定剂有效含量,最终确定Ⅱ号中定剂有效含量消耗最快的中心区域作为关重部位,以关重部位Ⅱ号中定剂有效含量消耗50%作为临界判据,获得温度71℃安定剂消耗的临界时间tn=102d,采用双基推进剂装药安定期预估方程(6),预估贮存环境温度30℃下的双基推进剂装药安定期为41a。
Claims (1)
1.一种双基推进剂装药安定期快速预估方法,其特征在于,该方法利用加速老化法制备双基推进剂装药Ⅱ号中定剂有效含量近红外定量模型建模样品,采用化学滴定法或高效液相色谱法获取Ⅱ号中定剂有效含量化学值,利用近红外漫反射光谱仪采集光谱,通过化学计量学建立Ⅱ号中定剂有效含量近红外定量模型;以Ⅱ号中定剂有效含量消耗50%作为安定期的临界判据,监测多种典型双基推进剂热加速老化过程安定剂含量变化获得温度系数,温度系数均值与标准偏差差值的0.9作为可靠温度系数r0,建立安定期预估方程;开展单一温度加速老化试验,近红外方法跟踪检测老化过程中双基推进剂装药样品Ⅱ号中定剂有效含量,获得Ⅱ号中定剂有效含量到达临界判据消耗的时间,根据安定期预估方程及可靠温度系数预估双基推进剂装药工况及贮存环境温度下的安定期,具体按下列步骤进行:
(一)加速老化法制备Ⅱ号中定剂有效含量近红外建模样品
双基推进剂药饼三明治式夹心结构试样,由上层药片(1)、夹心层药片(2)、下层药片(3)依次叠加组成药饼,药饼装入带有毛细排气管的磨口减量瓶中,于单一温度下加速老化试验,间隔时间取夹心层药片(2),获得不同Ⅱ号中定剂有效含量的双基推进剂样品,分为校正集样品和验证集样品;
(二)Ⅱ号中定剂有效含量近红外定量模型
1、采用化学滴定法或高效液相色谱法获取双基推进剂样品Ⅱ号中定剂有效含量化学值,确定含量范围;
2、采用近红外光谱仪,对上述选定的样品采集近红外光谱,采样方式为漫反射,特征波段为1600nm-2400nm,分辨率为8cm-1,扫描64次,每个样品重复测量5次,采集到的近红外光谱传输至计算机;
3、采用化学计量学方法对上述近红外光谱进行预处理和线性拟合,综合采用归一化、基线平滑和一阶求导三种方法对采集光谱进行预处理,继而采用偏最小二乘法即PLS法进行回归拟合,建立校正模型,用验证集样品对所建模型进行外部验证,依据模型的内部检验和外部验证相结合的评价体系,逐步优化模型;
根据建立的Ⅱ号中定剂有效含量近红外定量模型,采用近红外光谱仪对选定的样品进行有效安定剂的含量测定,获得近红外预测值,化学值与近红外预测值的最大误差应小于标准方法重复性误差;
(三)安定期预估方程
双基推进剂装药安定期预估方程
式中:r0表示温差为10℃的反应速率可靠温度系数;
t0表示贮存环境温度安定期,d;
tn表示高温加速老化时间,d;
Tn表示高温加速老化试验温度,℃;
T0表示工况或贮存环境温度,℃;
(四)可靠温度系数r0
可靠温度系数保证了安定期内化学安定性良好,以多种典型双基推进剂温度系数r的均值与标准偏差差值的0.9倍作为可靠温度系数r0,经式(2)进行一元线性回归,进而计算得到温度系数r
Tn=A+Blgtn 式(2)
式中A、B为系数,其中A=T0+10×lgt0/lgr、B=-10/lgr
温度系数r=10-10/B,采用近红外方法跟踪测试双基推进剂药饼不同温度Tni下加速老化试验不同老化时间的Ⅱ号中定剂有效含量,获得不同温度Ⅱ号中定剂有效含量随时间变化关系曲线,以Ⅱ号中定剂有效含量消耗50%作为临界判据,获得不同温度Tni对应的时间tni,经代入(2)式进行一元线性回归,采用最小二乘法获得系数B,进而计算得到温度系数r;
采用上述方法获得多种典型双基推进剂装药的温度系数,以其均值与标准偏差差值的0.9倍作为可靠温度系数r0;
(五)加速老化试验与安定期
将双基推进剂药柱沿径向等分为两段,将两段药柱沿轴向上下叠加制作成两段式药柱,于单一温度下Tn下加速老化试验,获得装药两段式老化试样,选取下段药柱叠合面(9)的中心区域(8)及药柱表面进行近红外光谱扫描,获得药柱在不同老化时间、不同部位的Ⅱ号中定剂有效含量,以Ⅱ号中定剂有效含量消耗最快的区域作为关重部位,以关重部位Ⅱ号中定剂有效含量消耗50%作为临界判据,获得温度Tn下安定剂消耗的时间tn,依公式(1)外推贮存环境温度下的装药安定期。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610614757.9A CN106290233B (zh) | 2016-07-29 | 2016-07-29 | 双基推进剂装药安定期快速预估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610614757.9A CN106290233B (zh) | 2016-07-29 | 2016-07-29 | 双基推进剂装药安定期快速预估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106290233A CN106290233A (zh) | 2017-01-04 |
CN106290233B true CN106290233B (zh) | 2019-04-19 |
Family
ID=57663179
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610614757.9A Active CN106290233B (zh) | 2016-07-29 | 2016-07-29 | 双基推进剂装药安定期快速预估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106290233B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110082384B (zh) * | 2019-05-15 | 2021-07-23 | 湖北航天化学技术研究所 | 高能固体推进剂药柱产气产生孔洞或开裂时间预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104792964B (zh) * | 2015-05-04 | 2016-12-07 | 西安近代化学研究所 | 一种基于自重诱发压力条件下nepe推进剂贮存寿命的预估方法 |
-
2016
- 2016-07-29 CN CN201610614757.9A patent/CN106290233B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106290233A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106290236B (zh) | 硝酸酯类火炸药装药安定期快速预估方法 | |
CN106290231B (zh) | 硝酸酯类火炸药装药安全贮存寿命预估方法 | |
CN103018195B (zh) | 近红外光谱测定pbx炸药中pctfe含量的方法 | |
Oliveira et al. | Adulteration of diesel/biodiesel blends by vegetable oil as determined by Fourier transform (FT) near infrared spectrometry and FT-Raman spectroscopy | |
Li et al. | Application of Fourier transform infrared spectroscopy for the quality and safety analysis of fats and oils: A review | |
CN106290237B (zh) | 复合改性双基推进剂装药安定期快速预估方法 | |
de Lima et al. | In-line monitoring of the transesterification reactions for biodiesel production using NIR spectroscopy | |
Nie et al. | Enhanced accuracy of near-infrared spectroscopy for traditional Chinese medicine with competitive adaptive reweighted sampling | |
De Géa Neves et al. | Monitoring of adulteration and purity in coconut oil using raman spectroscopy and multivariate curve resolution | |
CN106290233B (zh) | 双基推进剂装药安定期快速预估方法 | |
CN106290222B (zh) | 复合改性双基推进剂装药安全贮存寿命预估方法 | |
CN106290232B (zh) | 堆积贮存单基发射药安全贮存寿命预估方法 | |
CN103134770A (zh) | 消除水分对近红外光谱检测土壤全氮含量影响的方法 | |
Xie et al. | An approach to detecting diphenylamine content and assessing chemical stability of single-base propellants by near-infrared reflectance spectroscopy | |
Liu et al. | A feasibility quantitative analysis of free fatty acids in polished rice by fourier transform near‐infrared spectroscopy and chemometrics | |
CN106018333B (zh) | 堆积双基发射药安定期快速预估方法 | |
CN106290172B (zh) | 堆积三基发射药安定期快速预估方法 | |
Zhou et al. | Application of near infrared spectroscopy in fast assay of liquid components in single-base propellant intermediates | |
CN106093292B (zh) | 一种nepe推进剂装药安全贮存寿命预估方法 | |
Kartnaller et al. | Evaluating the kinetics of the esterification of oleic acid with homo and heterogeneous catalysts using in-line real-time infrared spectroscopy and partial least squares calibration | |
CN106290223B (zh) | 堆积贮存三基发射药安全贮存寿命预估方法 | |
CN106226266B (zh) | 一种双基推进剂装药安全贮存寿命预估方法 | |
CN106290235B (zh) | 堆积单基发射药安定期快速预估方法 | |
CN106290234B (zh) | 堆积双基发射药安全贮存寿命预估方法 | |
CN102954945B (zh) | 一种由红外光谱测定原油酸值的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |