CN110688771B - 一种基于渗透扩散的片状多层发射药燃烧性能计算方法 - Google Patents

一种基于渗透扩散的片状多层发射药燃烧性能计算方法 Download PDF

Info

Publication number
CN110688771B
CN110688771B CN201910969062.6A CN201910969062A CN110688771B CN 110688771 B CN110688771 B CN 110688771B CN 201910969062 A CN201910969062 A CN 201910969062A CN 110688771 B CN110688771 B CN 110688771B
Authority
CN
China
Prior art keywords
propellant
layer
combustion
propellant powder
multilayer
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
Application number
CN201910969062.6A
Other languages
English (en)
Other versions
CN110688771A (zh
Inventor
马中亮
何利明
贺增弟
代淑兰
梁泰鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North University of China
Original Assignee
North University of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by North University of China filed Critical North University of China
Priority to CN201910969062.6A priority Critical patent/CN110688771B/zh
Publication of CN110688771A publication Critical patent/CN110688771A/zh
Application granted granted Critical
Publication of CN110688771B publication Critical patent/CN110688771B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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的关系为
Figure BDA0002231473790000021
Figure BDA0002231473790000031
其中
Figure BDA0002231473790000032
Figure BDA0002231473790000033
ρ为发射药的密度,ρ1为第一层的密度,ρ2为过渡层的密度,ρ3为第三层的密度,u为发射药的燃速,u1为第一层的燃速,u2为过渡层的燃速,u3为第三层的燃速,x1为第一层结束时的坐标,x2为过渡层结束时的坐标,xt为发射药总的厚度的坐标。增加了过渡层密度ρ2、过渡层燃速u2和过渡层结束时的坐标x2,可使燃烧过程完整,具有燃烧的连续性。
再进一步地,所述步骤3中的计算发射药已燃百分数,具体方法为:
ψ=χZ(1+λZ)0≤Z≤1
Figure BDA0002231473790000034
其中,ψ为发射药的已燃百分数,χ、λ为发射药的药形系数,Z为燃烧层的已燃相对厚度。
更进一步地,所述步骤4中的计算发射药燃气生成猛度,具体方法为:
Figure BDA0002231473790000035
假设发射药的燃速采用正比式表示,即
u=kp
进而
Figure BDA0002231473790000036
其中,Γ为发射药燃气生成猛度,p为压力,k为燃速系数,t为燃烧时间。有效解决了Γ-ψ曲线的突越点,变得有连续性,更接近实际燃烧曲线。
与现有技术相比,本发明的有益效果体现在以下几个方面:
(1)采用由外层向内层方向的同一坐标方式,使得发射药在燃烧过程的特性参数全部转变为与该坐标为变量的单一关系,能够实现任意多层发射药的参数表征,并且随着层数的增加并不增加计算的复杂度。与现有技术相比,实现了发射药特性参数用同一参数表征,拓展了计算范围,降低了计算难度。
(2)对片状多层发射药的过渡层采用渐变函数的处理方法,能够更符合物理过程,较好的对片状多层发射药的燃烧渐增性进行表征。与现有技术相比,实现了燃烧渐增性的平滑过渡表征,为片状多层发射药的燃烧特性的表征方法提供了较好的技术手段,填补了现有技术空白。
(3)所建立的方法不需要指定缓燃层和速燃层,因此,能够使用于各类复杂结构的发射药。与现有技术相比,适用范围更广,能够应用于各类具有层状结构的发射药。
因此,本发明解决了片状多层发射药的燃烧特性的表征方法,提供了一种燃烧渐增性评估的技术手段。
附图说明
图1是片状多层发射药燃烧特性的计算流程图;
图2是计算获得的Γ-ψ曲线。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
本实施例中一种基于渗透扩散的片状多层发射药燃烧性能计算方法,包括以下步骤:
(1)确定片状多层发射药的层数
一种片状多层发射药由三层基体药组成,上层和下层为相同基体的材料,厚度0.5mm,中间为其他的基体材料,厚度为1mm,物理压实后的总厚度为2mm,呈现外直径为45mm的圆盘状,过渡层厚度为0.2mm。
(2)给出发射药不同层的特性与x关系
片状多层发射药由外层向内层的方向为x的正向,并且发射药不同部位的特性与x的关系为
Figure BDA0002231473790000051
Figure BDA0002231473790000052
其中
Figure BDA0002231473790000053
Figure BDA0002231473790000054
通过以上公式,即可获得在不同位置处发射药的密度及燃速。
(3)计算发射药的已燃百分数
根据已有的药形尺寸数据,可求得该片状多层发射药的药形函数χ=1.0442,λ=-0.0423。
发射药的已燃百分数可以用下式表示
ψ=1.0442×Z×(1-0.0423×Z)0≤Z≤1
Figure BDA0002231473790000055
发射药的燃速采用正比式的表达式表示,以上所给的基体材料燃速为在100MPa下的燃速数值,单位为cm/s。根据燃速的正比式计算公式
u=kp
进而获得
Figure BDA0002231473790000061
其中
Figure BDA0002231473790000062
(4)计算发射药的燃气生成猛度
发射药燃气的生成猛度可以用下式表示
Figure BDA0002231473790000063
将Z离散化后,步长取0.001。
(5)通过上述计算步骤可计算出在x点处的燃烧特性数值。
(6)重复步骤2、3和4,逐次求解,即可计算出所有点处的片状多层发射药的燃烧特性值。所计算的Γ-ψ曲线如图2所示。

Claims (3)

1.一种基于渗透扩散的片状多层发射药燃烧性能计算方法,其特征在于:包括以下步骤:
步骤1:确定片状多层发射药的层数;
步骤2:确定某时刻的x,x表示此时刻已燃烧的火药的弧厚;
步骤3:计算发射药的已燃百分数;
步骤4:计算发射药的燃气生成猛度;
步骤5:通过以上计算步骤,即可计算出x点处的片状多层发射药的燃烧特性值;
步骤6:重复步骤(2)、(3)和(4),获得所有点处的燃烧特性数值;
所述步骤2中的给出发射药不同层的特性与x的关系,具体方法为:
片状多层发射药由外层向内层的方向为x的正向,发射药不同层的特性与x的关系为
Figure FDA0003883257820000011
Figure FDA0003883257820000012
其中
Figure FDA0003883257820000013
Figure FDA0003883257820000014
ρ为发射药的密度,ρ1为第一层的密度,ρ2为过渡层的密度,ρ3为第三层的密度,u为发射药的燃速,u1为第一层的燃速,u2为过渡层的燃速,u3为第三层的燃速,x1为第一层结束时的坐标,x2为过渡层结束时的坐标,xt为发射药总的厚度的坐标。
2.根据权利要求1所述的一种基于渗透扩散的片状多层发射药燃烧性能计算方法,其特征在于:所述步骤3中的计算发射药已燃百分数,具体方法为:
ψ=χZ(1+λZ) 0≤Z≤1
Figure FDA0003883257820000021
其中,ψ为发射药的已燃百分数,χ、λ为发射药的药形系数,Z为燃烧层的已燃相对厚度。
3.根据权利要求1所述的一种基于渗透扩散的片状多层发射药燃烧性能计算方法,其特征在于:所述步骤4中的计算发射药燃气生成猛度,具体方法为:
Figure FDA0003883257820000022
假设发射药的燃速采用正比式表示,即
u=kp
进而
Figure FDA0003883257820000023
其中,Γ为发射药燃气生成猛度,p为压力,k为燃速系数,t为燃烧时间。
CN201910969062.6A 2019-10-12 2019-10-12 一种基于渗透扩散的片状多层发射药燃烧性能计算方法 Active CN110688771B (zh)

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 CN110688771A (zh) 2020-01-14
CN110688771B true 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 (8)

* Cited by examiner, † Cited by third party
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
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 西安近代化学研究所 一种底部多层等离子体点火装置
CN106202797A (zh) * 2016-07-21 2016-12-07 中国人民解放军装甲兵工程学院 一种火药燃烧时炮管内等离子体密度计算方法
CN107480388A (zh) * 2017-08-23 2017-12-15 西安近代化学研究所 一种预定设计尺寸发射药的压伸成型方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6736942B2 (en) * 2000-05-02 2004-05-18 Johns Hopkins University Freestanding reactive multilayer foils
JP5986234B2 (ja) * 2015-01-15 2016-09-06 日油技研工業株式会社 焼尽性容器
CN107353170B (zh) * 2017-08-23 2018-10-12 西安近代化学研究所 一种发射药压伸工艺参数的动态计算方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
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
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 西安近代化学研究所 一种底部多层等离子体点火装置
CN106202797A (zh) * 2016-07-21 2016-12-07 中国人民解放军装甲兵工程学院 一种火药燃烧时炮管内等离子体密度计算方法
CN107480388A (zh) * 2017-08-23 2017-12-15 西安近代化学研究所 一种预定设计尺寸发射药的压伸成型方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Dependence of density and burning rate of composite solid propellant on mixer size;SumitVerma 等;《Acta Astronautica》;20140131;第93卷;130-137 *
球形变燃速发射药的燃气生成规律;刘海情 等;《广州化工》;20160308;第44卷(第5期);29-31 *

Also Published As

Publication number Publication date
CN110688771A (zh) 2020-01-14

Similar Documents

Publication Publication Date Title
Huang et al. Study of processing parameters in fused deposition modeling based on mechanical properties of acrylonitrile‐butadiene‐styrene filament
CN110688771B (zh) 一种基于渗透扩散的片状多层发射药燃烧性能计算方法
Trantina et al. Elastic-plastic analysis of small defects—voids and inclusions
Yi et al. Study of the key technologies of LOM for functional metal parts
JPS5916831B2 (ja) 膜構造型酸素センサ−の製造方法
CN110685651A (zh) 一种多层合采气井产量劈分方法及系统
Chakrabarti et al. A new triangular element to model inter‐laminar shear stress continuous plate theory
Pilipović et al. Influence of laser sintering parameters on mechanical properties of polymer products
Herle et al. Hypoplastic description of the interface behaviour
CN112818469A (zh) 固体火箭发动机映射设计方法、装置和设备
CN110765407B (zh) 一种片状多层发射装药内弹道特性参数的计算方法
CN109448039B (zh) 一种基于深度卷积神经网络的单目视觉深度估计方法
CN106886615B (zh) 一种rdx基多组分含能复合物的模拟方法
CN111680434B (zh) 包含裂纹尖端非奇异应力项的裂缝扩展路径数值计算方法
FR2461948A1 (fr) Procede de production d'une couche plate d'un electrolyte solide pour un capteur d'oxygene du type en pellicule plate
JP2008280190A (ja) アルミナ材料の焼成プロファイル調整方法及びそれを用いたセラミック積層体の製造方法
CN115619930A (zh) 一种基于参数化体积布尔运算的燃面估计方法
JP2010181252A (ja) セラミックシートの製造方法及びセラミック積層体の製造方法
Maggi et al. Burn-rate measurement on small-scale rocket motors
CN105840926A (zh) 一种深海油气工程用复合柔性管及其制造方法
AU2016319724A1 (en) Propellant charge
CN112507469A (zh) 固体火箭发动机燃烧室绝热层设计方法
Barskov et al. The study of the properties of ABS plastics and polyamide products, made by 3D printing for use in turbine construction
CN109948269A (zh) 一种铠装脐带缆模型六面体网格划分方法
CN114996787B (zh) 梯度功能材料球结构的应力分布确定方法及装置

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