CN110688771A - 一种基于渗透扩散的片状多层发射药燃烧性能计算方法 - Google Patents
一种基于渗透扩散的片状多层发射药燃烧性能计算方法 Download PDFInfo
- Publication number
- CN110688771A CN110688771A CN201910969062.6A CN201910969062A CN110688771A CN 110688771 A CN110688771 A CN 110688771A CN 201910969062 A CN201910969062 A CN 201910969062A CN 110688771 A CN110688771 A CN 110688771A
- Authority
- CN
- China
- Prior art keywords
- propellant
- layer
- multilayer
- combustion
- sheet
- 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.)
- Granted
Links
- 239000003380 propellant Substances 0.000 title claims abstract description 92
- 238000002485 combustion reaction Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000009792 diffusion process Methods 0.000 title claims abstract description 11
- 230000003204 osmotic effect Effects 0.000 title abstract 2
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 230000007704 transition Effects 0.000 claims abstract description 16
- 230000035515 penetration Effects 0.000 claims description 7
- 239000003721 gunpowder Substances 0.000 claims description 2
- 230000008595 infiltration Effects 0.000 claims 1
- 238000001764 infiltration Methods 0.000 claims 1
- 230000001965 increasing effect Effects 0.000 abstract description 4
- 239000000843 powder Substances 0.000 abstract 7
- 238000011156 evaluation Methods 0.000 abstract 1
- 239000000463 material Substances 0.000 description 7
- 238000001125 extrusion Methods 0.000 description 4
- 238000012512 characterization method Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000012856 packing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及发射药技术领域,具体涉及一种基于渗透扩散的片状多层发射药燃烧性能计算方法。本发明解决了片状多层发射药的燃烧特性的表征方法,提供了一种燃烧渐增性评估的技术手段。采用由外层向内层方向的同一坐标方式,能够实现任意多层发射药的参数表征,并且随着层数的增加并不增加计算的复杂度。对片状多层发射药的过渡层采用渐变函数的处理方法,能够更符合物理过程,较好的对片状多层发射药的燃烧渐增性进行表征。所建立的方法不需要指定缓燃层和速燃层,因此,能够适用于各类复杂结构的发射药。
Description
技术领域
本发明涉及发射药技术领域,具体涉及一种基于渗透扩散的片状多层发射药燃烧性能计算方法。
背景技术
片状多层发射药是由至少三层不同组分的发射药燃烧层组成,每一层具有不同的燃烧速度,经过物理挤压的方式结合在一起,具有一定强度、几何外形的复合式发射药。该类发射药具有较高的燃烧渐增性,同时发射药具有较高的装填密度,对比粒状散装发射药,装填密度能够提高25%以上,最高可以达到1.25g/cm3以上。应用多层发射药的高渐增性燃烧特性,可以在保持最大膛压不变的条件下提高炮口动能15%-50%,因此,该类技术的研究具有重大意义。
现在的片状多层发射药的燃烧性能计算方法主要是认为缓燃层和速燃层的密度相同,并且不考虑燃烧过程中缓燃层与速燃层边界的影响。然而,现有计算方法存在如下问题:(1)认为不同分层的发射药密度均相同,且层与层之间为直接连接的方式,这显然与实际发射药的组成不同,由于发射药材料的密度范围从1.2g/cm3到1.8g/cm3,跨度较大,因此,在挤压成型后的不同层之间,物料为两种材料的混合物,现有方法无法较好的表征燃烧渐增性,尤其是层与层之间的过渡部分;(2)在发射药的燃烧过程中不考虑缓燃层与速燃层边界的影响,由于边界处在经过物理挤压和化学渗透后已经形成了一种具有渐变特性的新材料,相对于发射药的整体尺寸而言,所占比例较大,如弧厚为0.5mm的发射药,边界层厚度约0.1mm,因此,现有方法无法较好的表征片状多层发射药的燃烧特性。
从目前检索的技术资料中,尚未见具有能够考虑发射药层与层之间的过渡部分组成及其燃烧特性的片状多层发射药燃烧性能计算方法。
发明内容
为解决现有技术无法较好的表征片状多层发射药的燃烧渐增性、燃烧特性的问题,本发明提供一种基于渗透扩散的片状多层发射药燃烧性能计算方法。
为了解决上述问题,本发明采用了下列技术方案:
一种基于渗透扩散的片状多层发射药燃烧性能计算方法,包括以下步骤:
步骤1:确定片状多层发射药的层数;
步骤2:确定某时刻的x,x表示此时刻以燃烧的火药的弧厚;
步骤3:计算发射药已燃百分数;
步骤4:计算发射药燃气生成猛度;
步骤5:通过以上计算步骤,即可计算出x点处的片状多层发射药的燃烧特性值;
步骤6:重复步骤(2)、(3)和(4),获得所有点处的燃烧特性数值。
基于片状多层发射药在物理挤压中组分的渗透扩散,采用同一坐标尺度的材料特性分布,对过渡层采用渐变函数的处理方法,能够获得较符合实际过程的发射药燃烧特性数据,较好的对片状多层发射药的燃烧特性能进行表征。
进一步地,所述步骤2中的给出发射药不同层的特性与x的关系,具体方法为:
片状多层发射药由外层向内层的方向为x的正向,发射药不同层的特性与x的关系为
其中
ρ为发射药的密度,ρ1为第一层的密度,ρ2为过渡层的密度,ρ3为第三层的密度,u为发射药的燃速,u1为第一层的燃速,u2为过渡层的燃速,u3为第三层的燃速,x1为第一层结束时的坐标,x2为过渡层结束时的坐标,xt为发射药总的厚度的坐标。增加了过渡层密度ρ2、过渡层燃速u2和过渡层结束时的坐标x2,可使燃烧过程完整,具有燃烧的连续性。
再进一步地,所述步骤3中的计算发射药已燃百分数,具体方法为:
ψ=χZ(1+λZ)0≤Z≤1
其中,ψ为发射药的已燃百分数,χ、λ为发射药的药形系数,Z为燃烧层的已燃相对厚度。
更进一步地,所述步骤4中的计算发射药燃气生成猛度,具体方法为:
假设发射药的燃速采用正比式表示,即
u=kp
进而
其中,Γ为发射药燃气生成猛度,p为压力,k为燃速系数,t为燃烧时间。有效解决了Γ-ψ曲线的突越点,变得有连续性,更接近实际燃烧曲线。
与现有技术相比,本发明的有益效果体现在以下几个方面:
(1)采用由外层向内层方向的同一坐标方式,使得发射药在燃烧过程的特性参数全部转变为与该坐标为变量的单一关系,能够实现任意多层发射药的参数表征,并且随着层数的增加并不增加计算的复杂度。与现有技术相比,实现了发射药特性参数用同一参数表征,拓展了计算范围,降低了计算难度。
(2)对片状多层发射药的过渡层采用渐变函数的处理方法,能够更符合物理过程,较好的对片状多层发射药的燃烧渐增性进行表征。与现有技术相比,实现了燃烧渐增性的平滑过渡表征,为片状多层发射药的燃烧特性的表征方法提供了较好的技术手段,填补了现有技术空白。
(3)所建立的方法不需要指定缓燃层和速燃层,因此,能够使用于各类复杂结构的发射药。与现有技术相比,适用范围更广,能够应用于各类具有层状结构的发射药。
因此,本发明解决了片状多层发射药的燃烧特性的表征方法,提供了一种燃烧渐增性评估的技术手段。
附图说明
图1是片状多层发射药燃烧特性的计算流程图;
图2是计算获得的Γ-ψ曲线。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
本实施例中一种基于渗透扩散的片状多层发射药燃烧性能计算方法,包括以下步骤:
(1)确定片状多层发射药的层数
一种片状多层发射药由三层基体药组成,上层和下层为相同基体的材料,厚度0.5mm,中间为其他的基体材料,厚度为1mm,物理压实后的总厚度为2mm,呈现外直径为45mm的圆盘状,过渡层厚度为0.2mm。
(2)给出发射药不同层的特性与x关系
片状多层发射药由外层向内层的方向为x的正向,并且发射药不同部位的特性与x的关系为
其中
通过以上公式,即可获得在不同位置处发射药的密度及燃速。
(3)计算发射药的已燃百分数
根据已有的药形尺寸数据,可求得该片状多层发射药的药形函数χ=1.0442,λ=-0.0423。
发射药的已燃百分数可以用下式表示
ψ=1.0442×Z×(1-0.0423×Z)0≤Z≤1
发射药的燃速采用正比式的表达式表示,以上所给的基体材料燃速为在100MPa下的燃速数值,单位为cm/s。根据燃速的正比式计算公式
u=kp
进而获得
其中
(4)计算发射药的燃气生成猛度
发射药燃气的生成猛度可以用下式表示
将Z离散化后,步长取0.001。
(5)通过上述计算步骤可计算出在x点处的燃烧特性数值。
(6)重复步骤2、3和4,逐次求解,即可计算出所有点处的片状多层发射药的燃烧特性值。所计算的Γ-ψ曲线如图2所示。
Claims (4)
1.一种基于渗透扩散的片状多层发射药燃烧性能计算方法,其特征在于:包括以下步骤:
步骤1:确定片状多层发射药的层数;
步骤2:确定某时刻的x,x表示此时刻已燃烧的火药的弧厚;
步骤3:计算发射药的已燃百分数;
步骤4:计算发射药的燃气生成猛度;
步骤5:通过以上计算步骤,即可计算出x点处的片状多层发射药的燃烧特性值;
步骤6:重复步骤(2)、(3)和(4),获得所有点处的燃烧特性数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910969062.6A CN110688771B (zh) | 2019-10-12 | 2019-10-12 | 一种基于渗透扩散的片状多层发射药燃烧性能计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910969062.6A CN110688771B (zh) | 2019-10-12 | 2019-10-12 | 一种基于渗透扩散的片状多层发射药燃烧性能计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110688771A true CN110688771A (zh) | 2020-01-14 |
CN110688771B CN110688771B (zh) | 2022-12-02 |
Family
ID=69112503
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910969062.6A Active CN110688771B (zh) | 2019-10-12 | 2019-10-12 | 一种基于渗透扩散的片状多层发射药燃烧性能计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110688771B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6419129B1 (en) * | 1994-06-02 | 2002-07-16 | Robert Henry Abplanalp | Flexible barrier member useful in aerosol dispensers |
US20020182436A1 (en) * | 2000-05-02 | 2002-12-05 | Weihs Timothy P. | Freestanding reactive multilayer foils |
CN101907427A (zh) * | 2010-07-15 | 2010-12-08 | 西安近代化学研究所 | 发射药变容燃速测试装置 |
WO2011006805A2 (fr) * | 2009-07-16 | 2011-01-20 | Astrium Sas | Procédé de réalisation de liaisons mécaniques |
CN103708985A (zh) * | 2013-12-19 | 2014-04-09 | 中北大学 | 一种包覆多微孔结构渐增性燃烧的固体发射药 |
WO2014141196A1 (en) * | 2013-03-15 | 2014-09-18 | Glaxosmithkline Intellectual Property (No.2) Limited | Liquid dispensing devices and methods of controlling the same |
CN104180721A (zh) * | 2014-08-25 | 2014-12-03 | 西安近代化学研究所 | 一种底部多层等离子体点火装置 |
JP2016130617A (ja) * | 2015-01-15 | 2016-07-21 | 日油技研工業株式会社 | 焼尽性容器 |
CN106202797A (zh) * | 2016-07-21 | 2016-12-07 | 中国人民解放军装甲兵工程学院 | 一种火药燃烧时炮管内等离子体密度计算方法 |
CN107353170A (zh) * | 2017-08-23 | 2017-11-17 | 西安近代化学研究所 | 一种发射药压伸工艺参数的动态计算方法 |
CN107480388A (zh) * | 2017-08-23 | 2017-12-15 | 西安近代化学研究所 | 一种预定设计尺寸发射药的压伸成型方法 |
-
2019
- 2019-10-12 CN CN201910969062.6A patent/CN110688771B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6419129B1 (en) * | 1994-06-02 | 2002-07-16 | Robert Henry Abplanalp | Flexible barrier member useful in aerosol dispensers |
US20020182436A1 (en) * | 2000-05-02 | 2002-12-05 | Weihs Timothy P. | Freestanding reactive multilayer foils |
WO2011006805A2 (fr) * | 2009-07-16 | 2011-01-20 | Astrium Sas | Procédé de réalisation de liaisons mécaniques |
CN101907427A (zh) * | 2010-07-15 | 2010-12-08 | 西安近代化学研究所 | 发射药变容燃速测试装置 |
WO2014141196A1 (en) * | 2013-03-15 | 2014-09-18 | Glaxosmithkline Intellectual Property (No.2) Limited | Liquid dispensing devices and methods of controlling the same |
CN103708985A (zh) * | 2013-12-19 | 2014-04-09 | 中北大学 | 一种包覆多微孔结构渐增性燃烧的固体发射药 |
CN104180721A (zh) * | 2014-08-25 | 2014-12-03 | 西安近代化学研究所 | 一种底部多层等离子体点火装置 |
JP2016130617A (ja) * | 2015-01-15 | 2016-07-21 | 日油技研工業株式会社 | 焼尽性容器 |
CN106202797A (zh) * | 2016-07-21 | 2016-12-07 | 中国人民解放军装甲兵工程学院 | 一种火药燃烧时炮管内等离子体密度计算方法 |
CN107353170A (zh) * | 2017-08-23 | 2017-11-17 | 西安近代化学研究所 | 一种发射药压伸工艺参数的动态计算方法 |
CN107480388A (zh) * | 2017-08-23 | 2017-12-15 | 西安近代化学研究所 | 一种预定设计尺寸发射药的压伸成型方法 |
Non-Patent Citations (2)
Title |
---|
SUMITVERMA 等: "Dependence of density and burning rate of composite solid propellant on mixer size", 《ACTA ASTRONAUTICA》 * |
刘海情 等: "球形变燃速发射药的燃气生成规律", 《广州化工》 * |
Also Published As
Publication number | Publication date |
---|---|
CN110688771B (zh) | 2022-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5625672B2 (ja) | 変性ポリテトラフルオロエチレンファインパウダー及び変性ポリテトラフルオロエチレン成形体 | |
HA et al. | Molecular dynamics simulation of sintering and surface premelting of silver nanoparticles | |
JP2014169458A5 (zh) | ||
Benson et al. | Dynamic compaction of copper powder: Computation and experiment | |
CN110688771A (zh) | 一种基于渗透扩散的片状多层发射药燃烧性能计算方法 | |
Whitworth | CREST modelling of PBX 9502 corner turning experiments at different initial temperatures | |
CN106834807A (zh) | 一种石墨烯增强双相金属铝化物复合材料及其制备方法 | |
CN110765407B (zh) | 一种片状多层发射装药内弹道特性参数的计算方法 | |
US11976015B2 (en) | Nanoenergetic material composite-based solid propellant, method of preparing the same, and projectile using the same | |
CN107145621A (zh) | 一种提高SiCp/2024Al基复合材料挤压致密度均匀性的方法 | |
Reaugh | Computer simulations to study the explosive consolidation of powders into rods | |
Beckstead | An overview of aluminum agglomeration modeling | |
KR20200041115A (ko) | 나노고에너지물질과 폴리머 바인더의 복합체 제조 방법 및 이를 이용하여 제조된 복합체와 이의 연소 특성 측정 및 분석 방법 | |
CN116682511A (zh) | 基于近球形粉末的钛合金增材制造工艺参数的设计方法 | |
Xiao et al. | Interior ballistic prediction of gun propellants based on experimental pressure-apparent burning rate model in closed vessel | |
CN116796630A (zh) | 基于机器学习的高性能制件的增材制造制备方法 | |
CN115619930A (zh) | 一种基于参数化体积布尔运算的燃面估计方法 | |
CN110348072B (zh) | 一种提高电弧增材构件热力分析有限元数值计算效率方法 | |
CN111027222A (zh) | 一种航天器Lambert变轨的多脉冲发动机总冲修正方法及系统 | |
Jixiang et al. | Simulation of the neck growth of non-isometric biosphere during initial sintering | |
Lamy et al. | Modeling the porosity evolution of a powder under uniaxial compression | |
CN117852383A (zh) | 一种基于机器学习的压伸发射药能量性能计算方法 | |
CN107563110A (zh) | 一种rdx与石蜡缺陷晶体复合体系的模拟方法 | |
CN113744812A (zh) | 一种渗透工艺球扁发射药的燃速计算方法 | |
Borisovskaya et al. | Mechanical properties of powder titanium at different production stages. I. Densification curves for titanium powder billets |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |