CN107451344B - 一种固体推进剂燃速测试用靶线尺寸的确定方法 - Google Patents

一种固体推进剂燃速测试用靶线尺寸的确定方法 Download PDF

Info

Publication number
CN107451344B
CN107451344B CN201710595136.5A CN201710595136A CN107451344B CN 107451344 B CN107451344 B CN 107451344B CN 201710595136 A CN201710595136 A CN 201710595136A CN 107451344 B CN107451344 B CN 107451344B
Authority
CN
China
Prior art keywords
target line
target
tau
propellant
wire
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
CN201710595136.5A
Other languages
English (en)
Other versions
CN107451344A (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.)
Xian Modern Chemistry Research Institute
Original Assignee
Xian Modern Chemistry Research Institute
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 Xian Modern Chemistry Research Institute filed Critical Xian Modern Chemistry Research Institute
Priority to CN201710595136.5A priority Critical patent/CN107451344B/zh
Publication of CN107451344A publication Critical patent/CN107451344A/zh
Application granted granted Critical
Publication of CN107451344B publication Critical patent/CN107451344B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N31/00Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods
    • G01N31/12Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods using combustion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明公开了一种固体推进剂燃速测试用靶线尺寸的确定方法。在靶线法燃速测试过程中,靶线尺寸的选择直接关系到系统的计时误差和实验操作可靠性。本发明以燃速测试过程中靶线的熔断过程为基础,采用传热学公式描述靶线与燃气之间的对流传热过程,并通过公式求解获得满足计时准确性要求的最大靶线直径,从而形成燃速测试用靶线尺寸的确定方法。本发明弥补了当前靶线尺寸通过经验进行确定的不足,对提高靶线计时精度、增加实验操作可靠性、提高燃速测试标准化程度等皆有重要意义。

Description

一种固体推进剂燃速测试用靶线尺寸的确定方法
技术领域
本发明涉及一种固体推进剂燃速测试用靶线尺寸确定方法,尤其适用于低燃速、低燃烧温度推进剂燃速测试。
背景技术
固体推进剂是各类固体火箭发动机的推进能源,其性能参数对发动机性能水平的提高至关重要。燃速是固体推进剂的核心性能参数之一,燃速值的准确获得对火箭发动机的设计和推进剂配方调节均具有重要意义。固体推进剂的燃速只能通过实验方法准确获得,常用的燃速测试方法有恒压弹法、密闭燃烧器法和发动机法三种,其中恒压弹法由于具有测试结果可靠、成本低廉、数据处理方便等优点成为目前最为普及的测试方法。恒压弹法按计时方式的不同又可分为靶线法、声发射法、超声波法等,但考虑到信号的稳定性,靶线计时方法最为常用。靶线法计时原理是沿着推进剂药条的燃烧方向分别穿入两根低熔点金属丝,这些低熔点金属丝称为靶线。实验前准确测量出两靶线间的距离,由于推进剂的燃烧产生大量的热,因此当推进剂燃烧到靶线位置时会将其熔断。在实验过程中,可测得两靶线熔断的时间差,再根据靶线间距便可计算出推进剂的燃速。
燃速测试用的靶线一般采用铅锑合金丝,靶线的熔点为500~600K,而推进剂的燃面温度一般为700K以上,因此当推进剂药条烧到靶线位置时靶线很容易被烧断,从而实现准确计时。在国军标GJB 700B-2005《火药实验方法》方法706.1燃速-靶线法中对药条尺寸的规定为:双基推进剂为Φ5×150mm的圆柱体药条,复合推进剂为4×4×150mm的长方体药条,由于药条的细长特性,当推进剂的燃速较低时药条的燃烧时间较长,燃烧时存在较为严重的热散失。靶线法燃速测试过程中,要求靶线在燃面处能够被迅速熔断,若药条燃烧时的热损失较大,则会造成火焰对燃面热反馈的不足而引起的靶线不易烧断的现象,从而导致计时误差变大。靶线的熔断时间是计时系统准确性的重要标志,在其它条件相同的情况下,靶线的直径越小,靶线的升温越快,达到熔点的时间也就越短,因此减小靶线的直径有利于计时准确性的提高。但是铅锑合金材质的靶线强度较低,直径较小的靶线在装配过程中容易折断,从实验操作的角度上来看需要靶线的直径较大,因此需要一种最佳靶线尺寸的确定方法,使满足计时准确性的基础上靶线的尺寸尽量大。
发明内容
本发明提出一种固体推进剂燃速测试用靶线尺寸的确定方法,使靶线尺寸确定过程简单、准确性高、条件适应性强。
本发明的实现过程如下:
步骤1:设靶线体积为V,表面积为A,初始温度为T0,推进剂燃气温度为为T,T>T0,靶线与流体间的表面换热系数为h,靶线的比热容为c,靶线的密度为ρ,应用集中参数法可得到靶线温度随时间的变化关系;
初始时刻到某一瞬间为止的时间间隔内,靶线的瞬时热流量
Figure BDA0001355707940000021
引入过余温度θ=T-T,则式(a)变为
Figure BDA0001355707940000022
以过余温度表示的初始条件为
θ(0)=T-T0………………………………(c)
将式(b)分离变量,得
Figure BDA0001355707940000023
将式(d)对τ从0到τ积分,有
Figure BDA0001355707940000024
Figure BDA0001355707940000025
Figure BDA0001355707940000026
步骤2:将已穿靶线和点火线的推进剂药条安装至密封燃烧器中,开启恒压系统并点燃推进剂药条,当燃面退移至靶线处时,靶线与燃气之间通过对流传热进行热量交换,燃气沿着垂直于靶线轴线的方向流过靶线表面,通过外部强制对流传热的实验关联式
Figure BDA0001355707940000027
确定表面换热系数h,
其中Nu数通过邱吉尔与朋斯登准则式确定:
Figure BDA0001355707940000031
Re为雷诺数,Pr为普朗特数,λf为流体导热系数,d为靶线直径;
步骤3:在确定靶线的尺寸时,若计时准确度要求靶线的熔断时间不大于[τ],则先假设一个靶线直径,再根据药条尺寸确定埋入到药条内靶线的长度,根据数学模型计算出熔断时间τ,若τ=[τ],则靶线直径满足使用要求;若τ<[τ],则逐步加大靶线直径,最终使τ=[τ];若τ<[τ],则逐步减小靶线直径,最终使τ=[τ]。
本发明的优点:本发明通过传热学相关理论,对靶线在推进剂燃烧过程中的升温熔断过程进行计算,从而获得综合考虑计时精度和装配强度的最佳靶线尺寸。发明所涉及到的靶线尺寸确定过程简单、准确性高、条件适应性强,具有较好的经济效益。
具体实施方式
下面通过实施例对本发明做进一步解释说明。
本发明具体实施如下:
步骤1:根据所测推进剂的类型,制备与燃速测试样品具有相同截面尺寸的推进剂药条,并在药条中埋置直径为100μm的细钨铼热电偶丝,测试推进剂的燃烧波结。通过对药条燃烧波曲线的处理得到推进剂药条的燃面温度。
步骤2:根据燃速测试精度的要求,指定相应的靶线熔断时间[τ],并根据待测推进剂药条尺寸确定埋入到药条中的靶线长度。初步假定一个靶线直径d,进而可得到埋入推进剂药条中的靶线表面积A,然后通过数学模型确定靶线的熔断时间τ。
步骤3:若初算得到的靶线熔断时间τ=[τ],则可认为假设的靶线直径d为最佳直径;若τ<[τ],则逐步加大靶线直径,直至τ=[τ],此时的靶线直径便为最佳直径;若τ<[τ],则逐步减小靶线直径,直至τ=[τ],此时的靶线直径便为最佳直径。
实施例1
燃速待测推进剂为双基推进剂,将推进剂切成尺寸为Φ5×60mm的圆柱体药条,并测试药条在3MPa下的燃烧波曲线,通过燃烧波测试实验得到的推进剂燃面温度t=600K。燃速测试使用的靶线为熔点519K的铅锑合金丝,由推进剂药条尺寸可知埋入推进剂药条中的靶线长度为5mm。燃速测试要求靶线的最大熔断时间[τ]为1.7ms,由数学模型计算得到的该推进剂在3MPa燃速测试的最佳靶线直径为1.1mm。
实施例2
燃速待测推进剂为双基推进剂,将推进剂切成尺寸为Φ5×60mm的圆柱体药条,并测试药条在10MPa下的燃烧波曲线,通过燃烧波测试实验得到的推进剂燃面温度t=600K。燃速测试使用的靶线为熔点500K的铅锑合金丝,由推进剂药条尺寸可知埋入推进剂药条中的靶线长度为5mm。燃速测试要求靶线的最大熔断时间[τ]为4.1ms,由数学模型计算得到的该推进剂在10MPa燃速测试的最佳靶线直径为2.3mm。
实施例3
燃速待测推进剂为富燃料推进剂,将推进剂切成尺寸为4×4×60mm的长方体药条,并测试药条在8MPa下的燃烧波曲线,通过燃烧波测试实验得到的推进剂燃面温度t=700K。燃速测试使用的靶线为熔点500K的铅锑合金丝,由推进剂药条尺寸可知埋入推进剂药条中的靶线长度为4mm。燃速测试要求靶线的最大熔断时间[τ]为4.7ms,由数学模型计算得到的该推进剂在8MPa下燃速测试的最佳靶线直径为3.3mm。
实施例4
燃速待测推进剂为富燃料推进剂,将推进剂切成尺寸为4×4×60mm的长方体药条,并测试药条在0.2MPa下的燃烧波曲线,通过燃烧波测试实验得到的推进剂燃面温度t=711K。燃速测试使用的靶线为熔点500K的铅锑合金丝,由推进剂药条尺寸可知埋入推进剂药条中的靶线长度为4mm。燃速测试要求靶线的最大熔断时间[τ]为2.9ms,由数学模型计算得到的该推进剂在0.2MPa下燃速测试的最佳靶线直径为2.1mm。

Claims (1)

1.一种固体推进剂燃速测试用靶线尺寸的确定方法,其特征在于步骤如下:
步骤1:设靶线体积为V,表面积为A,初始温度为T0,推进剂燃气温度为T,T>T0,靶线与流体间的表面换热系数为h,靶线的比热容为c,靶线的密度为ρ,某一瞬时时刻τ靶线温度为T,应用集中参数法可得到靶线温度随时间的变化关系;
初始时刻到某一瞬间为止的时间间隔内,靶线的瞬时热流量
Figure FDA0002493735720000011
引入过余温度θ=T-T,则式(a)变为
Figure FDA0002493735720000012
以过余温度表示的初始条件为
θ0=T-T0……………………………………(c)
将式(b)分离变量,得
Figure FDA0002493735720000013
将式(d)对τ从0到τ积分,得
Figure FDA0002493735720000014
式(e)积分结果为
Figure FDA0002493735720000015
对式(f)等号两边取对数,得
Figure FDA0002493735720000016
步骤2:将已穿靶线和点火线的推进剂药条安装至密封燃烧器中,开启恒压系统并点燃推进剂药条,当燃面退移至靶线处时,靶线与燃气之间通过对流传热进行热量交换,燃气沿着垂直于靶线轴线的方向流过靶线表面,通过外部强制对流传热的实验关联式
Figure FDA0002493735720000017
确定表面换热系数h,
其中Nu数通过邱吉尔与朋斯登准则式确定:
Figure FDA0002493735720000021
Re为雷诺数,Pr为普朗特数,λf为流体导热系数,d为靶线直径;
步骤3:在确定靶线的尺寸时,若计时准确度要求靶线的熔断时间不大于[τ],则假设一个靶线的直径,再根据药条尺寸确定埋入到药条内靶线的长度,根据步骤1和步骤2建立的数学模型计算出熔断时间τm,若τm=[τ],则靶线直径满足使用要求;若τm<[τ],则逐步加大靶线直径,最终使τm=[τ];若τm<[τ],则逐步减小靶线直径,最终使τm=[τ]。
CN201710595136.5A 2017-07-20 2017-07-20 一种固体推进剂燃速测试用靶线尺寸的确定方法 Active CN107451344B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710595136.5A CN107451344B (zh) 2017-07-20 2017-07-20 一种固体推进剂燃速测试用靶线尺寸的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710595136.5A CN107451344B (zh) 2017-07-20 2017-07-20 一种固体推进剂燃速测试用靶线尺寸的确定方法

