计算机辅助双基系固体推进剂配方概念设计方法研究

董文静, 赵宏安, 耿国华, 李满荣

(西北大学信息科学与技术学院, 陕西 西安 710127)

1 引 言

随着国防科技技术的不断发展,现代火箭武器的产品用户需求也更加多元化,如何准确有效地明确用户需求,设计出符合客户要求的产品对设计者来说尤为重要。概念设计是产品设计的重要阶段,通常占用整个设计环节的75%以上[1-3]。它直接决定最终产品的能量、安全、经济等性能的优劣程度,因此设计出满足顾客需求的概念设计产品是双基系固体推进剂设计的关键。

传统火炸药配方设计大多是通过寻找满足配方某种需求(或设计指标)的功能组分,再对担任不同化学机理的相关成分进行组合搭配,最终得到配方方案。如Long Y[4]等为了研究低感度炸药,将石墨与高能炸药组合,在获得较大冲击能量的同时,又因为石墨能使炸药能量密度减小,热点形成的概率及灵敏度都降低,设计出了一种具有较大动态冲击性能的奥克托今-石墨混合炸药。William A.MacCrehan等[5]在设计无烟火药时,了解无烟火药的影响因素除了主要推进剂材料硝化纤维外还需要对它的添加剂硝化甘油(NG)、稳定剂二苯胺(DPA)、硝化产物N-nitrosodiphenylamine(NnDPA)以及乙基中定剂(N,N′-diethyl-N,N′-diphenylurea,EC)进行准确分析设计。陆明[6]在进行高威力乳化炸药的设计时,由于氧化剂对乳化炸药的爆炸性能起着关键作用,因此转而研究氧化剂(硝酸铵、硝酸钠)等物质的特性。这些文献都体现了火炸药产品也是以简单的功能(如: 低感度炸药中的降敏功能)-载体(如: 低感度炸药中的石墨成分)的概念设计结构模型的思想。目前,把双基系固体推进剂配方概念设计的思想单独提出并加以研究的报道鲜见,仅查阅到Шварц А Γ[7]在1987年莫斯科国际橡胶会议上报告了对橡胶结构和性能进行现代概念设计的思想; 随后傅圣勇[8]初步提出对水泥建立广义的配方概念,希望全面揭示能量、质量、成本的关系、提高水泥原料资源的配置效率。但传统的双基系固体推进剂配方设计方法中涉及的组分、性能以及功能结构之间的关系复杂,概念知识模糊、抽象,优化设计中常会出现因变量成分的改变、增加而无法求解出满足性能要求的成分含量最终解,出现返工,极大浪费成本。

质量功能展开(Quality Function Deployment,QFD)是一种由顾客需求驱动的产品开发方法[9],在产品的设计和制造过程中,围绕用户需求,令用户得到最大满意度。公理化设计(Axiomatic Design,AD)是由美国Nam P.Suh教授于1990年提出的指导设计过程的基本公理[10],其针对用户需求能进行功能-载体的层次映射推理,该过程十分符合现代化学类产品的设计特点。本研究结合QFD和AD的各自优势,根据双基系固体推进剂配方的设计特点,提出一种基于QFD和公理化设计的双基系固体推进剂配方概念设计模型。本研究提出的配方概念设计理论思想,可对配方进行前期的方案规划,保证了所有产品方案都会有解,该设想为双基系推进剂的配方设计提供新的思路。

2 双基系固体推进剂配方概念设计过程

2.1 传统公理化设计的分析

利用传统公理化设计中功能域到载体域的"Z"型映射分解对双基系固体推进剂配方进行设计,其过程为: 把收集到的顾客需求用质量屋技术转化为技术需求,再将技术需求作为功能域直接映射到载体域,按照自顶向下的循环往复的过程,直至找到最终解(配方成分)。但这种功能到载体(配方成分)的直接映射只能证明“载体(配方成分)实现了某种功能”,并不能说明“为什么载体(配方成分)可以实现此功能”。所以在功能域和载体域之间引入原理域(principle domain)、措施域(measure domain),把两个域作为功能域与载体域之间的转换桥梁,在符合人类思维过程的情况下,既实现对载体的求解,又增加产品的创新设计空间。本文即在QFD的基础上,提出“功能域-原理域-措施域-载体域”的概念设计系统模型。

2.2 双基系固体推进剂配方概念设计系统模型的建立

本文基于QFD和AD建立了概念设计产品模型,将产品在"功能域→原理域→措施域→载体域"模式下进行“Z”字形循环往复映射分解,直到映射结束,得到产品概念设计的方案解。该映射过程体现了公理化设计中的域(domains)、层级(hierarchies)、Z形映射(zigzagging)三大概念[10],采用自顶向下的方式依次向下分解。同时映射过程可分为横向映射和纵向映射,横向映射是指以功能域为起点,依次向原理域、措施域和载体域映射的过程,这是域间的映射; 而纵向映射是指域内的映射,每一个域中的元素都可不断地划分相应的子元素,直到可以被右边相邻域中元素实现为止,如功能域中的功能可划分为子功能,子功能再根据是否能被子原理实现决定是否继续划分。该设计模型以层次设计的方法应用于双基系固体推进剂配方设计中, (1)方案管理便捷化,便于更加灵活、方便的对每个层中的元素进行修改、删除; (2)方案准确化,对所需配方都能找到满足性能需求的最终解; (3)方案呈多元化,可根据不同的用户需求选取配方; (4)方案呈创新性,原理域、措施域的加入使产品方案设计的过程增加更多的创新性。基于QFD和AD的产品概念设计系统模型如图1所示。

图1 基于QFD和AD的产品概念设计系统模型

Fig.1 The system model of product conceptual design based on QFD and AD

2.3 改性双基推进剂域间映射关系结构模型图建立

以改性双基推进剂配方概念设计为例,对功能域、原理域、措施域、载体域之间进行映射推理,得到改性双基推进剂域间映射结构模型图如图2所示。

图2 改性双基推进剂域间映射结构模型图

Fig.2 Model of modified double-base propellant structure mapping between domains

3 双基系固体推进剂概念配方的评价方法及组分含量修正方法

完成对双基系固体推进剂配方的方案设计(即获得了概念配方)后,概念配方是将用户需求运用QFD和 AD的概念设计产品模型,按照"功能域→原理域→措施域→载体域"的层次结构进行“Z”形循环往复映射分解,最后得到满足用户性能需求的配方功能成分,形成的配方设计方案。为了选择最优的配方方案,需要对此概念配方的性能进行初步评价。以配方的性能、代价以及可持续发展为主要考虑因素[11-13],通过采用配方成分(组分)选取的推荐度对方案的设计结果进行评价,使推荐度成为衡量配方优劣程度的指标,为后续配方的详细设计提供指导; 采用含量修正方法对配方组分进行约束,保证以最有效的含量满足顾客需求功能。

3.1 配方组分选取推荐度评价方法

本文中双基系固体推进剂配方组分的推荐度用该配方在组分中功能的(强弱程度)、成本代价(大小)、可持续发展(程度)性能的加权平均值来表示。

功能强弱程度: 功能强弱程度是衡量配方推荐度的重要标准,是衡量组分在配方中担任的各项功能强弱大小的能力。如双基系固体推进剂配方中担任能量性能成分的功能强弱程度; 担任增塑性能成分的功能强弱程度。根据配方中组分担任功能能力的强弱程度,建立功能值定义表如表1所示。

表1 功能值定义

Table 1 The definition of functional numerical value

runctionalperformancedescriptionvalueintheformulaasthisfunctionindependently，andtheindispensable1．0asthefeatureofstrongcapacityintheformula0．5-0．9asthefeatureofweakcapacityintheformulation0．1-0．4withoutthecapacityofrealizingthefunction0