Publications (2)

Publication Number Publication Date
CN107451344A CN107451344A (zh) 2017-12-08
CN107451344B true CN107451344B (zh) 2020-07-17

Family

ID=60488747

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710595136.5A Active CN107451344B (zh) 2017-07-20 2017-07-20 一种固体推进剂燃速测试用靶线尺寸的确定方法

Country Status (1)

Country Link
CN (1) CN107451344B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112623270B (zh) * 2020-12-23 2022-08-09 北京空天技术研究所 一种基于高速摄像和靶线的高速分离测量装置及方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102175830A (zh) * 2011-02-17 2011-09-07 西安近代化学研究所 固体推进剂多靶线准动态燃烧性能测试方法
CN107132308A (zh) * 2017-05-15 2017-09-05 西北工业大学 一种基于燃面标识的固体推进剂燃速辨识方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102175830A (zh) * 2011-02-17 2011-09-07 西安近代化学研究所 固体推进剂多靶线准动态燃烧性能测试方法
CN107132308A (zh) * 2017-05-15 2017-09-05 西北工业大学 一种基于燃面标识的固体推进剂燃速辨识方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种小型熔断体熔体直径的计算;张叔韬;《低压电器》;19861028;全文 *
熔断电阻器瞬态特性仿真分析及其性能的研究;邱文文;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20140115;第32-36页 *

Also Published As

Publication number Publication date
CN107451344A (zh) 2017-12-08

Similar Documents

Publication Publication Date Title
CN107132308B (zh) 一种基于燃面标识的固体推进剂燃速辨识方法
CN101581649B (zh) 固体推进剂燃烧温度测试装置
CN103439360B (zh) 固体推进剂多热电偶动态燃烧性能测试系统和方法
CN108287030B (zh) 一种内埋式热电偶表面热流测量方法
CN107451344B (zh) 一种固体推进剂燃速测试用靶线尺寸的确定方法
CN113343495B (zh) 一种管壳式燃滑油散热器热性能模型修正方法
CN104655360A (zh) 内燃机的缸内压力推定装置
Demuynck A fuel independent heat transfer correlation for premixed spark ignition engines
CN108226222A (zh) 一种用于测试可燃气体燃爆参数的装置及方法
CN113047981B (zh) 冲量法固体推进剂燃速测试原始实验数据有效性判定方法
CN109738196B (zh) 一种基于地面直连试验的冲压发动机性能天地换算方法
Sieder et al. Development history and verification of the flight model of a 500 N Ethanol/LOX rocket engine
CN111006870A (zh) 一种定容燃烧传热装置
CN109737836A (zh) 一种弹体侵彻过载快速评价方法
CN203365362U (zh) 固体推进剂多热电偶动态燃烧性能测试的燃烧室装置
CN115221638A (zh) 亚燃冲压发动机非稳态过程的性能时间响应分析方法
CN109632613B (zh) 一种压装侵彻炸药JHLDu-1装药安全使用寿命的预估方法
CN113790899A (zh) 液体火箭发动机喷注器燃烧状态评估装置及方法
CN114034731A (zh) 含能材料爆热、爆压的双功能测定装置和方法
CN108182297B (zh) 火箭发动机燃烧室效率分析方法及系统
CN111046531A (zh) 一种htpe型固体发动机烤燃数值模拟计算方法
CN106093284B (zh) 一种金属燃料燃速的快速测试系统及方法
Zhang Piston temperature measurement using voltage recorder and numerical simulation of the temperature field
Petkovic et al. Experimental verification of mathematical model of the heat transfer in exhaust system
RU2701522C1 (ru) Манометрическая бомба высокого давления

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