成本代价大小: 代价性是指获取该组分所耗费物质、精力。如通常为了提高配方的能量,往往加入燃烧时单位质量放热量大的金属粉末。铝、铍为两种可添加的金属燃料,以添加铍的理论比冲最高,其次为铝,但因为铍的资源缺乏、价格昂贵、使用代价高,所以未得到实际应用; 铝的来源最广,价格便宜,所以在推进剂等火炸药产品中广泛应用。代价值定义表如表2所示。

表2 代价值定义

Table 2 The definition of cost numerical value

costperformancedescriptionvaluesourceiswide，andthepriceischeap1．0obtainingisrelativelydifficult，thepriceismoreexpensive0．5-0．9obtainingisverydifficultandveryexpensive0．1-0．4don′tget0

可持续发展程度: 可持续发展性是指该组分对于生物即环境的危害程度。如用氢化铝和氢化铍等金属氢化物来提高火炸药能量,虽然都能获得高的理论比冲,但由于铍化物为剧毒物,对环境污染严重,所以未得到实际应用。可持续发展值定义表如表3所示。

表3 可持续发展值定义

Table 3 The definition of sustainable numerical value

sustainableperformancedescriptionvaluedon′tcausepollutiontoecologicalenvironment，notobviouslyimpairpeople′shealth，canrecycleandreuse1．0lessdamagetotheecologicalenvironmentandhumanhealth，mostcanberecycled0．5-0．9harmisbiggerfortheecologicalenvironmentandhumanhealth，asmallnumberofrecyclable0．1-0．4seriousdamagetotheecologicalenvironmentandhumanhealth，cannotberecycled0

配方组分的推荐度加权平均值定义如下:

R(i)=Rf(i)×wf+Rc(i)×wc+Rs(i)×ws

(1)

式中,R(i)为配方组分i的加权总推荐度,Rf(i)、Rc(i)和Rs(i)分别为组分i的功能性数值、代价性数值和可持续发展性数值;wf、wc和ws为相应的权重值,且wf+wc+ws=1,经过验证,本文中的wf、wc、ws定为0.7、0.2、0.1; 权重值的大小由层次分析法[14](analysis hierarchy process,AHP)或德尔菲法[15]确定。

根据以上分析,可定义配方的推荐度计算公式如(2)式:

(2)

式中,H为配方的总推荐度;Nr为配方中组分的个数,wf、wc和ws分别指功能性、代价性、可持续发展性分别在组分中的比重。

表4给出常用配方中常见高能单质炸药各性能推荐值如下。

表4 常见高能单质炸药各性能推荐值

Table 4 Recommended values of performances for common high-energy single compound explosives

acronymfunctionalnumericalvaluecostnumericalvaluesustainablenumericalvalueRDX0．70．90．7HMX0．90．50．5PETN0．80．30．2Tetryl0．60．40．4TNT0．40．90．6

3.2 组分含量的修正

配方中每个子功能需求在用户方面的重要程度是各异的,而组分的含量对各子功能的影响十分重要,因此合理的利用、分配组分含量对最大限度满足用户需求来说是相当必要的。

设有m个组分的配方(1个或多个组分对应同一种功能),不出现一种组分对应不同的功能结构,组分之间满足独立性和完整性。

(1) 选取每种功能下,措施域映射的所有载体使用频率最高的20%。

(2) 对这些载体含量按照所对应功能的权重大小计算组分推荐的含量(这里的含量指的是整个配方总含量为100情况下,每种组分的推荐含量)。

计算方法如下:

C=100/(F1×β1+F2×β2+…+Fm×βm)

(3)

M(i)=C(Fi×βi)

(4)

式中,C表示含量指数;F1,F2,…,Fm表示分别表示第1到m组分的初始推荐最小值,β1,β2,…,βm表示第1到m组方在配方中的权值大小,M(i)表示第i个组分在原始配方总含量为100情况下的推荐含量。β1+β2+βm=1,权重值的大小由经验分析法[14]给出。

(3) 对新的组分含量进行修正。

如果M(i)≥Fi,则组分i的含量为Fi; 如果M(i)

4 双基系固体推进剂配方概念设计系统的实现

现以改性双基推进剂为例,对概念设计方案生成模型进行说明并给出实验结果。

4.1 改性双基推进剂配方概念设计及分析

采用上述概念设计模型、方案结果评价方法以及含量修正方法等设计了一套适用于改性双基推进剂计算机辅助配方概念设计的原型系统,该系统具备根据用户性能需求设计满足其的概念配方、根据类型查询相关配方、配方的具体查看以及用户管理等功能,其系统框架见图3,其配方概念设计系统主要操作流程和关键操作界面见图4和图5。

图3 改性双基推进剂概念设计系统框架图

Fig.3 The framework map of modified double-base propellant concept design system

图4 配方概念设计主要操作流程

Fig.4 The main operation processes of designing conceptual formulation

图5 配方概念设计关键操作界面

Fig.5 The key operation interface of conceptual design for formulation

为了更方便地展示该系统的工作原理和过程,现假设需要设计一个改性双基固体推进剂,要求: 能量性能比冲为2400～2450 N·s·kg-1; 低特征信号为无烟或微烟,力学性能要具有良好的低温性能,燃速可控。将用户以上需求指标输入到配方概念设计系统中,系统依据用户需求,根据功能-原理-措施-载体的概念设计过程模型进行映射推理,完成配方成分的选择及各成分的含量确定。最后利用性能预估软件对配方方案进行性能预估,如果性能满足用户需求,则设计结束; 若不满足,则重新调整原理、措施、组分进行再设计。其设计结果如图6所示,图中显示了配方的主要组分以及相应的含量。

图6 改性双基推进剂配方设计结果

Fig.6 The design result of formulation for modified double-base propellant

经分析,由概念设计模型设计的新配方与已测定性能配方库中的某配方Ⅰ[14]的能量比冲很接近。二者的成分组成、各成分含量的对比及性能对比见表5和表6。

由表5可知,两个配方的主要功能成分相同,硝化纤维(N=12%)、硝化甘油、奥克托今的主要能量成分十分接近。另外由新配方的成分组成也可以看出新配方属于改性双基推进剂,符合用户需求的配方类型(设计过程中用到图1的概念设计模型),由表5可知,由配方概念设计系统所设计出来的新配方经性能预估得出的比冲与配方Ⅰ中已测定比冲相差不大。

二者的差异如下: (1)新设计的配方没有铝粉、碳黑,由于在设计配方时依照如前所述的用户需求为无烟推进剂,根据产生烟气的原理,在设计过程中会避免选择铝、铅、高氯酸铵等物质,防止产生三氧化二铝的白色固体颗粒烟雾、产生氯化氢白色气体烟雾。(2)新设计的配方中加入了吉纳(DINA),由于在设计配方时要求配方具有低温力学性能,就是要有良好延伸率,因此在配方中加入相应的增塑剂(溶剂)可使大分子的作用力降低,玻璃化温度降低、模量降低、延伸率提高。所以在配方中加入吉纳,而且吉纳属于良好的含能增塑剂,对于配方能量的提升也具有很好的作用。(3)由于要求推进剂燃速可控,故氧化铜作为双基推进剂中十分有效的催化剂加入新配方中。(4)虽然新配方与配方Ⅰ在小含量成分上存在一定的差异,但不会特别影响配方的能量性能。因此可得出,概念配方设计系统在概念设计模型、配方成分推荐以及含量修正方法的运用下得出的配方是合理的。

表5 新设计配方与某配方Ⅰ的成分组成及成分含量的对比

Table 5 Comparsion of the component and its content of new design formulation and formulation Ⅰ

newdesignformulationnamecontent/%formulationⅠnamecontent/%NC(N=12．0%)50．09NC(N=12．0%)51．0NG26．92NG25．3HMX13．85HMX12．8DINA6．92aluminum42，4⁃dinitrodiphenylamine1．38carbonblack3oxideCu0．69C2-vaseline(VSL)0．14vaseline(VSL)-

表6 新设计的配方与某配方Ⅰ的性能对比

Table 6 Comparsion of performances for new design formulation and formulation Ⅰ

formulationspecialimpulse/N·s·kg-1signaturepropellanttypenewdesignformulation1)2444．06nosmokedoubleBaseformulationⅠ2)2445．86smokedoubleBase

Note: 1) specific impulse estimated; 2) specific impulse determined.

由系统设计的某两种双基推进剂配方A与配方B,其主要性能指标密度都为1.6 g·cm-3,爆热都为3640 kJ·kg-1,为了选择最优的配方,对两配方的功能值、代价值、可持续发展值进行确定,并计算配方推荐度,如表6所示,由于配方A的推荐度3.93大于配方B的推荐度3.914,因此配方A为较优配方。

双基系固体推进剂配方概念设计系统通过对用户关于配方的需求进行归纳、整理,继而总结为一般的技术功能需求,通过结合配方设计专业知识推理映射出相应实现此功能的原理、措施,进而得到配方的具体组分,并对其含量进行修正,利用综合评价方法进行方案优选,最终得到符合相应性能指标的概念设计方案。实验结果表明,采用本研究的配方概念设计模型和配方评价以及含量修正方法给出的概念配方能为后续配方详细设计工作起到指导作用。

表6 配方A和配方B的推荐度对比

Table 6 Comparsion of recommended degree of formulation A and formulation B

basicinformationcompositioncontent/%theamountofplaninformationfunctionalnumericalvaluecostnumericalvaluesustainablenumericalvaluerecommendeddegreeformulationANC⁃120540．90．850．8NG270．90．80．6DNT150．70．750．2VSL20．90．20．6oxideCu20．80．90．73．93formulationBNC⁃120570．90．90．8NG280．90．820．6centralite110．90．60．6DNT110．70．70．2VSL10．90．10．63．914

5 结 论

(1)将质量功能展开与公理化设计的优势相结合,根据双基系固体推进剂配方的产品需求,提出一种基于QFD和公理化设计的双基系固体推进剂配方概念设计系统模型,并利用配方概念设计系统设计的一则改性双基推进剂配方实例展示了该模型在推进剂配方概念设计系统中的运用效果。结果表明该模型能够有效地对用户需求实施映射分解,最终获取满足用户需求指标的配方最终解(配方成分)。

(2)提出了双基系固体推进剂配方组分的评价方法和组分含量修正方法。通过实验表明,在其支持下,可以设计出在最低成分含量下满足用户性能的配方,并根据评价指标对配方进行优选。两种方法的应用既避免了实际配方时的浪费,又对设计者的决策选择起到了指导作用。

(3)相较于传统的配方设计方法,前期概念设计的引入,能够通过由功能需求到载体成分的映射过程保证配方一定有方案解,避免了生产中的重复返工; 而且通过用户需求的科学规划设计能大量地减少详细设计的工作量,缩短了产品开发周期,为双基系固体推进剂配方设计领域提供了一种新的切实可行的方法。

参考文献:

[1] Wang Z, El-Gizawy A S, Billatos S B. A knowledge-based system methodology for conceptual design of mechanical systems[J].AmericanSocietyofMechanicalEngineersDesignEngineeringDivision, 1996: 1-12.

[2] Casakin H. Metaphors in design problem-solving: implications for creativity[J].InternationalJournalofDesign, 2007, 1(2): 23-35.

[3] Sas C, Whittaker S, Dow S, et al. Generating implications for design through design research[C]∥ACM Conference on Human Factors in Computing Systems. ACM, 2014: 1971-1980.

[4] Long Y, Liu Y G, Nie F D, et al. Theoretical study of impacting and desensitizing for HMX-graphite mixture explosive[J].ShockWaves, 2012, 22(6): 605-614.

[5] Maccrehan W A, Bedner M. Development of a smokeless powder reference material for propellant and explosives analysis[J].ForensicScienceInternational, 2006, 163(1): 119-124.

[6] 陆明. 乳化炸药氧化剂选择的理论研究[J]. 爆破器材, 2006(4): 4-7.

LU Ming.Theory study on the selection of oxidizer for emulsion explosive[J].ExplosiveMaterials, 2006(4): 4-7.

[7] Шварц А Γ, 陈根度. 按照橡胶结构和性能的现代概念设计配方[J]. 橡胶参考资料, 1986(9): 1-10.

Шварц А Γ, CHEN Gen-du. The modern concept design formula of the rubber structure and performance[J].RubberResources, 1986(9): 1-10.

[8] 傅圣勇, 秦至刚, 李启钧.建立广义的配方概念(下)——论水泥烧成的原料配制效率[J]. 四川水泥, 2003(1): 5-9.

FU Sheng-yong, QIN Zhi-gang, LI Qi-jun. Establish a generalized formula concept (below)——theory of cement raw material preparation efficiency of burning[J].SichuanCementMagazine, 2003(1): 5-9.

[9] Brown P G. QFD: Echoing the voice of the customer[J].BellLabsTechnicalJournal, 1991, 70(2): 18-32.

[10] SUH N P. The principles of design[M]. London: Oxford University Press, 1990.

[11] 蒋浩龙, 王晓峰, 陈松, 等 炸药混合技术的发展和应用[J]. 飞航导弹, 2014(12): 74-78.

JANG Hao-long, WANG Xiao-feng, CHEN Song, et al. The development and application of explosive mixture technology[J].AerodynamicMissileJournal, 2014 (12): 74-78.

[12] 李志峰, 迟淑萍, 宋仁峰. 国内外乳化炸药现状及技术发展建议[J]. 矿业工程, 2011(5): 42-44.

LI Zhi-feng, CHI Shu-ping, SONG Ren-feng. Recent status and development trend of home and abroad emulsive explosive[J].MiningEngineering, 2011(5): 42-44.

[13] 任晓雪. 国外工业炸药的研究与发展[J]. 火炸药学报, 2011(5): 50-53.

REN Xiao-xue. Study and development of foreign commercial explosives[J].ChineseJournalofExplosivesandPropellants, 2011(5): 50-53.

[14] 西北大学·西安近代化学研究所.“改性双基推进剂配方设计专家系统”技术总结报告[R]. 西安: 西北大学·西安近代化学研究所, 2013.

North West University·Xi′an Modern Chemistry Research Institute. "The modified double-base propellant design expert system" technology summary report[R]. Xi′an: Northwestern University·Xi′an Modern Chemistry Research Institute, 2013.

[15] 秦吉, 张翼鹏. 现代统计信息分析技术在安全工程方面的应用——层次分析法原理[J]. 工业安全与防尘, 1999 (5): 44-48.

QIN Ji, ZHANG Yi-peng. Application of contemporary statistical information analysis method in safety engineering-the principle of AHP[J].IndustrialSafetyandEnvironmentalProtection, 1999(5): 44-48.

[16] Delphi method [ EB/OL] .http:∥en.wikipedia.org/wiki/Delphi method [2011-04-11].

[17] 李志琴, 赵宏安, 赵凤起, 等.基于语义网络的推进剂配方设计知识表示[J]. 含能材料, 2014, 22(4): 542-547.

LI Zhi-qin, ZHAO Hong-an, ZHAO Feng-qi, et al.Knowledge representation of propellant formula design based on semantic network[J].ChineseJournalofEnergeticMaterials(HannengCailiao), 2014, 22(4): 542-547.